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

    利用EST-SSR、IT-ISJ分子標記研究四倍體與八倍體柳枝稷的遺傳多樣性

    2016-11-12 06:14:26鐘鳴周思凡張新全黃秀嚴海東張愛玲羅茜黃琳凱
    草業(yè)學報 2016年10期
    關(guān)鍵詞:四倍體柳枝多態(tài)性

    鐘鳴,周思凡,張新全,黃秀,嚴海東,張愛玲,羅茜,黃琳凱

    (四川農(nóng)業(yè)大學動物科技學院草業(yè)科學系,四川 成都 611130)

    ?

    利用EST-SSR、IT-ISJ分子標記研究四倍體與八倍體柳枝稷的遺傳多樣性

    鐘鳴**,周思凡**,張新全,黃秀,嚴海東,張愛玲,羅茜,黃琳凱*

    (四川農(nóng)業(yè)大學動物科技學院草業(yè)科學系,四川 成都 611130)

    柳枝稷是重要的能源植物之一,為研究四倍體Alamo與通過秋水仙素加倍后的八倍體之間的遺傳差異,本研究以65份四倍體Alamo與其誘變的八倍體為研究材料,通過EST-SSR與IT-ISJ兩種分子標記進行遺傳分析。結(jié)果表明,EST-SSR引物共擴增出245條帶,多態(tài)性比率為99.59%,觀測等位基因Na為1.9878,有效等位基因Ne為1.1718,基因多樣性指數(shù)H為0.1308,Shannon信息指數(shù)I為0.2388,其中四倍體Alamo的各遺傳參數(shù)(226條,88.93%,1.8816,1.1823,0.1349,0.2398)多高于八倍體Alamo(217條,92.24%,1.9184,1.1652,0.1228,0.2234)。IT-ISJ引物共擴增出165條條帶,多態(tài)性比率為96.36%,觀測等位基因Na為1.9878,有效等位基因Ne為1.5777,基因多樣性指數(shù)H為0.3354,Shannon信息指數(shù)I為0.4999,四倍體Alamo各項指數(shù)(145條,94.55%,1.9333,1.6366,0.3602,0.5280)均高于八倍體Alamo(139條,84.24%,1.7515,1.4377,0.2573,0.3864)。通過UPGMA分析表明,65個柳枝稷材料中四倍體與八倍體沒有明顯區(qū)分開,表明其沒有產(chǎn)生明顯的遺傳分化。通過AMOVA分析表明,EST-SSR和IT-ISJ的遺傳變異分別有94.78%和80.76%發(fā)生在倍性內(nèi),有5.22%和19.24%發(fā)生在倍性間,表明柳枝稷倍性間沒有明顯的差異。

    柳枝稷;EST-SSR;IT-ISJ;Alamo;分子標記;遺傳多樣性

    柳枝稷(Panicumvirgatum)是一種禾本科多年生C4植物,具有地域分布廣,基因型豐富,纖維素含量高等優(yōu)勢[1]。在美國,柳枝稷已被作為能源模式物種被廣泛研究;在中國,柳枝稷的開發(fā)利用也逐步得到重視和發(fā)展。柳枝稷作為異花授粉植物經(jīng)歷了漫長的進化演變后,其倍性材料主要有四倍體、六倍體和八倍體。其中,多倍體由于基因的劑量和互作效應(yīng)往往會導(dǎo)致植物產(chǎn)生一系列生理生化功能的變化[2],在柳枝稷種質(zhì)資源中八倍體較四倍體有新的表型和特性,比如器官增大,抗旱性、抗病蟲性、木質(zhì)纖維素構(gòu)成變化等。柳枝稷豐富的倍性引起的表型變異使得柳枝稷遺傳多樣性豐富,遺傳差異較大[3]。

    物種的遺傳多樣性是種質(zhì)資源開發(fā)利用與保護的基礎(chǔ)。分子標記是目前研究植物遺傳多樣性最準確、最主流的檢測方法之一,它不僅可以直接反映基因組DNA上的變化,而且也能反映基因序列上的單個核苷酸變異[4]。目前,分子標記已廣泛應(yīng)用于植物遺傳圖譜構(gòu)建、目的基因定位以及種質(zhì)資源研究等[5],如隨機擴增多態(tài)DNA(random amplification polymorphic DNA,RAPD)構(gòu)建苜蓿(Medicago)的遺傳連鎖圖譜[6],微衛(wèi)星DNA標記(simple sequence repeats,SSR)定位斯卑爾脫小麥(Triticumspelta)YrSp基因[7],隨機擴增多態(tài)DNA標記(random amplification polymorphic DNA,RAPD)研究高羊茅(Festucaarundinacea)種質(zhì)資源遺傳多樣性[8]。本試驗利用EST-SSR(expressed sequence tag-simple sequence repeat)和IT-ISJ(intron targeted intron-exon splice junction)分子標記進行四倍體Alamo與八倍體Alamo的遺傳多樣性分析。在眾多標記中,簡單重復(fù)序列(simple sequence repeat,SSR)具有多態(tài)性高、多等位性、共顯性、重復(fù)性高、易用于PCR檢測等優(yōu)點[9]。SSR可分為基因組SSR(g-SSR)和表達序列標簽SSR(EST-SSR)等,而EST-SSR與其他相比較經(jīng)濟,利用相對充分[10]。目前,EST-SSR用于柳枝稷種質(zhì)資源研究已有報道[11-13]。Weining等[14]根據(jù)內(nèi)含子-外顯子拼接位點的保守序列設(shè)計擴增內(nèi)含子或外顯子的ISJ標記引物,并應(yīng)用于大麥(Hordeumvulgare)和小麥(Triticumaestivum),發(fā)現(xiàn)ISJ單引物PCR擴增效果差,而ISJ引物與α-淀粉酶基因保守序列引物組合的聚合酶鏈反應(yīng)(polymerase chain reaction,PCR)擴增效果好。西南大學的張正圣實驗室根據(jù)植物基因內(nèi)含子剪接位點設(shè)計合成的IT-ISJ分子標記,其正反引物可進行自由的組合,成本低,重現(xiàn)性好,已成功用于陸地棉(Gosspiumhirsutum)重組近交系群體構(gòu)建遺傳圖譜[15],也可用于比較基因組學的研究。目前在棉花[16]、花生(Arachishypogaea)[17]等物種上都得到了運用,但尚未應(yīng)用于能源植物中。兩種分子標記的結(jié)合,不僅可以驗證兩種分子標記在柳枝稷多態(tài)性上的可行性,還可對比驗證柳枝稷不同倍性間的遺傳多樣性[18-20]。

    本研究采用EST-SSR和IT-ISJ對柳枝稷Alamo四倍體及通過秋水仙素離體誘導(dǎo)獲得的柳枝稷八倍體Alamo進行了遺傳多樣性分析,以了解這些多倍體材料的遺傳背景和DNA遺傳多樣性分析,有助于全面深入了解柳枝稷四倍體Alamo與八倍體Alamo的遺傳背景信息,為柳枝稷品種的定向選育和其他禾本科多倍化遺傳變異研究提供一定參考。

    1 材料與方法

    1.1材料

    本試驗采用65份柳枝稷品種Alamo不同倍性的材料(表1),其中28份為柳枝稷四倍體Alamo材料(編號1-22為隨機選取的四倍體單株,編號23-28為從四倍體中選擇的出愈率和轉(zhuǎn)化率高的種子),37份為柳枝稷八倍體Alamo材料(其中編號29為四倍體Alamo材料經(jīng)秋水仙素加倍得來[3],編號30-65為編號29自交產(chǎn)生的子一代)。實驗材料于2012年4月種植于四川農(nóng)業(yè)大學雅安草學試驗基地,其新鮮幼嫩葉片于2014年7月采摘,放入裝有足量變色硅膠的收集袋中干燥保存?zhèn)溆谩?/p>

    表1 柳枝稷供試材料

    1.2基因組DNA的提取和濃度測定

    每份材料選取長勢良好的一個單株,隨機取5片幼嫩葉片混合,采用植物基因組DNA提取試劑盒(北京天根生化科技有限公司)進行基因組DNA的提取,對提取出的DNA采用2.0%瓊脂糖凝膠電泳和超微量分光光度儀檢測其純度和濃度,經(jīng)提純合格后稀釋至20 ng/μL的工作濃度,-20 ℃冷藏備用。

    1.3EST-SSR和IT-ISJ標記擴增

    選取柳枝稷四倍體4Alamo、HR6和八倍體8Alamo、8Alamo-17四個供試材料作為DNA模板,分別對24對禾本科通用EST-SSR引物[11]及參考鄭靚等[15]實驗中11條上游IT-ISJ引物和34條下游IT-ISJ引物隨機組合產(chǎn)生的40對IT-ISJ引物進行篩選,從中選出條帶清晰,多態(tài)性好的引物進行后續(xù)擴增。所篩選出的EST-SSR引物和IT-ISJ引物詳細內(nèi)容見表2。

    表2 EST-SSR和IT-ISJ引物信息

    注:a表示EST-SSR引物,b表示IT-ISJ引物。

    Note: a means EST-SSR primers, while b means IT-ISJ primers.

    EST-SSR的擴增反應(yīng)體系為15 μL:DNA 1 μL (20 ng/μL),上下游引物各1 μL(10 pmol/μL),Taq DNA酶0.5 μL(2.5 U/μL),Mix混合液7.5 μL(包含Mg2+、dNTPs、10×PCR buffer),ddH2O 3.5 μL。EST-SSR擴增反應(yīng)程序為:95 ℃預(yù)變性10 min;9 ℃變性30 s,55~60 ℃退火30 s,72 ℃變性30 s,共35個循環(huán);最后72 ℃延伸7 min,4 ℃保存。

    IT-ISJ 擴增采用15 μL反應(yīng)體系:DNA 2 μL (20 ng/μL),上下游引物0.75 μL(10 pmol/μL),Taq DNA酶0.5 μL(2.5 U/μL),Mix 混合液7.5 μL(包含Mg2+、dNTPs、10×PCR buffer),ddH2O 3.5 μL。IT-ISJ擴增反應(yīng)程序為:94 ℃預(yù)變性5 min;94 ℃變性45 s,50~55 ℃退火45 s,72 ℃變性1 min,共40個循環(huán);最后72 ℃延伸10 min,4 ℃保存。

    EST-SSR和IT-ISJ擴增產(chǎn)物均用6%變性聚丙烯酰胺凝膠電泳進行檢測,400 V電壓下電泳90 min,以DL2000 DNA Marker作為標準分子量對照,電泳結(jié)束后沖洗膠板,1 g/L AgNO3染色14 min,顯影后拍照保存(圖1,圖2)。

    圖1 EST-SSR引物51擴增圖譜Fig.1 Result of PCR amplification of No.51 EST-SSR primer

    圖2 IT-ISJ引物5F5R擴增圖譜Fig.2 Result of PCR amplification of 5F5R IT-ISJ primer

    1.4擴增譜帶統(tǒng)計分析

    對EST-SSR和IT-ISJ擴增產(chǎn)物的電泳結(jié)果分別比對與校正后,根據(jù)標準分子質(zhì)量標記,對擴增條帶在電泳圖譜上的位置,從大到小依次編號和統(tǒng)計,對有擴增條帶(包括強帶或清晰弱帶)的記為“1”,無帶的記為“0”。將建立的EST-SSR和IT-ISJ 標記0,1原始矩陣采用Excel 2010和POPGENE V1.32軟件[21]分別計算65個供試柳枝稷材料的各個遺傳多樣性參數(shù),包括總擴增條帶數(shù)(TNB)、多態(tài)性條帶數(shù)(NPB)、多態(tài)性條帶百分比(PP)、基因多樣性指數(shù)(H)、Shannon指數(shù)(I)、觀測等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、群體總基因多樣性(Ht)、群體內(nèi)基因多樣性(Hs)。利用UPGMA法(unweighted pair-group method with arithmetic mean,算術(shù)平均非加權(quán)配組法)進行聚類分析。并通過AMOVE 1.55軟件對柳枝稷四倍體Alamo和八倍體Alamo間與內(nèi)的遺傳距離進行分子變異方差分析(analysis of molecular variance)。

    2 結(jié)果與分析

    2.1EST-SSR和IT-ISJ擴增多樣性分析

    從18對EST-SSR引物和40對IT-ISJ引物中分別篩選得到12和15對具有較好擴增結(jié)果的引物,其擴增條帶大小分別在100~250 bp和100~500 bp之間。EST-SSR和IT-ISJ標記的多態(tài)性水平和信息量如表3所示。從中可以看出EST-SSR和IT-ISJ分別擴增出條帶245和165條,多態(tài)性比率(PPB)分別為99.59%和96.36%。兩種標記的多態(tài)性比例都比較高,且差異不大。但兩種標記效率的比較分析可見,EST-SSR的有效多元比率(EMR)和標記指數(shù)(MI)均高于IT-ISJ,且高達一倍左右。表明EST-SSR標記具有較強的品種辨別能力和引物效率。

    對EST-SSR和IT-ISJ分子標記的遺傳參數(shù)計算結(jié)果顯示觀測等位基因Na分別為1.9878和1.9455,有效等位基因Ne分別為1.1718和1.5777,基因多樣性指數(shù)H分別為0.1308和0.3354,Shannon信息指數(shù)I分別為0.2388和0.4999(表4,表5)。兩分子標記測得的觀測等位基因結(jié)果基本一致,但IT-ISJ測得的有效等位基因、基因多樣性指數(shù)與Shannon指數(shù)分別都高于EST-SSR。除EST-SSR的四倍體Alamo的Na低于八倍體Alamo,EST-SSR和IT-ISJ標記的四倍體Alamo相應(yīng)各遺傳參數(shù)(Na、Ne、H和I)均高于八倍體Alamo。

    表3 EST-SSR和IT-ISJ標記效率比較分析

    表4 2種倍性柳枝稷EST-SSR分析結(jié)果

    表5 2種倍性柳枝稷IT-ISJ分析結(jié)果

    2.265份柳枝稷不同倍性Alamo材料的資源遺傳聚類分析

    利用UPGMA法對供試柳枝稷材料EST-SSR、IT-ISJ以及EST-SSR與IT-ISJ數(shù)據(jù)進行了聚類分析(圖3,圖4,圖5)。EST-SSR數(shù)據(jù)表明四倍體Alamo與八倍體Alamo沒有明顯區(qū)分開,總體來看,65份供試材料主要可分為三類。類群1為2個四倍體Alamo和22個八倍體Alamo,類群2為9個四倍體Alamo和15個八倍體Alamo,類群3為17個四倍體Alamo。根據(jù)IT-ISJ聚類圖顯示,四倍體Alamo與八倍體Alamo大體分開,65份材料可分為兩大類。其中第1類為6個四倍體Alamo和36個八倍體Alamo,第2類為22個四倍體Alamo和1個八倍體Alamo。將兩種分子標記聯(lián)合分析可得出,65個材料可分為三大類,且四倍體Alamo與八倍體Alamo沒有完全分開。其中第1類包含6個四倍體Alamo和35個八倍體Alamo,第2類包含12個四倍體Alamo和2個八倍體Alamo,第3類包含10個四倍體Alamo。

    圖3 EST-SSR標記對65個柳枝稷植株親緣關(guān)系聚類圖Fig.3 Dendrogram of the relationship UPGMA cluster of 65 switchgrass based on EST-SSR markers

    圖4 IT-ISJ標記對65個柳枝稷植株親緣關(guān)系聚類圖Fig.4 Dendrogram of the relationship unweighted pair-group method with arithmetic means (UPGMA) cluster of 65 switchgrass based on IT-ISJ markers

    圖5 EST-SSR和IT-ISJ標記對65個柳枝稷植株親緣關(guān)系聚類圖Fig.5 Dendrogram of the relationship UPGMA cluster of 65 switchgrass based on EST-SSR and IT-ISJ markers

    2.365份柳枝稷不同倍性Alamo材料遺傳分化分析

    對柳枝稷兩個倍性Alamo材料建立的EST-SSR和IT-ISJ原始矩陣,運用AMOVA分析了兩個倍性群體的遺傳結(jié)構(gòu)。結(jié)果表明:兩個標記總的遺傳變異中有95%(EST-SSR)和81%(IT-ISJ)發(fā)生在柳枝稷Alamo兩倍性之內(nèi),而Alamo兩倍性之間的遺傳變異較低,只有5%(EST-SSR)和19%(IT-ISJ),兩個標記Alamo倍性群體之間和之內(nèi)的差異均極顯著(P<0.001)(表6和表7)。此外,兩標記得到的Alamo倍性間的遺傳分化系數(shù)Фst值0.052(EST-SSR)和0.192(IT-ISJ),也說明遺傳結(jié)構(gòu)的變異主要存在于倍性之內(nèi)。

    表6 EST-SSR 柳枝稷群體分子變異方差分析

    表7 IT-ISJ 柳枝稷群體分子變異方差分析

    3 討論

    目前,關(guān)于其他植物多倍體的遺傳多樣性研究有較多報道,所得到結(jié)果基本呈兩種分化,一種為倍性材料DNA水平間沒有明顯的差異[22-23];而另一種為在倍性材料間發(fā)現(xiàn)了明顯的基因差異[24-25]。根據(jù)柳枝稷材料間遺傳關(guān)系的聚類圖來看(圖3,4,5),各材料的聚類分布與其倍性間沒有存在明顯的差異。本研究中EST-SSR和IT-ISJ標記在柳枝稷Alamo兩倍性間,四倍體Alamo相應(yīng)各遺傳參數(shù)(Na、Ne、H和I)多數(shù)略高于八倍體Alamo(表3,4,5),根據(jù)相關(guān)柳枝稷倍性的研究表明,柳枝稷兩個生態(tài)型之間遺傳差異不大,但四倍體具有較高水平的多態(tài)性,且具有異位顯性基因和互補的基因作用,通常將四倍體作為雜種優(yōu)勢群體[26]。Brunken等[27]的研究表明,整倍體倍性變異和非整倍體倍性變異能夠增加位點的等位基因數(shù),但這種情況在等位頻率極低的條件下不會發(fā)生或者可能發(fā)生,這與群體內(nèi)染色體數(shù)目變異的狀況一致。本試驗所得結(jié)果與之基本相似,分析導(dǎo)致這一實驗結(jié)果的原因:其一,可能由于本實驗中八倍體Alamo是由四倍體Alamo經(jīng)過秋水仙素人工加倍而來,所以基因組總體差異不大,但是在秋水仙素誘變過程中的影響可能導(dǎo)致兩群體間產(chǎn)生差異[28];其二,通過基因加倍后,八倍體在某些基因位點經(jīng)過基因重復(fù)后某些特異性條帶變成了非特異性條帶,基因更趨于穩(wěn)定,從而減少了八倍體的多態(tài)性。柳枝稷有低地型和高地型兩種生態(tài)型,其中低地型為四倍體,植株高大生物產(chǎn)量高;高地型多為八倍體,抗寒能力強[29]。采用雜交育種方式聚合兩種生態(tài)型一直是育種家致力于解決的問題,但兩種倍性間雜交相對較難,低地型八倍體的創(chuàng)制及遺傳多樣性的信息可為兩種生態(tài)型雜交育種提供寶貴的物質(zhì)和數(shù)據(jù)支撐。

    本試驗采用的兩個分子標記分別從12對EST-SSR和15對IT-ISJ引物中擴增出245和165條清晰可辨條帶(TNB),分子標記檢測效率的評價指標很多,目前廣泛應(yīng)用的有MI、PIC、DI、EMR等[30]。本研究對EST-SSR和IT-ISJ的PPB、PIC、EMR、MI指標進行了比較研究。結(jié)果表明,對于EST-SSR和IT-ISJ來說,都能產(chǎn)生各自有效的多態(tài)性條帶,但多態(tài)性水平及檢測效率各不相同。EST-SSR的基本重復(fù)單元是由幾十個核苷酸組成的,重復(fù)次數(shù)一般為10~50,同一類微衛(wèi)星DNA可分布在基因組的不同位置上[19],基因組中EST-SSR多態(tài)性較高,并且等位基因數(shù)目多,呈共顯性遺傳,所以可提供的多態(tài)性信息量較高,結(jié)果穩(wěn)定可靠[20],因此,被廣泛應(yīng)用于遺傳多樣性研究。隨著EST-SSR標記位點在圖譜上的豐富,它以高度的多態(tài)性信息量和PCR技術(shù)的方便性使其表現(xiàn)出更大的優(yōu)勢。IT-ISJ作為一種新興的分子標記方式,在技術(shù)上可能略顯不成熟,還需要進一步完善。IT-ISJ作為一種新型的DNA分子標記,迄今尚未在能源草上應(yīng)用,本研究探討了其在柳枝稷上的應(yīng)用效果?,F(xiàn)今報道的植物遺傳多樣性相關(guān)研究多是單一分子標記鑒定的結(jié)果,本研究運用EST-SSR和IT-ISJ標記技術(shù)同時對柳枝稷遺傳多樣性進行了分析,不僅驗證了兩種分子標記的可行性,還對比了柳枝稷不同倍性間的遺傳多樣性差異。試驗表明,EST-SSR和IT-ISJ標記對柳枝稷同一群體多樣性的評價存在一定的差異,四倍體Alamo和八倍體Alamo沒有明顯的遺傳分化,但四倍體Alamo的遺傳多樣性高于八倍體Alamo。本試驗結(jié)果可為柳枝稷種質(zhì)資源的開發(fā)利用、遺傳變異等方面的研究提供一定的參考。

    [1]Hou X C, Fan X F, Wu J Y,etal. Evaluation of economic benefits and ecological values of cellulosic bioenergy grasses in Beijing suburban areas. Acta Prataculturae Sinica, 2011, 20(6): 12-17.

    [2]Broz A K, Manter D K, Bowman G,etal. Plant origin and ploidy influence gene expression and life cycle characteristics in an invasive weed. BMC Plant Biology, 2009, 9(2): 33.

    [3]Yang Z, Shen Z, Tetreault H,etal. Production of autopolyploid lowland switchgrass lines throughinvitrochromosome doubling. Bioenergy Research, 2014, 7(1): 232-242.

    [4]Yang D J, Ma R J. The reviews of application of microsatellite DNA marker in crop genetic diversity research. Gansu Science and Technology, 2007, 8(8): 99-102.

    [5]Li N, Jiao Z, Qin G Y. DNA molecular marker techniques and their application in wheat breeding and genetic study. Acta Agriculturae Nucleatae Sinica, 2005, 4: 322-326.

    [6]Liu S N, Yu L Q, Zhou Y L,etal. The construction of genetic linkage frame map in tetraploidMedicagousing RAPD markers. Acta Prataculturae Sinica, 2012, 21(1): 170-175.

    [7]Wei Y L, Ni Z F, Jie C J,etal. Molecular marker location of a novel yellow rust resistance gene YrSP derived from Spelt Wheat. Journal of Agricultural Biotechnology, 2003, 1(1): 30-33.

    [8]Li H Y, Lou Y H, Hu T,etal. Genetic diversity analysis of Chinese natural tall fescue (Festucaarundinacea) germplasms by RAPD. Acta Prataculturae Sinica, 2010, 19(6): 208-214.

    [9]He L, Wang Q. Review on the PCR-Based SSR markers isolation methods. Genomics and Applied Biology, 2010, 29(4): 775-782.

    [10]Li L, He W M, Ma L P,etal. Construction Chinese Cabbage (BrassicarapaL.) core collection and its EST-SSR fingerprint database by EST-SSR molecular markers. Genomics and Applied Biology, 2009, 28(1): 76-88.

    [11]Huang L K, Bughrara S S, Zhang X Q,etal. Genetic diversity of switchgrass and its relative species inPanicumgenus using molecular markers. Biochemical Systematics & Ecology, 2011, 39(4): 685-693.

    [12]Jiang X M, Huang L K, Zhang X Q. Genetic diversity of switchgrass (PanicumvirgatumL.) germplasms detected by sCoT markers. Journal of Plant Genetic Resources, 2014, 15(1): 89-95.

    [13]Zhang Y, Huang L K, Zhang X A,etal. Analysis on genetic diversity of germplasm resources of switchgrass (Panicumvirgatum) as shown by ISSR. Acta Prataculturae Sinica, 2013, 22(6): 213-222.

    [14]Weining S, Langridge P. Identification and mapping of polymorphisms in cereals based on the polymerase chain reaction. Theoretical and Applied Genetics, 1991, 92: 209-216.

    [15]Zheng L, Zhang Z S, Chen L,etal. IT-ISJ marker and its application in construction of upland cotton linkage map. China Agriculture Science, 2008, 8(8): 2241-2248.

    [16]Lin Z, He D, Zhang X,etal. Linkage map construction and mapping QTL for cotton fiber quality using SRAP, SSR and RAPD. Plant Breeding, 2008, 124(2): 180-187.

    [17]Yu S T, Wang C T, Wang X Z,etal. Application of IT-ISJ markers in Peanut. Journal of Peanut Science, 2009, 3(3): 10-14.

    [18]Zhang J, Huang L K, Zhao B Y,etal. Comparison of SRAP and SSR markers for hybrid identification of switchgrass (Panicumvirgatum). Acta Prataculturae Sinica, 2012, 21(5): 128-133.

    [19]Lu Y F, Jiang J X, Yi Z L. Study on the transferability of maize SSR and sugarcane EST-SSR markers toMiscanthus(Poaceae). Acta Prataculturae Sinica, 2012, 21(5): 86-95.

    [20]Huang X X, Tang T, Jiang Y L,etal. Genetic diversity of wild tea plant in different altitude in Qianijazhai. Journal of Tea Science, 2015, (4): 347-353.

    [21]Yeh F C, Yang R C, Boyle T. Popgene Version 1. 31 Quick User Guide[M]. Canada: University of Alberta, and Centre for International Forestry Research, 1999.

    [22]Jiao F, Lou C F, Zhang Y Z,etal. DNA amplified polymorphism among mulberry mutants. Acta Sericologica Sinica, 2001, 27(3): 165-169.

    [23]Li X H, Wang C Y, Zhang S P,etal. Genetic polymorphism among diploid and autopolyploidy watermelons revealed by SRAP. Molecular Plant Breeding, 2007, 11: 61-63.

    [24]Wang Z W, Yu M D, Lu C,etal.AFLP molecular marker technique and its application in the study of mulberry genetics and breeding.Acta Sericologica Sinica, 2001, 27(3): 170-176.

    [25]Jiang D S. Study on Genetic Diversity and New Multiploid Germplasm of Dioscorea Zingiberensis[D]. Hunan: Hunan Agricultural University, 2006.

    [26]Du F, Yang F Y, Casler M D,etal. The research progress of switchgrass-energy grass in America. Journal of Anhui Agricultural Sciences, 2010, 38(35): 20334-20339.

    [27]Brunken J N, Estes J R. Cytological and morphological variation inPanicumvirgatumL. Southwestern Naturalist, 1975, 19(4): 379-385.

    [28]Zhong X X, Liu Z W, Liu W G,etal. Analysis of characteristics of a hexaploid somatic of hybridPennisetum. Acta Prataculturae Sinica, 2014, 23(4): 107-113.

    [29]Evans J, Crisovan E, Barry K,etal. Diversity and population structure of northern switchgrass as revealed through exome capture sequencing. Plant Journal, 2015, 4: 800-815.

    [30]Archak S, Gaikwad A B, Gautam D,etal. Comparative assessment of DNA finger printing techniques (RAPD, ISSR and AFLP) for genetic analysis of cashew (AnacardiumoccidentaleL.) accessions of India. Genomics, 2003, 46(3): 362-369.

    [1]侯新村, 范希峰, 武菊英, 等. 纖維素類能源草在京郊地區(qū)的經(jīng)濟效益與生態(tài)價值評價. 草業(yè)學報, 2011, 20(6): 12-17.

    [4]楊東娟, 馬瑞君. SSR分子標記在作物遺傳多樣性研究中的應(yīng)用現(xiàn)狀. 甘肅科技, 2007, 8(8): 99-102.

    [5]李娜, 焦湞, 秦廣雍. DNA分子標記技術(shù)及其在小麥育種及遺傳研究中的應(yīng)用. 核農(nóng)學報, 2005, 4: 322-326.

    [6]劉曙娜, 于林清, 周延林, 等. 利用RAPD技術(shù)構(gòu)建四倍體苜蓿遺傳連鎖圖譜. 草業(yè)學報, 2012, 21(1): 170-175.

    [7]魏艷玲, 倪中福, 解超杰, 等. 來自斯卑爾脫小麥新的抗條銹病基因YrSp的分子標記定位. 農(nóng)業(yè)生物技術(shù)學報, 2003, 1(1): 30-33.

    [8]李惠英, 婁燕宏, 胡濤, 等. 中國高羊茅種質(zhì)資源遺傳多樣性的RAPD分析. 草業(yè)學報, 2010, 19(6): 208-214.

    [9]何琳, 王群. 基于PCR的SSR標記分離方法綜述. 基因組學與應(yīng)用生物學, 2010, 29(4): 775-782.

    [10]李麗, 何偉明, 馬連平, 等. 用EST-SSR分子標記技術(shù)構(gòu)建大白菜核心種質(zhì)及其指紋圖譜庫. 基因組學與應(yīng)用生物學, 2009, 28(1): 76-88.

    [12]蔣曉梅, 黃琳凱, 張新全, 等. 利用SCoT分析柳枝稷遺傳資源的多樣性. 植物遺傳資源學報, 2014, 15(1): 89-95.

    [13]張瑜, 黃琳凱, 張新全, 等. 柳枝稷種質(zhì)資源遺傳多樣性的ISSR分析. 草業(yè)學報, 2013, 22(6): 213-222.

    [15]鄭靚, 張正圣, 陳利, 等. IT-ISJ標記及其在陸地棉遺傳圖譜構(gòu)建中的應(yīng)用. 中國農(nóng)業(yè)科學, 2008, 8(8): 2241-2248.

    [17]于樹濤, 王傳堂, 王秀貞, 等. IT-ISJ標記在花生上的應(yīng)用研究. 花生學報, 2009, 3(3): 10-14.

    [18]張婧, 黃琳凱, Zhao B Y, 等. SRAP與SSR分子標記在柳枝稷雜種鑒定中的比較分析. 草業(yè)學報, 2012, 21(5): 128-133.

    [19]盧玉飛, 蔣建雄, 易自力. 玉米SSR引物和甘蔗EST-SSR引物在芒屬中的通用性研究. 草業(yè)學報, 2012, 21(5): 86-95.

    [20]黃曉霞, 唐探, 姜永雷, 等. 千家寨不同海拔野生茶樹的 EST-SSR 遺傳多樣性研究. 茶葉科學, 2015, (4): 347-353.

    [22]焦鋒, 樓程富, 張有做, 等. 桑樹變異株系基因組DNA擴增多態(tài)性(RAPD)研究. 蠶業(yè)科學, 2001, 27(3): 165-169.

    [23]李曉慧, 王從彥, 張四普, 等. 西瓜二倍體及同源多倍體SRAP多態(tài)性分析. 分子植物育種, 2007, 11: 61-63.

    [24]王卓偉,余茂德,魯成, 等.桑樹多倍體育種材料遺傳背景的AFLP分析.蠶業(yè)科學, 2001, 27(3): 170-176.

    [25]蔣道松. 盾葉薯蕷遺傳多樣性及多倍體新種質(zhì)培育的研究[D]. 長沙: 湖南農(nóng)業(yè)大學, 2006.

    [26]杜菲, 楊富裕, Casler M D, 等. 美國能源草柳枝稷的研究進展. 安徽農(nóng)業(yè)科學, 2010, 38(35): 20334-20339.

    [28]鐘小仙, 劉智微, 劉偉國, 等. 六倍體雜交狼尾草體細胞突變體特異性分析. 草業(yè)學報, 2014, 23(4): 107-113.

    Genetic diversity assessment of tetraploid and octoploid switchgrass using EST-SSR and IT-ISJ molecular markers

    ZHONG Ming**, ZHOU Si-Fan**, ZHANG Xin-Quan, HUANG Xiu, YAN Hai-Dong, ZHANG Ai-Ling, LUO-Qian, HUANG Lin-Kai*

    Department of Grassland Science, Sichuan Agricultural University, Chengdu 611130, China

    Switchgrass is an important bioenergy plant. In this study, the genetic diversity of 65 samples of tetraploid Alamo and mutated octoploid Alamo (using colchicine) was assessed using expressed sequence tag-simple sequence repeat (EST-SSR) and intron-targeted intron-exon splice junction (IT-ISJ) markers. The results showed that EST-SSR primers amplified 245 fragments; the percentage of polymorphic loci was 99.59%, the number of alleles (Na) was 1.9878, the effective number of alleles (Ne) was 1.1718, Nei’s gene diversity (H) was 0.1308, and Shannon’s information index (I) was 0.2388. These genetic parameters for the tetraploid Alamo were as follows: 226 fragments; 88.93% polymorphic loci; Na, 1.8816; Ne, 1.1823; H, 0.1349; and I, 0.2398. These values were higher than those obtained for octoploid Alamo: 217 fragments; 92.24% polymorphic loci; Na, 1.9184; Ne, 1.1652; H, 0.1228; and I, 0.2234. The IT-ISJ primers amplified 165 fragments and the percentage of polymorphic loci was 96.36%, Na was 1.9878, Ne was 1.5777, H was 0.3354, and I was 0.4999. These parameters for tetraploid Alamo were as follows: 145 fragments; 94.55% polymorphic loci; Na, 1.9333; Ne, 1.6366; H, 0.3602; and I, 0.5280. These values were higher than those obtained for octoploid Alamo: 139 fragments; 84.24% polymorphic loci; Na, 1.7515; Ne, 1.4337; H, 0.2573; and I, 0.3864. A clustering analysis with the unweighted pair-group method with arithmetic means did not clearly separate 68 samples of tetraploid Alamo and octoploid Alamo, which had a low level of genetic differentiation. An analysis of molecular variance showed that 94.78% and 80.76% of genetic diversity occurred within ploidy levels, while 5.22% and 19.24% of genetic diversity occurred between ploidy levels.

    switchgrass (Panicumvirgatum); EST-SSR (expressed sequence tag-simple sequence repeat); IT-ISJ (intron targeted intron-exon splice junction); Alamo; molecular marker; genetic diversity

    10.11686/cyxb2015554

    2015-12-07;改回日期:2016-01-26

    科技部“十二五”863計劃(2012AA101801-2)和國家自然基金(31201845)資助。

    鐘鳴(1987-),女,四川富順人,碩士。E-mail:zmzl7758@163.com。周思凡(1993-),女,重慶人,在讀碩士。E-mail:swanchou93@163.com。**共同第一作者These authors contributed equally to this work.

    Corresponding author. E-mail:huanglingkai@sicau.edu.cn

    http://cyxb.lzu.edu.cn

    鐘鳴, 周思凡, 張新全, 黃秀, 嚴海東, 張愛玲,羅茜,黃琳凱. 利用EST-SSR、IT-ISJ分子標記研究四倍體與八倍體柳枝稷的遺傳多樣性. 草業(yè)學報, 2016, 25(10): 113-123.

    ZHONG Ming, ZHOU Si-Fan, ZHANG Xin-Quan, HUANG Xiu, YAN Hai-Dong, ZHANG Ai-Ling, LUO-Qian, HUANG Lin-Kai. Genetic diversity assessment of tetraploid and octoploid switchgrass using EST-SSR and IT-ISJ molecular markers. Acta Prataculturae Sinica, 2016, 25(10): 113-123.

    猜你喜歡
    四倍體柳枝多態(tài)性
    小果型西瓜四倍體誘變及其鑒定分析
    柳枝上的春天
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進展
    柳枝柳枝 變變變
    生命力
    文苑(2020年2期)2020-11-22 01:45:32
    生命力
    四倍體泡桐育苗高效管理技術(shù)綜述
    四倍體巴戟天根的結(jié)構(gòu)與其蒽醌類化合物的關(guān)系
    中成藥(2016年4期)2016-05-17 06:07:54
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    GlobalFiler~? PCR擴增試劑盒驗證及其STR遺傳多態(tài)性
    香蕉丝袜av| 国产一级毛片七仙女欲春2| 可以在线观看的亚洲视频| 黄色视频,在线免费观看| 午夜福利免费观看在线| 免费看美女性在线毛片视频| 午夜福利免费观看在线| 国内毛片毛片毛片毛片毛片| 舔av片在线| 最近视频中文字幕2019在线8| 免费观看人在逋| 欧美性感艳星| 90打野战视频偷拍视频| 91在线精品国自产拍蜜月 | 日本在线视频免费播放| 日本熟妇午夜| 成人国产一区最新在线观看| 97超视频在线观看视频| 日韩欧美精品免费久久 | 国产高清三级在线| 亚洲第一电影网av| 日本黄色片子视频| 51午夜福利影视在线观看| 十八禁网站免费在线| 婷婷六月久久综合丁香| 三级国产精品欧美在线观看| 亚洲国产精品合色在线| 白带黄色成豆腐渣| 亚洲专区国产一区二区| 亚洲成人免费电影在线观看| 日本黄大片高清| 精品电影一区二区在线| 欧美又色又爽又黄视频| 2021天堂中文幕一二区在线观| 午夜福利在线在线| 精品国产三级普通话版| 日本 av在线| 最近视频中文字幕2019在线8| av女优亚洲男人天堂| 色尼玛亚洲综合影院| 国产精品亚洲美女久久久| 国产视频一区二区在线看| 在线天堂最新版资源| 三级男女做爰猛烈吃奶摸视频| 熟妇人妻久久中文字幕3abv| 国产激情欧美一区二区| 亚洲欧美日韩高清在线视频| 国产69精品久久久久777片| 99久久精品国产亚洲精品| 日本撒尿小便嘘嘘汇集6| 99久国产av精品| 一本精品99久久精品77| 一区二区三区国产精品乱码| 伊人久久精品亚洲午夜| 国产精品久久久久久久电影 | 成人高潮视频无遮挡免费网站| 久久久久久久久中文| 老熟妇乱子伦视频在线观看| 草草在线视频免费看| 亚洲av电影在线进入| 熟女人妻精品中文字幕| 成人永久免费在线观看视频| 日韩欧美 国产精品| av国产免费在线观看| 欧美又色又爽又黄视频| 亚洲 国产 在线| 一个人看视频在线观看www免费 | 丰满的人妻完整版| 欧美一区二区精品小视频在线| 国内精品美女久久久久久| 日韩人妻高清精品专区| 亚洲成人精品中文字幕电影| 免费电影在线观看免费观看| 欧美日韩中文字幕国产精品一区二区三区| 国产黄色小视频在线观看| www.www免费av| 非洲黑人性xxxx精品又粗又长| 老司机午夜福利在线观看视频| 最新美女视频免费是黄的| 国产伦精品一区二区三区四那| 国产精品乱码一区二三区的特点| 午夜福利在线观看免费完整高清在 | 手机成人av网站| 无遮挡黄片免费观看| 在线免费观看的www视频| 色综合站精品国产| a级毛片a级免费在线| 久久这里只有精品中国| 不卡一级毛片| 精品一区二区三区人妻视频| 久久国产精品影院| 色哟哟哟哟哟哟| 波多野结衣高清无吗| 成人三级黄色视频| 久久久久久久久大av| 麻豆国产97在线/欧美| 久久久精品欧美日韩精品| 久久久久久久精品吃奶| 日韩亚洲欧美综合| 夜夜夜夜夜久久久久| 国产精品美女特级片免费视频播放器| 欧美中文综合在线视频| 亚洲国产日韩欧美精品在线观看 | 久久久色成人| 麻豆一二三区av精品| 欧美一区二区精品小视频在线| 欧美性猛交黑人性爽| 免费在线观看亚洲国产| 欧美日韩精品网址| 日本精品一区二区三区蜜桃| 成人av在线播放网站| 午夜日韩欧美国产| 国产一区二区在线观看日韩 | 中文字幕高清在线视频| 三级国产精品欧美在线观看| 九九久久精品国产亚洲av麻豆| 国产一区在线观看成人免费| 亚洲国产精品sss在线观看| ponron亚洲| 变态另类成人亚洲欧美熟女| 国语自产精品视频在线第100页| 99久久99久久久精品蜜桃| 亚洲av二区三区四区| 一区福利在线观看| 午夜福利视频1000在线观看| 欧美高清成人免费视频www| 国产av在哪里看| 免费观看人在逋| 一级黄片播放器| 国产在线精品亚洲第一网站| 国产精品影院久久| 久久久久久久精品吃奶| 久久久久精品国产欧美久久久| 在线天堂最新版资源| 亚洲国产欧洲综合997久久,| 一级黄片播放器| 两性午夜刺激爽爽歪歪视频在线观看| 欧美bdsm另类| 国产淫片久久久久久久久 | 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网| 在线观看一区二区三区| 中文字幕人成人乱码亚洲影| 亚洲真实伦在线观看| 夜夜看夜夜爽夜夜摸| 欧美三级亚洲精品| 身体一侧抽搐| 欧美日本亚洲视频在线播放| 一个人免费在线观看的高清视频| 久久草成人影院| 亚洲欧美日韩东京热| 国产一区在线观看成人免费| 男人和女人高潮做爰伦理| 亚洲在线观看片| 亚洲专区中文字幕在线| а√天堂www在线а√下载| 熟女少妇亚洲综合色aaa.| 国产伦在线观看视频一区| 亚洲 国产 在线| 免费看美女性在线毛片视频| 99久久成人亚洲精品观看| 欧美三级亚洲精品| 我要搜黄色片| 亚洲国产精品久久男人天堂| 欧美一级a爱片免费观看看| 两个人看的免费小视频| 制服人妻中文乱码| 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av涩爱 | 色视频www国产| 天天一区二区日本电影三级| 亚洲精品美女久久久久99蜜臀| 精品熟女少妇八av免费久了| 白带黄色成豆腐渣| АⅤ资源中文在线天堂| 色综合欧美亚洲国产小说| 一级黄色大片毛片| 欧美黑人欧美精品刺激| 老司机午夜福利在线观看视频| 老熟妇乱子伦视频在线观看| 久久久久久久久大av| 国产真实乱freesex| 9191精品国产免费久久| 又粗又爽又猛毛片免费看| 亚洲精品国产精品久久久不卡| 久久久久性生活片| 亚洲av二区三区四区| 88av欧美| 中出人妻视频一区二区| 身体一侧抽搐| 色播亚洲综合网| 日韩大尺度精品在线看网址| 国产精品99久久99久久久不卡| 成人国产一区最新在线观看| 精品久久久久久,| 中文字幕人成人乱码亚洲影| 少妇丰满av| 高清日韩中文字幕在线| 俄罗斯特黄特色一大片| 黄片大片在线免费观看| 精品日产1卡2卡| 一级毛片高清免费大全| 不卡一级毛片| 久久精品国产亚洲av香蕉五月| 香蕉丝袜av| aaaaa片日本免费| 亚洲在线观看片| 日韩中文字幕欧美一区二区| 男人舔女人下体高潮全视频| 超碰av人人做人人爽久久 | 黄色丝袜av网址大全| 欧美不卡视频在线免费观看| 亚洲精品日韩av片在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 日韩人妻高清精品专区| 亚洲精品久久国产高清桃花| 一级毛片高清免费大全| 午夜免费观看网址| 毛片女人毛片| 国产高清三级在线| 特级一级黄色大片| 在线观看66精品国产| 国产一区二区三区视频了| 日韩欧美精品v在线| 91在线精品国自产拍蜜月 | 欧美性猛交╳xxx乱大交人| 国产高清videossex| 国产精品99久久99久久久不卡| 中文字幕人妻丝袜一区二区| 精品99又大又爽又粗少妇毛片 | 亚洲av五月六月丁香网| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久成人av| 亚洲精品在线观看二区| 十八禁人妻一区二区| av专区在线播放| 亚洲成人中文字幕在线播放| 给我免费播放毛片高清在线观看| 日韩国内少妇激情av| 国产91精品成人一区二区三区| 亚洲成人久久爱视频| 网址你懂的国产日韩在线| 亚洲国产欧美人成| 亚洲人成网站高清观看| 免费av毛片视频| 很黄的视频免费| 国产精品亚洲一级av第二区| 免费人成视频x8x8入口观看| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 欧美大码av| 亚洲成人精品中文字幕电影| 18禁黄网站禁片午夜丰满| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩中文字幕国产精品一区二区三区| 久久久精品大字幕| 国产欧美日韩一区二区三| 久久久久性生活片| 日韩人妻高清精品专区| 欧美色视频一区免费| 欧美不卡视频在线免费观看| 一夜夜www| 国产精品电影一区二区三区| 在线免费观看的www视频| 精品日产1卡2卡| 欧美又色又爽又黄视频| 一进一出抽搐动态| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月婷婷丁香| 国产精品乱码一区二三区的特点| 深爱激情五月婷婷| 少妇的逼水好多| 免费电影在线观看免费观看| 久久久久精品国产欧美久久久| 在线视频色国产色| 欧美色欧美亚洲另类二区| 我的老师免费观看完整版| 夜夜爽天天搞| 天天一区二区日本电影三级| 观看美女的网站| 看片在线看免费视频| 美女黄网站色视频| 亚洲精品乱码久久久v下载方式 | 香蕉丝袜av| 性欧美人与动物交配| 人妻夜夜爽99麻豆av| 最近最新中文字幕大全免费视频| 又黄又粗又硬又大视频| 叶爱在线成人免费视频播放| 中文字幕高清在线视频| 深夜精品福利| 嫩草影院入口| 一个人看的www免费观看视频| 乱人视频在线观看| 9191精品国产免费久久| 中文字幕熟女人妻在线| 99热精品在线国产| 99在线人妻在线中文字幕| 丰满人妻熟妇乱又伦精品不卡| 搡老岳熟女国产| 国产日本99.免费观看| 人妻夜夜爽99麻豆av| 18禁裸乳无遮挡免费网站照片| 伊人久久大香线蕉亚洲五| 高清毛片免费观看视频网站| 亚洲va日本ⅴa欧美va伊人久久| 少妇的逼好多水| 一本精品99久久精品77| 欧美又色又爽又黄视频| 一本一本综合久久| a级一级毛片免费在线观看| 香蕉av资源在线| 亚洲国产中文字幕在线视频| 五月伊人婷婷丁香| 老司机午夜十八禁免费视频| 日韩欧美 国产精品| 夜夜爽天天搞| 国产精品98久久久久久宅男小说| 99国产精品一区二区蜜桃av| 特级一级黄色大片| 国产精品自产拍在线观看55亚洲| 最近最新中文字幕大全电影3| 男人舔女人下体高潮全视频| 色综合欧美亚洲国产小说| 亚洲成人久久性| 亚洲最大成人手机在线| 岛国在线观看网站| 天美传媒精品一区二区| 午夜视频国产福利| 欧美黑人欧美精品刺激| 欧美中文综合在线视频| 亚洲av成人不卡在线观看播放网| 亚洲成人久久爱视频| h日本视频在线播放| av视频在线观看入口| 日韩大尺度精品在线看网址| 搡老岳熟女国产| 成人鲁丝片一二三区免费| 久久精品91蜜桃| 啦啦啦免费观看视频1| 欧美日韩一级在线毛片| av天堂中文字幕网| 中国美女看黄片| 人人妻人人看人人澡| 中文字幕人成人乱码亚洲影| 人妻久久中文字幕网| 成人av一区二区三区在线看| 日本 欧美在线| 国产精品乱码一区二三区的特点| 国产精品爽爽va在线观看网站| 日日摸夜夜添夜夜添小说| 久久久久久久久大av| 最新在线观看一区二区三区| 久久久久亚洲av毛片大全| 嫁个100分男人电影在线观看| 日日摸夜夜添夜夜添小说| 久久久久久久久大av| 男人的好看免费观看在线视频| 最新在线观看一区二区三区| 91九色精品人成在线观看| 色在线成人网| 久久香蕉国产精品| 黄色片一级片一级黄色片| 国产午夜精品论理片| 99热只有精品国产| 国产高清有码在线观看视频| 欧美极品一区二区三区四区| 久久久久免费精品人妻一区二区| 免费在线观看影片大全网站| 精品国产美女av久久久久小说| 中国美女看黄片| 高清在线国产一区| 成人国产一区最新在线观看| 极品教师在线免费播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美日本亚洲视频在线播放| 欧美性猛交╳xxx乱大交人| 在线看三级毛片| av福利片在线观看| 露出奶头的视频| 国产精品久久久久久久久免 | av天堂中文字幕网| 俄罗斯特黄特色一大片| www日本在线高清视频| 国产一区二区在线av高清观看| 精品久久久久久成人av| 国产蜜桃级精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 精品不卡国产一区二区三区| 亚洲欧美日韩卡通动漫| 国产熟女xx| 欧美午夜高清在线| 国产精品久久电影中文字幕| 美女免费视频网站| 成熟少妇高潮喷水视频| 99国产极品粉嫩在线观看| 日韩高清综合在线| 国产高潮美女av| 欧美黑人巨大hd| 国产精品98久久久久久宅男小说| 18禁美女被吸乳视频| 内射极品少妇av片p| 久久久国产精品麻豆| 真实男女啪啪啪动态图| 在线免费观看的www视频| 亚洲av不卡在线观看| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频| 欧美日韩精品网址| 午夜免费激情av| 手机成人av网站| 国模一区二区三区四区视频| 日韩精品中文字幕看吧| 女人十人毛片免费观看3o分钟| 99久久成人亚洲精品观看| 波多野结衣巨乳人妻| 日本一本二区三区精品| 在线天堂最新版资源| 99精品欧美一区二区三区四区| 18禁黄网站禁片午夜丰满| 日韩国内少妇激情av| 天天添夜夜摸| 99久久成人亚洲精品观看| 日韩精品青青久久久久久| 少妇裸体淫交视频免费看高清| 又爽又黄无遮挡网站| 又粗又爽又猛毛片免费看| 国产高清有码在线观看视频| 人妻丰满熟妇av一区二区三区| 亚洲激情在线av| 一级作爱视频免费观看| 香蕉久久夜色| 九九久久精品国产亚洲av麻豆| 51国产日韩欧美| 白带黄色成豆腐渣| 亚洲精品色激情综合| 国产成人欧美在线观看| 97碰自拍视频| 国产一区二区激情短视频| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 一进一出抽搐gif免费好疼| 国内毛片毛片毛片毛片毛片| 日韩欧美国产在线观看| 99久久久亚洲精品蜜臀av| 国产精品久久视频播放| 人妻丰满熟妇av一区二区三区| 国产三级在线视频| 亚洲七黄色美女视频| 九色国产91popny在线| 亚洲真实伦在线观看| av在线天堂中文字幕| 又爽又黄无遮挡网站| 久久国产精品影院| 成人午夜高清在线视频| 国产精品 欧美亚洲| 国产精品综合久久久久久久免费| 一个人免费在线观看的高清视频| 首页视频小说图片口味搜索| 精品一区二区三区人妻视频| 九九热线精品视视频播放| 性欧美人与动物交配| 熟女电影av网| 一二三四社区在线视频社区8| 一边摸一边抽搐一进一小说| 99久久精品一区二区三区| 欧美又色又爽又黄视频| 性色av乱码一区二区三区2| 法律面前人人平等表现在哪些方面| 久久精品亚洲精品国产色婷小说| 午夜免费成人在线视频| 欧美zozozo另类| 欧美另类亚洲清纯唯美| 欧美区成人在线视频| 国产精品av视频在线免费观看| 黄色片一级片一级黄色片| 欧美日韩瑟瑟在线播放| 色吧在线观看| 狂野欧美白嫩少妇大欣赏| 国产视频一区二区在线看| 欧美一区二区国产精品久久精品| 特大巨黑吊av在线直播| 国产午夜福利久久久久久| 最近最新免费中文字幕在线| 日日摸夜夜添夜夜添小说| www日本黄色视频网| 真人一进一出gif抽搐免费| 一级毛片女人18水好多| 97超视频在线观看视频| 国产亚洲欧美98| 亚洲专区中文字幕在线| 99久久综合精品五月天人人| 日韩成人在线观看一区二区三区| 日本黄色片子视频| 丰满的人妻完整版| 国产在视频线在精品| 亚洲精品国产精品久久久不卡| 亚洲在线观看片| 久久精品国产清高在天天线| 在线国产一区二区在线| 90打野战视频偷拍视频| 国产乱人伦免费视频| 成人午夜高清在线视频| 少妇熟女aⅴ在线视频| 久久久久九九精品影院| 九九在线视频观看精品| 久久久久性生活片| 亚洲真实伦在线观看| 亚洲第一电影网av| 国产蜜桃级精品一区二区三区| 色av中文字幕| 国产成人福利小说| 精品国内亚洲2022精品成人| 午夜福利免费观看在线| 中文字幕久久专区| 精品久久久久久久毛片微露脸| 老司机深夜福利视频在线观看| 两个人视频免费观看高清| 成人永久免费在线观看视频| 国产精华一区二区三区| 精品一区二区三区视频在线观看免费| 国产精品亚洲av一区麻豆| 午夜福利18| 久久人妻av系列| 亚洲,欧美精品.| 国产亚洲精品一区二区www| 天天添夜夜摸| 国产亚洲精品一区二区www| a在线观看视频网站| 国产亚洲精品久久久久久毛片| 欧美色欧美亚洲另类二区| 亚洲乱码一区二区免费版| 在线免费观看不下载黄p国产 | 亚洲欧美精品综合久久99| av国产免费在线观看| 亚洲人成网站高清观看| 十八禁网站免费在线| 少妇熟女aⅴ在线视频| 亚洲欧美一区二区三区黑人| 亚洲成人久久性| 亚洲成人久久爱视频| 丁香欧美五月| 亚洲avbb在线观看| 精品人妻1区二区| 制服人妻中文乱码| 亚洲av成人av| 手机成人av网站| 成人无遮挡网站| 国产精品影院久久| 国产亚洲欧美98| 深夜精品福利| 国产精品野战在线观看| 亚洲黑人精品在线| 亚洲 欧美 日韩 在线 免费| 国产单亲对白刺激| 五月伊人婷婷丁香| 校园春色视频在线观看| 亚洲内射少妇av| 天堂av国产一区二区熟女人妻| 国产精品 国内视频| 欧美日韩综合久久久久久 | 国产真人三级小视频在线观看| 午夜福利成人在线免费观看| 女生性感内裤真人,穿戴方法视频| 日韩欧美在线乱码| 男插女下体视频免费在线播放| 动漫黄色视频在线观看| 欧美成狂野欧美在线观看| www日本黄色视频网| 国内久久婷婷六月综合欲色啪| 中文字幕熟女人妻在线| 又黄又粗又硬又大视频| av女优亚洲男人天堂| 2021天堂中文幕一二区在线观| 99热这里只有是精品50| 床上黄色一级片| 免费av观看视频| 国产精品自产拍在线观看55亚洲| 伊人久久精品亚洲午夜| 亚洲精华国产精华精| 人妻夜夜爽99麻豆av| 最新在线观看一区二区三区| 天堂√8在线中文| 日日夜夜操网爽| 久久精品91无色码中文字幕| 国产一区二区三区视频了| 精品99又大又爽又粗少妇毛片 | 韩国av一区二区三区四区| 亚洲五月婷婷丁香| 亚洲专区中文字幕在线| 制服人妻中文乱码| 亚洲18禁久久av| 欧美日韩精品网址| 禁无遮挡网站| 最新中文字幕久久久久| 日韩欧美国产在线观看| 亚洲av第一区精品v没综合| 亚洲av二区三区四区| 欧美3d第一页| 五月伊人婷婷丁香| 中文字幕人妻熟人妻熟丝袜美 | 看免费av毛片| 国产三级黄色录像| 又紧又爽又黄一区二区| 青草久久国产| 热99在线观看视频| 天堂网av新在线| 中亚洲国语对白在线视频| 亚洲一区二区三区色噜噜| 久久伊人香网站| 长腿黑丝高跟| 国产老妇女一区| 麻豆国产av国片精品| 搡女人真爽免费视频火全软件 | 手机成人av网站|