2021年3月18日，由河南農(nóng)業(yè)大學王道文研究員團隊聯(lián)合多個單位共同破譯黑麥基因組。相關研究成果“A high-quality genome assembly highlights rye genomic characteristics and agronomically important genes”見刊于國際知名期刊Nature Genetics。河南農(nóng)業(yè)大學李廣偉博士、王立建博士、靳懷冰博士，北京大學現(xiàn)代農(nóng)業(yè)研究院何航研究員，四川農(nóng)業(yè)大學任天恒副教授以及百邁客生物科技有限公司李緒明作為共同一作共同參與本項研究。河南農(nóng)業(yè)大學楊青華教授、楊建平研究員，中國科學院遺傳與發(fā)育生物學研究所張坤普副研究員以及河南農(nóng)業(yè)大學王道文研究員為并列通訊作者。

此次論文的發(fā)表引起了學術界和媒體界的高度關注，百邁客作為合作參與方及測序技術服務企業(yè)，為促進科研成果更廣泛更科學的傳播，本次特邀請到國內(nèi)外兩篇論文團隊一起解析黑麥基因組的奧秘。歡迎各位專家學者報名參加本次國際交流會,會議定于北京時間2021年3月24日（周三）下午15:45—18:00，詳細會議流程請點擊：黑麥基因組學研究成果《Nature Genetics》見刊國際交流會。

（或點擊二維碼直接報名）

?

PART ONE

論文內(nèi)容詳解

研究背景

黑麥（Secale cereale, 2n = 2x = 14, RR）屬于禾本科小麥族黑麥屬，雖然與普通小麥（Triticum aestivum，Ta；2n=6x=42；AABBDD）和大麥（Hordeum vulgare，Hv）有親緣關系，但黑麥具有獨特的農(nóng)藝性狀和基因組特性。黑麥1RS染色體臂攜帶的抗病基因通過遠緣雜交轉入小麥基因組，為小麥生產(chǎn)中白粉病和條銹病的防治做出了巨大貢獻。此外，黑麥也可以與普通小麥遠緣雜交，染色體加倍后，人工合成八倍體小黑麥，具有比黑麥更高的生物量和產(chǎn)量。因此，黑麥是許多國家重要的糧食和飼料作物，也是全球小麥和小黑麥改良的重要遺傳資源。

威寧黑麥是我國的栽培黑麥早抽穗優(yōu)良品種，具有抗白粉病和條銹病的能力。為了解析黑麥優(yōu)良性狀的遺傳和分子基礎，促進黑麥及相關作物的基因組和育種研究，作者對威寧黑麥進行了基因組測序和分析。

?

材料與方法

基因組：

材料：威寧黑麥自交系，自交18代；熒光原位雜交（FISH）確認基因組染色體對=7組；

測序：270 bp小片段文庫（430Gb, ~54X），Illumina 平臺；20 Kb大片段文庫（497Gb,~63X），Pacbio平臺；Hi-C 文庫（560 Gb, ~71X），Illumina平臺；光學圖譜（779.55 Gb，~99X），Bionano Irys平臺；

基因注釋轉錄組測序：正?；蛎{迫（冷或干旱）條件下栽培的植物中的葉、莖、根和穗樣品，以及在開花10、20、30、40天后收獲的發(fā)育中的樣品, 采用Illumina及Pacbio平臺進行RNA-seq及Iso-Seq；

遺傳圖譜QTL: 295份威寧黑麥×荊州黑麥雜交的F2代樣品，采用SLAF-seq進行標記開發(fā)；

抽穗基因表達：威寧黑麥和荊州黑麥樣品，播種后4、7和10天采集葉片樣品，每個時間點使用3個生物重復，使用Illumina平臺進行測序；

選擇進化分析：已發(fā)表的101份家養(yǎng)黑麥和野生黑麥品種材料的公共數(shù)據(jù)，包括81份S. cereala、5份S. vavilovii、11份S. strictum和4份S. sylvestre。

?

研究內(nèi)容

1. 基因組裝與注釋
2. 轉座子分析
3. 黑麥基因組進化和染色體共線性研究
4. 基因重復及其對淀粉生物合成基因的影響分析
5. 黑麥種子貯藏蛋白(SSP)基因位點的解析
6. 轉錄因子(TF)及抗病基因的注釋和分析
7. 早抽穗性狀相關基因表達特征的研究
8. 黑麥馴化過程中染色體區(qū)域和基因座的挖掘

主要研究結果

1.?基因組裝與注釋

流式預估黑麥基因組大小約為7.86 Gb，結合PacBio、Illumina、Hi-C、遺傳圖譜、BioNano光學圖譜等技術進行基因組組裝。最終組裝7.74 Gb基因組（為預估基因組大小的98.47%），scaffold N50 為1.04 Gb，并將93.67%的序列掛載至7條染色體上（圖1）。黑麥中每條染色體基因組大小均在1G左右（2R、3R、4R、6R、7R（~1Gb）；1R（0.94097 Gb）、5R（0.99891 Gb) ），比烏拉爾圖小麥（T. urartu；Tu；AA型）、節(jié)節(jié)麥（Aegilops tauschii；Aet；DD型）、野生二粒小麥（T. turgidum?ssp. dicoccoides；WEW；AABB型）、普通小麥（Ta）和大麥（Hv）等復雜的麥類中的單條染色體基因組均要大。超大的基因組及染色體，對黑麥基因組組裝，特別是染色體掛載帶來了巨大的挑戰(zhàn)。

