中国农业大学农学院宋伟彬团队解锁黍稷基因宝藏为谷物育种提供种质和基因资源

　　中国农大新闻网讯 近日，中国农业大学玉米生物育种全国重点实验室、国家玉米改良中心宋伟彬团队在《自然·遗传学》（ Nature Genetics ） 上发表论文《黍稷优异基因资源挖掘》（ Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits ）。

 

　　玉米，作为我国的第一大粮食作物，对于保障国家粮食安全具有至关重要的意义。长期的人工选择造成了优良等位基因的缺失，导致种质基础趋于狭窄，严重限制了玉米育种创新。因此，开展玉米近缘物种组学研究，通过比较基因组学技术手段挖掘玉米近缘种特有、特优基因资源并应用于玉米育种，对玉米种质基础的拓展具有重要的理论和实践指导意义。玉米近缘种黍稷起源于中国而且是最为古老的农作物之一，具有生育期短、耐旱耐瘠薄、适应盐碱胁迫等特点，暗示着黍稷基因组中蕴藏着丰富的优异基因资源。宋伟彬前期与合作者组装了黍稷 （陇糜 4 号） 基因组，共注释出 63,671 个编码基因，进一步与玉米自交系 Mo17 （39,756个编码基因） 的比较基因组学分析发现，63,671 个编码基因中 11521 个是黍稷特有 （Nature Communications, 2019）；为了从 11,521 个黍稷特有基因中发掘出与农艺性状和育种选择相关的基因，团队开展了育种过程中黍稷基因组变异图谱构建工作。

 

　　该研究团队系统地完成了中国野生黍稷种质资源的收集和基因资源挖掘工作。团队从内蒙古、陕西、甘肃、辽宁、吉林、河北、山西、宁夏、黑龙江的等10个省份的路边、丘陵、山地、农田等地收集了 967 份野生黍稷，结合 937 份栽培黍稷组成了一个 1,904 份的黍稷大群体。利用双末端 150 bp 测序获得黍稷群体的全基因组重测序数据，产生了 65.1 Tb 的数据，平均基因组覆盖度为 40.3 倍。根据基因组覆盖度和 IBS，对野生黍稷和栽培黍稷进行筛选，过滤低覆盖度和遗传背景相似的品种。最后保留了 1,094 份黍稷进行后续分析，包含 316 份野生黍稷和 778 份栽培黍稷。

 

　　研究人员利用鉴定的 5,603,840 个 SNPs 解析了黍稷的群体结构，将黍稷群体分为四个亚群：LR （Liao River，辽河野生黍稷亚群）、YY （Yin Mountain and Yan Mountain，阴山-燕山野生黍稷亚群）、LP （Loess Plateau，黄土高原野生黍稷亚群）、Cultivated （栽培黍稷亚群），发现野生黍稷和栽培黍稷具有明显的遗传差异。进一步分析发现，黄土高原野生黍稷亚群与栽培黍稷分化时间最为接近 （4,600 年前），这与考古学的时间相吻合，也证实了现代栽培黍稷起源于黄土高原地区的单一起源理论，为黍稷的起源中心提供了理论支持。这些结果填补了中国黍稷遗传多样性研究的空白，并明确了现代栽培黍稷在中国的精确起源地里位置，有助于我们更加全面地了解中国黍稷的遗传特征、起源及其与栽培黍稷的关系 。

 

　　分析发现，在驯化过程中黍稷基因组不仅发生了单核苷酸变异和小片段的序列变异，也存在大片段的结构变异。每个野生黍稷材料平均含有 1,431 个 SVs，每个栽培黍稷平均含有 538 个 SVs，每个野生黍稷的 SVs 数量接近是栽培黍稷的 3 倍。为了分析 SVs 在驯化过程中的动态变化，该研究分析了 SVs 和 TE 的关系。发现 SVs 主要发生在 Gypsy 类型的 TE 上，另外 Gypsy 类型的 SVs 在黍稷驯化过程中受到选择，可能在驯化过程中发挥重要作用。为了探究 SVs 对基因的影响，分析了 SVs 和基因的位置关系，结果显示结构变异使基因组 10% 的基因为 PAV。其中，栽培黍稷特有 2,149 个 PAV 基因，野生黍稷特有 2,892 个 PAV 基因，这些基因更多表现在影响植株的抗性和驯化性状。比如，栽培黍稷群体中的一个长度为 41 kb 的 deletion 引起影响籽粒大小基因 PmTGW6 的缺失、另外一个长度为 10.3 kb 的 deletion 引起影响落粒性基因 PmSh1 的缺失。

 

