• <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ī)律
    免费搜索国产男女视频| 久久精品影院6| 亚洲片人在线观看| 欧美一区二区精品小视频在线| 国产免费一级a男人的天堂| 两个人的视频大全免费| 国产野战对白在线观看| 亚洲成人中文字幕在线播放| 啦啦啦观看免费观看视频高清| 无限看片的www在线观看| 琪琪午夜伦伦电影理论片6080| 麻豆久久精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 动漫黄色视频在线观看| 韩国av一区二区三区四区| 在线观看免费午夜福利视频| www日本黄色视频网| 成年版毛片免费区| 国产欧美日韩精品一区二区| 国产v大片淫在线免费观看| 别揉我奶头~嗯~啊~动态视频| 国产精品 国内视频| 怎么达到女性高潮| 国产伦一二天堂av在线观看| 日韩中文字幕欧美一区二区| av天堂在线播放| 变态另类丝袜制服| 久久久久国内视频| 女警被强在线播放| 中出人妻视频一区二区| 波多野结衣巨乳人妻| 老司机深夜福利视频在线观看| 亚洲国产精品合色在线| 国产高清视频在线播放一区| eeuss影院久久| 草草在线视频免费看| 国产一区二区激情短视频| 2021天堂中文幕一二区在线观| 国产精品香港三级国产av潘金莲| 亚洲av电影在线进入| 综合色av麻豆| 国产av不卡久久| www日本黄色视频网| 无遮挡黄片免费观看| 亚洲片人在线观看| 最近视频中文字幕2019在线8| 亚洲国产欧洲综合997久久,| 欧美一区二区亚洲| 欧美成人性av电影在线观看| 18美女黄网站色大片免费观看| 久久久久久国产a免费观看| 国产精品久久久久久精品电影| 午夜福利高清视频| 欧美+亚洲+日韩+国产| 欧美bdsm另类| 欧美一区二区国产精品久久精品| 桃红色精品国产亚洲av| 51国产日韩欧美| 欧美极品一区二区三区四区| 韩国av一区二区三区四区| 亚洲专区国产一区二区| 国产精品影院久久| 久久久久久久久大av| 国产亚洲欧美在线一区二区| 一二三四社区在线视频社区8| 变态另类丝袜制服| 亚洲精品在线观看二区| 精品国产亚洲在线| 岛国在线观看网站| 精品电影一区二区在线| 日本精品一区二区三区蜜桃| 观看美女的网站| 免费看十八禁软件| 成年女人毛片免费观看观看9| 成熟少妇高潮喷水视频| 波多野结衣巨乳人妻| 在线观看免费视频日本深夜| eeuss影院久久| 久久国产精品影院| aaaaa片日本免费| 真人做人爱边吃奶动态| 国产高清视频在线观看网站| 久久精品91蜜桃| 色视频www国产| 国产一区二区三区视频了| 国产v大片淫在线免费观看| a级毛片a级免费在线| 国产精品影院久久| 淫秽高清视频在线观看| 欧美又色又爽又黄视频| 免费观看人在逋| 欧美高清成人免费视频www| 亚洲人成网站在线播放欧美日韩| 嫩草影院入口| 黄色日韩在线| 欧美丝袜亚洲另类 | www.www免费av| 免费一级毛片在线播放高清视频| 99热6这里只有精品| 免费高清视频大片| 亚洲精品在线美女| 天堂影院成人在线观看| 国产成年人精品一区二区| 日韩国内少妇激情av| 美女高潮的动态| 我要搜黄色片| 丰满人妻熟妇乱又伦精品不卡| 在线观看一区二区三区| 欧美国产日韩亚洲一区| 久久精品91无色码中文字幕| 欧美bdsm另类| 亚洲在线观看片| 久久这里只有精品中国| 中文字幕高清在线视频| 亚洲七黄色美女视频| eeuss影院久久| 岛国在线观看网站| 国产一区二区亚洲精品在线观看| 色哟哟哟哟哟哟| 俺也久久电影网| 国产蜜桃级精品一区二区三区| 亚洲国产精品久久男人天堂| 一个人看视频在线观看www免费 | 不卡一级毛片| 亚洲国产精品合色在线| 两个人视频免费观看高清| 一个人免费在线观看电影| 91av网一区二区| 成人特级黄色片久久久久久久| 亚洲欧美精品综合久久99| 久久人人精品亚洲av| 精品免费久久久久久久清纯| 欧美色视频一区免费| 日韩免费av在线播放| 99视频精品全部免费 在线| 国产私拍福利视频在线观看| 久久久久国产精品人妻aⅴ院| 久久久精品大字幕| 给我免费播放毛片高清在线观看| 欧美成人免费av一区二区三区| 日韩成人在线观看一区二区三区| 一进一出抽搐动态| 成人无遮挡网站| 亚洲人成电影免费在线| 99国产精品一区二区蜜桃av| 日韩成人在线观看一区二区三区| 尤物成人国产欧美一区二区三区| 日日夜夜操网爽| 99在线人妻在线中文字幕| 亚洲av成人不卡在线观看播放网| 免费看十八禁软件| www.999成人在线观看| 欧美黄色片欧美黄色片| 五月玫瑰六月丁香| 天堂动漫精品| 欧美一区二区国产精品久久精品| 深爱激情五月婷婷| 99国产精品一区二区蜜桃av| 免费搜索国产男女视频| 18禁裸乳无遮挡免费网站照片| 午夜免费激情av| 久99久视频精品免费| 少妇的丰满在线观看| 啦啦啦韩国在线观看视频| 又黄又爽又免费观看的视频| 可以在线观看毛片的网站| 精品人妻一区二区三区麻豆 | 又黄又粗又硬又大视频| 色精品久久人妻99蜜桃| 两个人的视频大全免费| 亚洲国产中文字幕在线视频| 久久久久久九九精品二区国产| 搞女人的毛片| 看黄色毛片网站| svipshipincom国产片| 级片在线观看| 精品无人区乱码1区二区| 国产午夜福利久久久久久| 特大巨黑吊av在线直播| 国产成年人精品一区二区| 亚洲av免费高清在线观看| 亚洲av第一区精品v没综合| 91av网一区二区| 免费看a级黄色片| 99国产精品一区二区三区| 高潮久久久久久久久久久不卡| 国产v大片淫在线免费观看| 国产亚洲欧美在线一区二区| 精品一区二区三区人妻视频| 99久久成人亚洲精品观看| 一进一出抽搐动态| 色综合亚洲欧美另类图片| 在线观看一区二区三区| 欧美+亚洲+日韩+国产| 免费av不卡在线播放| 精品国产亚洲在线| 午夜日韩欧美国产| 日日夜夜操网爽| 高潮久久久久久久久久久不卡| 国产成人啪精品午夜网站| 免费在线观看成人毛片| 琪琪午夜伦伦电影理论片6080| 久久久久久久久久黄片| 久久欧美精品欧美久久欧美| 欧美一区二区精品小视频在线| 99国产极品粉嫩在线观看| 亚洲天堂国产精品一区在线| 男人和女人高潮做爰伦理| 在线观看日韩欧美| 国产三级黄色录像| 男人舔女人下体高潮全视频| 此物有八面人人有两片| 2021天堂中文幕一二区在线观| 老熟妇乱子伦视频在线观看| 欧美不卡视频在线免费观看| 国产毛片a区久久久久| 亚洲av五月六月丁香网| 少妇高潮的动态图| 亚洲色图av天堂| 成人永久免费在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 中出人妻视频一区二区| 精品乱码久久久久久99久播| 人人妻,人人澡人人爽秒播| 99热精品在线国产| 亚洲无线观看免费| 国产淫片久久久久久久久 | 国内精品久久久久精免费| 国产黄a三级三级三级人| 国产亚洲精品久久久久久毛片| 丰满人妻熟妇乱又伦精品不卡| 丝袜美腿在线中文| 亚洲精品色激情综合| 搡女人真爽免费视频火全软件 | 搡老妇女老女人老熟妇| 国产一区二区三区视频了| 欧美色欧美亚洲另类二区| 日本撒尿小便嘘嘘汇集6| 久久久国产成人免费| 国产91精品成人一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 岛国在线观看网站| 日韩国内少妇激情av| 特大巨黑吊av在线直播| 看免费av毛片| av女优亚洲男人天堂| 国产欧美日韩一区二区精品| 亚洲无线在线观看| 97人妻精品一区二区三区麻豆| 窝窝影院91人妻| 天堂av国产一区二区熟女人妻| 国产精品 欧美亚洲| av黄色大香蕉| 精品一区二区三区av网在线观看| 亚洲第一电影网av| 欧美极品一区二区三区四区| 国产aⅴ精品一区二区三区波| 波多野结衣巨乳人妻| 亚洲精品国产精品久久久不卡| 日韩欧美免费精品| 欧美zozozo另类| 精品熟女少妇八av免费久了| xxxwww97欧美| 亚洲男人的天堂狠狠| 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| 天天添夜夜摸| 久久午夜亚洲精品久久| 欧洲精品卡2卡3卡4卡5卡区| 国产综合懂色| 国产欧美日韩精品一区二区| 国内精品久久久久精免费| 可以在线观看的亚洲视频| 美女高潮的动态| 日本免费a在线| 好男人在线观看高清免费视频| 亚洲精品乱码久久久v下载方式 | 精品无人区乱码1区二区| 午夜福利成人在线免费观看| av专区在线播放| 国产毛片a区久久久久| 伊人久久大香线蕉亚洲五| 日本 av在线| 久久这里只有精品中国| 热99在线观看视频| 国产一区在线观看成人免费| 哪里可以看免费的av片| 变态另类丝袜制服| 长腿黑丝高跟| 成人特级av手机在线观看| 此物有八面人人有两片| 国产午夜精品久久久久久一区二区三区 | 成人欧美大片| 九九在线视频观看精品| 老司机深夜福利视频在线观看| 国产激情欧美一区二区| 女同久久另类99精品国产91| 婷婷亚洲欧美| 美女被艹到高潮喷水动态| 欧美激情久久久久久爽电影| 精品一区二区三区av网在线观看| 国产激情偷乱视频一区二区| 日本成人三级电影网站| 亚洲精品日韩av片在线观看 | 别揉我奶头~嗯~啊~动态视频| 日本免费a在线| 少妇人妻精品综合一区二区 | 我的老师免费观看完整版| 人妻夜夜爽99麻豆av| 欧美一区二区精品小视频在线| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av涩爱 | 天堂影院成人在线观看| 国产成人欧美在线观看| 日韩精品中文字幕看吧| 日韩免费av在线播放| 中文字幕熟女人妻在线| 成年人黄色毛片网站| 精品一区二区三区av网在线观看| 97碰自拍视频| 亚洲国产精品成人综合色| 亚洲人成伊人成综合网2020| 免费在线观看成人毛片| 久久性视频一级片| 麻豆国产av国片精品| 18禁裸乳无遮挡免费网站照片| 国产成人影院久久av| 午夜精品久久久久久毛片777| 亚洲内射少妇av| 亚洲av一区综合| 亚洲五月婷婷丁香| 国产97色在线日韩免费| 日韩欧美在线二视频| 欧美一级a爱片免费观看看| 伊人久久精品亚洲午夜| 全区人妻精品视频| 伊人久久大香线蕉亚洲五| 97人妻精品一区二区三区麻豆| 亚洲美女黄片视频| 国产一区二区激情短视频| 精品免费久久久久久久清纯| 在线播放国产精品三级| 久久久久久九九精品二区国产| 欧美性猛交黑人性爽| xxxwww97欧美| 夜夜看夜夜爽夜夜摸| 麻豆国产97在线/欧美| 亚洲一区高清亚洲精品| 夜夜看夜夜爽夜夜摸| 欧美日韩中文字幕国产精品一区二区三区| 亚洲男人的天堂狠狠| 亚洲av中文字字幕乱码综合| 国产一区二区在线观看日韩 | 美女被艹到高潮喷水动态| 亚洲,欧美精品.| 欧美一级毛片孕妇| 少妇丰满av| 欧美丝袜亚洲另类 | 国产精品国产高清国产av| 99国产综合亚洲精品| 免费一级毛片在线播放高清视频| 亚洲熟妇中文字幕五十中出| 国内精品久久久久精免费| 日本 av在线| 美女高潮的动态| ponron亚洲| 国产淫片久久久久久久久 | 国产精品影院久久| 中文字幕精品亚洲无线码一区| 日韩欧美 国产精品| 国产欧美日韩一区二区精品| 亚洲电影在线观看av| 国产精品 国内视频| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放| 18+在线观看网站| 天堂动漫精品| 一级毛片女人18水好多| 亚洲国产欧洲综合997久久,| 黑人欧美特级aaaaaa片| 亚洲18禁久久av| 亚洲人成网站高清观看| 18禁国产床啪视频网站| 国产av一区在线观看免费| 亚洲美女视频黄频| 亚洲不卡免费看| 欧美在线黄色| 欧美一级毛片孕妇| 色噜噜av男人的天堂激情| 国产中年淑女户外野战色| 成人av一区二区三区在线看| 国产淫片久久久久久久久 | 日韩人妻高清精品专区| 亚洲欧美日韩高清专用| 身体一侧抽搐| 一级作爱视频免费观看| 国产高清视频在线播放一区| 亚洲国产色片| 亚洲 国产 在线| 真人一进一出gif抽搐免费| 亚洲中文字幕日韩| 亚洲五月婷婷丁香| 熟女电影av网| 久久性视频一级片| 不卡一级毛片| 18禁裸乳无遮挡免费网站照片| 小说图片视频综合网站| 美女cb高潮喷水在线观看| av福利片在线观看| 国产精品自产拍在线观看55亚洲| netflix在线观看网站| 欧美性猛交黑人性爽| 观看免费一级毛片| 国产精品99久久99久久久不卡| 国产精品乱码一区二三区的特点| 亚洲在线观看片| 国产午夜精品论理片| 动漫黄色视频在线观看| 最新中文字幕久久久久| 舔av片在线| 久久久国产成人免费| 国产精品永久免费网站| 好男人电影高清在线观看| 在线观看免费午夜福利视频| 成人国产一区最新在线观看| 美女 人体艺术 gogo| 免费看a级黄色片| 日本撒尿小便嘘嘘汇集6| 午夜福利在线在线| 成年免费大片在线观看| 日本免费一区二区三区高清不卡| 51国产日韩欧美| 午夜福利在线观看免费完整高清在 | 成人永久免费在线观看视频| 舔av片在线| 国产色婷婷99| 精品熟女少妇八av免费久了| 日本一本二区三区精品| 熟女电影av网| 91av网一区二区| 国产精品女同一区二区软件 | 成年女人看的毛片在线观看| 草草在线视频免费看| 午夜精品久久久久久毛片777| 波多野结衣巨乳人妻| 在线视频色国产色| 欧美成人性av电影在线观看| 日本在线视频免费播放| 国产精品久久久久久久电影 | 99在线视频只有这里精品首页| 亚洲在线观看片| 欧美又色又爽又黄视频| 一二三四社区在线视频社区8| 最后的刺客免费高清国语| 午夜精品一区二区三区免费看| 波多野结衣高清无吗| 69人妻影院| 成熟少妇高潮喷水视频| 久久精品亚洲精品国产色婷小说| 在线观看一区二区三区| 国产精品嫩草影院av在线观看 | 五月伊人婷婷丁香| 欧美成狂野欧美在线观看| 18+在线观看网站| av欧美777| 欧美成人一区二区免费高清观看| 啦啦啦韩国在线观看视频| 免费看a级黄色片| 18美女黄网站色大片免费观看| 久久香蕉精品热| 好男人在线观看高清免费视频| 国产免费一级a男人的天堂| 99精品久久久久人妻精品| 久久久久久人人人人人| 成人av在线播放网站| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久末码| 老熟妇乱子伦视频在线观看| 午夜福利视频1000在线观看| 成人av一区二区三区在线看| 午夜免费激情av| 成年版毛片免费区| 俺也久久电影网| 日韩人妻高清精品专区| 色视频www国产| 最新中文字幕久久久久| 日韩大尺度精品在线看网址| 国产精品爽爽va在线观看网站| 亚洲一区二区三区不卡视频| 国模一区二区三区四区视频| 亚洲成人久久性| 国产97色在线日韩免费| 日韩欧美一区二区三区在线观看| 国内精品久久久久精免费| 18禁在线播放成人免费| 亚洲色图av天堂| 久久久久国内视频| 91字幕亚洲| 成人亚洲精品av一区二区| 香蕉av资源在线| 少妇的逼水好多| 两个人视频免费观看高清| 中文亚洲av片在线观看爽| 非洲黑人性xxxx精品又粗又长| 亚洲精品日韩av片在线观看 | 男女午夜视频在线观看| 亚洲男人的天堂狠狠| 久久久久久久久久黄片| 18+在线观看网站| 欧美成人a在线观看| 国产午夜福利久久久久久| 国产亚洲欧美98| 日韩欧美一区二区三区在线观看| 99热精品在线国产| 成人午夜高清在线视频| 亚洲中文字幕日韩| 无人区码免费观看不卡| 国内揄拍国产精品人妻在线| 精品国内亚洲2022精品成人| 性色av乱码一区二区三区2| 日本黄色视频三级网站网址| 非洲黑人性xxxx精品又粗又长| 色噜噜av男人的天堂激情| 亚洲在线观看片| or卡值多少钱| 婷婷精品国产亚洲av| 一进一出抽搐动态| 美女大奶头视频| 亚洲在线观看片| 香蕉丝袜av| 99视频精品全部免费 在线| 亚洲欧美激情综合另类| 久久精品夜夜夜夜夜久久蜜豆| 国产色婷婷99| 51午夜福利影视在线观看| 一级黄色大片毛片| aaaaa片日本免费| 久久久成人免费电影| 国产精品电影一区二区三区| 欧美在线一区亚洲| 精品一区二区三区人妻视频| 美女高潮的动态| 国产精品一及| 午夜免费男女啪啪视频观看 | 国内精品美女久久久久久| 亚洲欧美激情综合另类| 丰满人妻一区二区三区视频av | 一二三四社区在线视频社区8| 国产午夜精品论理片| 精品久久久久久久末码| 欧美+亚洲+日韩+国产| 99久久成人亚洲精品观看| 亚洲精品粉嫩美女一区| 又黄又爽又免费观看的视频| 欧美乱妇无乱码| 国产真实伦视频高清在线观看 | 国内精品久久久久久久电影| 国产成人系列免费观看| 亚洲男人的天堂狠狠| 在线免费观看的www视频| 脱女人内裤的视频| 亚洲无线在线观看| 搞女人的毛片| 国产精品 国内视频| 国产午夜福利久久久久久| 欧美日韩一级在线毛片| 国产男靠女视频免费网站| 国产v大片淫在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美国产在线观看| 亚洲av美国av| 免费观看精品视频网站| 琪琪午夜伦伦电影理论片6080| 中文字幕av在线有码专区| 中文字幕熟女人妻在线| 又黄又粗又硬又大视频| 搡老妇女老女人老熟妇| 精品福利观看| tocl精华| 又粗又爽又猛毛片免费看| 在线看三级毛片| 日韩 欧美 亚洲 中文字幕| 欧美国产日韩亚洲一区| 一二三四社区在线视频社区8| 天堂av国产一区二区熟女人妻| or卡值多少钱| 欧美精品啪啪一区二区三区| 嫩草影院入口| 我的老师免费观看完整版| 成人性生交大片免费视频hd| 99国产极品粉嫩在线观看| 成人亚洲精品av一区二区| 国产欧美日韩一区二区三| 天堂网av新在线| 嫩草影院精品99| 精品福利观看| 免费在线观看成人毛片| 中文字幕人妻丝袜一区二区| 国产精品乱码一区二三区的特点| 国产美女午夜福利| 国产欧美日韩一区二区三| 日韩国内少妇激情av| 欧美日韩福利视频一区二区| 国产免费av片在线观看野外av| 91九色精品人成在线观看| 久9热在线精品视频| 久久精品国产自在天天线| 欧美高清成人免费视频www| 男插女下体视频免费在线播放| 3wmmmm亚洲av在线观看| 午夜视频国产福利| 久久久久九九精品影院| 91av网一区二区| 人人妻人人看人人澡| 亚洲av免费高清在线观看|