將組裝結果與兩個冬季黑麥品種（Lo7和Lo225）構建的染色體連鎖圖相比，威寧黑麥1R至7R物理圖具有較高的一致性。在先前報道的Lo7的pyro-sequencing reads中97.45%可以定位到威寧基因組，平均序列同源性為97.71%，平均序列覆蓋率為97.27%。LAI值為18.42，遠高于先前發(fā)表的小麥和大麥基因組的LAI值。BUSCO評估結果為96.74%。并注釋了86,991個蛋白編碼基因。盡管黑麥基因組非常復雜，通過以上評估結果說明，本次研究構建了一個高質量的威寧黑麥基因組。

圖1 黑麥基因組特征

2.?轉座子分析

威寧黑麥基因組中90.31%被注釋為轉座子（TE），共包含537個家族的2,671,941個成員，這些TE的含量明顯比普通小麥 (84.70%), 烏拉爾圖小麥 (81.42%), 節(jié)節(jié)麥 (84.40%), 野生二粒小麥 (82.20%)以及大麥 (80.80%)更高。其中長末端重復反轉錄轉座子（LTR-RTs）是主要的轉座子，在注釋的TEs中占 84.49%。

與烏拉爾圖小麥、節(jié)節(jié)麥及大麥的LTR-RTs進行比較發(fā)現(xiàn)：（1）Gypsy是威寧黑麥基因組擴張的主要原因之一(Fig. 2a)，并且有3個LTR-RT家族（Daniela, Sumaya, Sumana）在威寧黑麥中特異擴張，其中Daniela的占比最高 (Fig. 2b)。（2）威寧黑麥中完整的LTR-RTs插入時間存在獨特的雙峰分布(Fig. 2c)，最近一個擴增峰出現(xiàn)在50萬年前，另出現(xiàn)在1.7 百萬年（MYA）左右（與大麥相同） (Fig. 2c)。（3）這種雙峰分布模式由Gypsy?RTs的擴增主導(Fig. 2d)。

圖2 黑麥中轉座子分析

3.?黑麥基因組進化和染色體共線性研究

通過比較威寧黑麥、烏拉爾圖小麥、節(jié)節(jié)麥、普通小麥（亞基因組TaA、TaB和TaD）、大麥、水稻Os（Oryza sativa ssp.?japonica）、二穗短柄草Bd（Brachypodium distachyon）、玉米Zm（Zea mays）、高粱Sb（Sorghum bicolor）、谷子Si（Setaria italica）基因組，共找到2517個單拷貝同源基因。通過對單拷貝基因組構建進化樹和計算分化時間發(fā)現(xiàn)，在大麥和小麥分化后（15MYA）黑麥和二倍體小麥發(fā)生了分化（9.6 MYA） (Fig. 3a)。

作者以水稻為祖先參考基因組，研究了威寧黑麥的染色體進化。威寧黑麥與水稻共鑒定出23個大的共線性區(qū)，包含10949對同源基因，可推斷出祖先染色體片段在1R到7R之間的排列（圖3b）：（1）3R來源于一條古老的染色體AGK1/Os1，該染色體的一段易位到6RL；（2）1R和2R分別由兩條祖先染色體組成，1R與AGK10/Os10嵌套插入AGK5/Os5有關，2R與AGK7/Os7嵌套插入AGK4/Os4有關；（3）4R, 5R, 6R,7R是通過復雜易位從至少三條祖先染色體上獲得的（圖3b）。

