參考案例
玉米株高性狀QTG-Seq精細等位案例分享
標題:
QTG-Seq Accelerates QTL Fine Mapping through QTL Partitioning and Whole-Genome Sequencing of Bulked Segregant Samples
期刊:
Molecular Plant (IF: 10.812)
研究背景:
農藝性狀是作物改良育種的直接目標,這一過程對全球糧食安全至關重要。現代作物改良依賴于對農藝性狀的遺傳和分子機制的剖析。圖位克隆是一種有效的分離QTGs的方法,然而其費時費力,對檢測罕見等位基因的能力是有限的。數量性狀位點測序(quantitative trait locus sequencing, QTLseq)分辨率太低,難以識別候選基因,特別是在基因組規模較大的物種中。因此,仍然需要一種精細的策略來快速識別潛在數量性狀的候選基因。玉米是一種被廣泛研究的遺傳模型,是世界上種植最廣泛的作物之一,其株高性狀可以作為QTG-seq方法驗證的性狀模型。
主要研究結果:
QTG-seq策略步驟如下:(1)利用兩個自交系培育F1、F2和BC1F1群體。(2)采用F2群體進行QTL定位,采用F2:3家系確定QTL定位可靠性。(3)使用與每個QTL緊密連鎖的分子標記進行QTL分區,在BC1F1群體中選擇感興趣的QTL的雜合位點,同時選擇其他QTL的純合位點。(4)BC1F1群體自交構建BC1F1:2衍生系,選擇表現出極端表型。(5)構建極端表型混池。(6)混池測序,遺傳變異識別。(7)計算混池的等位基因頻率,并使用新的統計量smooth LOD值,確定目標位點的峰值位置。

QTG策略的研究步驟
候選基因定位:
通過玉米自交系HZS和1462,構建F2分離群體,混池測序,使用LOD閾值3.7,檢測了位于染色體1、3、6和7上的4個QTL,分別解釋了總表型方差的11.22%、7.49%、9.39%和16.60%。衍生F2:3家系的分析,驗證了4個QTL位點的可靠性。同時,選擇LOD值最高的QTL(在F2群體中為11.99),位于7號染色體上(以下簡稱qPH7),采用QTG-seq策略進行精細等位。使用12個分子標記,篩選出813個在qPH7處為雜合子,在其他QTL處為純合子的BC1F1個體。選擇15株,自交,構建BC1F1:2家系,表型分析進行極端混池構建,混池測序(>280×)檢測到197,021個高質量SNP位點。smooth LOD及ED4分析7號染色體上只有一個峰,位置在135,216,475 bp處,屬于Zm00001d020874基因區域,其編碼一個含有NF-YC-結構域的蛋白,其擬南芥同源物已被報道影響開花時間。ED分析表明,qPH7候選基因位于峰值位置周圍的300-kb區間,該區間包含13個基因,也包括Zm00001d020874基因。

qPH7的精細定位和候選基因挖掘
候選基因驗證:
使用RNA-seq,評估兩個親本中300-kb區間內的任何候選基因是否存在表達的顯著差異。在莖頂端分生組織(SAM)和節間分生組織中,Zm00001d020874的表達量均具有顯著的差異。根據其擬南芥同系物在控制開花時間方面的已知功能和由smoothLOD算法識別出的強精細定位信號,我們認為Zm00001d020874是qPH7的強候選信號。為了核實Zm00001d020874在qPH7位點上的功能,構建了一個由RNA指導的CRISPR-Cas9表達載體,以Zm00001d020874外顯子上的兩個區域為目標,對野生型玉米自交系CAL1進行轉化,獲得兩棵轉基因陽性植株T13和T16。轉基因植株株高顯著低于對照組,支持Zm00001d020874是qPH7功能基因的假設。對親本HZS和1462中的Zm00001d020874序列分析發現6個變異位點,對181個從Chinese association mapping panel上獲得的自交系進行RNA測序分析,表明3個基因型之間Zm00001d020874的表達量存在顯著差異,且植株高度和Zm00001d020874的表達量存在顯著相關性。類似的現象在穗高和發育時間參數—花粉脫落日、吐絲日和抽穗日等性狀上發現,表明Zm00001d020874基因在株高上的影響可能與生殖發育有關。選擇足跡分析揭示Zm00001d020874在玉米改良過程中受到選擇

候選基因的驗證
研究意義:
QTG-seq只需要4代就可以達到基因的精細定位,其速度可與質量性狀的基因相媲美。因此,與其他方法相比,QTG-seq是一種更強大的快速QTL精細定位的方法。
參考文獻:
1. Adhikari K, Mendoza-Revilla J, Sohail A, et al. A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia. Nature communications, 2019, 10(1): 1-16.
2. Du X, Huang G, He S, et al. Resequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits. Nature genetics, 2018, 50(6): 796-802.
3. Kaya H B, Akdemir D, Lozano R, et al. Genome wide association study of 5 agronomic traits in olive (Olea europaea L.). Scientific Reports, 2019, 9(1): 1-14.
4. Styrkarsdottir U, Stefansson O A, Gunnarsdottir K, et al. GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures. Nature communications, 2019, 10(1): 1-13.
5. Tamisier L, Szadkowski M, Nemouchi G, et al. Genome‐wide association mapping of QTLs implied in potato virus Y population sizes in pepper: evidence for widespread resistance QTL pyramiding. Molecular plant pathology, 2020, 21(1): 3-16.
6. Wu D, Liang Z, Yan T, et al. Whole-genome resequencing of a worldwide collection of rapeseed accessions reveals the genetic basis of ecotype divergence. Molecular plant, 2019, 12(1): 30-43.