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

    基于SSR分子標(biāo)記的部分烤煙種質(zhì)資源遺傳多樣性研究

    2020-05-25 02:47丁永亮陳仁霄苑舉民管成偉何寬信張興偉李卓張啟明
    中國煙草科學(xué) 2020年2期
    關(guān)鍵詞:遺傳多樣性種質(zhì)資源煙草

    丁永亮 陳仁霄 苑舉民 管成偉 何寬信 張興偉 李卓 張啟明

    摘 ?要:為準(zhǔn)確評價部分煙草種質(zhì)資源的遺傳多樣性,利用42對已發(fā)表多態(tài)性較好的煙草SSR標(biāo)記引物對38份烤煙種質(zhì)資源進(jìn)行分析,篩選出22對引物在這些種質(zhì)間存在多態(tài)性,平均等位位點(diǎn)數(shù)為3.68個,Neis多樣性指數(shù)平均為0.451,Shannon信息指數(shù)平均達(dá)0.789。UPGMA聚類分析表明,紅花大金元有明顯特異性。遺傳結(jié)構(gòu)分析顯示,參試材料可分為2個類群,類群1和類群2分別包含16份和22份種質(zhì)材料,且類群1變異范圍大于類群2,類群2材料多表現(xiàn)為自然株高在210 cm以下,葉數(shù)20~25,莖圍10~11 cm,田間赤星病較嚴(yán)重,青枯病發(fā)病相對較輕。研究結(jié)果為提升SSR分子標(biāo)記篩選效率和部分烤煙種質(zhì)創(chuàng)新應(yīng)用提供理論依據(jù)。

    關(guān)鍵詞:煙草;種質(zhì)資源;SSR分子標(biāo)記;遺傳多樣性;遺傳結(jié)構(gòu)

    Genetic Diversity of Flue-cured Tobacco Germplasm Resources Based on SSR Molecular Markers

    DING Yongliang1, CHEN Renxiao1, YUAN Jumin1, GUAN Chengwei1, HE Kuanxin1,

    ZHANG Xingwei2, LI Zhuo1, ZHANG Qiming1*

    (1. Tobacco Science Institute of Jiangxi Province, Nanchang 330009, China; 2. Institute of Tobacco Research of CAAS, Qingdao 266101, China)

    Abstract: In order to accurately evaluate the genetic diversity of introduced tobacco germplasm resources, 38 tobacco germplasm resources were analyzed with 42 SSR primers with good polymorphism. A total of 22 pairs of polymorphic primers were selected based on higher identification ability. The average number of alleles was 3.68. The average diversity index of Nei was 0.451, and the average Shannon information index was 0.789. The UPGMA cluster analysis showed the specificity of Honghuadajinyuan. The genetic structure analysis showed that these 38 accessions were classified into two groups (group 1 and group 2), which included 16 and 22 accessions, respectively. The range of variation of group 1 was larger than that of group 2. Most accessions in group 2 showed less than 210 cm in plant height, more than 20 in the number of leaves, 10-11 cm in the stem circumference. Furmore,the accessions in group 2 were susceptible to brown spot disease but resistant to bacterial wilt. The results can provide theoretical basis for improving SSR molecular marker screening efficiency and innovative utilization of some flue-cured tobacco germplasm.

    Keywords: tobacco; germplasm resources; SSR molecular marker; genetic diversity; genetic structure

    煙草是我國重要經(jīng)濟(jì)作物,在國民經(jīng)濟(jì)發(fā)展中發(fā)揮著重要作用。優(yōu)質(zhì)煙葉是卷煙工業(yè)的重要原料,而品種是彰顯煙葉質(zhì)量特色的關(guān)鍵因素之一,對煙葉品質(zhì)起主導(dǎo)作用[1]。種質(zhì)資源是選育優(yōu)良煙草品種的重要基礎(chǔ)[2],烤煙遺傳背景狹窄是煙草育種發(fā)展中的突出問題[3-4]。因此,建立遺傳資源豐富的種質(zhì)庫尤為重要。分子標(biāo)記技術(shù)目前已廣泛應(yīng)用于煙草的遺傳多樣性和親緣關(guān)系研究,而且在煙草化學(xué)成分、普通農(nóng)藝性狀和抗病性等性狀連鎖基因定位中也取得了一定進(jìn)展[5]。在基于PCR技術(shù)的分子標(biāo)記方法中,SSR(Simple Sequence Repeats)標(biāo)記具有多態(tài)性高、操作簡單、穩(wěn)定性好等諸多優(yōu)點(diǎn),在許多作物的遺傳多樣性研究中廣泛應(yīng)用[6]。對于煙草,BINDLER等[7-8]首次報(bào)道了其SSR標(biāo)記,構(gòu)建了包括2317個SSR標(biāo)記和2363個位點(diǎn),并覆蓋24個連鎖群的高密度遺傳圖譜,總遺傳長度為3267 cM,標(biāo)記間的平均距離小于1.5 cM。此后,SSR分子標(biāo)記在煙草種質(zhì)間遺傳多樣性研究中大量應(yīng)用。葉蘭欽等[9]利用92對煙草SSR標(biāo)記分析了云南省13個煙草主栽品種的遺傳多樣性,檢測到20對SSR多態(tài)性引物。陳夏曄等[10]從960對引物中篩選出34對SSR多態(tài)性引物,分析73份煙草抗黑脛病種質(zhì),認(rèn)為其遺傳背景較狹窄。尹國英等[6]從278對引物中篩選出32對多態(tài)性較好的SSR引物并對140份煙草種質(zhì)進(jìn)行分析,結(jié)果表明晾曬煙具有豐富的遺傳多樣性。本研究直接利用前人研究中多態(tài)性較高的SSR分子標(biāo)記,對江西省自存和引進(jìn)的烤煙種質(zhì)資源開展遺傳多樣性分析,旨在為煙草新品種的創(chuàng)育和分子標(biāo)記輔助育種提供理論指導(dǎo)。

    1 ?材料與方法

    1.1 ?試驗(yàn)材料

    試驗(yàn)所用烤煙種質(zhì)資源共計(jì)38份(表1),其中18份由國家農(nóng)作物種質(zhì)資源平臺煙草種質(zhì)資源子平臺提供;4份為江西省內(nèi)自育品系和收集的突變體,其余均從國內(nèi)各省引進(jìn)。

    1.2 ?試驗(yàn)地點(diǎn)和方法

    1.2.1 ?試驗(yàn)地點(diǎn) ?所有材料于2018年3月10日移栽到江西省煙草科學(xué)研究所樂安試驗(yàn)田內(nèi),土壤為

    砂壤土,無嚴(yán)重病蟲害發(fā)生。土壤肥力狀況:pH值5.3,有機(jī)質(zhì)含量19.1 g/kg,速效氮含量90.3 mg/kg,速效磷含量24.5 mg/kg,速效鉀含量99.5 mg/kg,氯離子含量0.12 mg/kg。

    1.2.2 ?農(nóng)藝性狀及病害調(diào)查 ?采用隨機(jī)區(qū)組設(shè)計(jì),每份材料種植40株,設(shè)3次重復(fù),行株距1.20 m×

    0.50 m,田間管理和調(diào)制按照當(dāng)?shù)亟y(tǒng)一模式進(jìn)行。按照YC/T 142—2010煙草農(nóng)藝性狀調(diào)查測量方法調(diào)查各種質(zhì)單株著生葉數(shù)、自然株高、莖圍和節(jié)距等農(nóng)藝性狀,每個小區(qū)調(diào)查10株。赤星病和青枯病參照GB/T 23222—2008煙草病蟲害分級及調(diào)查方法進(jìn)行調(diào)查,其中赤星病每個小區(qū)調(diào)查50片葉,青枯病調(diào)查全部煙株。

    1.2.3 ?煙草DNA提取 ?以每份種質(zhì)嫩葉為材料,利用磁珠法植物DNA提取試劑盒KFR-804496(Thermo Fisher公司)提取總DNA,用超微量分光光度計(jì)NANODROP 2000(Thermo Fisher公司)檢測DNA濃度和純度,將其稀釋至約50 ng/μL,?20 ℃保存?zhèn)溆谩?/p>

    1.2.4 ?PCR擴(kuò)增與檢測 ?采用已發(fā)表多態(tài)性較好的42對煙草SSR標(biāo)記引物[6,11-15],引物由英濰捷基(上海)貿(mào)易有限公司合成。PCR擴(kuò)增體系參照Ex Taq試劑盒(購自TAKARA公司)說明配制。擴(kuò)增程序:94 ℃預(yù)變性5 min,PCR擴(kuò)增35個循環(huán),每循環(huán)94 ℃變性30 s,55 ℃(視引物而變)退火30 s,72 ℃延伸30 s,最后72 ℃延伸10 min。利用10%非變性聚丙烯酰胺凝膠電泳檢測PCR產(chǎn)物,在180 V恒壓下電泳3 h,用銀染法顯影。

    1.2.5 ?數(shù)據(jù)處理 ?利用Excel 2013軟件統(tǒng)計(jì)農(nóng)藝性狀和病害病情指數(shù)。對于SSR帶型,在相同遷移率的位置上,有帶記錄為“1”,無帶記錄為“0”。利用Dataformater軟件對0、1矩陣進(jìn)行格式轉(zhuǎn)換,并過濾無多態(tài)位點(diǎn)。利用Popgene軟件計(jì)算基因多樣性和Shannon信息指數(shù)。利用NTSYS-pc 2.10軟件構(gòu)建UPGMA(非加權(quán)組平均法)聚類圖,并利用DCENTER和EIGEN模塊對供試煙草材料進(jìn)行主成分分析(PCA)。利用Structure 2.3軟件對煙草材料進(jìn)行群體結(jié)構(gòu)分析,潛在的分組數(shù)值K的范圍設(shè)置在1~10。根據(jù)原始結(jié)果,每一個K值將獨(dú)立運(yùn)行10次再取平均值,其中burn-in period 10 000 步,MCMC重復(fù)100 000,ΔK參考EVANNO等[16]的算法用來決定最佳組群數(shù)量。

    2 ?結(jié) ?果

    2.1 ?SSR引物多態(tài)性分析

    PCR擴(kuò)增產(chǎn)物通過10%非變性聚丙烯酰胺凝膠檢測結(jié)果表明,篩選的42對SSR標(biāo)記中有27對檢測出等位變異。利用Dataformater軟件對0、1矩陣進(jìn)行格式轉(zhuǎn)換,并過濾無多態(tài)位點(diǎn),篩選出22對具有多態(tài)位點(diǎn)的SSR標(biāo)記(表2)。觀測等位基因數(shù)變幅在2~5之間,平均為3.68個。有效等位基因數(shù)變幅在1.174~2.885之間,平均為1.94個,Shannon信息指數(shù)平均為0.789,Neis多樣性指數(shù)變幅在0.148~0.653之間,平均為0.451,其中PT20213多樣性指數(shù)最高,表明具有較高的鑒別能力(圖1)。

    2.2 ?種質(zhì)資源聚類分析

    為確定38份種質(zhì)材料在基因組上的遺傳關(guān)系,利用NTSYS軟件按UPGMA方法進(jìn)行聚類分析,構(gòu)建了參試材料的親緣關(guān)系圖。由圖2可知,22對多態(tài)性引物可將36份種質(zhì)材料區(qū)分歸類,聚類結(jié)果與已知種質(zhì)親緣關(guān)系大致相近,美國引進(jìn)種質(zhì)和部分突變體材料遺傳距離較近,國內(nèi)各種質(zhì)間距離較近。利用DCENTER和EIGEN對供試材料進(jìn)行主成分分析(圖3),結(jié)果表明,38份種質(zhì)材料大致可分為2個類群,且兩個類群的離散程度存在較大差異。

    2.3 ?群體遺傳結(jié)構(gòu)分析

    利用Structure軟件對煙草種質(zhì)材料進(jìn)行群體結(jié)構(gòu)分析(圖4)得出,群體等位變異頻率特征類型數(shù)K呈逐漸增大的趨勢,通過EVANNO等[16]ΔK值來確定K值。ΔK值在K=2時值最大,表明等位變異頻率特征類型數(shù)K=2時其模型后驗(yàn)概率最大,由此可將38份種質(zhì)材料分為2個類群。煙草38個種質(zhì)材料的群體遺傳結(jié)構(gòu)如圖5,其中類群1(POP 1)包括G80、K326、Coker86、CV70、CV87、中煙98、中煙102、云煙97、云煙99、云煙110、HY06、粵煙98、FJ109、閩煙9號、紅花大金元、新豐變等16個種質(zhì)材料,類群2(POP2)包括剩余的22份種質(zhì)材料。遺傳結(jié)構(gòu)分析與聚類分析在云煙110、云煙97和Coker86等材料的群體劃分結(jié)果上存在差異,與主成分分析在云煙110、NC2514、豫煙10號等材料的群體劃分結(jié)果上存在差異。

    2.4 ?兩個遺傳群體農(nóng)藝性狀與抗病性變異分析

    如圖6所示,群體遺傳結(jié)構(gòu)分析所得兩個群體的農(nóng)藝性狀和田間主要病害發(fā)病情況存在差異。POP1在各農(nóng)藝性狀統(tǒng)計(jì)范圍內(nèi)均有分布,POP2自然株高均在210 cm以下,99.91%的種質(zhì)材料自然株高不超過185 cm,著生葉片17.6片以上,莖圍在10.1~11.0 cm范圍內(nèi)種質(zhì)數(shù)占比高達(dá)86.36%,節(jié)距集中在3.6~4.9 cm。在赤星病方面,POP1有12.5%的個體發(fā)病指數(shù)在10以內(nèi),POP2發(fā)病指數(shù)均在10.1以上。在青枯病方面,兩個類群病情指數(shù)主要處于0~5.0之間,POP1和POP2在此范圍內(nèi)的株數(shù)分別占比56.25%和81.82%,其中POP2發(fā)病指數(shù)均在15.0以下。由此可見,POP1的農(nóng)藝性狀變異范圍大于POP2,POP2材料多表現(xiàn)為株高在210 cm以下,葉片數(shù)20~25,莖圍10~11 cm,田間赤星病較嚴(yán)重,青枯病發(fā)病相對較輕。

    3 ?討 ?論

    篩選出多態(tài)性較好的SSR分子標(biāo)記是開展煙草SSR相關(guān)研究的基礎(chǔ)。張銘真等[17]從250對引物中篩選出64對SSR引物,分析了81份煙草種質(zhì)資源的遺傳多樣性,平均等位位點(diǎn)數(shù)為4.2個。馬冰等[18]從2317對引物中篩選出24對核心引物對20份烤煙種質(zhì)進(jìn)行遺傳多樣性分析,平均等位位點(diǎn)數(shù)為3.54個。鞠馥竹等[19]從580對引物中篩選出43對多態(tài)性引物,平均等位位點(diǎn)數(shù)為3.39個。為充分利用前人研究成果,本研究利用42對已發(fā)表的多態(tài)性較好的SSR標(biāo)記引物對38份煙草種質(zhì)材料進(jìn)行遺傳多樣性分析,篩選到22對引物在這些種質(zhì)間存在多態(tài)性,檢測到81條多態(tài)性擴(kuò)增條帶,平均等位位點(diǎn)數(shù)為3.68個,與前人研究相當(dāng),但較大程度減少了標(biāo)記篩選工作量,提升了效率。Neis多樣性指數(shù)平均為0.451,Shannon信息指數(shù)平均達(dá)0.789,說明這些引物總體上有較好的鑒別能力,種質(zhì)材料具有較高的遺傳多樣性。

    本研究所用種質(zhì)材料包括國外引進(jìn)種質(zhì)、國內(nèi)地方種質(zhì)和本省自有突變體等,聚類分析和主成分分析發(fā)現(xiàn)紅花大金元與大部分種質(zhì)材料遺傳距離大,可以考慮將其作為親本材料拓寬遺傳背景,發(fā)掘優(yōu)異性狀。聚類及遺傳結(jié)構(gòu)分析顯示美國引進(jìn)的種質(zhì)材料遺傳距離近,遺傳背景相似,與MOON等[20]研究結(jié)果相一致。聚類分析未能將參試的38份種質(zhì)材料劃分成明顯的類群,主成分分析初步顯示38份材料可劃分為2個類群,與類群2相比,類群1種質(zhì)間變異和離散程度更大,這與類群遺傳結(jié)構(gòu)分析的類群間農(nóng)藝性狀及抗病性變異分析相一致。但是兩個分類結(jié)果在云煙110、云煙97和Coker86 3個材料的類群劃分上存在差異,造成這種差異的主要原因是兩者分類依據(jù)不同。聚類分析依據(jù)的是烤煙材料間的遺傳距離,反映其親緣關(guān)系的遠(yuǎn)近;而群體遺傳結(jié)構(gòu)分類依據(jù)服從Hardy-Weinberg平衡的群體數(shù)目,基于數(shù)學(xué)模型之上,并計(jì)算各供試材料間相應(yīng)的Q值[21]。因此,本試驗(yàn)利用Structure軟件將供試烤煙材料分為2個類群,具有較高的準(zhǔn)確性。此外,關(guān)聯(lián)分析已廣泛應(yīng)用于大麥[22]、玉米[23]、小麥[24]、水稻[25]、花生[26]等作物的遺傳育種研究,煙草與農(nóng)藝性狀[21]、抗病性[27]、致香物質(zhì)[28]等方面的關(guān)聯(lián)分析也大量出現(xiàn),本研究可進(jìn)一步在農(nóng)藝性狀、抗病性、化學(xué)成分等方面開展多年重復(fù)試驗(yàn),尋找與性狀緊密相關(guān)的標(biāo)記,以用于分子標(biāo)記輔助育種研究。

    4 ?結(jié) ?論

    利用已發(fā)表多態(tài)性較好的42對煙草SSR標(biāo)記引物對38份煙草種質(zhì)材料進(jìn)行分析,篩選到22對引物在這些種質(zhì)間存在多態(tài)性,平均等位位點(diǎn)數(shù)為3.68個。參試的38份煙草種質(zhì)具有較高的遺傳多樣性,UPGMA聚類分析顯示紅花大金元具有明顯特異性。遺傳結(jié)構(gòu)分析顯示,參試材料可分為2個類群,在農(nóng)藝性狀和抗病性方面,類群1變異范圍大于類群2,類群2材料多表現(xiàn)為株高在210 cm以下,葉數(shù)20~25,莖圍10~11 cm,田間赤星病較嚴(yán)重,青枯病發(fā)病相對較輕。本研究結(jié)果為提升SSR分子標(biāo)記篩選效率和部分烤煙種質(zhì)創(chuàng)新應(yīng)用提供理論依據(jù),并為進(jìn)一步與農(nóng)藝性狀、抗病性、化學(xué)成分等方面開展關(guān)聯(lián)分析奠定了基礎(chǔ)。

    參考文獻(xiàn)

    [1]中國農(nóng)業(yè)科學(xué)院煙草研究所. 中國煙草栽培學(xué)[M]. 上海:上??茖W(xué)技術(shù)出版社,2005.

    Institute of Tobacco Research of CAAS. Tobacco cultivation in China[M]. Shanghai: Shanghai Science And Technology Press, 2005.

    [2]張艇. 煙草種質(zhì)資源的遺傳多樣性研究[J]. 農(nóng)技服務(wù),2010,27(1):101-103.

    ZHANG T. Research on genetic diversity of tobacco germplasm resource[J]. Agricultural Technology Services, 2010, 27(1): 101-103.

    [3]孫計(jì)平,吳照輝,李雪君,等. 21世紀(jì)中國烤煙種植區(qū)域及主栽品種變化分析[J]. 中國煙草科學(xué),2016,37(3):86-92.

    SUN J P, WU Z H, LI X J, et al. Analysis of regional variation and major varieties of flue-cured tobacco planted in China in the 21 century[J]. Chinese Tobacco Science, 2016, 37(3): 86-92.

    [4]陳雅瓊,王衛(wèi)峰,丁安明,等. 基于SSR熒光標(biāo)記分析我國主要審(認(rèn))定烤煙及白肋煙品種遺傳多樣性[J]. 中國煙草學(xué)報(bào),2013,19(2):98-105.

    CHEN Y Q, WANG W F, DING A M, et al. Analysis of genetic diversity in national listed tobacco varieties using SSR fluorescent marker data[J]. Acta Tabacaria Sinica, 2013, 19(2): 98-105.

    [5]羅昭標(biāo),陳歡,毛華,等. 分子標(biāo)記在煙草性狀連鎖基因定位中的應(yīng)用進(jìn)展[J]. 生物技術(shù)通訊,2014,25(1):139-142.

    LUO Z B, CHEN H, MAO H, et al. Progress of the application of molecular markers in gene mapping of tobacco characters [J]. Letters In Biotechnology, 2014, 25(1): 139-142.

    [6]尹國英,楊小燕,何其波,等. 鑒定煙草種質(zhì)資源SSR核心引物篩選和驗(yàn)證[J]. 植物遺傳資源學(xué)報(bào),2013,14(5):960-965.

    YIN G Y, YANG X Y, HE Q B, et al. Screen and identification of SSR core primers for tobacco germplasm[J]. Journal of Plant Genetic Resources, 2013, 14(5): 960-965.

    [7]BINDLER G, VAN DER HOEVEN R, GUNDUZ I, et al. A microsatellite marker based linkage map of tobacco[J]. Theoretical and Applied Genetics, 2007, 114(2): 341-349.

    [8]BINDLER G, PLIESKE J, BAKAHER N, et al. A high density genetic map of tobacco (Nicotiana tabacum L.) obtained from large scale microsatellite marker development[J]. Theoretical and Applied Genetics, 2011, 123(2): 219.

    [9]葉蘭欽,辛明明,杜金昆,等. SSR標(biāo)記應(yīng)用于煙草品種遺傳多樣性研究[J]. 中國農(nóng)學(xué)通報(bào),2009,25(1):56-62.

    YE L Q, XIN M M, DU J K, et al. Genetic diversity revealed by SSR markers of the tobacco cultivars[J]. Chinese Agricultural Science Bulletin, 2009, 25(1): 56-62.

    [10]陳夏曄,劉艷華,張興偉,等. 抗黑脛病煙草種質(zhì)的遺傳多樣性分析[J]. 中國農(nóng)學(xué)通報(bào),2014,30(7):58-63

    CHEN X Y, LIU Y H, ZHANG X W, et al. Genetic diversity analysis on resistant tobacco germplasm to black shank disease[J]. Chinese Agricultural Science Bulletin, 2014, 30(7): 58-63.

    [11]賴瑞強(qiáng),李榮華,夏巖石,等. 煙草種質(zhì)的SSR標(biāo)記遺傳多樣性及青枯病抗性的關(guān)聯(lián)分析[J]. 中國煙草學(xué)報(bào),2018,24(6):67-77.

    LAI R Q, LI R H, XIA Y S, et al. SSR marker-based genetic diversity analysis of tobacco germplasm and association analysis with resistance to bacterial wilt[J]. Acta Tabacaria Sinica, 2018, 24(6): 67-77.

    [12]張興偉,任民,劉艷華,等. 煙草核心引物的篩選[C]//中國煙草2013年學(xué)術(shù)年會論文集. 北京:中國煙草學(xué)會,2013:109-114.

    ZHANG X W, REN M, LIU Y H, et al. Screening of SSR core primers with polymorphism on 12 tobacco accessions[C]//Proceedings of the 2013 China tobacco annual conference. Beijing: China Tobacco Society, 2013: 109-114.

    [13]向世鵬,胡日生,周向平,等. 煙草種質(zhì)資源黑脛病抗性鑒定及親緣關(guān)系SSR分析[J]. 華北農(nóng)學(xué)報(bào),2016,31(1):156-161.

    XIANG S P, HU R S, ZHOU X P, et al. Identification of resistance to black shank disease of tobacco germplasm resources and analysis of genetic relationship of SSR[J]. Acta Agriculturae Boreali-sinica, 2016, 31(1): 156-161.

    [14]陳芳,徐世曉,李曉輝,等. 基于SSR標(biāo)記的80份煙草種質(zhì)指紋圖譜的構(gòu)建及遺傳多樣性分析[J]. 作物雜志,2019(1):22-31.

    CHEN F, XU S X, LI X H, et al. Construction of molecular fingerprinting and analysis of genetic diversity for 80 tobacco (Nicotiana tabacum) germplasms[J]. Crops, 2019(1): 22-31.

    [15]張銘真,邢雪霞,李曉輝,等. 煙草種質(zhì)資源遺傳多樣性與群體結(jié)構(gòu)分析[J]. 中國煙草科學(xué),2016,37(4):42-47.

    ZHANG M Z, XING X X, LI X H, et al. Analysis of genetic diversity and population structure of germplasm resources in Nicotiana tobacum[J]. Chinese Tobacco Science, 2016, 37(4): 42-47.

    [16]EVANNO G, REGNAUT S, GOUDET J. Detecting the number of clusters of individuals using the software structure: a simulation study[J]. Molecular Ecology, 2005, 14(8): 2611-2620.

    [17]張銘真,李曉輝,王袁,等. 81份煙草種質(zhì)資源遺傳多樣性分析[J]. 江西農(nóng)業(yè)學(xué)報(bào),2017,29(1):62-68.

    ZHANG M Z, LI X H, WANG Y, et al. Analysis of genetic diversity of 81 tobacco germplasm resources[J], Acta Agricultural Jiangxi, 2017, 29(1): 62-68.

    [18]馬冰,代帥帥,程亞增,等. 烤煙種質(zhì)資源SSR核心引物的篩選及驗(yàn)證[J]. 中國煙草科學(xué),2016,37(5):1-5,9.

    MA B, DAI S S, CHENG Y Z, et al. Screen and identification of SSR core primers for flue-cured tobacco germplasm[J]. Chinese Tobacco Science, 2016, 37(5): 1-5, 9.

    [19]鞠馥竹,趙文濤,劉元德,等. 不同產(chǎn)區(qū)晾曬煙資源多樣性的鑒定與評價[J]. 中國煙草科學(xué),2019,40(2):8-15.

    JU F Z, ZHAO W T, LIU Y D, et al. Identification and evaluation of the diversity of air/sun-cured tobacco germplasm resources in different producing areas[J]. Chinese Tobacco Science, 2019, 40(2): 8-15.

    [20]MOON H S, NIFONG J M, NICHOLSON J S, et al. Microsatellite-based analysis of tobacco (Nicotiana tabacum L.) genetic resources[J]. Crop Science, 2009, 49: 2149-2159.

    [21]童治軍,陳學(xué)軍,方敦煌,等. 231份烤煙種質(zhì)資源SSR標(biāo)記遺傳多樣性及其與農(nóng)藝性狀和化學(xué)成分的關(guān)聯(lián)分析[J]. 中國煙草學(xué)報(bào),2017,23(5):31-61.

    TONG Z J, CHEN X J, FANG D H, et al. SSR marker-based analyses on genetic diversity and relevant variations of agronomic traits and chemical composition of 231 flue-cured tobacco germplasm resource[J]. Acta Tabacaria Sinica, 2017, 23(5): 31-61.

    [22]IVANDIC V, HACKETT C A, NEVO E, et al. Analysis of simple sequence repeats (SSR) in wild barley from the fertile crescent: associations with ecology, geography and flowering time[J]. Plant Molecular Biology, 2002, 48: 511-527.

    [23]FLINT-GARCIA S A, THUILLET A C, YU J, et al. Maize association population: a high-resolution platform for quantitative trait locus dissection[J]. The Plant Journal, 2005, 44(6): 1054-1064.

    [24]ROY J K, BANDOPADHYAY R, RUSTGI S, et al. Association analysis of agronomical important traits using SSR, SAMPL and AFLP markers in bread wheat[J]. Current Science, 2006, 90(5): 683-689.

    [25]AGRAMA H A, EIZENGA G C, YAN W. Association mapping of yield and its components in rice cultivars[J]. Molecular Breeding, 2007, 19(4): 341-356.

    [26]李東霞,石鵬,楊偉波,等. 利用SSR分子標(biāo)記分析海南與國內(nèi)外花生種質(zhì)資源的遺傳多樣性[J]. 廣東農(nóng)業(yè)科學(xué),2015,42(20):118-124.

    LI D X, SHI P, YANG W B, et al. Analysis of genetic diversity in peanut(Arachis hypogaea L.)germplasm resources from Hainan and other regions by using SSR markers[J]. Guangdong Agricultural Sciences, 2015, 42(20): 118-124.

    [27]叢鑫,劉艷華,戴培剛,等. 抗PVY煙草種質(zhì)的遺傳多樣性分析[J]. 植物遺傳資源學(xué)報(bào),2014,15(3):679-684.

    CONG X, LIU Y H, DAI P G, et al. Genetic diversity analysis of tobacco germplasm resistance to PVY[J]. Journal of Plant Genetic Resources, 2014, 15(3): 679-684.

    [28]任民,張長靜,蔣彩虹,等. 基于高密度SSR連鎖群的煙草致香物質(zhì)關(guān)聯(lián)分析[J]. 中國煙草學(xué)報(bào),2014,20(4):88-93.

    REN M, ZHANG C J, JIANG C H, et al. Association analysis of tobacco aroma constituents based on high density SSR linkage group[J]. Acta Tabacaria Sinica, 2014, 20(4): 88-93.

    基金項(xiàng)目:江西省煙草公司科技項(xiàng)目“煙草種質(zhì)資源的收集引進(jìn)與種質(zhì)庫建立”(201401001)

    作者簡介:丁永亮(1990-),男,農(nóng)藝師,碩士,從事煙草育種及生物技術(shù)研究。E-mail:yld9001@163.com

    *通信作者,E-mail:zhangqiming1979@163.com

    收稿日期:2019-10-23 ? ? ? ? ? ? ? ? ? ? 修回日期:2020-02-14

    猜你喜歡
    遺傳多樣性種質(zhì)資源煙草
    煙草產(chǎn)品展廳設(shè)計(jì)
    綠肥作物紫云英研究進(jìn)展
    大白菜種質(zhì)資源抗根腫病基因CRa和CRb的分子標(biāo)記鑒定與分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    玉米種質(zhì)資源抗旱性鑒定研究進(jìn)展
    一地?zé)煵?,半世孤?dú)
    全國首次青少年煙草調(diào)查
    天堂动漫精品| 桃色一区二区三区在线观看| 久久久精品大字幕| av.在线天堂| 最新中文字幕久久久久| 久久久久久久午夜电影| 九九热线精品视视频播放| 亚洲av美国av| 亚洲av中文av极速乱 | 啪啪无遮挡十八禁网站| 欧美激情在线99| 女人十人毛片免费观看3o分钟| 美女xxoo啪啪120秒动态图| 亚洲内射少妇av| 天堂动漫精品| 免费人成视频x8x8入口观看| 男插女下体视频免费在线播放| 精品人妻熟女av久视频| 一本一本综合久久| 午夜激情欧美在线| 午夜激情欧美在线| 桃色一区二区三区在线观看| 女的被弄到高潮叫床怎么办 | 成人国产麻豆网| 深夜精品福利| videossex国产| 在线免费观看不下载黄p国产 | 精品久久久久久久久亚洲 | 色av中文字幕| 亚洲欧美精品综合久久99| 亚洲男人的天堂狠狠| 午夜日韩欧美国产| 成人国产一区最新在线观看| 成人特级av手机在线观看| 欧美色视频一区免费| 中文在线观看免费www的网站| 久久久久国产精品人妻aⅴ院| 午夜影院日韩av| a级毛片a级免费在线| 别揉我奶头~嗯~啊~动态视频| 九色国产91popny在线| 国产精品98久久久久久宅男小说| 久久久成人免费电影| 免费人成视频x8x8入口观看| 久久久久精品国产欧美久久久| 亚洲色图av天堂| www日本黄色视频网| 欧美性猛交╳xxx乱大交人| 少妇丰满av| 色综合站精品国产| 久9热在线精品视频| 级片在线观看| 小说图片视频综合网站| 美女黄网站色视频| 精品欧美国产一区二区三| 国产又黄又爽又无遮挡在线| 九九在线视频观看精品| 91麻豆精品激情在线观看国产| 两个人视频免费观看高清| 国产69精品久久久久777片| 国产伦一二天堂av在线观看| 亚洲成人中文字幕在线播放| 女人十人毛片免费观看3o分钟| 91麻豆av在线| 成人国产综合亚洲| 久久99热这里只有精品18| 国产高清有码在线观看视频| 亚洲久久久久久中文字幕| 欧美日韩瑟瑟在线播放| 婷婷色综合大香蕉| 日本a在线网址| 国产精品伦人一区二区| 白带黄色成豆腐渣| 亚洲美女视频黄频| 欧美xxxx黑人xx丫x性爽| 亚洲成a人片在线一区二区| 精品久久久久久久久久久久久| 又黄又爽又刺激的免费视频.| 小蜜桃在线观看免费完整版高清| 国产免费男女视频| 亚洲一区二区三区色噜噜| 日日啪夜夜撸| 日韩欧美国产一区二区入口| 18禁黄网站禁片免费观看直播| 99热只有精品国产| 国产伦人伦偷精品视频| 好男人在线观看高清免费视频| 国内精品久久久久久久电影| 18禁黄网站禁片午夜丰满| 国产69精品久久久久777片| 中文字幕免费在线视频6| 亚洲人成网站高清观看| 中国美白少妇内射xxxbb| 国产精品久久久久久久电影| 精华霜和精华液先用哪个| 国国产精品蜜臀av免费| 欧美成人一区二区免费高清观看| 乱系列少妇在线播放| 91麻豆av在线| 成人精品一区二区免费| 免费av不卡在线播放| 午夜影院日韩av| 国产成人福利小说| 免费观看的影片在线观看| 亚洲精品456在线播放app | 午夜福利在线观看免费完整高清在 | 国产一区二区亚洲精品在线观看| 久久久久国内视频| 可以在线观看的亚洲视频| 毛片女人毛片| 国产日本99.免费观看| a级一级毛片免费在线观看| 久久精品人妻少妇| 一夜夜www| 国产人妻一区二区三区在| 国产探花在线观看一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲无线在线观看| 成人二区视频| 亚洲av二区三区四区| 黄色女人牲交| 亚洲av免费高清在线观看| 窝窝影院91人妻| 日韩精品青青久久久久久| 欧美性感艳星| 别揉我奶头 嗯啊视频| 国产熟女欧美一区二区| 久久久久久国产a免费观看| 国产精品嫩草影院av在线观看 | 国产一区二区三区av在线 | 精品乱码久久久久久99久播| 99精品久久久久人妻精品| 欧美日本视频| 午夜影院日韩av| 国产一区二区激情短视频| 18禁裸乳无遮挡免费网站照片| 日本撒尿小便嘘嘘汇集6| 免费在线观看日本一区| 成年版毛片免费区| www日本黄色视频网| 18禁在线播放成人免费| 一个人观看的视频www高清免费观看| 极品教师在线视频| 国产精品1区2区在线观看.| 亚洲国产欧洲综合997久久,| 久久精品国产99精品国产亚洲性色| 99久久九九国产精品国产免费| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕精品亚洲无线码一区| 中国美白少妇内射xxxbb| 色噜噜av男人的天堂激情| 国产精华一区二区三区| 国产午夜精品论理片| 亚洲美女视频黄频| 午夜激情欧美在线| av天堂中文字幕网| 亚洲国产欧美人成| 色av中文字幕| 久久热精品热| 免费搜索国产男女视频| 十八禁国产超污无遮挡网站| 国产男靠女视频免费网站| 91在线精品国自产拍蜜月| 亚洲图色成人| 干丝袜人妻中文字幕| 在线观看舔阴道视频| 3wmmmm亚洲av在线观看| 国产欧美日韩一区二区精品| 国产亚洲精品综合一区在线观看| 国产免费av片在线观看野外av| 午夜精品久久久久久毛片777| 亚洲午夜理论影院| 自拍偷自拍亚洲精品老妇| 久久久久久国产a免费观看| 毛片女人毛片| 欧美日韩亚洲国产一区二区在线观看| 亚洲中文日韩欧美视频| 亚洲一区二区三区色噜噜| 97人妻精品一区二区三区麻豆| 日韩精品有码人妻一区| 国产一区二区激情短视频| 免费搜索国产男女视频| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 亚洲 国产 在线| 在线免费观看的www视频| 男女视频在线观看网站免费| 99久国产av精品| 一个人看视频在线观看www免费| 日韩精品青青久久久久久| 欧美性猛交╳xxx乱大交人| 免费在线观看成人毛片| 免费观看精品视频网站| 91久久精品电影网| 亚洲av免费在线观看| 久久久久九九精品影院| 色综合婷婷激情| 国产不卡一卡二| 亚洲精品影视一区二区三区av| 欧美国产日韩亚洲一区| 日韩精品中文字幕看吧| eeuss影院久久| 国国产精品蜜臀av免费| 亚洲av免费高清在线观看| 蜜桃亚洲精品一区二区三区| 麻豆久久精品国产亚洲av| 亚洲一级一片aⅴ在线观看| 欧美xxxx性猛交bbbb| 免费一级毛片在线播放高清视频| 久久久久久九九精品二区国产| 亚洲内射少妇av| 99久久精品国产国产毛片| 欧美一级a爱片免费观看看| 亚洲精品亚洲一区二区| 国产精品久久久久久av不卡| 尤物成人国产欧美一区二区三区| 欧美最新免费一区二区三区| 亚洲成人中文字幕在线播放| 免费观看人在逋| 99热6这里只有精品| 国产精品电影一区二区三区| 老熟妇仑乱视频hdxx| 欧美激情在线99| av专区在线播放| 91av网一区二区| 1024手机看黄色片| 亚洲国产精品合色在线| 乱人视频在线观看| 国产又黄又爽又无遮挡在线| 国产不卡一卡二| 18禁黄网站禁片午夜丰满| 亚洲精品一卡2卡三卡4卡5卡| 五月玫瑰六月丁香| 高清在线国产一区| 亚洲真实伦在线观看| 午夜免费激情av| 老司机深夜福利视频在线观看| 亚洲一级一片aⅴ在线观看| 成人二区视频| 久久精品人妻少妇| 久久久色成人| 亚洲av二区三区四区| 日韩强制内射视频| 国产 一区 欧美 日韩| 欧美日本亚洲视频在线播放| 国产精品美女特级片免费视频播放器| 99视频精品全部免费 在线| 日韩欧美免费精品| 亚洲人成网站高清观看| 久久久成人免费电影| 国产亚洲欧美98| 欧美日韩黄片免| 亚洲aⅴ乱码一区二区在线播放| 波多野结衣巨乳人妻| 99久久中文字幕三级久久日本| 男人狂女人下面高潮的视频| 欧美+亚洲+日韩+国产| 国产成人影院久久av| 国产伦一二天堂av在线观看| 露出奶头的视频| 午夜久久久久精精品| 久久国产乱子免费精品| 中文字幕精品亚洲无线码一区| 99视频精品全部免费 在线| 国产伦人伦偷精品视频| 亚洲精华国产精华液的使用体验 | 欧美日本亚洲视频在线播放| 午夜爱爱视频在线播放| av在线蜜桃| 精品久久久久久久末码| 欧美日韩瑟瑟在线播放| 午夜亚洲福利在线播放| 亚洲无线观看免费| 天堂√8在线中文| 99久久成人亚洲精品观看| 午夜影院日韩av| 国产人妻一区二区三区在| 女同久久另类99精品国产91| 日韩欧美免费精品| 中国美白少妇内射xxxbb| 亚洲熟妇熟女久久| 日韩欧美国产在线观看| 日韩精品有码人妻一区| 无人区码免费观看不卡| 国产精品国产高清国产av| 神马国产精品三级电影在线观看| 国产高清三级在线| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 国产淫片久久久久久久久| 熟妇人妻久久中文字幕3abv| 久久婷婷人人爽人人干人人爱| 99久久精品国产国产毛片| 国产精品野战在线观看| 在线天堂最新版资源| 哪里可以看免费的av片| 国产美女午夜福利| 国产av麻豆久久久久久久| 人人妻人人澡欧美一区二区| 深爱激情五月婷婷| 亚洲久久久久久中文字幕| 在线观看免费视频日本深夜| 免费看美女性在线毛片视频| 亚洲无线观看免费| 深爱激情五月婷婷| 亚洲,欧美,日韩| 亚洲国产欧美人成| 看片在线看免费视频| 日韩欧美三级三区| 国产精品一区二区免费欧美| 久久精品国产清高在天天线| 天堂影院成人在线观看| 日日撸夜夜添| 亚洲avbb在线观看| 国产主播在线观看一区二区| 老师上课跳d突然被开到最大视频| 搡老岳熟女国产| 久久99热6这里只有精品| 999久久久精品免费观看国产| 毛片一级片免费看久久久久 | 五月伊人婷婷丁香| 国产单亲对白刺激| 一级a爱片免费观看的视频| 热99re8久久精品国产| 我的老师免费观看完整版| 成人无遮挡网站| 免费无遮挡裸体视频| 午夜亚洲福利在线播放| 日日摸夜夜添夜夜添av毛片 | 日韩,欧美,国产一区二区三区 | 国产黄片美女视频| 亚洲黑人精品在线| 亚洲av五月六月丁香网| 亚洲国产精品成人综合色| 国产精品99久久久久久久久| 美女xxoo啪啪120秒动态图| 看黄色毛片网站| 变态另类丝袜制服| 99热这里只有是精品50| 久久久久久久久久黄片| 尤物成人国产欧美一区二区三区| 国产精品一及| 老司机午夜福利在线观看视频| 免费观看在线日韩| 禁无遮挡网站| 亚洲成人久久性| 成年人黄色毛片网站| 国产老妇女一区| 69av精品久久久久久| 日日摸夜夜添夜夜添小说| 亚洲综合色惰| 美女高潮的动态| 国产乱人伦免费视频| 国产国拍精品亚洲av在线观看| 久久久久性生活片| 在线播放无遮挡| 久久这里只有精品中国| 少妇人妻精品综合一区二区 | 听说在线观看完整版免费高清| 午夜久久久久精精品| 自拍偷自拍亚洲精品老妇| 日韩高清综合在线| 赤兔流量卡办理| 久久亚洲精品不卡| 国产成人影院久久av| 国产伦一二天堂av在线观看| 97超视频在线观看视频| 亚洲欧美清纯卡通| 真人做人爱边吃奶动态| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 久久久久免费精品人妻一区二区| 亚州av有码| 直男gayav资源| www.色视频.com| 直男gayav资源| 国产精品久久久久久精品电影| 极品教师在线视频| 国产伦精品一区二区三区视频9| 国产久久久一区二区三区| 婷婷色综合大香蕉| 在线播放国产精品三级| 国产精品久久久久久精品电影| 校园春色视频在线观看| 欧美日韩精品成人综合77777| 极品教师在线视频| 九九久久精品国产亚洲av麻豆| 性欧美人与动物交配| 最好的美女福利视频网| 成人鲁丝片一二三区免费| 91av网一区二区| 亚洲av电影不卡..在线观看| 婷婷六月久久综合丁香| 精品人妻一区二区三区麻豆 | 国产男人的电影天堂91| 女的被弄到高潮叫床怎么办 | 女人被狂操c到高潮| 国产精品久久久久久久电影| 亚洲黑人精品在线| 亚洲中文日韩欧美视频| 极品教师在线免费播放| 成人永久免费在线观看视频| 免费搜索国产男女视频| 一区二区三区激情视频| 国内精品宾馆在线| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩高清专用| 神马国产精品三级电影在线观看| 91精品国产九色| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美精品综合久久99| 亚洲美女视频黄频| 精品一区二区三区av网在线观看| 色综合站精品国产| 一个人看视频在线观看www免费| 欧美激情在线99| av天堂中文字幕网| 1024手机看黄色片| av中文乱码字幕在线| 91精品国产九色| 搡女人真爽免费视频火全软件 | netflix在线观看网站| 极品教师在线免费播放| 亚洲av五月六月丁香网| 毛片女人毛片| 国产精品久久久久久亚洲av鲁大| 极品教师在线视频| 色吧在线观看| 成人一区二区视频在线观看| av在线天堂中文字幕| 成人高潮视频无遮挡免费网站| 亚洲五月天丁香| 99视频精品全部免费 在线| 好男人在线观看高清免费视频| 亚洲 国产 在线| 美女高潮的动态| 男女下面进入的视频免费午夜| 一级av片app| 男女做爰动态图高潮gif福利片| 久久精品国产99精品国产亚洲性色| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 欧美精品国产亚洲| 尾随美女入室| 精品久久久久久久末码| 一进一出抽搐gif免费好疼| 真人一进一出gif抽搐免费| 日韩精品青青久久久久久| 韩国av在线不卡| 亚洲精品在线观看二区| 日韩欧美精品免费久久| 国产亚洲精品av在线| 看免费成人av毛片| 长腿黑丝高跟| 十八禁国产超污无遮挡网站| 欧美bdsm另类| 国产麻豆成人av免费视频| 麻豆一二三区av精品| eeuss影院久久| 国产中年淑女户外野战色| 欧美日韩瑟瑟在线播放| 亚洲av日韩精品久久久久久密| 级片在线观看| 男人舔奶头视频| 一本一本综合久久| 一个人看视频在线观看www免费| a级毛片a级免费在线| 女同久久另类99精品国产91| 国产精品精品国产色婷婷| 真实男女啪啪啪动态图| eeuss影院久久| 午夜老司机福利剧场| 老司机午夜福利在线观看视频| 精品福利观看| 美女xxoo啪啪120秒动态图| 亚洲午夜理论影院| 免费高清视频大片| 淫妇啪啪啪对白视频| 欧美一区二区精品小视频在线| 男插女下体视频免费在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av免费在线观看| 九九爱精品视频在线观看| 日本与韩国留学比较| 最后的刺客免费高清国语| 看免费成人av毛片| 亚洲成人中文字幕在线播放| 国产麻豆成人av免费视频| 免费无遮挡裸体视频| 噜噜噜噜噜久久久久久91| 亚洲av成人精品一区久久| 欧美bdsm另类| 最新在线观看一区二区三区| 大型黄色视频在线免费观看| 可以在线观看的亚洲视频| 国产精品综合久久久久久久免费| 蜜桃亚洲精品一区二区三区| 国产精品国产三级国产av玫瑰| 精品久久久久久久人妻蜜臀av| 亚洲七黄色美女视频| 国产精品1区2区在线观看.| 国产真实乱freesex| 欧美又色又爽又黄视频| 日本三级黄在线观看| 亚洲熟妇中文字幕五十中出| 国产免费男女视频| 毛片女人毛片| 日韩欧美精品v在线| 色av中文字幕| 精品午夜福利在线看| 色综合婷婷激情| 69av精品久久久久久| 欧美三级亚洲精品| 一进一出抽搐动态| 久久国产乱子免费精品| 国产成人a区在线观看| 听说在线观看完整版免费高清| av福利片在线观看| 一进一出抽搐动态| 久久亚洲精品不卡| 真人做人爱边吃奶动态| 啪啪无遮挡十八禁网站| 男女啪啪激烈高潮av片| 成人鲁丝片一二三区免费| 色吧在线观看| 十八禁国产超污无遮挡网站| 最近视频中文字幕2019在线8| 国产精品国产高清国产av| 别揉我奶头~嗯~啊~动态视频| 3wmmmm亚洲av在线观看| 免费人成在线观看视频色| 成人午夜高清在线视频| 国产爱豆传媒在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品粉嫩美女一区| h日本视频在线播放| 日韩 亚洲 欧美在线| 欧美bdsm另类| 国产激情偷乱视频一区二区| 一区二区三区四区激情视频 | 久久久精品大字幕| 岛国在线免费视频观看| 人妻夜夜爽99麻豆av| 男人狂女人下面高潮的视频| 十八禁网站免费在线| 内射极品少妇av片p| 午夜福利视频1000在线观看| 嫩草影视91久久| 精品日产1卡2卡| 免费高清视频大片| 亚洲av五月六月丁香网| 久久国内精品自在自线图片| 亚洲精品一区av在线观看| 久久国产精品人妻蜜桃| 亚洲av第一区精品v没综合| 麻豆久久精品国产亚洲av| 最新在线观看一区二区三区| 久久精品国产自在天天线| 日韩欧美国产在线观看| 国产成人a区在线观看| 成人国产麻豆网| 色尼玛亚洲综合影院| 热99在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 国产成人福利小说| 成人欧美大片| 黄色视频,在线免费观看| 一个人看视频在线观看www免费| 一边摸一边抽搐一进一小说| 又紧又爽又黄一区二区| 丝袜美腿在线中文| 午夜福利在线观看免费完整高清在 | 乱人视频在线观看| 亚洲av中文字字幕乱码综合| 天堂动漫精品| av在线老鸭窝| 免费电影在线观看免费观看| 人妻少妇偷人精品九色| 内射极品少妇av片p| 精品久久久久久久末码| 精品不卡国产一区二区三区| 欧美色视频一区免费| 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| 免费电影在线观看免费观看| a级一级毛片免费在线观看| 岛国在线免费视频观看| 国内精品美女久久久久久| 久久久久性生活片| 3wmmmm亚洲av在线观看| 天堂网av新在线| 免费看美女性在线毛片视频| 亚洲色图av天堂| 国产精品嫩草影院av在线观看 | 深夜精品福利| 日日夜夜操网爽| 国产黄色小视频在线观看| 国产精品亚洲一级av第二区| 久久午夜亚洲精品久久| 日韩一本色道免费dvd| 亚洲av电影不卡..在线观看| 99热只有精品国产| 午夜免费男女啪啪视频观看 | 久久精品夜夜夜夜夜久久蜜豆| 99热精品在线国产| 国产精品精品国产色婷婷| .国产精品久久| 非洲黑人性xxxx精品又粗又长| 少妇丰满av| 大又大粗又爽又黄少妇毛片口| 免费看美女性在线毛片视频| 一级黄色大片毛片| 国内精品宾馆在线| 午夜福利18| 亚洲国产精品久久男人天堂|