在威寧黑麥基因組與普通小麥3個亞基因組的比較中，1R、2R和3R分別與小麥1、2和3組染色體完全共線。在4R中發(fā)現(xiàn)與4A/4B/4D、7A/7B/7D或6A/6B/6D部分共線性的三個區(qū)域。5R與5A完全共線，與5B、5D部分共線是由于易位的4B或4D片段在5BL或5DL的長臂融合（圖3c）。在6R中，觀察到3個區(qū)域與6A/6B/6D、3A/3B/3D或7A/7B/7D部分共線。這些數(shù)據(jù)將有助于黑麥在禾本科比較基因組學研究以及黑麥與普通小麥雜交研究中的應用。

圖3 黑麥基因組進化及染色體共線性分析

4.?基因重復及其對淀粉生物合成基因的影響分析

作者在威寧黑麥中檢測到4217個單拷貝基因、23753個分散重復基因 (DDGs)、6659個近端重復基因 (PDGs)、7077個串聯(lián)重復基因(TDGs)和1866個片段重復基因。轉座重復基因（TrDGs）由TE活性誘導的，它們是DDGs的主要組成部分。作者以大麥為參考，在威寧黑麥中鑒定出10357個TrDG，遠遠大于以相同方式計算的烏拉爾圖小麥（7145）和節(jié)節(jié)麥（7351）中TrDG的數(shù)量（圖4b）。威寧黑麥所特有的TrDG（5926）也比烏拉爾圖小麥（3513）和節(jié)節(jié)麥（3327）所特有的TrDG更多（圖4b）。

接下來作者研究了黑麥淀粉生物合成相關基因（SBRGs）中的基因復制。研究發(fā)現(xiàn)這些重復類型在黑麥淀粉生物合成相關基因（SBRGs）中普遍存在，并且不同的SBRG的復制基因之間往往表現(xiàn)出表達差異，說明不同類型的基因復制可以豐富黑麥基因在重要生物過程中的遺傳多樣性（圖4c），這些黑麥SBRGs的新變化可能為調(diào)控植物淀粉生物合成和性質提供新的酶活性，因此解析全套黑麥SBRGs有利于提高黑麥的產(chǎn)量潛力和營養(yǎng)品質。

?圖4 黑麥基因重復及其對淀粉合成相關基因多樣性的影響

5.?黑麥種子貯藏蛋白(SSP)基因位點的分離

與小麥和大麥相似，黑麥在胚乳組織中積累了豐富的儲存蛋白SSPs。作者利用威寧基因組組裝技術對黑麥SSP基因座進行了分析。

在威寧黑麥中未發(fā)現(xiàn)小麥低分子量麥谷蛋白亞基（LMW-GSs）或大麥B-hordein的同源序列（圖5b），表明在黑麥進化過程中攜帶這些基因的染色體片段缺失，這可能是1BL1RS易位系小麥品種中，品質受影響的主要原因。并且在威寧黑麥基因組中未發(fā)現(xiàn)α-醇溶蛋白基因，這說明小麥及其近緣種的α-醇溶蛋白（α-gliadin）基因可能在小麥和黑麥分化之后進化產(chǎn)生的。這些SSP分析結果闡明了黑麥堿基因座的結構和組成，這將有助于進一步研究黑麥、小黑麥和小麥的加工和營養(yǎng)品質。

圖5 黑麥堿位點分析

6.?轉錄因子(TF)及抗病基因的注釋和分析

作者利用iTAK預測了威寧黑麥和其他8種禾本科植物的TF基因，在注釋的65個轉錄因子基因家族中，威寧黑麥有28個家族成員增加，其中AP2/ERF TF基因家族成員增加幅度較大。威寧黑麥的抗病相關基因（DRA）數(shù)量（1989）比烏拉爾圖小麥（1621）、節(jié)節(jié)麥（1758）、大麥（1508）、二穗短柄草（1178）、水稻（1575）以及普通小麥的A（1836）、B（1728）和D（1888）亞基因組多。鑒于AP2/ERF TFs和DRA基因在植物對非生物逆境和生物逆境的反應中的重要作用，上述發(fā)現(xiàn)可能有助于黑麥和相關作物的有效遺傳研究和分子改良。

 

7.?早抽穗性狀相關基因表達特征的研究

