被廣泛應用于遺傳基因定位的另一種方法，BSA，也可有效解析質(zhì)量與數(shù)量性狀。高通量測序技術(shù)的快速發(fā)展使BSA在基因定位中發(fā)揮著巨大潛力，例如，MutMap，SHOREmap與MMAPPR在質(zhì)量性狀定位中可取代傳統(tǒng)的圖位克隆。對于數(shù)量性狀，可利用QTL-seq定位QTLs，但定位區(qū)間較大很難識別候選基因，所以非常需要一種新的方法，用于快速定位數(shù)量性狀候選基因。
介于上述種種定位方法的局限性，本文介紹一種新的方法——QTG-Seq，即通過QTL分離（即，將數(shù)量性狀轉(zhuǎn)化為質(zhì)量性狀），極端表型混池，高通量測序以及新算法挖掘候選基因。這一方法是由華中農(nóng)業(yè)大學李林課題組、章元明課題組與中國農(nóng)科院王國英課題組開發(fā)并應用成功，最后以“QTG-seq accelerates QTL fine mapping through QTL partitioning and whole-genome sequencing on bulked segregant samples”為題，在線發(fā)表在Molecular Plant上。
QTG-Seq基本原理（流程如下）——
（1）群體構(gòu)建：?F1，?F2，?BC1F1群體
（2）QTL定位：F2群體定位QTL，F(xiàn)2:3家系證明QTL定位結(jié)果
（3）QTL分離：分子標記篩選目標QTL雜合，其它QTL純合的BC1F1單株（BC1F1數(shù)目要充足）；然后自交得到BC1F1:2家系
（4）極端表型選擇：從BC1F1:2中選2個極端表型組（前20%，后20%，至少1000個單株）
（5）極端表型混池：提取DNA并混池，形成2個極端池
（6）測序并變異calling：對兩個混池測序，與參考基因組比對分析變異
（7）關(guān)聯(lián)分析：新方法smoothLOD（準確關(guān)聯(lián)靶標位點），結(jié)合其他方法，如ED（Hill et al., 2013）與G’（Magwene et al., 2011）（關(guān)聯(lián)peak位置）
?圖1.?QTG-Seq基本原理
應用——QTG-Seq快速定位玉米株高QTL，qPH7
以玉米株高為模式性狀展開QTG-seq研究：利用4代玉米群體材料，通過QTL-seq精細定位株高主效QTL，這與傳統(tǒng)QTL精細定位用到許多代創(chuàng)制高世代回交群體明顯不同。首先，利用株高差異明顯的兩個自交系HZS與1462，構(gòu)建F2分離群體，個體株高變化范圍為199cm~307 cm，呈超親分離的現(xiàn)象，暗示株高是由多個QTLs控制；然后，利用1028個多態(tài)性標記定位主效QTLs，結(jié)果定位到4個QTL（位于Chr.1,3,6,7），PVE介于7%-17%，又利用F2:3家系驗證QTL定位結(jié)果可靠；最后，選擇Chr7.上的QTL?qPH7利用QTL-seq策略進行精細定位。
利用12個標記篩選出813個BC1F1單株（qPH7位點雜合，另外3個QTL位點純合），從中選擇15個材料，自交獲得BC1F1:2家系，并考察株高。共獲得3120個BC1F1:2單株，從中劃分出2個極端表型組，低值組580個單株，高值組567個單株，并分別組成混池（低池與高池）；接著，對混池展開全基因組測序（＞280×），檢測到197021個高質(zhì)量SNPs。利用smoothLOD與ED4算法分析出Chr.7只有一個峰，并且峰的位置在135.3Mb。KASP標記分型也證實峰的位置處于135.1與135.2之間。有趣的是，smoothLOD法將qPH7定位在chr：135216475bp，此位點正好處在Zm00001d020874上。Zm00001d020874編碼一個NF-YC結(jié)構(gòu)域蛋白，它的擬南芥同源蛋白與開花時間有關(guān)。
圖2.?qPH7精細定位與候選基因挖掘
慮到分離群體中重組blocks的存在，利用增加滑窗大小的方法識別目標區(qū)域的重組blocks。由于smoothLOD與G’為連續(xù)變量，所以采用ED/SNP的方法分析重組blocks，結(jié)果顯示候選基因與peak間距離＜150kb。因此，候選基因的區(qū)間在300?kb，其中包含13個基因，也包括基因Zm00001d020874。
因為交換單株數(shù)量與定位區(qū)間大小密切相關(guān)，所以確保較大的混池規(guī)模至關(guān)重要。根據(jù)模擬數(shù)據(jù)，QTG-seq相對低的測序深度（＜100×）足夠檢測出所有的重組blocks，因為每個重組block包含許多SNPs，沒必要對每個重組block的所有reads測序。因此，為了平衡定位精度與測序成本，在使用QTG-seq方法時，建議增大混池規(guī)模與相對低的測序深度。
為進一步確定候選基因，本研究又通過RNA-seq分析發(fā)現(xiàn)，Zm00001d020874表達量在雙親的莖頂端分生組織（SAM）與幼嫩節(jié)間組織中存在顯著性差異。基于擬南芥同源蛋白功能與smoothLOD定位信號，認為Zm00001d020874是qPH7的候選基因，并且后續(xù)的CRISPR基因敲除，蛋白互作驗證等實驗均暗示Zm00001d020874是控制株高的基因。
圖3.?候選基因的驗證
Monte Carlo模擬支持QTG-Seq的可靠性
4個Monte Carlo模擬實驗如下：
1.?PVE與QTG位置偏差大?。?%時，接近80kb；≥10%時，30-60kb
2.?取樣率與QTG位置偏差大?。寒斎勇蕿?0%時，位置偏差略大；較大的取樣比例會增加群體大小，但有效重組事件數(shù)量也會增加
3.?群體大小與QTG位置偏差大?。寒斎后w大于4000時，位置偏差小于60kb
4.?測序深度與檢測力：當測序深度大于100×時（混池超過1000個個體），檢測力達到100%
圖4.?SmoothLOD算法影響定位準確性的各種因素
QTG-Seq特點總結(jié)：
1.?省時省力省錢，不需要構(gòu)建高世代復雜群體，只需F1，?F2，?BC1F1群體即可實現(xiàn)；
2.?將多個QTL分解為單個QTL分析（數(shù)量性狀→質(zhì)量性狀），即，在每個QTL區(qū)域選擇3個（或多個）多態(tài)性標記篩選目標QTL雜合，其它QTL純合的交換單株，確保交換單株后續(xù)自交發(fā)生表型分離以對目標QTL進行解析；
3.?測序策略：高的混池個體數(shù)目（數(shù)百至上千個個體）與相對低的測序深度（~100×）；相對低的測序深度，是因為每條重組block上有多個SNPs，檢測其中1個SNP就可代表；
4.?關(guān)聯(lián)分析采用新方法smoothLOD，并結(jié)合ED與G’，定位精度高，可識別候選基因；smoothLOD可準確關(guān)聯(lián)靶標位點，而ED與G’關(guān)聯(lián)peak位置；
5.?可結(jié)合KASP標記分型尋找重組斷點等，識別候選基因。
參考文獻:
Zhang H., Wang X., Pan Q., Li P., Liu Y., Lu X., Zhong W., Li M., Han L., Li J., Wang P., Li D., Liu Y., Li Q., Yang F., Zhang Y.-M., Wang G., and Li L. QTG-Seq Accelerates QTL Fine Mapping through QTL Partitioning and Whole-Genome Sequencing of Bulked Segregant Samples.?Mol Plant.?2019 Mar 4;12(3):426-437.
如果您的科研項目有問題，歡迎點擊下方按鈕咨詢我們，我們將免費為您設(shè)計文章方案。
]]>
在這40余篇文章中，遺傳圖譜文章就有21篇，占據(jù)了優(yōu)勢，GWAS有8篇，遺傳進化也有6篇，BSA文章也有8篇之多。
所研究的物種也有27個，包括：
作物類：大豆、大麥、甘薯、高粱、花生、棉花、水稻、油菜、芝麻、甜菜、鴨茅、黃麻等；林木類：桂花、旱柳、梨、蘋果、葡萄、銀杏等；
水產(chǎn)類：中華絨螯蟹、白對蝦、三疣梭子蟹、沙塘鱧、螠蟶等；
果蔬類：甜瓜、辣椒、菠菜等；
家禽類：豬等。
所研究的方向：抗病性、產(chǎn)量、品質(zhì)、開花期、花色、株型、性別決定、雄性不育等相關(guān)功能基因的定位，另外還有多個物種的遺傳多樣性和進化分析。
百邁客自2009年成立，8年來深耕群體遺傳研究領(lǐng)域，目前已完成上百個物種的測序分析項目，項目經(jīng)驗十分豐富，共發(fā)表文章200+篇，累計影響因子達650+， 其中包含多篇Nature Genetics、Nature Communications及PNAS等國際等級期刊文章。百邁客致力于群體遺傳研究，為您提供專業(yè)的測序分析服務，助力您的科學研究！
百邁客2017年文章列表
