提到性狀定位，我們通常會(huì)想到大群體、復(fù)雜算法、耗資巨大，例如全基因組關(guān)聯(lián)分析。今天向大家介紹的 BSA 則不同，這種方法不僅材料容易準(zhǔn)備、算法原理清晰明了，更重要的是不需要花費(fèi)高額經(jīng)費(fèi)，十分適合性狀初定位。
 

BSA 混池分組分析法簡(jiǎn)介

BSA（bulked segregant analysis）混池分組分析法，根據(jù)目標(biāo)性狀表型對(duì)分離群體中的個(gè)體分別進(jìn)行分組混合，依據(jù)目標(biāo)性狀的相對(duì)差異選擇表型極端的個(gè)體分成兩組，然后將兩組的個(gè)體或株系 DNA 混合，形成相對(duì)的 DNA 混合池。通過(guò)在親本和子代混池之間的多態(tài)性標(biāo)記篩選即可完成對(duì)目標(biāo)性狀的定位。

BSA 適用群體

關(guān)于性狀定位常用的群體類(lèi)型有兩種，一種是自然群體（個(gè)體間遺傳背景差異較大，通常來(lái)自不同地域、不同品種、甚至于不同進(jìn)化程度的材料），另一種是家系群體（群體遺傳背景較為相似，通常全部來(lái)自于兩個(gè)親本個(gè)體，在目標(biāo)性狀上有明顯表型分離）。

BSA 采用的是家系群體，常用的家系群體有 F2、BC（n）、RIL、DH 等等。對(duì)于家系群體的構(gòu)建，親本選擇建議：同一物種前提下，目標(biāo)性狀盡可能差異大的兩個(gè)個(gè)體。

BSA 原理及分析邏輯
BSA 中通常需要測(cè)序 4 個(gè)樣本：兩個(gè)親本+兩個(gè)子代混池

測(cè)序拿到原始數(shù)據(jù)后，可通過(guò)質(zhì)控、與參考基因組比對(duì)、SNP calling 等步驟得到各樣本 SNP 標(biāo)記分型信息。對(duì)于重測(cè)序而言，無(wú)論是 BSA 還是其他性狀定位方法，SNP 分型信息是后續(xù)分析的基礎(chǔ)。

接下來(lái)如何從海量的 SNP 標(biāo)記中篩選到與性狀性格的那一部分，進(jìn)而定位到性狀相關(guān)基因是關(guān)鍵。對(duì)于“雙親健全”（因?yàn)閷?shí)際情況可能只有兩個(gè)子代混池）的 BSA 項(xiàng)目來(lái)說(shuō)，是一個(gè)三步走策略：

Step 1. 篩選兩親本間純合且差異的標(biāo)記

家系群體構(gòu)建時(shí)我們知道，首先選擇在目標(biāo)性狀上差異明顯的個(gè)體為親本，表型有差異，跟表型相關(guān)的 SNP 標(biāo)記就一定有差異。在親本之間無(wú)差異的SNP可以先篩掉，剩余的親本間純合且差異的這部分再借助子代混池進(jìn)一步聚焦。
 

Step 2. 子代混池篩選SNP標(biāo)記并獲得候選區(qū)間

標(biāo)記與性狀進(jìn)行關(guān)聯(lián)是性狀定位最核心的部分，對(duì)于 BSA 而言，有兩種關(guān)聯(lián)方法，一種是 ΔSNP-Index，一種是歐式距離（Euclidean distance, ED）。

ΔSNP-Index

SNP-Index可以理解成突變位點(diǎn)遺傳物質(zhì)來(lái)源于某一個(gè)親本的頻率，對(duì)于混池的某一個(gè)位點(diǎn)來(lái)說(shuō)，SNP-Index即為該位點(diǎn)與某親本序列相同的reads占總reads的比例。而Δ SNP-Index即為2個(gè)混池之間SNP-Index的差值�；�窗方式計(jì)算各區(qū)域的Δ SNP-Index，通過(guò)置信區(qū)間的設(shè)置可以判斷區(qū)域是否與目的性狀存在關(guān)聯(lián)。

歐式距離（Euclidean distance, ED）

歐式距離計(jì)算將混池中每一個(gè)位點(diǎn)的突變信息抽象成坐標(biāo)上的一個(gè)點(diǎn)，并計(jì)算混池之間每個(gè)相同位點(diǎn)之間的歐氏距離：

ED²=(A1-A2)²+(T1-T2)²+(C1-C2)²+(G1-G2)²

An、Tn、Cn、Gn 分別指 ATCG 四種堿基在混池 n 中的頻率

取所有位點(diǎn)擬合值的 median+3SD 作為分析的關(guān)聯(lián)閾值，根據(jù)關(guān)聯(lián)閾值判定候選區(qū)間。

Step 3. 區(qū)間內(nèi)候選基因注釋及功能富集

