最近到了提交國自然的最后時期，小編猜今年一定還有老師想要申請——解析一個性狀的調(diào)控機制——這類課題，那這么大的課題，第一步我們要干什么呢？那肯定是獲得群體材料了，一般為全世界收集的核心種質(zhì)資源材料組成的自然群體和雙親構(gòu)建遺傳群體，那這些群體材料怎么用呢，這是老師們所關(guān)注的重點。
想要研究一個性狀的調(diào)控機制或一個生物學(xué)現(xiàn)象的原因，通常是先從基因定位開始，而其主流的三種方法是遺傳圖譜、混池分離分析（BSA）和全基因組關(guān)聯(lián)分析（GWAS）。通常三者之間可以聯(lián)合分析，又或是和轉(zhuǎn)錄組進行聯(lián)合分析，將定位區(qū)域縮小或候選基因的數(shù)目減少。
思路1、遺傳圖譜+GWAS
傳統(tǒng)的QTL分析是利用雙親得到的遺傳群體構(gòu)建遺傳圖譜結(jié)合表型對目標(biāo)性狀進行定位，這種方法主要依賴于雙親的遺傳多樣性，且QTL的檢測效果在不同群體間存在差異，同時有的定位結(jié)果可能比較大，包含了較多的基因，無法找到潛在的候選基因進行研究。而利用自然群體的全基因組關(guān)聯(lián)研究(GWAS)可以克服QTL分析的局限性，縮小候選區(qū)域，但GWAS假陽性比QTL分析高，因此兩者聯(lián)合可以在一定程度上彌補彼此的不足。
接下來就跟大家分享一遍辣椒遺傳圖譜和GWAS聯(lián)合分析的文章。
研究意義
材料方法
材料：2個RIL群體（PD群體：120株；TH群體：85株）+1個自然群體（208份）
表型調(diào)查：利用高效液相色譜法（HPLC）測定辣椒堿類物質(zhì)含量（辣椒堿、二氫辣椒堿和總辣椒堿）
1）‘PD’RIL和‘TH’RIL群體：分別進行3個和2個環(huán)境表型鑒定，為了減小果實大小的影響，均為取其果實的胎盤組織用于辣椒堿類物質(zhì)提取。
2）GWAS群體：一年表型鑒定，冷凍干燥整個果實用于辣椒堿類物質(zhì)提取，三次生物學(xué)重復(fù)。
測序手段：簡化測序（TH群體和GWAS材料）和重測序數(shù)據(jù)（PD群體）
技術(shù)路線如下：
主要結(jié)果
QTL定位結(jié)果：利用120個PD單株構(gòu)建的重測序遺傳圖譜和3個環(huán)境下的56個單株表型數(shù)據(jù)定位到22個QTL位點。其中LOD值最高的是PD-dicap1.1 （LOD=8.7），表型變異率（R2）最高的是PD-dicap10.2 (28.8%)。同樣地，85個TH單株構(gòu)建的簡化圖譜和2個環(huán)境的表型數(shù)據(jù)定位到22個QTL位點。其中LOD值（LOD=9.9）和R2（18.8%）最高的都是TH-cap2.2。兩者分析時均用CM334作為參考基因組，因此可以對2個群體的QTL定位結(jié)果進行比較，發(fā)現(xiàn)PD群體的9個QTL位點與TH群體的7個QTL位點的物理位置重疊。
圖1、QTL定位結(jié)果
GWAS結(jié)果：208份材料的簡化數(shù)據(jù)與一年表型數(shù)據(jù)進行GWAS分析定位到99個與辣椒堿顯著相關(guān)的SNP，分別對應(yīng)69個block區(qū)，含有213個基因。
聯(lián)合分析結(jié)果：基于相同的參考基因組，通過比較QTL分析與GWAS分析的定位結(jié)果的物理位置發(fā)現(xiàn)GWAS結(jié)果中10個區(qū)域與本研究QTL定位結(jié)果重疊，查詢注釋結(jié)果，發(fā)現(xiàn)5個基因與辣椒素生物合成途徑相關(guān)基因。
圖2、GWAS定位結(jié)果
思路2、遺傳圖譜+轉(zhuǎn)錄組
轉(zhuǎn)錄組是通過研究多個樣品在不同時期和或不同組織部位之間基因表達的差異，從而找到與性狀相關(guān)的基因，但這種通常是全基因組范圍內(nèi)且不能鎖定性狀相關(guān)區(qū)域。而遺傳圖譜可以鎖定性狀相關(guān)區(qū)域，但不能確定具體是哪個基因。結(jié)合轉(zhuǎn)錄組和遺傳圖譜，可以鎖定在某個區(qū)域內(nèi)，某些基因或某個基因與性狀相關(guān)，用于后續(xù)驗證。
接下來就跟大家分享一篇高粱遺傳圖譜和轉(zhuǎn)錄組聯(lián)合分析的文章。
研究意義
材料與方法
作圖材料：CK60（氮敏感）和San Chi San（耐氮）雜交得到RIL群體，208個單株。
RNA-seq：雙親及2個混池（RIL中高、低NUE池，分別混5株），處理后3周時取根部
表型評價：2年2點重復(fù)試驗
田間條件：分別將雙親及RIL群體種植于低N田（LN）和正常田（NN）；低N田供給的是N含量為0 kg. ha?1的合成肥，且用燕麥和玉米輪作耗盡其中的N素；正常田的供給的是N含量為100 kg. ha?1的無水氨肥，與大豆輪作補充N素。
表型鑒定：葉綠素含量、形態(tài)和產(chǎn)量性狀，共11個性狀
利用便攜式葉綠素儀測量營養(yǎng)期、開花期時第三片葉以及成熟期植株頂部向下6片葉的葉綠素含量（Chl1、Chl2和Chl3），均進行3次生物學(xué)重復(fù)；株高（PH）：在生理成熟時從植株基部到穗尖的測量；花期（AD）：統(tǒng)計50%的植株開花的天數(shù)；穗頭水分含量（MC2）：收獲時測量鮮重，熱風(fēng)機10-14d后測量干重，3株單穗鮮重和干重的差異（%）為MC2；秸稈含水率（MC1）：方法穗部測量；生物量（BY）：3株地上組織干重的均值；產(chǎn)量（GY）：3株脫粒后粒重的均值；千粒重（TGW）；粒稈比（GS，%）：產(chǎn)量/生物量。
方法：Mapmaker/EXP 3.0和?IciMapping用于遺傳圖譜構(gòu)建，WinQTLcart2.5用于QTL定位
主要結(jié)果
可靠QTL位點：9個主效QTL與本實驗室之前CK60/China17群體的的定位結(jié)果重疊，其中qAD-9在2個環(huán)境（LN和NN）下的定位結(jié)果一致。比較本實驗多個環(huán)境下定位結(jié)果，找到的2個穩(wěn)定QTL，分別為qAD-9和qTGW-3，R2分別為9.2%和15%。這些均側(cè)面驗證了哪些是可靠QTL位點。
RNA-seq結(jié)果：雙親RNA-seq找到486個DEG，兩個RIL池的RNA-seq找到131個DEG，兩者取交集得到54個差異基因，認(rèn)為這些基因可能與耐氮性狀相關(guān)的基因。
QTL定位與RNA-seq所用的參考基因組相同，從而兩者結(jié)果聯(lián)合得到2個基因（此處，作者寫的很含蓄~~~）。這篇文章就寫到了這里，但是我們可以看出后續(xù)作者會對2個主效QTL及這2個基因進行進一步驗證。
圖3、QTL定位結(jié)果
以上2篇文章都通過聯(lián)合分析找到了個位數(shù)的基因，并沒有做后續(xù)的驗證就發(fā)文章了，是不是很心動~~~
你也許會有疑問，怎么不見BSA的身影，因為篇幅的限制，這里并沒有寫，但是遺傳圖+BSA，BSA+GWAS和BSA+轉(zhuǎn)錄組這種思路當(dāng)然也是可以的，但需要注意的是BSA最適合的是質(zhì)量性狀（數(shù)量性狀的主效位點呢，可以是可以，但是表型調(diào)查要細心，誰讓數(shù)量性狀就是這么不讓人省心呢），而遺傳圖譜就不用注意這一點了，無論質(zhì)量性狀還是數(shù)量性狀，它都能定位。
綜上給大家介紹的是聯(lián)合分析的策略來縮小區(qū)間，那是否還有其他的方法呢？當(dāng)然有，想了解嘛？敬請期待下期…….
如果您的項目有任何問題，歡迎點擊下方按鈕咨詢我們，我們將免費為您設(shè)計文章方案。
京公網(wǎng)安備 11011302003368號