在長日照條件下，威寧黑麥比荊州黑麥提前抽穗10-12天（圖6a），這與威寧黑麥莖尖分生組織發(fā)育更快有關（圖6b）。在威寧黑麥基因組中，注釋到在長日照條件下高表達的兩個開花位點（FT）基因ScFT1（ScWN4R01G446100）和ScFT2（ScWN3R01G192500）。播種后7天和10天，ScFT1和ScFT2在威寧黑麥中的表達水平顯著高于荊州黑麥（圖6c），且ScFT蛋白在黑麥中積累到相對較高水平，而荊州黑麥中幾乎沒有（圖6d）。檢測到的ScFT蛋白的大小（~29 kDa）比預測出的ScFT1和ScFT2的分子量大（~19 kDa）（圖6d），表明ScFT蛋白具有潛在的翻譯后修飾。用高效檢測磷蛋白的磷酸標記SDS-PAGE分析表明，威寧黑麥中ScFT確實發(fā)生了翻譯后磷酸化修飾。

作者突變了ScFT2磷酸化相關的兩個殘基（S76和T132），并為ScFT2構建了一系列去磷模擬位點（S76A、T132A和S76A/T132A）和磷模擬位點（S76D、T132D和S76D/T132D）。當使用馬鈴薯X病毒載體在煙草中進行外源表達時，ScFT2和去磷雙突變體ScFT2^S76A/T132A相對于游離GFP（對照）和其他ScFT2突變體，表現(xiàn)出持續(xù)促進煙草生長（圖6e）。與GFP相比，ScFT2和三個去磷突變體（ScFT2^S76A、ScFT2^T132A及ScFT2^S76A/T132A）的異位表達提高了開花植株的百分比，ScFT2^S76A/T132A尤其明顯，但在表達三個擬磷突變體（ScFT2^S76D, ScFT2^T132D, or ScFT2^S76D/T132D）中沒有觀察到這種促進作用（圖6f）。免疫印跡分析表明，ScFT2、ScFT2^S76A、ScFT2^T132A和ScFT2^S76A/T132A在煙草植株中的積累量相當高，但ScFT2^S76D、ScFT2^T132D和ScFT2^S76D/T132D在煙草植株中的積累量卻很低（圖6g）。因此，保守的S76和T132殘基的改變影響了ScFT2控制植物開花的功能，這與ScFT2蛋白穩(wěn)定性的改變有關。本研究首次發(fā)現(xiàn)FT磷酸化對開花時間控制的影響，為更全面地探索FT蛋白控制植物開花的分子和生化機制提供了新的途徑。

作者進一步研究了光周期Photoperiod（Ppd）基因的表達，該基因在長日照條件下正調(diào)控FT的表達。在威寧和荊州黑麥的轉錄組分析中發(fā)現(xiàn)了一個表達Ppd的基因ScPpd1（ScWN2R01G043000）。該基因在威寧黑麥內(nèi)的表達非常早，在播種2 天后達到表達高峰；而荊州黑麥在播種4天后才達到高峰（圖6h）。根據(jù)ScPpd1對黑麥抽穗期的調(diào)控作用，研究者利用威寧×荊州F2代群體，檢測到與前期研究一致的三個主效抽穗期QTL（Hd2R、Hd5R和Hd6R）。

?圖6 威寧黑麥早期抽穗性狀的發(fā)育及基因表達特征

8.?黑麥馴化過程中染色體區(qū)域和基因座的挖掘。

對馴化基因的分析可以促進對作物性狀的理解和改良，但在黑麥中這類基因的分子分析方面進展甚微。作者通過全基因組選擇清除分析，利用在栽培黑麥和瓦維洛夫黑麥（S. vavilovii）之間鑒定的123647個SNPs，挖掘黑麥馴化相關染色體區(qū)域和基因座。DRI、F_ST、XP-CLR分析中，共同識別到11個選擇信號（圖7a-c）。通過與水稻和大麥的共線性比較，發(fā)現(xiàn)了一些可能的選擇清除位點，包括已在水稻或大麥中已被功能分析的ScBC1、ScBtr、ScGW2、ScMOC1、ScID1和ScWx的同源基因（圖7a-c）。