篩選 SNP 位點(diǎn)不是目的，關(guān)于結(jié)果，候選區(qū)域更具有生物學(xué)意義。確定候選區(qū)間后，注釋區(qū)間內(nèi)所包含的基因，并對(duì)候選基因集做 KEGG、GO 富集分析。

BSA方案設(shè)計(jì)
說(shuō)到方案，主要涉及三個(gè)方面：親本測(cè)序深度、子代混池樣本數(shù)、子代混池測(cè)序深度。方案設(shè)計(jì)是否合理直接關(guān)系到結(jié)果定位的精度和項(xiàng)目費(fèi)用。
 

極端混池的樣本數(shù)量越多，混池測(cè)序深度越高，定位結(jié)果精度越高，但同時(shí)項(xiàng)目費(fèi)用也會(huì)越高。有學(xué)者借助模式植物擬南芥專(zhuān)門(mén)針對(duì)子代混池樣本數(shù)和測(cè)序深度對(duì)定位結(jié)果的影響進(jìn)行了研究。結(jié)果發(fā)現(xiàn)，當(dāng)混池樣本數(shù)量達(dá)到 20 個(gè)，同時(shí)測(cè)試深度在 20X 以上的情況下，能夠得到相對(duì)理想的定位結(jié)果。超過(guò)這個(gè)數(shù)量，隨著成本消耗的增加，定位精度的提升效果有限。綜合考慮，通常推薦方案如下：

注：推薦方案不能適用于所有項(xiàng)目，具體方案要根據(jù)項(xiàng)目具體情況而定

BSA 結(jié)果展示

高通量測(cè)序項(xiàng)目的結(jié)題報(bào)告內(nèi)容向來(lái)比較多，包括各種形式的圖表展示，抓住重點(diǎn)才能讓我們對(duì)結(jié)果的解讀準(zhǔn)確而高效。對(duì)于 BSA 而言，要重點(diǎn)關(guān)注以下幾個(gè)方面：

• 候選區(qū)域

• 候選區(qū)域內(nèi)基因注釋

• 基因功能富集

結(jié)果細(xì)節(jié)

< BSA-demo報(bào)告 >

BSA 常見(jiàn)問(wèn)題

Q1：沒(méi)有參考基因組的物種可以做 BSA 嗎？

BSA 性狀定位的方法適用于有參考基因組的物種。參考基因組組裝質(zhì)量好壞以及注釋是否完全對(duì) BSA 結(jié)果有一定的影響。建議盡量選擇組裝至染色體水平，且基因組注釋較完全的參考基因組。

Q2：混池樣本如何選擇，一定要到 20 個(gè)嗎？

混池的樣本數(shù)通常選取極端性狀（群體總樣本數(shù)的 5～10% 內(nèi)）的個(gè)體進(jìn)行混池�；斐刂腥魳O端性狀個(gè)體數(shù)太少，可能造成樣品數(shù)不夠，不具備代表性；若樣本數(shù)過(guò)高可能會(huì)引入雜合個(gè)體，產(chǎn)生干擾。需要根據(jù)群體的實(shí)際情況，首先保證表型足夠極端，其次考慮樣本數(shù)量。不一定必須到 20。例如，某兔子毛色性狀定位研究，通過(guò)兩親本雜交、F1 個(gè)體混交，得到 35 只 F2 個(gè)體，將極端表型控制在 20% 以?xún)?nèi)，分別選擇 5 只、6 只進(jìn)行混合，測(cè)序深度各為 10X，也得到了較為理想的定位結(jié)果。

Q3：影響定位結(jié)果的因素？

評(píng)估定位結(jié)果主要參考候選區(qū)域大小和候選基因個(gè)數(shù)，影響它們的因素包括但不限于兩親本材料的差異大小、家系群體大小及性狀分離情況、子代混池個(gè)數(shù)、親本及子代混池測(cè)序深度、參考基因組水平、目標(biāo)性狀特點(diǎn)等。

Q4：子代混池構(gòu)建是先提取在混合還是先混合再提�。�

先提 DNA 再等量混合，可以減少系統(tǒng)誤差，近年發(fā)表的多數(shù)文獻(xiàn)都是先提 DNA 再等量混合。

Q5：BSA性狀定位是否可以采用簡(jiǎn)化基因組測(cè)序？

簡(jiǎn)化基因組捕獲的基因組區(qū)域有限，一般僅能捕獲 3%～30%，如果變異的區(qū)域正好沒(méi)有捕獲到就不能找到目標(biāo)性狀的基因了。因此用簡(jiǎn)化基因組的風(fēng)險(xiǎn)很大，不建議做簡(jiǎn)化基因組的BSA性狀定位。

當(dāng)然要做性狀定位研究，除了 BSA，還有遺傳連鎖圖譜、全基因組關(guān)聯(lián)分析等經(jīng)典方案，有任何疑問(wèn)可聯(lián)系技術(shù)專(zhuān)家為您免費(fèi)答疑解惑哦。