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

    樟樹(shù)5種化學(xué)類(lèi)型葉片轉(zhuǎn)錄組分析

    2014-11-17 06:42:20江香梅伍艷芳肖復(fù)明熊振宇徐海寧
    遺傳 2014年1期
    關(guān)鍵詞:芳樟醇樟樹(shù)測(cè)序

    江香梅,伍艷芳,肖復(fù)明,熊振宇,徐海寧

    江西省林業(yè)科學(xué)院,國(guó)家林業(yè)局樟樹(shù)工程技術(shù)研究中心,南昌 330032

    樟樹(shù)(Cinnamomum camphora(L.) Presl)是我國(guó)特有國(guó)家Ⅱ級(jí)保護(hù)樹(shù)種,是集材用、藥用、香料、油用、化工、觀賞、生態(tài)環(huán)境和生態(tài)文化建設(shè)等于一體的多用途樹(shù)種,可作為樟科的代表種,具有極重要的開(kāi)發(fā)利用價(jià)值[1]。樟樹(shù)的化學(xué)成分復(fù)雜多樣,其根、莖、葉、花、果中均富含精油,從其精油中已鑒定出 150余種化合物,且新化合物還在不斷被鑒定出來(lái)。我國(guó)是世界上生產(chǎn)樟油最多的國(guó)家,其產(chǎn)量占世界的 80%,產(chǎn)品質(zhì)量享譽(yù)全球。樟樹(shù)按葉精油主要化學(xué)成分種類(lèi)及含量的不同,可分為芳樟(精油中主成分為芳樟醇,下同)、腦樟(樟腦)、油樟(1,8-桉葉油素)、異樟(異-橙花叔醇)和龍腦樟(右旋龍腦)5種主要化學(xué)類(lèi)型[2]。迄今為止,對(duì)樟樹(shù)精油化學(xué)成分分析與利用的研究較多[3~5],而對(duì)化學(xué)成分代謝途徑的研究則較少[6,7]。目前,以樟樹(shù)為代表的樟科植物尚未完成全基因組測(cè)序,EST文庫(kù)尚未建立,GenBank上樟樹(shù)EST序列也較少,相關(guān)基因的功能注釋亦不完善,給樟樹(shù)萜類(lèi)化合物次生代謝途徑的研究帶來(lái)極大困難,亟需更多的基因組或轉(zhuǎn)錄組信息來(lái)解決這一問(wèn)題。

    新一代測(cè)序技術(shù)(Next-generation sequencing technologies,NGSTs)對(duì)分子生物學(xué)的發(fā)展起到了巨大的推動(dòng)作用。其中用于轉(zhuǎn)錄組測(cè)序和數(shù)字基因表達(dá)譜(Digital gene expression profiling,DGE)研究的RNA-seq技術(shù)不僅能夠廣泛應(yīng)用于有參考基因組序列的物種研究,如水稻(Oryza sativa)[8]、葡萄(Vitis vinifera)[9]、黃瓜(Cucumis sativus)[10]等,也能應(yīng)用于無(wú)參考基因組序列的物種,如青蒿(Artemisia annua)[11]、人參(Panax ginseng)[12]、紅豆杉(Taxus mairei)[13]、羅漢果(Siraitia grosvenorii)[14]等,應(yīng)用較為廣泛[15]?;谵D(zhuǎn)錄組測(cè)序得到的基因功能注釋、蛋白編碼區(qū)注釋等大量的信息,可以從基因的表達(dá)水平[9,13]、SNP鑒定[11,16]、SSR分子標(biāo)記篩選[17,18]、候選基因的挖掘[19,20]、融合轉(zhuǎn)錄本的表達(dá)[8,21]、可變剪接[8,22]等方面展開(kāi)相關(guān)研究,并建立相關(guān)轉(zhuǎn)錄組數(shù)據(jù)庫(kù)[17],為進(jìn)一步研究提供重要基礎(chǔ)。本研究采用 Illumina HiSeq? 2000新一代高通量測(cè)序技術(shù)對(duì)樟樹(shù)5種化學(xué)類(lèi)型葉片轉(zhuǎn)錄組進(jìn)行測(cè)序,對(duì)測(cè)序得到的大量Unigene進(jìn)行GO、COG和KEGG分類(lèi)統(tǒng)計(jì),給出功能注釋和Pathway注釋,并預(yù)測(cè)Unigene蛋白CDS。這些注釋信息的完成將為樟樹(shù)精油主要成分合成、功能基因及相關(guān)候選基因的發(fā)掘提供基礎(chǔ)數(shù)據(jù),同時(shí)也為進(jìn)一步克隆功能基因全長(zhǎng)、研究基因功能奠定了基礎(chǔ),對(duì)推動(dòng)我國(guó)樟科植物分子生物學(xué)研究將起到極大促進(jìn)作用。

    1 材料和方法

    1.1 材 料

    試驗(yàn)材料采自江西省林業(yè)科學(xué)院內(nèi)年齡約 30年的樟樹(shù)成年植株。于4月底統(tǒng)一采集腦樟、芳樟、油樟、異樟、龍腦樟 5種化學(xué)類(lèi)型幼嫩葉片,液氮速凍后-70℃低溫保存,用于提取RNA,測(cè)序。

    1.2 方法

    1.2.1 RNA提取、純化和文庫(kù)構(gòu)建

    采用通用植物總RNA提取試劑盒RNeasy Plant Mini Kit提取樟樹(shù)5種化學(xué)類(lèi)型葉片總RNA,瓊脂糖凝膠電泳檢測(cè)RNA完整性,Agilent 2100 Bioanalyzer檢測(cè)總 RNA 濃度。用 Oligotex? mRNA kit(TaKaRa)分離純化 mRNA。得到的 mRNA采用fragmentation buffer打斷成小片段,經(jīng)過(guò)PCR擴(kuò)增,建立小片段測(cè)序文庫(kù)。文庫(kù)測(cè)序采用 Illumina Hi-Seq? 2000完成。

    1.2.2 序列拼接

    測(cè)序后得到的Raw reads,去除含有帶接頭的、重復(fù)的(N>5%)、測(cè)序質(zhì)量很低的reads(質(zhì)量值Q≤10的堿基數(shù)占整個(gè) reads的 20%以上),獲得 Clean reads。采用短 reads組裝軟件 Trinity[23]做轉(zhuǎn)錄組從頭組裝。先將具有一定長(zhǎng)度overlap的reads連成更長(zhǎng)的片段,得到Contig組裝片段,再將reads比對(duì)回Contig,通過(guò) paired-end reads來(lái)確定來(lái)自同一轉(zhuǎn)錄本的不同Contig以及這些Contig之間的距離,將這些 Contig連在一起,得到兩端不能再延長(zhǎng)的序列,即為Unigene。

    1.2.3 功能注釋和CDS預(yù)測(cè)

    功能注釋信息給出Unigene的蛋白功能注釋、COG功能注釋。先通過(guò)BlastX將Unigene序列比對(duì)到蛋白數(shù)據(jù)庫(kù) Nr(NCBI非冗余蛋白庫(kù))、SwissProt(去冗余的蛋白數(shù)據(jù)庫(kù))、KEGG(系統(tǒng)分析基因產(chǎn)物在細(xì)胞中的代謝途徑以及這些基因產(chǎn)物功能的數(shù)據(jù)庫(kù))和 COG(對(duì)基因產(chǎn)物進(jìn)行直系同源分類(lèi)的數(shù)據(jù)庫(kù))(E值<1e-5),再通過(guò) BlastN將Unigene比對(duì)到核酸數(shù)據(jù)庫(kù) Nt(NCBI非冗余核酸數(shù)據(jù)庫(kù))(E值<1e-5),得到跟給定 Unigene具有最高序列相似性的蛋白,從而得到該 Unigene的蛋白功能注釋信息。采用 Blast2GO軟件[24],根據(jù)Nr注釋信息得到 GO注釋信息; 采用 WEGO軟件[25]對(duì) All-Unigene(按分子功能、細(xì)胞組分、生物學(xué)過(guò)程)進(jìn)行GO功能分類(lèi)統(tǒng)計(jì),從宏觀上認(rèn)識(shí)樟樹(shù)的基因功能分布特征。具體測(cè)序數(shù)據(jù)分析參照林萍等[26]的轉(zhuǎn)錄組注釋分類(lèi)方法。將 Unigene序列按 Nr、SwissProt、KEGG和 COG的優(yōu)先級(jí)順序做BlastX比對(duì)(E值<1e-5)。取比對(duì)結(jié)果中rank最高的蛋白確定該 Unigene的編碼區(qū)序列,根據(jù)標(biāo)準(zhǔn)密碼子表將編碼區(qū)序列翻譯成氨基酸序列,從而得到該Unigene編碼區(qū)的核酸序列(序列方向5′→3′)和氨基酸序列。與 4個(gè)數(shù)據(jù)庫(kù)皆比對(duì)不上的Unigene用ESTscan[27]軟件預(yù)測(cè)其編碼區(qū),得到其編碼區(qū)的核酸序列(序列方向5'→3')和氨基酸序列。

    2 結(jié)果與分析

    2.1 樟樹(shù)5種化學(xué)類(lèi)型葉片轉(zhuǎn)錄組測(cè)序產(chǎn)量

    樟樹(shù) 5種化學(xué)類(lèi)型葉片轉(zhuǎn)錄組測(cè)序后得到的Raw reads及去除雜質(zhì)之后的Clean reads結(jié)果列于表1。由表1可以看出,質(zhì)量參數(shù)Q20最低為95.08%(異樟),最高達(dá)到 98.51%(腦樟); 過(guò)濾后不確定的堿基比例N值≤0.01%; 5種化學(xué)類(lèi)型的GC含量在46.55%~49.29%之間。

    2.2 樟樹(shù)5種化學(xué)類(lèi)型葉片轉(zhuǎn)錄組組裝質(zhì)量

    2.2.1 數(shù)據(jù)統(tǒng)計(jì)

    采用短reads組裝軟件Trinity從頭組裝的樟樹(shù)5種化學(xué)類(lèi)型序列重疊群Contig和Unigene數(shù)目列于表 2。由表 2可知,經(jīng)過(guò)拼接最終獲得長(zhǎng)度在 100~2 000 bp之間的Unigene 156 278個(gè),總長(zhǎng)87.09 Mb,序列平均長(zhǎng)度584 bp,N50為1 023 bp。

    表1 測(cè)序產(chǎn)量統(tǒng)計(jì)表

    表2 組裝質(zhì)量統(tǒng)計(jì)表

    2.2.2 Contig和Unigene的長(zhǎng)度分布

    Contig和Unigene的長(zhǎng)度分布分別見(jiàn)圖1和圖2。樟樹(shù)All-Unigene長(zhǎng)度在100~500 bp的有114 115個(gè),占73.02%; 500~1 000的有18 316個(gè),占11.72%;1 000~1 500的有9 144個(gè),占5.85%; 1 500~2 000的有6 365個(gè),占4.07%; ≥2 000的有8 338個(gè),占5.34%。

    2.3 Unigene的功能注釋

    2.3.1 注釋結(jié)果統(tǒng)計(jì)

    分別將 Unigene注釋到 Nr、Nt、SwissProt、KEGG、COG、GO庫(kù),并分別對(duì)注釋到每個(gè)庫(kù)以及所有注釋上的Unigene數(shù)目進(jìn)行統(tǒng)計(jì),結(jié)果見(jiàn)表3。通過(guò)Blast搜索比對(duì),共有55 955條Unigene獲得了基因注釋,占All-Unigene的35.80%; 有100 323條Unigene(64.20%)未被注釋。Nr數(shù)據(jù)庫(kù)比對(duì)注釋的信息最多,注釋了55 257條Unigene,COG注釋的信息最少,僅21 806條Unigene得到了注釋。

    2.3.2 Unigene的COG分類(lèi)

    為了進(jìn)一步評(píng)價(jià)轉(zhuǎn)錄組文庫(kù)的完整性和注釋的有效性,對(duì)Unigene進(jìn)行了COG分類(lèi)(圖3)。在COG 25個(gè)類(lèi)別中,“一般功能基因”為最大 5的組(7 466,34.24%),其次是“轉(zhuǎn)錄”(4 955,22.72%)和“復(fù)制、重組和修復(fù)”(3 803,17.44%)。3個(gè)最小的組別分別為“核苷酸結(jié)構(gòu)”(2,0.01%)、“細(xì)胞外結(jié)構(gòu)”(11,0.05%)和“RNA 加工與修飾”(358,1.64%)。此外,參與次生代謝生物合成、運(yùn)輸和分解的有 1 296個(gè)Unigene,占5.94%。

    2.4 Unigene的GO分類(lèi)

    在已經(jīng)得到的 Nr注釋信息基礎(chǔ)上,采用Blast2GO獲得樟樹(shù)Unigene的GO分類(lèi)信息,共有24 717條Unigene得到GO注釋。在GO分類(lèi)體系中,生物學(xué)過(guò)程、細(xì)胞組分和分子功能 3個(gè)大的類(lèi)別被劃分為詳細(xì)的 44個(gè)小的類(lèi)別,其中“細(xì)胞”(16 006,14.52%)、“細(xì)胞要素”(14 515,13.17%)和“細(xì)胞器”(11 547,10.48%)3個(gè)類(lèi)群占了主要部分,隨后是“催化活性”(11 074,10.05%)、“結(jié)合活性”(10 642,9.65%)和“代謝過(guò)程”(9 930,9.01%)3 個(gè)類(lèi)群,而“細(xì)胞殺傷”(1,0.001%)、“律動(dòng)過(guò)程”(2,0.002%)和“氮素利用”(3,0.003%)僅有非常少的基因歸入,這一分類(lèi)結(jié)果顯示了樟樹(shù)葉基因表達(dá)譜的總體情況(圖4)。

    2.5 KEGG pathways分析

    圖1 Contig的長(zhǎng)度分布統(tǒng)計(jì)圖

    利用KEGG數(shù)據(jù)庫(kù)進(jìn)行了功能分類(lèi)和Pathway注釋[28]。首先,將156 278條All-Unigene采用BlastX比對(duì)到KEGG數(shù)據(jù)庫(kù),結(jié)果有48 875條序列能夠比對(duì)上(E值<1e-5),共有32 885條Unigene能夠注釋到217個(gè)KEGG標(biāo)準(zhǔn)Pathway,另外123 393條序列則沒(méi)有相應(yīng)的生物學(xué) Pathway注釋(Condition:expect≤1e-5; rank≤5),選取注釋基因比例大于1%(占所有注釋基因)的Pathway列于表4。由表4可知,注釋到“代謝途徑”中的Unigene最多,有7 808條,占23.74%; 有3 350條基因(10.19%)注釋到次生代謝生物合成途徑,其中,參與單萜[PATH:ko00902](66,0.2%)、二萜[PATH:ko00904](113,0.34%)、倍半萜[PATH:ko00909](34,0.1%)和萜類(lèi)骨架合成[PATH:ko00900](211,0.24%)的Unigene共424個(gè)。參與不飽和脂肪酸生物合成的 Unigene有 191條,占0.58%。這些有代表性的注釋為研究樟樹(shù)特殊生物學(xué)進(jìn)程、功能和代謝提供了重要依據(jù)。

    圖2 Unigene的長(zhǎng)度分布統(tǒng)計(jì)圖

    2.6 Unigene的編碼蛋白框(CDS)預(yù)測(cè)

    將 Unigene序列按 Nr、SwissProt、KEGG和COG數(shù)據(jù)庫(kù)的優(yōu)先級(jí)順序分別做BlastX比對(duì)(E值<1e-5),確定該 Unigene的編碼區(qū)序列,然后根據(jù)標(biāo)準(zhǔn)密碼子表將編碼區(qū)序列翻譯成氨基酸序列,從而得到該Unigene編碼區(qū)的核酸序列(序列方向 5′→3′)和氨基酸序列。最后,跟以上 4個(gè)數(shù)據(jù)庫(kù)皆比對(duì)不上的Unigene用ESTscan軟件預(yù)測(cè)其編碼區(qū),得到其編碼區(qū)核酸序列(序列方向?yàn)?5′→3′)和氨基酸序列。圖 5分別顯示All-Unigene與數(shù)據(jù)庫(kù)Blast的CDS核酸和氨基酸序列分布,ESTscan預(yù)測(cè)編碼區(qū)的核酸和氨基酸序列的長(zhǎng)度分布。

    表3 Unigene注釋結(jié)果統(tǒng)計(jì)表

    2.7 樟樹(shù)芳樟醇合酶基因在不同化學(xué)類(lèi)型中的表達(dá)模式分析

    樟樹(shù)5種化學(xué)類(lèi)型精油的主要成分為萜類(lèi)物質(zhì),且多數(shù)為單萜和倍半萜類(lèi)。在 Pathway單萜合成的代謝通路中,找到9條Unigene可能編碼樟樹(shù)芳樟醇合成途徑的關(guān)鍵酶芳樟醇合酶基因,它們?cè)?5種化學(xué)類(lèi)型中的表達(dá)水平(Reads Per kb per Million reads,RPKM值)見(jiàn)表5。表5顯示,芳樟醇合酶基因在芳樟中優(yōu)勢(shì)表達(dá),而在油樟中表達(dá)水平較低。這些Unigene的注釋信息將為進(jìn)一步克隆功能基因的全長(zhǎng)、研究其功能提供基礎(chǔ)數(shù)據(jù),也為樟樹(shù)精油的代謝調(diào)控研究奠定了基礎(chǔ)。

    圖3 Unigene的COG分類(lèi)圖

    圖4 Unigene的GO分類(lèi)圖

    表4 KEGG pathway注釋結(jié)果統(tǒng)計(jì)表

    3 討 論

    轉(zhuǎn)錄組是功能基因組研究的一個(gè)重要功能指標(biāo)[29]。本研究首次采用 Illumina高通量測(cè)序技術(shù)對(duì)樟樹(shù) 5種化學(xué)類(lèi)型葉片轉(zhuǎn)錄組進(jìn)行測(cè)序,共拼接得到156 278條Unigene。其中,未獲得鑒定的新Unigene 100 323條,占總數(shù)的64.20%,為進(jìn)一步挖掘并鑒定新的功能基因提供了豐富的數(shù)據(jù)信息。本研究測(cè)序所得到的Unigene序列平均長(zhǎng)度584 bp,N50為1 023 bp,完全滿足轉(zhuǎn)錄組測(cè)序的要求。傳統(tǒng)的基因表達(dá)序列標(biāo)簽(Expressed sequence tags,ESTs)技術(shù)被認(rèn)為是一種研究轉(zhuǎn)錄組的有效方法,廣泛應(yīng)用于新基因發(fā)現(xiàn)、基因表達(dá)分析和蛋白質(zhì)組學(xué)?;蛐酒?Gene chip)技術(shù)也廣泛用于大量核酸分子的檢測(cè)分析,為研究不同層次多基因協(xié)同作用提供手段。與EST 技術(shù)及基因芯片技術(shù)相比,基于 Illumina HiSeqTM2000 高通量測(cè)序?qū)D(zhuǎn)錄組進(jìn)行比較和分析,所需的RNA量較少,背景噪音比基因芯片低,繪制轉(zhuǎn)錄組遺傳圖譜所需的費(fèi)用更低。該平臺(tái)可同時(shí)用于有(無(wú))參考基因組的轉(zhuǎn)錄組測(cè)序,通過(guò)短序列組裝軟件Trinity獲得的信息量完全可滿足研究需求。

    木本植物由于多數(shù)為異交物種,雜合性較強(qiáng),基因組相對(duì)較大且較為復(fù)雜,從而導(dǎo)致遺傳背景研究相對(duì)滯后。而對(duì)于遺傳背景不清晰的木本植物,可先采用高通量測(cè)序技術(shù)進(jìn)行轉(zhuǎn)錄組測(cè)序,以獲得的大量 Unigene信息構(gòu)建遺傳和物理圖譜,為待測(cè)序的物種提供遺傳背景信息[31]。林萍等[27]曾采用Illumina的Solexa技術(shù)對(duì)普通油茶種子4個(gè)發(fā)育時(shí)期的轉(zhuǎn)錄組進(jìn)行測(cè)序,經(jīng)組裝分析獲得 80 310 條Unigene,其中確定編碼蛋白功能的有 21 789 條,占All-Unigene 的27.13%。本研究中,Unigene注釋結(jié)果顯示共有 55 955條 Unigene獲得了基因注釋,占 All-Unigene的 35.80%,略高于油茶轉(zhuǎn)錄組的注釋,但是相對(duì)于草本植物來(lái)說(shuō)注釋結(jié)果偏低,表明現(xiàn)有數(shù)據(jù)庫(kù)中,樟科植物基因注釋信息量極為缺乏。樟樹(shù)作為樟科植物的代表樹(shù)種,其轉(zhuǎn)錄組測(cè)序及相關(guān)注釋具有極為重要的意義,將為樟科其他物種提供注釋信息參考。高通量測(cè)序技術(shù)可以大規(guī)模的對(duì)生物體組織樣本進(jìn)行測(cè)序分析,有利于建立特定時(shí)空條件下的物質(zhì)代謝途徑[29]。Sun等[19]對(duì)西洋參(Panax quinquefoliusL.)轉(zhuǎn)錄組進(jìn)行測(cè)序,KEGG分析確定有 4 097條序列被定位到特定的代謝途徑中,并且初步確定了從 acetyl CoA開(kāi)始經(jīng)過(guò)類(lèi)異戊二烯途徑的所有參與人參皂苷骨架合成的酶。此外,東北紅豆杉(Taxus cuspidataSieb.et Zucc.)中的紫杉醇[32]、蛇足石杉(Huperzia serrata(Thunb.ex Murray)Trev.)和龍骨馬尾杉(Phlegmariurus carinatus(Desv.ex Poir.) Ching.)中的石松生物堿[33]等物質(zhì)也通過(guò)高通量測(cè)序進(jìn)一步確定其次生代謝途徑。樟樹(shù)Unigene的GO分類(lèi)結(jié)果顯示了樟樹(shù)5種化學(xué)類(lèi)型葉基因表達(dá)譜的總體情況,共有24 717條Unigene得到GO注釋,在 44個(gè)小類(lèi)別中,歸入“細(xì)胞”和“細(xì)胞器”類(lèi)別的分別為 16 006(14.52%)和 11 547(10.48%)條Unigene,說(shuō)明有較多的基因與細(xì)胞和細(xì)胞器中的生物代謝相關(guān)。KEGG pathways分析結(jié)果表明,共有3 350(10.19%)條基因注釋到次生代謝生物合成途徑,其中,參與萜類(lèi)代謝的Unigene共424條。在此基礎(chǔ)上進(jìn)行芳樟醇合成途徑分析,結(jié)果顯示芳樟醇合酶基因在芳樟中優(yōu)勢(shì)表達(dá),而在油樟中表達(dá)水平較低,這一結(jié)果與樟樹(shù) 5種化學(xué)型中芳樟醇的含量高低相一致[2]。這些有代表性的注釋為研究樟樹(shù)次生代謝過(guò)程、脂肪酸合成途徑及其他特殊生物學(xué)進(jìn)程提供了重要依據(jù)。

    圖5 CDS的長(zhǎng)度分布統(tǒng)計(jì)圖

    表5 編碼芳樟醇合酶的Unigene注釋結(jié)果統(tǒng)計(jì)表

    [1]戴寶合.野生植物資源學(xué).北京:農(nóng)業(yè)出版社,1993.

    [2]石皖陽(yáng),何偉,文光裕,郭德選,龍光遠(yuǎn),劉銀茍.樟精油成分和類(lèi)型劃分.植物學(xué)報(bào),1989,31(3):209–214.

    [3]彭東輝.樟樹(shù)優(yōu)良單株選擇與組培研究[學(xué)位論文].福建農(nóng)林大學(xué),2004.

    [4]Lee HJ,Hyuna EA,Yoon WJ,Kim BH,Rhee MH,Kang HK,Cho JY,Yoo ES.In vitroanti-inflammatory and anti-oxidative effects ofCinnamomum camphoraextracts.J Ethnopharmacol,2006,103(2):208–216.

    [5]Liu CH,Mishra AK,Tan RX,Tang C,Yang H,Shen YF.Repellent and insecticidal activities of essential oils fromArtemisia princepsandCinnamomum camphoraand their effect on seed germination of wheat and broad bean.Bioresour Technol,2006,97(15):1969–1973.

    [6]Yang T,Li J,Wang HX,Zeng Y.A geraniol-synthase gene fromCinnamomum tenuipilum.Phytochemistry,2005,66(3):285–293.

    [7]陳美蘭.藥用植物樟化學(xué)型形成機(jī)理的基礎(chǔ)研究[學(xué)位論文].中國(guó)中醫(yī)科學(xué)院,2007.

    [8]Zhang GJ,Guo GW,Hu XD,Zhang Y,Li QY,Li RQ,Zhang RH,Lu ZK,He ZQ,Fang XD,Chen L,Tian W,Tao Y,Kristiansen K,Zhang XQ,Li SG,Yang HM,Wang J.Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome.Genome Res,2010,20(5):646–654.

    [9]Wu J,Zhang YL,Zhang HQ,Huang H,Folta KM,Lu J.Whole genome wide expression profiles ofVitis amurensisgrape responding to downy mildew by using Solexa sequencing technology.BMC Plant Biol,2010,10:234.

    [10]Guo SG,Zheng Y,Joung JG,Liu SQ,Zhang ZH,Crasta OR,Sobral BW,Xu Y,Huang SW,Fei ZJ.Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types.BMC Genomics,2010,11:384.

    [11]Graham IA,Besser K,Blumer S,Branigan CA,Czechowski T,Elias L,Guterman I,Harvey D,Isaac PG,Khan AM,Larson TR,Li Y,Pawson T,Penfield T,Rae AM,Rathbone DA,Reid S,Ross J,Swallwood MF,Segura V,Townsend T,Vyas D,Winzer T,Bowles D.The genetic map ofArtemisia annuaL.identifies loci affecting yield of the antimalarial drug artemisinin.Science,2010,327(5963):328–331.

    [12]Chen S,Luo H,Li Y,Sun Y,Wu Q,Niu Y,Song J,Lv A,Zhu Y,Sun C,Steinmetz A,Qian Z.454 EST analysis detects genes putatively involved in ginsenoside biosynthesis inPanax ginseng.Plant Cell Rep,2011,30(9):1593–1601.

    [13]Hao DC,Ge GB,Xiao PG,Zhang YY,Yang L.The first insight into the tissue specificTaxustranscriptome via Illumina second generation sequencing.PLoS ONE,2011,6(6):e21220.

    [14]Tang Q,Ma XJ,Mo CM,Wilson IW,Song C,Zhao H,Yang YF,Fu W,Qiu DY.An efficient approach to findingSiraitia grosvenoriitriterpene biosynthetic genes by RNA-seq and digital gene expression analysis.BMC Genomics,2011,12(1):343.

    [15]Wilhelm BT,Landry JR.RNA-Seq quantitative measurement of expression through massively parallel RNA Sequencing.Methods,2009,48(3):249–257.

    [16]Zenoni S,Ferrarini A,Giacomelli E,Xumerle L,Fasoli M,Malerba G,Bellin D,Pezzotti M,Delledonne M.Characterization of transcriptional complexity during berry development inVitis viniferausing RNA-Seq.Plant Physiol,2010,152(4):1787–1795.

    [17]Wang ZY,Fang BP,Chen JY,Zhang XJ,Luo ZX,Huang LF,Chen XL,Li YJ.De novoassembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweetpotato(Ipomoea batatas).BMC Genomics,2010,11(1):726.

    [18]李瀅,孫超,羅紅梅,李西文,牛云云,陳士林.基于高通量測(cè)序454 GS FLX的丹參轉(zhuǎn)錄組學(xué)研究.藥學(xué)學(xué)報(bào),2010,45(4):524–529.

    [19]Sun C,Li Y,Wu Q,Luo HM,Sun YZ,Song JY,Lui E MK,Chen SL.De novosequencing and analysis of the American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis.BMC Genomics,2010,11(1):262.

    [20]Logacheva MD,Kasianov AS,Vinogradov DV,Samigullin TH,Gelfand MS,Makeev VJ,Penin AA.De novosequencing and characterization of floral transcriptome in two species of buckwheat(Fagopyrum).BMC Genomics,2011,12(1):30.

    [21]Francis RW,Thompson-Wicking K,Carter KW,Anderson D,Kees UR,Beesley AH.FusionFinder:A software tool to identify expressed gene fusion candidates from RNA-Seq data.PLoS ONE,2012,7(6):e39987.

    [22]Li PH,Ponnala L,Gandotra N,Wang L,Si YQ,Tausta SL,Kebrom TH,Provart N,Patel R,Myers CR,Reidel EJ,Turgeon R,Liu P,Sun Q,Nelson T,Brutnell TP.The developmental dynamics of the maize leaf transcriptome.Nat Genet,2010,42(12):1060–1067.

    [23]Grabherr MG,Haas BJ,Yassour M,Levin JZ,Thompson DA,Amit I,Adiconis X,Fan L,Raychowdhury R,Zeng Q,Chen Z,Mauceli E,Hacohen N,Gnirke A,Rhind N,Palma F,Birren BW,Nusbaum C,Lindblad-Toh K,Friedma N,Regev A.Full-length transcriptome assembly from RNA-Seq data without a reference genome.Nat Biotechnol,2011,29(7):644–652.

    [24]Conesa A,G?tz S,García-Gómez JM,Terol J,Talón M,Robles M.Blast2GO:a universal tool for annotation,visualization and analysis in functional genomics research.Bioinformatics,2005,21(18):3674–3676.

    [25]Ye J,Fang L,Zheng HK,Zhang Y,Chen J,Zhang ZG,Wang J,Li ST,Li RQ,Bolund L,Wang J.WEGO:a web tool for plotting GO annotations.Nucleic Acids Res,2006,34(1):W293–297.

    [26]林萍,曹永慶,姚小華,王開(kāi)良,滕建華.普通油茶種子 4個(gè)發(fā)育時(shí)期的轉(zhuǎn)錄組分析.分子植物育種,2011,9(4):498–505.

    [27]Iseli C,Jongeneel CV,Bucher P.ESTScan:a program for detecting,evaluating,and reconstructing potential coding regions in EST sequences.Proc Int Conf Intell Syst Mol Biol,1999,138–148.

    [28]Kanehisa M,Goto S,Kawashima S,Okuno Y,Hattori M.The KEGG resource for deciphering the genome.Nucleic Acids Res,2004,32(1):D277–D280.

    [29]梁燁,陳雙燕,劉公社.新一代測(cè)序技術(shù)在植物轉(zhuǎn)錄組研究中的應(yīng)用.遺傳,2011,33(12):1317–1326.

    [30]施季森,王占軍,陳金慧.木本植物全基因組測(cè)序研究進(jìn)展.遺傳,2012,34(2):145–156.

    [31]Li X,Chen GH,Zhang WY,Zhang X.Genome-wide transcriptional analysis of maize endosperm in response to ae wx double mutations.J Genet Genomics,2010,37(11):749–762.

    [32]Wu Q,Sun C,Luo H,Li Y,Niu Y,Sun Y,Lu A,Chen S.Transcriptome analysis ofTaxus cuspidataneedles based on 454 pyrosequencing.Planta Med,2010,77(4):394–400.

    [33]Luo HM,Li Y,Sun C,Wu Q,Song JY,Sun YZ,Steinmetz A,Chen SL.Comparison of 454-ESTs fromHuperzia serrataandPhlegmariurus carinatusreveals putative genes involved in lycopodium alkaloid biosynthesis and developmental regulation.BMC Plant Biol,2010,10(1):209.

    猜你喜歡
    芳樟醇樟樹(shù)測(cè)序
    香樟樹(shù)
    杰 Sir 帶你認(rèn)識(shí)宏基因二代測(cè)序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    氧化芳樟醇合成及應(yīng)用的研究進(jìn)展*
    廣州化工(2021年19期)2021-10-25 14:03:02
    二代測(cè)序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    香樟樹(shù),樟樹(shù)香
    樟樹(shù)扦插繁殖研究進(jìn)展
    花椒酒中檸檬烯和芳樟醇的測(cè)定
    從八角茴香油前餾分中單離芳樟醇和草蒿腦工藝研究
    芳樟型樟樹(shù)葉精油減壓連續(xù)精餾分離芳樟醇的工藝模擬
    基因捕獲測(cè)序診斷血癌
    欧美日韩av久久| 天天影视国产精品| 亚洲欧美清纯卡通| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产欧美日韩av| 99久久99久久久精品蜜桃| 久久ye,这里只有精品| 国产精品国产三级专区第一集| 成年动漫av网址| 亚洲精品第二区| 亚洲 欧美一区二区三区| 各种免费的搞黄视频| 国产成人精品久久久久久| 国产不卡av网站在线观看| 女的被弄到高潮叫床怎么办| 亚洲情色 制服丝袜| 最黄视频免费看| 色94色欧美一区二区| 国产不卡av网站在线观看| 欧美精品高潮呻吟av久久| 国产 精品1| 午夜免费鲁丝| 国产片内射在线| 欧美中文综合在线视频| 中文精品一卡2卡3卡4更新| 一区二区三区精品91| 午夜福利,免费看| 国产精品国产三级国产专区5o| 中文字幕色久视频| 美女中出高潮动态图| 秋霞在线观看毛片| 一区在线观看完整版| 欧美日韩综合久久久久久| 下体分泌物呈黄色| 免费观看性生交大片5| 一区在线观看完整版| 色94色欧美一区二区| 日韩,欧美,国产一区二区三区| 一区福利在线观看| 免费黄网站久久成人精品| 亚洲自偷自拍图片 自拍| 日本av免费视频播放| 成人毛片60女人毛片免费| 久久精品aⅴ一区二区三区四区| 欧美成人精品欧美一级黄| 黄网站色视频无遮挡免费观看| 啦啦啦啦在线视频资源| 黄色视频在线播放观看不卡| 国产黄色视频一区二区在线观看| 少妇人妻精品综合一区二区| 丰满迷人的少妇在线观看| 在线天堂中文资源库| 免费日韩欧美在线观看| 最近中文字幕高清免费大全6| 日本wwww免费看| 亚洲av国产av综合av卡| 成年av动漫网址| 伦理电影免费视频| 亚洲欧美成人综合另类久久久| 少妇人妻 视频| 自拍欧美九色日韩亚洲蝌蚪91| 最近最新中文字幕免费大全7| 欧美日韩亚洲国产一区二区在线观看 | 这个男人来自地球电影免费观看 | 成人影院久久| 国产一区亚洲一区在线观看| 大香蕉久久成人网| 国产免费福利视频在线观看| 日本91视频免费播放| 国产精品熟女久久久久浪| 中文乱码字字幕精品一区二区三区| 又黄又粗又硬又大视频| 国产在视频线精品| 成年人午夜在线观看视频| 中文精品一卡2卡3卡4更新| 亚洲熟女精品中文字幕| 中文字幕人妻熟女乱码| 亚洲av成人精品一二三区| 51午夜福利影视在线观看| 成人亚洲精品一区在线观看| 亚洲国产欧美在线一区| av网站免费在线观看视频| av在线观看视频网站免费| 欧美日韩福利视频一区二区| av福利片在线| 中文字幕色久视频| 精品福利永久在线观看| 久久久国产欧美日韩av| 亚洲熟女毛片儿| 五月天丁香电影| xxxhd国产人妻xxx| 欧美另类一区| 波多野结衣一区麻豆| 母亲3免费完整高清在线观看| 国产在线视频一区二区| 不卡视频在线观看欧美| 久久婷婷青草| 999精品在线视频| 中文乱码字字幕精品一区二区三区| 日日爽夜夜爽网站| 黄色怎么调成土黄色| 精品久久蜜臀av无| 成人国产av品久久久| 精品人妻一区二区三区麻豆| 国产日韩欧美视频二区| 国产一级毛片在线| 麻豆乱淫一区二区| 日本爱情动作片www.在线观看| 人人澡人人妻人| 七月丁香在线播放| 国产成人免费观看mmmm| 成人毛片60女人毛片免费| 日韩av不卡免费在线播放| 夫妻性生交免费视频一级片| 深夜精品福利| 巨乳人妻的诱惑在线观看| av在线老鸭窝| 日韩制服丝袜自拍偷拍| 男人舔女人的私密视频| 亚洲精品在线美女| 欧美成人午夜精品| 一边摸一边做爽爽视频免费| 在现免费观看毛片| 建设人人有责人人尽责人人享有的| 一级a爱视频在线免费观看| 日本欧美国产在线视频| 韩国av在线不卡| 国产欧美日韩一区二区三区在线| 欧美av亚洲av综合av国产av | 亚洲美女视频黄频| 黑人猛操日本美女一级片| 日韩成人av中文字幕在线观看| 亚洲av男天堂| 国产精品 欧美亚洲| 国产成人av激情在线播放| 超碰成人久久| 亚洲欧美一区二区三区黑人| 80岁老熟妇乱子伦牲交| 久久av网站| 黄色视频不卡| 精品国产露脸久久av麻豆| 国产成人精品久久二区二区91 | 大香蕉久久成人网| 欧美精品人与动牲交sv欧美| 人体艺术视频欧美日本| 亚洲成色77777| 性高湖久久久久久久久免费观看| 99国产精品免费福利视频| 少妇精品久久久久久久| av在线老鸭窝| avwww免费| 亚洲国产精品999| 欧美激情高清一区二区三区 | 黄色毛片三级朝国网站| 超碰成人久久| 熟妇人妻不卡中文字幕| 久久国产精品男人的天堂亚洲| 亚洲成人免费av在线播放| av天堂久久9| 五月天丁香电影| 日本欧美视频一区| 成人免费观看视频高清| 欧美日韩一级在线毛片| 男男h啪啪无遮挡| 在线观看www视频免费| 久久久精品区二区三区| 宅男免费午夜| 日韩一本色道免费dvd| 国产精品 欧美亚洲| 亚洲欧美激情在线| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 日本午夜av视频| 啦啦啦在线免费观看视频4| 水蜜桃什么品种好| 搡老岳熟女国产| 久久久久久久久久久免费av| 黄片无遮挡物在线观看| 国产成人欧美在线观看 | 婷婷色av中文字幕| 十八禁高潮呻吟视频| 少妇被粗大的猛进出69影院| 欧美激情极品国产一区二区三区| 在线看a的网站| 欧美日韩视频精品一区| 国产成人欧美在线观看 | 99国产精品免费福利视频| 亚洲第一区二区三区不卡| 亚洲av福利一区| 欧美日韩国产mv在线观看视频| 一区二区三区激情视频| 国产又色又爽无遮挡免| 伊人亚洲综合成人网| 色吧在线观看| 啦啦啦中文免费视频观看日本| 天堂俺去俺来也www色官网| 国产一区二区 视频在线| 久久精品久久精品一区二区三区| 91老司机精品| 啦啦啦在线免费观看视频4| av电影中文网址| 2021少妇久久久久久久久久久| 五月天丁香电影| 国产精品免费视频内射| 黑人巨大精品欧美一区二区蜜桃| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品男人的天堂亚洲| kizo精华| 99re6热这里在线精品视频| 久久久久精品人妻al黑| 精品国产乱码久久久久久男人| 亚洲一码二码三码区别大吗| 观看av在线不卡| 欧美精品av麻豆av| 一本久久精品| 最近的中文字幕免费完整| 精品亚洲乱码少妇综合久久| 成人亚洲精品一区在线观看| 亚洲一区二区三区欧美精品| 日本wwww免费看| 黄色a级毛片大全视频| 曰老女人黄片| av免费在线观看网站| 亚洲国产高清在线一区二区三 | 亚洲人成伊人成综合网2020| 十分钟在线观看高清视频www| 亚洲av电影不卡..在线观看| 久久 成人 亚洲| 国产精品1区2区在线观看.| 女人被躁到高潮嗷嗷叫费观| 少妇被粗大的猛进出69影院| 国产精品自产拍在线观看55亚洲| 国产伦人伦偷精品视频| 少妇裸体淫交视频免费看高清 | 亚洲国产中文字幕在线视频| 国产成人啪精品午夜网站| 91大片在线观看| 美国免费a级毛片| 国产精品一区二区精品视频观看| 亚洲五月色婷婷综合| 乱人伦中国视频| 美女免费视频网站| 亚洲av日韩精品久久久久久密| 黄色视频,在线免费观看| 日韩av在线大香蕉| 又大又爽又粗| 日韩有码中文字幕| 女人被狂操c到高潮| svipshipincom国产片| 91老司机精品| 国产亚洲精品av在线| 欧美日韩瑟瑟在线播放| 一边摸一边做爽爽视频免费| 韩国精品一区二区三区| 一进一出好大好爽视频| 深夜精品福利| 嫩草影院精品99| 日韩一卡2卡3卡4卡2021年| 精品国产乱子伦一区二区三区| 欧美日韩精品网址| 免费看a级黄色片| 国产一区二区激情短视频| 在线免费观看的www视频| 精品久久久久久久毛片微露脸| 亚洲aⅴ乱码一区二区在线播放 | 又黄又粗又硬又大视频| 动漫黄色视频在线观看| 女人高潮潮喷娇喘18禁视频| 国产私拍福利视频在线观看| 日本a在线网址| 热99re8久久精品国产| 国产成人av激情在线播放| 在线观看66精品国产| 国产又爽黄色视频| 精品久久久久久久毛片微露脸| 国产成人精品久久二区二区91| 女警被强在线播放| 黄色片一级片一级黄色片| 美女扒开内裤让男人捅视频| 欧美日本视频| 精品一区二区三区视频在线观看免费| 波多野结衣高清无吗| 宅男免费午夜| 免费搜索国产男女视频| 国产亚洲欧美98| 国产精品野战在线观看| 美女国产高潮福利片在线看| 黄网站色视频无遮挡免费观看| 韩国精品一区二区三区| 欧美黑人精品巨大| 日韩三级视频一区二区三区| 精品国产乱码久久久久久男人| 1024视频免费在线观看| 精品国产一区二区久久| 亚洲国产中文字幕在线视频| 日日摸夜夜添夜夜添小说| 国产伦一二天堂av在线观看| 激情视频va一区二区三区| 午夜影院日韩av| 亚洲第一青青草原| 香蕉丝袜av| 亚洲成av人片免费观看| 狂野欧美激情性xxxx| 日韩大尺度精品在线看网址 | 亚洲av美国av| 欧美日韩亚洲综合一区二区三区_| 丝袜人妻中文字幕| 女人精品久久久久毛片| 欧美久久黑人一区二区| 欧美乱色亚洲激情| 激情视频va一区二区三区| 亚洲美女黄片视频| 国产一区二区在线av高清观看| 黄片大片在线免费观看| 欧美色视频一区免费| 在线播放国产精品三级| 欧美最黄视频在线播放免费| 岛国在线观看网站| 神马国产精品三级电影在线观看 | 麻豆久久精品国产亚洲av| 日韩成人在线观看一区二区三区| 亚洲av日韩精品久久久久久密| 中亚洲国语对白在线视频| 亚洲国产精品久久男人天堂| av电影中文网址| 国产又色又爽无遮挡免费看| 大陆偷拍与自拍| 国产精品美女特级片免费视频播放器 | 国产一区二区在线av高清观看| 给我免费播放毛片高清在线观看| 美女扒开内裤让男人捅视频| 亚洲国产日韩欧美精品在线观看 | 女人精品久久久久毛片| 欧美乱色亚洲激情| 日本欧美视频一区| 99国产精品一区二区蜜桃av| 女警被强在线播放| 欧美大码av| 亚洲av片天天在线观看| 国产成人av教育| 级片在线观看| 在线十欧美十亚洲十日本专区| 91成人精品电影| 99在线人妻在线中文字幕| 1024视频免费在线观看| 伦理电影免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 婷婷精品国产亚洲av在线| 琪琪午夜伦伦电影理论片6080| 十八禁人妻一区二区| 午夜福利成人在线免费观看| 亚洲av成人一区二区三| 国产精品综合久久久久久久免费 | 中文字幕另类日韩欧美亚洲嫩草| 在线av久久热| 国产成+人综合+亚洲专区| 国产精品香港三级国产av潘金莲| 人成视频在线观看免费观看| www.精华液| 精品国产亚洲在线| 亚洲精品美女久久av网站| 一本久久中文字幕| 欧美成人免费av一区二区三区| 欧美性长视频在线观看| 国产精品久久久久久亚洲av鲁大| 久久精品人人爽人人爽视色| 一级作爱视频免费观看| 夜夜躁狠狠躁天天躁| 亚洲国产高清在线一区二区三 | 国产精品国产高清国产av| 亚洲国产精品999在线| 亚洲九九香蕉| www国产在线视频色| 成人国产综合亚洲| 国产精品免费一区二区三区在线| 日韩高清综合在线| 亚洲国产精品久久男人天堂| 午夜激情av网站| 免费在线观看完整版高清| 好男人在线观看高清免费视频 | 国产av一区在线观看免费| 日韩中文字幕欧美一区二区| 国产精品免费一区二区三区在线| 麻豆av在线久日| 最新美女视频免费是黄的| 亚洲国产欧美日韩在线播放| 十八禁网站免费在线| 亚洲欧美激情综合另类| 中文字幕高清在线视频| 国产野战对白在线观看| 亚洲国产精品成人综合色| 亚洲国产中文字幕在线视频| 日韩中文字幕欧美一区二区| 国产精品爽爽va在线观看网站 | 久久精品人人爽人人爽视色| 一级作爱视频免费观看| 亚洲欧洲精品一区二区精品久久久| 国产成人欧美| 人人妻人人爽人人添夜夜欢视频| 青草久久国产| 国产av一区在线观看免费| 变态另类成人亚洲欧美熟女 | 亚洲成a人片在线一区二区| 亚洲国产欧美网| 无遮挡黄片免费观看| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 黄色成人免费大全| 中文字幕人成人乱码亚洲影| 欧美在线一区亚洲| 黄色丝袜av网址大全| 亚洲一区二区三区色噜噜| 国产成年人精品一区二区| 热99re8久久精品国产| 一边摸一边抽搐一进一出视频| 久久中文看片网| 国产区一区二久久| 日韩中文字幕欧美一区二区| 午夜福利免费观看在线| 无人区码免费观看不卡| 日日夜夜操网爽| 精品一品国产午夜福利视频| 给我免费播放毛片高清在线观看| 韩国av一区二区三区四区| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| 女人精品久久久久毛片| 午夜两性在线视频| 日本精品一区二区三区蜜桃| 国产熟女午夜一区二区三区| 国产成人免费无遮挡视频| 国产欧美日韩一区二区三| 亚洲一卡2卡3卡4卡5卡精品中文| 好看av亚洲va欧美ⅴa在| 一区二区日韩欧美中文字幕| 99精品久久久久人妻精品| 国产免费男女视频| 日韩欧美三级三区| 中出人妻视频一区二区| 国产又爽黄色视频| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 成人三级做爰电影| 成年版毛片免费区| 亚洲激情在线av| 久99久视频精品免费| 国产成人精品久久二区二区免费| 国产av一区二区精品久久| 国产成人影院久久av| 淫妇啪啪啪对白视频| 久久欧美精品欧美久久欧美| 国产伦人伦偷精品视频| 国产激情欧美一区二区| 色在线成人网| 亚洲精华国产精华精| 国产精品野战在线观看| 91成人精品电影| 午夜精品在线福利| 国产成人精品久久二区二区免费| 久久国产乱子伦精品免费另类| www.www免费av| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久亚洲av鲁大| 亚洲精华国产精华精| 色综合婷婷激情| 色综合亚洲欧美另类图片| 女警被强在线播放| 手机成人av网站| 淫秽高清视频在线观看| 久久国产精品人妻蜜桃| 村上凉子中文字幕在线| 午夜福利视频1000在线观看 | 国产精品二区激情视频| 亚洲第一欧美日韩一区二区三区| 亚洲精品在线美女| 亚洲av熟女| 精品高清国产在线一区| 亚洲第一av免费看| 一区在线观看完整版| 一边摸一边抽搐一进一出视频| 日本精品一区二区三区蜜桃| 极品人妻少妇av视频| 成年版毛片免费区| 国产欧美日韩精品亚洲av| 亚洲免费av在线视频| 久热爱精品视频在线9| 久久久国产精品麻豆| 多毛熟女@视频| 此物有八面人人有两片| 电影成人av| 国产精品乱码一区二三区的特点 | 91麻豆av在线| 亚洲欧美日韩另类电影网站| 国产黄a三级三级三级人| 人妻久久中文字幕网| 午夜视频精品福利| 免费搜索国产男女视频| 亚洲中文日韩欧美视频| 亚洲片人在线观看| 美女 人体艺术 gogo| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 亚洲 欧美一区二区三区| 久热这里只有精品99| x7x7x7水蜜桃| 麻豆成人av在线观看| 日韩欧美免费精品| 免费观看精品视频网站| 丰满的人妻完整版| 亚洲av电影不卡..在线观看| 少妇粗大呻吟视频| 99国产精品一区二区蜜桃av| 一级毛片精品| 国产91精品成人一区二区三区| 可以免费在线观看a视频的电影网站| 麻豆久久精品国产亚洲av| 人成视频在线观看免费观看| 黑人操中国人逼视频| 亚洲成人国产一区在线观看| 99久久精品国产亚洲精品| 一级黄色大片毛片| 国产av精品麻豆| av天堂久久9| 久久国产乱子伦精品免费另类| 如日韩欧美国产精品一区二区三区| 超碰成人久久| 19禁男女啪啪无遮挡网站| 亚洲国产精品久久男人天堂| 亚洲伊人色综图| 男人的好看免费观看在线视频 | 亚洲视频免费观看视频| 最新在线观看一区二区三区| 亚洲精品国产区一区二| 69av精品久久久久久| 波多野结衣巨乳人妻| 日韩欧美一区二区三区在线观看| 99久久精品国产亚洲精品| 国产一区二区三区视频了| cao死你这个sao货| 9191精品国产免费久久| 国产精品一区二区三区四区久久 | 国产三级黄色录像| 黑人欧美特级aaaaaa片| 欧美在线黄色| 国产麻豆69| 国产成人欧美| 久久性视频一级片| 精品不卡国产一区二区三区| 免费搜索国产男女视频| 亚洲一区中文字幕在线| 亚洲精品中文字幕一二三四区| 午夜福利影视在线免费观看| 国产精品日韩av在线免费观看 | 国产私拍福利视频在线观看| 亚洲国产中文字幕在线视频| 999久久久精品免费观看国产| 日韩av在线大香蕉| 国产成+人综合+亚洲专区| 国产精品精品国产色婷婷| 国产成+人综合+亚洲专区| 精品国产乱码久久久久久男人| 亚洲精品av麻豆狂野| 国产三级黄色录像| 50天的宝宝边吃奶边哭怎么回事| 久久亚洲精品不卡| 成年女人毛片免费观看观看9| 国产三级在线视频| 性少妇av在线| 久久人人97超碰香蕉20202| 国产av又大| 欧美日韩黄片免| 国产成人影院久久av| 亚洲精品一卡2卡三卡4卡5卡| 狠狠狠狠99中文字幕| 国产精品一区二区精品视频观看| 国产成人精品久久二区二区91| 国产伦一二天堂av在线观看| 91成年电影在线观看| 亚洲中文日韩欧美视频| 日韩视频一区二区在线观看| 国产1区2区3区精品| 国产一区二区三区综合在线观看| 人妻丰满熟妇av一区二区三区| 天堂√8在线中文| 国产1区2区3区精品| videosex国产| 成人欧美大片| 国产亚洲欧美在线一区二区| 亚洲午夜理论影院| 久久中文字幕人妻熟女| 91九色精品人成在线观看| 免费女性裸体啪啪无遮挡网站| 天堂√8在线中文| 久久影院123| 欧美黑人欧美精品刺激| 亚洲 欧美一区二区三区| 国产黄a三级三级三级人| 久久人人爽av亚洲精品天堂| 国产91精品成人一区二区三区| svipshipincom国产片| 国产精品久久久久久精品电影 | 精品国产一区二区久久| 亚洲全国av大片| 99久久国产精品久久久| 99久久久亚洲精品蜜臀av| 99精品久久久久人妻精品| 嫩草影视91久久| 国产麻豆69| 国产真人三级小视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 巨乳人妻的诱惑在线观看| 精品久久久久久成人av| 国产欧美日韩一区二区三| 久久国产精品人妻蜜桃| 91精品三级在线观看|