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

    利用2b-RAD測序結合HRM分析技術開發(fā)與梨矮生性狀相關的DNA分子標記

    2017-09-11 09:04:23肖玉雄王彩虹田義軻楊紹蘭李鼎立張海月
    中國農業(yè)科學 2017年15期
    關鍵詞:普通型分型基因組

    肖玉雄,王彩虹,田義軻,楊紹蘭,李鼎立,張海月

    ?

    利用2b-RAD測序結合HRM分析技術開發(fā)與梨矮生性狀相關的DNA分子標記

    肖玉雄,王彩虹,田義軻,楊紹蘭,李鼎立,張海月

    (青島農業(yè)大學園藝學院/青島市園藝植物遺傳改良與育種重點實驗室,山東青島 266109)

    【目的】果樹的矮生性狀是一種重要的農藝性狀,對果樹的集約化栽培具有重要意義。來自西洋梨實生變異品種‘Le Nain Vert’的矮生性狀受控于一個單顯性基因,目前關于該基因的序列信息等還不清楚。本研究的目的是開發(fā)與其緊密連鎖的DNA分子標記,為鑒定該基因提供依據(jù)。【方法】根據(jù)分離群體分組分析的原理,以‘矮生梨’ב茌梨’和‘2-3’ב綠寶石’2個F1雜交分離群體為試材,應用IIB型限制性內切酶的RAD技術(Restriction association site DNA,2b-RAD)對2對矮生型/普通型對比基因池進行基因組測序分析,在對比基因池間篩選出單核苷酸多態(tài)性(Single nucleotide polymorphism,SNP)位點,從中篩查出位于定位染色體上的SNPs,用高分辨率熔解曲線分析技術(High-resolution melting analysis,HRM)在群體上進一步檢測驗證,以確定其與位點的連鎖關系?!窘Y果】對4個樣本(即4個對比基因池)的2b-RAD標簽測序文庫的測序結果共產生67 186 260條reads,平均每個樣本測序reads數(shù)為16 796 565。將原始reads進行質量過濾后的統(tǒng)計結果表明,每個樣本獲得平均unique標簽數(shù)目為86 810,平均測序深度為77×,該測序深度能夠達到準確分型的標準。SOAP軟件定位結果表明,4個測序文庫中含有酶切位點的高質量reads占測序原始reads的70%以上,表明測序質量較好。在來自2個不同群體的矮生型基因池與普通型基因池間進行對比分析,初步篩選出SNP位點1 317個,其中有8個位于的定位染色體scaffold00074上。用HRM技術對這8個SNP標記在群體上的進一步檢測結果表明,在‘矮生梨’ב茌梨’群體上有2個SNP標記、在‘2-3’ב綠寶石’群體上有4個SNP標記表現(xiàn)出與位點共分離的特性,根據(jù)其擴增子的熔解曲線形狀差異,可有效區(qū)分矮生型和普通型表型。在來自‘矮生梨’ב茌梨’群體的215個雜種后代和來自‘2-3’ב綠寶石’群體的168個雜種后代中,未發(fā)現(xiàn)有標記與性狀的重組類型?!窘Y論】2b-RAD測序技術與HRM分析技術相結合,進行果樹重要農藝性狀分子標記的開發(fā),是一種行之有效的方法。基于這一策略,本研究鑒定獲得了4個與西洋梨矮生性狀單顯性基因連鎖的SNP標記。

    梨;矮生性狀;;SNP標記;2b-RAD;HRM

    0 引言

    【研究意義】果樹的矮化密植栽培因具有早果、高產、品質優(yōu)良、易于管理、能夠有效節(jié)約成本和增加利潤等優(yōu)點,已發(fā)展成為國內外極為重要的栽培模式。目前,實現(xiàn)矮密栽培的主要途徑是利用矮化砧木和矮生品種,因此,發(fā)掘和利用優(yōu)良的基因資源,進行矮化砧木和矮生品種的培育對果樹生產具有重要意義?!厩叭搜芯窟M展】西洋梨矮生型實生變異品種‘Le Nain Vert’具有樹體緊湊、矮小的特點,非常適合矮密栽培。研究發(fā)現(xiàn),矮生梨在莖、葉的解剖結構上有區(qū)別于普通型梨的顯著特征[1]。由于‘Le Nain Vert’的矮生性狀受控于顯性單基因[2],因此,該基因資源在梨矮生品種的培育方面具有重要應用價值。Wang等[3]將該基因命名為,并通過SSR標記將其定位于西洋梨基因組第16連鎖群上。最近,李煒等[4]依據(jù)蘋果與梨基因組間高度的保守性和共線性關系,以二者的參考基因組序列為基礎,又篩選鑒定了一組與連鎖的DNA分子標記,將該基因定位到了西洋梨的染色體片段scaffold00074上。利用簡化基因組測序技術進行大規(guī)模高通量SNP分型是當前國際上動植物基因組學研究熱點之一。美國Oregon大學發(fā)明的RAD(Restriction Association site DNA)技術是廣泛應用的簡化基因組測序技術[5]。2007年,Miller等[6]發(fā)表了應用芯片雜交的RAD技術。2008年,Baird等[7]將RAD和NGS(Next-generation sequencing)結合,建立了RAD-seq技術。2012年,Peterson等[8]對RAD技術進行改進,提出了ddRAD(double digest RAD)技術。針對RAD技術流程復雜的問題,Elshire等[9]曾提出GBS(genotyping- by-sequencing)技術。隨后,Wang等[10]推出基于llB型限制性內切酶的2b-RAD技術,該技術流程簡便,重復性和標簽代表性高,基于混合泊松分布模型的de novo SNP分型新算法(iML),可使假陽性率得到較大幅度降低[11]。目前,RAD測序技術在許多物種上得到了應用,例如大麥(L. )[12]、黑麥草()[13]等,但在果樹上還未見報道?!颈狙芯壳腥朦c】為了最終達到分離鑒定的目的,還需要大量的遺傳標記進一步對該基因進行更為精細的染色體定位。因此,在以往研究的基礎上,繼續(xù)高效開發(fā)更多與其緊密連鎖的遺傳標記是十分必要的。將2b-RAD測序技術與HRM分析技術相結合,進行果樹重要農藝性狀SNP標記的開發(fā),是一個值得探索的途徑?!緮M解決的關鍵問題】隨著西洋梨全基因組序列的公布[14],使大量獲得SNP標記成為可能。本研究利用分離群體分組分析的原理[15],采用2b-RAD技術在矮生型/普通型梨對比基因池間進行比較分析,初步篩查可能與目標性狀相關的SNP位點,然后利用高效的高分辨率熔解曲線(HRM)分析技術進行大規(guī)模的群體分析,以確認與梨矮生性狀決定基因共分離的SNP標記,為進一步實現(xiàn)該基因的精細定位及分離鑒定提供重要依據(jù)。

    1 材料與方法

    試驗于2016年4—11月在青島農業(yè)大學園藝學院植物遺傳改良與育種重點實驗室進行。

    1.1 基因組DNA 的提取與對比基因池的構建

    以‘矮生梨’(L.)ב茌梨’(Rehd.)及‘2-3’(L.×Rehd.)ב綠寶石’(Nakai.)2個F1雜交分離群體(即群體1和群體2)為試材?!妗汀?-3’為矮生型親本;‘茌梨’和‘綠寶石’為普通型親本。群體1共有215株,其中矮生型和普通型分別有115株和100株。群體2共有168株,其中矮生型和普通型各有89株和79株。2016年4月取梨樹的幼葉,用天根植物基因組DNA提取試劑盒提取并純化基因組DNA,用瓊脂糖凝膠電泳和NanoDrop 2000超微量分光光度計(Thermo Fisher Scienti?c, USA)對DNA的質量和濃度進行檢測,并將其濃度稀釋到10 ng·L-1。

    從每個群體中分別抽取15個矮生型雜種和15個普通型雜種,將其DNA等量混合,形成矮生型基因池D1和D2以及普通型基因池S1和S2。

    1.2 2b-RAD測序分析與SNP位點篩選

    利用2b-RAD技術構建梨4個基因池的標簽測序文庫,并在Hiseq2500 v2平臺進行single-end測序(此項工作由上海歐易生物醫(yī)學科技有限公司完成)。將原始reads進行質量過濾,剔除不含有XⅠ酶切識別位點的序列,剔除低質量序列(即大于10個堿基的質量分數(shù)小于20的序列),剔除有10個以上連續(xù)相同堿基的序列。從梨參考基因組序列(https:// www.rosaceae.org/species/pyrus/pyrus_communis/genome_v1.0)中提取包含XⅠ酶切位點的標簽作為參考序列。將每個樣本的高質量reads利用SOAP軟件[16](參數(shù)設置為:-M4–v2–r0)定位到參考序列上。獲得可用于分型的unique標簽數(shù)目及深度。利用最大似然法(maximum likelihood,ML)進行SNP位點的分型(分型條件:只留取標簽內最多有3個SNP的標簽位點;分型標簽深度設置為500以下。選取在樣本間能分型的SNP位點,采用種群結構分析軟件對樣本間及組間差異SNP位點篩選數(shù)據(jù)進行分析。

    1.3 與連鎖的SNP標記的鑒定

    對在4個對比基因池上產生的SNP標記進行分析,篩選出矮生型基因池與普通型基因池間的多態(tài)性SNP標記,然后查找位于西洋梨染色體片段scaffold00074上的SNP位點,在其兩側合適位置設計引物(Primer 4.0 software),用HRM技術分析這些SNP標記在群體上的分離行為,以確定是否為與緊密連鎖的遺傳標記。

    1.4 高分辨率熔解曲線(HRM)分析

    HRM分析在LightCycler?480Ⅱ熒光定量PCR儀(Roche)上進行。反應試劑來自LightCycler?480 High Resolution Melting Master 試劑盒。反應體系為10 μL,內含10 ng基因組DNA 1 μL,1×Master Mix 5 μL,2.0 mmol·L-1MgCl21 μL,上、下游引物(0.2 μmol·L-1)各1 μL,ddH2O 1 μL。PCR擴增程序為95℃預變性10 min,然后按95℃變性10 min、55℃退火15 s、72℃延伸10 s的程序進行45個循環(huán)。在PCR循環(huán)結束后,立即對擴增產物進行HRM檢測,程序為:95℃ 1 min、40℃ 1 min、65℃1 s;在65℃升溫至95℃的過程中,以25次/℃的頻率收集熒光信息,最后降溫至40℃。

    HRM分析用LightCycler?480的Gene Scanning軟件(1.5 version)進行。

    2 結果

    2.1 2b-RAD測序結果

    對4個樣本(即4個近等基因池)的2b-RAD標簽測序文庫的測序結果共產生67 186 260條reads,平均每個樣本測序reads數(shù)為16 796 565。將原始reads進行質量過濾后的統(tǒng)計結果表明,每個樣本獲得平均unique標簽數(shù)目為86 810,平均測序深度為77×。測序深度能夠達到準確分型的標準,可以進行下一步SNP標記分型分析。SOAP軟件定位結果表明,4個測序文庫中含有酶切位點的高質量reads占測序原始reads的70%以上(表1),表明梨文庫的測序質量較好。

    2.2 矮生型/普通型對比基因池間SNP位點篩選

    選取在樣本間能分型的SNP位點,分別在來自群體1和群體2的4個對比基因池間進行比較分析,共篩查出可在矮生型/普通型表型間進行分型的SNP標記1 317個,進一步的標簽篩查結果表明,其中有8個SNP位于西洋梨染色體片段scaffold00074(即的定位染色體片段,長度441 741 bp)上,其染色體位置及樣品分型結果見表2。

    表1 2b-RAD測序結果

    D1、D2:分別來自群體1和群體2的矮生型基因池;S1、S2:分別來自群體1和群體2的普通型基因池。下同

    D1 and D2 represent the dwarf bulk from population 1 and population 2, respectively. S1 and S2 represent the standard bulk from population 1 and population 2, respectively. The same as below

    表2 定位到西洋梨染色體片段scaffold00074上的8個SNP標記及其分型信息

    方框為堿基差異位點位置The box represents the position of the base difference

    2.3 與位點共分離的SNP標記鑒定

    針對以上8個定位于scaffold00074上的SNP位點,分別設計8對PCR引物(表3),并分別在群體1和群體2上進行HRM分析,最終在群體1上檢測到SNP1和SNP7兩個與位點共分離的SNP標記,在群體2上檢測到SNP1、SNP6、SNP7和SNP8 4個可與位點共分離的SNP標記。根據(jù)其擴增子的熔解曲線形狀差異,可有效地區(qū)分矮生型和普通型表型(圖1、圖2)。在群體1的215個雜種后代和群體2的168個雜種后代中,未發(fā)現(xiàn)有標記與性狀的重組類型。

    表3 8個SNP位點群體分型檢測所用的HRM 引物

    圖 1 SNP1和SNP7在‘矮生梨’ב茌梨’群體中的HRM分析結果

    圖 2 SNP1、SNP6、SNP7和SNP8在‘2-3’ב綠寶石’群體中的HRM 分析結果

    3 討論

    單核苷酸多態(tài)性(SNP)是基因組中最常見的遺傳多態(tài)性。SNP標記已在果樹的遺傳圖譜構建、基因標記與定位等研究中發(fā)揮了重要作用[17-20]。與RFLP、AFLP、SSR等傳統(tǒng)基因分型技術相比,SNP分子具有分布廣泛、代表性強、穩(wěn)定性好、易于分析等優(yōu)勢。SNP分子標記的獲得往往依賴于大規(guī)模、高通量的基因組數(shù)據(jù)分析[21]。對果樹來說,其復雜的基因組是造成測序分析成本高的主要原因。

    簡化基因組測序技術2b-RAD是改進的RAD技術,使測序更加簡易、精確,已在高密度遺傳圖譜的構建、群體的遺傳學研究、群體進化分析、標記開發(fā)等領域得到了應用[22-28]。高分辨率熔解曲線分析是一種高效的遺傳變異檢測方法,被廣泛的應用于突變掃描、基因分型和標記開發(fā)等方面。目前,HRM技術已發(fā)展成為SNP標記分型的重要手段[29-31]。李煒等[4]曾依據(jù)蘋果與梨基因組間高度的保守性和共線性關系,以蘋果基因組序列為參考設計了大量引物,通過HRM分析僅篩查到1個與梨矮生性狀決定基因緊密連鎖的SNP標記。為了提高SNP位點的檢測效率,本研究利用2b-RAD測序技術通過對4個對比基因池間共有標簽的序列分析,篩查到了大量的可能與目標性狀相關的SNP位點,對于初步篩查到的SNP標記,再利用高通量的HRM分析技術進行群體上的大規(guī)模驗證,以確定其是否為與目標性狀緊密連鎖的遺傳標記。這樣,可避免了在群體上繼續(xù)用2b-RAD技術進行驗證分析造成的高成本問題。本研究基于這一策略,鑒定獲得了一組與梨矮生性狀決定基因位點連鎖的新標記。可見,將2b-RAD測序技術與HRM分析技術相結合,進行果樹重要農藝性狀分子標記的開發(fā),是一種行之有效的方法。

    根據(jù)不同親本構建的雜交群體的遺傳背景是有差異的,因此,依據(jù)不同群體開發(fā)的遺傳標記是不可能完全相同的。本研究結果顯示,在群體1上,鑒定出了2個與連鎖的SNP標記;在群體2上,鑒定出了4個與連鎖的SNP標記。另外,對在不同群體上鑒定到的相同的SNP標記來說,其HRM分析熔解曲線形狀也可能會出現(xiàn)差異。出現(xiàn)這一現(xiàn)象的原因,也與群體的遺傳基礎差別有關。本研究從2個不同的雜種群體上共檢測到了4個新的與梨緊密連鎖的DNA分子標記。這些標記,不僅可以作為目標基因定位的參考依據(jù),也可成為目標性狀標記輔助選擇的有力工具。

    4 結論

    應用2b-RAD技術對2對矮生型/普通型對比基因池進行基因組測序分析,在對比基因池間篩選出單核苷酸多態(tài)性(SNP)位點1 317個,其中有8個SNP位點位于的定位染色體scaffold0007上。用HRM分析技術對這8個SNP位點進行群體上的進一步檢測,共鑒定出4個與位點共分離的SNP標記,其中SNP1和SNP7為2對群體共有的SNP標記。這一研究結果對的精細定位及其候選基因的篩查具有重要參考。

    References

    [1] Chen B Y, Wang C H, Tian Y K, Chu Q, Hu C. Anatomical characteristics of young stems and mature leaves of dwarf pear., 2015, 186: 172-179.

    [2] Rivalta L, Dradi M, Rosati C. Thirty years of pear breeding activity at ISF Forli, Italy., 2002, 596(1): 233-238.

    [3] Wang C H, Tian Y K, Buck E J, Gardiner S E, Dai H Y, Jia Y. Genetic mapping of PcDw determining pear dwarf trait., 2011, 136(1): 48-53.

    [4] 李煒, 田義軻, 王彩虹, 白牡丹, 侯董亮. 通過HRM技術篩查與梨矮生性狀決定位點緊密連鎖的SNP標記. 園藝學報, 2015, 42(2): 214-220.

    Li W, Tian Y K, Wang C H, Bai M D, Hou D L. Screening of SNP markers tightly linked tolocus determining pear dwarf trait using HRM technology., 2015, 42(2): 214-220. (in Chinese)

    [5] Davey J W, Blaxter M L. RADSeq: Next-generation population genetics., 2011, 9(5): 416-423.

    [6] Miller M R, Dunham J P, Amores A, Cresko W A, Johnson E A. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers., 2007, 17(2): 240-248.

    [7] Baird N A, Etter P D, Atwood T S, Currey M C, Shiver A L, Lewis Z A, Selker E U, Cresko W A, Johnson E A. Rapid SNP discovery and genetic mapping using sequenced RAD markers., 2008, 3(10): e3376.

    [8] Peterson B K, Weber J N, Kay E H, Fisher H S, Hoekstra H E. Double digest RADseq: An inexpensive method for de novo SNP discovery and genotyping in model and non-model species., 2012, 7(5): e37135.

    [9] Elshire R J, Glaubitz J C, Sun Q, Poland J A, Kawamoto K, Buckler E S, Mitchell S E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species., 2011, 6(5): e19379.

    [10] Wang S, Meyer E, Mckay J K, Matz M V. 2b-RAD: A simple and flexible method for genome-wide genotyping., 2012, 9(8): 808-810.

    [11] Dou J, Zhao X, Fu X, Jiao W, Wang N, Zhang L, Hu X, Wang S, Bao Z. Reference-free SNP calling: improved accuracy by preventing incorrect calls from repetitive genomic regions., 2012, 7(1): 1-9.

    [12] Chutimanitsakun Y, Nipper R W, Cuesta-Marcos A, Cistué L, Corey A, Filichkina T, Johnson E A, Hayes P M. Construction and application for QTL analysis of a restriction site associated DNA (RAD) linkage map in barley., 2011, 12(1): 4.

    [13] Pfender W F, Saha M C, Johnson E A, Slabaugh M B. Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in., 2011, 122(8): 1467-1480.

    [14] Chagné D, Crowhurst R N, Pindo M, Thrimawithana A, Deng C, Ireland H, Fiers M, Dzierzon H, Cestaro A, Fontana P, Bianco L, Lu A, Storey R, Kn?bell M, Saeedl M, Montanari1 S, Kim Y K, Nicolini D, Larger S, Stefani E, Allan A C, Bowen J, Harvey I, Johnston J, Malnoy M, Troggio M, Perchepied L, Sawyer G, Wiedowl C, Won K, Viola R, Hellens R P, Brewer L, Bus V G , Schaffer R J, Gardiner S E, Velasco R. The draft genome sequence of european pear (L. ‘Bartlett’)., 2014, 9: e92644.

    [15] Michelmore R W, Paran I, Kesseli R V. Identification of markers linked to disease-resistance genes by bulked segregant analysis: A rapid method to detect markers in specific genomic regions by using segregating populations., 1991, 88: 9828-9832.

    [16] Li R, Yu C, Li Y, Lam T W, Yiu S M, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment., 2009, 25 (15): 1966-1967.

    [17] Baumgartner I O, Kellerhals M, Costa F, Dondini L, Pagliarani G, Gregori R, Tartarini S, Leumann L, Laurens F, Patocchi A. Development of SNP-based assays for disease resistance and fruit quality traits in apple (Borkh.) and validation in breeding pilot studies., 2016, 12(3): 1-21.

    [18] Zeballos J L, Abidi W, Giménez R, Monforte A J, Moreno M á, Gogorcena Y. Mapping QTLs associated with fruit quality traits in peach [(L.) Batsch] using SNP maps., 2016, 12(3): 1-17.

    [19] Wu J, Li L T, Li M, Khan M A, Li X G, Chen H, Yin H, Zhang S L. High-density genetic linkage map construction and identification of fruit-related QTLs in pear using SNP and SSR markers., 2014, 65(20): 5771-5781.

    [20] SUN R, CHANG Y, YANG F, WANG Y, LI H, ZHAO Y, CHEN D, WU T, ZHANG X, HAN Z. A dense SNP genetic map constructed using restriction site-associated DNA sequencing enables detection of QTLs controlling apple fruit quality., 2015, 16(1): 747.

    [21] 李全林, 王義發(fā), 韓晴, 袁政, 沈雪芳. 糯玉米‘滬五彩花糯1號’品系特異性2b-RAD分子標記的開發(fā)及應用. 植物生理學報, 2016, 52(5): 669-677.

    Li Q L, Wang Y F, Han Q, Yuan Z, Shen X F. Development and application of molecular markers for event-speci?c identi?cation of waxy corn ‘Huwucaihuanuo 1’ based on 2b-RAD technique., 2016, 52(5): 669-677. (in Chinese)

    [22] Jiao W, Fu X, Dou J, Li H, Su H, Mao J, Yu Q, Zhang L, Hu X, Huang X, Wang Y, Wang S, Bao Z. High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc., 2013, 21(1): 85-101.

    [23] Seetharam A S, Stuart G W. Whole genome phylogeny for 21species using predicted 2b-RAD fragments., 2013, 1: e226.

    [24] Fletcher R S, Mullen J L, Yoder S, Bauerle W L, Reuning G, Sen S, Meyer E, Juenger T E, McKay J K. Development of a next-generation NIL library infor dissecting complex traits., 2013, 14(1): 655-655.

    [25] Guo Y, Yuan H, Fang D, Song L, Liu Y, Liu Y, Wu L, Yu J, Li Z, Xu X. An improved 2b-RAD approach (2b-RAD) offering genotyping tested by a rice (L.) F2 population., 2014, 15 (1): 956.

    [26] Li C, Li Y, Bradbury P J, Wu X, Shi Y, Song Y, Zhang D, Rodgers-Melnick E, Buckler E S, Zhang Z. Construction of high-quality recombination maps with low-coverage genomic sequencing for joint linkage analysis in maize., 2015, 13: 78.

    [27] Pecoraro C, Babbucci M, Villamor A, Franch R, Papetti C, Leroy B, Ortega-Garcia S, Muir J, Rooker J, Arocha F, Murua H, Zudaire I, Chassot E, Bodin N, Tinti F, Bargelloni L, Cariani A. Methodological assessment of 2b-RAD genotyping technique for population structure inferences in yellowfin tuna ()., 2015, 25: 43-48.

    [28] Pauletto M, Carraro L, Babbucci M, Lucchini R, Bargelloni L, Cardazzo B. Extending RAD tag analysis to microbial ecology: A comparison between multilocus sequence typing and 2b-RAD to investigate listeriamonocytogenes genetic structure., 2016, 16(3): 823-835.

    [29] Lu Z, Niu L, Chagné D, Cui G, Pan L, Foster T, Zhang R P, Zeng W F, Wang Z Q. Fine mapping of the temperature-sensitive semi-dwarf (Tssd) locus regulating the internode length in peach ()., 2016, 36(2): 1-11.

    [30] Wang C H, Li W, Tian Y K, Hou D L, Bai M D. Development of molecular markers for genetic and physical mapping of thelocus in pear (L.)., 2016, 91(3): 299-307.

    [31] 趙俊生, 楊曉燕, 曾祥有, 鐘聲, 方靜, 羅劍斌, 曾運友, 向旭. 利用SNP 分子標記分析化橘紅種質資源. 分子植物育種, 2016, 14(5): 1203-1211.

    Zhao J S, Yang X Y, Zeng X Y, Zhong S, Fang J, Luo J B, Zeng Y Y, Xiang X. Analysis on germplasm resources ofusing SNP molecular markers., 2016, 14(5): 1203-1211.

    (責任編輯 趙伶俐)

    Development of DNA Molecular Markers for the Dwarf Trait in Pear through the Method of 2b-RAD Sequencing and HRM Analysis

    XIAO YuXiong, WANG CaiHong, TIAN YiKe, Yang ShaoLan, Li DingLi, ZHANG HaiYue

    (College of Horticulture, Qingdao Agricultural University/Qingdao Key Laboratory of Genetic Improvement and Breeding in Horticultural Plants, Qingdao 266109, Shandong)

    【Objective】As an important agronomic trait, the dwarf character of tree is significant to fruit intensive culture. The dwarf traits of pear originated from variety ‘Le Nain Vert’, a chance seedling of, is determined by a dominant gene, which sequence information is still unknown. Development of DNA molecular markers tightly linked to thelocus could provide important information for identifying of this gene.【Method】Using 2 different F1populations obtained from ‘Aishengli’בChili’ and ‘2-3’בLvbaoshi’, respectively, as plant materials, 2 pairs of bulks for dwarf type/standard type were sequenced through IIB restriction association site DNA (2b-RAD) technology, and then single nucleotide polymorphism (SNP) markers were detected between the dwarf and standard bulks based on the principle of bulked segregant analysis. After that, those SNPs located on themapped chromosome were selected and tested on the whole population to confirm their linkage relationship to thegene by high-resolution melting (HRM) analysis.【Result】A total of 67 186 260 reads were produced from the 2b-RAD sequencing of the 4 samples (2 pairs of bulks). That means the average reads per sample was 16 796 565. Statistical analysis of the filtered raw reads showed that the average unique tags per sample were 86 810, and the average sequencing depth was 77× which was enough for accurate genotyping. For the 4 libraries, mapping through SOAP software indicated that the high quality reads containing the restriction site accounted for more than 70% of the raw reads, which means the high quality of sequencing. Totally 1 317 polymorphic SNP markers were screened between the dwarf bulks and standard bulks derived from the 2 different populations, and 8 of them were mapped on the-localized chromosome scaffold00074. Further detection of the 8 SNPs through HRM analysis showed that 2 and 4 SNP markers co-segregated with thelocus were identified from population ‘Aishengli’בChili’ and population ‘2-3’בLvbaoshi’, respectively. According to the different shapes of the melting curves of amplicons, the dwarf/normal phenotype could be distinguished effectively. The recombinants for each marker and the target trait were not found in both of the two populations, which contained 215 and 168 progenies, respectively. 【Conclusion】The combination of 2b-RAD sequencing and HRM analysis technology is an efficient method for exploiting molecular markers of important agronomic traits in fruit trees. In this study, a total of 4 SNP markers co-segregated with thelocus were identified based on this strategy.

    pear; dwarf traits;; SNP marker; 2b-RAD; HRM

    2017-01-20;接受日期:2017-04-13

    國家自然科學基金(31372049)、山東省現(xiàn)代農業(yè)產業(yè)技術體系果品創(chuàng)新團隊(SDAIT-06-06)

    肖玉雄,E-mail:1427449415@qq.com。通信作者王彩虹,E-mail:chwang@qau.edu.cn

    猜你喜歡
    普通型分型基因組
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    調峰保供型和普通型LNG接收站罐容計算
    煤氣與熱力(2021年5期)2021-07-22 09:02:14
    周大勇從濕論治普通型新型冠狀病毒肺炎經(jīng)驗
    失眠可調養(yǎng),食補需分型
    普通型新型冠狀病毒肺炎的CT影像特點及其與臨床表現(xiàn)的關系
    便秘有多種 治療須分型
    基于分型線驅動的分型面設計研究
    康復新液治療普通型手足口病112例臨床觀察
    頸椎病中醫(yī)辨證分型與影像探討
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    亚洲av成人精品一区久久| 亚洲成人精品中文字幕电影| 国产成人一区二区在线| 美女 人体艺术 gogo| 少妇的逼好多水| 日韩人妻高清精品专区| 又黄又爽又免费观看的视频| 一级av片app| 嫩草影院新地址| 国产乱人视频| 国产精品久久久久久久久免| 国产精品久久久久久亚洲av鲁大| 国产精品美女特级片免费视频播放器| 成人国产一区最新在线观看| 成年人黄色毛片网站| 日韩中文字幕欧美一区二区| 久久久久九九精品影院| 精品一区二区三区视频在线观看免费| 欧美黑人巨大hd| 特级一级黄色大片| 男女视频在线观看网站免费| 18+在线观看网站| 成人性生交大片免费视频hd| 99热只有精品国产| 一级黄色大片毛片| 大又大粗又爽又黄少妇毛片口| 国产69精品久久久久777片| 国产一区二区亚洲精品在线观看| 亚洲av.av天堂| 长腿黑丝高跟| 免费观看人在逋| 亚洲精华国产精华精| 久久精品影院6| 2021天堂中文幕一二区在线观| 国产精品乱码一区二三区的特点| 99热精品在线国产| 一级av片app| 日本欧美国产在线视频| 国产一区二区在线av高清观看| 免费看a级黄色片| 全区人妻精品视频| 欧美高清成人免费视频www| 中出人妻视频一区二区| 毛片女人毛片| 国产精品爽爽va在线观看网站| 亚洲成人中文字幕在线播放| 熟女电影av网| 国产精品野战在线观看| 精品久久久久久成人av| 国产美女午夜福利| 中亚洲国语对白在线视频| 极品教师在线免费播放| 国内揄拍国产精品人妻在线| 亚洲一区二区三区色噜噜| 久久久久免费精品人妻一区二区| 日本三级黄在线观看| 欧美色视频一区免费| 俄罗斯特黄特色一大片| 国产男人的电影天堂91| 18禁黄网站禁片免费观看直播| 两个人的视频大全免费| 亚洲成人久久性| 亚洲欧美日韩东京热| 男女那种视频在线观看| 亚洲一区高清亚洲精品| 色播亚洲综合网| 综合色av麻豆| 国产成人福利小说| 精品午夜福利视频在线观看一区| 久久久久九九精品影院| 亚洲av一区综合| 国产一区二区亚洲精品在线观看| 国产精品久久久久久精品电影| 国产乱人视频| 永久网站在线| 国语自产精品视频在线第100页| 免费人成在线观看视频色| 国产精品久久久久久久久免| 1024手机看黄色片| 最后的刺客免费高清国语| 我要搜黄色片| 有码 亚洲区| 国产欧美日韩精品亚洲av| 国产精品av视频在线免费观看| 又粗又爽又猛毛片免费看| 麻豆久久精品国产亚洲av| 国产午夜福利久久久久久| 一区二区三区激情视频| 全区人妻精品视频| 婷婷精品国产亚洲av| 国产免费男女视频| 国产精品亚洲美女久久久| 搡老熟女国产l中国老女人| 国产激情偷乱视频一区二区| 精品久久久久久久久久免费视频| 亚洲第一电影网av| 欧美3d第一页| 人妻少妇偷人精品九色| 校园人妻丝袜中文字幕| 久久久久久久午夜电影| 日本精品一区二区三区蜜桃| 深夜a级毛片| 欧美+亚洲+日韩+国产| 成人特级黄色片久久久久久久| 中国美白少妇内射xxxbb| 十八禁国产超污无遮挡网站| 成年女人毛片免费观看观看9| 淫秽高清视频在线观看| 无遮挡黄片免费观看| 最好的美女福利视频网| 免费看光身美女| 韩国av一区二区三区四区| 国产大屁股一区二区在线视频| 亚洲欧美日韩高清专用| 毛片一级片免费看久久久久 | 国产精品98久久久久久宅男小说| 色5月婷婷丁香| 少妇的逼水好多| 国产又黄又爽又无遮挡在线| 内地一区二区视频在线| 美女 人体艺术 gogo| 可以在线观看的亚洲视频| 精品欧美国产一区二区三| 一个人观看的视频www高清免费观看| 在线观看美女被高潮喷水网站| 国产一区二区亚洲精品在线观看| 长腿黑丝高跟| 国产精品一区二区三区四区久久| 久久婷婷人人爽人人干人人爱| www日本黄色视频网| 亚洲真实伦在线观看| 久久久久久九九精品二区国产| 一个人看的www免费观看视频| netflix在线观看网站| 精品国内亚洲2022精品成人| 久久久久久久久久久丰满 | 又爽又黄a免费视频| 国产大屁股一区二区在线视频| 又爽又黄a免费视频| 成熟少妇高潮喷水视频| 老司机福利观看| 夜夜夜夜夜久久久久| 日本撒尿小便嘘嘘汇集6| 麻豆国产av国片精品| 亚洲精品久久国产高清桃花| 久久久久久久精品吃奶| 啪啪无遮挡十八禁网站| 欧美成人性av电影在线观看| 亚洲欧美日韩卡通动漫| 狠狠狠狠99中文字幕| 国产一区二区三区在线臀色熟女| 亚洲18禁久久av| 长腿黑丝高跟| 欧美一区二区国产精品久久精品| 黄色女人牲交| 精品乱码久久久久久99久播| 久久久久久伊人网av| 久久这里只有精品中国| 搡女人真爽免费视频火全软件 | 日本免费一区二区三区高清不卡| 免费av不卡在线播放| 久久久久国产精品人妻aⅴ院| 少妇高潮的动态图| 亚洲精品一卡2卡三卡4卡5卡| av.在线天堂| 亚洲一区二区三区色噜噜| 国产精品久久久久久久电影| 久久精品91蜜桃| 12—13女人毛片做爰片一| 精品一区二区三区视频在线观看免费| 欧美日韩黄片免| 国产爱豆传媒在线观看| av在线观看视频网站免费| 亚洲自拍偷在线| 美女大奶头视频| 婷婷色综合大香蕉| 日韩精品有码人妻一区| 一a级毛片在线观看| 国产精品无大码| 草草在线视频免费看| 日本 av在线| 亚洲美女视频黄频| 国产蜜桃级精品一区二区三区| avwww免费| 色综合亚洲欧美另类图片| 久久久久久久久大av| 国产av不卡久久| 成年女人永久免费观看视频| a级毛片免费高清观看在线播放| 少妇裸体淫交视频免费看高清| 免费搜索国产男女视频| 色精品久久人妻99蜜桃| 桃红色精品国产亚洲av| 国产亚洲av嫩草精品影院| 欧美xxxx性猛交bbbb| 成人二区视频| 久久精品国产亚洲av天美| 久久久精品大字幕| av在线老鸭窝| 午夜免费成人在线视频| 日本-黄色视频高清免费观看| 三级国产精品欧美在线观看| 性色avwww在线观看| 欧美精品国产亚洲| 嫩草影院精品99| 亚洲七黄色美女视频| 久久午夜亚洲精品久久| 97超级碰碰碰精品色视频在线观看| 两人在一起打扑克的视频| 九九爱精品视频在线观看| 亚洲在线自拍视频| 免费av观看视频| 少妇熟女aⅴ在线视频| 欧美区成人在线视频| 国内少妇人妻偷人精品xxx网站| 日韩欧美国产一区二区入口| 亚洲精品日韩av片在线观看| 精品乱码久久久久久99久播| 有码 亚洲区| 91在线精品国自产拍蜜月| 久久国内精品自在自线图片| 国产精品一区www在线观看 | 国产午夜福利久久久久久| 神马国产精品三级电影在线观看| 国内揄拍国产精品人妻在线| 长腿黑丝高跟| 69人妻影院| 国内精品宾馆在线| 一级av片app| 亚洲欧美日韩高清在线视频| 国产探花极品一区二区| 国产亚洲精品久久久久久毛片| av在线蜜桃| 亚洲成人中文字幕在线播放| 亚洲18禁久久av| 91在线精品国自产拍蜜月| 一本一本综合久久| 97超级碰碰碰精品色视频在线观看| 亚洲精品乱码久久久v下载方式| 精品国内亚洲2022精品成人| 99久久九九国产精品国产免费| 国产免费一级a男人的天堂| 日本与韩国留学比较| 国产私拍福利视频在线观看| 嫩草影院精品99| 最后的刺客免费高清国语| 精品无人区乱码1区二区| 久久精品国产亚洲av天美| 国产激情偷乱视频一区二区| 国产探花极品一区二区| 国产免费一级a男人的天堂| av国产免费在线观看| 高清日韩中文字幕在线| 欧美性感艳星| 男人的好看免费观看在线视频| 日本黄色片子视频| 亚洲图色成人| 真人做人爱边吃奶动态| 婷婷亚洲欧美| 成人国产一区最新在线观看| 婷婷六月久久综合丁香| 少妇丰满av| 久久久久久伊人网av| 亚洲av日韩精品久久久久久密| 两个人视频免费观看高清| 国产精品一区www在线观看 | 亚洲第一区二区三区不卡| 91久久精品电影网| 亚洲国产精品sss在线观看| 高清日韩中文字幕在线| 国产一区二区三区在线臀色熟女| 狠狠狠狠99中文字幕| 麻豆国产av国片精品| 亚洲图色成人| 亚洲中文日韩欧美视频| 婷婷亚洲欧美| 97碰自拍视频| 99久久精品热视频| 欧美性猛交╳xxx乱大交人| a在线观看视频网站| 欧美丝袜亚洲另类 | 免费看a级黄色片| 深夜精品福利| 男女视频在线观看网站免费| 真人一进一出gif抽搐免费| 如何舔出高潮| 中文字幕久久专区| 国产毛片a区久久久久| 久久亚洲精品不卡| 国产精品嫩草影院av在线观看 | 日韩一本色道免费dvd| 欧美性感艳星| 乱系列少妇在线播放| 高清在线国产一区| 看黄色毛片网站| 成年女人看的毛片在线观看| 3wmmmm亚洲av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美精品免费久久| 成人亚洲精品av一区二区| 亚洲精华国产精华液的使用体验 | 亚洲av一区综合| 久久久色成人| 亚洲av不卡在线观看| 在线国产一区二区在线| 在线观看66精品国产| 成人三级黄色视频| 色综合站精品国产| 99热6这里只有精品| 一本精品99久久精品77| 久久精品国产亚洲av涩爱 | 黄色视频,在线免费观看| 女人被狂操c到高潮| 搞女人的毛片| 身体一侧抽搐| 国产主播在线观看一区二区| 国产一区二区三区在线臀色熟女| 俺也久久电影网| 热99re8久久精品国产| 久久国内精品自在自线图片| 亚洲欧美日韩高清专用| 老司机午夜福利在线观看视频| 少妇人妻精品综合一区二区 | 免费观看人在逋| 中亚洲国语对白在线视频| 精品人妻熟女av久视频| 国产精品伦人一区二区| 人妻丰满熟妇av一区二区三区| 国内精品一区二区在线观看| 国产一级毛片七仙女欲春2| 欧美日韩精品成人综合77777| 香蕉av资源在线| 制服丝袜大香蕉在线| 狂野欧美白嫩少妇大欣赏| 日韩欧美国产在线观看| 欧美精品啪啪一区二区三区| 亚洲专区中文字幕在线| 欧美bdsm另类| 男女之事视频高清在线观看| 成熟少妇高潮喷水视频| 国产av麻豆久久久久久久| 国产一级毛片七仙女欲春2| 日本精品一区二区三区蜜桃| 97热精品久久久久久| 日韩欧美在线二视频| 少妇高潮的动态图| 国产不卡一卡二| 日韩精品有码人妻一区| 久久精品国产亚洲av香蕉五月| 亚洲七黄色美女视频| 国产视频内射| 亚洲性久久影院| av天堂在线播放| 99九九线精品视频在线观看视频| 麻豆精品久久久久久蜜桃| 国产精品不卡视频一区二区| 免费看a级黄色片| 极品教师在线免费播放| 国产视频内射| 特级一级黄色大片| 亚洲精品日韩av片在线观看| 最近中文字幕高清免费大全6 | 中文字幕人妻熟人妻熟丝袜美| 黄色配什么色好看| 国产主播在线观看一区二区| 女人十人毛片免费观看3o分钟| 99热这里只有精品一区| 变态另类成人亚洲欧美熟女| 在线观看一区二区三区| 亚洲精品乱码久久久v下载方式| 国产精品99久久久久久久久| 亚洲 国产 在线| 久久人人爽人人爽人人片va| 最后的刺客免费高清国语| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品乱码久久久v下载方式| 黄色日韩在线| 亚洲av二区三区四区| av国产免费在线观看| 女人十人毛片免费观看3o分钟| 亚洲欧美精品综合久久99| 少妇的逼水好多| 日本黄大片高清| 长腿黑丝高跟| 简卡轻食公司| 大又大粗又爽又黄少妇毛片口| 国产极品精品免费视频能看的| 国产私拍福利视频在线观看| 一区二区三区免费毛片| 亚洲男人的天堂狠狠| 18+在线观看网站| 1000部很黄的大片| 国产一区二区激情短视频| 亚洲国产色片| 国产欧美日韩精品一区二区| 久久久久久九九精品二区国产| 国产亚洲91精品色在线| 真人一进一出gif抽搐免费| 国产精品电影一区二区三区| 中文字幕av成人在线电影| 噜噜噜噜噜久久久久久91| 97超级碰碰碰精品色视频在线观看| 亚洲精华国产精华液的使用体验 | 久久久久国内视频| 国产色爽女视频免费观看| 白带黄色成豆腐渣| 一级a爱片免费观看的视频| 18禁裸乳无遮挡免费网站照片| 久久久久精品国产欧美久久久| 亚洲中文字幕日韩| 美女高潮喷水抽搐中文字幕| 麻豆成人午夜福利视频| 99久久精品一区二区三区| 欧美+亚洲+日韩+国产| 搡老妇女老女人老熟妇| 九九爱精品视频在线观看| 黄色日韩在线| 亚洲av成人av| 国产精品一区www在线观看 | 日韩欧美三级三区| 熟女人妻精品中文字幕| 亚洲精华国产精华液的使用体验 | 午夜福利成人在线免费观看| 日本成人三级电影网站| 国产精品精品国产色婷婷| 人妻丰满熟妇av一区二区三区| 欧美日韩精品成人综合77777| 中国美女看黄片| 精品久久久噜噜| 免费在线观看影片大全网站| 亚洲黑人精品在线| 亚洲第一电影网av| 亚洲人成网站在线播放欧美日韩| 精华霜和精华液先用哪个| 啪啪无遮挡十八禁网站| 亚洲不卡免费看| 午夜福利高清视频| 狂野欧美激情性xxxx在线观看| 精品日产1卡2卡| 久久精品91蜜桃| 露出奶头的视频| 成人二区视频| 国产精品亚洲一级av第二区| 长腿黑丝高跟| 91精品国产九色| 国产精品久久久久久久电影| 99在线人妻在线中文字幕| 简卡轻食公司| 国产精品国产高清国产av| 久久久久久久久大av| 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 国产精品无大码| 国产精品无大码| 亚洲va在线va天堂va国产| 麻豆一二三区av精品| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放| 国产又黄又爽又无遮挡在线| 久久久久久九九精品二区国产| 免费大片18禁| 老司机福利观看| 99久久九九国产精品国产免费| 最新中文字幕久久久久| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看| 久久精品国产自在天天线| 毛片女人毛片| 午夜免费男女啪啪视频观看 | 亚洲经典国产精华液单| 天堂√8在线中文| 最新在线观看一区二区三区| 国语自产精品视频在线第100页| 美女高潮喷水抽搐中文字幕| 久久精品国产亚洲网站| 成人高潮视频无遮挡免费网站| 99热网站在线观看| 99久久精品国产国产毛片| 国产精品永久免费网站| 成人精品一区二区免费| 精品久久久久久久久久久久久| 亚洲最大成人av| 男人舔奶头视频| 91麻豆精品激情在线观看国产| 亚洲av熟女| 亚洲在线自拍视频| 深夜a级毛片| 俺也久久电影网| 美女被艹到高潮喷水动态| 一级a爱片免费观看的视频| 色精品久久人妻99蜜桃| 亚洲三级黄色毛片| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| 小蜜桃在线观看免费完整版高清| 亚洲自拍偷在线| 亚洲欧美日韩高清专用| 日韩国内少妇激情av| 国产精品,欧美在线| 国产久久久一区二区三区| 久久99热这里只有精品18| 最新在线观看一区二区三区| 看十八女毛片水多多多| www.色视频.com| 国产69精品久久久久777片| 国产午夜福利久久久久久| 搡老妇女老女人老熟妇| 免费一级毛片在线播放高清视频| 美女大奶头视频| 熟女人妻精品中文字幕| h日本视频在线播放| 18禁在线播放成人免费| 久久久久久久亚洲中文字幕| 日本色播在线视频| 亚洲人成伊人成综合网2020| 日韩中文字幕欧美一区二区| 91在线精品国自产拍蜜月| 窝窝影院91人妻| 国产精品,欧美在线| 亚洲av电影不卡..在线观看| 欧美又色又爽又黄视频| 欧美在线一区亚洲| 国产在线男女| 久久中文看片网| aaaaa片日本免费| 三级国产精品欧美在线观看| 啪啪无遮挡十八禁网站| 欧美日韩国产亚洲二区| 精品久久久久久久久av| 久久人妻av系列| 草草在线视频免费看| 成年女人毛片免费观看观看9| 97碰自拍视频| 最近在线观看免费完整版| 女同久久另类99精品国产91| 成人午夜高清在线视频| 国产亚洲精品综合一区在线观看| 亚洲一级一片aⅴ在线观看| 老司机福利观看| 亚洲 国产 在线| 国产精品嫩草影院av在线观看 | 亚洲人成网站在线播放欧美日韩| 亚洲aⅴ乱码一区二区在线播放| 久久久久免费精品人妻一区二区| 国产成人福利小说| 国产大屁股一区二区在线视频| 午夜福利18| 免费在线观看成人毛片| 91久久精品电影网| 久久久久九九精品影院| 熟女电影av网| 国产精品久久视频播放| 日韩,欧美,国产一区二区三区 | a级毛片a级免费在线| 国产探花在线观看一区二区| 国产精品三级大全| 色5月婷婷丁香| 日本a在线网址| 尾随美女入室| 亚洲精品粉嫩美女一区| 黄色欧美视频在线观看| 久久亚洲精品不卡| 麻豆av噜噜一区二区三区| 两个人的视频大全免费| 日本一二三区视频观看| 精品不卡国产一区二区三区| 日韩av在线大香蕉| 搞女人的毛片| 极品教师在线免费播放| 亚洲国产日韩欧美精品在线观看| 欧美一区二区国产精品久久精品| 99热精品在线国产| 久久人人爽人人爽人人片va| 国产成人a区在线观看| 国产成年人精品一区二区| 国产在线精品亚洲第一网站| 他把我摸到了高潮在线观看| 免费黄网站久久成人精品| 男女下面进入的视频免费午夜| 又爽又黄无遮挡网站| 欧美日本亚洲视频在线播放| 久99久视频精品免费| 老熟妇仑乱视频hdxx| 国产美女午夜福利| 91久久精品电影网| 亚洲成人久久爱视频| 男女啪啪激烈高潮av片| 日韩欧美精品免费久久| 韩国av一区二区三区四区| 国产精品国产三级国产av玫瑰| 一本精品99久久精品77| 久久久久久久久久黄片| 又黄又爽又刺激的免费视频.| 亚洲精品成人久久久久久| 国产精品嫩草影院av在线观看 | 精品一区二区三区人妻视频| 最近视频中文字幕2019在线8| 最近最新中文字幕大全电影3| 亚洲综合色惰| 级片在线观看| 亚洲av.av天堂| 久久人妻av系列| 亚洲熟妇中文字幕五十中出| 久久人妻av系列| 欧美精品国产亚洲| 日日夜夜操网爽| 久久人人爽人人爽人人片va| 午夜精品一区二区三区免费看| 亚洲一级一片aⅴ在线观看| 黄色丝袜av网址大全| 久久久久久久久大av| av在线老鸭窝| 国产欧美日韩精品亚洲av| 亚洲三级黄色毛片| 色哟哟·www|