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

    利用BSA-Seq方法快速定位作物農(nóng)藝性狀QTL/基因概述

    2022-05-10 23:32:52周文期劉忠祥王曉娟何海軍周玉乾楊彥忠連曉榮李永生
    甘肅農(nóng)業(yè)科技 2022年4期
    關(guān)鍵詞:基因克隆農(nóng)藝性狀

    周文期 劉忠祥 王曉娟 何海軍 周玉乾 楊彥忠 連曉榮 李永生

    摘要:為了利用生物學(xué)技術(shù)和手段來(lái)解決農(nóng)業(yè)問(wèn)題,讓基因組測(cè)序技術(shù)、轉(zhuǎn)基因及基因編輯等現(xiàn)代分子生物學(xué)技術(shù)為農(nóng)業(yè)遺傳育種服務(wù),實(shí)現(xiàn)常規(guī)育種與分子設(shè)計(jì)育種的緊密結(jié)合,加速作物精準(zhǔn)育種進(jìn)程,近年頗受學(xué)術(shù)界關(guān)注。現(xiàn)梳理了一些分子生物學(xué)專(zhuān)用名詞概念,概述了目前在正向遺傳學(xué)研究中如何利用極端表型材料基于全基因組測(cè)序的BSA-Seq方法(MutMap法、MutMap+法、MutMap-Gap法、QTL-seq法)快速定位重要農(nóng)藝性狀的QTL/基因,為快速定位和克隆候選基因提供新思路。

    關(guān)鍵詞:QTL定位; BSA-Seq方法;遺傳群體;基因克隆;農(nóng)藝性狀;分子育種

    中圖分類(lèi)號(hào):S336;Q78 ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼:A ? ? ? ? ? ?文章編號(hào):1001-1463(2022)04-0001-10

    doi:10.3969/j.issn.1001-1463.2022.04.001

    Rapid Mapping of QTL/gene for Agronomic Traits in Crops Using BSA-Seq Method

    ZHOU Wenqi, LIU Zhongxiang, WANG Xiaojuan, HE Haijun, ZHOU Yuqian, YANG Yanzhong, LIAN Xiaorong, LI Yongsheng

    (Institute of Crops, Gansu Academy of Agricultural Sciences, Lanzhou Gansu 730070, China)

    Abstract:In order to resolve the agriculture problem using biological technologies and tools, modern molecular biology technologies such as genome sequencing, transgenic technology and gene editing could serve agricultural genetics and breeding, to realize the close intergration of traditional breeding and molecular design breeding, and accelerate the process of precision breeding of crops, which has attracted academic attention in recent years. In this study, we sorted out some molecular biological terminology concepts, and summarized how to rapidly locate QTL/gene for important agronomic traits using extreme phenotypic materials and BSA-Seq (Bulked segregate analysis-sequencing) method (MutMap, MutMap+, MutMap-Gap and QTL-seq method) based on the whole genome sequencing in forward genetics, providing new ideas for rapid mapping and cloning of candidate genes.

    Key words:QTL mapping;BSA-seqmethod;Genetical population;Map-based cloning;Agronomic traits;Molecular breeding

    過(guò)去人們依靠傳統(tǒng)育種方法獲取農(nóng)作物新品種,根據(jù)所需新品種的特性,輪回選擇并純化其親本,然后通過(guò)雜交、回交或直系篩選程序來(lái)完成。傳統(tǒng)育種在整個(gè)農(nóng)業(yè)發(fā)展中發(fā)揮了至關(guān)重要的作用,但同時(shí)也存在著育種上的盲目性、經(jīng)驗(yàn)性、長(zhǎng)周期和不確定因素等諸多問(wèn)題。目前,遺傳育種已經(jīng)歷了從原始馴化選育階段由耕作者選育具有廣泛表型變異的地方農(nóng)家品種(育種1.0版);常規(guī)育種階段由職業(yè)育種家通過(guò)預(yù)先設(shè)計(jì)雜交育種試驗(yàn)選育現(xiàn)代栽培品種(育種2.0版);進(jìn)入依靠生物技術(shù)的分子育種階段,采用分子標(biāo)記輔助選擇及轉(zhuǎn)基因技術(shù)的手段實(shí)現(xiàn)對(duì)品種單一目標(biāo)性狀的導(dǎo)入與修飾(育種3.0版)。隨著基因工程和現(xiàn)代信息技術(shù)的迅速崛起,育種也逐漸步入互聯(lián)網(wǎng)+大數(shù)據(jù)+人工智能“三位一體”的設(shè)計(jì)育種時(shí)代,利用分子生物學(xué)技術(shù)和遺傳學(xué)方法控制動(dòng)植物群體數(shù)量和質(zhì)量均得以實(shí)現(xiàn),從而將全面改寫(xiě)生物育種的策略和方法,推動(dòng)育種技術(shù)向智能化階段(育種4.0版)發(fā)展[1 - 2 ]。生物技術(shù)向人類(lèi)展示了巨大發(fā)展?jié)摿?,?chuàng)造著農(nóng)業(yè)革命的未來(lái)。利用分子設(shè)計(jì)育種彰顯出比傳統(tǒng)育種更為突出的優(yōu)越性,結(jié)合基因組學(xué)、表型組學(xué)、蛋白組學(xué)及代謝組學(xué)等多組學(xué)數(shù)據(jù)進(jìn)行生物信息學(xué)的整合和分析、優(yōu)化和篩選,獲取最佳育種目標(biāo)基因型,通過(guò)轉(zhuǎn)基因技術(shù)或基因編輯手段,可高效精準(zhǔn)地培育出目標(biāo)新品種,縮短育種周期2~5 a,提高育種效率,因此智能化育種是未來(lái)作物育種的必然選擇,其精準(zhǔn)性、科學(xué)性和高效性都將帶領(lǐng)作物育種進(jìn)入一個(gè)新的時(shí)代[2 - 3 ]。智能化育種或分子設(shè)計(jì)育種4.0版的關(guān)鍵方案就是如何合理設(shè)計(jì)具有理想株型(玉米、小麥和水稻等有各自不同的理想株型,如玉米合理株高、穗上葉夾角減小、增強(qiáng)耐密性、株型性狀之間的協(xié)同調(diào)控機(jī)制等)、高產(chǎn)、高抗、綠色、優(yōu)質(zhì)等綜合農(nóng)藝性狀的優(yōu)良品種。為了實(shí)現(xiàn)設(shè)計(jì)目標(biāo),定位和克隆控制農(nóng)藝性狀的關(guān)鍵基因,及解析其遺傳調(diào)控機(jī)理至關(guān)重要[4 - 7 ]。通過(guò)正向遺傳方法的圖位克隆技術(shù)是傳統(tǒng)經(jīng)典的基因克隆方法,但需要構(gòu)建永久遺傳群體或F2 ∶ 3分離群體,周期很長(zhǎng),存在高密度分子標(biāo)記篩選工作量大、標(biāo)記難以獲得等缺點(diǎn)。我們通過(guò)介紹如何利用作物極端表型材料基于混合分組測(cè)序(Bulked segregate analysis-sequencing BSA-Seq)的MutMap方法快速獲得重要農(nóng)藝性狀的QTL/基因,旨在為定位和克隆候選基因提供思路和方法。

    1 ? 分子生物學(xué)一些專(zhuān)用名詞介紹

    1.1 ? 基因測(cè)序及克隆相關(guān)概念

    1.1.1 ? 表型與基因型(Phenotype andgenotype) ? 又稱(chēng)性狀,是基因型(Genetype)和環(huán)境(Environment)共同作用的結(jié)果,即P=G+E;基因型是指某一生物個(gè)體全部基因組合的總稱(chēng)[8 ]。

    1.1.2 ? 質(zhì)量性狀(Qualitativecharacter) ? 指屬性性狀,是指同一種性狀的不同表現(xiàn)型之間不存在連續(xù)性的數(shù)量變化,而呈現(xiàn)質(zhì)的中斷性變化的那些性狀,組間差異顯著。質(zhì)量性狀通常受一個(gè)或少數(shù)幾個(gè)主效基因調(diào)控,不易受環(huán)境因素影響。比如花色、抗病性、抗蟲(chóng)性、育性等表現(xiàn)為質(zhì)量性狀遺傳的特點(diǎn)[8 ]。

    1.1.3 ? 數(shù)量性狀(Quantitative character) ? 指某一性狀在群體中變異呈連續(xù)性正態(tài)分布,無(wú)法明確分組,如動(dòng)植物的高度、長(zhǎng)度、成熟期、產(chǎn)量等經(jīng)濟(jì)特性,這類(lèi)性狀無(wú)法通過(guò)表型來(lái)推斷其基因型,易受環(huán)境影響和遺傳背景影響[8 ]。

    1.1.4 ? ?數(shù)量性狀基因座(Quantitative trait locus,QTL) ? 指控制數(shù)量性狀的基因在基因組中的位置。把能控制數(shù)量性狀的位點(diǎn)就叫做QTL[8 ]。

    1.1.5 ? 遺傳標(biāo)記(Genetic marker) ? ?指可追蹤染色體、染色體某一片段或某個(gè)基因座在家系中傳遞的任何一種遺傳特性。具有兩個(gè)基本特征,即可遺傳性和可識(shí)別性。某種生物的任何有差異表型的基因突變型均可作為遺傳標(biāo)記[8 ]。

    1.1.6 ? 遺傳連鎖圖譜(Genetic linkage map) ? 指基因或DNA標(biāo)記在染色體上的相對(duì)位置與遺傳距離。單位由基因或DNA片段在染色體交換過(guò)程中分離的頻率厘摩(cM)來(lái)表示。1 cM表示每次減數(shù)分裂的重組頻率為1%[8 ]。

    1.1.7 ? ?圖位克?。∕ap-based cloning) ? ?又稱(chēng)定位克隆,1986年首先由劍橋大學(xué)的Alan coulson提出[9 ],可用于分離大部分性狀的相關(guān)基因,它是通過(guò)分析突變位點(diǎn)與已知分子標(biāo)記的連鎖關(guān)系來(lái)確定候選基因,根據(jù)目的基因在染色體上的位置進(jìn)行定位[8 ]。群體中的個(gè)體在某一表型性狀上有遺傳差異,導(dǎo)致這些差異的基因可以被已知的分子標(biāo)記定位在特定染色體位置。隨著各種分子標(biāo)記技術(shù)的發(fā)展和高密度連鎖圖譜構(gòu)建完成,圖位克隆技術(shù)已經(jīng)成為分離生物基因最常規(guī)的手段之一。

    1.1.8 ? 全基因組重測(cè)序(Genome-wide resequencing)

    全基因組重測(cè)序是對(duì)基因組序列已知的個(gè)體進(jìn)行基因組水平測(cè)序,并在個(gè)體或群體水平上進(jìn)行差異性分析的方法。如擬南芥兩個(gè)已測(cè)序亞種Columbia和Landsberg、水稻粳稻和秈稻兩個(gè)亞種代表日本晴和9311、玉米R(shí)eid 和Lancaster類(lèi)群代表自交系B73和Mo17等,均已經(jīng)有參考的基因組序列。對(duì)有優(yōu)異農(nóng)藝性狀的個(gè)體或者群體進(jìn)行基因組重測(cè)序,將結(jié)果與已有的基因序列進(jìn)行比對(duì)分析,找出控制重要農(nóng)業(yè)性狀的QTLs或目標(biāo)基因,為育種精準(zhǔn)控制表型所用[8 ]。

    1.1.9 ? 簡(jiǎn)化基因組測(cè)序(Reduced-representation genome sequencing,RRGS) ? 指利用限制性?xún)?nèi)切酶分割基因組DNA,選擇性地回收一定長(zhǎng)度范圍內(nèi)的酶切片段進(jìn)行高通量測(cè)序,獲得海量遺傳多態(tài)性標(biāo)簽序列來(lái)充分代表目標(biāo)物種全基因組信息的測(cè)序策略[8 ]。

    1.1.10 ? 混合分組測(cè)序分析法(Bulked segregate analysis-sequencing,BSA-Seq) ? 又稱(chēng)分離體分組混合分析法或集團(tuán)分離分析法[10 ],常用于單個(gè)QTL、突變位點(diǎn)定位等。將目標(biāo)性狀在F2后代或RILs群體中的極端表型,如株高或株矮、葉寬或葉窄、抗逆性狀表現(xiàn)高抗或敏感等中2組個(gè)體的DNA分別混合成2個(gè)DNA池,測(cè)序比較兩組群體在多態(tài)位點(diǎn)(SNPs)的等位基因頻率(AF)是否具有顯著差異,然后利用高密度分子標(biāo)記在兩池中進(jìn)行標(biāo)記與性狀間的共分離分析,通過(guò)SNP值計(jì)算,篩選與性狀相關(guān)聯(lián)的QTL/基因[8 ]。

    1.1.11 ? RNA的混合分組測(cè)序法(Bulked segregant RNA-sequencing,BSR-Seq) ? 指將轉(zhuǎn)錄組測(cè)序與集群分離分析相結(jié)合,在轉(zhuǎn)錄組范圍內(nèi)開(kāi)發(fā)SNPs,篩選與性狀緊密連鎖的SNPs,進(jìn)行功能基因的定位,同時(shí)進(jìn)行基因差異表達(dá)分析等轉(zhuǎn)錄組常規(guī)分析的技術(shù)?;赗NA表達(dá)水平的BSR,更多的是轉(zhuǎn)錄本的比對(duì)[8 ]。

    1.1.12 ? 全基因組關(guān)聯(lián)分析(Genome-wide association study,GWAS) ? 是一種對(duì)全基因組范圍內(nèi)的常見(jiàn)遺傳變異(單核苷酸多態(tài)性和拷貝數(shù))多態(tài)性總體關(guān)聯(lián)分析的方法[11 ],在全基因組范圍內(nèi)進(jìn)行整體研究,能夠一次性對(duì)目標(biāo)性狀進(jìn)行輪廓性概覽,適用于復(fù)雜性狀的研究[8 ]。

    1.2 ? 遺傳群體定義及構(gòu)建方法

    1.2.1 ? 作圖群體(Mapping population) ? 指雙親通過(guò)自交和雜交產(chǎn)生后代,后代中性狀和標(biāo)記位點(diǎn)基因型均是分離的群體,也叫分離群體,常被用于遺傳圖譜構(gòu)建,因此也叫作圖群體[8 ]。作圖群體被分為初級(jí)作圖群體和次級(jí)作圖群體2種。初級(jí)作圖群體根據(jù)遺傳穩(wěn)定性又可分為臨時(shí)性分離群體(F2、F3、BC1)和永久性分離群體(單倍體群體、重組自交系群體、永久性F2分離群體IF2等),個(gè)體間遺傳背景差異較大。次級(jí)作圖群體的群體基因組背景高度一致且與輪回親本相同,個(gè)體僅含少量供體親本片段,主要包括近等基因系類(lèi)群體(NILs)、單片段代換系(SSSL)群體和染色體片段置換系(CSSL)群體。

    1.2.2 ? 重組自交系(Recombinant inbred line,RIL)

    群體是雜種后代經(jīng)過(guò)多代自交而產(chǎn)生的一種作圖群體,即重組近交系群體,是由重組近交系組成的分離群體[8 ],通常從F2代開(kāi)始采用單粒傳的方法來(lái)建立。由于自交的作用是使基因型純合化,RIL群體中每個(gè)株系都是純合的,因而RIL群體是一種可以長(zhǎng)期使用的永久性分離群體。目前RIL群體被廣泛應(yīng)用,能保持單個(gè)純合基因型個(gè)體。

    1.2.3 ? ?近等基因系(Near isogenic line,NIL) ? ?群體將多個(gè)QTL位點(diǎn)分解成單個(gè)孟德?tīng)栠z傳因子,將數(shù)量性狀轉(zhuǎn)化為質(zhì)量形狀,從而可以對(duì)主效QTL進(jìn)行精細(xì)定位和圖位克隆[8 ]。

    1.2.4 ? ?深度雜交系(Advanced intercross lines,AIL)

    群體AIL開(kāi)始于F2群體,F(xiàn)1雜交后裔繼續(xù)雜交一定數(shù)目的世代,F(xiàn)2 ∶ 3繼續(xù)雜交后代(與RIL近似,但是遠(yuǎn)交,而不是近交)[8 ]。

    1.2.5 ? 雙單倍體(Doubled haploid,DH) ? ?群體指高等植物的單倍體經(jīng)過(guò)染色體加倍形成的二倍體。DH群體產(chǎn)生的途徑很多,亦因物種不同而異,最常見(jiàn)的方法是通過(guò)花藥培養(yǎng),即取F1 植株的花藥進(jìn)行離體培養(yǎng)誘導(dǎo)產(chǎn)生單倍體植株,然后對(duì)染色體進(jìn)行加倍產(chǎn)生DH植株。也可通過(guò)單倍體誘導(dǎo)系誘導(dǎo),通過(guò)自然加倍或化學(xué)藥劑處理加倍產(chǎn)生DH群體[8 ]。

    1.2.6 ? ?回交(Back Cross,BC)群體 ? 也是常用的作圖群體,如BC1為A與B的雜交F1代與A或者B回交產(chǎn)生的子代,直接反映了F1代配子的分離比例[8 ]。

    1.2.7 ? ?F2群體 ? ?為A與B的雜交F1代經(jīng)自交產(chǎn)生的子代,是一種暫時(shí)性分離群體。特點(diǎn)是群體內(nèi)個(gè)體間基因型不同,各個(gè)體的基因型雜合。這類(lèi)群體的優(yōu)點(diǎn)是不但可提供豐富的遺傳信息,而且可以用來(lái)估算加性效應(yīng)及顯性效應(yīng)。F2群體或衍生的F2 :3群體是最常用的作圖群體[8 ]。

    1.2.8 ? ?連鎖群體(Linkage group) ? 標(biāo)記組成的序列,染色體是堿基組成的序列,連鎖群上的標(biāo)記只是存在多態(tài)性的堿基位點(diǎn),即SNP,所以來(lái)自于同一條染色體的標(biāo)記一定構(gòu)建到同一連鎖群上[8 ]。

    2 ? 基于遺傳群體進(jìn)行基因定位方法

    目前,定位質(zhì)量性狀控制基因主要利用近等基因系分析法、連鎖分析法和混合分組分析法等途徑,在不同作物基因定位及遺傳改良中均有成功案例。

    2.1 ? 近等基因系分析法(Near isogenic line analysis)

    將遺傳背景相近或相同、多個(gè)農(nóng)藝性狀相似、個(gè)別染色體位置存在差異的材料稱(chēng)之為近等基因系。其主要獲得手段是將兩個(gè)具有不同目標(biāo)表型的親本P1和P2雜交,再與P1或P2多代回交后篩選得到在目標(biāo)表型上有差異的品系,這樣,品系以及品系間與輪回親本間就構(gòu)成了近等基因系。在育種中,親本輪回改良或品種特性改良應(yīng)用最多的就是近等基因系分析法[12 - 13 ],當(dāng)市場(chǎng)推廣品種缺少某個(gè)優(yōu)良性狀時(shí)(如不抗蟲(chóng)、易感?。?,常采用輪回轉(zhuǎn)育方法將外源品種中含有目標(biāo)性狀的抗性基因?qū)氲浆F(xiàn)有品種中,提高該品種的綜合抗性;若是雜交種,只改變存在缺陷性狀的父本或母本之一,以保證雜交種其他農(nóng)藝性狀的穩(wěn)定性。因此將用于多次回交的親本稱(chēng)為受體親本或輪回親本,是目標(biāo)表型的接受者,僅用于首次雜交時(shí)的親本稱(chēng)為供體親本或非輪回親本,是目標(biāo)表型的提供者。如此,多代回交的結(jié)果就是不斷提高后代中輪回親本的遺傳基因,減少供體材料的遺傳成分,逐漸向輪回親本遺傳背景純合,7~8代后理論上除了含有需改良的目標(biāo)性狀基因片段外,其他染色體遺傳信息與輪回親本幾乎相同(育種上講就是擁有了99%以上的輪回親本血緣)。為了加快回交后代基因組恢復(fù)成輪回親本的速度,在每代選擇繼續(xù)回交的植株時(shí),除了要保證含有供體目標(biāo)基因外,應(yīng)盡量選擇形態(tài)上與輪回親本接近的植株。因此,改良的品系與輪回親本間實(shí)際上構(gòu)成了一對(duì)近等基因系,通過(guò)持續(xù)的自交保持目標(biāo)性狀位點(diǎn)的雜合可以建立成對(duì)的近等基因系,利用近等基因系分析法可以克隆控制重要農(nóng)藝性狀的候選基因。

    2.2 ? 連鎖分析法(Linkage analysis)

    用連鎖分析法定位QTL的基本原理是經(jīng)典遺傳學(xué)中的連鎖分析(兩點(diǎn)測(cè)交、三點(diǎn)測(cè)交),利用與QTL連鎖的分子標(biāo)記進(jìn)行推斷[14 - 15 ]。因?yàn)楫?dāng)一個(gè)分子標(biāo)記與QTL連鎖時(shí),傾向于一起遺傳,這樣就可以借助分子標(biāo)記信息近似估計(jì)QTL在群體中的傳遞情況。還可利用區(qū)間定位方法。與單標(biāo)記分析法比較,區(qū)間定位法能大大提高QTL的檢測(cè)效率,并能較準(zhǔn)確的估計(jì)出QTL的位置和效應(yīng)值。區(qū)間定位已成為目前QTL檢測(cè)中的標(biāo)準(zhǔn)方法。QTL連鎖分析常用的統(tǒng)計(jì)分析方法有區(qū)間作圖法(IM)、復(fù)合區(qū)間作圖法(CIM)、混合線(xiàn)性模型(MLM)、完備區(qū)間作圖法(ICIM)等。關(guān)聯(lián)分析需要依賴(lài)作圖群體,構(gòu)建遺傳圖譜周期長(zhǎng),成本較高。

    2.3 ? 關(guān)聯(lián)分析法(Correlation analysis)

    關(guān)聯(lián)分析可以解決連鎖分析在挖掘種質(zhì)資源等位變異時(shí)效率低、成本高、在群體數(shù)量和研究目標(biāo)性狀方面依賴(lài)作圖群體等缺陷。GWAS分析是對(duì)多個(gè)個(gè)體全基因組范圍的遺傳變異(標(biāo)記)進(jìn)行檢測(cè),獲得基因型,進(jìn)而將基因型與可觀測(cè)的性狀即表型進(jìn)行群體水平的統(tǒng)計(jì)學(xué)分析,根據(jù)統(tǒng)計(jì)量或顯著性P值篩選出最有可能影響該性狀的遺傳標(biāo)記,定位與性狀相關(guān)的遺傳位點(diǎn)或挖掘與性狀變異相關(guān)的候選基因,在群體水平上解析性狀遺傳基礎(chǔ)[10, 16 ]。關(guān)聯(lián)分析具有許多優(yōu)點(diǎn),如可以直接使用自然群體或種質(zhì)資源,不需要專(zhuān)門(mén)構(gòu)建作圖群體;檢測(cè)效率高,可同時(shí)檢測(cè)同一座位的多個(gè)等位基因;節(jié)省研究時(shí)間且分辨率高。但相應(yīng)地,GWAS 也存在一定的缺點(diǎn),如當(dāng)群體結(jié)構(gòu)分化明顯時(shí)容易造成假陽(yáng)性,需要足夠多的樣本才能保證P值足夠低,來(lái)避免假陰性和假陽(yáng)性。同時(shí),精確的表型鑒定是GWAS成功定位的重要保障。

    2.4 ? 離體分組混合分析法(Bulked segregate analysis, BSA)

    離體分組混合分析法是1991年由R. W. MICHELMORE[10 ]在萵苣上首次應(yīng)用的一種快速定位控制目標(biāo)性狀基因的方法。取F2群體中具有極端表型的15~100個(gè)單株等量混合其DNA形成2個(gè)DNA池,然后在親本和兩個(gè)池之間進(jìn)行標(biāo)記多態(tài)篩選。如果某個(gè)標(biāo)記在親本和混池之間具有一致的多態(tài),則該標(biāo)記很可能和性狀連鎖,通過(guò)F2群體對(duì)這些篩選到的多態(tài)標(biāo)記進(jìn)行基因型分析即可完成對(duì)目標(biāo)基因的定位,而不需要對(duì)每個(gè)標(biāo)記在群體里進(jìn)行基因型分析。BSA的方法常用于單個(gè)QTL、突變位點(diǎn)定位等,BSA建立的混池因?yàn)榻?jīng)過(guò)了目標(biāo)性狀篩選,所以保證了2個(gè)混池之間除了目標(biāo)性狀外其他遺傳背景基本相同,僅目標(biāo)性狀關(guān)聯(lián)的DNA區(qū)段存在差異,所以?xún)蓚€(gè)混池又被稱(chēng)為近等基因池。該方法主要用于植物性狀的初定位,適用于單基因的質(zhì)量性狀和主效基因的數(shù)量性狀上,常用分離群體為F2、BC、RIL、DH群體等,常用的分子標(biāo)記有RAPD、RFLP、AFLP、SSR、SNP、InDel等。

    BSA技術(shù)應(yīng)用廣泛的主要有以下4種方法,分別是MutMap[17 - 18 ]、MutMap+[19 ]、MutMap-Gap[20 ]、QTL-seq[21 - 22 ],他們的基本原理相同,即針對(duì)研究的目標(biāo)性狀,選擇表型差異顯著的親本構(gòu)建出分離群體(或家系群體),再?gòu)姆蛛x群體中選取目標(biāo)性狀表型極端的一定數(shù)量的單體,混合構(gòu)建2個(gè)DNA池(DNA pools)。通常以雙親的DNA作為對(duì)照(參考基因組以利于對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行正確的分析和判斷),比較2個(gè)DNA池之間的差異,2個(gè)DNA池之間的差異片段即為候選區(qū)域,所關(guān)注的基因或者QTL可能存在于該候選區(qū)域中[23 ]。

    3 ? 不同MutMap方法及應(yīng)用

    3.1 ? MutMap方法

    3.1.1 ? ?MutMap定義及適宜條件 ? ?MutMap是基于高通量第二代測(cè)序技術(shù)的全基因組測(cè)序(Whole genome sequencing,WGS)發(fā)展起來(lái)的新的正向遺傳學(xué)基因定位和遺傳分析方法[17 ]。2012年Abe ? 等[17 ]開(kāi)發(fā)的MutMap方法應(yīng)用于水稻性狀篩選改良,利用甲基磺酸乙酯(ethylmethanesulfonate,EMS)誘變處理野生稻后獲得穩(wěn)定遺傳突變體,且目標(biāo)性狀由隱性單基因控制。取突變表型個(gè)體DNA等量混合獲得突變體DNA池(選擇30~200株不等)進(jìn)行全基因組重測(cè)序,混池測(cè)序數(shù)據(jù)比對(duì)誘變親本基因組獲取SNP位點(diǎn),計(jì)算SNP指數(shù),指數(shù)峰作為目的性狀關(guān)聯(lián)的目的基因所在的候選位置。利用該技術(shù)成功定位到4個(gè)水稻半矮稈、1個(gè)雄性不育、2個(gè)葉色淡綠性狀關(guān)聯(lián)基因在染色體上的遺傳位置[17 ]。

    3.1.2 ? MutMap原理及分析 ? 利用二代測(cè)序技術(shù),根據(jù)BSA極端性狀混池測(cè)序原理[17,23 ],分別抽提具有極端表型性狀植株的親本及子代植株DNA或RNA建庫(kù)測(cè)序。在分離群體搭建過(guò)程中,子代會(huì)根據(jù)表型進(jìn)行選擇,篩選出突變型子代池和野生型子代池。根據(jù)遺傳連鎖交換定律,子代池的基因型會(huì)和表型產(chǎn)生共分離,反應(yīng)在物理圖譜層面,與表型連鎖的染色體區(qū)段會(huì)和不連鎖的染色體區(qū)間產(chǎn)生穩(wěn)定的SNP-指數(shù)差異(SNP-index)。SNP指數(shù)的概念:SNP-index =突變SNP的reads數(shù)/對(duì)應(yīng)SNP位點(diǎn)的總reads數(shù)。差異大小范圍為0~1。如果SNP-index=0,說(shuō)明該池中僅包含來(lái)自參考基因序列的親本基因組片段;如果SNP-index=1,說(shuō)明該池中僅僅包含來(lái)自于另一個(gè)親本(突變體)的基因組片段;如果SNP-index=0.5,說(shuō)明該池中存在2個(gè)親本的基因組片段,此子代混池中SNP來(lái)自2個(gè)親本的基因組的頻率一致(圖1)。操作及分析步驟:①篩選親本間純合差異的SNP位點(diǎn);②子代SNP根據(jù)親本的基因型,確定突變親本來(lái)源的Reads,計(jì)算SNP-index;③只有突變子代池(Highbulk)的群體采用SNP-index定位,同時(shí)存在突變子代池和野生子代池(Lowbulk)的群體,采用ΔSNP-index定位,即將突變池的SNP指數(shù)減去野生池的SNP指數(shù)以獲得每個(gè)位點(diǎn)的SNP指數(shù)差值,即Δ(SNP-index);④定位時(shí)以每15個(gè)點(diǎn)為窗口,5個(gè)點(diǎn)步移統(tǒng)計(jì)該窗口內(nèi)所有SNP-index平均值,做出紅色擬合線(xiàn),排除部分噪音影響,使連鎖區(qū)域峰值更加明顯;⑤ΔSNP-index的群體可以利用二項(xiàng)分布檢驗(yàn)Highbulk和Lowbulk的覆蓋深度,考慮其混池?cái)?shù)量以及群體類(lèi)型,算出其統(tǒng)計(jì)學(xué)上95%的和99%的置信線(xiàn);⑥ΔSNP-index定位時(shí),擬合線(xiàn)超出置信線(xiàn)的染色體區(qū)域?yàn)榭赡艿谋硇瓦B鎖區(qū)域。根據(jù)擬合線(xiàn)超出染色體置信區(qū)間,判斷候選基因位置[17 ]。測(cè)序完成后,要對(duì)測(cè)序結(jié)果進(jìn)行系統(tǒng)分析,常用的且最重要的分析參數(shù)有以下6種:一是SNP,單核苷酸多態(tài)性(Single nucleotide polymorphism),主要是指在基因組水平上由單個(gè)核苷酸的變異所引起的DNA序列多態(tài)性;二是InDel,插入/缺失(Insertion/Deletion)基因組中小片段的插入和缺失序列;三是SV,基因組結(jié)構(gòu)變異(Structure variation),染色體結(jié)構(gòu)變異是指在染色體上發(fā)生了大片段的插入和缺失、重復(fù)復(fù)制、翻轉(zhuǎn)顛換、易位等;四是CNV,基因組拷貝數(shù)變異(Copy number variation),是基因組變異的一種形式,通常使基因組中大片段的DNA形成非正常的拷貝數(shù)量;五是基因組測(cè)序深度,測(cè)序得到的總堿基數(shù)與待測(cè)基因組大小的比值;六是基因組覆蓋率,由于大片段拼接的gap、測(cè)序讀長(zhǎng)有限、重復(fù)序列等問(wèn)題的存在,測(cè)序分析后組裝得到的基因組序列通常無(wú)法完全覆蓋所有區(qū)域,覆蓋度就是最終得到的結(jié)果占整個(gè)基因組的比例。

    3.1.3 ? MutMap應(yīng)用實(shí)例 ? 2011年日本地震導(dǎo)致海嘯爆發(fā)、海水倒灌農(nóng)田,使當(dāng)?shù)卦痉饰值耐恋刈兂闪他}堿地,水稻產(chǎn)量大幅下降??茖W(xué)家采用EMS誘變的方法對(duì)當(dāng)?shù)厮酒贩N“Hitomebore” 進(jìn)行了處理,F(xiàn)1代自交得到的 F2群體,獲得了耐鹽型單基因控制耐鹽突變株,野生型和突變型表型個(gè)體數(shù)的比值約為 3∶1,符合孟德?tīng)栯[性單基因控制遺傳定律。選取帶有目標(biāo)性狀的子代個(gè)體構(gòu)建混合池并測(cè)序,將測(cè)序結(jié)果與野生親本基因組比對(duì),得到SNP位點(diǎn),通過(guò)Δ(SNP_index)的方法將目標(biāo)性狀定位在6號(hào)染色體2.76~8.57 Mb之間。該候選區(qū)域其中1個(gè)SNP的突變使一個(gè)無(wú)義色氨酸密碼子突變?yōu)榻K止子,在篩選到與耐鹽性狀相關(guān)的基因OsRR22,導(dǎo)致突變體hst1產(chǎn)生耐鹽表型。該研究縮短了耐鹽品種水稻的育種進(jìn)程,將對(duì)當(dāng)?shù)厮井a(chǎn)量恢復(fù)做出巨大貢獻(xiàn)[18 ]。

    3.2 ? MutMap+方法

    3.2.1 ? MutMap+適宜條件及原理 ? 常規(guī)的MutMap方法是將 M3~M5世代的突變體與野生型親本株系進(jìn)行雜交構(gòu)建后代群體,在此基礎(chǔ)上進(jìn)行BSA分析。對(duì)于人工雜交困難的突變體,以及發(fā)育早期死亡或不育的突變體,因?yàn)闊o(wú)法直接構(gòu)建子代群體,MutMap方法就不適用了。為了解決這個(gè)問(wèn)題,F(xiàn)ekih等[19 ]開(kāi)發(fā)了MutMap+方法:隱性純合M1代株系不育或死亡,因此利用雜合的M1代自交產(chǎn)生的雜合后代來(lái)保存變異基因,并在構(gòu)建群體時(shí)重點(diǎn)關(guān)注那些自交后代表型發(fā)生分離的株系。在導(dǎo)致表型的突變位點(diǎn)上野生型株系理論上有2/3是雜合的。以水稻為例,簡(jiǎn)單闡述MutMap+方法的原理(圖2)。

    3.2.2 ? MutMap+應(yīng)用實(shí)例 ? ?應(yīng)用該方法不雜交只自交,在雜合的株系自交獲得的后代群體中,野生型表型和隱性純合突變表型分離比例約為3∶1。將野生型 M1后代個(gè)體自交,收集雜合的M1株系自交獲得的M2代群體。同樣收集雜合的M2代自交獲得的M3群體。根據(jù)表型將M3株系分為兩組,構(gòu)建2組混合DNA池,并進(jìn)行全基因組重測(cè)序?;斐販y(cè)序數(shù)據(jù)比對(duì)野生型親本基因組獲取SNP位點(diǎn),2個(gè)混池分別計(jì)算SNP-index,Δ(SNP-index)峰作為目的性狀關(guān)聯(lián)的目的基因所在的候選位置。該方法主要用來(lái)定位早衰、致死、不育、白化等性狀的定位分析,也用來(lái)分析花蕊太小或去雄難度大而尚未建立有效的人工雜交物種。例如谷子:利用MutMap+方法,F(xiàn)ekin等找到了表型為早衰和過(guò)早死亡的突變體Hit9188的候選基因OsNAP6,并且找到了表型為苗白化和過(guò)早死亡的突變體Hit11440 的候選基因Os08g0139100[19 ]。

    3.3 ? Mutmap-Gap方法

    3.3.1 ? Mutmap-Gap適宜條件 ? Mutmap或Mutmap+都需要根據(jù)參考基因組序列先構(gòu)建一個(gè)誘變親本的Pseudo-genome。如果野生型親本的基因組和參考基因組的序列有結(jié)構(gòu)性差異,而導(dǎo)致性狀改變的區(qū)域位點(diǎn)恰巧不在參考基因組序列上,相對(duì)參考基因組是一個(gè)大的插入片段,是參考基因組上的缺乏的區(qū)域(Gap),直接做Mutmap分析將得不到關(guān)聯(lián)位點(diǎn)的全部信息。為了鑒別在參考基因組序列Gap內(nèi)的變異,Takagi等[20 ]開(kāi)發(fā)了Map-Gap方法。該方法將MutMap方法與對(duì)基因組gap區(qū)域的從頭(de novo)組裝相結(jié)合,對(duì)這類(lèi)特異的突變基因進(jìn)行定位和分離。MutMap-Gap將MutMap方法得出的靶基因組區(qū)域的gap通過(guò)de novo組裝補(bǔ)齊,以判斷不在參考基因組上變異候選區(qū)段位置(圖3)。

    3.3.2 ? Mutmap-Gap原理及應(yīng)用 ? 首先將比對(duì)到候選區(qū)域的臨近區(qū)域和未對(duì)應(yīng)參考序列的reads序列合并起來(lái)做de novo組裝,將組裝得到scaffold序列和pseudo-genome序列一起當(dāng)作參考序列(圖3),重新比對(duì)計(jì)算SNP指數(shù),關(guān)聯(lián)定位區(qū)域。如突變位于P特異性基因組區(qū)域內(nèi),僅通過(guò)MutMap分析無(wú)法識(shí)別。因?yàn)槎ㄎ粎^(qū)域相對(duì)參考基因組是一個(gè)新的發(fā)現(xiàn),參考序列里面沒(méi)有基因相關(guān)的信息,通過(guò)對(duì)P的全基因組測(cè)序及與參考基因組的重新組裝pseudo-genome序列,需要對(duì)該區(qū)域進(jìn)行基因預(yù)測(cè)和注釋?zhuān)缘玫侥康幕?。一般情況下Mutmap或Mutmap+沒(méi)定位到結(jié)果,可以嘗試用Mutmap-Gap。

    3.4 ? QTL-Seq方法

    3.4.1 ? QTL-Seq方法原理及操作 ? MutMap系列僅適用于質(zhì)量性狀,對(duì)于數(shù)量性狀,產(chǎn)生了QTL-seq法。QTL-seq是一種在極端表型混池高通量測(cè)序的基礎(chǔ)上,根據(jù)極端表型池間SNP頻率的差異進(jìn)行QTL/基因定位的方法,是一種結(jié)合BSA與NGS 來(lái)快速定位單個(gè)數(shù)量性狀的方法[21 - 22 ]。具體操作:選取目標(biāo)性狀差異大的兩個(gè)親本雜交,產(chǎn)生分離群體如F2、RILs、DH 等,目標(biāo)性狀理論上將會(huì)呈現(xiàn)正態(tài)分布選取目標(biāo)性狀表型極端的20%的比例的個(gè)體分別混合成2個(gè)表型極端池,進(jìn)行重測(cè)序。以對(duì)照親本的基因型為參照,計(jì)算子代極端池中的SNP-index[17 ]。若該參數(shù)為0,代表子代所有測(cè)到的Reads都來(lái)自野生親本;若該參數(shù)為1,代表子代所有Reads都來(lái)自突變親本;該參數(shù)為0.5,說(shuō)明該池中存在2個(gè)親本的基因組片段。一般而言,大部分位點(diǎn)的SNP index值在0.5左右的區(qū)域,但如果某SNP與目標(biāo)QTL連鎖,導(dǎo)致SNP-index的值偏離0.5。比較ΔSNP-index,候選區(qū)間內(nèi)的ΔSNP-index應(yīng)該接近于1,從而將QTL定位到了基因組上的某一區(qū)域[18 - 19 ](圖4)。

    3.4.2 ? QTL-Seq方法應(yīng)用實(shí)例 ? QTL-seq定位黃瓜早花農(nóng)藝性狀:采用BSA混樣策略對(duì)10株F2子代群體極端性狀(早花和晚花)的樣品混合的DNA池,及其親本進(jìn)行基因組重測(cè)序。通過(guò)全基因組掃描SNP,分析頻率差異,檢測(cè)F2群體早花性狀的QTL,找到了一個(gè)位于早花QTL Ef1.1中的候選基因。利用F2群體,SSR標(biāo)記構(gòu)建的遺傳圖譜進(jìn)行QTL定位,兩種策略結(jié)合,將Ef1.1縮小到980 Kb,其中包含基因Csa1G651710,與FT基因同源[22 ]。

    3.5 ? 4種不同MutMap方法的比較

    從表1可以看出,MutMap、MutMap+、Mut Map-Gap、QTL-seq這4種方法在樣本材料、親本類(lèi)型、雜交方式、所需測(cè)序樣本、適用范圍等不盡相同,但均可快速獲得與目標(biāo)基因連鎖的分子標(biāo)記,且遺傳定位的周期短、效率高。基于誘變突變體的MutMap方法、自交的MutMap+方法、用于識(shí)別基因組缺失區(qū)間變異的MutMap-Gap方法以及用于定位數(shù)量性狀基因座的QTL-seq方法等均不需要建立繁瑣的后代定位群體,加快了對(duì)重要表型變異位點(diǎn)的識(shí)別過(guò)程,縮短了育種時(shí)間,大大加快了新品種的開(kāi)發(fā)速度。

    4 ? 小結(jié)與展望

    由于MutMap方法的應(yīng)用使得基因得以快速定位,這些基因資源將有可能促進(jìn)玉米、水稻等作物的遺傳改良[1 - 2 ]?;赪GS的MutMap方法的正向遺傳學(xué)研究,不僅有利于QTL/基因的定位,也利于新的基因資源在作物育種中的應(yīng)用。早在1991年,Michelmore等[10 ]就建立了BSA法,克服了許多作物沒(méi)有或難以創(chuàng)建相應(yīng)的NIL群體的限制,在自交和異交作物中均有廣泛的應(yīng)用前景。對(duì)于尚無(wú)連鎖圖或連鎖圖飽和程度較低的植物,BSA法也是快速獲得與目標(biāo)基因連鎖的分子標(biāo)記的有效方法,該方法的優(yōu)點(diǎn)是遺傳定位的周期短且效率高。在此基礎(chǔ)上擴(kuò)展的新方法也不斷出現(xiàn),如基于自交的MutMap+、用于識(shí)別基因組缺失區(qū)間變異的MutMap-Gap,以及用于定位數(shù)量性狀基因座的QTL-seq方法等。這些方法的優(yōu)點(diǎn)是不需要建立繁瑣的后代定位群體,甚至不依賴(lài)于遺傳雜交和任何連鎖信息,加快了對(duì)重要表型變異位點(diǎn)的識(shí)別過(guò)程。Mutmap采用EMS誘變技術(shù)將育種時(shí)間縮短至傳統(tǒng)育種的1/5左右,大大加快了新品種的開(kāi)發(fā),通過(guò)利用MutMap對(duì)突變體進(jìn)行遺傳定位研究,可以快速鑒別出靶基因[23 - 24 ]。

    過(guò)去,人們依靠傳統(tǒng)育種方法獲取作物新品種,根據(jù)選擇的特性確定植物親本,然后通過(guò)雜交、回交或者直系篩選程序來(lái)完成。傳統(tǒng)育種方法存在著極大的盲目性、經(jīng)驗(yàn)性、不確定因素和長(zhǎng)周期等問(wèn)題。為了實(shí)現(xiàn)育種4.0的目標(biāo),需要為分子育種制定一個(gè)育種路線(xiàn)圖。應(yīng)用新的生物技術(shù),依托人工智能、基因組測(cè)序、基因編輯等相關(guān)技術(shù),實(shí)現(xiàn)玉米組學(xué)基因型與表型大數(shù)據(jù)的快速積累。通過(guò)遺傳變異等數(shù)據(jù)的整合,實(shí)現(xiàn)作物性狀調(diào)控基因的快速挖掘與表型的精準(zhǔn)預(yù)測(cè)。通過(guò)人工改造基因元器件與人工合成基因回路,使作物具備新的抗逆、高效等生物學(xué)性狀,創(chuàng)建智能組合優(yōu)良等位基因的自然變異、人工變異、數(shù)量性狀位點(diǎn)的育種設(shè)計(jì)方案,最終實(shí)現(xiàn)智能、高效、定向培育新品種[25 ]。分子育種的精準(zhǔn)實(shí)施可幫助育種家在最短的時(shí)間里實(shí)現(xiàn)以下目標(biāo)。①大大縮短育種時(shí)間,加速育種進(jìn)程。比如玉米自交系選育中,通過(guò)常規(guī)育種回交6代及以上才能達(dá)到99%純合的后代,而分子標(biāo)記選育只需要3代就可以達(dá)到,利用單倍體誘導(dǎo)技術(shù)只需要2代就能達(dá)到100%純合。②最大限度地繞過(guò)物種生殖隔離的障礙,實(shí)現(xiàn)生物界遺傳物質(zhì)的自由交流。利用基因重組操作技術(shù),將不同物種中特定的基因與載體結(jié)合,可以導(dǎo)入到其他物種中表達(dá),發(fā)揮其基因功能,利用轉(zhuǎn)基因技術(shù)讓現(xiàn)存的作物品種快速適應(yīng)新環(huán)境[26 ]。如抗棉鈴蟲(chóng)轉(zhuǎn)基因棉花的普及和大面積種植推廣,以及抗旱耐鹽堿、抗除草劑等玉米品種的獲得,都大大降低了人工打藥、鋤草、灌溉等成本,提高了產(chǎn)量和收入,且轉(zhuǎn)基因觀賞園藝植物在休閑觀光農(nóng)業(yè)中被廣泛應(yīng)用。③可幫助育種從多個(gè)目標(biāo)性狀帶有經(jīng)驗(yàn)性的多代重組選擇,轉(zhuǎn)變成針對(duì)目標(biāo)性狀進(jìn)行基因編輯的精準(zhǔn)性狀選擇和改良,可將主效基因及其優(yōu)異等位基因集于一個(gè)骨干自交系或品種中,大大提高了育種效率。如孟山都公司利用分子標(biāo)記將多個(gè)不同的抗旱位點(diǎn)聚合育成高抗旱玉米品種(gene stacking),是復(fù)雜性狀通過(guò)分子標(biāo)記獲得育成品種的成功案例[3 ]。劉忠祥[27 ]系統(tǒng)闡述玉米株高主效QTL定位研究進(jìn)展及與株高相關(guān)基因的功能與響應(yīng)途徑, 對(duì)玉米分子育種的選育工作有積極借鑒作用。④隨著人們對(duì)身體健康,美容保健等方面的重視,育種家通過(guò)分子育種可選育功能保健型(低脂低糖、高蛋白、高氨基酸等)品種,以滿(mǎn)足人們通過(guò)飲食對(duì)特殊營(yíng)養(yǎng)成分需求,如“黃金大米”,其β-胡蘿卜素的含量是普通大米的23倍。⑤基因組編輯技術(shù)作為一項(xiàng)變革性的新興技術(shù),尤其是CRISPR/Cas 9基因編輯技術(shù)的興起和廣泛應(yīng)用,將推動(dòng)新技術(shù)產(chǎn)品產(chǎn)業(yè)化的發(fā)展[28 - 29 ]。中國(guó)科學(xué)院高彩霞團(tuán)隊(duì)等聯(lián)合攻關(guān),闡明了小麥新型mlo突變體既抗白粉病又高產(chǎn)的分子機(jī)制,并通過(guò)CRISPR/Cas9多重基因組編輯,使主栽小麥品種快速獲得廣譜抗白粉病的優(yōu)異性狀[30 ]。然而,沒(méi)有前期的對(duì)關(guān)鍵基因進(jìn)行克隆和功能的深入剖析,根本談不上DNA水平的分子標(biāo)記輔助選擇育種。因此,加強(qiáng)農(nóng)業(yè)種質(zhì)資源保護(hù)和開(kāi)發(fā)利用,創(chuàng)新農(nóng)作物分子改良與技術(shù)應(yīng)用體系,加強(qiáng)作物分子育種方面的應(yīng)用基礎(chǔ)性研究是未來(lái)提高中國(guó)種業(yè)發(fā)展水平最為重要的步驟和方向[31 - 35 ]。

    參考文獻(xiàn):

    [1] 種 ? 康,李家洋. ?科學(xué)發(fā)展催生新一輪育種技術(shù)革命[J]. ?中國(guó)科學(xué):生命科學(xué),2021,51(10):1353-1355.

    [2] 景海春,田志喜,種 ? 康,等. ?分子設(shè)計(jì)育種的科技問(wèn)題及其展望概論[J]. ?中國(guó)科學(xué):生命科學(xué),2021,

    51(10):1356-1365.

    [3] 薛勇彪,種 ? 康,韓 ? 斌,等. ?創(chuàng)新分子育種科技支撐我國(guó)種業(yè)發(fā)展[J]. ?中國(guó)科學(xué)院院刊,2018,33(9):893-899.

    [4] PINGALI P L. ?Green revolution:impacts,limits,and the path ahead[J]. ?Proceedings of the National Academy of Sciences of the United States of America,USA,2012,109:12302-12308.

    [5] 周文期,連曉榮,劉忠祥,等. ?玉米株高和穗位高的調(diào)控機(jī)理研究[J]. ?分子植物育種,2021,19(23):7965-7976.

    [6] 周文期,連曉榮,周玉乾,等. ?EMS誘變玉米自交系種質(zhì)創(chuàng)新應(yīng)用[J]. ?玉米科學(xué),2020,28(6):31-38.

    [7] 周文期,王曉娟,寇思榮,等. ?玉米葉夾角形成的分子調(diào)控機(jī)理研究[J]. ?土壤與作物,2019,8(3):339-348.

    [8] 朱玉賢,李 ? 毅,鄭曉峰,等. ?現(xiàn)代分子生物學(xué)[M]. ?5版. ?北京:高等教育出版社,2019.

    [9] COULSON A,SULSTON J,BRENNER S,KARN J. ?Toward a physical map of the genome of the nematode Caenorhabditis elegans[J]. ?Proceedings of the National Academy of Sciences of the United States of America,USA,1986,83:7821-7825.

    [10] 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[J]. ?Proceedings of the National Academy of Sciences, 1991,88(21):9828-9832.

    [11] VISSCHER P M,WRAY N R,ZHANG Q,et al. ?10 years of GWAS discovery:Biology,function,and translation[J]. ?American Journal of Human Genetics, 2017,101(1):5-22.

    [12] 張 ? 葉,王梓鈺,邢躍先,等. ?利用近等基因系定位玉米無(wú)葉舌基因的研究[J]. ?玉米科學(xué),2021,29(3):29-35

    [13] 傅大雄,阮仁武,劉大軍,等. ?近等基因系法對(duì)小麥顯性矮源的研究[J]. ?中國(guó)農(nóng)業(yè)科學(xué),2007,40(4):655-664

    [14] 王 ? 蕓. ?關(guān)聯(lián)分析和連鎖分析定位控制水稻有效穗數(shù)QTL [D]. ?北京:中國(guó)農(nóng)業(yè)科學(xué)院,2018.

    [15] 田 ? 玉,馬春紅,宋麗華,等. ?數(shù)量遺傳性狀基因定位方法研究進(jìn)展[J]. ?河北農(nóng)業(yè)科學(xué),2021,25(5):88-91.

    [16] 涂雨辰,田 ? 云,盧向陽(yáng). ?全基因組關(guān)聯(lián)分析在植物中的應(yīng)用[J]. ?化學(xué)與生物工程,2013,30(6):1672-5425.

    [17] ABE A,KOSUGI S,YOSHIDA K,et al. ?Genome sequencing reveals agronomically important loci in rice using MutMap[J]. ?Nature Biotechnology,2012,30(2):174-178.

    [18]TAKAGI,H.,TAMIRU,M.,ABE,A. et al. ?MutMap accelerates breeding of a salt-tolerant rice cultivar[J]. ? ?Nature Biotechnology,2015,33,445-449.

    [19] FEKIH R,TAKAGI H,TAMIRU M,et al. ?MutMap+:Genetic mapping and mutant identification without crossing in rice[J]. ?PLoS One,2013,8(7):e68529.

    [20] TAKAGI H,UEMURA A,YAEGASHI H,et al.,MutMap-Gap:whole-genome resequencing of mutant F2 progeny bulk combined with de novo assembly of gap regions identifies the rice blast resistance gene Pii[J]. ?New Phytologist,2013,200(1):276-283.

    [21] TAKAGI H,ABE A,YOSHIDA K,et al. ?QTL-seq:rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations[J]. ?Plant Journal,2013,74(1):174-183.

    [22] LU H,LIN T,KLEIN J,et al. ?QTL-seq identifies an early flowering QTL located near Flowering Locus T. in cucumber[J]. ?Theoretical and Applied Genetics,2014,

    127:1491-1499

    [23] 袁金紅,李俊華,袁嬌嬌,等. 基于全基因組測(cè)序的MutMap方法在正向遺傳學(xué)研究中的應(yīng)用[J]. ?遺傳, 2017,39(12):1168-1177

    [24] 陸才瑞,鄒長(zhǎng)松,宋國(guó)立. ?高通量測(cè)序技術(shù)結(jié)合正向遺傳學(xué)手段在基因定位研究中的應(yīng)用[J]. ?遺傳,2015,37(8):765-776.

    [25] 王向峰,才 ? 卓. ?中國(guó)種業(yè)科技創(chuàng)新的智能時(shí)代—“玉米育種4.0”[J]. ?玉米科學(xué) 2019,27(1):1-9.

    [26] 王紅梅,陳玉梁,石有太,等. ?中國(guó)作物分子育種現(xiàn)狀與展望[J]. ?分子植物育種,2020,18(2):507-513.

    [27] 劉忠祥. ?玉米株高主效QTL定位研究綜述[J]. ?甘肅農(nóng)業(yè)科技,2018(9):62-69.

    [28] 盧俊南,褚 ? 鑫,潘燕平,等. ?基因編輯技術(shù):進(jìn)展與挑戰(zhàn)[J]. ?中國(guó)科學(xué)院院刊,2018,33(11):1184-1192.

    [29] 陳云偉,陶 ? 誠(chéng),周海晨,等. ?基因編輯技術(shù)研究進(jìn)展與挑戰(zhàn)[J]. ?世界科技研究與發(fā)展,2021,43(1):8-23.

    [30] LI S N,LIU D X,ZHANG Y W,et al. ?Genome-edited powdery mildew resistance in wheat without growth penalties[J]. ?Nature,2022,602:455-460.

    [31] 鄭懷國(guó),趙靜娟,秦曉婧,等. ?全球作物種業(yè)發(fā)展概況及對(duì)我國(guó)種業(yè)發(fā)展的戰(zhàn)略思考[J]. ?中國(guó)工程科學(xué),2021,23(4):45-55.

    [32] 周喜旺,劉鴻燕,王 ? 娜,等. ?DNA分子標(biāo)記技術(shù)在小麥遺傳育種中的應(yīng)用綜述[J]. ?甘肅農(nóng)業(yè)科技,2017(5):64-68.

    [33] 李聞娟,齊燕妮,王利民,等. ?運(yùn)用SRAP分子標(biāo)記對(duì)胡麻雜交種純度的鑒定研究[J]. ?甘肅農(nóng)業(yè)科技,2019(9):59-62.

    [34] 陳 ? 琛,羅俊杰,陳衛(wèi)國(guó). ?利用SSR技術(shù)快速鑒定2個(gè)辣椒雜交品種純度[J]. ?甘肅農(nóng)業(yè)科技,2020(4):53-58.

    [35] 劉新星,歐巧明,羅俊杰,等. ?利用SSR標(biāo)記鑒定當(dāng)歸的真實(shí)性[J]. ?甘肅農(nóng)業(yè)科技,2020(7):43-48.

    收稿日期:2022 - 03 - 09

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(32160490、31860384);2020年甘肅省科協(xié)青年科技人才托舉工程項(xiàng)目;甘肅省農(nóng)業(yè)科學(xué)院創(chuàng)新專(zhuān)項(xiàng)-博士基金(2020GAAS34)。

    作者簡(jiǎn)介:周文期(1985 — ),男,甘肅靜寧人,副研究員,博士,主要從事玉米遺傳育種及基因功能研究工作。Email: zhouwenqi850202@163.com。

    通信作者:周玉乾(1979 — ),男,甘肅靖遠(yuǎn)人,研究員,主要從事玉米育種工作。Email:yuqianzhou2008@163.com。

    猜你喜歡
    基因克隆農(nóng)藝性狀
    種植密度對(duì)伊犁河谷移栽橡膠草農(nóng)藝性狀及產(chǎn)量的影響
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    三個(gè)小麥防御素基因的克隆及序列分析
    新型緩釋肥在全膜雙壟溝播玉米上的應(yīng)用效果研究
    武運(yùn)粳24號(hào)水稻機(jī)插精確定量高產(chǎn)栽培技術(shù)研究
    不同夾心肥料對(duì)玉米生長(zhǎng)的影響
    從農(nóng)藝性狀及化學(xué)成分測(cè)定分析遠(yuǎn)志藥材商品品規(guī)和良種選育的合理性
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達(dá)分析
    山桃醇腈酶基因PdHNL1的克隆、序列分析與原核表達(dá)
    鵝PRL基因克隆及在繁殖相關(guān)組織中的表達(dá)規(guī)律
    露出奶头的视频| 特大巨黑吊av在线直播| 亚洲成av人片在线播放无| 精品一区二区三区视频在线观看免费| 男女做爰动态图高潮gif福利片| 香蕉丝袜av| 亚洲国产色片| 一级毛片高清免费大全| 嫩草影视91久久| 五月玫瑰六月丁香| 精品欧美国产一区二区三| 法律面前人人平等表现在哪些方面| 一本综合久久免费| 久久人人精品亚洲av| 国产一区二区在线av高清观看| 男人的好看免费观看在线视频| 日本三级黄在线观看| www日本在线高清视频| 午夜免费成人在线视频| 免费看美女性在线毛片视频| 99久国产av精品| 国产亚洲精品av在线| 久久精品人妻少妇| 久久久精品欧美日韩精品| 久99久视频精品免费| 最新美女视频免费是黄的| 精品福利观看| 国产精品,欧美在线| 欧美成人一区二区免费高清观看 | 我的老师免费观看完整版| 一进一出好大好爽视频| 两性夫妻黄色片| 中文字幕最新亚洲高清| 国产精品亚洲一级av第二区| 国产精品一及| 午夜福利视频1000在线观看| 99热精品在线国产| 国产精品久久久久久亚洲av鲁大| 久久中文字幕一级| 91麻豆av在线| 成人高潮视频无遮挡免费网站| 人人妻人人澡欧美一区二区| 亚洲欧美日韩卡通动漫| 天堂影院成人在线观看| 国产精品亚洲一级av第二区| 久久中文看片网| 69av精品久久久久久| 99久久综合精品五月天人人| 国产视频一区二区在线看| 男女做爰动态图高潮gif福利片| tocl精华| 一级毛片高清免费大全| 亚洲国产精品999在线| 国产伦人伦偷精品视频| 国产精品 欧美亚洲| 成人高潮视频无遮挡免费网站| 美女免费视频网站| 夜夜爽天天搞| 免费观看精品视频网站| 我的老师免费观看完整版| 九色成人免费人妻av| 国产淫片久久久久久久久 | 亚洲欧美日韩高清专用| 天堂网av新在线| 成年女人看的毛片在线观看| 色综合站精品国产| 欧美中文日本在线观看视频| av天堂中文字幕网| 久久久水蜜桃国产精品网| 麻豆国产av国片精品| 国内精品久久久久精免费| 国产精品99久久久久久久久| 国产精品av视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 校园春色视频在线观看| 国产精品久久久久久久电影 | 欧美成狂野欧美在线观看| 亚洲在线观看片| 动漫黄色视频在线观看| 一个人免费在线观看的高清视频| 成熟少妇高潮喷水视频| 成年女人看的毛片在线观看| 亚洲狠狠婷婷综合久久图片| 欧美中文综合在线视频| 国产午夜福利久久久久久| 五月玫瑰六月丁香| 国产成人一区二区三区免费视频网站| 少妇熟女aⅴ在线视频| 国产一区二区在线观看日韩 | 成人性生交大片免费视频hd| 国产午夜精品论理片| 亚洲乱码一区二区免费版| 我要搜黄色片| 国产 一区 欧美 日韩| 午夜两性在线视频| 中文字幕熟女人妻在线| 丰满人妻熟妇乱又伦精品不卡| 中文亚洲av片在线观看爽| 五月玫瑰六月丁香| 欧美3d第一页| 国产毛片a区久久久久| 色噜噜av男人的天堂激情| 男女视频在线观看网站免费| 级片在线观看| 午夜福利在线在线| 日本 欧美在线| 国产 一区 欧美 日韩| 两个人的视频大全免费| 久久中文字幕人妻熟女| 欧美大码av| 全区人妻精品视频| 亚洲成人中文字幕在线播放| 综合色av麻豆| 免费观看的影片在线观看| svipshipincom国产片| 一区二区三区国产精品乱码| 成人特级黄色片久久久久久久| 嫩草影视91久久| 久久亚洲精品不卡| 少妇裸体淫交视频免费看高清| 久久久久国产一级毛片高清牌| 精品99又大又爽又粗少妇毛片 | 亚洲 国产 在线| 国产一区在线观看成人免费| 欧美一区二区国产精品久久精品| 欧美在线一区亚洲| e午夜精品久久久久久久| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av| 男女午夜视频在线观看| 亚洲精品粉嫩美女一区| 欧美成狂野欧美在线观看| 午夜久久久久精精品| 白带黄色成豆腐渣| 老汉色av国产亚洲站长工具| 99久国产av精品| 又紧又爽又黄一区二区| 国产一区二区在线av高清观看| 夜夜夜夜夜久久久久| 天天一区二区日本电影三级| 中文字幕人妻丝袜一区二区| 亚洲第一欧美日韩一区二区三区| 三级国产精品欧美在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 日本 av在线| 一本一本综合久久| 日韩国内少妇激情av| 欧美午夜高清在线| xxxwww97欧美| 动漫黄色视频在线观看| 亚洲七黄色美女视频| 国产人伦9x9x在线观看| 日本与韩国留学比较| 啪啪无遮挡十八禁网站| 中文字幕高清在线视频| 巨乳人妻的诱惑在线观看| 久久久久国内视频| 中文资源天堂在线| 亚洲av美国av| 亚洲第一欧美日韩一区二区三区| 国产视频一区二区在线看| 国产一区二区激情短视频| 性欧美人与动物交配| 在线观看66精品国产| 国产av不卡久久| 欧美黄色淫秽网站| 身体一侧抽搐| cao死你这个sao货| 国产成人精品久久二区二区免费| 亚洲激情在线av| 欧美日韩一级在线毛片| 老司机在亚洲福利影院| 嫩草影院精品99| 国内久久婷婷六月综合欲色啪| 12—13女人毛片做爰片一| 偷拍熟女少妇极品色| 午夜福利成人在线免费观看| 久久精品人妻少妇| 欧美xxxx黑人xx丫x性爽| 免费一级毛片在线播放高清视频| 精品国产美女av久久久久小说| 欧美日韩福利视频一区二区| 在线a可以看的网站| 床上黄色一级片| 老司机午夜福利在线观看视频| 午夜两性在线视频| 最近最新中文字幕大全免费视频| 可以在线观看毛片的网站| 美女免费视频网站| av女优亚洲男人天堂 | 男插女下体视频免费在线播放| 色综合婷婷激情| 国产精品av视频在线免费观看| 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女| 午夜亚洲福利在线播放| 亚洲 欧美一区二区三区| 欧美av亚洲av综合av国产av| 成人av在线播放网站| www日本在线高清视频| 亚洲18禁久久av| 久久精品国产99精品国产亚洲性色| 小说图片视频综合网站| 我要搜黄色片| 国产激情偷乱视频一区二区| 草草在线视频免费看| 一级作爱视频免费观看| 欧美日韩综合久久久久久 | 国产男靠女视频免费网站| 欧美黑人巨大hd| 99热精品在线国产| 久久亚洲精品不卡| 性欧美人与动物交配| 色哟哟哟哟哟哟| 一夜夜www| 亚洲在线自拍视频| 久久精品国产亚洲av香蕉五月| 精品熟女少妇八av免费久了| 日韩中文字幕欧美一区二区| 欧美zozozo另类| 国产伦精品一区二区三区四那| 色噜噜av男人的天堂激情| 99精品久久久久人妻精品| 欧美日韩瑟瑟在线播放| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 天堂动漫精品| 级片在线观看| 一区福利在线观看| 首页视频小说图片口味搜索| 国产爱豆传媒在线观看| 欧美成人一区二区免费高清观看 | 中国美女看黄片| 亚洲国产日韩欧美精品在线观看 | 久久九九热精品免费| 国产精品亚洲美女久久久| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 三级国产精品欧美在线观看 | 悠悠久久av| 色播亚洲综合网| 国产激情偷乱视频一区二区| 亚洲国产欧洲综合997久久,| 中文字幕最新亚洲高清| 精品国产三级普通话版| www国产在线视频色| 99久久久亚洲精品蜜臀av| 免费搜索国产男女视频| 日韩国内少妇激情av| 国产成人福利小说| 亚洲av成人精品一区久久| 美女午夜性视频免费| 黄色片一级片一级黄色片| 国产视频一区二区在线看| 国产黄色小视频在线观看| 欧美最黄视频在线播放免费| 国产91精品成人一区二区三区| 首页视频小说图片口味搜索| 人人妻,人人澡人人爽秒播| 久久性视频一级片| 欧美不卡视频在线免费观看| 观看美女的网站| 丁香六月欧美| 亚洲av美国av| 制服人妻中文乱码| 婷婷六月久久综合丁香| 中出人妻视频一区二区| 免费看a级黄色片| 国产真实乱freesex| 黄色成人免费大全| 免费人成视频x8x8入口观看| 又黄又粗又硬又大视频| 国产精品av久久久久免费| 久久久精品大字幕| 1024香蕉在线观看| 99热这里只有精品一区 | 午夜精品久久久久久毛片777| 国语自产精品视频在线第100页| 国产一区二区三区在线臀色熟女| 99国产精品99久久久久| 床上黄色一级片| 亚洲天堂国产精品一区在线| 久久久成人免费电影| www国产在线视频色| 日韩精品青青久久久久久| 欧美乱妇无乱码| 国产乱人伦免费视频| 国产精品av久久久久免费| 欧美最黄视频在线播放免费| 长腿黑丝高跟| 久久久精品大字幕| 深夜精品福利| 日本一本二区三区精品| 丰满人妻一区二区三区视频av | 免费在线观看日本一区| 国产高清视频在线观看网站| 国产精品一区二区免费欧美| 在线免费观看的www视频| 亚洲国产中文字幕在线视频| 丁香欧美五月| 亚洲成a人片在线一区二区| 麻豆久久精品国产亚洲av| 18美女黄网站色大片免费观看| 老司机深夜福利视频在线观看| 久久久水蜜桃国产精品网| 久久国产精品人妻蜜桃| 色吧在线观看| 国产成人啪精品午夜网站| 2021天堂中文幕一二区在线观| 熟女人妻精品中文字幕| ponron亚洲| 免费观看的影片在线观看| 久久国产精品影院| 亚洲精品久久国产高清桃花| 在线观看免费视频日本深夜| 丝袜人妻中文字幕| 一二三四在线观看免费中文在| 狠狠狠狠99中文字幕| 国产熟女xx| 国产亚洲欧美在线一区二区| 99热精品在线国产| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 国产欧美日韩精品亚洲av| 国产aⅴ精品一区二区三区波| 国产又色又爽无遮挡免费看| 亚洲av电影不卡..在线观看| 国产精品久久久av美女十八| 91麻豆av在线| 国产一区二区在线观看日韩 | 国产在线精品亚洲第一网站| 亚洲18禁久久av| 久久国产精品人妻蜜桃| 国产精品,欧美在线| 在线播放国产精品三级| 一进一出好大好爽视频| 丰满人妻熟妇乱又伦精品不卡| 一a级毛片在线观看| 亚洲精品456在线播放app | 法律面前人人平等表现在哪些方面| 亚洲成人中文字幕在线播放| 午夜福利成人在线免费观看| 精品国产三级普通话版| 丰满人妻熟妇乱又伦精品不卡| 亚洲片人在线观看| 成人18禁在线播放| 免费电影在线观看免费观看| 男女那种视频在线观看| 长腿黑丝高跟| 国产1区2区3区精品| 日本精品一区二区三区蜜桃| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 成人一区二区视频在线观看| 亚洲精品国产精品久久久不卡| 99热精品在线国产| 国产成人福利小说| 午夜日韩欧美国产| 国产精品女同一区二区软件 | 高清毛片免费观看视频网站| 国产午夜福利久久久久久| 国产精华一区二区三区| 精品国产超薄肉色丝袜足j| 桃红色精品国产亚洲av| 黄色 视频免费看| 香蕉国产在线看| 身体一侧抽搐| bbb黄色大片| 久久精品人妻少妇| 日本黄色片子视频| 欧美在线黄色| 日本三级黄在线观看| 国产伦精品一区二区三区视频9 | 小蜜桃在线观看免费完整版高清| 听说在线观看完整版免费高清| 久久精品国产综合久久久| 噜噜噜噜噜久久久久久91| 亚洲成av人片在线播放无| 黄色视频,在线免费观看| 好看av亚洲va欧美ⅴa在| 免费看美女性在线毛片视频| 久久精品亚洲精品国产色婷小说| 亚洲国产精品999在线| 成人性生交大片免费视频hd| 国产精品久久久久久亚洲av鲁大| 一级作爱视频免费观看| 午夜两性在线视频| 嫁个100分男人电影在线观看| 日韩中文字幕欧美一区二区| 成人性生交大片免费视频hd| 精品人妻1区二区| 欧美色欧美亚洲另类二区| 日本与韩国留学比较| 国产单亲对白刺激| 亚洲美女视频黄频| 99热这里只有是精品50| 亚洲中文av在线| 一区二区三区国产精品乱码| 嫩草影院精品99| 成人永久免费在线观看视频| 69av精品久久久久久| 黄色日韩在线| 中国美女看黄片| 久久精品91蜜桃| 真人一进一出gif抽搐免费| 成人av一区二区三区在线看| bbb黄色大片| 亚洲精品一区av在线观看| 少妇的丰满在线观看| 搡老岳熟女国产| 美女扒开内裤让男人捅视频| 国产91精品成人一区二区三区| 亚洲av成人精品一区久久| 床上黄色一级片| 午夜日韩欧美国产| 日日干狠狠操夜夜爽| 激情在线观看视频在线高清| 亚洲五月婷婷丁香| 又粗又爽又猛毛片免费看| 99热6这里只有精品| 午夜影院日韩av| 成人三级做爰电影| 在线永久观看黄色视频| 欧美最黄视频在线播放免费| av在线天堂中文字幕| 成人高潮视频无遮挡免费网站| 亚洲国产欧美一区二区综合| 夜夜看夜夜爽夜夜摸| 别揉我奶头~嗯~啊~动态视频| 此物有八面人人有两片| 婷婷丁香在线五月| 美女黄网站色视频| 成人精品一区二区免费| 日本a在线网址| 一二三四在线观看免费中文在| 国产亚洲av嫩草精品影院| 国产精品综合久久久久久久免费| 欧美丝袜亚洲另类 | 中文字幕最新亚洲高清| 19禁男女啪啪无遮挡网站| 午夜日韩欧美国产| 亚洲aⅴ乱码一区二区在线播放| 又紧又爽又黄一区二区| 亚洲人成电影免费在线| 中文字幕熟女人妻在线| 欧美日韩中文字幕国产精品一区二区三区| 在线观看免费午夜福利视频| 免费看光身美女| 俄罗斯特黄特色一大片| 亚洲第一电影网av| 人人妻人人看人人澡| 久久久久久国产a免费观看| 国产伦精品一区二区三区四那| 国产黄a三级三级三级人| 精品国内亚洲2022精品成人| 男女床上黄色一级片免费看| www.熟女人妻精品国产| 亚洲国产精品sss在线观看| 午夜精品久久久久久毛片777| 欧美三级亚洲精品| 欧美xxxx黑人xx丫x性爽| 亚洲av五月六月丁香网| 狠狠狠狠99中文字幕| 无限看片的www在线观看| 成在线人永久免费视频| 国产av麻豆久久久久久久| 亚洲av成人av| 国产 一区 欧美 日韩| 国产精品久久久久久精品电影| 亚洲avbb在线观看| 啦啦啦免费观看视频1| www国产在线视频色| 午夜福利18| 老司机午夜十八禁免费视频| 日韩免费av在线播放| 午夜精品在线福利| 亚洲精品在线美女| 午夜精品在线福利| 色哟哟哟哟哟哟| 一边摸一边抽搐一进一小说| 亚洲自偷自拍图片 自拍| 国产午夜精品久久久久久| 制服人妻中文乱码| 在线永久观看黄色视频| www.精华液| 精品国产三级普通话版| 精品99又大又爽又粗少妇毛片 | 美女高潮喷水抽搐中文字幕| 丰满人妻熟妇乱又伦精品不卡| 男女做爰动态图高潮gif福利片| 美女扒开内裤让男人捅视频| 国产精品香港三级国产av潘金莲| 91麻豆av在线| 成人三级黄色视频| 老司机午夜福利在线观看视频| 国产精品久久久av美女十八| 少妇的丰满在线观看| 日韩欧美精品v在线| 女人被狂操c到高潮| 亚洲欧美日韩无卡精品| 一卡2卡三卡四卡精品乱码亚洲| 免费电影在线观看免费观看| 中文资源天堂在线| 俺也久久电影网| 人人妻,人人澡人人爽秒播| 亚洲精品乱码久久久v下载方式 | 精品无人区乱码1区二区| 在线观看一区二区三区| 怎么达到女性高潮| 最近最新中文字幕大全免费视频| 人妻夜夜爽99麻豆av| 老司机午夜福利在线观看视频| 日韩欧美三级三区| 欧美日韩国产亚洲二区| 美女cb高潮喷水在线观看 | 亚洲欧洲精品一区二区精品久久久| 女生性感内裤真人,穿戴方法视频| 国产精品99久久久久久久久| 久久九九热精品免费| 无遮挡黄片免费观看| 中文字幕最新亚洲高清| 极品教师在线免费播放| 成人特级黄色片久久久久久久| 99久久精品国产亚洲精品| 国产精品久久久久久精品电影| 色综合婷婷激情| 欧美xxxx黑人xx丫x性爽| 日韩精品青青久久久久久| 国产精品 欧美亚洲| 色综合站精品国产| 久久国产精品影院| 97超级碰碰碰精品色视频在线观看| 国产精品久久久人人做人人爽| 在线国产一区二区在线| www.www免费av| 2021天堂中文幕一二区在线观| 亚洲天堂国产精品一区在线| 亚洲九九香蕉| 国产一级毛片七仙女欲春2| 少妇的丰满在线观看| www.999成人在线观看| 国产成+人综合+亚洲专区| 青草久久国产| www日本黄色视频网| 免费在线观看日本一区| 91久久精品国产一区二区成人 | 国内揄拍国产精品人妻在线| 99国产精品99久久久久| 成人无遮挡网站| 国产精品影院久久| 最近最新中文字幕大全电影3| 成人特级av手机在线观看| 亚洲专区国产一区二区| 日本撒尿小便嘘嘘汇集6| 丰满人妻一区二区三区视频av | 国产极品精品免费视频能看的| 久久草成人影院| 国产成人福利小说| 一二三四在线观看免费中文在| 免费在线观看成人毛片| 99热只有精品国产| 欧美日韩瑟瑟在线播放| 18禁黄网站禁片午夜丰满| 黄色女人牲交| 一本综合久久免费| 99在线视频只有这里精品首页| 最新中文字幕久久久久 | 免费高清视频大片| 白带黄色成豆腐渣| aaaaa片日本免费| 18禁国产床啪视频网站| 一个人看视频在线观看www免费 | 看片在线看免费视频| 国产成人影院久久av| 亚洲真实伦在线观看| 国产高清三级在线| 亚洲av第一区精品v没综合| 久久久久久人人人人人| 欧美3d第一页| 国产免费av片在线观看野外av| 人妻丰满熟妇av一区二区三区| 亚洲无线在线观看| 91九色精品人成在线观看| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 国产黄色小视频在线观看| 婷婷亚洲欧美| 天天躁日日操中文字幕| 在线观看66精品国产| 午夜福利视频1000在线观看| 日韩免费av在线播放| 欧美中文日本在线观看视频| 黄色女人牲交| 99热精品在线国产| 香蕉丝袜av| 一二三四在线观看免费中文在| tocl精华| 午夜影院日韩av| www.熟女人妻精品国产| 国产av不卡久久| 成人永久免费在线观看视频| 精品久久久久久久末码| 99久久久亚洲精品蜜臀av| 最新在线观看一区二区三区| 99久久精品热视频| 美女被艹到高潮喷水动态| 日本免费a在线| 欧美日韩综合久久久久久 | 国产成人一区二区三区免费视频网站| 丁香欧美五月| 真实男女啪啪啪动态图| 视频区欧美日本亚洲| 日韩国内少妇激情av| 免费看十八禁软件| 熟女电影av网|