檢測到的ScID1基因座包含一對具有相同編碼序列的ScID1同源序列（ScWN6R01G057200和ScWN6R01G057300，下稱ScID1.1和ScID1.2）（圖7d）。ScID1.1和ScID1.2蛋白與玉米ID1（63.19%）和水稻RID1（65.34%）具有很強的同源性，這兩個蛋白在玉米和水稻中都被發(fā)現(xiàn)調(diào)控著從營養(yǎng)體向花發(fā)育的轉換。ScID1.1和ScID1.2在威寧黑麥幼葉中的表達水平高于荊州黑麥（圖7e）。在威寧×荊州分離的F2群體中，ScID1^JZ/JZ純合植株的平均抽穗期顯著晚于ScID1^JZ/WN或ScID1^WN/WN個體（圖7f），這與荊州黑麥相對于威寧黑麥的晚花表型相一致（圖6a）。以上結果表明ScID1可能參與了抽穗期的調(diào)控，并可能通過黑麥馴化進行選擇，使作物成熟度得到適當?shù)恼{(diào)整，以更好地適應生長環(huán)境。

?圖7 黑麥馴化相關染色體區(qū)域和位點的鑒定與分析

?

總結

本研究對我國栽培的優(yōu)良品種威寧黑麥進行了基因組測序，基因組組裝大小為7.74 Gb，其中93.67%被掛載到7條染色體上，重復序列占總基因組的90.31%。并基于高質量基因組揭示了全基因組基因復制及其對淀粉生物合成基因的影響，解析了復雜儲存蛋白基因座位點、早抽穗性狀的基因表達特征以及黑麥中與馴化相關的染色體區(qū)域和基因座。本次研究結果對黑麥的基因組特性及其農(nóng)藝性狀調(diào)控基因有了新的認識，獲得了對進一步研究黑麥馴化遺傳基礎可能有用的染色體區(qū)域和基因座。威寧黑麥基因組組裝對于破解黑麥基因組生物學，深化比較谷類基因組學研究，加速黑麥及相關谷類作物的遺傳改良具有重要價值。

?

PART TWO 團隊介紹

河南農(nóng)業(yè)大學“植物基因組學和分子育種” 研究團隊組建于2017年4月，其主要研究方向是深入解析植物重要生命過程以及小麥等作物關鍵農(nóng)藝性狀控制的分子機理，發(fā)掘具有科學意義和育種實用價值的基因與優(yōu)良變異，創(chuàng)制具有實用價值的新種質和新品種，為促進小麥等作物的遺傳改良作出基礎性、戰(zhàn)略性、前瞻性和實用性貢獻。

團隊現(xiàn)有5位教授（王道文、楊建平、茍明月、鄭旭、姬祥），其中國家杰出青年基金獲得者1名。團隊擁有實驗室面積1440平方米，建有基因組學、分子生物學、分子育種等設施。團隊主要學術帶頭人為王道文研究員，目前為河南農(nóng)業(yè)大學農(nóng)學院擔任特聘教授，主要研究植物分子生物學和小麥遺傳改良。王道文研究員以第一或通訊作者身份在SCI刊物 (Nature、Science、Nature genetics、PNAS、Plant Cell、Plant Journal、New Phytologist、Molecular Breeding等) 發(fā)表論文50余篇，獲國家杰出青年基金支持，獲得中華農(nóng)業(yè)科技獎一等獎1項 (第二完成單位和第二完成人)?，F(xiàn)主持國家重點研發(fā)計劃項目“主要農(nóng)作物品質性狀形成的分子基礎”的研究，并擔任核心期刊“作物學報”常務編委、SCI刊物Scientific Reports編委。

團隊成立后取得了一系列重大科研進展及成果，包括黑麥基因組測序和深度分析工作。在黑麥基因組解析工作中，該團隊利用三代DNA技術，結合Hi-C染色質構象捕捉技術的應用，獲得了黑麥全基因組范圍內(nèi)整個染色質 DNA 在空間位置上的關系，在此基礎上實現(xiàn)了黑麥基因組序列的組裝、排序和定向，以及整個基因組的注釋分析工作，成果最終發(fā)表在國際科研期刊《Nature genetics》上。

百邁客生物科技有限公司(Biomarker?Technologies)成立于2009年5月，是一家高新技術企業(yè)，公司基于功能基因組學的BT+IT(高通量基因測序技術、生物信息技術、高性能生物云計算技術)技術，面向廣大科研工作者提供以綜合科技服務、BMKCloud生物云分析、二代及三代高通量測序分析，儀器，試劑等為主的科技服務。公司目前擁有Nanopore、Pacbio、Illumina、Waters、10X genomics等測序服務平臺，構建了基于云架構的自主生物云計算分析平臺—百邁客云，提供涵蓋人重外顯子、三維基因組、單細胞與空間轉錄組、基因組de?novo、轉錄調(diào)控、微生物、群體遺傳、質譜及表觀遺傳等研究方向的技術服務。自公司成立起，服務客戶在《Cell》、《Nature》、《Nature Genetics》、《Nature Communications》、《Plant Cell》、《PNAS》等學術刊物發(fā)表論文數(shù)千篇。