• <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| 成人亚洲精品一区在线观看| 欧美日韩国产mv在线观看视频| 国产一区二区三区综合在线观看| 精品久久久精品久久久| 国产精品美女特级片免费视频播放器 | 中文亚洲av片在线观看爽| 久久九九热精品免费| 在线观看免费视频网站a站| 麻豆一二三区av精品| 欧美成狂野欧美在线观看| 国产主播在线观看一区二区| netflix在线观看网站| 亚洲国产欧美一区二区综合| 超碰97精品在线观看| 午夜成年电影在线免费观看| 日韩视频一区二区在线观看| 精品国产亚洲在线| 国产精品一区二区在线不卡| 91av网站免费观看| 黄色成人免费大全| 久久精品国产99精品国产亚洲性色 | 91精品国产国语对白视频| 亚洲精华国产精华精| 免费看a级黄色片| av天堂久久9| 一级,二级,三级黄色视频| 亚洲伊人色综图| 99久久99久久久精品蜜桃| 亚洲欧美激情综合另类| 国产精品九九99| 午夜激情av网站| 欧美久久黑人一区二区| 在线免费观看的www视频| 精品卡一卡二卡四卡免费| a级毛片在线看网站| 成人三级黄色视频| 免费观看精品视频网站| 国产有黄有色有爽视频| 久久精品aⅴ一区二区三区四区| 精品卡一卡二卡四卡免费| 亚洲人成电影观看| 亚洲av成人不卡在线观看播放网| 热99re8久久精品国产| 久久久久亚洲av毛片大全| 黄色丝袜av网址大全| 亚洲av电影在线进入| 国产精品98久久久久久宅男小说| 国产单亲对白刺激| 香蕉丝袜av| 最好的美女福利视频网| 午夜91福利影院| bbb黄色大片| 午夜福利欧美成人| 在线观看舔阴道视频| 两人在一起打扑克的视频| 欧美激情极品国产一区二区三区| 午夜视频精品福利| 日韩欧美在线二视频| 亚洲精品美女久久久久99蜜臀| 老司机亚洲免费影院| 日本撒尿小便嘘嘘汇集6| 免费人成视频x8x8入口观看| 婷婷六月久久综合丁香| 91九色精品人成在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利免费观看在线| 激情在线观看视频在线高清| 亚洲精品国产一区二区精华液| 亚洲av熟女| 国产成人系列免费观看| 欧美大码av| 男女之事视频高清在线观看| 久久国产乱子伦精品免费另类| 色播在线永久视频| 女警被强在线播放| 五月开心婷婷网| 日本五十路高清| 欧美国产精品va在线观看不卡| 久久久国产一区二区| 一级,二级,三级黄色视频| 19禁男女啪啪无遮挡网站| 欧美亚洲日本最大视频资源| 看片在线看免费视频| 亚洲人成77777在线视频| 亚洲精品成人av观看孕妇| 在线观看午夜福利视频| 亚洲成人免费av在线播放| 制服诱惑二区| 精品国产亚洲在线| 最近最新中文字幕大全电影3 | 国产精品秋霞免费鲁丝片| 亚洲狠狠婷婷综合久久图片| 精品国产亚洲在线| 亚洲第一欧美日韩一区二区三区| 99re在线观看精品视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av日韩精品久久久久久密| 亚洲午夜理论影院| 日韩精品免费视频一区二区三区| av电影中文网址| 久久久国产成人免费| 国产aⅴ精品一区二区三区波| 宅男免费午夜| 成年女人毛片免费观看观看9| 男人舔女人的私密视频| 黑人欧美特级aaaaaa片| 操美女的视频在线观看| 视频区欧美日本亚洲| 亚洲第一av免费看| 久久精品亚洲精品国产色婷小说| 久久久久九九精品影院| 黄色视频不卡| 丰满的人妻完整版| 黄色怎么调成土黄色| 一级作爱视频免费观看| 黑丝袜美女国产一区| 女人被躁到高潮嗷嗷叫费观| 91老司机精品| 一进一出抽搐动态| 另类亚洲欧美激情| 国产精品国产高清国产av| 精品午夜福利视频在线观看一区| 亚洲熟妇中文字幕五十中出 | 亚洲aⅴ乱码一区二区在线播放 | 欧美人与性动交α欧美软件| 亚洲aⅴ乱码一区二区在线播放 | 色精品久久人妻99蜜桃| 欧美一区二区精品小视频在线| 怎么达到女性高潮| 久99久视频精品免费| 亚洲av日韩精品久久久久久密| 国产一区在线观看成人免费| 97超级碰碰碰精品色视频在线观看| 亚洲狠狠婷婷综合久久图片| 三上悠亚av全集在线观看| 丰满迷人的少妇在线观看| 男女之事视频高清在线观看| 欧美黑人精品巨大| a级毛片黄视频| 韩国av一区二区三区四区| 搡老乐熟女国产| 日韩欧美免费精品| 亚洲中文av在线| 国产成人精品久久二区二区免费| a在线观看视频网站| 美女高潮到喷水免费观看| 巨乳人妻的诱惑在线观看| 国产视频一区二区在线看| 50天的宝宝边吃奶边哭怎么回事| 狂野欧美激情性xxxx| 一边摸一边抽搐一进一小说| 琪琪午夜伦伦电影理论片6080| 午夜a级毛片| 老司机靠b影院| 激情在线观看视频在线高清| 精品久久久久久久久久免费视频 | 一级a爱片免费观看的视频| 亚洲av熟女| 精品日产1卡2卡| 十八禁网站免费在线| 女性被躁到高潮视频| 国产成人精品无人区| 日日爽夜夜爽网站| 免费看十八禁软件| 日韩欧美三级三区| 激情在线观看视频在线高清| 校园春色视频在线观看| 精品第一国产精品| 日本a在线网址| 另类亚洲欧美激情| 日韩免费av在线播放| 欧美乱色亚洲激情| 正在播放国产对白刺激| 午夜精品久久久久久毛片777| 夜夜爽天天搞| 国产精品免费一区二区三区在线| 国产av在哪里看| 欧美最黄视频在线播放免费 | 搡老乐熟女国产| 中文字幕人妻丝袜一区二区| 成人国产一区最新在线观看| 老司机在亚洲福利影院| 成熟少妇高潮喷水视频| 国产亚洲欧美在线一区二区| 日韩成人在线观看一区二区三区| 又黄又爽又免费观看的视频| 丝袜人妻中文字幕| 手机成人av网站| 又黄又粗又硬又大视频| 一级毛片精品| 久久 成人 亚洲| 欧美日韩亚洲国产一区二区在线观看| 午夜免费激情av| 不卡一级毛片| 中文字幕最新亚洲高清| а√天堂www在线а√下载| 俄罗斯特黄特色一大片| 啦啦啦免费观看视频1| 亚洲专区中文字幕在线| 免费观看精品视频网站| 99精国产麻豆久久婷婷| 99re在线观看精品视频| 国产99白浆流出| 黄片播放在线免费| 性少妇av在线| www.www免费av| 亚洲精品一二三| 国产精品影院久久| 久久香蕉精品热| 一边摸一边抽搐一进一出视频| 免费在线观看日本一区| 久久久国产成人免费| 久久天躁狠狠躁夜夜2o2o| 91av网站免费观看| 久久久精品欧美日韩精品| 久久香蕉国产精品| 国产黄a三级三级三级人| 无遮挡黄片免费观看| 国产深夜福利视频在线观看| 日韩欧美国产一区二区入口| 丝袜在线中文字幕| 丰满饥渴人妻一区二区三| 在线观看免费视频日本深夜| 亚洲色图av天堂| 免费日韩欧美在线观看| 国产伦一二天堂av在线观看| 91精品国产国语对白视频| 午夜免费鲁丝| 窝窝影院91人妻| 成人18禁高潮啪啪吃奶动态图| 法律面前人人平等表现在哪些方面| 国产黄色免费在线视频| 亚洲 欧美 日韩 在线 免费| 777久久人妻少妇嫩草av网站| 国产成+人综合+亚洲专区| 中文字幕高清在线视频| 亚洲午夜理论影院| 久久精品国产综合久久久| 女警被强在线播放| 免费女性裸体啪啪无遮挡网站| 欧美激情 高清一区二区三区| 欧美黑人欧美精品刺激| 三上悠亚av全集在线观看| 一本大道久久a久久精品| 国产伦一二天堂av在线观看| 国产单亲对白刺激| 国产单亲对白刺激| 99香蕉大伊视频| 一区在线观看完整版| 男人操女人黄网站| 免费看十八禁软件| 19禁男女啪啪无遮挡网站| 亚洲自拍偷在线| 午夜福利影视在线免费观看| 午夜精品国产一区二区电影| 久久 成人 亚洲| 午夜久久久在线观看| 露出奶头的视频| 涩涩av久久男人的天堂| 亚洲国产毛片av蜜桃av| 欧美 亚洲 国产 日韩一| 满18在线观看网站| 免费高清视频大片| 一区二区三区国产精品乱码| 亚洲男人天堂网一区| 午夜福利,免费看| 国产三级黄色录像| 日本欧美视频一区| 韩国av一区二区三区四区| 可以在线观看毛片的网站| 脱女人内裤的视频| 免费观看精品视频网站| 欧美久久黑人一区二区| 久久热在线av| 久久午夜综合久久蜜桃| 亚洲国产中文字幕在线视频| 在线观看免费视频网站a站| 一区二区日韩欧美中文字幕| 国产91精品成人一区二区三区| 国产不卡一卡二| 欧美不卡视频在线免费观看 | 欧美日韩亚洲高清精品| 精品高清国产在线一区| 精品国内亚洲2022精品成人| 久久草成人影院| 免费不卡黄色视频| 国内久久婷婷六月综合欲色啪| 国产精品一区二区在线不卡| av在线天堂中文字幕 | 又大又爽又粗| 中文字幕人妻丝袜制服| 满18在线观看网站| 久久精品成人免费网站| 久久九九热精品免费| 免费人成视频x8x8入口观看| 99国产精品一区二区蜜桃av| 亚洲成人国产一区在线观看| 国产精品一区二区三区四区久久 | 操美女的视频在线观看| 国产高清激情床上av| 妹子高潮喷水视频| 两个人看的免费小视频| av国产精品久久久久影院| 最新美女视频免费是黄的| 日本a在线网址| 国产精品久久久av美女十八| 母亲3免费完整高清在线观看| 欧美激情极品国产一区二区三区| 母亲3免费完整高清在线观看| 欧美日韩av久久| 伦理电影免费视频| 日本一区二区免费在线视频| 国产色视频综合| 日韩有码中文字幕| 满18在线观看网站| 99国产综合亚洲精品| 精品免费久久久久久久清纯| 女生性感内裤真人,穿戴方法视频| 成人亚洲精品一区在线观看| 国产精品自产拍在线观看55亚洲| 岛国视频午夜一区免费看| 久久伊人香网站| 久久青草综合色| 9191精品国产免费久久| 精品欧美一区二区三区在线| 日本vs欧美在线观看视频| 欧美中文日本在线观看视频| 成年人黄色毛片网站| 在线观看免费午夜福利视频| 免费在线观看日本一区| 亚洲av第一区精品v没综合| 欧美乱码精品一区二区三区| 久久精品影院6| 亚洲av日韩精品久久久久久密| 欧美日韩亚洲综合一区二区三区_| 夜夜爽天天搞| 午夜精品久久久久久毛片777| 成在线人永久免费视频| 女生性感内裤真人,穿戴方法视频| 国产精品乱码一区二三区的特点 | 欧美最黄视频在线播放免费 | 亚洲精品国产精品久久久不卡| 久久久水蜜桃国产精品网| 老汉色av国产亚洲站长工具| 国产伦人伦偷精品视频| 啦啦啦 在线观看视频| 久久人妻熟女aⅴ| 日韩高清综合在线| 久久久久九九精品影院| 亚洲一区高清亚洲精品| 制服诱惑二区| 国产极品粉嫩免费观看在线| 色婷婷av一区二区三区视频| 久久久久精品国产欧美久久久| 精品久久蜜臀av无| 多毛熟女@视频| 亚洲精品久久午夜乱码| 久久久久久久久中文| 757午夜福利合集在线观看| 日本黄色视频三级网站网址| 午夜激情av网站| 一进一出抽搐gif免费好疼 | 久久久久久久久中文| 中文字幕精品免费在线观看视频| 久久精品亚洲熟妇少妇任你| 人人澡人人妻人| 欧美日韩福利视频一区二区| av网站免费在线观看视频| 亚洲精品国产色婷婷电影| 午夜福利影视在线免费观看| 高清在线国产一区| 日韩中文字幕欧美一区二区| 欧美日本亚洲视频在线播放| 午夜福利一区二区在线看| 欧美黑人欧美精品刺激| 1024香蕉在线观看| 极品教师在线免费播放| 又大又爽又粗| 精品国产一区二区久久| 欧美中文日本在线观看视频| 国产精品永久免费网站| av福利片在线| 日本 av在线| 1024香蕉在线观看| 国产精品av久久久久免费| 美女高潮到喷水免费观看| 51午夜福利影视在线观看| 亚洲国产精品一区二区三区在线| 亚洲性夜色夜夜综合| 亚洲五月天丁香| 国产av精品麻豆| 亚洲五月天丁香| 欧美中文日本在线观看视频| 黑人欧美特级aaaaaa片| 国产伦一二天堂av在线观看| 国产成人欧美| 人人妻,人人澡人人爽秒播| 午夜福利欧美成人| 午夜福利免费观看在线| 亚洲av电影在线进入| 又大又爽又粗| 母亲3免费完整高清在线观看| 18禁黄网站禁片午夜丰满| 国产精品一区二区免费欧美| 波多野结衣一区麻豆| 欧美另类亚洲清纯唯美| 国产aⅴ精品一区二区三区波| 热re99久久国产66热| 久久精品91蜜桃| 美女午夜性视频免费| 老汉色av国产亚洲站长工具| 国产成人精品无人区| 最近最新免费中文字幕在线| 亚洲欧美激情综合另类| 欧美国产精品va在线观看不卡| 日韩国内少妇激情av| 国产精品 欧美亚洲| 精品国产国语对白av| 人人澡人人妻人| 国产乱人伦免费视频| 欧美激情久久久久久爽电影 | 欧美日韩瑟瑟在线播放| 一级黄色大片毛片| 国产精品国产av在线观看| а√天堂www在线а√下载| 午夜福利在线观看吧| 亚洲,欧美精品.| 亚洲第一av免费看| 国产欧美日韩精品亚洲av| 日韩欧美免费精品| 性少妇av在线| 亚洲国产中文字幕在线视频| 99精品久久久久人妻精品| 亚洲欧洲精品一区二区精品久久久| 丰满迷人的少妇在线观看| 色精品久久人妻99蜜桃| 精品久久蜜臀av无| 国产视频一区二区在线看| 露出奶头的视频| 法律面前人人平等表现在哪些方面| 老汉色∧v一级毛片| 国产精品一区二区三区四区久久 | 国产欧美日韩一区二区三| 国产午夜精品久久久久久| 午夜免费激情av| 成熟少妇高潮喷水视频| 麻豆一二三区av精品| 欧美性长视频在线观看| 亚洲精品国产区一区二| av在线播放免费不卡| 久久国产精品男人的天堂亚洲| 国产精品免费一区二区三区在线| 俄罗斯特黄特色一大片| 欧美日韩一级在线毛片| 亚洲黑人精品在线| 欧美成狂野欧美在线观看| 麻豆国产av国片精品| 色婷婷av一区二区三区视频| 夜夜夜夜夜久久久久| 免费观看精品视频网站| 色婷婷久久久亚洲欧美| 国产精品影院久久| 十八禁网站免费在线| 国产色视频综合| 一区二区三区精品91| 老司机亚洲免费影院| 久久精品国产综合久久久| 美女国产高潮福利片在线看| 国产日韩一区二区三区精品不卡| www日本在线高清视频| 韩国精品一区二区三区| 欧美老熟妇乱子伦牲交| 天堂动漫精品| 久久99一区二区三区| 欧美黄色淫秽网站| 18禁黄网站禁片午夜丰满| 中文字幕色久视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品一卡2卡三卡4卡5卡| 免费看a级黄色片| 天堂√8在线中文| 久久午夜亚洲精品久久| 黄色怎么调成土黄色| 精品国产一区二区久久| 欧美乱码精品一区二区三区| 久热这里只有精品99| 午夜老司机福利片| 亚洲av片天天在线观看| 亚洲精品成人av观看孕妇| 国产精品香港三级国产av潘金莲| 在线看a的网站| 可以免费在线观看a视频的电影网站| 黄色毛片三级朝国网站| 国产男靠女视频免费网站| 久久精品国产清高在天天线| 啪啪无遮挡十八禁网站| av网站在线播放免费| 在线十欧美十亚洲十日本专区| 老司机靠b影院| 三上悠亚av全集在线观看| 精品一品国产午夜福利视频| 亚洲激情在线av| 夜夜爽天天搞| 熟女少妇亚洲综合色aaa.| 我的亚洲天堂| 国产av一区在线观看免费| 久久影院123| 日日夜夜操网爽| 日韩精品青青久久久久久| 国产一区二区激情短视频| 久久久久久久精品吃奶| 在线观看午夜福利视频| 夫妻午夜视频| 国产亚洲精品久久久久久毛片| 亚洲一区二区三区不卡视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日本黄色日本黄色录像| 成人手机av| 亚洲午夜精品一区,二区,三区| 日韩 欧美 亚洲 中文字幕| 日本a在线网址| 欧洲精品卡2卡3卡4卡5卡区| 欧美性长视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 成人三级做爰电影| 在线av久久热| 亚洲欧洲精品一区二区精品久久久| 国产精品99久久99久久久不卡| 啦啦啦免费观看视频1| 中文字幕人妻丝袜制服| aaaaa片日本免费| 99国产精品99久久久久| 欧美+亚洲+日韩+国产| 免费高清视频大片| 88av欧美| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美精品济南到| 色精品久久人妻99蜜桃| 成人亚洲精品一区在线观看| 久久久国产欧美日韩av| 亚洲一区高清亚洲精品| e午夜精品久久久久久久| 国产亚洲精品综合一区在线观看 | 老熟妇乱子伦视频在线观看| 国产一区二区三区视频了| 日韩精品免费视频一区二区三区| 国产精品 国内视频| 十分钟在线观看高清视频www| 亚洲av五月六月丁香网| 久久久国产欧美日韩av| 夜夜爽天天搞| 人人澡人人妻人| 国产伦一二天堂av在线观看| 老汉色∧v一级毛片| 亚洲美女黄片视频| 久久人人精品亚洲av| 亚洲人成77777在线视频| 午夜福利免费观看在线| 国产麻豆69| 国产成+人综合+亚洲专区| 国产野战对白在线观看| 欧美日韩亚洲综合一区二区三区_| 中文字幕人妻熟女乱码| 色播在线永久视频| 欧美av亚洲av综合av国产av| 51午夜福利影视在线观看| 另类亚洲欧美激情| 欧美人与性动交α欧美精品济南到| 久久精品aⅴ一区二区三区四区| 国产精品日韩av在线免费观看 | а√天堂www在线а√下载| 欧美中文日本在线观看视频| 老司机在亚洲福利影院| 99国产精品免费福利视频| 一边摸一边抽搐一进一小说| 亚洲男人的天堂狠狠| 女同久久另类99精品国产91| 日韩欧美一区二区三区在线观看| 男人操女人黄网站| 看片在线看免费视频| 久久精品亚洲熟妇少妇任你| 99久久人妻综合| 天堂中文最新版在线下载| 级片在线观看| 国产成人欧美| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美一区二区综合| 欧美日韩中文字幕国产精品一区二区三区 | 桃色一区二区三区在线观看| 亚洲精品久久午夜乱码| 欧美av亚洲av综合av国产av| 久久久久九九精品影院| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 日韩高清综合在线| 夜夜躁狠狠躁天天躁| 亚洲精品一区av在线观看| 日本免费a在线| 正在播放国产对白刺激| 97碰自拍视频| 18禁裸乳无遮挡免费网站照片 | 日韩免费高清中文字幕av| 亚洲中文日韩欧美视频| 村上凉子中文字幕在线| 性色av乱码一区二区三区2| 男女高潮啪啪啪动态图| 精品电影一区二区在线| 天堂中文最新版在线下载| 久久午夜亚洲精品久久| 久久草成人影院| 中文字幕高清在线视频| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 1024香蕉在线观看|