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

    小麥類病斑突變基因 lm3的定位

    2016-12-30 05:56:28楊佳秀耿皆飛李倩倩劉錄祥謝彥周王成社
    麥類作物學(xué)報(bào) 2016年12期
    關(guān)鍵詞:親本突變體病斑

    楊佳秀,耿皆飛,李倩倩,劉錄祥,謝彥周,王成社

    (1.西北農(nóng)林科技大學(xué)農(nóng)學(xué)院,陜西楊陵 712100; 2.中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所,北京 100081)

    小麥類病斑突變基因 lm3的定位

    楊佳秀1,耿皆飛1,李倩倩1,劉錄祥2,謝彥周1,王成社1

    (1.西北農(nóng)林科技大學(xué)農(nóng)學(xué)院,陜西楊陵 712100; 2.中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所,北京 100081)

    為進(jìn)一步對(duì)類病斑突變基因lm3進(jìn)行克隆和功能研究,以EMS誘變普通小麥花培品系H261所得后代中選育出的一個(gè)穩(wěn)定遺傳的小麥類病斑突變體LF2010為研究材料,通過(guò)BSA(bulked segregant analysis)、MBSA(multiple bulked segregant analysis)、小麥90K基因芯片結(jié)合BSA三種不同的基因定位策略對(duì)該突變基因進(jìn)行定位。結(jié)果表明,lm3基因位于小麥6B染色體的長(zhǎng)臂,與其兩側(cè)標(biāo)記SWES2和Xbarc134的遺傳距離分別為13.1 cM和3.8 cM,為一個(gè)新的小麥類病斑突變基因。

    小麥;類病斑突變;基因定位;lm3

    在小麥中,目前已報(bào)道的僅有7個(gè)類病斑突變體,分別為M66[11]、C591(M8)[12]、AIM9[13]、HLP[8]、Ning 7840[14]、Yanzhan 1/Zaosui 30分離群體中的一個(gè)穩(wěn)定斑點(diǎn)后代[15]、LF2010[16]。不同類病斑突變體的表型存在不同程度的差異,并且不同類病斑表型對(duì)各自農(nóng)藝性狀、抗病性等的影響也各有不同,如突變體HLP與野生型相比產(chǎn)量性狀差異不明顯,對(duì)葉銹菌成株期抗病性顯著增強(qiáng)[8];突變體M66雖然產(chǎn)量較親本明顯下降,但對(duì)白粉病[11]、條銹病和葉銹病[9]的抗病性增強(qiáng)。迄今為止,小麥中僅有3個(gè)類病斑突變基因被定位,只有不斷創(chuàng)造、發(fā)掘、定位更多的小麥類病斑基因,才能推動(dòng)小麥類病斑突變機(jī)理、基因克隆等研究地進(jìn)行。

    通常情況下,基因定位是一個(gè)費(fèi)時(shí)費(fèi)力的過(guò)程,尤其是在小麥這種基因組復(fù)雜且龐大的植物中。為了盡快確定目標(biāo)染色體,不同的基因定位策略隨著生物技術(shù)的發(fā)展相繼被嘗試。BSA(bulked segregant analysis)法是快速定位目標(biāo)基因的有效策略[17]。但試驗(yàn)中我們發(fā)現(xiàn),當(dāng)目標(biāo)基因附近標(biāo)記密度較低時(shí),運(yùn)用BSA法則可能會(huì)漏掉距離目標(biāo)基因相對(duì)較遠(yuǎn)的分子標(biāo)記,從而導(dǎo)致篩選到連鎖標(biāo)記的可能性大大降低,有時(shí)甚至?xí)Y不到連鎖的分子標(biāo)記,導(dǎo)致定位失敗。2012年,Hiebert等[18]提出了新的基因定位策略MBSA(multiple bulked segregant analysis)來(lái)提高BSA的靈敏度,Ghazvini等[19]對(duì)其做了具體地改進(jìn)及驗(yàn)證。此外,隨著基因芯片技術(shù)的出現(xiàn),使得高通量SNP基因分型得以實(shí)現(xiàn),進(jìn)而為在小麥這樣的復(fù)雜基因組中構(gòu)建高分辨率的遺傳連鎖圖譜提供了快速而有效的策略。2014年,Wang等[20]利用小麥90K芯片對(duì)8個(gè)不同作圖群體進(jìn)行了遺傳變異分析,將46 977對(duì)SNP分子標(biāo)記進(jìn)行遺傳定位,并利用其中的6個(gè)群體構(gòu)建了包含40 267個(gè)遺傳位點(diǎn)的整合圖譜。這些標(biāo)記數(shù)量可觀的高密度遺傳圖譜為小麥基因組的研究提供強(qiáng)大的支持。在小麥中,9K、90K、660K芯片相繼被開(kāi)發(fā),結(jié)合BSA,可以利用大量的SNP對(duì)一個(gè)較小的群體進(jìn)行富集分析,為小麥基因定位研究開(kāi)辟了一條全新、快速、高效的新策略。

    LF2010是由EMS誘變普通小麥花培品系H261得到的一個(gè)穩(wěn)定遺傳的類病斑突變體,其斑點(diǎn)性狀受1對(duì)隱性核基因控制,且表達(dá)受光照和溫度的影響,推測(cè)LF2010為一類新的突變類型[16]。由于已定位的3個(gè)小麥類病斑突變基因分別為lm[14]、lm1[15]和lm2[15],因此,這個(gè)類病斑突變基因暫時(shí)命名為lm3。本研究以突變體LF2010為研究材料,利用BSA、MBSA、小麥90K基因芯片結(jié)合BSA三種不同的基因定位策略,對(duì)突變基因lm3進(jìn)行定位研究,以期為lm3基因的克隆和功能研究以及最終揭示小麥類病斑突變發(fā)生機(jī)制奠定基礎(chǔ)。

    1 材料與方法

    1.1 材料及其種植

    LF2010是普通小麥花培品系H261經(jīng)過(guò)EMS(ethane methyl sulfonate) 誘變所得后代中選育出的小麥類病斑突變體,經(jīng)過(guò)連續(xù)5代自交已遺傳穩(wěn)定。2013年5月,將LF2010與普通六倍體小麥品種中國(guó)春雜交,隨后分別播種親本、F1、F2及F2∶3,從三葉期開(kāi)始,以葉片是否出現(xiàn)黃色斑點(diǎn)為標(biāo)準(zhǔn),鑒定植株的表型。

    1.2 DNA的提取

    一是構(gòu)建事業(yè)平臺(tái)。智庫(kù)要有事情做,有研究的問(wèn)題、項(xiàng)目、課題等,要有相關(guān)的科學(xué)有效的理論、方法、信息、手段等的支持,也要有專業(yè)、有水平團(tuán)隊(duì)支撐;

    將LF2010與中國(guó)春雜交后自交得到的190個(gè)F2∶3代株系作為定位群體。按CTAB[21]法提取LF2010、中國(guó)春及定位群體各單株的全基因組DNA。

    1.3 突變基因 lm3的初定位

    1.3.1 BSA策略

    參照Michelmore等[22]文章中所述BSA方法,隨機(jī)選取F2∶3定位群體中8株綠色野生型單株的DNA等量混合構(gòu)建正常性狀池(BW),8株斑點(diǎn)型單株的DNA等量混合構(gòu)建突變性狀池(BM),然后將親本LF2010和中國(guó)春間有差異的SSR分子標(biāo)記在BW和BM之間進(jìn)行篩選,在不同性狀池間存在多態(tài)性的標(biāo)記即為連鎖標(biāo)記。

    1.3.2 MBSA策略

    參照Ghazvini等[19]文章中所述MBSA方法,隨機(jī)選取F2∶3定位群體中2株斑點(diǎn)型單株的DNA(含隱性純合的目標(biāo)基因)等量混合組建一個(gè)迷你池(mini bulk),一共組建14個(gè)這樣的mini bulk。將親本間存在多態(tài)性的所有SSR標(biāo)記在這14個(gè)mini bulk中進(jìn)行篩選,在5個(gè)以上mini-bulk中與隱性親本帶型一致的標(biāo)記即為連鎖標(biāo)記。

    1.3.3 小麥90K芯片結(jié)合BSA策略

    參照Michelmore等[22]文章中所述BSA方法,選取F2∶3定位群體中6株綠色野生型單株的DNA等量混合構(gòu)建正常性狀池(C-BW),6株斑點(diǎn)型單株DNA等量混合構(gòu)建突變性狀池(C-BM),送至北京康普森生物公司進(jìn)行小麥90K芯片檢測(cè),2次樣本重復(fù)。根據(jù)Wang等[20]提供的SNP位點(diǎn)信息確定基因池間的差異SNP在小麥各個(gè)染色體上的富集情況,最終根據(jù)這個(gè)SNP富集結(jié)果即可初步確定目標(biāo)基因所在的染色體。

    將90K芯片初步確定染色體上的多態(tài)性SNP對(duì)照Wang等[20]文章中小麥與二穗短柄草、水稻、高粱的對(duì)應(yīng)關(guān)系表,可以得到SNP對(duì)應(yīng)的二穗短柄草、水稻、高粱各自的基因區(qū)間,從而得到兩基因池之間SNP的富集區(qū)。利用富集區(qū)SNP的序列設(shè)計(jì)基于小麥基因組序列的SSR引物并在基因池間檢測(cè)其是否與目標(biāo)基因連鎖。

    1.4 突變基因 lm3的精細(xì)定位

    從F2∶3定位群體中隨機(jī)挑選一定單株用于驗(yàn)證得到的連鎖標(biāo)記及檢測(cè)初定位所得染色體上其他在親本間有多態(tài)性的分子標(biāo)記的多態(tài)性。將突變體LF2010特有的帶型記作A,中國(guó)春特有的帶型記作B,雜合帶型記作H;若連鎖標(biāo)記為顯性標(biāo)記,則單株帶型記作A/C或者B/D。利用JoinMap 4.0軟件分析,并繪制相應(yīng)的遺傳連鎖圖。

    2 結(jié)果與分析

    2.1 遺傳分析

    類病斑突變體LF2010與中國(guó)春雜交,F(xiàn)1代單株全部為野生型的綠色表型,說(shuō)明該基因?yàn)殡[性基因;經(jīng)卡方測(cè)驗(yàn),F(xiàn)2群體中綠色表型單株與斑點(diǎn)表型單株的比例符合3∶1的分離比例,F(xiàn)2∶3群體符合1∶2∶1的分離比例,表明在這個(gè)群體中的目標(biāo)性狀受1對(duì)核基因控制。以上說(shuō)明類病斑突變體LF2010的斑點(diǎn)性狀在LF2010與中國(guó)春雜交的F2群體中為單基因隱性遺傳(表1)。這與本實(shí)驗(yàn)室以前研究的結(jié)果相同。

    2.2 突變基因lm3的BSA策略定位

    杜麗芬[23]在京花一號(hào)與LF2010雜交的F2∶3定位群體中篩選到2對(duì)分子標(biāo)記Xwmc388和cfa2187,本研究基于此,將這2對(duì)分子標(biāo)記在親本LF2010、中國(guó)春及其F2∶3群體單株構(gòu)建的兩基因池間進(jìn)行篩選,結(jié)果顯示,cfa2187在親本間無(wú)差異,Xwmc388在親本及基因池之間均存在多態(tài)性,初步表明Xwmc388與類病斑突變基因lm3連鎖,但由于Xwmc388有多個(gè)位點(diǎn),分別位于3A、4D、5A、6A、6B、7A染色體上,所以并不能根據(jù)分子標(biāo)記Xwmc388初步判斷目標(biāo)基因lm3所在的染色體。

    2.3 突變基因lm3的MBSA策略定位

    選取Ghazvini等[19]列出的均勻分布在小麥21條染色體上帶型單一的423對(duì)引物中標(biāo)記Xwmc388所在的6條染色體上的110對(duì)SSR引物,及本實(shí)驗(yàn)室已有的位于這6條染色體上的27對(duì)SSR引物(表2),進(jìn)行親本間多態(tài)性的篩選。在這137對(duì)引物中,68對(duì)在雙親間表現(xiàn)出多態(tài)性,作為進(jìn)一步篩選的候選引物。將候選引物在14個(gè)mini-bulk間進(jìn)行PCR擴(kuò)增,經(jīng)檢測(cè)發(fā)現(xiàn),分子標(biāo)記Xbarc134 在14個(gè)池中均與隱性親本LF2010帶型一致(圖1A),Xwmc417在6個(gè)mini-bulk中與隱性親本LF2010帶型一致(圖1B),初步表明Xwmc417、Xbarc134與類病斑突變基因lm3連鎖。Xwmc417有2個(gè)位點(diǎn),分別位于6A、6B染色體上,Xbarc134僅有1個(gè)位點(diǎn),位于小麥6B染色體上,因此可以初步確定突變基因lm3位于小麥6B染色體上。

    2.4 突變基因lm3的小麥90K芯片結(jié)合BSA策略定位

    根據(jù)整合的高密度90K SNP遺傳連鎖圖譜[20],可以得到基因池間差異SNP所在的染色體位置。從圖2A中可以看出,6B染色體上多態(tài)性SNP數(shù)目最多,為95個(gè),占已知位置的總SNP數(shù)目的36.7%,與圖2B結(jié)果基本一致,說(shuō)明控制斑點(diǎn)的基因最可能在6B染色體上,但仍需進(jìn)一步地驗(yàn)證分析。

    表1 突變體LF2010與中國(guó)春雜交后代性狀分離比

    Table 1 Segregation for spotted mutation

    群體Population總株(家系)數(shù)Totalnumberofplants(lines)綠色株(家系)數(shù)Numberofgreenplants(lines)分離家系數(shù)Numberofsegregationlines斑點(diǎn)株(家系)數(shù)Numberofspottedplants(lines)分離比Segregationratioχ2P值PvalueF11111F2190139513∶10.0250.84F2∶31804393441∶2∶10.2110.90

    M:DL2000;P1:LF2010;P2:中國(guó)春;1~14:分別由2個(gè)突變型單株所構(gòu)建的14個(gè)迷你池。

    M:DL2000; P1:LF2010; P2:Chinese Spring; 1~14:Fourteen mini bulks consisting of two mutant plants,respectively.

    圖1 SSR分子標(biāo)記Xbarc134(A)和Xwmc417(B)的電泳結(jié)果

    Fig.1 Electrophoresis of fragments amplified with SSR markers Xbarc134 (A) and Xwmc417 (B)

    A和B為2個(gè)重復(fù)。 A and B indicate two replicates.

    綜合2次重復(fù)結(jié)果,得到多態(tài)性SNP在6B染色體上的富集區(qū),對(duì)照Wang等[20]文中各物種之間的線性關(guān)系,得到6B染色體上SNP富集區(qū)對(duì)應(yīng)的二穗短柄草、水稻、高粱各自的基因區(qū)間(表3)。利用所得區(qū)間內(nèi)的二穗短柄草的基因?qū)?yīng)的SNP序列信息設(shè)計(jì)了69對(duì)基于小麥基因組的SSR引物(表2),其中37對(duì)在親本間有多態(tài)性。由于染色體已經(jīng)初步確定,試驗(yàn)重心已由確定目標(biāo)基因所在染色體轉(zhuǎn)為得到與目標(biāo)基因緊密連鎖的分子標(biāo)記,因此利用BSA法對(duì)37對(duì)親本間有差異的引物進(jìn)行篩選,結(jié)果得到6BLm15、6BLm22、6BLm59在池之間存在差異,初步表明這3對(duì)引物與lm3連鎖,同時(shí)驗(yàn)證了芯片混池快速確定染色體的結(jié)果,基本確定了突變基因lm3位于6B染色體上。

    2.5 連鎖標(biāo)記的驗(yàn)證及連鎖圖的構(gòu)建

    隨機(jī)選取F2∶3定位群體中110個(gè)單株進(jìn)一步驗(yàn)證引物Xwmc388、Xwmc417、Xbarc134、6BLm15、6BLm22、6BLm59是否與類病斑突變基因lm3連鎖。數(shù)據(jù)分析結(jié)果表明,該突變基因位于SSR標(biāo)記Xbarc134和6BLm22之間。由于Xbarc134和6BLm22均位于6BL染色體上,因此搜索6BL上的SSR及EST-SSR引物,將得到的6BL染色體上親本間有差異且?guī)颓逦囊锞贔2∶3群體中進(jìn)行分析,最終得到如圖3的遺傳連鎖圖譜。該圖譜包括了11對(duì)SSR引物、1對(duì)EST-SSR引物、3對(duì)基于芯片數(shù)據(jù)開(kāi)發(fā)的SSR引物,其中距離目標(biāo)基因最近的兩側(cè)引物為SWES2和Xbarc134,其遺傳距離分別為13.1 cM 和3.8 cM。

    表2 本實(shí)驗(yàn)室已有的SSR引物(No.1-No.27)及本研究基于SNP序列重新合成的SSR引物(No.28-No.96)

    Table 2 SSR primers kept in our laboratory(No.1-No.27) and resynthesized based on SNP sequences in this study(No.28-No.96)

    編號(hào)No.SSR引物名稱NameofSSRprimer正向引物序列(5′?3′)Forwardprimersequence(5′?3′)反向引物序列(5′?3′)Reverseprimersequence(5′?3′)1barc310GGGCGGCGCATGTGCACCTAGCGTGGAAGCGACTAAATCAACT2barc356CGCTAGAGCTGTTTGAGGGGAGGAGCGCATGTAGGGGGAGGCTTCTTTT3gwm480TGCTGCTACTTGTACAGAGGACCCGAATTGTCCGCCATAG4cfa2234AATCTGACCGAACAAAATCACATCGGAGAGTATTAGAACAGTGCC5gwm155CAATCATTTCCCCCTCCCAATCATTGGAAATCCATATGCC6wmc388TGTGCGGAATGATTCAATCTGTGGCCATTAGACTGCAATGGTTT7wmc331CCTGTTGCATACTTGACCTTTTTGGAGTTCAATCTTTCATCACCA8barc334ATCCGCGTGTCAAACTTCTTCCGGGCTGGCTGGGCTAAATG9cfd54TCGTTCCAAAATGCATGAAAAAGGGCCAGAAATCTGTGTG10gwm205CGACCCGGTTCACTTCAGAGTCGCCGTTGTATAGTGCC11barc117TCATGCGTGCTAAGTGCTAAGAGGGCAGGAAAAAGTGACT12barc141GGCCCATGGATAATTTTTGAAATGCAATTCGGCCAAAGAAGAAGTCA13barc330GCACTAAGCGCTCTTTATTTACCCTGCATCTGGTATGGAGA14barc319GCAGAGCTACGGCAATGTGCGTAAGTCCCGGAAGTAACAGAA15barc171GCGGGGTCATCTTAGTAACTCAAATAACTGTCAACGTTGGTTCACATTCA16barc113GCGCACAACAACGGACACTTAACAATTGGGACTCATTTAGCTTCTACTCGCCATTA17gwm169ACCACTGCAGAGAACACATACGGTGCTCTGCTCTAAGTGTGGG18gwm617GATCTTGGCGCTGAGAGAGACTCCGATGGATTACTCGCAC19barc79GCGTTGGAAAGGAGGTAATGTTAGATAGTCGTGGGTTACAAGTTTGGGAGGTCA20barc24CGCCTCTTATGGACCAGCCTATGCGGTGAGCCATCGGGTTACAAAG21barc178GCGTATTAGCAAAACAGAAGTGAGGCGACTAGTACGAACACCACAAAA22barc134CCGTGCTGCAAATGAACACAGTTGCCGGTTCCCATTGTCA23wmc168AACACAAAAGATCCAACGACACCAGTATAGAAGGATTTTGAGAG24gwm60TGTCCTACACGGACCACGTGCATTGACAGATGCACACG25cfa2110TCACTACCCGCATGAACAAATTCTGCACAAACCGTTCTGA26cfa2123CGGTCTTTGTTTGCTCTAAACCACCGGCCATCTATGATGAAG27gwm332AGCCAGCAAGTCACCAAAACAGTGCTGGAAAGAGTAGTGAAGC286BLm01TTTCTAACCCATTCGTGACTAGAATAGGGGTTATCTTCCAGA296BLm02GGAAGAAGGGAGGAAATAAATGCGTCATTTCATTGAATTTT306BLm03TTGGACAACAAGTAGAGCAGTTCGCTAGTGACATACAGTGAA316BLm04GGAAAATTCAATGAAATGACGGAGAATTATACTCCCCAAACG326BLm05TCCCAAATGATGCATATAGACCTTTCAGGAGCAATTAGTCAA336BLm06TGATGGAGCAGTCAATATTTCGTAGATCCTGATCCCAAAGAC346BLm07CCAAATGTTCATCTTGTTTGTTGGTTGGTGATCAGAACTTAC356BLm08TTGCTTATTTGCTTACCTCACACTTCGACTGTACAACACCAT366BLm09GGTTTCATAGGCCTGTTCTCCATCACATGAAACAGATCACA376BLm10CAGATCATGTCTCCAAGGTCATCAGAAACAAGTGCTAATCG386BLm11AGGTTTGGACATGTGGATAATTTTACCTTTTATATGCCATCG396BLm12CTAAAAATTGACTGCACCATTTTTGTTTGGTAGTCTGTTCGT

    (續(xù)表3 Continued table 3)

    編號(hào)No.SSR引物名稱NameofSSRprimer正向引物序列(5′?3′)Forwardprimersequence(5′?3′)反向引物序列(5′?3′)Reverseprimersequence(5′?3′)406BLm13ATACTCGAAACCTTCTGAACCGTGAAATTTCGATCCCAAAG416BLm14CAAAGTCTGCTAGCTTTAGGAAGCAGCTATAGGCGTTCTAGT426BLm15GCACAGCCTTTAATCATCTACACACGAAGCTCTCTATCAATG436BLm16ACGTATGCCATCAGTAACCTACACCAAGCGTAAAGTGAATAA446BLm17GAGATCTTATGGGTTTCACCTCAGTTTCCAGAAGGAACATAA456BLm18CGCAAGTTGTATGTTGACTTCCCATCAGATAACACCAATAGC466BLm19TACTCCCTAATCTGCAGTAGCCATGTATGTAGCAGCAGAACA476BLm20TCATTGTACTGCGATGTTGTAGGTGGATGAAGCTAGGTAGAG486BLm21ATCGACATCAACTCGATCCTTTGTGGTTGCTAAAAGTGAC496BLm22ACCTTCATCTTGGACTTCTTGGTAGCGGTTGAGAGAAAGAG506BLm23AGAGCTTTGTCTTCACCTTCTTGAATTTCGATCTGAACAAAC516BLm24AAAAAGGTTCACATTGTTTCCAACTAATCAATGCATCCTCAA526BLm25TCTCTTCTTTTTAGGGCTTTCCAGTAGAAATCGAACTGAGGA536BLm26AAATGCCATGGTTAAACTACAGTAGGCTGAATGATGCTACC546BLm27ACCACCAAAAGATTTAGAAGCGCACCGCTAAGCAACACTA556BLm28AATGTTTTAACTCAGGCACAACCCAAGACATCAAAGTGTATC566BLm29GAATCGTGTGTTGATTTGGTCAGAAACAAGGAACGGTTAAT576BLm30CTGCTTATTGGAAATTGTTTCTTTAACCGATCTTCTGTTTTCA586BLm31ACACTCTTACATTGCTTTCCATATTCTGGTGGCAGATTTGTA596BLm32AGGTGTGACATGGAAACAAGTTATTAGTGGCCTCTCTCTCC606BLm33ACAGAGGGAGTAGTGAGTGGTCATTGATGTGTGTGTGTTCAT616BLm34ATGGAGTTGTACGAGCAATAACTTTGCAACTATGCTAATCGT626BLm35TCTGAATCACCCAGAAATAAATGATTAAGGCGTCTCCTCT636BLm36TTGTGAAAGGACAGAGTTCATTCCAGGTAAGATAACCAATCA646BLm37TCGAGATGAGTTTGAGCAGTGAGACGTACAGGTGATTAGA656BLm38GTGTCGAGGTCTTGGATGTAAGATAGAGCAGGTTTTCACC666BLm39GCTCAACATCATTTGAAAAACTGTTAAAAGTGGTGTCCAAAG676BLm40GACGTAGGCTGTTGGATCTCTAATTAATCGCAGCTGTACG686BLm41CTGGAACAGGAGGCAAAGCAGTGTAAATAACGTCGCTTC696BLm42TCTAGATCAAACTCCGAAACCGTTTTGCTTGAGCAGACG706BLm43AGAAGAAGAGAAGGACGTGACGTCGTTGTCTCCAGTGTAAAT716BLm44GTCAGAGGGCTAGTGGTAAATTTGTGCTGATAGTGTCAAAAA726BLm45TCAGCTAAAAGTGATGGCTACCATTTCGATTCCCAATCTT736BLm46GTTCAATGTATTGTTCCTTGGAACGATGACCCTTTCACTATT746BLm47ATAATTTCTATGCGGAAGAGCTCCCTAAATCCTGAAAGAAAC756BLm48TTCTTTCAGGATTTAGGGAAGGGATATTTTCCTTTTTCATCG766BLm49CAGACGATGATGACTTTATGCTCTGTTCTGATCTCTTTTGCT776BLm50CAAGCTGAAAACAGATCAAAGAGCTTGCACTACCCTAGAATC786BLm51CATGATACGCTACAATGTTGACCTCCCATCTCCATATTATTT796BLm52CGGTAAAATATTCTTCAGCAATGATAGATATCAATGGGGATG806BLm53CTGTTTTTAATGCCTCTTTTGCTGGTCCATATTAAAGCAGTG816BLm54TATATAGTGCGATGGTTTTGGGATAAAAACATCCATGTCACG826BLm55GCGATGTCAATGAAACACTACGGCTACTTGGTTGGAAACTAC836BLm56AGGGTAAACCAATTAAGAAGCACTTCACTGAGGAGGAGGAG846BLm57TGCATCTTAGCTAATTGTGGTGTCAATCGAATCCATAGAGGT856BLm58CCACCATTGCTTGTTAGTTTAGCTGCTACTGCTGCTATTATC866BLm59TTGATAGAGGCGTAGAAATGAGGCTTCTTATTTGAGATGGTT876BLm60AAGATAGATAACGTTCGTTGCGAGAGTGATGGAGTCAAGGTC886BLm61TAATTTTCATCGGTCATTCTGAGTAGTGGAGGAGTGATGGAT896BLm62ATCCATCACTCCTCCACTACTGGGAATGGGAGTTAGTACATT906BLm63ATTGCCCTTTCATACTTTTTCTGGTAGGGTACGTACATTTTG916BLm64TGAGAAGCAGAACCCAGTATACAATCATCCTCCCCTAGATA926BLm65TTTTACATGGAGAAACCACACAGTGAACAAGAAATCCATGTG936BLm66CACATGGATTTCTTGTTCACTGTTGTAGTGGTTGTGGTTTGT946BLm67GATGAAGACGATGAAGATGAGGTGCTACAACAGACATTAGGG956BLm68ATAAAAACAGTAGCGGTGGTTAATGTCTCACAGAAAATCGAA966BLm69TATAATACGCAGGCAACAGACATCCTTCCTCTGCTTGTTAAT

    表3 小麥、水稻、二穗短柄草和高粱的共線性關(guān)系

    Table 3 Macro-colinearity relationship among wheat,rice,Brachypodiumdistachyonand sorghum

    90KSNP標(biāo)記90KSNPmarker二穗短柄草基因Brachypodiumdis?tachyongene水稻基因Ricegene高粱基因Sorghumgenewsnp_CAP11_c3599_1741800Bradi3g50010.1LOC_Os02g42890.1Sb04g033650.1_PACid_1968195Ex_c20946_574Bradi3g51617.1wsnp_Ku_c5744_10166561Bradi3g51910.1LOC_Os02g45920.1Sb04g031640.1_PACid_1967964RAC875_c4420_371Bradi3g51940.1LOC_Os02g45950.1Sb04g031610.1_PACid_1967960Ra_c23874_338Bradi3g52367.2LOC_Os02g46990.1Sb04g030980.1_PACid_1967873wsnp_Ra_c7889_13482164Bradi3g55120.1Tdurum_contig7986_704Bradi3g55950.1Jagger_c1231_85Bradi3g55960.1LOC_Os02g57830.2Sb04g037890.1_PACid_1968698RAC875_c45987_132Bradi3g56110.1LOC_Os02g57990.1Sb04g037940.1_PACid_1968703CAP8_c7260_77Bradi3g56460.1LOC_Os02g49320.1Excalibur_c28014_99Bradi3g56460.1LOC_Os02g49320.1Sb04g029450.1_PACid_1967686Tdurum_contig12648_732Bradi3g60430.1LOC_Os02g58230.1Sb04g038220.1_PACid_1968733RAC875_c86296_99Bradi3g60430.1LOC_Os02g58230.1Sb04g038220.1_PACid_1968733BobWhite_c447_897Bradi3g60446.1RAC875_c28848_330Bradi3g60477.1LOC_Os02g58280.1Sb04g038270.1_PACid_1968738

    圖3 突變基因 lm3在小麥6BL染色體上的遺傳連鎖圖

    3 討 論

    在小麥中,已定位的類病斑突變基因有3個(gè),Li和Bai[14]報(bào)道的小麥品系Ning 7840為自然突變產(chǎn)生的類病斑材料,利用Ning 7840/Chokwang的重組自交系將控制Ning 7840黃色斑點(diǎn)性狀的隱性基因lm定位在小麥1BL染色體上;Yao等[15]在Yanzhan 1/Zaosui 30的分離群體中發(fā)現(xiàn)了受光照誘導(dǎo)的類病斑突變性狀,其實(shí)是兩個(gè)隱性基因互作產(chǎn)生的斑點(diǎn)癥狀,對(duì)控制斑點(diǎn)性狀的這兩個(gè)互作基因進(jìn)行定位,其中來(lái)自親本Zaosui 30的lm1被定位于小麥3BS染色體上Xwmc674與Xbarc133/Xbarc147之間,遺傳距離分別為1.2 cM、3.8 cM,另一個(gè)來(lái)自親本Yanzhan 1的lm2被定位到小麥4B染色體的長(zhǎng)臂,距離兩側(cè)標(biāo)記Xgwm513、Xksum154的遺傳距離分別為1.5 cM、3.0 cM。本研究的類病斑突變體LF2010的突變基因雖然其表達(dá)也受光照誘導(dǎo)[16],但被定位于小麥6BL染色體上,不同于之前定位的任何一個(gè)類病斑基因。因此,為一個(gè)新的類病斑突變基因,將其命名為lm3。

    在基因定位中,為了找到目標(biāo)染色體區(qū)段,人們提出了多種策略,其中應(yīng)用最廣泛的是Michelmore等[24]于1991年提出的BSA法,該方法簡(jiǎn)單方便,能夠快速檢測(cè)到目標(biāo)基因相關(guān)的區(qū)域。但由于小麥SSR遺傳圖譜標(biāo)記密度不足夠高,有些染色體區(qū)段還存在較大的gap(>30 cM),因此該方法在結(jié)合傳統(tǒng)SSR分子標(biāo)記時(shí)也存在一定的弊端,就是靈敏度不高。例如,小麥桿銹基因SrCad的定位,一開(kāi)始采用BSA法定位時(shí)失敗,最終由全基因組定位(whole genome mapping,WGM)法定位成功[25]。但全基因組定位需要消耗大量的人力物力,由此Hiebert等[18]提出了MBSA法,同時(shí)兼?zhèn)淞巳蚪M定位的靈敏度與BSA法的操作高效性。Hiebert等[26]運(yùn)用MBSA法將小麥葉銹抗性基因 Lr70成功定位在小麥5D染色體短臂。突變基因lm3的定位工作,一開(kāi)始單純依靠BSA法進(jìn)行基因定位時(shí),656對(duì)SSR引物在京花一號(hào)與L2010的F2∶3定位群體中只篩選到兩個(gè)連鎖的分子標(biāo)記Xwmc388和cfa2187[23],而在中國(guó)春與LF2010的F2∶3定位群體中cfa2187在雙親間沒(méi)有多態(tài)性,Xwmc388在染色體上有多個(gè)位點(diǎn),由于本實(shí)驗(yàn)室沒(méi)有缺體材料,所以并不能確定目標(biāo)帶型是位于哪一條染色體上?;诖?,我們初步將目標(biāo)染色體鎖定在Xwmc388所在的染色體,運(yùn)用MBSA提高檢測(cè)的靈敏度進(jìn)行連鎖標(biāo)記的篩選,得到兩個(gè)連鎖標(biāo)記,最終確定突變基因的目標(biāo)染色體。說(shuō)明MBSA在小麥這樣基因組龐大、遺傳連鎖圖密度不夠高的作物中進(jìn)行基因定位是切實(shí)可行的。

    隨著生物芯片的不斷發(fā)展,新的快速定位基因的策略應(yīng)運(yùn)而生。相對(duì)于傳統(tǒng)的利用全基因組SSR標(biāo)記篩選目標(biāo)染色體區(qū)段的方法,更加省時(shí)省力、快速精確。同時(shí),90K芯片是通過(guò)不同小麥品系的轉(zhuǎn)錄組數(shù)據(jù)開(kāi)發(fā)而來(lái),檢測(cè)到的SNP大多數(shù)位于編碼區(qū),可以提高定位精度和準(zhǔn)確性。本研究即是利用90K芯片與BSA的結(jié)合,對(duì)兩個(gè)相對(duì)性狀基因池進(jìn)行SNP檢測(cè),根據(jù)SNP的富集結(jié)果快速地發(fā)現(xiàn)了多態(tài)性SNP較多的的集中在6B染色體上,經(jīng)6B富集區(qū)設(shè)計(jì)引物的篩選及F2∶3群體的驗(yàn)證分析確定突變基因lm3在6B染色體上,與MBSA結(jié)果一致。說(shuō)明利用90K SNP芯片快速進(jìn)行小麥基因定位的策略是切實(shí)可行、簡(jiǎn)單高效。需要說(shuō)明的是,本研究中90K芯片在池之間的檢測(cè)結(jié)果在6B以外的個(gè)別染色體上SNPs的數(shù)量也相對(duì)較多,并且得到的SNP富集區(qū)間也相對(duì)較大,推測(cè)是由于建池時(shí)所選取的F2∶3單株數(shù)目較少,使得池之間除目標(biāo)區(qū)域外其他背景值干擾性太強(qiáng)??梢栽诖_保F2單株正確的基礎(chǔ)上選取較多的單株進(jìn)行混池檢測(cè),所得的富集區(qū)間會(huì)大大縮小,更加富有針對(duì)性,有利于開(kāi)發(fā)標(biāo)記加密目標(biāo)基因。后續(xù)將選用lm3兩側(cè)的共顯性標(biāo)記對(duì)更大的F2∶3群體進(jìn)行篩選,得到足夠數(shù)量的重組單株作為篩選加密標(biāo)記的群體,并選取更多數(shù)量的單株混池進(jìn)行基因芯片(90K或660K)的檢測(cè),根據(jù)富集區(qū)的信息開(kāi)發(fā)標(biāo)記對(duì)類病斑突變基因lm3進(jìn)一步加密。

    4 結(jié) 論

    小麥類病斑突變體LF2010的黃色斑點(diǎn)性狀為單基因隱性遺傳。通過(guò)運(yùn)用BSA、MBSA、小麥90K芯片3種不同的基因定位策略將小麥類病斑突變基因lm3定位于小麥6B染色體的長(zhǎng)臂,距兩側(cè)標(biāo)記SWES2、Xbarc134的遺傳距離分別為13.1 cM和3.8 cM,為一個(gè)新的小麥類病斑突變基因。

    [1]FAMBRINI M,CASTAGNA A,DALLA VECCHIA F,etal.Characterization of a pigment-deficient mutant of sunflower (HelianthusannuusL.) with abnormal chloroplast biogenesis,reduced PS II activity and low endogenous level of abscisic acid [J].PlantScience, 2004,167(1):79.

    [2]JUNG K H,HUR J,RYU C H,etal.Characterization of a rice chlorophyll-deficient mutant using the T-DNA gene-trap system [J].PlantandCellPhysiology, 2003,44(5):463.

    [3]SCHWARTZ S H,QIN X,ZEEVAART J A D.Elucidation of the indirect pathway of abscisic acid biosynthesis by mutants,genes,and enzymes [J].PlantPhysiology,2003,131(4):1591.

    [4]JOHAL G S,HULBERT S H,BRIGGS S P.Disease lesion mimics of maize:a model for cell death in plants [J].Bioessays, 1995,17(8):685.

    [5] 孫惠敏,張春嬌,李保同,等.水稻類病斑突變體的研究進(jìn)展[J].上海農(nóng)業(yè)學(xué)報(bào),2014,30(3):142.

    SUN H M,ZHANG C J,LI B T,etal.Advances of study on rice lesion mimic mutants [J].ActaAgricuturaeShanghai, 2014,30(3):142.

    [6] 楊雪靜,劉燕娟,高必達(dá),等.植物類病斑壞死突變體研究進(jìn)展[J].湖南農(nóng)業(yè)科學(xué),2010(23):92.

    YANG X J,LIU Y J,GAO B D,etal.Progresses of the plant lesion mimic mutation [J].HunanAgriculturalSciences, 2010(23):92.

    [7]YIN Z,CHEN J,ZENG L,etal.Characterizing rice lesion mimic mutants and identifying a mutant with broad-spectrum resistance to rice blast and bacterial blight [J].MolecularPlant-MicrobeInteractions, 2000,13(8):869.

    [8]KAMLOFSKI C A,ANTONELLI E,BENDER C,etal.A lesion-mimic mutant of wheat with enhanced resistance to leaf rust [J].PlantPathology, 2007,56(1):46.

    [9]BOYD L A,MINCHIN P N.Wheat mutants showing altered adult plant disease resistance [J].Euphytica, 2001,122(2):361.

    [10]WU C,BORDEOS A,MADAMBA M R S,etal.Rice lesion mimic mutants with enhanced resistance to diseases [J].MolecularGeneticsandGenomics, 2008,279(6):605.

    [11]KINANE J T,JONES P W.Isolation of wheat mutants with increased resistance to powdery mildew from small induced variant populations [J].Euphytica,2001,117(3):251.

    [12]NAIR S K,TOMAR S M S.Genetical and anatomical analyses of a leaf flecking mutant inTriticumaestivumL.[J].Euphytica, 2001,121(1):53.

    [13]LUO P G,REN Z L.Wheat leaf chlorosis controlled by a single recessive gene [J].JournalofPlantPhysiologyandMolecularBiology, 2006,32(3):330.

    [14]LI T,BAI G.Lesion mimic associates with adult plant resistance to leaf rust infection in wheat [J].TheoreticalandAppliedGenetics,2009,119(1):13.

    [15]YAO Q,ZHOU R,FU T,etal.Characterization and mapping of complementary lesion-mimic geneslm1 andlm2 in common wheat [J].TheoreticalandAppliedGenetics, 2009,119(6):1005.

    [16] 杜麗芬,李明飛,劉錄祥,等.一個(gè)化學(xué)誘變的小麥斑點(diǎn)葉突變體的生理和遺傳分析[J].作物學(xué)報(bào),2014,40(6):1020.

    DU L F,LI M F,LIU L X,etal.Physiological characteristics and genetic analysis on a spotted-leaf wheat derived from chemical mutation [J].2014,40(6):1020.

    [17] 白鳳虎,李德芳,陳安國(guó),等.基于BSA分析法的分子標(biāo)記基因定位技術(shù)在農(nóng)作物中的應(yīng)用[J].中國(guó)麻業(yè)科學(xué),2006,28(6):282.

    BAI F H,LI D F,CHEN A G.etal.The application of the molecular gene localization technique based on BSA in crop [J].PlantFiberScuencesinChina, 2006,28(6):282.

    [18]HIEBERT C.An alternative approach for identifying the chromosome location of new genes [J].CanadianJournalofPlantScience,2012,92:597.

    [19]GHAZVINI H,HIEBERT C W,THOMAS J B,etal.Development of a multiple bulked segregant analysis (MBSA) method used to locate a new stem rust resistance gene( Sr54) in the winter wheat cultivar Norin 40 [J].TheoreticalandAppliedGenetics, 2013,126(2):443.

    [20]WANG S,WONG D,FORREST K,etal.Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array [J].PlantBiotechnologyJournal, 2014,12(6):787.

    [21]ROGERS S O,BENDICH A J.Extraction of DNA from milligram amounts of fresh,herbarium and mummified plant tissues [J].PlantMolecularBiology, 1985,5(2):69.

    [22]MICHELMORE R W,PARAN I,KESSELI R.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 [J].ProceedingsoftheNationalAcademyofSciences, 1991,88(21):9828.

    [23] 杜麗芬.小麥斑點(diǎn)葉突變體LF2010的生理與遺傳定位研究[D] .楊凌:西北農(nóng)林科技大學(xué),2014:25.

    DU L F.Study on physiological characteristics and genetic mapping of a spotted-leaf wheat mutant LF2010 [D].Yangling:Northwest A&F University,2014:25.

    [24]SOURDILLE P,SINGH S,CADALEN T,etal.Microsatellite-based deletion bin system for the establishment of genetic-physical map relationships in wheat(TriticumaestivumL.) [J].Functional&IntegrativeGenomics,2004,4(1):12.

    [25]HIEBERT C W,FETCH T G,ZEGEYE T,etal.Genetics and mapping of seedling resistance to Ug99 stem rust in Canadian wheat cultivars ‘Peace’ and ‘AC Cadillac’ [J].TheoreticalandAppliedGenetics, 2011,122(1):143.

    [26]HIEBERT C W,MCCALLUM B D,THOMAS J B. Lr70,a new gene for leaf rust resistance mapped in common wheat accession KU3198 [J].TheoreticalandAppliedGenetics, 2014,127(9):2005.

    Mapping of Lesion Mimic Mutant Genelm3 in Common Wheat

    YANG Jiaxiu1,GENG Jiefei1,LI Qianqian1,LIU Luxiang2,XIE Yanzhou1,WANG Chengshe1

    (1.College of Agronomy,Northwest A&F University,Yangling,Shaanxi 712100,China; 2.Institute of Crop Sciences,Chinese Academy of Agricultural Sciences,Beijing 100081,China)

    In order to clone the lesion mimic mutant genelm3 and analyse its function further,a lesion mimic wheat mutant LF2010,obtained from an anther culture line H261 which was mutated with EMS,was chosen as material,and three kinds of different gene mapping strategies were taken to map it,including BSA(bulked segregant analysis),MBSA(multiple bulked segregant analysis),and wheat 90K gene chip combined with BSA.The result showed that a genetic linkage map containing 15 pairs markers linked withlm3 was obtained. Finally,lm3 was mapped on 6BL of wheat chromosome,13.8 cM away from marker SWES2 and 3.8 cM away from Xbarc134.lm3 is a new gene of wheat lesion mimic mutant.

    Wheat; Lesion mimic mutant; Gene localization;lm3

    時(shí)間:2016-12-07

    2016-05-26

    2016-11-13

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFD0102101);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2014YB081);陜西省科技統(tǒng)籌創(chuàng)新工程計(jì)劃項(xiàng)目(2014KTZB02-01-01,2015KTZDNY01-01-02);西北農(nóng)林科技大學(xué)南陽(yáng)小麥試驗(yàn)示范站建設(shè)項(xiàng)目;西北農(nóng)林科技大學(xué)唐仲英育種基金項(xiàng)目

    E-mail:jiaxiu_0906@126.com(楊佳秀);18700895579@163.com(耿皆飛,與第一作者同等貢獻(xiàn))

    謝彥周(E-mail:xieyanzhou@gmail.com); 王成社(E-mail:wangcs2008@126.com)

    S512.1;S330

    A

    1009-1041(2016)12-1578-09

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1359.S.20161207.1748.006.html

    猜你喜歡
    親本突變體病斑
    甘蔗親本農(nóng)藝性狀評(píng)價(jià)與分析
    改進(jìn)自適應(yīng)分水嶺方法分割棉花葉部粘連病斑
    幾種蘋(píng)果砧木實(shí)生后代與親本性狀的相關(guān)性
    一種桑輪紋病病葉發(fā)病程度的評(píng)價(jià)指數(shù)
    9種藥劑防治山核桃干腐病試驗(yàn)研究
    CLIC1及其點(diǎn)突變體與Sedlin蛋白的共定位研究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    云瑞10系列生產(chǎn)性創(chuàng)新親本2種方法評(píng)價(jià)
    Survivin D53A突變體對(duì)宮頸癌細(xì)胞增殖和凋亡的影響
    油用向日葵親本繁殖及雜交制種技術(shù)
    别揉我奶头~嗯~啊~动态视频 | 久久鲁丝午夜福利片| 久久精品国产亚洲av高清一级| 亚洲天堂av无毛| 夫妻午夜视频| 久久人人爽人人片av| 在线观看人妻少妇| 亚洲中文av在线| 制服诱惑二区| 久久青草综合色| 建设人人有责人人尽责人人享有的| 国产野战对白在线观看| 欧美激情极品国产一区二区三区| 日韩不卡一区二区三区视频在线| 亚洲色图综合在线观看| 晚上一个人看的免费电影| 欧美黄色片欧美黄色片| 日本色播在线视频| 国产片内射在线| 9热在线视频观看99| 熟妇人妻不卡中文字幕| 最近中文字幕2019免费版| 日本一区二区免费在线视频| 9191精品国产免费久久| 一级爰片在线观看| 欧美精品一区二区大全| 久久久精品94久久精品| 最近2019中文字幕mv第一页| 免费看不卡的av| 母亲3免费完整高清在线观看| 国产精品人妻久久久影院| 女人爽到高潮嗷嗷叫在线视频| 老熟女久久久| 国产熟女欧美一区二区| 中文精品一卡2卡3卡4更新| 国产在视频线精品| avwww免费| 国产福利在线免费观看视频| 亚洲精品国产色婷婷电影| 国产黄色视频一区二区在线观看| 黄片小视频在线播放| 日韩一本色道免费dvd| 亚洲精品国产一区二区精华液| 国产精品熟女久久久久浪| 免费黄频网站在线观看国产| 亚洲第一区二区三区不卡| 成人国产av品久久久| 99国产精品免费福利视频| 国产精品偷伦视频观看了| 制服人妻中文乱码| 丝袜在线中文字幕| 亚洲人成77777在线视频| 亚洲欧美成人综合另类久久久| 高清不卡的av网站| videos熟女内射| 天天添夜夜摸| 亚洲伊人久久精品综合| 久久精品亚洲熟妇少妇任你| 国产色婷婷99| 亚洲精品久久久久久婷婷小说| 成人影院久久| 在线亚洲精品国产二区图片欧美| 黑丝袜美女国产一区| 久久久精品国产亚洲av高清涩受| 无限看片的www在线观看| 少妇精品久久久久久久| 日韩人妻精品一区2区三区| 国产精品 欧美亚洲| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 老鸭窝网址在线观看| 亚洲国产中文字幕在线视频| 国产无遮挡羞羞视频在线观看| 亚洲,欧美,日韩| 欧美日韩视频精品一区| 好男人视频免费观看在线| 男女午夜视频在线观看| a级毛片在线看网站| 男女高潮啪啪啪动态图| 国产精品国产三级国产专区5o| kizo精华| 欧美人与性动交α欧美精品济南到| 国产av国产精品国产| 18禁观看日本| 中文字幕人妻丝袜制服| 亚洲欧洲精品一区二区精品久久久 | 亚洲成人免费av在线播放| 叶爱在线成人免费视频播放| 亚洲一级一片aⅴ在线观看| 我要看黄色一级片免费的| 日本欧美视频一区| 日韩大码丰满熟妇| 亚洲,一卡二卡三卡| 伊人久久大香线蕉亚洲五| 国产精品一区二区在线观看99| 桃花免费在线播放| 久久久久精品性色| 91精品三级在线观看| 亚洲天堂av无毛| 一二三四中文在线观看免费高清| 丝瓜视频免费看黄片| 黑丝袜美女国产一区| 久久久久视频综合| 国产亚洲av高清不卡| 国产伦理片在线播放av一区| 精品一区二区三区四区五区乱码 | 一二三四在线观看免费中文在| 欧美老熟妇乱子伦牲交| 人人妻人人澡人人看| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美色中文字幕在线| 成人毛片60女人毛片免费| 国产精品 欧美亚洲| 建设人人有责人人尽责人人享有的| 亚洲精品美女久久久久99蜜臀 | 热re99久久国产66热| 高清在线视频一区二区三区| 永久免费av网站大全| 男女高潮啪啪啪动态图| 自线自在国产av| 国产av码专区亚洲av| 中文精品一卡2卡3卡4更新| 狂野欧美激情性bbbbbb| 国产精品.久久久| 国产成人91sexporn| 日本欧美国产在线视频| 欧美人与善性xxx| 国产国语露脸激情在线看| 一本久久精品| 熟女av电影| 亚洲精华国产精华液的使用体验| 久久精品久久久久久噜噜老黄| netflix在线观看网站| 免费少妇av软件| 视频在线观看一区二区三区| 男男h啪啪无遮挡| 在线观看人妻少妇| 在线观看免费日韩欧美大片| 日韩欧美精品免费久久| 热99久久久久精品小说推荐| 精品一区二区三区四区五区乱码 | 亚洲成人一二三区av| 精品酒店卫生间| 十分钟在线观看高清视频www| 欧美日韩av久久| 欧美激情 高清一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线 av 中文字幕| 国产精品二区激情视频| 久久久久视频综合| 亚洲熟女毛片儿| 国产在线免费精品| www.精华液| 国产一级毛片在线| 国产成人精品久久二区二区91 | 在线观看人妻少妇| 国产精品一国产av| 久久精品亚洲熟妇少妇任你| 日韩,欧美,国产一区二区三区| 国产一区二区在线观看av| 久久久久网色| 五月开心婷婷网| 成年女人毛片免费观看观看9 | 国产在线一区二区三区精| 久久99一区二区三区| 久久影院123| 日韩制服丝袜自拍偷拍| 汤姆久久久久久久影院中文字幕| 一区二区三区乱码不卡18| 人妻 亚洲 视频| 99精品久久久久人妻精品| 国语对白做爰xxxⅹ性视频网站| 午夜免费观看性视频| 国产 一区精品| 久久精品国产综合久久久| h视频一区二区三区| 欧美精品一区二区免费开放| 99久国产av精品国产电影| 宅男免费午夜| 成年女人毛片免费观看观看9 | 国产亚洲av片在线观看秒播厂| 母亲3免费完整高清在线观看| 久久青草综合色| 午夜影院在线不卡| 欧美精品一区二区免费开放| 国产精品久久久久久精品电影小说| 亚洲久久久国产精品| 国产成人免费观看mmmm| 丝袜在线中文字幕| 麻豆av在线久日| 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 久久精品久久久久久噜噜老黄| 亚洲精品日本国产第一区| 亚洲精品美女久久久久99蜜臀 | 久久久久网色| 久久国产精品大桥未久av| 免费av中文字幕在线| 中文字幕制服av| 亚洲欧美一区二区三区黑人| 青青草视频在线视频观看| tube8黄色片| 免费高清在线观看视频在线观看| 老司机影院成人| 不卡视频在线观看欧美| 男女边吃奶边做爰视频| 午夜福利视频在线观看免费| 国产欧美日韩综合在线一区二区| 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 中文精品一卡2卡3卡4更新| 蜜桃国产av成人99| 色精品久久人妻99蜜桃| 免费在线观看视频国产中文字幕亚洲 | 国产精品国产三级国产专区5o| 天天添夜夜摸| 亚洲av电影在线观看一区二区三区| www.自偷自拍.com| 日本av免费视频播放| 夜夜骑夜夜射夜夜干| 欧美黑人欧美精品刺激| 欧美国产精品一级二级三级| 91成人精品电影| 两性夫妻黄色片| 一本大道久久a久久精品| 人妻一区二区av| 亚洲av男天堂| 国产欧美日韩一区二区三区在线| 婷婷色麻豆天堂久久| 亚洲欧美一区二区三区久久| 熟女av电影| 亚洲精品乱久久久久久| 丝瓜视频免费看黄片| 午夜日本视频在线| av电影中文网址| 亚洲成av片中文字幕在线观看| 又粗又硬又长又爽又黄的视频| 丝袜在线中文字幕| 亚洲精品视频女| 欧美日韩国产mv在线观看视频| 免费看不卡的av| 一本色道久久久久久精品综合| 国产精品久久久久久精品古装| 在线看a的网站| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 成人国产av品久久久| 午夜av观看不卡| 国产精品熟女久久久久浪| 国产探花极品一区二区| 黑人猛操日本美女一级片| 亚洲图色成人| 69精品国产乱码久久久| 国产精品一二三区在线看| 免费观看av网站的网址| 国产成人a∨麻豆精品| 两个人免费观看高清视频| 亚洲欧洲国产日韩| 国产亚洲av高清不卡| 免费观看av网站的网址| 咕卡用的链子| 久久久久精品人妻al黑| 国产人伦9x9x在线观看| 久久精品aⅴ一区二区三区四区| 极品人妻少妇av视频| 一边摸一边做爽爽视频免费| 人成视频在线观看免费观看| 夜夜骑夜夜射夜夜干| 岛国毛片在线播放| 综合色丁香网| 久久久欧美国产精品| 一本色道久久久久久精品综合| 中文字幕制服av| 汤姆久久久久久久影院中文字幕| 国产黄色视频一区二区在线观看| 欧美xxⅹ黑人| 咕卡用的链子| 夫妻性生交免费视频一级片| 国精品久久久久久国模美| 国产成人午夜福利电影在线观看| 亚洲欧美激情在线| 纵有疾风起免费观看全集完整版| 国产精品久久久久成人av| 看免费成人av毛片| 久久午夜综合久久蜜桃| 久久婷婷青草| 国产97色在线日韩免费| 国产爽快片一区二区三区| 亚洲人成电影观看| 一级黄片播放器| 欧美日韩一区二区视频在线观看视频在线| 久久97久久精品| 国产一卡二卡三卡精品 | 伊人久久国产一区二区| bbb黄色大片| 亚洲av综合色区一区| 国产精品 国内视频| 99久国产av精品国产电影| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 国产av一区二区精品久久| 国产精品亚洲av一区麻豆 | 大片电影免费在线观看免费| tube8黄色片| 国产精品秋霞免费鲁丝片| 999精品在线视频| 亚洲精品视频女| 美女国产高潮福利片在线看| 国产成人精品无人区| 国产在视频线精品| 欧美日韩一区二区视频在线观看视频在线| 日韩免费高清中文字幕av| 午夜福利免费观看在线| 美女福利国产在线| 伦理电影大哥的女人| xxxhd国产人妻xxx| 亚洲久久久国产精品| 国产视频首页在线观看| 国产精品久久久av美女十八| 狂野欧美激情性xxxx| 女人高潮潮喷娇喘18禁视频| av片东京热男人的天堂| 精品人妻一区二区三区麻豆| 女人精品久久久久毛片| 无遮挡黄片免费观看| 性少妇av在线| 成年av动漫网址| 哪个播放器可以免费观看大片| 国产精品国产三级国产专区5o| 亚洲国产欧美一区二区综合| 黄色视频在线播放观看不卡| 久久免费观看电影| 久热爱精品视频在线9| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 1024视频免费在线观看| 国产精品国产三级国产专区5o| 日韩av不卡免费在线播放| 国产片特级美女逼逼视频| 亚洲第一区二区三区不卡| 一区二区三区乱码不卡18| 免费高清在线观看日韩| 一本久久精品| 制服丝袜香蕉在线| 一本大道久久a久久精品| 叶爱在线成人免费视频播放| 亚洲自偷自拍图片 自拍| 最近中文字幕2019免费版| 精品国产国语对白av| 黄色一级大片看看| 男男h啪啪无遮挡| 成人免费观看视频高清| 91精品伊人久久大香线蕉| av有码第一页| 亚洲欧美一区二区三区久久| 精品国产乱码久久久久久男人| 99九九在线精品视频| 久久国产精品大桥未久av| 亚洲婷婷狠狠爱综合网| 一级毛片黄色毛片免费观看视频| 中文字幕精品免费在线观看视频| 满18在线观看网站| 97人妻天天添夜夜摸| 伊人亚洲综合成人网| 欧美国产精品va在线观看不卡| 欧美日韩亚洲综合一区二区三区_| 亚洲成人手机| 国产av精品麻豆| 美国免费a级毛片| 叶爱在线成人免费视频播放| 国产精品熟女久久久久浪| 久久精品国产亚洲av高清一级| 日本欧美国产在线视频| 久久久久久久久久久久大奶| 国产精品国产三级国产专区5o| 欧美亚洲 丝袜 人妻 在线| 国产乱来视频区| 亚洲精品久久久久久婷婷小说| 国产乱来视频区| 国产精品成人在线| 在线亚洲精品国产二区图片欧美| 精品一区二区免费观看| 午夜福利视频在线观看免费| 亚洲精品国产一区二区精华液| 日韩一区二区视频免费看| 国产成人啪精品午夜网站| 欧美黑人精品巨大| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| 日日爽夜夜爽网站| 最近中文字幕2019免费版| 一本久久精品| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 国产精品偷伦视频观看了| 婷婷成人精品国产| 少妇人妻 视频| 97在线人人人人妻| 午夜激情久久久久久久| 成人午夜精彩视频在线观看| 久久午夜综合久久蜜桃| 国产探花极品一区二区| 精品卡一卡二卡四卡免费| 国产高清不卡午夜福利| 日本欧美国产在线视频| 亚洲精品自拍成人| 国产一级毛片在线| 亚洲伊人久久精品综合| 搡老岳熟女国产| 久久久国产一区二区| 久热这里只有精品99| 男女高潮啪啪啪动态图| 十八禁网站网址无遮挡| 男男h啪啪无遮挡| 精品国产一区二区三区久久久樱花| 美国免费a级毛片| 在线天堂最新版资源| 亚洲自偷自拍图片 自拍| 交换朋友夫妻互换小说| 国产精品一区二区精品视频观看| 久久性视频一级片| 熟妇人妻不卡中文字幕| 日韩中文字幕视频在线看片| 赤兔流量卡办理| 91老司机精品| 伦理电影大哥的女人| 精品久久久精品久久久| 天天操日日干夜夜撸| 亚洲欧美精品综合一区二区三区| 亚洲美女黄色视频免费看| 丝袜在线中文字幕| 桃花免费在线播放| 人妻 亚洲 视频| 久久久久久久久免费视频了| 亚洲av电影在线进入| 欧美人与性动交α欧美软件| 亚洲精品aⅴ在线观看| 欧美日韩av久久| av在线观看视频网站免费| 亚洲精品国产区一区二| 日韩一卡2卡3卡4卡2021年| 在线观看人妻少妇| 久久久精品94久久精品| 捣出白浆h1v1| 国产在视频线精品| 亚洲精品自拍成人| 欧美 日韩 精品 国产| av片东京热男人的天堂| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品久久久久久婷婷小说| 一级黄片播放器| 99精品久久久久人妻精品| 国产日韩欧美亚洲二区| 一级毛片电影观看| 看免费成人av毛片| 国产精品三级大全| 国产午夜精品一二区理论片| 黄色怎么调成土黄色| 久久国产精品大桥未久av| www日本在线高清视频| 亚洲精品日韩在线中文字幕| 狠狠精品人妻久久久久久综合| 国产成人91sexporn| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区三区| 免费高清在线观看日韩| 亚洲一区中文字幕在线| 亚洲美女黄色视频免费看| 亚洲av成人不卡在线观看播放网 | 日韩精品有码人妻一区| 日韩视频在线欧美| 久久久久人妻精品一区果冻| 久久久久视频综合| 欧美激情 高清一区二区三区| 少妇猛男粗大的猛烈进出视频| 在线天堂最新版资源| 国产精品免费视频内射| 看非洲黑人一级黄片| av在线老鸭窝| 美女扒开内裤让男人捅视频| 捣出白浆h1v1| 亚洲四区av| 久久99精品国语久久久| av线在线观看网站| 美女扒开内裤让男人捅视频| 亚洲国产日韩一区二区| 无限看片的www在线观看| 亚洲av国产av综合av卡| 亚洲精品国产av蜜桃| xxx大片免费视频| 亚洲国产精品999| 国产免费一区二区三区四区乱码| 老鸭窝网址在线观看| 欧美成人精品欧美一级黄| 成人国产av品久久久| 久久ye,这里只有精品| 伦理电影大哥的女人| 女性被躁到高潮视频| 91国产中文字幕| 十八禁高潮呻吟视频| svipshipincom国产片| 九草在线视频观看| 亚洲精品日本国产第一区| 香蕉国产在线看| 日本91视频免费播放| 男女国产视频网站| 国产在线免费精品| 久久精品亚洲熟妇少妇任你| 制服诱惑二区| 亚洲第一av免费看| 综合色丁香网| 亚洲欧美精品自产自拍| 亚洲av综合色区一区| 18禁裸乳无遮挡动漫免费视频| 亚洲成人国产一区在线观看 | videos熟女内射| 成人影院久久| 国产av精品麻豆| 亚洲五月色婷婷综合| 久久久久精品久久久久真实原创| 秋霞伦理黄片| 91精品国产国语对白视频| 精品福利永久在线观看| 999久久久国产精品视频| 一边亲一边摸免费视频| 19禁男女啪啪无遮挡网站| 久久久久网色| 在线观看免费午夜福利视频| 久久婷婷青草| 国产成人a∨麻豆精品| 亚洲国产欧美一区二区综合| 亚洲美女搞黄在线观看| 一级毛片电影观看| 国产高清不卡午夜福利| 国产精品国产三级国产专区5o| 亚洲av男天堂| 国产成人一区二区在线| 国产精品久久久久成人av| 香蕉国产在线看| 日本欧美视频一区| 国产av国产精品国产| 中文天堂在线官网| 免费黄网站久久成人精品| 免费av中文字幕在线| 国产成人精品在线电影| 亚洲精品第二区| av片东京热男人的天堂| 人妻一区二区av| 不卡av一区二区三区| 国产熟女午夜一区二区三区| 久久婷婷青草| 午夜日韩欧美国产| 国产极品天堂在线| 少妇精品久久久久久久| av国产精品久久久久影院| 男女下面插进去视频免费观看| 超碰97精品在线观看| 中文精品一卡2卡3卡4更新| www.av在线官网国产| 国产 精品1| 久久这里只有精品19| 亚洲欧美成人精品一区二区| 中文字幕精品免费在线观看视频| 国产 一区精品| avwww免费| 午夜日韩欧美国产| 亚洲精品,欧美精品| 亚洲第一av免费看| 黄色毛片三级朝国网站| 午夜日本视频在线| 午夜福利视频在线观看免费| 久久青草综合色| 国产精品 国内视频| 97在线人人人人妻| 综合色丁香网| av网站免费在线观看视频| 久久久久久人妻| 国产高清国产精品国产三级| 中文字幕人妻丝袜制服| 久久久久久人人人人人| 亚洲欧洲日产国产| 大陆偷拍与自拍| 婷婷成人精品国产| 欧美中文综合在线视频| 欧美日韩综合久久久久久| 国产精品一二三区在线看| 国产男人的电影天堂91| 国产成人精品福利久久| 精品少妇一区二区三区视频日本电影 | netflix在线观看网站| 两性夫妻黄色片| avwww免费| 欧美日韩国产mv在线观看视频| 香蕉国产在线看| h视频一区二区三区| 欧美乱码精品一区二区三区| 国产精品成人在线| 一区福利在线观看| 日韩欧美精品免费久久| 日本午夜av视频| av视频免费观看在线观看| 无限看片的www在线观看| 欧美日韩福利视频一区二区| 热99久久久久精品小说推荐| av不卡在线播放| 国产欧美日韩综合在线一区二区| 1024香蕉在线观看| 在线 av 中文字幕| 日韩一卡2卡3卡4卡2021年| 精品人妻熟女毛片av久久网站| 国产精品免费视频内射| 少妇的丰满在线观看| 久久久国产精品麻豆| 在线观看www视频免费|