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

    基于葉綠體基因rbcL和trnL-F序列的河北省野生苔草屬親緣關(guān)系分析

    2024-12-31 00:00:00周欣瑩董彥娜邱玉鵬孫浩男劉亞偉李艷劉冬云
    草地學(xué)報 2024年9期

    摘要:為研究河北省野生苔草屬(Carex)植物的親緣關(guān)系,本研究以43種河北省野生苔草及國外引種的27份苔草種質(zhì)為試驗(yàn)材料,克隆了葉綠體基因rbcL和trnL-F序列,進(jìn)行遺傳多樣性及親緣關(guān)系分析,并構(gòu)建了系統(tǒng)發(fā)育樹。結(jié)果表明,從70份苔草屬植物中成功克隆到的rbcL和trnL-F序列數(shù)量分別為67條、60條,測序成功率分別為95.7%、85.7%,經(jīng)正反序列拼接對齊獲得序列長度分別為716 bp、976 bp,并上傳至NCBI獲得登錄號。二者合并得到59條序列,長度為1692 bp。rbcL,trnL-F以及合并序列分別定義了43,56,58個單倍型?;诤喜⑿蛄械?9份種質(zhì)間遺傳距離為0~0.0421,平均遺傳距離為0.0142,整體劃分為3大類,與形態(tài)學(xué)分類基本相符。trnL-F序列有望開發(fā)為DNA條形碼,以期為苔草屬的分類鑒定、系統(tǒng)發(fā)育等研究提供依據(jù)。

    關(guān)鍵詞:苔草屬;rbcL;trnL-F;親緣關(guān)系;遺傳多樣性

    中圖分類號:Q37""" 文獻(xiàn)標(biāo)識碼:A""""" 文章編號:1007-0435(2024)09-2728-09

    收稿日期:2024-01-15;修回日期:2024-04-06

    基金項(xiàng)目:河北省苔草屬植物種質(zhì)資源收集與新種質(zhì)創(chuàng)制項(xiàng)目(21326348D)資助

    作者簡介:

    周欣瑩(1999-),女,漢族,河北保定人,碩士研究生,主要從事園林植物遺傳育種與栽培研究,E-mail: 1240260550@qq.com;*通信作者Author for correspondence, E-mail: 123360961@qq.com; liudongyun0505@163.com

    doi:10.11733/j.issn.1007-0435.2024.09.006

    引用格式:

    周欣瑩, 董彥娜, 邱玉鵬,等.基于葉綠體基因rbcL和trnL-F序列的河北省野生苔草屬親緣關(guān)系分析[J].草地學(xué)報,2024,32(9):2728-2736

    ZHOU Xin-ying, DONG Yan-na, QIU Yu-peng,et al.Genetic Relationship Analysis of Wild Carex Species in Hebei Province Based on Sequences of Chloroplast Genes rbcL and trnL-F[J].Acta Agrestia Sinica,2024,32(9):2728-2736

    Genetic Relationship Analysis of Wild Carex Species in Hebei Province

    Based on Sequences of Chloroplast Genes rbcL and trnL-F

    ZHOU Xin-ying1, DONG Yan-na2, QIU Yu-peng1, SUN Hao-nan1, LIU Ya-wei1,

    LI Yan1*, LIU Dong-yun1*

    (1. Collage of Landscape and Travel, Hebei Agricultural University, Baoding, Hebei Province 071000, China;

    2. Shijiazhuang Information Engineering Vocational College, Shijiazhuang, Hebei Province 050000, China)

    Abstract:To study the genetic relationships of wild Carex plants in Hebei Province,43 species of wild Carex plants in Hebei Province and 27 germplasms of Carex introduced from abroad were used as experimental materials. The chloroplast genes rbcL and trnL-F sequences were cloned for genetic diversity and phylogenetic analyses. And a phylogenetic tree was constructed. The results showed that the number of rbcL and trnL-F sequences successfully cloned from 70 plants of the genus Carex was 67 and 60,with sequencing success rates of 95.7% and 85.7%,respectively. The sequence lengths obtained after forward and backward sequence splicing and alignment were 716 bp and 976 bp, respectively and were uploaded to NCBI to obtain the GenBank accession number. The rbcL and trnL-F sequences combination obtained 59 sequences with a length of 1 692 bp. The rbcL,trnL-F and combined sequence define 43,56 and 58 haplotypes. Based on the combined sequence,the genetic distance between 59 germplasms is 0~0.0421,with an average genetic distance of 0.0142. Based on NJ phylogenetic tree of combined sequence,59 germplasms could be divided into 3 groups,which was basically consistent with morphological classification. The trnL-F sequence is expected to be developed as DNA barcodes,in order to provide basis for the classification,identification and phylogenetic research of Carex.

    Key words:Carex; rbcL;trnL-F;Genetic Relationship;Genetic Diversity

    莎草科(Cyperaceae)苔草屬(Carex)植物種類豐富,在世界范圍分布廣泛,全球約有2000種,我國約有500種,其中,河北省內(nèi)約有45種[1-3]。苔草屬植物春季返青早、生長期長,具有耐寒、耐旱、耐瘠薄、耐踐踏等特點(diǎn),且自繁能力及適應(yīng)性強(qiáng),有一定觀賞價值,能夠作為替代傳統(tǒng)草坪的優(yōu)質(zhì)地被植物[4]。同時,苔草屬植物也在生態(tài)修復(fù)方面發(fā)揮著重要的作用[5]。根據(jù)形態(tài)學(xué)分類,苔草屬植物可分為二柱苔草亞屬(Carex Subg. Vignea)、復(fù)序苔草亞屬(Carex Subg. Indocarex)和苔草亞屬(Carex Subg. Carex),二柱苔草亞屬和苔草亞屬的分布較為廣泛,而復(fù)序苔草亞屬主要分布在長江以南地區(qū)[7-8],河北省野生苔草屬于苔草亞屬和二柱苔草亞屬。眾多學(xué)者發(fā)現(xiàn),苔草屬具有形態(tài)復(fù)雜、高度變異等特點(diǎn),因此,苔草屬內(nèi)的分類工作是一大難題[9-10]。近年來,隨著對苔草屬研究的深入,國內(nèi)外不斷發(fā)現(xiàn)苔草屬新物種,此外,苔草屬的遺傳背景復(fù)雜、親緣關(guān)系不明確等原因也限制了引種和育種工作的進(jìn)行,阻礙了進(jìn)一步推廣和開發(fā)利用。

    葉綠體有一套完整的遺傳系統(tǒng),其獨(dú)立的遺傳信息被稱為葉綠體基因組(cpDNA),大多數(shù)植物遵循母系遺傳,總體保守性較高,多用于植物種間遺傳信息分析[11]。cpDNA編碼區(qū)和非編碼區(qū)在進(jìn)化速率上存在差異,非編碼區(qū)的進(jìn)化速率更快,系統(tǒng)發(fā)育信息含量更多[12]。因此,在進(jìn)行系統(tǒng)發(fā)育研究時,需要考慮所選基因序列的進(jìn)化速率,篩選合適的基因序列[13]。常用于植物系統(tǒng)發(fā)育研究的葉綠體基因序列有rbcL,rpl16,rpoC1,trnL-F,trnT-L,psbA-trnH等[14-16],其不僅可以用于構(gòu)建系統(tǒng)發(fā)育樹還可進(jìn)一步開發(fā)為DNA條形碼。Oda等[17]基于多個葉綠體基因和核基因序列(ITS)數(shù)據(jù)進(jìn)行了分子系統(tǒng)發(fā)育分析,了解日本苔草屬Rarae組種間關(guān)系,并結(jié)合其形態(tài)特征進(jìn)一步建立了該組的分類檢索表。Starr等[18]利用葉綠體基因開發(fā)了DNA條形碼,發(fā)現(xiàn)matK條形碼可區(qū)分47%的試驗(yàn)材料,并通過該片段發(fā)現(xiàn)新物種。Yano等[19]基于atpB-rbcL,trnT-trnL,trnL內(nèi)含子和trnL-trnF序列,分析了C. conica及其近緣種的分子地理學(xué)關(guān)系,在258份材料中得到了16個單倍型,將其分為日本東部和西部兩大地理結(jié)構(gòu)類型。

    野生苔草屬植物經(jīng)過長期自然進(jìn)化,其遺傳背景復(fù)雜,極大地影響研究和開發(fā)利用。河北省有著豐富的苔草屬種質(zhì)資源,對其進(jìn)行系統(tǒng)性的調(diào)查與引種后,發(fā)現(xiàn)部分種質(zhì)的適應(yīng)性強(qiáng),且株型、葉色觀賞價值較高,有待進(jìn)一步開發(fā)為優(yōu)質(zhì)的草坪地被或觀賞草。本研究以43種河北省野生苔草及國外引種的27份種質(zhì)為試材,克隆和分析了其葉綠體基因rbcL和trnL-F序列,以此進(jìn)行了親緣關(guān)系和遺傳多樣性分析,并進(jìn)一步鑒定了種質(zhì)資源,旨在從分子角度上為苔草屬的分類鑒定、DNA條形碼等研究提供理論依據(jù),推動苔草屬植物的開發(fā)利用。

    1" 材料方法

    1.1" 試驗(yàn)材料

    供試材料包括43種河北省野生苔草(根據(jù)《河北省植物志》及《內(nèi)蒙古植物志》鑒定),采自保定市洪崖山、承德市霧靈山、張家口市小五臺山等地,以及27份國外引種的苔草種質(zhì),共計70份種質(zhì)。引種收集后統(tǒng)一栽植于河北農(nóng)業(yè)大學(xué)試驗(yàn)基地,詳見表1。

    1.2" 試驗(yàn)方法

    1.2.1" DNA提取和PCR擴(kuò)增" 使用DNA提取試劑盒(CW0531S,康為世紀(jì))提取70份苔草樣本的DNA用Nano Drop微量分光光度計進(jìn)行DNA純度檢測并使用1.2%的瓊脂糖凝膠進(jìn)行電泳檢測。使用ddH2O將DNA母液稀釋至80 ng·μL-1,-20℃保存。

    rbcL序列使用通用引物rbcL-1f、rbcL-724r[20],trnL-F間隔區(qū)序列使用通用引物trnL-c、trnF-f[21],詳情見表2。

    PCR擴(kuò)增體系采用20 μL體系:1 μL DNA模板,1 μL上下游引物(10 μmol·L-1),10 μL 2×Super Pfx Master Mix,7 μL ddH2O。

    擴(kuò)增程序:98℃預(yù)變性3 min,98℃變性10 s,退火30 s,72℃延伸1 min,35個循環(huán),72℃延伸10 min,4℃保存。

    1.2.2" 擴(kuò)增產(chǎn)物的回收與純化" PCR產(chǎn)物采用1.2%的瓊脂糖凝膠電泳進(jìn)行檢測,待條帶分離后使用膠回收試劑盒(DC3511,BIOMIGA)進(jìn)行回收、純化,送上海生物公司進(jìn)行雙向測序。

    1.2.3" 數(shù)據(jù)分析" 利用BioEdite 7對雙向峰圖進(jìn)行校對、拼接,去除兩端峰圖雜亂、質(zhì)量較低的序列,將拼接后的序列作為查詢序列輸入到GenBank數(shù)據(jù)庫中,利用NCBI-BLAST相似性比對工具進(jìn)行序列比對。利用MEGA 7軟件對序列的基本特點(diǎn)進(jìn)行統(tǒng)計分析,包括了位點(diǎn)長度、簡約信息位點(diǎn)(Parsimony information sites,Ps)、變異位點(diǎn)(Variable sites,Vs)、保守位點(diǎn)(Conservedsites,C),堿基組成(Nucleotide composition)等,采用鄰接法(Neighbor-Joining,NJ)進(jìn)行建樹?;贙2P(Kimura 2-parameter model)雙參數(shù)模型計算遺傳距離。利用DnaSP 5.0統(tǒng)計序列的特征信息和變異位點(diǎn)信息,Network 10.0基于MJ算法構(gòu)建單倍型網(wǎng)絡(luò)圖。最終以同為莎草科的具芒碎米莎草(Cyperus microiria登錄號MN205054、KR869862)為外類群,構(gòu)建系統(tǒng)發(fā)育樹。

    2" 結(jié)果與分析

    2.1" PCR擴(kuò)增結(jié)果

    PCR擴(kuò)增產(chǎn)物經(jīng)1.2%瓊脂糖凝膠電泳檢測,所有樣品的rbcL、trnL-F序列目的條帶分別在750~1000 bp、1000~1500 bp,擴(kuò)增出的條帶清晰、明亮,無明顯雜帶,可用于后續(xù)的測序試驗(yàn)。部分條帶結(jié)果見圖1,2。

    2.2" rbcL、trnL-F序列特征分析

    將rbcL和trnL-F序列測序結(jié)果進(jìn)行分析,得到了67條、60條序列信息,測序成功率分別為95.7%,85.7%。正反序列經(jīng)拼接并去除兩端不準(zhǔn)確側(cè)翼序列后,分別獲得長度為716 bp、976 bp。

    序列基本信息見表3。獲得的rbcL和trnL-F序列合并后得到59條序列,長度為1692 bp。核苷酸構(gòu)成分析結(jié)果顯示,rbcL和trnL-F序列堿基A+T平均含量均高于G+C含量,rbcL、trnL-F及合并序列的平均堿基頻率為A=30.5%,39.6%,35.8%,T=28.6%,35.2%,32.3%;C=20.0%,12.4%,15.6%;G=20.9%,12.8%,16.3%。trnL-F及合并序列的平均GC含量低于rbcL序列,但變異位點(diǎn)、單一突變位點(diǎn)及簡約信息位點(diǎn)均高于rbcL序列。將所測得苔草屬植物rbcL和trnL-F序列提交到NCBI,并獲得登錄號為PP531050~PP531176(表1)。

    2.3" rbcL、trnL-F單倍型多樣性分析

    對rbcL、trnL-F以及合并序列進(jìn)行多態(tài)性分析,詳見表4。rbcL、trnL-F序列的單倍型數(shù)目(h)為43,56個,合并序列的單倍型數(shù)目為58個,多于rbcL和trnL-F序列的單倍型數(shù)目。苔草屬植物基于單一序列和合并序列分析的單倍型多態(tài)性(Hd)均大于0.5,核苷酸多態(tài)性(π)均高于0.005,合并序列的單倍型多態(tài)性高于單一序列。

    通過Tajima’ s D和Fu and Li’ s值來檢驗(yàn)rbcL、trnL-F序列在物種水平的進(jìn)化過程中是否遵循中性進(jìn)化模型,得到的Tajima’ s D值分別為-2.305,-2.137,-2.192,F(xiàn)u and Li’ s D*值為-4.025,-2.583,-2.783,F(xiàn)u and Li’ s F*值為-4.023,-2.887,-3.058,結(jié)果均顯著,不遵循中性進(jìn)化模型,合并序列分析的35種河北省野生苔草屬植物Tajima’ s D和Fu and Li’ s值均為負(fù)值且結(jié)果顯著,也不遵循中性進(jìn)化模型。

    通過DnaSP軟件對苔草rbcL、trnL-F兩個序列聯(lián)合分析,共得到58個單倍型,利用Network軟件構(gòu)建其單倍型網(wǎng)絡(luò)圖,見圖3。大多單倍型僅包含一個種質(zhì),其中矮生苔草和直穗苔草共享一個單倍型,同為H47。

    2.4" 遺傳距離分析

    基于rbcL、trnL-F以及合并序列分析的遺傳距離分別為0~0.027,0~0.053,0~0.042,平均遺傳距離分別為0.006,0.021,0.014。rbcL序列中鴨綠苔草、白頭山苔草、紫鱗苔草、紅穗苔草、華北苔草等的遺傳距離最小。trnL-F序列中異鱗苔草和灰脈苔草、藍(lán)葉苔草和‘藍(lán)辛格’、矮生苔草和直穗苔草、索蘭德苔草和三裂苔草遺傳距離最小?;趓bcL和trnL-F序列結(jié)合分析,矮生苔草和直穗苔草遺傳距離最小。

    2.5" 基于rbcL、trnL-F序列的苔草屬植物親緣關(guān)系分析

    使用MEGA軟件對苔草屬植物rbcL、trnL-F合并序列建樹,以具芒碎米莎草(C. microiria)為外類群,構(gòu)建NJ系統(tǒng)發(fā)育樹(圖4),59份苔草屬種質(zhì)可劃分為3類,第Ⅰ類包含寬葉苔草和毛緣寬葉苔草‘寶藏島’,二者葉片較寬,通常為披針形至寬條形,小穗均為兩性,4~10個,雄雌順序,苞片成佛焰苞狀,柱頭3個,屬于苔草亞屬寬葉苔草組(Sect. Siderostictae)。第Ⅱ類包含10種苔草,其中白穎苔草、中亞苔草、尖嘴苔草、砂鉆苔草、疏穗苔草、卵囊苔草屬于二柱苔草亞屬,棕櫚葉苔草、棕櫚葉苔草‘小矮人’、阿巴拉契亞苔草和多刺苔草分類地位暫不明確。二柱苔草亞屬小穗多數(shù),全部為兩性,無柄,雌雄順序或雄雌順序,聚集在小穗頂端,柱頭通常2個。其中尖嘴苔草、卵囊苔草、中亞苔草、白穎苔草小穗通常為雄雌順序,且密集為圓柱狀花序;疏穗苔草小穗為雌雄順序;砂鉆苔草為小穗單性,雌雄異株,花序長圓狀卵形。其余47個材料均屬于第Ⅲ類,包含苔草亞屬和分類地位不明確的種質(zhì),如藍(lán)葉苔草、全裂苔草等國外引種的種質(zhì)。其中形態(tài)分類學(xué)黑穗苔草組(Sect. Atratae)、指狀苔草組(Sect. Digitatae)下材料聚在一起,親緣關(guān)系更近。除未觀察到小穗的種質(zhì),其余材料小穗1個至多數(shù),單性或單性與兩性兼有,小穗單個或多個生于苞片腋內(nèi),具柄,少數(shù)柄很短或近無柄,柱頭通常3個。

    3" 討論

    3.1" rbcL、trnL-F序列特征

    在實(shí)際應(yīng)用中,常將cpDNA中進(jìn)化速率不一的序列結(jié)合來區(qū)分屬下不同物種[22]。本研究發(fā)現(xiàn),在兩個葉綠體序列中,trnL-F的突變總數(shù)、變異位點(diǎn)、單一突變位點(diǎn)、平均核苷酸多樣性等均高于rbcL。合并序列中247個變異位點(diǎn),109個單一突變位點(diǎn),138個簡約信息位點(diǎn),共58個單倍型,單倍型間聯(lián)系密切,通過研究各種單倍型間的關(guān)系能夠了解其種群遺傳變異及進(jìn)化關(guān)系。合并序列中變異位點(diǎn)、單倍型數(shù)目等指標(biāo)高于單一序列,證明多個序列拼接后具有豐富的遺傳信息,與Kim等[23-24]的研究結(jié)果一致。張仕蘭等[25]利用rbcL和trnL-F等5個葉綠體基因?qū)H馑懖荩邸?Bolboschoenoplectus mariqueter (Tang amp; F. T. Wang) Tatanov]及其近緣種的分類進(jìn)行研究,trnL-F的變異位點(diǎn)少于rbcL,變異程度較低。本研究中trnL-F的變異位點(diǎn)多于rbcL,變異程度較高,表現(xiàn)為較高的遺傳多樣性,這與rbcL和trnL-F序列在鳶尾(Iris tectorum)[26]中表現(xiàn)相似,因此,葉綠體基因序列的變異程度可能與物種密切相關(guān)。

    rbcL序列相對保守,種間變異小,因其強(qiáng)大的通用性被廣泛用于植物系統(tǒng)發(fā)育研究中,trnL-F為cpDNA的非編碼區(qū),其進(jìn)化速率快,片段兩端存在保守序列,適用于分類階元較低的類群間的系統(tǒng)發(fā)育研究[27-28]。Starr等[18]在使用引物rbcL-a_f和rbcL-a_r對苔草屬的rbcL序列進(jìn)行擴(kuò)增和測序,發(fā)現(xiàn)該序列的獲得率僅為18%,本研究采引物rbcL-1f和rbcL-724r對rbcL序列進(jìn)行擴(kuò)增和測序,序列的獲得率為95.8%,但序列變異率低,不能對苔草屬種質(zhì)進(jìn)行有效區(qū)分,與Clerc-Blain等[29]的研究結(jié)果一致。無論是單一序列還是多序列組合對苔草屬植物的鑒別率均不超過60%[18],本研究將rbcL和trnL-F片段組合,對物種的鑒別率較高,將進(jìn)化速率不一的片度進(jìn)行結(jié)合分析,可以提供更多的有效信息點(diǎn),從而提高研究結(jié)果的準(zhǔn)確性和可靠性[30]。

    3.2" rbcL、trnL-F單倍型多樣性

    rbcL和trnL-F合并序列將59份種質(zhì)定義為58個單倍型,較古老的單倍型通常位于聚類簇內(nèi)部,而稀有單倍型通常位于網(wǎng)絡(luò)分支末端,H1位于單倍型網(wǎng)絡(luò)分支的中心位置,則認(rèn)為該單倍型可能為原始單倍型。中性檢驗(yàn)的結(jié)果為負(fù)值,說明研究的群體可能在歷史上經(jīng)歷過突然擴(kuò)張或者受到了定向選擇,在分子水平導(dǎo)致進(jìn)化的原因可能是突變壓力或遺傳漂變[31-32]。葉綠體單倍型網(wǎng)絡(luò)圖中的各單倍體呈現(xiàn)的星狀輻射也能夠說明這一現(xiàn)象。矮生苔草和直穗苔草均屬于H47,二者分別屬于苔草亞屬中沼生苔草組、苔草組,其形態(tài)較相似,葉片較硬,扁平,小穗3~5個,上部2~3個為雄小穗,接近生,其余為雌小穗,但矮生苔草有粗壯而堅韌的匍匐根狀莖。矮生苔草分布于沿海地區(qū)的海邊沙地,采集于河北省秦皇島市黃金海岸,直穗苔草分布于沼澤地、河邊及河邊潮濕地,采集于河北省唐山市曹妃甸湖區(qū)域,二者的采集地環(huán)境氣候較為相似,地理距離近,可能導(dǎo)致其被定義為同一種葉綠體單倍型,有待結(jié)合核基因序列進(jìn)一步研究。

    3.3" 遺傳多樣性與親緣關(guān)系

    遺傳多樣性的高低可以通過遺傳距離的大小來判斷,遺傳距離越大說明親緣關(guān)系越遠(yuǎn),表明其遺傳多樣性越高,反之,遺傳距離越小,說明遺傳多樣性越低。59份苔草屬材料的遺傳距離為0~0.0421,平均遺傳距離為0.0142,親緣關(guān)系較近,遺傳多樣性較低,與白羊草(Bothriochloa ischaemum)[33]種質(zhì)遺傳多樣性的研究結(jié)果一致,可能與自身繁育系統(tǒng)的特性有關(guān)。無性繁殖、長期的自交或近交會降低種群的遺傳多樣性,根據(jù)現(xiàn)有報道來看,苔草屬植物的繁殖方式為營養(yǎng)繁殖和種子繁殖,由于苔草屬小穗頂生雄小穗、雄雌順序或雌雄順序的結(jié)構(gòu)特征以及開花習(xí)性,導(dǎo)致散粉后花粉幾乎全部落在雌花柱頭上,極大的增加了自交的概率[34-35]。此外,環(huán)境、氣候、人為干擾等外部因素也會間接影響植物的遺傳多樣性[36]。

    除分類不明確的材料外,本研究中苔草材料屬于苔草亞屬和二柱苔草亞屬,59份苔草屬種質(zhì)被劃分為3大類,聚類結(jié)果與形態(tài)學(xué)分類基本相符,與Villaverde等[37]對苔草屬重新分類的結(jié)果一致,將寬葉苔草組從苔草亞屬中獨(dú)立出來,成為寬葉亞屬,但目前我國還是基于形態(tài)學(xué)分類將苔草屬植物分為3個亞屬:苔草亞屬、二柱苔草亞屬、復(fù)序苔草亞屬。而基于ISSR[38]對河北省野生苔草屬植物進(jìn)行親緣關(guān)系分析,將53份供試材料分為6大類,并不能將亞屬級別的材料進(jìn)行區(qū)分,且與傳統(tǒng)形態(tài)學(xué)分類存在一定分歧,也與本研究基于葉綠體基因rbcL和trnL-F序列的聚類結(jié)果存在不同,推測可能是存在核質(zhì)沖突,但也可能是該種標(biāo)記的缺陷或人工讀帶誤差等導(dǎo)致的。從自展值來看,苔草亞屬中黑穗苔草組短鱗苔草、紫喙苔草與華北苔草、白頭山苔草、紫鱗苔草的自展值低于50%,有待結(jié)合更多序列來進(jìn)行種間的區(qū)分與鑒定,而部分苔草品種間的自展值高,可以進(jìn)一步區(qū)分同一種的不同品種,如芒髯苔草品種‘香草冰’ ‘冰舞’ ‘永恒光輝’和大島苔草品種‘金芯’ ‘銀邊’ ‘黃金發(fā)絲’,有望作為區(qū)分苔草品種的DNA條形碼。綜合來看,應(yīng)采用形態(tài)學(xué)和分子標(biāo)記相結(jié)合的方式對苔草屬植物進(jìn)行分類和鑒定。

    目前,苔草屬中應(yīng)用較多的為苔草亞屬和二柱苔草亞屬中的種質(zhì),已有不同種質(zhì)被開發(fā)利用,如青綠苔草、澇峪苔草、白穎苔草等。除此之外,還發(fā)現(xiàn)這兩個亞屬中披針葉苔草、早春苔草、矮叢苔草、疏穗苔草等植物葉長葉寬適中、株型圓整,觀賞效果好,而寬葉亞屬中寬葉苔草、‘寶藏島’葉片更寬大,株型整齊,整體較美觀。這些苔草適應(yīng)性強(qiáng),在華北地區(qū)可以露地越冬越夏,維護(hù)成本低,但目前應(yīng)用的較少,可以開發(fā)為園林地被或觀賞草。利用誘變育種或依據(jù)其親緣關(guān)系進(jìn)行雜交育種等多種方法,進(jìn)行品種改良和新品種選育,進(jìn)一步選育出具有良好觀賞性的種質(zhì)。

    4" 結(jié)論

    葉綠體基因rbcL、trnL-F序列以及二者合并序列分析均顯示苔草屬植物具有一定的遺傳多樣性。經(jīng)裁剪拼接得到的苔草的rbcL和trnL-F序列分別為716 bp和976 bp,上傳至NCBI獲得登錄號。rbcL序列較為穩(wěn)定,trnL-F序列具有更豐富的變異,兩條序列拼接長度為1692 bp,序列合并后定義了58種單倍型。59份種質(zhì)劃分為3類,與形態(tài)學(xué)分類有一定相關(guān)性,首次在河北省苔草屬中驗(yàn)證了寬葉亞屬的地位。trnL-F序列定義了56個單倍型,有望作為苔草屬植物鑒定的DNA條形碼。本研究為苔草屬植物在分類鑒定、系統(tǒng)發(fā)育及開發(fā)利用等方面提供了理論依據(jù)。

    參考文獻(xiàn)

    [1]" 張玉,王麗媛,王繼豐,等. 基于文獻(xiàn)計量學(xué)的苔草屬植物研究進(jìn)展[J]. 黑龍江科學(xué),2018,9(20):17-19

    [2]" TENG K,TENG W J,WEN H F,et al. PacBio single-molecule long-read sequencing shed new light on the complexity of the Carex breviculmis transcriptome[J]. BMC Genomics,2019,20(1):789

    [3]" 河北省植物智編輯委員會. 河北植物志III[M]. 石家莊:河北科學(xué)技術(shù)出版社,1991:408-439

    [4]" 張全剛,孫建華,郝愛,等. 淺談苔草在園林綠地中的應(yīng)用[J]. 農(nóng)業(yè)與技術(shù),2022,42(24):127-129

    [5]" 豆豆,李紅. 黃河源區(qū)退化草地褐鱗苔草種植技術(shù)及對土壤改良效果研究[J]. 甘肅農(nóng)業(yè)大學(xué)學(xué)報,2022,57(5):171-180

    [6]" 中國科學(xué)院中國植物委員會. 中國植物志.第十二卷[M]. 北京:科學(xué)出版社,2000:56-519

    [7]" 劉凌云,范希峰,滕珂,等. 苔草屬植物遺傳多樣性研究進(jìn)展[J]. 草地學(xué)報,2021,29(6):1148-1157

    [8]" 冷建紅,樓爐煥,魏琦,等. 薹草屬植物研究進(jìn)展及園林應(yīng)用前景研究[J]. 北方園藝,2011(6):196-201

    [9]" SHEKHOVTSOV S V,SHEKHOVTSOV I N,PELTEK S E. Phylogeny of Siberian species of Carex sect. Vesicariae based on nuclear and plastid markers[J]. Nordic Journal of Botany,2012,30(3):343-351

    [10]LVEILL-BOURRET,TIENNE,STARR J R,F(xiàn)ORD B A. Why are there so many sedges? Sumatroscirpeae,a missing piece in the evolutionary puzzle of the giant genus Carex (Cyperaceae)[J]. Molecular Phylogenetics and Evolution,2018,119:93-104

    [11]PERDEREAU A,KLAAS M,BARTH S,et al. Plastid genome sequencing reveals biogeographical structure and extensive population genetic variation in wild populations of Phalaris arundinacea L. in north-western Europe.[J].Global Change Biology Bioenergy,2017,9(1):46-56

    [12]姜棟柱. 基于葉綠體基因組和ITS序列的姜屬植物系統(tǒng)發(fā)育研究[D]. 荊州:長江大學(xué),2023:3-4

    [13]AMAR M H,MAGDY M,WANG L,et al. Peach chloroplast genome variation architecture and phylogenomic signatures of cpDNA introgression in Prunus[J]. Canadian Journal of Plant Science,2019,99(6):885-896

    [14]苗藝明,黃松殿,覃永華,等. 基于三個葉綠體基因序列的廣西盾翅藤遺傳多樣性分析[J/OL]. http://kns.cnki.net/kcms/detail/46.1068.S.20231020.1334.006.html,2023-10-20/2024-03-19

    [15]龍炳宏,蔣向輝,宋榮,等. DNA條形碼在黃精屬藥用植物鑒定與遺傳多樣性分析中的應(yīng)用[J]. 植物科學(xué)學(xué)報,2022,40(4):533-543

    [16]錢志瑤,饒海,曹旭林,等. 黔產(chǎn)寬葉纈草葉綠體trnL-F,matK基因遺傳關(guān)系的分析[J]. 中國實(shí)驗(yàn)方劑學(xué)雜志,2018,24(5):36-42

    [17]ODA J,F(xiàn)USE S,YAMASHITA J,et al. Phylogeny and Taxonomy of Carex (Cyperaceae) in Japan I. C. sect. Rarae[J]. Acta Phytotaxonomica et Geobotanica:APG,2019(2):69-85

    [18]STARR J R,NACZI R F C,CHOUINARD B N. Plant DNA barcodes and species resolution in sedges (Carex,Cyperaceae)[J]. Molecular Ecology Resources,2010,9(s1):151-163

    [19]YANO O,IKEDA H,HOSHINO T. Phylogeography of the Japanese common sedge,Carex conica complex (Cyperaceae),based on chloroplast DNA sequence data and chromosomal variation[J]. American Journal of Botany,2010,97(8):1365-1376

    [20]KRESS W J,ERICKSON D L,SHIU S H. A Two-Locus Global DNA Barcode for Land Plants:The Coding rbcL Gene Complements the Non-Coding trnH-psbA Spacer Region[J]. Plos One,2007,2(6):e508

    [21]TABERLET P,GIELLY L,PAUTOU G,et al. Universal primers for amplification of three non-coding regions of chloroplast DNA[J]. Plant Molecular Biology,1991,17(5):1105-1109

    [22]孫浩男,李蓉,王鑫,等. 基于三種分子標(biāo)記的栽培金雞菊遺傳多樣性及親緣關(guān)系分析[J]. 草地學(xué)報,2022,30(11):2922-2936

    [23]KIM Y,SHIN J,OH D R,et al. Comparative Analysis of Complete Chloroplast Genome Sequences and Insertion-Deletion (Indel) Polymorphisms to distinguish five vaccinium species[J]. Multidisciplinary Digital Publishing Institute,2020,11(9):927

    [24]王崇,焦春海,楊新筍,等. 基于葉綠體基因trnL-trnF、trnH-psbA和trnT-trnL序列的甘薯種質(zhì)遺傳多樣性分析[J]. 南方農(nóng)業(yè)學(xué)報,2021,52(6):1536-154

    [25]張仕蘭,王嶼岑,劉文亮. 海三棱藨草及其近緣種的DNA條形碼研究[J]. 廣西植物,2024,44(4):699-709

    [26]蔣喻林. 基于trnL-F、rbcL基因的四川省鳶尾屬植物系統(tǒng)分類進(jìn)化研究[D]. 成都:四川農(nóng)業(yè)大學(xué),2016:31-32

    [27]朱斌,胡利娟,田忠靜,等. 基于葉綠體基因組序列對50種菊科植物親緣關(guān)系分析[J]. 貴州師范大學(xué)學(xué)報(自然科學(xué)版),2021,39(2):15-25

    [28]王霞,劉艷,賈秀秀,等. 矩鐮莢苜蓿與近緣種花苜蓿的DNA條形碼篩選[J]. 草地學(xué)報,2023,31(8):2323-2333

    [29]CLERC-BLAIN J L,STARR J R,BULL R D,et al. A regional approach to plant DNA barcoding provides high species resolution of sedges (Carex and Kobresia,Cyperaceae) in the Canadian Arctic Archipelago[J]. Molecular Ecology Resources,2010,10(1):69-91

    [30]FAZEKAS A J,BURGESS K S,KESANAKURTI P R,et al. Multiple multilocus DNA barcodes from the plastid genome discriminate plant species equally well[J]. Plos One,2008,3(7):e2802

    [31]萬誠,李蕾嘉,陳勝,等. 貴州老鷹基于葉綠體DNA的遺傳多樣性分析[J]. 農(nóng)業(yè)與技術(shù),2022,42(15):22-25

    [32]趙字陳,趙玉娟,龔洵,等. 龍血樹柴胡的遺傳多樣性及其分布特征研究[J]. 植物科學(xué)學(xué)報,2022,40(1):1-10

    [33]李鈺瑩,董寬虎. 山西白羊草種質(zhì)資源遺傳多樣性的ISSR分析[J]. 草業(yè)學(xué)報,2014,23(1):217-222

    [34]文亞峰,韓文軍,吳順. 植物遺傳多樣性及其影響因素[J]. 中南林業(yè)科技大學(xué)學(xué)報,2010,30(12):80-87

    [35]林鵬,楊虎彪. 吊羅山薹草開花生物學(xué)特征及繁育系統(tǒng)研究[J]. 熱帶作物學(xué)報,2021,42(3):689-694

    [36]BOLTON P E,ROLLINS L A,BRAZILL-BOAST J,et al. Genetic diversity through time and space:diversity and demographic history from natural history specimens and serially sampled contemporary populations of the threatened Gouldian finch (Erythrura gouldiae)[J]. Conservation Genetics,2018,19(3):737-754

    [37]VILLAVERDE T,JIMNEZ-MEJAS P,LUCENO M,et al. A new classification of Carex (Cyperaceae) subgenera supported by a HybSeq backbone phylogenetic tree[J]. Botanical Journal of the Linnean Society,2020,194(2):141-163

    [38]周欣瑩,耿宇航,賈艷艷,等. 河北省野生苔草屬植物ISSR分子標(biāo)記開發(fā)及親緣關(guān)系分析[J]. 草地學(xué)報,2023,31(10):2949-2959

    (責(zé)任編輯" 彭露茜)

    日韩av在线免费看完整版不卡| 国产精品福利在线免费观看| 又爽又黄a免费视频| 精品久久久久久久末码| 国产成年人精品一区二区| 亚洲精品aⅴ在线观看| 一个人看视频在线观看www免费| 久热久热在线精品观看| 日韩欧美精品免费久久| 国产男人的电影天堂91| www.色视频.com| 国产高潮美女av| 熟女电影av网| 久久精品熟女亚洲av麻豆精品 | 黄片wwwwww| 精品熟女少妇av免费看| 日韩,欧美,国产一区二区三区 | 日本午夜av视频| 免费观看在线日韩| 亚洲性久久影院| 看黄色毛片网站| 乱系列少妇在线播放| 青春草亚洲视频在线观看| 亚洲自偷自拍三级| 亚洲综合色惰| 床上黄色一级片| 国产精品av视频在线免费观看| 天堂中文最新版在线下载 | 日日撸夜夜添| 精品久久久久久电影网 | 国产精品.久久久| 国产综合懂色| 建设人人有责人人尽责人人享有的 | 禁无遮挡网站| 国产一区二区在线av高清观看| 成人二区视频| 亚洲aⅴ乱码一区二区在线播放| 视频中文字幕在线观看| 在线观看一区二区三区| 免费看a级黄色片| 中文天堂在线官网| 欧美区成人在线视频| 久久久久久久久久成人| 少妇的逼水好多| av黄色大香蕉| 成人鲁丝片一二三区免费| 久久99热6这里只有精品| 一本一本综合久久| 亚洲国产精品成人综合色| 乱人视频在线观看| 久久精品综合一区二区三区| 成人国产麻豆网| 国内揄拍国产精品人妻在线| ponron亚洲| 国产高潮美女av| 亚洲欧美精品自产自拍| 22中文网久久字幕| 两性午夜刺激爽爽歪歪视频在线观看| 国产在视频线在精品| 久久精品国产99精品国产亚洲性色| 久久人人爽人人片av| 久久久色成人| 亚洲国产最新在线播放| 日本一二三区视频观看| 日日干狠狠操夜夜爽| АⅤ资源中文在线天堂| 亚洲一级一片aⅴ在线观看| 亚洲国产欧洲综合997久久,| 中文欧美无线码| 男的添女的下面高潮视频| 国产黄片美女视频| 赤兔流量卡办理| 色尼玛亚洲综合影院| 亚洲一区高清亚洲精品| a级毛色黄片| 亚洲国产精品成人久久小说| 色播亚洲综合网| 国产一区二区在线av高清观看| 国产精品久久久久久精品电影| 亚洲怡红院男人天堂| 日本wwww免费看| 日韩人妻高清精品专区| 国产精品麻豆人妻色哟哟久久 | 亚洲内射少妇av| 亚洲熟妇中文字幕五十中出| 久久亚洲精品不卡| 日日撸夜夜添| 日韩精品有码人妻一区| 久久热精品热| 夜夜看夜夜爽夜夜摸| 女人久久www免费人成看片 | 国产亚洲91精品色在线| 精华霜和精华液先用哪个| 一区二区三区高清视频在线| 精品国产一区二区三区久久久樱花 | 久久久久久久久久成人| 亚洲丝袜综合中文字幕| 又黄又爽又刺激的免费视频.| 亚洲不卡免费看| 国产在线一区二区三区精 | 99久久精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久久久丰满| 国产日韩欧美在线精品| 日韩成人伦理影院| 国产成人午夜福利电影在线观看| 成人美女网站在线观看视频| 国产 一区 欧美 日韩| av黄色大香蕉| 国产成人a区在线观看| av在线老鸭窝| 九九爱精品视频在线观看| 男人狂女人下面高潮的视频| 干丝袜人妻中文字幕| 久久99热6这里只有精品| 国产精品一区www在线观看| 晚上一个人看的免费电影| 99在线视频只有这里精品首页| 免费大片18禁| 婷婷六月久久综合丁香| 99热网站在线观看| 久久久亚洲精品成人影院| 亚洲中文字幕一区二区三区有码在线看| 在线天堂最新版资源| 成人亚洲欧美一区二区av| 日韩国内少妇激情av| 亚洲精华国产精华液的使用体验| 啦啦啦观看免费观看视频高清| 国产大屁股一区二区在线视频| 久热久热在线精品观看| 人妻系列 视频| 韩国高清视频一区二区三区| 午夜激情欧美在线| 亚洲成人久久爱视频| 男插女下体视频免费在线播放| 91精品伊人久久大香线蕉| 天美传媒精品一区二区| 久久久精品欧美日韩精品| 国产欧美另类精品又又久久亚洲欧美| 国产精品,欧美在线| 精品久久久久久电影网 | 日日干狠狠操夜夜爽| 最近视频中文字幕2019在线8| 日韩强制内射视频| 国产精品爽爽va在线观看网站| 变态另类丝袜制服| 欧美性猛交╳xxx乱大交人| 精品免费久久久久久久清纯| 久久99蜜桃精品久久| 国产精品99久久久久久久久| 欧美激情国产日韩精品一区| 一个人看视频在线观看www免费| 岛国毛片在线播放| 亚洲成色77777| 欧美激情久久久久久爽电影| 国产精品99久久久久久久久| 人妻制服诱惑在线中文字幕| 欧美不卡视频在线免费观看| 国语自产精品视频在线第100页| 国产黄片美女视频| 国产乱人视频| 日本五十路高清| 中文欧美无线码| 少妇丰满av| 亚洲一级一片aⅴ在线观看| 欧美日韩在线观看h| 国产精品av视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲熟妇中文字幕五十中出| 国产精品一区二区三区四区久久| 国产精品一区二区三区四区久久| 狂野欧美激情性xxxx在线观看| 九九热线精品视视频播放| 美女xxoo啪啪120秒动态图| 亚洲怡红院男人天堂| 亚洲第一区二区三区不卡| 边亲边吃奶的免费视频| 国产三级中文精品| 免费看a级黄色片| 亚洲精品成人久久久久久| 看非洲黑人一级黄片| 精品久久久久久久久亚洲| 99久久精品热视频| 亚洲国产欧美在线一区| 97在线视频观看| 狠狠狠狠99中文字幕| 男女视频在线观看网站免费| 禁无遮挡网站| 午夜激情欧美在线| 午夜福利成人在线免费观看| 国产在视频线在精品| 国产色婷婷99| 熟妇人妻久久中文字幕3abv| .国产精品久久| 精品一区二区三区人妻视频| 男人舔奶头视频| 久久99蜜桃精品久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | www日本黄色视频网| 免费观看a级毛片全部| 97超视频在线观看视频| 麻豆国产97在线/欧美| 国产真实乱freesex| 午夜视频国产福利| 亚洲成人精品中文字幕电影| 日韩强制内射视频| 中文乱码字字幕精品一区二区三区 | 亚洲欧美日韩高清专用| 日本一本二区三区精品| 亚洲中文字幕一区二区三区有码在线看| 黄片wwwwww| 欧美zozozo另类| 极品教师在线视频| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 99国产精品一区二区蜜桃av| 秋霞在线观看毛片| av在线老鸭窝| 国产av码专区亚洲av| 丰满人妻一区二区三区视频av| 成人无遮挡网站| 级片在线观看| 国产av码专区亚洲av| 国产精品一区二区三区四区久久| 中文乱码字字幕精品一区二区三区 | 亚洲欧美日韩卡通动漫| 亚洲国产成人一精品久久久| 日韩欧美在线乱码| 两个人的视频大全免费| 长腿黑丝高跟| 简卡轻食公司| 国产在线男女| 又粗又爽又猛毛片免费看| 国内揄拍国产精品人妻在线| 一级毛片电影观看 | 最近2019中文字幕mv第一页| 女的被弄到高潮叫床怎么办| 观看免费一级毛片| 亚洲国产日韩欧美精品在线观看| 老司机影院毛片| 欧美又色又爽又黄视频| 日韩制服骚丝袜av| 中国国产av一级| 插阴视频在线观看视频| 免费看日本二区| 国产成人精品一,二区| 国产精品乱码一区二三区的特点| 你懂的网址亚洲精品在线观看 | 亚洲va在线va天堂va国产| 最近的中文字幕免费完整| 中文字幕免费在线视频6| 亚洲av中文字字幕乱码综合| 国产精品一区二区在线观看99 | 大又大粗又爽又黄少妇毛片口| 国产毛片a区久久久久| 久久精品人妻少妇| 99久久精品热视频| 亚洲最大成人av| 亚洲国产最新在线播放| 久久精品国产亚洲网站| 真实男女啪啪啪动态图| 国产伦一二天堂av在线观看| 在线天堂最新版资源| 2021天堂中文幕一二区在线观| 日日撸夜夜添| 丝袜美腿在线中文| 一级黄片播放器| 久久精品国产鲁丝片午夜精品| 亚洲人成网站在线观看播放| 22中文网久久字幕| av国产免费在线观看| 长腿黑丝高跟| 色播亚洲综合网| 99热6这里只有精品| 在线播放国产精品三级| 少妇的逼水好多| 国产av在哪里看| 男女啪啪激烈高潮av片| 少妇被粗大猛烈的视频| 国产成人精品婷婷| 天天躁日日操中文字幕| 韩国高清视频一区二区三区| 99热精品在线国产| 国产白丝娇喘喷水9色精品| 国产极品天堂在线| av又黄又爽大尺度在线免费看 | 亚洲av电影不卡..在线观看| 国产午夜精品论理片| 国产高潮美女av| 国产精品久久久久久精品电影| 国产精品女同一区二区软件| 亚洲欧美一区二区三区国产| 亚洲天堂国产精品一区在线| 亚洲性久久影院| av在线蜜桃| 久久99热这里只有精品18| av在线播放精品| 干丝袜人妻中文字幕| 久久精品熟女亚洲av麻豆精品 | 日本黄色片子视频| 日韩一区二区视频免费看| .国产精品久久| 亚洲精品日韩av片在线观看| 一边亲一边摸免费视频| 久久久久网色| 亚洲久久久久久中文字幕| 日韩欧美在线乱码| 天堂中文最新版在线下载 | 国产精品久久久久久精品电影小说 | 久久精品综合一区二区三区| 亚洲国产最新在线播放| 亚洲精品456在线播放app| 欧美一级a爱片免费观看看| 国产精品国产高清国产av| 久久精品国产自在天天线| 综合色av麻豆| 麻豆成人午夜福利视频| 国产色爽女视频免费观看| 亚洲色图av天堂| 亚洲真实伦在线观看| 亚洲av成人av| av国产久精品久网站免费入址| 超碰97精品在线观看| 国产午夜精品一二区理论片| 寂寞人妻少妇视频99o| 久久久久久久久久黄片| 一个人观看的视频www高清免费观看| 麻豆乱淫一区二区| 精品国内亚洲2022精品成人| 日本-黄色视频高清免费观看| 色噜噜av男人的天堂激情| 观看美女的网站| 最近的中文字幕免费完整| 亚洲三级黄色毛片| av.在线天堂| 99热网站在线观看| 久久韩国三级中文字幕| 午夜精品国产一区二区电影 | 人人妻人人澡欧美一区二区| 成年女人永久免费观看视频| 99久国产av精品国产电影| 久久精品国产亚洲av天美| 九草在线视频观看| 日本色播在线视频| 99九九线精品视频在线观看视频| 日韩 亚洲 欧美在线| 国产探花在线观看一区二区| 亚洲国产高清在线一区二区三| 热99在线观看视频| av视频在线观看入口| 久久久国产成人精品二区| 亚洲人成网站在线观看播放| 永久免费av网站大全| 丰满人妻一区二区三区视频av| 免费搜索国产男女视频| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 波多野结衣高清无吗| 有码 亚洲区| 亚洲国产色片| 精品久久久噜噜| 夜夜看夜夜爽夜夜摸| 高清日韩中文字幕在线| 美女内射精品一级片tv| 欧美极品一区二区三区四区| 国产精品乱码一区二三区的特点| 国产精品福利在线免费观看| 国产精品女同一区二区软件| 国产精品国产三级国产专区5o | 亚洲伊人久久精品综合 | 国产单亲对白刺激| 国产精品99久久久久久久久| 亚洲伊人久久精品综合 | 亚洲欧美精品综合久久99| 久久久久国产网址| 国产不卡一卡二| 免费观看的影片在线观看| 国产高清三级在线| 人妻夜夜爽99麻豆av| 99在线人妻在线中文字幕| 简卡轻食公司| 少妇高潮的动态图| 搡女人真爽免费视频火全软件| 好男人视频免费观看在线| 成年av动漫网址| 免费av不卡在线播放| 又粗又硬又长又爽又黄的视频| 1024手机看黄色片| 国产在线男女| 久久热精品热| 亚洲无线观看免费| 亚洲欧美日韩无卡精品| 美女xxoo啪啪120秒动态图| 国产精品一二三区在线看| 欧美人与善性xxx| av国产免费在线观看| 精品久久久久久久末码| 国产一区有黄有色的免费视频 | 中文字幕人妻熟人妻熟丝袜美| 美女脱内裤让男人舔精品视频| 爱豆传媒免费全集在线观看| 男人和女人高潮做爰伦理| 免费人成在线观看视频色| 日韩av不卡免费在线播放| 欧美性猛交黑人性爽| 色视频www国产| 女人久久www免费人成看片 | 美女高潮的动态| 精品少妇黑人巨大在线播放 | 欧美3d第一页| 欧美日本视频| 亚洲成人中文字幕在线播放| 国产免费视频播放在线视频 | 韩国av在线不卡| 久久久久久久久久成人| av国产免费在线观看| 亚洲婷婷狠狠爱综合网| 国产精品麻豆人妻色哟哟久久 | 国产一级毛片七仙女欲春2| 亚洲国产欧美人成| 超碰97精品在线观看| 日本-黄色视频高清免费观看| 欧美zozozo另类| av在线老鸭窝| 日本五十路高清| 一二三四中文在线观看免费高清| 日韩视频在线欧美| 我的老师免费观看完整版| 九九热线精品视视频播放| 又爽又黄无遮挡网站| 网址你懂的国产日韩在线| 美女黄网站色视频| 精品无人区乱码1区二区| www.av在线官网国产| 国产精品久久久久久精品电影| 亚洲欧美清纯卡通| 成人美女网站在线观看视频| 日日啪夜夜撸| 天天躁日日操中文字幕| 亚洲av成人精品一区久久| 亚洲精品日韩在线中文字幕| 欧美性猛交╳xxx乱大交人| 如何舔出高潮| 最近中文字幕2019免费版| 国产国拍精品亚洲av在线观看| 人妻系列 视频| 一个人看的www免费观看视频| 精品午夜福利在线看| 日韩一区二区三区影片| 亚洲经典国产精华液单| 日韩欧美三级三区| 热99re8久久精品国产| 寂寞人妻少妇视频99o| 日韩三级伦理在线观看| 亚洲欧美中文字幕日韩二区| 国产伦一二天堂av在线观看| 亚洲av成人精品一二三区| 国产人妻一区二区三区在| 欧美+日韩+精品| 免费av毛片视频| 亚洲成人久久爱视频| 两性午夜刺激爽爽歪歪视频在线观看| av免费在线看不卡| 非洲黑人性xxxx精品又粗又长| 天天一区二区日本电影三级| 久久久久久久久久黄片| 少妇人妻精品综合一区二区| 插逼视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲人成网站在线播| 国产成人精品久久久久久| 美女xxoo啪啪120秒动态图| 在线播放国产精品三级| 午夜日本视频在线| 亚洲精品影视一区二区三区av| 国产探花在线观看一区二区| 精品人妻视频免费看| 最近手机中文字幕大全| 久久亚洲国产成人精品v| www.色视频.com| 美女高潮的动态| 亚洲av日韩在线播放| 久久6这里有精品| 国产在线男女| 99热网站在线观看| 建设人人有责人人尽责人人享有的 | 麻豆国产97在线/欧美| 中文字幕亚洲精品专区| 国产乱人视频| 18禁在线无遮挡免费观看视频| 赤兔流量卡办理| 七月丁香在线播放| 少妇人妻精品综合一区二区| 白带黄色成豆腐渣| 嫩草影院精品99| 亚洲人成网站高清观看| av免费观看日本| 欧美xxxx性猛交bbbb| 欧美区成人在线视频| 99九九线精品视频在线观看视频| 观看免费一级毛片| 五月伊人婷婷丁香| 精品一区二区三区视频在线| 在线a可以看的网站| 亚洲av一区综合| 精品久久久噜噜| 热99在线观看视频| 日本黄色片子视频| 午夜福利网站1000一区二区三区| 日韩av在线大香蕉| 欧美日本视频| 久久精品影院6| 国产大屁股一区二区在线视频| 午夜a级毛片| 欧美丝袜亚洲另类| 国产老妇女一区| av专区在线播放| 国产一区二区在线观看日韩| 禁无遮挡网站| av福利片在线观看| 少妇裸体淫交视频免费看高清| 联通29元200g的流量卡| 国产一区二区在线观看日韩| 国产中年淑女户外野战色| 亚洲成色77777| 青春草亚洲视频在线观看| 超碰97精品在线观看| 小说图片视频综合网站| 青青草视频在线视频观看| 欧美成人午夜免费资源| 国内精品美女久久久久久| 毛片一级片免费看久久久久| 国产精品人妻久久久久久| 中文精品一卡2卡3卡4更新| 国产久久久一区二区三区| 亚洲欧美精品专区久久| 能在线免费观看的黄片| 国产精品久久久久久久电影| 色播亚洲综合网| 欧美日韩国产亚洲二区| 国产精品三级大全| 国产成人a区在线观看| 亚洲欧美成人综合另类久久久 | 亚洲第一区二区三区不卡| 国产成人a∨麻豆精品| 国产免费男女视频| 一级二级三级毛片免费看| av专区在线播放| 国产精品一区二区三区四区久久| 国产精品永久免费网站| eeuss影院久久| 99热这里只有是精品在线观看| 午夜福利高清视频| 插阴视频在线观看视频| 亚洲av二区三区四区| 两个人视频免费观看高清| av在线亚洲专区| 九九热线精品视视频播放| 午夜精品一区二区三区免费看| 久久精品夜夜夜夜夜久久蜜豆| 婷婷色av中文字幕| 日产精品乱码卡一卡2卡三| 最近2019中文字幕mv第一页| 三级国产精品欧美在线观看| 午夜福利视频1000在线观看| 国产成人精品久久久久久| 能在线免费观看的黄片| 国产女主播在线喷水免费视频网站 | 日韩欧美 国产精品| 久久久国产成人免费| 午夜福利在线观看免费完整高清在| 欧美xxxx黑人xx丫x性爽| 精品酒店卫生间| 亚洲国产精品成人综合色| 免费看光身美女| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验| 久久久久久久午夜电影| 亚洲丝袜综合中文字幕| 床上黄色一级片| 免费人成在线观看视频色| 有码 亚洲区| 看片在线看免费视频| 国产精品三级大全| 精品久久久久久久人妻蜜臀av| 久久久久久国产a免费观看| 看片在线看免费视频| 丰满乱子伦码专区| 国产又色又爽无遮挡免| 中文字幕熟女人妻在线| 欧美日韩国产亚洲二区| 免费大片18禁| 国产成人freesex在线| 两个人视频免费观看高清| 国语自产精品视频在线第100页| 天天躁日日操中文字幕| .国产精品久久| 国语对白做爰xxxⅹ性视频网站| 小说图片视频综合网站| 久久久久九九精品影院| 亚洲欧美精品自产自拍| 熟女人妻精品中文字幕| 亚洲精品久久久久久婷婷小说 | 精品久久久久久久人妻蜜臀av| 舔av片在线| 免费看光身美女| 性插视频无遮挡在线免费观看| 最近中文字幕2019免费版| 亚洲国产最新在线播放| 夜夜看夜夜爽夜夜摸| 久久久色成人| 国产精品综合久久久久久久免费| 亚洲婷婷狠狠爱综合网| 不卡视频在线观看欧美|