2021年4月，中國熱帶農業(yè)科學(xué)院李開(kāi)綿研究員團隊牽頭完成了高質(zhì)量單倍型木薯基因組，相關(guān)研究成果“Allele-defined Genome Reveals Biallelic Differentiation during Cassava Evolution”發(fā)表于國際期刊Molecular Plant（IF：13.162）。近日，該團隊再次牽頭完成了木薯群體研究，該研究以“Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity”為題見(jiàn)刊Genome Biology（IF：13.584）。百邁客有幸作為共同作者參與這兩項研究。接下來(lái)，隨小編一起來(lái)看看兩篇文章如何克服物種高雜合的困難，并以此為突破點(diǎn)開(kāi)展系列研究。
木薯基因組
文章亮點(diǎn)
本研究首先通過(guò)PacBio平臺構建具有“冗余”特征的基因組，再進(jìn)一步結合Hi-C互作信號，將二倍體兩套的基因組進(jìn)行拆分，最終獲得具有18對同源染色體的雙單倍型基因組。前期研究結合本次同源對基因組的差異比較，進(jìn)一步說(shuō)明了同源染色體之間不同遺傳位點(diǎn)的組合對生物表型有重要影響。該研究為廣大科研工作者提供了一個(gè)重要的研究思路：無(wú)雙親本背景下，構建高雜合物種的單倍型基因組。通過(guò)單倍型的分析，為解析物種重要性狀形成的遺傳機制奠定了重要基礎。
研究背景
通過(guò)種間或種內雜交產(chǎn)生雜合子基因組有助于維持植物多樣性并成為新物種的潛在來(lái)源。由于缺乏高質(zhì)量的參考基因組和有限的基因組資源，木薯基因組學(xué)和育種研究相對滯后。木薯（Manihot esculenta Crantz；2n=36）是非洲和其他熱帶地區的主要作物，其具有很高的雜合性，是研究雙等位基因分化的理想系統。
材料方法
基因組：木薯栽培品種SC205；~98.5X PacBio；~35X Hi-C
轉錄組：種下后100-340天（每隔40天取樣）的塊根；以及根莖葉和塊根組織
研究結果
木薯基因組組裝
作者通過(guò)PacBio和Hi-C技術(shù)進(jìn)行基因組的初始組裝，組裝大小為1.1Gb（流式預估單套大小770.3 Mb），表明在基因組的雜合區域組裝了兩個(gè)單倍型。作者通過(guò)自對z確定了冗余區域，并將非冗余木薯基因組作為參考，稱(chēng)為SC205-ref。組裝基因組的連續性指標ContigN50=1.1 Mb，Scaffold N50=34.5 Mb，并經(jīng)過(guò)BUSCO（92.8%）評估、LAI評估（18.46）、高密度遺傳圖譜（相關(guān)系數ρ>0.99）等多種評估方式，表明了本次構建了一個(gè)參考級的高質(zhì)量基因組。
進(jìn)一步使用ALLHiC算法進(jìn)行等位基因定義的基因組組裝，組裝1.52Gb木薯基因組，構建了18對同源染色體高質(zhì)量單體型基因組圖譜，鑒定了24128個(gè)雙等位基因（Bialleles）。并通過(guò)Hi-C信號、同源染色體間共線(xiàn)性比對、遺傳圖譜比對、二代/三代Reads回比、同源對BUSCO評估、雙等位基因鑒定等方式對本次單倍型基因組進(jìn)行驗證。
圖1 木薯雙單倍型基因組特征
不對稱(chēng)演化
通過(guò)同源染色體比對，鑒定到3,107,171 個(gè)SNPs、 500,018 個(gè)indels、57,643 個(gè)PAVs以及632 個(gè)倒位。高雜合度區域的進(jìn)化速度比其他區域更快。
蛋白質(zhì)編碼基因、不同等位基因和顯性表達的雙等位基因的密度，發(fā)現它們都傾向于富集在高度雜合的區域，通常位于染色體的一端或兩端附近。通過(guò)對不同種植時(shí)間和組織部位的轉錄組分析發(fā)現除了第1、15和17對染色體外，幾乎所有同源染色體對中，雙等位基因間廣泛發(fā)生表達不平衡，并且淀粉和蔗糖代謝途徑上的雙等位基因表現出顯著(zhù)的表達分化；
等位基因進(jìn)化
顯性表達的等位基因的啟動(dòng)子和CDS的序列差異顯著(zhù)大于同等表達的雙等位基因。根據非同義和同義核苷酸替換的比率（Ka/Ks）將等位基因對分為7個(gè)亞群。等位基因組進(jìn)化分析揭示了基因組快速進(jìn)化可能是木薯高雜合性形成的重要驅動(dòng)力，且基因組方向性選擇驅動(dòng)等位基因表達分化。
圖2 木薯等位基因組表達與進(jìn)化分析
總結
本研究通過(guò)構建高質(zhì)量的木薯單倍型參考基因組，深化了對木薯雙等位基因變異遺傳基礎的理解，有助于探索雙等位基因的分化和表達優(yōu)勢及其潛在的進(jìn)化驅動(dòng)力，為木薯改良的戰略性開(kāi)發(fā)提供了可能。這些基因組資源將為創(chuàng )新木薯和其他高雜合作物的育種策略提供新的見(jiàn)解。
木薯群體進(jìn)化
文章亮點(diǎn)
本研究突破前期純合子等位基因變異與物種關(guān)鍵農藝性狀的相關(guān)性研究，進(jìn)一步聚焦高雜合物種的等位基因變異。發(fā)現不同基因位點(diǎn)的純合與雜合變異涉及到塊根重量、淀粉含量、抗病性等多種木薯關(guān)鍵的農藝性狀。闡明了雜合位點(diǎn)的
保持對于優(yōu)異性狀的形成具有重要貢獻，為高雜合物種研究提供了一個(gè)新穎的研究思路。
研究背景
許多野生物種通過(guò)人工選擇馴化為栽培作物，隨著(zhù)人類(lèi)遷移、定向選擇和進(jìn)一步的性狀改良以滿(mǎn)足人類(lèi)需求，栽培作物在性狀上獲得了區域特異性差異。利用基因組測序技術(shù)，我們可以追蹤作物馴化和育種歷史，從而更好地了解人類(lèi)選擇如何形成作物基因組。群體中基因型和表型之間的關(guān)聯(lián)顯示純合子等位基因變異與許多作物的關(guān)鍵農藝性狀顯著(zhù)相關(guān)，從而加速了育種過(guò)程。然而，許多作物是高度雜合的，雜合性變異對農藝性狀的影響尚不清楚。
栽培木薯是6000多年前從亞馬遜流域的野生祖先Manihot esculenta ssp. flabellifolia馴化而來(lái)。木薯的馴化產(chǎn)生了包括年生長(cháng)周期、高初始生長(cháng)率、產(chǎn)量增加和高淀粉含量在內的特征，這使得木薯成為人類(lèi)消費和工業(yè)生物燃料應用的理想能源。但由于基因組資源有限，木薯基因組學(xué)和育種研究相對滯后。群體基因組資源的開(kāi)發(fā)將加速這一重要作物的遺傳改良進(jìn)程。
材料方法
代表來(lái)自全球15個(gè)國家主要木薯生產(chǎn)區的388份木薯種質(zhì)（14份野生、38份地方品種和336份栽培品種）；51份材料來(lái)自發(fā)表重測序數據，337份在本研究中測序?？偣伯a(chǎn)生3.35 Tb數據，平均測序深度為9.45×。對本次測序的337份木薯，進(jìn)行三年的33個(gè)形態(tài)性狀GWAS分析。
研究結果
基因組變異圖譜與系統發(fā)育關(guān)系
通過(guò)與參考基因組比對，共鑒定1,344,463個(gè)SNPs以及1,018,832個(gè)InDels。系統發(fā)育和主成分分析（PCA）將388份木薯材料分為兩大類(lèi)（圖1c，d），第一組包括13個(gè)野生祖先，第二組包含所有地方品種、栽培品種以及自巴西的野生祖先FLA433-2（具類(lèi)似木薯的塊根）。結果支持了FLA433-2是木薯野生祖先種的假說(shuō)。系統進(jìn)化及FST、p的分析結果與考古證據相符，表明南美和非洲木薯種群之間以及非洲和亞洲種群之間存在著(zhù)密切的系統發(fā)育關(guān)系，提出了木薯從南美洲到非洲再到亞洲的傳播馴化路徑。
圖1 木薯全基因組變異圖譜及其親緣關(guān)系
全基因組關(guān)聯(lián)分析（GWAS）
GWAS分析鑒定出與23個(gè)重要農藝性狀關(guān)聯(lián)的52遺傳標記。表皮類(lèi)型（光滑或粗糙）和疤痕嚴重影響木薯塊根的外觀(guān)質(zhì)量。10號染色體上Sc10g01040中的三個(gè)和Sc10g02050中的七個(gè)非同義SNP與塊根表皮型相關(guān)。在10號染色體Sc10g01040、Sc10g01050以及3號染色體Sc03g001750、13號染色體Sc13g000920中攜帶雜合等位基因和雜合單倍型的木薯材料比攜帶純合等位基因和純合單倍型的木薯材料具有更高的塊根平滑表皮頻率。木薯地上部分的重量反映了植物的生長(cháng)活力，并成為其用作飼料的一個(gè)重要因素；塊根支鏈淀粉含量決定木薯的工業(yè)價(jià)值。5號染色體上Sc05g013530、11號染色體上Sc11g000910攜帶雜合子等位基因的材料表現出顯著(zhù)更高的地上重量和支鏈淀粉含量。2號染色體Sc02g008280攜帶AT或TT等位基因的材料具更高紅色內皮顏色頻率，18號染色體Sc18g013220攜帶雜合子AC等位基因的材料比攜帶純合子CC或AA等位基因的材料顯示出顯著(zhù)更低的螨感染指數。對這些候選基因的表達水平進(jìn)行檢測，發(fā)現都在塊根發(fā)育的不同組織或階段大量表達。種種結果表明，雜合性等位基因變異在關(guān)鍵農藝性狀中的潛在作用，其中每個(gè)候選基因中攜帶雜合等位基因的木薯材料比攜帶相應純合等位基因的木薯材料具有更理想的表型。
圖2 GWAS鑒定與關(guān)鍵農藝性狀相關(guān)的雜合性變異候選基因
作者鑒定了79個(gè)與塊根中淀粉積累相關(guān)的高度雜合子片段，其中20個(gè)雜合塊與6個(gè)GWAS信號重疊。此外，在1號染色體上攜帶高度雜合子片段的木薯材料比攜帶相應純合子片段的木薯材料表現出更高的塊根紅色內皮細胞頻率。這些結果進(jìn)一步支持雜合性變異與關(guān)鍵農藝性狀相關(guān)的假設。
木薯雜合性的選擇特征
為了追蹤雜合性的可能進(jìn)化場(chǎng)景，作者根據參考基因組的雜合SNP將木薯基因組分為雜合區和純合區。然后比較了374個(gè)品種的雜合區和純合區的基因組多樣性。在大多數（18條染色體中的16條）染色體中，雜合區具有顯著(zhù)高于純合區的遺傳多樣性（π值）和正Tajima’s D值，表明平衡選擇有助于維持基因組雜合性。在品種和野生祖先之間的大多數染色體（18個(gè)染色體中的14個(gè)）中，雜合區的FST值低于純合區的FST值也支持了這一發(fā)現。檢測到81個(gè)雜合度和核苷酸多樣性降低的選擇性位點(diǎn)，包含548個(gè)基因，這些基因在多種生物過(guò)程中富集，包括生長(cháng)、發(fā)育、激素代謝和反應以及免疫相關(guān)過(guò)程。
圖3 木薯雜合性的選擇特征
馴化使木薯塊根的產(chǎn)量和淀粉含量顯著(zhù)增加。栽培品種的塊根重和淀粉含量分別比野生后代高約60倍和3倍（圖4a，b）。選擇性清除分析結合轉錄組分析發(fā)現2號染色體MeT
R1周?chē)s合性降低的選擇信號，并且MeTIR1在攜帶GG等位基因的材料中的表達顯著(zhù)高于攜帶CG和CC等位基因的材料（圖4i），說(shuō)明MeTIR1純合變異促使塊根淀粉含量提升。并通過(guò)瞬時(shí)過(guò)表達和沉默技術(shù)進(jìn)一步驗證了MeTIR1在木薯中的功能。進(jìn)一步結合關(guān)聯(lián)分析發(fā)現，對2號染色體攜帶MeAHL17 AA等位基因的材料，經(jīng)選擇馴化促使塊根產(chǎn)量提升的同時(shí)導致木薯細菌性枯萎病抗性丟失。
圖4 木薯高淀粉含量塊根馴化相關(guān)雜合性降低的選擇作用
圖5 木薯細菌性枯萎病候選基因MeAHL17的鑒定
MeAHL17和MeTIR1位于2號染色體上相同的選擇性清除區內。連鎖不平衡分析發(fā)現，兩個(gè)基因緊密連鎖，攜帶AG（MeAHL17）和GG（MeTIR1）等位基因組合的品種顯示出較高的淀粉含量和對木薯細菌性枯萎病的高抗性頻率。
圖6 MeAHL17和MeTIR1等位基因組合對性狀變異的影響
總結
本研究提供了388份木薯材料的變異圖譜，鑒定了23個(gè)農藝性狀的52個(gè)位點(diǎn)，揭示了與關(guān)鍵農藝性狀和木薯馴化相關(guān)的雜合性等位基因變異。檢測到81個(gè)雜合度和核苷酸多樣性降低的選擇性位點(diǎn)，相關(guān)基因在多種生物過(guò)程中富集，包括生長(cháng)、發(fā)育、激素代謝和反應以及免疫相關(guān)過(guò)程。人工選擇MeTIR1和MeAHL17中的純合等位基因有助于大淀粉塊根的馴化。而在MeAHL17中選擇純合子等位基因與增加塊根重和木薯細菌性枯萎病易感性相關(guān)。本研究將有助于闡明與關(guān)鍵農藝性狀和木薯馴化相關(guān)的雜合性變異的遺傳基礎，并對木薯馴化過(guò)程中雜合性的變異提供了見(jiàn)解，為木薯和其他高雜合物種育種改良的戰略發(fā)展提供依據。
寫(xiě)在最后：
大多數二倍體基因組組裝都忽略了同源染色體之間的差異，將基因組組裝成一個(gè)假的單倍體序列（嵌合型）。但隨著(zhù)研究的深入，發(fā)現僅單套的基因組數據難以完全演示該物種的全面信息，例如在木薯、馬鈴薯、茶樹(shù)等物種中已有相關(guān)研究（Zhou Q et al., Nature Genetics.2020;Zhang X et al., Nature Genetics.2021; Wang P et al.,Horticulture Research）。同源多倍體中問(wèn)題尤其明顯，同源染色體之間不同遺傳位點(diǎn)的組合對生物表型有重要影響，如動(dòng)植物中的雜種優(yōu)勢、某些物種雜交不育現象等。
單倍型等位基因間差異對基因表達、功能及其表型都有著(zhù)重要影響，大多數雜交品種的優(yōu)勢表型都受等位基因調控。通過(guò)單倍型的分析，可為解析物種重要性狀形成的遺傳機制奠定重要基礎。
無(wú)論是單倍型基因組的構建還是泛基因組的構建，單個(gè)的樣本往往不具有性狀分析的普適性，需要進(jìn)一步通過(guò)群體規模的大樣本檢測，挖掘與自然人工馴化、性狀功能關(guān)聯(lián)的關(guān)鍵位點(diǎn)，為物種的育種改良與相關(guān)基因研究奠定重要理論基礎。
百邁客自2009年成立以來(lái)深耕于群體遺傳研究，同時(shí)具有近10年基因組分析經(jīng)驗?，F已在基因組、Hi-C、遺傳圖譜等技術(shù)上擁有實(shí)驗+生信分析雙保障。并與國內外70余所科研單位在國際期刊合作發(fā)表500余篇高質(zhì)量文章，累計影響因子3000+。
為了答謝大家長(cháng)期以來(lái)的信任與支持，百邁客特在年末之際，推出多項優(yōu)惠活動(dòng)，助力科研！ 詳情請咨詢(xún)當地銷(xiāo)售經(jīng)理或撥打電話(huà)：400-600-3186
京公網(wǎng)安備 11011302003368號