　　为了鉴定与产量相关的 QTL 位点，研究人员利用鉴定到的 SNPs、Indels、SVs 对 12 个农艺性状进行全基因组关联分析，并鉴定到 186 个 QTL。其中与籽粒大小显著关联的 PmGW8 外显子上的 T 碱基缺失导致该基因发生移码突变，并且该基因在黍稷驯化过程中受到了强烈的选择；鉴定到两个与抽穗期显著关联的基因 PmFT 和 PmPRR37 ，通过在拟南芥中过表达验证了 PmFT 能够正向调控花期，可以显著缩短开花时间；鉴定到了与株高显著关联的水稻 sd1-1 的同源基因 Pmsd1 基因，并且通过在黍稷中敲除 Pmsd1，证实了 Pmsd1 参与了调控黍稷的株高和总穗长。综合分析发现，野生黍稷能鉴定更多籽粒相关的 QTL，暗示着野生黍稷资源在黍稷籽粒性状的改良育种过程中可发挥重要作用。

 

　　野生黍稷和栽培黍稷在我国北方地区具有广泛的分布，因而需要适应不同的生态和地理环境。开花期是植物适应当地环境的重要性状，解析花期的分子进化机制是了解黍稷适应性的关键。研究人员利用鉴定到的两个抽穗期 QTL，qHD2 和 qHD5 来探究野生黍稷和栽培黍稷的适应机制。结果表明野生黍稷通过选择 qHD5 的不同单倍型调节生长周期来适应不同纬度的环境，栽培黍稷通过选择 qHD2 的早花单倍型来适应高纬度的光周期条件 。

 

　　综上所述，通过野生黍稷和栽培黍稷基因组多态性分析，结合比较基因组学技术手段共鉴定出 3,398 个育种过程中受选择基因，其中 324 个是黍稷特有基因；调控 12 个农艺性状的 1698 个候选基因中 262 个是黍稷特有，去除掉 7 个重叠基因，共计鉴定到 579 个黍稷特有基因，这为利用生物育种技术拓展玉米种质基础提供了重要基因资源。

 

　　中国农业大学玉米生物育种全国重点实验室、国家玉米改良中心、农学院 宋伟彬 教授为该论文通讯作者。中国农业大学博士生 陆琼、中国农业大学农学院、分子设计育种前沿科学中心 赵海楠 研究员、博士生 张政权、博士生 白雨禾、中国农业大学 赵海铭 副教授、河北省农科院 刘国庆 研究员为该论文共同第一作者。中国农业大学博士生 郑云霄、博士生 龚慧慧、博士生 洪祥德、博士生 李百川、博士生 白岩、已毕业硕士 陈令玮、邓喜真、刘天翔、中国农业科学院 陆平 研究员、刘敏轩 助理研究员、通辽市农牧科学研究所 文峰 助理研究员、山西省农科院 王伦 研究员、黑龙江省农科院种质资源所 李志江 副研究员、山西农业大学 李海 研究员、河北省农科院 李海权 研究员、全国农业技术推广中心 张力科、马文慧、刘春青、中国农业大学科研助理 赵海月、中国农业大学 辛蓓蓓 副教授、陈建 副教授、鄂立柱 副教授在黍稷资源收集方面和基因型分析方面做出了重要贡献。赖锦盛 教授在该研究的实验设计和分析提供了重要帮助和指导。

 

　　该研究得到了国家重点研发计划 （2021YFD1200701）、国家自然科学基金委 （32271541, 32272143）、科技创新2030 （2022ZD04020），种质遗传改良 （2022ZD04020, 2022010202, B21HJ0509）、基本科研业务费 （2023TC019）、中央高校基本科研业务费专项资金资助 （2023TC019） 等项目的资助，通州于家务国际种业园区教授工作站给予了大力支持。

 

　　原文链接：https://www.nature.com/articles/s41588-024-01718-6#Sec2