• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    鏈霉菌Streptomyces sp. FXP04全基因組測序分析

    2023-11-23 09:11:40唐碧瑤付學(xué)鵬
    生物技術(shù)通報 2023年10期
    關(guān)鍵詞:基因簇霉菌基因組

    唐碧瑤 付學(xué)鵬

    (齊齊哈爾大學(xué)生命科學(xué)與農(nóng)林學(xué)院,齊齊哈爾 161006)

    鏈霉菌是具有復(fù)雜生長周期的有機化能、好氧型革蘭氏陽性放線菌(G+),在系統(tǒng)分類上歸為放線菌目,是放線菌中最為高等的一類菌群。鏈霉菌種類繁多,在陸地、海洋、極端環(huán)境以及一些生物體內(nèi)都有廣泛分布,其中主要分布在土壤環(huán)境中。其孢子抗逆性強、存活時間長,具有良好的分枝菌絲,能產(chǎn)生水溶性或脂溶性的色素。因為其復(fù)雜的生長發(fā)育周期,產(chǎn)生了多種具有復(fù)雜結(jié)構(gòu)的天然代謝產(chǎn)物,其中產(chǎn)生大量的各式各樣的抗生素,目前在臨床應(yīng)用的抗生素約 2/3 來源于鏈霉菌屬[1],包括井崗霉素、萬古霉素等。除抗生素外還包括抗腫瘤劑、殺蟲劑、酶抑制劑、色素及酶類,如幾丁質(zhì)酶、果膠酶、木聚糖酶和蛋白酶等[2]。例如Arasu 等[3]從海洋中分離出Streptomyces sp. AP?123鏈霉菌產(chǎn)生的聚酮類化合物對綠猴腎細(xì)胞和喉癌細(xì)胞有顯著毒性。Kim等[4]研究發(fā)現(xiàn)沙門氏鏈霉菌菌株KP1404能產(chǎn)生多烯類抗生素,抑制番茄枯萎病菌菌絲的生長,其拮抗效果與苯菌靈相當(dāng)。胡磊等[5]發(fā)現(xiàn)產(chǎn)水鏈霉菌和球孢鏈霉菌可明顯抑制由核盤菌引起的油菜菌核病,球孢鏈霉菌的防治機理可能是產(chǎn)生幾丁質(zhì)酶部分分解油菜菌核。這些產(chǎn)物在醫(yī)療、農(nóng)業(yè)、食品、化工、環(huán)保等領(lǐng)域都有巨大的發(fā)展?jié)摿Γ瑢θ祟惖纳a(chǎn)生活至關(guān)重要。因此,近幾十年來,關(guān)于鏈霉菌遺傳發(fā)育及代謝調(diào)控的研究備受重視。

    Streptomyces sp. FXP04是本課題組在健康馬鈴薯根際土壤中分離的鏈霉菌菌株。該菌對致病疫霉P.infestans有很強的拮抗效果,對于作物生長具有良好的防病和促生[6-7]作用,表現(xiàn)出了較好的應(yīng)用前景,有作為生物農(nóng)藥的潛力。近年來,對于鏈霉菌的研究多集中在抗菌、促生、抗腫瘤等生物活性以及分離純化的表層上,而利用基因組測序技術(shù)研究和預(yù)測微生物代謝產(chǎn)物的種類、結(jié)構(gòu)和合成途徑成為了發(fā)掘鏈霉菌屬代謝產(chǎn)物的一種新方向。在自然情況下,鏈霉菌次級代謝產(chǎn)物的實際表達(dá)率并不高,應(yīng)用效果大大縮減。所以利用全基因組分析有助于人們對這種具有多種次級代謝物生產(chǎn)能力的鏈霉菌進行更深層次的研究。全基因組測序是針對未知基因組序列的物種進行個體全部基因組測序,具有測序覆蓋面廣、準(zhǔn)確性極高的技術(shù)優(yōu)勢[8]。在本研究中,利用Illumina HiSeq高通量測序平臺對FXP04的全基因組進行了測序。同時,進行了基因預(yù)測和基因序列的功能注釋以及次級代謝產(chǎn)物合成基因簇預(yù)測,利用生物信息學(xué)分析強調(diào)了天然產(chǎn)物生物合成在菌株應(yīng)用中的重要性,比較基因組學(xué)為FXP04的進化提供了見解,為揭示FXP04生防作用機制、功能基因的挖掘與利用奠定基礎(chǔ)。

    1 材料與方法

    1.1 Streptomyces sp. FXP04培養(yǎng)和基因組DNA的提取

    FXP04菌株接種至Luria?Bertani(LB)培養(yǎng)基活化,28℃培養(yǎng) 4 d,挑取菌絲轉(zhuǎn)接至三角瓶液體培養(yǎng)基于 28℃、140 r/min 培養(yǎng) 4 d后,12 000 r/min 離心20 min 收集菌體送至華大基因公司進行de novo 測序分析。FXP04已于中國普通微生物菌種保藏管理中心進行保藏,保藏號為CGMCC No.16826。

    1.2 全基因組測序

    測序采用第三代 PacBio 平臺與第二代 Illumina平臺相結(jié)合的測序技術(shù),得到的測序結(jié)果采用Glim?mer軟件對基因組進行編碼基因預(yù)測;TRF(Tandem Repeat Finder)軟件對串聯(lián)重復(fù)序列預(yù)測;通過軟件PHAST預(yù)測前噬菌體(Prophage);利用相關(guān)軟件完成對非編碼RNA的預(yù)測,包括RNAmmer軟件預(yù)測rRNA、tRNAscan軟件預(yù)測tRNA區(qū)域和tRNA的二級結(jié)構(gòu)以及利用Infernal軟件與Rfam數(shù)據(jù)庫比對得到sRNA。測序結(jié)果再與GO(Gene Ontology)、KEGG(Kyoto Encyclopedia of Genes and Genomes)、COG(Cluster of Orthologous Groups of Proteins)、VFDB(Virulence Factor Database)、ARDB(Antibiotic Resistance Genes Database)、CAZy(Carbohydrate?Active EnZymes Database)數(shù)據(jù)庫進行比對,完成蛋白序列功能注釋[9]。

    1.3 比較基因組學(xué)分析

    先通過FXP04的16S RNA序列比對,從中選取5種與其近緣關(guān)系相近且不同類型的模式菌株,從NCBI Genome上下載S. varsoviensis NRRL ISP?5346(GCF_000718635.1),S. diacarni LHW51701(GCF_003323715.1),S. celluloflavus NRRL B?2493(GCF_000720995.1),S. kasugaensis BCRC 12349(GCF_002261115.1),S. platensis DSM 40041(GCF_00211 9195.1)共5種模式菌株的全基因組序列,與FXP04全基因組序列一起進行以下比較基因組學(xué)分析。其中包括共有基因(Core Gene)和特有基因(Specif?ic Gene)分析[10]、基因家族(Gene Family)分析[11-13],通過相關(guān)軟件進行聚類分析,根據(jù)基因集分布情況,繪制熱圖和構(gòu)建系統(tǒng)進化樹。

    1.4 次生代謝產(chǎn)物合成基因簇分析

    將菌株FXP04的全基因序列采用antiSMASH(V6.1.1)軟件(https://antismash.secondarymetabolites.org/#!/ start)在線預(yù)測次級代謝產(chǎn)物合成基因簇,并結(jié)合 NCBI BLAST 比對分析結(jié)果[14]。

    2 結(jié)果

    2.1 菌株FXP04全基因組概況

    本研究利用PacBio與Illumina 平臺相結(jié)合測序技術(shù)對FXP04基因組測序,共得到1 367 Mb的原始數(shù)據(jù)量,對原始數(shù)據(jù)過濾后得到1 141 Mb的數(shù)據(jù)量,K?mer預(yù)估其基因組大小為4.89 Mb(K?mer值為15,深度為37)。FXP04基因組大小為4 535 201 bp(圖1)[6],共有5 037個基因,基因的總長度為3 674 928 bp,GC含量為72.95%,基因長度占基因組總長的81.03%。在菌株 FXP04的核基因組中,利用RNAmmer軟件[15]和tRNAscan軟件[16]對FXP04基因組的ncRNA進行預(yù)測分析,其中 rRNA有 3 種類型:3個5S rRNA、3個16S rRNA、3個23S rRNA,33個tRNA以及23個sRNA;利用 TRF軟件[17]對FXP04基因組的重復(fù)序列進行搜尋,結(jié)果顯示FXP04共有1 699個串聯(lián)重復(fù)序列,總長度為77 244 bp,大小為4-672 bp,其中含1 284個小衛(wèi)星 DNA,含188個微衛(wèi)星DNA;利用 PHAST軟件[18]對基因組進行前噬菌體預(yù)測的結(jié)果僅得到一個可疑的前噬菌體,其長度為34 291 bp。

    圖1 菌株FXP04的全基因組Fig. 1 Whole genome of strain FXP04

    2.2 基因功能注釋

    2.2.1 VFDB數(shù)據(jù)庫注釋 將菌株FXP04的基因組與VFDB數(shù)據(jù)庫進行比對分析,在FXP04的全基因組中共發(fā)現(xiàn)139個毒力相關(guān)基因,其中最大基因序列為1 999 bp,最小基因序列為41 bp(表1)。通過比對發(fā)現(xiàn)基因序列主要在PDIM(二霉酚)和PGL(酚類糖脂)的生物合成和運輸(ddrA、pks15、ppsC基因)、ABC轉(zhuǎn)運蛋白(irtA、fbpC、fagA/B/C/D基因)、DevR/S雙組分系統(tǒng)(devR/dosR基因)、PhoP/R雙組分系統(tǒng)(phoP/R基因)4種功能中有較高的覆蓋度。ddrA基因編碼ABC轉(zhuǎn)運蛋白參與抗腫瘤藥物柔紅霉素合成,pks15基因編碼聚酮合酶,ppsC基因編碼β-酮?;厦福籪bpC、fagA/B/C/D基因與鐵轉(zhuǎn)運系統(tǒng)相關(guān),irtA基因編碼跨膜ATP結(jié)合蛋白幫助藥物轉(zhuǎn)運進入細(xì)胞;devR/dosR基因和phoP/R基因都為雙組分系統(tǒng)調(diào)控因子,綜上說明菌株FXP04所注釋到的139個基因序列并非真正的毒力基因,且identity(一致性)≤70%,所以我們推測FXP04在與病原菌對抗中通過控制某些基因調(diào)節(jié)自身或共生生物體的生物學(xué)過程來發(fā)揮作用,也表明了FXP04被用作生物防治時是安全無害的。

    2.2.2 ARDB數(shù)據(jù)庫注釋 通過ARDB數(shù)據(jù)庫的注釋可了解FXP04基因的潛在的抗生素耐受性,通過注釋結(jié)果找出了14個耐藥基因,基因長度介于100-480 bp之間,耐受藥物嘌呤霉素(Puromycin)、桿菌肽(Bacitracin)、氨基糖苷(Aminoglycoside)、四烯霉素c(Tetracenomycin_c)、抗菌物(Na_Antimicrobials)、氯霉素(Chloramphenicol)、鏈霉素(Streptomycin)、萬古霉素(Vancomycin)共8種,具體預(yù)測結(jié)果見表2。

    表2 菌株FXP04在ARDB數(shù)據(jù)庫注釋結(jié)果Table 2 Annotation results of strain FXP04 in ARDB database

    2.2.3 CAZy 數(shù)據(jù)庫注釋 將基因組序列與 CAZy 數(shù)據(jù)庫進行比對,發(fā)現(xiàn)菌株 FXP04的基因組中共有73個編碼蛋白結(jié)構(gòu)域歸于 CAZy 家族,其蛋白編碼基因分布狀況如下:糖苷水解酶(glycoside hydrolases,GHs)基因27個,糖類酯解酶(car bohydrate esterases, CEs)基因6個,糖基轉(zhuǎn)移酶(glycosyltransferases, GTs)基因12個,碳水化合物酶相關(guān)模塊(carbohydrate?binding modules, CBMs)基因18個,輔助功能酶(auxiliary activities, AAs)基因7個,多糖裂合酶(polysaccharide lyases, PLs)基因3個。在糖苷水解酶中,覆蓋度最高的主要是葡聚糖酶(EC 3.2.1.70)。在碳水化合物酶相關(guān)模塊中,幾丁質(zhì)酶(EC 3.2.1.14)、β-1,4?木聚糖內(nèi)切酶(EC 3.2.1.8)等有較高的覆蓋度。在糖基轉(zhuǎn)移酶中,磷壁酸酶(EC 2.4.1.129)有著較高覆蓋度。在糖類酯解酶中D?葡萄糖醛酸酶(EC 3.1.1.?)也有比較高覆蓋度。

    2.2.4 GO數(shù)據(jù)庫注釋 FXP04通過GO注釋,共有2 088個基因在GO數(shù)據(jù)庫中被注釋到,按照GO數(shù)據(jù)庫的分類方式,F(xiàn)XP04的基因組功能被分為31個分支(圖2)。其中有9個跟細(xì)胞組分有關(guān),12個與生物學(xué)過程有關(guān),10個與分子功能有關(guān)。在細(xì)胞組分中共有731個基因被注釋,其中與細(xì)胞膜、細(xì)胞膜部分有關(guān)的基因表現(xiàn)出最高相關(guān)性,各有306、160個;在生物學(xué)過程中共3 867個基因被注釋,與其細(xì)胞過程、代謝過程以及單有機體過程的基因表現(xiàn)出很高的相關(guān)性,各有831、1 193、752個;分子功能分支共2 531個注釋結(jié)果,其中涉及最多的基因與結(jié)合(binding)、催化活性(catalytic activity)有關(guān),分別為925、1 261個。

    圖2 菌株FXP04基因組的GO功能分類圖Fig. 2 GO functional classification map of strain FXP04 genome

    2.2.5 COG數(shù)據(jù)庫注釋 將FXP04基因組中能編碼蛋白的基因與COG進行比對,得到FXP04的COG 功能數(shù)據(jù)分析結(jié)果如圖3所示。在 FXP04菌株有2 319個可以在 COG 數(shù)據(jù)庫中找到分類信息,包含COG數(shù)據(jù)庫功能的23類。主要涉及 General function prediction only(通用功能預(yù)測)345個,占功能注釋基因的14.88%;Transcription(轉(zhuǎn)錄)371個,占功能注釋基因的16%;Amino acid transport and metabolism(氨基酸的運輸與代謝)229個,占功能注釋基因的9.87%;Carbohydrate transport and metabolism(碳水化合物的運輸與代謝)205個,占功能注釋基因的8.84%;Signal transduction mechanisms(信號轉(zhuǎn)導(dǎo)結(jié)構(gòu))198個,占功能注釋基因的8.54%,不同類型的基因執(zhí)行功能有所側(cè)重。

    圖3 菌株FXP04基因組的COG功能分類圖Fig. 3 COG functional classification diagram of the genome of strain FXP04

    2.2.6 KEGG數(shù)據(jù)庫注釋 FXP04菌株共有1 530個基因在KEGG中被注釋,被分為6個分類31個功能條目(圖4)?;驍?shù)量占比最多的是與新陳代謝相關(guān)的。在KEGG注釋的代謝中,參與概述的基因有535個,輔因子和維生素的代謝130個,氨基酸代謝182個,碳水化合物代謝175個。

    圖4 菌株FXP04基因組的KEGG功能分類圖Fig. 4 KEGG functional classification diagram of FXP04 genome

    2.3 比較基因組學(xué)分析

    2.3.1 共有基因和特有基因分析 通過 Core/Pan分析發(fā)現(xiàn),F(xiàn)XP04與5個模式菌株的所有非冗余的 Pan Gene 數(shù)量為18 057個,長度為 5 710 958 bp;其中共有的Core Gene 數(shù)量為532 個,長度為180 719 bp;非共有的Dispensable Gene 數(shù)量為6 778個(圖5?B),長度為2 315 398 bp。通過Dispensable Gene同源關(guān)系熱圖(圖5?A)和Core Gene系統(tǒng)發(fā)育樹(圖5?C)可看出,F(xiàn)XP04的基因組與S. varsoviensis NRRLISP?5346最為相似,其次是S. diacarni LHW51701和S. platensis DSM40041,與另外2 株模式菌株具有明顯的遺傳距離。

    圖5 6個菌株之間基因組Core/Pan 分析結(jié)果Fig. 5 Genome Core/Pan analysis results among 6 strains

    2.3.2 基因家族分析 基于6個菌株全基因組序列的 Gene Family分析發(fā)現(xiàn),F(xiàn)XP04與S. varsoviensis NRRLISP?5346、S. diacarni LHW51701、S. platensis DSM40041等 6個菌株能夠聚類到Gene Family 的基因數(shù)量分別為 3 010 個、5 827個、5 287 個、5 974個、6 692個和6 376個,Gene Family 數(shù)量分別為 1 437個、2 469 個、2 294個、2 613 個、3 101個和3 019個;其中共有的 Gene Family 數(shù)量為 871個,各菌株特有的 Gene Family 數(shù)量分別為 25 個、25 個、24個、40個、9個、3個(圖6?A、圖6?B)。通過 Gene Fam?ily系統(tǒng)發(fā)育樹可看出,F(xiàn)XP04 菌株與 S. varsoviensis NRRLISP?5346聚類到同一分支,其次是S. diacarni LHW51701和S. platensis DSM40041,與另外2株菌株之間的遺傳距離約為0.45(圖6?C)。通過Gene Family和Core Gene系統(tǒng)發(fā)育樹發(fā)現(xiàn)FXP04 菌株與S. varsoviensis NRRLISP?5346模式菌株的相似性最高,選取S. varsoviensis NRRLISP?5346的基因序列與FXP04進行共線性分析(圖7),結(jié)果與系統(tǒng)發(fā)育樹的結(jié)果相同,但也存在倒位、缺失和重排等現(xiàn)象。

    圖6 6個菌株之間基因組Gene family分析結(jié)果Fig. 6 Genome Gene family analysis results among 6 strains

    圖7 FXP04與S. varsoviensis NRRLISP-5346共線性分析Fig. 7 Collinearity analysis between FXP04 and S.varsoviensis NRRLISP-5346

    2.3.3 次生代謝產(chǎn)物基因簇分析 利用 anti SMASH(V6.1.1)[14]對FXP04基因組的次級代謝產(chǎn)物合成基因簇進行預(yù)測分析,結(jié)果(表3)表明該菌株基因組中預(yù)測到了13個次級代謝產(chǎn)物合成相關(guān)基因簇,分別為聚酮合酶(polyketide synthase,PKS)基因簇3種、非核糖體多肽合成酶(non?ribosomal peptide synthase,NRPS)基因簇3種和萜烯類簇(Terpene)2種、芳基聚烯烴(arylpolyene)基因簇2種、套索肽(lassopeptide)基因簇1種、四氫嘧啶(ectoine)基因簇1種和尚未被定義出類型的基因簇1種。進一步將FXP04預(yù)測的13個基因簇與其中已知合成基因簇序列進行 BLAST 比對后發(fā)現(xiàn),其中有些與參考基因簇表現(xiàn)為完全相似或高度相似的合成基因簇。Cluster 3與BGC0001742(Candidatus Streptomyces philanthi bv. Triangulum)來源的殺粉蝶菌素(piericidin A1)基因簇相似度達(dá)到100%;Cluster 8與BGC0002052(Streptomyces sp. ID38640)來源的四氫嘧啶(ectoine)基因簇相似度也達(dá)到100%; Cluster 13與BGC0002547來源的Youssoufene類合成基因簇的相似度達(dá)到了88%。據(jù)相關(guān)的文獻(xiàn)報道,基因簇分析中確定已知化合物的標(biāo)準(zhǔn)是基因簇與已報道基因簇之間的相似度≥85%。所以,上述基因簇與已知基因簇序列的相似度≥85% 時,由該基因編碼形成的化合物是已知的;即低于85% 時,可以推測出該基因中存在獨特的序列。如Cluster 5與BGC0000663(Streptomyces coelicolor A3(2))來源的何帕烯(Hopene)基因簇的相似度達(dá)到61%;Cluster 7與BGC0000368來源的Griseobactin基因簇相似度達(dá)到53%;Cluster 4與BGC0001110(Uncultured bacterium psy1)來源的羊海綿抑素(Psymberin/Irciniastatin B)基因簇的相似度達(dá)到31%;Cluster 1與BGC0002662來源的二聚化Piperazyl環(huán)肽Petrichorin A /Petrichorin B基因簇的相似度達(dá)到24%;Cluster 6與 BGC0000431來源的Stenothricin基因簇和Cluster 11與BGC0002591來源的Aurachin類基因簇的相似度都達(dá)到13%;Cluster 9與BGC0001764來源的s56?p1基因簇相似度達(dá)到11%,這些基因簇可能會合成新的抑菌物質(zhì),有待于后續(xù)研究。除了上述的這些基因簇和一個尚未被定義出類型的基因簇,還有2個功能未知的基因簇(Cluster2、10)存在,說明菌株FXP04可能存在能合成新型抑菌物質(zhì)的潛力,對農(nóng)業(yè)的發(fā)展需求有著重要意義。

    表3 次生代謝產(chǎn)物合成區(qū)域鑒定結(jié)果Table 3 Identification results of synthetic regions of secondary metabolites

    3 討論

    Streptomyces sp. FXP04是從健康馬鈴薯根際土壤中分離的一株對致病疫霉具有明顯拮抗作用的鏈霉菌菌株[6],且前期研究也發(fā)現(xiàn)其對水稻等作物有促生效果[7]。日常生活中常見的抗生素大多來自鏈霉菌,而傳統(tǒng)實驗和鑒定方法很難全面地分析鏈霉菌的作用機制,對此,深入研究FXP04在基因組層面的內(nèi)在原因有著重要意義。所以我們通過全基因組測序和生物信息學(xué)分析得到了FXP04的全基因組序列,確定了FXP04基因組大小為4 535 201 bp,共編碼5 037個基因,通過GO、KEGG、COG等數(shù)據(jù)庫比對分析,完成了FXP04基因組的注釋結(jié)果和數(shù)據(jù)整理。結(jié)果表明無論在GO、COG還是KEGG 數(shù)據(jù)庫注釋中,菌株FXP04的基因功能多集中于氨基酸的轉(zhuǎn)運及代謝和碳水化合物的轉(zhuǎn)運及代謝。此外,在KEGG注釋中還發(fā)現(xiàn)53個基因與萜類和聚酮化合物的代謝相關(guān),可能參與抗菌活性物質(zhì)的生物合成,例如,ncsB1(Gene636)、ncsB3(Gene2 206)為烯二炔類抗生素(biosynthesis of enediyne antibiotics)的生物合成基因[19];tktA(Gene751)為芳香族氨基酸生物合成中心途徑的關(guān)鍵酶基因,糖肽類抗生素(glycopeptide antibiotics)幾乎都含有芳香族氨基酸,如萬古霉素等[20-21]。

    碳水化合物在動物生長生產(chǎn)方面具有重要作用。在碳水化合物功能注釋結(jié)果中顯示,菌株FXP04基因組中含有與葡聚糖酶、幾丁質(zhì)酶、β-1,4?木聚糖內(nèi)切酶、D?葡萄糖醛酸酶等相關(guān)的基因。β-1,4?木聚糖內(nèi)切酶、D?葡萄糖醛酸酶屬于木聚糖降解酶家族,木聚糖是僅次于纖維素存在于植物細(xì)胞壁的一種重要的異源多糖。葡聚糖酶屬纖維素酶家族,對纖維素的分解有著重要意義,幾丁質(zhì)廣泛存在于甲殼類動物的外殼、真菌的胞壁以及一些高等植物的細(xì)胞壁等,而菌株FXP04中存在的幾丁質(zhì)酶可以降解和破壞病原菌的外壁,從而起到防治作用,推測可能為其抑菌機理之一。

    利用antiSMASH工具對菌株FXP04的次級代謝產(chǎn)物基因簇進行預(yù)測分析,發(fā)現(xiàn)了10種已知的抗菌活性物質(zhì)合成基因簇,Hopene是合成具有巨大商業(yè)價值的生物活性物質(zhì)Hopanoids(何帕烷類)的重要前體,據(jù)相關(guān)報道何帕烷類化合物具有抗炎、抗菌等生理活性[22-24];Psymberin也稱為Irciniastatin A,是一種聚酮酸的生物活性天然產(chǎn)物。Psymberin于2004年由Crews和Pettit的研究小組從海綿Psammocinia sp.中獨立發(fā)現(xiàn),是一種有效的抗癌劑[25];Petrichorin A是一種具有抗癌、抗腫瘤活性的二聚化Piperazyl環(huán)肽[26];Griseobactin是一種具有兒茶酚肽結(jié)構(gòu)的鐵載體,在2010年P(guān)atzer等[27]發(fā)現(xiàn)由鏈霉菌屬ATCC 700974和幾種灰鏈霉菌株合成的,參與灰霉菌素的合成需求;s56?p1是一種具有獨特的肼單元的二肽天然產(chǎn)物,s56?p1生物合成基因簇是第一個被確定為負(fù)責(zé)含肼類天然產(chǎn)物生物合成的基因簇[28];Aurachin是各種細(xì)胞色素復(fù)合體的抑制劑,具有強大的抗菌、抗真菌和抗瘧原蟲藥物[29-30]。上述7種活性物質(zhì)與已知基因簇相似度低,推測可能合成新的活性物質(zhì),對于鏈霉菌屬的潛在的天然產(chǎn)物開采提供理論,對于農(nóng)業(yè)、醫(yī)學(xué)方面的應(yīng)用有著重要的研究意義。同時,通過預(yù)測得到的已確定的3種活性物質(zhì),包括Youssoufene、四氫嘧啶、殺粉蝶菌素,四氫嘧啶(Ectoine)可以幫助微生物適應(yīng)在高鹽、高滲透壓和紫外輻射等逆境中正常生長[31],使得促生菌株在惡劣氣候下能在根系或根際穩(wěn)定繁殖。Youssoufene化合物是含有肉桂酰基團的天然產(chǎn)物,對于Youssoufene化合物的報道較少,但肉桂?;惢衔镉兄鴱V泛的生理活性,包括抗菌和抗結(jié)核活性等[32]。殺粉蝶菌素(Piericidin A1)最早由日本科學(xué)家發(fā)現(xiàn),是一種具有良好生物活性的 α-吡啶酮類天然產(chǎn)物,后來在多個鏈霉菌中都有發(fā)現(xiàn)[33-34]。其可能是線粒體呼吸作用抑制劑,對多種生物具有較強的生物活性。首先,殺粉蝶菌素 A1 對蠶、菜青蟲有很顯著的殺蟲效果[35],對多種癌細(xì)胞和腫瘤細(xì)胞也有很好的生長抑制作用[36-37]。其次,殺粉蝶菌素對多種病原真菌和粘紅酵母等也有很強的抑制活性[38]。最后,何亞文等[39]發(fā)現(xiàn)殺粉蝶菌素A1 對水稻白葉枯菌的生長有較強的抑制作用,在溫室條件下噴施殺粉蝶菌素A1能夠減輕白葉枯癥狀。而且在我們前期研究中已經(jīng)對FXP04的次生代謝物進行純化,與殺粉蝶菌素標(biāo)品對比出的抑菌效果相同[6],進一步說明菌株FXP04有抗菌、抗蟲、抗腫瘤以及促生長的潛力。

    通過以上的數(shù)據(jù)結(jié)果分析,我們了解了菌株FXP04的全部基因組序列信息,發(fā)現(xiàn)其屬于天然小基因組,其中存在有3種重要的次級代謝產(chǎn)物基因簇和7種有望合成新型抑菌物質(zhì)的抗生素基因簇,極可能挖掘出結(jié)構(gòu)新穎、活性更強、靶向性更好的天然藥物用于臨床研究;而且據(jù)相關(guān)報道表明鏈霉菌基因組中含有大量的冗余基因,其中就包括次級代謝合成基因簇、可移動元件等[40],它們是菌株正常生長發(fā)育中非必需的,可以通過大規(guī)模刪除冗余基因用于構(gòu)建高版本工業(yè)底盤細(xì)胞,用作天然產(chǎn)物基因的異源表達(dá);其次對于殺粉蝶菌素的抑菌效果有了初步了解,后期將通過預(yù)測到的殺粉蝶菌素的合成途徑進行超表達(dá)和構(gòu)建缺失突變體,以驗證其在FXP04中的地位,同時也為通過遺傳技術(shù)獲取高產(chǎn)殺粉蝶菌素的菌株提供理論基礎(chǔ),有待于對Streptomyces sp. FXP04的次級代謝產(chǎn)物進行后續(xù)研究與應(yīng)用。

    4 結(jié)論

    通過Illumina和PacBio平臺相結(jié)合的測序技術(shù),對鏈霉菌Streptomyces sp. FXP04進行全基因組測序,獲得了菌株FXP04的全部基因組序列信息。分析結(jié)果可知FXP04基因組大小僅為4.5 M,屬天然的小基因組,后續(xù)可以構(gòu)建成高版本工業(yè)底盤細(xì)胞用作天然產(chǎn)物基因的異源表達(dá),滿足微生物藥物工業(yè)化生產(chǎn)需求,發(fā)揮其應(yīng)用潛力。

    猜你喜歡
    基因簇霉菌基因組
    揭示水霉菌繁殖和侵染過程
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    冬瓜高通量轉(zhuǎn)錄組測序及分析
    霉菌的新朋友—地衣
    地衣和霉菌
    腸球菌萬古霉素耐藥基因簇遺傳特性
    遺傳(2015年5期)2015-02-04 03:06:55
    海洋稀有放線菌 Salinispora arenicola CNP193 基因組新穎PKS 和NRPS基因簇的發(fā)掘
    青霉菌柚苷酶的分離純化及酶學(xué)性質(zhì)
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    亚洲国产精品999在线| 在线观看免费视频网站a站| 久久狼人影院| 国产蜜桃级精品一区二区三区| 另类亚洲欧美激情| 亚洲男人天堂网一区| 人人妻人人爽人人添夜夜欢视频| 亚洲久久久国产精品| 久久久国产精品麻豆| 亚洲性夜色夜夜综合| 久久精品亚洲av国产电影网| 天堂影院成人在线观看| 激情视频va一区二区三区| 国产不卡一卡二| 丝袜美腿诱惑在线| 中文字幕色久视频| 久久精品国产综合久久久| 国产高清视频在线播放一区| 琪琪午夜伦伦电影理论片6080| 成人手机av| 久久久久精品国产欧美久久久| 亚洲一区中文字幕在线| 天天影视国产精品| 电影成人av| 99久久综合精品五月天人人| 热99国产精品久久久久久7| 欧美日本亚洲视频在线播放| 日本a在线网址| 亚洲成a人片在线一区二区| 美女福利国产在线| 日本五十路高清| 国产精品日韩av在线免费观看 | 91大片在线观看| 国产精品1区2区在线观看.| 久久狼人影院| 电影成人av| 免费久久久久久久精品成人欧美视频| 亚洲av五月六月丁香网| 精品一品国产午夜福利视频| 日韩高清综合在线| 50天的宝宝边吃奶边哭怎么回事| 日韩免费av在线播放| 亚洲av成人av| 国产国语露脸激情在线看| 国产精品电影一区二区三区| 亚洲 国产 在线| 19禁男女啪啪无遮挡网站| 欧美色视频一区免费| 夜夜躁狠狠躁天天躁| 成人特级黄色片久久久久久久| 免费av毛片视频| 久久天躁狠狠躁夜夜2o2o| 欧美另类亚洲清纯唯美| 热99国产精品久久久久久7| 搡老岳熟女国产| 天堂中文最新版在线下载| 国产精品亚洲av一区麻豆| 9191精品国产免费久久| 天天影视国产精品| 嫁个100分男人电影在线观看| 正在播放国产对白刺激| 国产高清视频在线播放一区| 欧美日韩亚洲国产一区二区在线观看| 国产成+人综合+亚洲专区| 国产野战对白在线观看| 久久中文字幕一级| 精品久久蜜臀av无| 午夜激情av网站| 男女做爰动态图高潮gif福利片 | 国产极品粉嫩免费观看在线| 看片在线看免费视频| www.自偷自拍.com| 精品免费久久久久久久清纯| 美女扒开内裤让男人捅视频| 一边摸一边做爽爽视频免费| 久久久水蜜桃国产精品网| 69精品国产乱码久久久| 欧美最黄视频在线播放免费 | 在线观看午夜福利视频| 精品国内亚洲2022精品成人| 国产深夜福利视频在线观看| 久久99一区二区三区| 校园春色视频在线观看| 国产精品98久久久久久宅男小说| 亚洲av成人av| 女人被狂操c到高潮| 国产成+人综合+亚洲专区| 搡老岳熟女国产| 母亲3免费完整高清在线观看| 午夜成年电影在线免费观看| 不卡一级毛片| 午夜福利在线免费观看网站| 性少妇av在线| 老司机在亚洲福利影院| 99国产精品免费福利视频| 麻豆国产av国片精品| 日日摸夜夜添夜夜添小说| 韩国精品一区二区三区| 叶爱在线成人免费视频播放| 久久草成人影院| 美女国产高潮福利片在线看| 免费高清视频大片| 久久久久久亚洲精品国产蜜桃av| 老司机深夜福利视频在线观看| 精品一区二区三卡| 制服诱惑二区| 国产伦一二天堂av在线观看| a级片在线免费高清观看视频| 91国产中文字幕| 丁香六月欧美| 国产不卡一卡二| 日韩高清综合在线| 亚洲欧美一区二区三区久久| 国产精品国产av在线观看| 黄片大片在线免费观看| 9191精品国产免费久久| 国产精品国产av在线观看| 中文欧美无线码| 欧美最黄视频在线播放免费 | 成人黄色视频免费在线看| 中文字幕av电影在线播放| 国产av在哪里看| 中文字幕人妻丝袜制服| 久久热在线av| 国产男靠女视频免费网站| 亚洲av第一区精品v没综合| 亚洲avbb在线观看| 99久久久亚洲精品蜜臀av| 午夜福利在线免费观看网站| 国产片内射在线| 激情视频va一区二区三区| 国产一区二区激情短视频| 国产97色在线日韩免费| 不卡av一区二区三区| 午夜a级毛片| 麻豆av在线久日| 夜夜躁狠狠躁天天躁| 嫩草影视91久久| 国产亚洲精品久久久久5区| 午夜成年电影在线免费观看| 精品国产亚洲在线| 香蕉丝袜av| 免费久久久久久久精品成人欧美视频| 一区二区三区激情视频| 亚洲中文av在线| 精品一区二区三区四区五区乱码| 午夜免费鲁丝| 一级黄色大片毛片| 老司机靠b影院| 99久久人妻综合| 亚洲熟妇中文字幕五十中出 | 美女福利国产在线| 老司机在亚洲福利影院| 在线观看www视频免费| 在线视频色国产色| 人人妻人人爽人人添夜夜欢视频| 视频区图区小说| 久久久久久久久中文| 天天躁夜夜躁狠狠躁躁| 1024香蕉在线观看| cao死你这个sao货| 女人精品久久久久毛片| 淫妇啪啪啪对白视频| 国产高清视频在线播放一区| 精品福利永久在线观看| 久久 成人 亚洲| 亚洲国产精品999在线| 久久久久国内视频| 香蕉久久夜色| 国产精品免费一区二区三区在线| 亚洲免费av在线视频| 欧美日韩精品网址| 久久久国产精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人免费电影在线观看| 欧美日韩乱码在线| 91大片在线观看| 波多野结衣高清无吗| 夜夜爽天天搞| 免费观看精品视频网站| 天天影视国产精品| bbb黄色大片| 欧美激情久久久久久爽电影 | 一区二区三区精品91| 精品卡一卡二卡四卡免费| av在线播放免费不卡| 青草久久国产| 成年版毛片免费区| 国产精品爽爽va在线观看网站 | 麻豆av在线久日| 日日摸夜夜添夜夜添小说| 国产xxxxx性猛交| 亚洲五月婷婷丁香| 在线十欧美十亚洲十日本专区| 欧美激情 高清一区二区三区| 国产麻豆69| 亚洲av五月六月丁香网| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 黄色 视频免费看| 人成视频在线观看免费观看| 国产精品野战在线观看 | 最新美女视频免费是黄的| 婷婷精品国产亚洲av在线| 一级片免费观看大全| 可以在线观看毛片的网站| 午夜福利欧美成人| 18禁黄网站禁片午夜丰满| 一边摸一边抽搐一进一出视频| 亚洲第一青青草原| 两人在一起打扑克的视频| 亚洲人成网站在线播放欧美日韩| 满18在线观看网站| 看黄色毛片网站| 一区二区日韩欧美中文字幕| 十分钟在线观看高清视频www| 亚洲人成网站在线播放欧美日韩| 亚洲国产精品一区二区三区在线| 美女福利国产在线| 国产精品一区二区三区四区久久 | 在线看a的网站| 91老司机精品| 国产主播在线观看一区二区| 热99re8久久精品国产| 水蜜桃什么品种好| 天天添夜夜摸| 免费在线观看日本一区| 国产精品国产av在线观看| 亚洲成a人片在线一区二区| 亚洲中文字幕日韩| 一级片免费观看大全| 精品国产国语对白av| 久久天堂一区二区三区四区| 欧美丝袜亚洲另类 | 午夜福利,免费看| 国产97色在线日韩免费| 自线自在国产av| 一边摸一边抽搐一进一小说| 精品久久久精品久久久| 五月开心婷婷网| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 中国美女看黄片| 欧美乱色亚洲激情| 免费高清在线观看日韩| 人人妻人人添人人爽欧美一区卜| 国产av一区二区精品久久| 免费不卡黄色视频| 亚洲人成网站在线播放欧美日韩| 夜夜爽天天搞| 女性被躁到高潮视频| 久久久久久久午夜电影 | 成人亚洲精品一区在线观看| 亚洲avbb在线观看| 亚洲精品国产区一区二| 国产精品久久久久久人妻精品电影| 日本三级黄在线观看| 亚洲欧美日韩另类电影网站| av网站免费在线观看视频| 看黄色毛片网站| 欧美黄色片欧美黄色片| www.熟女人妻精品国产| 精品一区二区三卡| 可以免费在线观看a视频的电影网站| 欧美日韩亚洲高清精品| 大香蕉久久成人网| 在线观看66精品国产| 少妇裸体淫交视频免费看高清 | 精品乱码久久久久久99久播| 成人18禁在线播放| 男人舔女人下体高潮全视频| 在线av久久热| 涩涩av久久男人的天堂| 99久久国产精品久久久| 欧美日韩亚洲高清精品| 男人操女人黄网站| 精品国产亚洲在线| 一个人观看的视频www高清免费观看 | 久久精品国产99精品国产亚洲性色 | 亚洲国产精品合色在线| 香蕉丝袜av| 亚洲精品av麻豆狂野| 精品国产美女av久久久久小说| 欧美中文日本在线观看视频| 午夜免费成人在线视频| 亚洲人成电影免费在线| 热re99久久精品国产66热6| 交换朋友夫妻互换小说| 亚洲自偷自拍图片 自拍| 琪琪午夜伦伦电影理论片6080| 久久热在线av| 欧美激情 高清一区二区三区| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av在线| 69精品国产乱码久久久| 麻豆久久精品国产亚洲av | 三上悠亚av全集在线观看| 视频在线观看一区二区三区| 欧美激情久久久久久爽电影 | www国产在线视频色| www.自偷自拍.com| 三级毛片av免费| 欧美日韩黄片免| av超薄肉色丝袜交足视频| 欧美人与性动交α欧美软件| 欧美激情高清一区二区三区| av电影中文网址| 黄色视频不卡| 亚洲 欧美一区二区三区| 精品乱码久久久久久99久播| 久久热在线av| 美女高潮喷水抽搐中文字幕| 国产精品一区二区在线不卡| 亚洲视频免费观看视频| av天堂久久9| 久热这里只有精品99| 91麻豆av在线| 日韩 欧美 亚洲 中文字幕| e午夜精品久久久久久久| 999久久久国产精品视频| 两性夫妻黄色片| 亚洲精品成人av观看孕妇| 高清黄色对白视频在线免费看| 一本综合久久免费| 免费观看人在逋| 午夜91福利影院| 国产一区二区在线av高清观看| 久久热在线av| 18禁国产床啪视频网站| 久久热在线av| 亚洲成av片中文字幕在线观看| 国产人伦9x9x在线观看| 亚洲激情在线av| 色尼玛亚洲综合影院| 国产又色又爽无遮挡免费看| 色综合欧美亚洲国产小说| 日韩成人在线观看一区二区三区| 欧美丝袜亚洲另类 | 丝袜美腿诱惑在线| 中文字幕精品免费在线观看视频| 黑丝袜美女国产一区| 中文字幕高清在线视频| 久久久精品欧美日韩精品| 美女福利国产在线| 日韩视频一区二区在线观看| 韩国av一区二区三区四区| 91成年电影在线观看| 看片在线看免费视频| 啦啦啦 在线观看视频| 日韩大尺度精品在线看网址 | 成人黄色视频免费在线看| 精品一区二区三区视频在线观看免费 | svipshipincom国产片| 别揉我奶头~嗯~啊~动态视频| 亚洲三区欧美一区| 欧美精品一区二区免费开放| 夜夜看夜夜爽夜夜摸 | 757午夜福利合集在线观看| 久久人妻熟女aⅴ| 好看av亚洲va欧美ⅴa在| 91成年电影在线观看| 美女福利国产在线| 极品教师在线免费播放| 日韩免费高清中文字幕av| 极品教师在线免费播放| 老司机深夜福利视频在线观看| 免费观看人在逋| 亚洲专区字幕在线| 欧美成人性av电影在线观看| 又黄又爽又免费观看的视频| 国产精品 国内视频| av免费在线观看网站| 18禁美女被吸乳视频| 成人永久免费在线观看视频| 久久99一区二区三区| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 一二三四在线观看免费中文在| 黑丝袜美女国产一区| 在线观看舔阴道视频| 黄色视频,在线免费观看| 精品久久久精品久久久| 亚洲国产精品一区二区三区在线| 欧美精品亚洲一区二区| 免费看十八禁软件| 身体一侧抽搐| 国产欧美日韩综合在线一区二区| 韩国精品一区二区三区| 女人被狂操c到高潮| 丁香欧美五月| 亚洲国产中文字幕在线视频| 亚洲视频免费观看视频| 黄色怎么调成土黄色| 亚洲熟妇熟女久久| 91在线观看av| 黄频高清免费视频| 激情在线观看视频在线高清| 国产精品久久视频播放| 国产无遮挡羞羞视频在线观看| 麻豆国产av国片精品| 午夜精品国产一区二区电影| 两人在一起打扑克的视频| 精品久久久久久,| 欧美中文日本在线观看视频| 国产高清激情床上av| 亚洲精品国产色婷婷电影| 亚洲第一青青草原| 中文字幕av电影在线播放| 久久久久久大精品| 不卡av一区二区三区| 日韩av在线大香蕉| 在线观看免费视频日本深夜| 如日韩欧美国产精品一区二区三区| 亚洲一区二区三区色噜噜 | 亚洲熟女毛片儿| 国产精品成人在线| 国产欧美日韩一区二区三区在线| 老司机午夜十八禁免费视频| 国产成人系列免费观看| 亚洲男人的天堂狠狠| 另类亚洲欧美激情| 脱女人内裤的视频| 极品教师在线免费播放| 激情在线观看视频在线高清| 亚洲五月婷婷丁香| 亚洲专区中文字幕在线| 国产成人欧美在线观看| 日韩欧美一区二区三区在线观看| 夜夜躁狠狠躁天天躁| 琪琪午夜伦伦电影理论片6080| 咕卡用的链子| 久久久久久久久久久久大奶| 超碰97精品在线观看| 99香蕉大伊视频| 韩国精品一区二区三区| 色综合欧美亚洲国产小说| av视频免费观看在线观看| 搡老乐熟女国产| 亚洲熟妇中文字幕五十中出 | 免费搜索国产男女视频| 亚洲精品在线美女| 精品国产一区二区久久| 精品国产一区二区三区四区第35| av天堂久久9| 在线天堂中文资源库| 9191精品国产免费久久| 国产精品久久久久久人妻精品电影| 一边摸一边做爽爽视频免费| 一级作爱视频免费观看| 亚洲欧美日韩无卡精品| 国产人伦9x9x在线观看| 国产亚洲欧美在线一区二区| 成人永久免费在线观看视频| 1024视频免费在线观看| 国产精品综合久久久久久久免费 | cao死你这个sao货| 黄色成人免费大全| 麻豆国产av国片精品| 久久久久亚洲av毛片大全| 首页视频小说图片口味搜索| 黄色女人牲交| 久久青草综合色| 极品教师在线免费播放| 九色亚洲精品在线播放| tocl精华| 亚洲成人精品中文字幕电影 | 午夜影院日韩av| 精品久久久久久成人av| 亚洲精品美女久久久久99蜜臀| 在线观看免费高清a一片| 日本免费a在线| 黄色成人免费大全| 欧美日韩福利视频一区二区| 日韩成人在线观看一区二区三区| av天堂久久9| 琪琪午夜伦伦电影理论片6080| 欧美人与性动交α欧美软件| 亚洲人成电影观看| 亚洲专区中文字幕在线| 精品日产1卡2卡| 午夜免费观看网址| 色综合婷婷激情| 人成视频在线观看免费观看| 国产成人欧美| 亚洲成国产人片在线观看| 少妇被粗大的猛进出69影院| 国产精品电影一区二区三区| 日本a在线网址| 日韩中文字幕欧美一区二区| 亚洲人成伊人成综合网2020| 久热这里只有精品99| 脱女人内裤的视频| 亚洲欧美一区二区三区黑人| 久久久久久免费高清国产稀缺| 性少妇av在线| 久久精品国产亚洲av香蕉五月| 黄色怎么调成土黄色| 丝袜美腿诱惑在线| 亚洲熟女毛片儿| x7x7x7水蜜桃| 在线观看免费视频网站a站| 人人澡人人妻人| 国产精品自产拍在线观看55亚洲| 亚洲色图av天堂| 麻豆国产av国片精品| 亚洲自拍偷在线| 久久久久精品国产欧美久久久| 五月开心婷婷网| a级毛片黄视频| 真人一进一出gif抽搐免费| 亚洲自拍偷在线| 91大片在线观看| 天天添夜夜摸| 国产视频一区二区在线看| 一进一出抽搐动态| 一区二区日韩欧美中文字幕| 久久 成人 亚洲| 人成视频在线观看免费观看| ponron亚洲| 亚洲avbb在线观看| 国产三级黄色录像| 97人妻天天添夜夜摸| 久9热在线精品视频| 午夜影院日韩av| 国产日韩一区二区三区精品不卡| 亚洲九九香蕉| 久久久久久免费高清国产稀缺| 亚洲一区高清亚洲精品| 88av欧美| 国产精品乱码一区二三区的特点 | 午夜免费鲁丝| 亚洲欧美一区二区三区黑人| 国产成人欧美| 欧美日韩亚洲国产一区二区在线观看| 黄色视频不卡| 水蜜桃什么品种好| 久久久精品国产亚洲av高清涩受| 另类亚洲欧美激情| 国产熟女xx| 午夜免费激情av| 亚洲美女黄片视频| 精品一区二区三区av网在线观看| 精品免费久久久久久久清纯| 黄色 视频免费看| 成人三级做爰电影| 女生性感内裤真人,穿戴方法视频| 日韩精品中文字幕看吧| 操出白浆在线播放| 男女做爰动态图高潮gif福利片 | 精品福利观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产三级在线视频| 国产精品一区二区免费欧美| 中文字幕另类日韩欧美亚洲嫩草| 不卡一级毛片| 国产亚洲欧美精品永久| 成人免费观看视频高清| 一级a爱视频在线免费观看| 无人区码免费观看不卡| 亚洲视频免费观看视频| 国产成+人综合+亚洲专区| 中文欧美无线码| 日本欧美视频一区| 亚洲精品av麻豆狂野| av超薄肉色丝袜交足视频| 乱人伦中国视频| 国产高清国产精品国产三级| 丝袜美腿诱惑在线| 女人被狂操c到高潮| 亚洲熟妇中文字幕五十中出 | 欧美色视频一区免费| 精品国产一区二区三区四区第35| 国产1区2区3区精品| 免费人成视频x8x8入口观看| 波多野结衣高清无吗| 最新美女视频免费是黄的| 久久人妻熟女aⅴ| 午夜视频精品福利| 在线播放国产精品三级| 香蕉丝袜av| 亚洲国产欧美日韩在线播放| 欧洲精品卡2卡3卡4卡5卡区| 成在线人永久免费视频| 一区二区三区激情视频| 精品日产1卡2卡| videosex国产| 亚洲精品成人av观看孕妇| 久久久国产成人精品二区 | 国产精品野战在线观看 | 99久久人妻综合| xxx96com| 人成视频在线观看免费观看| 久久久久亚洲av毛片大全| 国产精品一区二区免费欧美| 大型av网站在线播放| 99久久人妻综合| 精品久久久久久久久久免费视频 | 搡老熟女国产l中国老女人| 国产精华一区二区三区| 国产真人三级小视频在线观看| 成人精品一区二区免费| 一区二区三区精品91| 日韩 欧美 亚洲 中文字幕| 亚洲 欧美一区二区三区| 18美女黄网站色大片免费观看| 国产男靠女视频免费网站| 夜夜躁狠狠躁天天躁| a级毛片在线看网站| 亚洲自拍偷在线| 免费一级毛片在线播放高清视频 | 9色porny在线观看| 欧美在线一区亚洲| 女人高潮潮喷娇喘18禁视频| 50天的宝宝边吃奶边哭怎么回事| 电影成人av|