谷子基因组学研究进展

Research Advances in Foxtail Millet Genomics

摘要

作为古老的驯化作物之一，谷子在中国北方的干旱和半干旱地区广泛种植。谷子是重要的杂粮作物，其籽粒富含蛋白质、糖类、维生素等营养物质，具有C4光合、小的二倍体基因组、抗逆、适应性强等特点。自2012年谷子基因组测序项目完成以来，谷子已经进入了基因组学时代。谷子基因组学的发展可以分为3个阶段：结构基因组、功能基因组、比较基因组。二代测序技术的出现，使得基于测序的基因分型和全基因组关联分析极大地推动了谷子遗传学的研究。在健康中国战略的时代背景下，谷子已成为改善居民膳食结构和推动有机旱作高质量发展的优势作物。近年来，谷子产业发展呈上升趋势，但品种选育方面缺乏突破，不能满足现代有机旱作需求，主要原因在于谷子育种仍然依托传统育种手段。谷子基因组学的发展为分子育种奠定了基础，文章综述了谷子全基因组测序历史与现状、遗传转化体系的研究进展以及基于谷子全基因组序列发掘基因资源取得的成果。谷子功能基因组学与传统育种技术的有机结合必将促进高效、精准谷子育种体系的构建，加快突破性谷子品种的培育，满足产业发展和消费者的需求。

Abstract

Foxtail millet(Setaria italica(L.) P. Beauv.), a venerable domesticated crop, is widely cultivated in arid and semi-arid regions of northern China. As an important orphan crop, foxtail millet boasts traits such as C4 photosynthesis, a compact diploid genome, robust stress tolerance, and wide adaptability. Its grains are rich in protein, sugar, vitamins, and miscellaneous nutrients. The completion of the foxtail millet genome sequencing project in 2012 marked the dawn of genomics for this ancient crop. The research of the foxtail millet genome has unfolded through three key stages: structural genomics, functional genomics, and comparative genomics. The advent of next-generation high-throughput DNA sequencing technologies has facilitated sequencing-based genotyping and genome-wide association studies(GWAS), propelling genetic research on foxtail millet. In the context of “Healthy China”, foxtail millet has emerged as a pivotal crop for improving dietary patterns and driving the sustainable development of organic dryland farming. While the foxtail millet industry has seen remarkable growth in recent years, the paucity of breakthroughs in breeding has hindered its alignment with the demands of modern organic dryland farming, largely due to the constant reliance on traditional breeding methods. Significant strides in foxtail millet genomics have been witnessed, laying the groundwork for molecular breeding. This review provided a comprehensive overview of the historical and current landscape of whole-genome sequencing in foxtail millet, the advancements in genetic transformation systems, and the achievements on gene discovery stemming from the whole-genome sequence. The effective integration of functional genomics with traditional breeding methods will guarantee to enhance the establishment of an efficient and precise breeding system for foxtail millet, hastening breakthroughs in breeding and aligning with the requirements of related industry and consumers.

关键词

谷子 / 基因组 / 全基因组测序 / 遗传转化 / 基因挖掘

Key words

foxtail millet / genome / whole-genome sequencing / genetic transformation / gene discovery

Author summay

侯思宇（1980-），男 ，山西陵川人 ，教授，博士，主要从事杂粮品质及分子育种研究工作。

随着2000年世界上第一个测序完成的模式植物——拟南芥（Arabidopsis thaliana）全基因组序列的公布^[1]，推动了基于基因组学、分子生物学等相关技术开展的植物生长发育、抗逆、抗病、养分利用等大量的模式植物功能基因组学研究^[2]，同时也极大地促进了作物科学领域的发展，其中水稻（Oryza sativa）基因组学研究最为深入，2002年水稻基因组测序完成^[3]，利用多组学、全基因组关联研究揭示了水稻养分利用、生物和非生物响应、生殖发育等重要分子机制，同时高通量表型平台、精确基因组编辑工具等新技术也广泛应用于水稻研究^[4]。随着相关技术的逐渐成熟，基因组学研究在多种主粮和杂粮作物中迅速发展。

谷子（Setaria italica）起源于我国，是最早完成基因组测序的杂粮作物^[5]。谷子基因组较小，且生育期短，具备模式植物的特点。同时，谷子作为C4植物，可以用于揭示C4植物高光效、抗逆耐瘠薄的特殊机理^[6]。近年来，我国科学家利用EMS诱变获得了株高和生育期与拟南芥极为接近的“模式化”谷子——xiaomi^[7]，同时，基于110份谷子种质的高质量泛基因组的完成^[8]，为推动谷子成为C4模式植物、深入开展基因组学研究奠定了坚实的基础。本文概述了自谷子基因组测序研究工作开展以来，国内外科学家围绕谷子开展的功能基因组学研究，包括谷子基因组的情况，以及利用谷子基因组开展的遗传转化体系的开发，调控特色营养物质积累和抗逆性调控等重要功能基因的发掘工作，这既是对谷子分子水平研究的一次梳理，也是对谷子功能基因组学研究的未来研究方向的一次展望。

1 谷子全基因组测序历史与现状

华大基因与河北张家口市农业科学院等单位的科研人员于2012年利用全基因组鸟枪法结合二代测序技术对中国北方广泛种植的谷子品系“张谷”进行了全基因组测序和组装，获得了大小约为423 Mb（Contig N50为25.4 kb，Scaffold N50为1.0 Mb）的谷子全基因组序列图谱^[5]。通过基因组注释和分析发现，重复序列约占整个基因组的46%，大约含有38 801个蛋白质编码基因，其中有81%可以被表达。研究人员通过谷子、水稻和高粱之间的共线性分析鉴定到了关键的染色体重组事件，明确了谷子9条染色体是在3次染色体重组事件后形成的（在谷子和水稻分化后2次，谷子的2号和9号染色体分别由水稻的7号和9号，3号和10号染色体融合而成；在谷子和高粱分化后一次：谷子的3号染色体由水稻的5号和12号染色体融合而成），为禾本科的进化提供了见解。此外，研究人员使用Illumina GA II平台对谷子光热敏的雄性不育系“A2”进行了10×深度的重测序，并利用“张谷”和“A2”的F₂群体绘制了高密度的遗传图谱，将该图谱与谷子基因组序列及后代表型等信息相结合定位到了一个谷子抗除草剂相关基因。

美国国家能源部所属联合基因组研究所的科研人员采用全基因组鸟枪法结合一代测序技术（Sanger测序法）对谷子品种豫谷1号进行了全基因组测序和组装，并与上述“张谷”基因组序列在同一期刊同时发表相关论文。他们获得了大小约为400 Mb（Contig N50为126.3 kb，Scaffold N50为47.3 Mb）的谷子参考基因组^[9]，覆盖了谷子近80%的基因组，超过95%的基因区域。为促进基因表达分析，利用不同组织、不同发育时期的豫谷1号的mRNA序列生成了超过130万个表达序列标签（EST）文库和超过5.8亿个RNA-Seq reads。同时，该研究还使用Illumina GA II平台对谷子的野生近缘种狗尾草A10进行了重测序，利用谷子自交系B100和狗尾草A10杂交后代构建了重组自交系（RIL）群体，该RIL群体具有992个SNP标记。该研究确定了差异单核苷酸多态性密度、转座子分布、小RNA含量、染色体重排和偏分离区域，通过比较5个已测序的禾本科植物基因组，研究了狗尾草属的广泛适应性的遗传基础。

美国唐纳德丹佛斯植物科学中心研究人员于2020年使用PacBio和Illumina HISeq测序平台更新了狗尾草A10.1的参考基因组，完成了端粒到端粒的染色体组装，获得了大小为395.1 Mb的高质量参考基因组（2.0版本）^[10]。2.0版本的狗尾草参考基因组的组装质量更高，Contig N50为11.2 Mb，在9条染色体上基因组组装的完整性为99.95%，共注释得到38 334个编码蛋白基因和14 125个可供选择的转录本。此外，该团队基于二代测序使用Illumina 2×150PE文库从头组装了598份狗尾草高质量基因组，构建了狗尾草泛基因组数据库，进一步利用高质量基因组、泛基因组鉴定到一个调控狗尾草种子落粒性基因Sh1并通过CRISPR-Cas9基因编辑技术对其功能进行验证，同时也证实了该基因在谷子中也是影响落粒性的关键位点。

同年，山西农业大学杂粮分子育种团队对名优谷子品种晋谷21号EMS诱变获得的超早熟谷子突变体xiaomi构建了功能基因组学研究体系。利用PacBio平台和二代、三代基因测序技术对xiaomi进行了基因组测序和组装，获得了大小为429.94 Mb（Contig N50高达18.8 Mb，Scaffold N50为42.41 Mb）的高质量参考基因组^[7]。xiaomi基因组序列中重复序列约占55.19%，根据从头预测、同源基因预测以及Iso-seq和RNA-seq分析结果，共注释到33 789个蛋白质编码基因。xiaomi基因组仅包含48个Gap且组装的错误率约为0.001%，在覆盖率、Contig N50等方面的表现均优于张谷、豫谷1号参考基因组。此外，该研究对迷你谷子xiaomi不同发育阶段的11个组织进行了转录组测序，结合高质量参考基因组，构建了首个谷子全生育期动态基因表达图谱和谷子多组学数据库，不仅促进了谷子功能基因组学的研究，还使得xiaomi成为C4植物功能研究的理想模型系统。

中国农业科学院作物科学研究所特色农作物优异种质资源发掘与创新利用团队于2023年建立了首个谷子图形结构泛基因组，在此基础上鉴定了多个关键基因，展示了泛基因组无与伦比的优势^[8]。他们从1 844份不同亚组、区域和生态型的谷子和狗尾草品系中选取了35个谷子野生种、40个农家品种和35个现代育成品种。利用PacBio、Illumina和Bionano测序技术从头组装（de novo assembly）代表谷子和青狗尾草最广泛的多样性的110份核心参考水平基因组。核心种质包含了对谷子育种或研究作出重大贡献的种质，如选育骨干亲本（60日和矮88）、高食用和蒸煮品质品种（晋谷21号和黄金苗）、强耐旱性品种（中谷2号）、广泛的气候适应性品种（豫谷18号），以及易于遗传转化用于基因功能分析的材料（Ci846），这些材料涵盖了株型结构、穗形和产量等性状多样性。这些材料包含了1 844份狗尾草属中85%以上的单核苷酸多态性（SNP）变异。谷子泛基因组中包含73 528个基因家族的全基因组。通过整合112份谷子和青狗尾草泛基因组中的107 151个插入、76 915个缺失和363个倒位变异，并将其整合到豫谷1号参考基因组序列中，构建了首个C4作物的泛基因组。泛基因组为探索谷子群体进化、驯化与改良、功能基因组学和广泛的适应性等基础研究，及其农艺学应用方面提供了推动力。利用泛基因组，鉴定了4 582个与驯化相关、152个与育种改良相关的结构变异，同时共鉴定出680个在驯化改良过程中被持续选择的基因；鉴定出有关落粒基因Sh1（855 bp的缺失）和籽粒大小基因SiGW3（366 bp的缺失）等多个结构变异基因；确定了共1 084个与表型显著相关的位点。泛基因组为未来的谷子分子育种提供了新平台，同时也为其他杂粮作物精准育种选择奠定了理论基础。

2 谷子转化体系的构建

谷子是C4作物，且基因组小、抗逆性强，具有成为禾本科功能基因组学模式作物的潜力，因此，近年来关于其转化体系的研究得到了广泛关注。

谷子遗传转化方法包括基因枪法、农杆菌介导法、PEG介导法、超声波介导法、花粉管通道法、子房注射法、基因编辑技术等。基于狗尾草愈伤诱导的成功，近些年来农杆菌介导法被应用于谷子的遗传转化相关研究。早在20世纪70年代日本学者就开始了谷子的组培研究，后续各国学者在此基础上进行了更为深入的研究。2007年中国农业大学LIU等^[11]探究了多种影响谷子遗传转化效率的因素，并确定了最优的转化条件，为农杆菌侵染愈伤组织的转化系统打下了基础。2020年山西农业大学YANG等^[7]以迷你谷子突变体xiaomi开发了高效的农杆菌介导的谷子遗传转化方案，由于其生命周期短、植株小，与拟南芥相似，一年可以在培养室中生长5~6代，更有利于谷子功能基因组学的研究^[7]。2020年PRIYANKA等^[12]优化了谷子简单有效的再生和转化方案，同年，SANTOS等^[13]利用成熟干燥的种子作为外植体进行农杆菌转化，大大提高了转化率和再生率。大大提高了转化率和再生率。2021年杨澜等^[14]利用基于农杆菌介导的CRISPR/Cas9系统转化谷子，胞嘧啶和腺嘌呤碱基编辑系统进行单、多基因敲除和单碱基替换，创造了一种纯合抗除草剂种质。除愈伤组织外，2021年杨澜等^[14]还建立了稳定的谷子体外茎尖遗传转化体系。

青狗尾草是谷子的近缘野生种，二者核型基本相同、基因组大小相近，是研究谷子驯化、发掘基因资源的重要材料。2016年SAHA等^[15]通过蘸穗法将目的基因成功整合到狗尾草基因组中，且在后代中成功表达。随后，2017年赵辉等^[16]通过高效胚性愈伤诱导技术，大大提高了狗尾草的转化效率。

谷子和狗尾草转化体系的建立为谷子功能基因组学研究和遗传改良奠定了基础。

3 基于谷子全基因组序列发掘基因资源

3.1 基于参考基因组，利用转录组探究谷子基因表达模式与调控机制

随着谷子参考基因组的公布发表，谷子重要性状相关基因表达模式和调控机制的研究成果呈“井喷”式增长。当前转录组分析主要聚焦在谷子非生物和生物胁迫、穗发育、初级代谢和次生代谢物等机制解析方面。尤其在转录因子家族成员鉴定及表达模式方面已有大量的研究文献。在谷子基因组中，鉴定出147个NAC、122个C2H2类锌指蛋白、225个WD40蛋白、209个MYB、171个AP/ERF、35个Dof、47个HD-zip、72个MADS-box、39个NF-Y及37个CCT，不同非生物胁迫（干旱、盐、冷等）和激素（水杨酸、脱落酸、茉莉酸甲酯、乙烯等）处理条件下，上述转录因子家族成员多样性的表达模式暗示了潜在的生物学功能分化，从中鉴定到一系列胁迫响应特异的候选基因，如SiNAC128、SiAP2/ERF-069、SiAP2/ERF-103、SiAP2/ERF-120、SiNF-YA1、SiNF-YB8和SiMADS51等^[17-26]。谷子在低氮胁迫下，利用转录组分析找到75个差异表达的转录因子，其中SiMYB3基因能够调控下游生长素合成相关基因TAR2表达，进而促进根系发育^[27]。结合干旱胁迫的时序动态转录组分析，鉴定出13 294个响应干旱胁迫基因，包括842个昼夜节律基因，初步解析了干旱胁迫与昼夜节律交叉调控谷子发育的分子机制^[28]。谷子籽粒灌浆期动态转录组揭示了11 399个差异表达基因，包括902个转录因子，涉及到淀粉合成、细胞壁活性、激素信号转导及多胺代谢途径^[29]。此外，涉及到非生物胁迫响应的73个谷子谷胱甘肽转移酶及331个细胞色素P450单加氧酶基因家族成员，也通过全基因组生物信息学及转录组分析，分别鉴定出37个组织特异性表达和21个响应非生物胁迫的谷胱甘肽转移酶基因家族成员，以及6个细胞色素P450单加氧酶基因表达水平响应不同的除草剂处理（硝磺草酮、双氟磺草胺、烟嘧磺隆、氟草定及拿普净）^[30-31]。丛枝菌根真菌定植显著增强谷子籽粒产量，比较转录组分析揭示了超过2 000个基因受到丛枝菌根真菌的诱导响应，GO富集分析表明大部分基因富集在氮同化和转运、细胞壁重组和木质化途径^[32]。低单宁谷子品种可以吸引红蜘蛛取食，从而保护板栗树免受虫害，比较转录组揭示了335个共有的差异表达基因，KEGG^[33]富集发现PTAL、CCR及POX基因涉及到苯丙烷类合成途径的单宁生物合成分支途径，此外，预测了20个转录因子，包括bHLH、WRKY、FAR1、ERF、C2H2等参与调控单宁生物合成。板栗树与谷子套种系统有效地提高园林土地利用率及增加其产量和品质，遮荫条件下不同基因型谷子转录组分析鉴定出9个差异表达基因，功能注释为HSP70-8、HsfA2、SPDS、SPMS等基因涉及光系统和热胁迫响应途径^[34]。谷子紫叶基因型具有较高的胁迫容忍性，比较转录组分析，鉴定出9个差异表达基因与花青素合成相关，可作为分子育种的潜在功能标记^[35]。在黄米和白米谷子品种穗发育阶段，基于比较转录组和类胡萝卜素靶向代谢组联合分析，构建了基因共表达网络，鉴定出54个类胡萝卜素代谢途径相关基因及3个转录因子协同调节类胡萝卜素代谢流^[36]。联合脂质组和转录组分析，鉴定出2 633个脂质分子，归属于13种甘油磷脂、11种甘油糖脂、4种鞘脂类及脂肪酰和甾醇类，同时鉴定出9个差异表达基因与甘油二脂代谢途径相关^[37]；谷子籽粒不饱和脂肪酸含量远大于饱和脂肪酸含量，152个差异表达基因与脂肪酸代谢和植物甾醇合成相关^[38]。禾生指梗霉是引起谷子白发病的一类真菌，其病害发生使晋谷21号造成大量减产，转录组分析表明，禾生指梗霉侵染谷子后导致萜类生物合成酶（CPS和KS）高表达，赤霉素、脱落酸及生长素合成途径相关基因表达量上调，初步揭示了内源激素调控谷子白发病表型的分子机制^[39]。环境CO₂浓度升高可提升谷子高光效、生物量及产量，转录组分析鉴定出一系列差异表达基因，涉及到细胞重组、茎段发育、气孔导度、碳固定、糖酵解和糖异生途径^[40]。龙谷25属于低钾胁迫敏感基因型，转录组分析发现了1 982个差异表达基因响应低钾胁迫，其中一个响应低钾胁迫的转录因子SiMYB3功能分析表明，超表达拟南芥可促进根系延长及应对钾素缺乏能力^[41]。在谷子noncoding RNA（LncRNA、siRNAs、microRNA）测序分析方面，也初步明确了其种类、染色体分布特征、表达模式、共表达网络及谷物类保守进化特征，预测了响应干旱和脱水胁迫、水杨酸处理的noncoding RNA及其调控的靶基因^[42-47]。此外，利用基因组和转录组分析，还鉴定出一些重要的基因家族，如SWEET、SET-domain、HAT（Histone acetyltransferase）、AATs（Amino acid transporters）和CDPKs（Calcium-dependent protein kinases）等，找到与产量、蛋白品质、抗非生物胁迫及甲基化相关的一些候选基因^[26，48-52]。

综上所述，利用基因组和转录组分析，针对谷子的生长发育调节、生物胁迫与非生物胁迫及代谢调控相关的基因发掘方面已做了较为重要的创新性工作，但仍存在不同基因型不同处理所获得的结果不尽相同，其相关共性的规律尚需进一步总结归纳，挖掘的候选基因多数尚未验证其功能，需要在生信分析的基础上，展开实质性的基因功能验证和机理深度解析，为将来谷子分子育种提供可利用、主效和可靠的基因资源。

3.2 基于GWAS、转录组定位重要性状基因

谷子作为起源于我国的重要粮食作物，随着谷子基因组的公布，基于全基因组关联分析、QTL定位等手段，已有很多调控谷子重要性状的关键基因或QTL被发掘。JIA等^[53]基于全基因组重测序对916份谷子材料的47个农艺性状进行全基因组关联分析，在5个环境下共鉴定到512个显著关联的SNP位点；LI等^[54]基于312个谷子材料的全基因组关联分析，挖掘到影响谷子高海拔适应性的关键基因SiPRR37；VANDANA等^[55]和LIU等^[56]分别对142个和407个谷子品种的农艺性状及穗型性状进行全基因组关联分析，挖掘了大量的显著关联SNP。大规模mGWAS研究也被报道。基于广泛靶向代谢组学与GWAS的协同分析，挖掘了控制谷子米色及类胡萝卜素含量的关键基因SiPSY^[57]。通过谷子泛基因组图谱，鉴定了多个与谷子驯化及育种改良的染色体结构变异及重要候选基因^[8]。

除全基因组关联分析外，基于遗传作图，谷子中也发现了多个重要QTL。DOUST等^[58]通过QTL分析在来源于谷子和青狗尾草的双亲群体中检测到14个与穗部性状相关的QTL。在另外2个QTL定位研究中，也分别鉴定出12、32个QTL与穗部性状相关，如PWP、GWP、TGW等^[59-60]。这些研究大部分是基于少量的多态性位点，并且很少有QTL被鉴定。最近一项基于12个环境的QTL定位研究，确定了159个与穗性状相关的QTL，深化了对谷子穗部和产量形成的遗传基础的理解^[61]。然而，基于双亲群体的QTL结果严重依赖于群体和环境的遗传背景，这极大地限制了这些QTL在谷子育种计划中的广泛适应性和稳定性。总的来说，尽管谷子上已有部分基因或QTL被发掘，但是与小麦、水稻、玉米等大作物相比，还远远不足。进一步充分发掘利用优异种质资源及其富含的基因资源是谷子现代化分子育种研究的基础。

3.3 基于参考基因组对基因家族的鉴定

谷子参考基因组测序的完成对于谷子功能基因组的研究尤其是谷子基因家族的鉴定提供了便利，近年来涉及到抗逆胁迫、生长发育、以及营养吸收相关的基因家族相继被报道，为谷子功能基因的深入研究提供了思路与参考。

3.3.1 调控植物生长发育过程的基因家族的分离与鉴定

MADS-box基因家族部分成员参与谷子穗分生组织发育^[62]；SPL基因参与组织分化以及种子发育^[63]；bHLH转录因子（187个家族成员）参与谷子根、茎、叶、花以及果实的发育调控^[64-65]；AAT基因家族成员在谷子籽粒品质形成过程中具有重要作用^[49]；GRAS基因家族成员调控谷子株高，并影响谷子的穗质量^[66]；GRF基因家族倾向于在发芽的种子、芽、花蕾和嫩叶等活跃生长和发育的组织中高表达，暗示其可能参与谷子生长发育的调控^[67-70]；LBD家族成员在谷子根系中特异表达^[71]，CCT转录因子（39个成员）中有15个SiCCT基因与谷子抽穗期具有重要关系^[72]。以上工作为谷子基因功能的深入研究提供了资源与参考。

3.3.2 参与谷子逆境胁迫的基因家族的分离与鉴定

参与谷子干旱、盐碱以及低氮、低磷胁迫等非生物胁迫的基因家族相继被报道。其中包括通过氧化各种脂质参与植物的生理活动的谷子脂氧合酶（LOXs）基因（12个成员），调控脂肪酸和蜡质的合成与积累的酮脂酰辅酶A合酶（β-ketoacyl-CoA synthase，KCS，33个成员）^[73]，LBD转录因子（33个家族成员）^[71]，糖转运蛋白（Sugar transporter proteins，STPs）家族成员^[74]，SPL（18个家族成员）基因家族^[75]，植物氨基酸转运载体（Amino acid transporter，AAT）^[49]，组蛋白乙酰转移酶（HAT）基因^[52]。另外，WRKY转录因子^[76]、MYB^[20]、CCCH类锌指蛋白（27个成员）^[77]以及GRF^[78]家族成员也参与谷子多种非生物胁迫响应；ANK家族基因成员参与低温、干旱、盐等胁迫响应^[79]；谷子APX家族成员受干旱胁迫和高氮诱导后，表达显著上调^[80]。

3.3.3 调控营养元素转运吸收功能的基因家族的鉴定

参与谷子硒、叶酸等微量元素吸收以及氮、磷、钾等营养元素吸收转运的基因家族近年来也被报道，其中谷子SULTR基因家族成员SiSULTRC3可以提高植株转运SeO₄^2-的能力^[81]。FPGS基因家族的SiFPGS1-3、SiFPGS1-4、SiFPGS-2基因表达与谷子叶酸含量呈正相关^[82]。Shaker K⁺通道基因家族参与植物钾的摄取和分布^[83]。ZIP家族中的基因编码转录本可以储存和运输二价金属微量营养素，特别是铁（Fe）和锌（Zn）。NAC、bZIP和bHLH是SiZIP基因中存在的主要Fe和Zn响应转录因子^[84]。PHT1家族转运蛋白基因在低磷胁迫下，可促进植株对磷的吸收^[85]。MYB-CC基因已经被证明可调控植物中磷酸盐吸收^[86]。NRT基因家族由NRT1/PTR，NRT2和NRT3亚家族组成，其在硝酸盐从土壤到植物的吸收和运输中起关键作用^{[87- 88]}。

4 展望

2012年全基因组测序的完成为谷子功能基因组学揭开了序幕，2023年泛基因组的发布吹响了谷子分子育种的集结号。在丰富种质资源的基础上，谷子基因资源的发掘将走上快车道，越来越多的重要性状相关基因及其功能将会被发现并用于分子育种，分子育种与传统育种的紧密结合将极大地提升谷子育种效率，推动突破性品种的培育。

参考文献

[1] INITIATIVE T A G. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana [J]. Nature，2000，408：796-815.

[2] PROVART N J， ALONSO J， ASSMANN S M，et al. 50 years of Arabidopsis research：highlights and future directions[J]. The New Phytologist，2016，209（3）：921-944.

[3] GOFF S A， RICKE D， LAN T H，et al. A draft sequence of the rice genome（Oryza sativa L.ssp.japonica）[J]. Science，2002，296（5565）：92-100.

[4] CHEN R Z， DENG Y W， DING Y L，et al. Rice functional genomics：decades’ efforts and roads ahead[J]. Science China Life Sciences，2022，65（1）：33-92.

[5] ZHANG G Y， LIU X， QUAN Z W，et al. Genome sequence of foxtail millet（Setaria italica） provides insights into grass evolution and biofuel potential[J]. Nature Biotechnology，2012，30：549-554.

[6] PENG R H， ZHANG B H. Foxtail millet：a new model for C4 plants[J]. Trends in Plant Science，2021，26（3）：199-201.

[7] YANG Z R， ZHANG H S， LI X K，et al. A mini foxtail millet with an Arabidopsis-like life cycle as a C4 model system[J]. Nature Plants，2020，6：1167-1178.

[8] HE Q， TANG S， ZHI H，et al. A graph-based genome and pan-genome variation of the model plant Setaria [J]. Nature Genetics，2023，55：1232-1242.

[9] BENNETZEN J L， SCHMUTZ J， WANG H，et al. Reference genome sequence of the model plant Setaria [J]. Nature Biotechnology，2012，30：555-561.

[10] MAMIDI S， HEALEY A， HUANG P，et al. A genome resource for green millet Setaria viridis enables discovery of agronomically valuable loci[J]. Nature Biotechnology，2020，38：1203-1210.

[11] LIU Y H， YU J J， AO G M，et al. Factors affecting Agrobacterium tumefaciens-mediated transformation of foxtail millet[J]. Chinese Journal of Biochemistry and Molecular Biology，2007（7）：531-536.

[12] SOOD P， SINGH R K， PRASAD M. An efficient Agrobacter- ium-mediated genetic transformation method for foxtail millet（Setaria italica L.）[J]. Plant Cell Reports，2020，39（4）：511-525.

[13] SANTOS C M， ROMEIRO D， SILVA J P，et al. An improved protocol for efficient transformation and regeneration of Setaria italica [J]. Plant Cell Reports，2020，39（4）：501-510.

[14] 杨澜，刘雅，项阳，等. 谷子茎尖体外遗传转化体系的建立与优化[J]. 植物学报，2021，56（1）：71-79.

[15] YANG L， LIU Y， XIANG Y，et al. Establishment and optimization of a shoot tip-based genetic transformation system for foxtail millet[J]. Chinese Bulletin of Botany，2021，56（1）：71-79.

[16] SAHA P， BLUMWALD E. Spike-dip transformation of Setaria viridis [J]. The Plant Journal：for Cell and Molecular Biology，2016，86（1）：89-101.

[17] 赵辉，郭静远，孔华，等. 狗尾草胚性愈伤高效农杆菌转化体系的建立[J]. 热带作物学报，2017，38（6）：1106-1112.

[18] ZHAO H， GUO J Y， KONG H，et al. The establishment of high efficient Agrobacterium transformation system of Setaria viridis [J]. Chinese Journal of Tropical Crops，2017，38（6）：1106-1112.

[19] PURANIK S， SAHU P P， MANDAL S N，et al. Comprehensive genome-wide survey，genomic constitution and expression profiling of the NAC transcription factor family in foxtail millet（Setaria italica L.）[J]. PLoS One，2013，8（5）：e64594.

[20] MISHRA A K， MUTHAMILARASAN M， KHAN Y，et al. Genome-wide investigation and expression analyses of WD40 protein family in the model plant foxtail millet（Setaria italica L.）[J]. PLoS One，2014，9（1）：e86852.

[21] MUTHAMILARASAN M， BONTHALA V S， MISHRA A K，et al. C₂H₂ type of zinc finger transcription factors in foxtail millet define response to abiotic stresses[J]. Functional & Integrative Genomics，2014，14（3）：531-543.

[22] MUTHAMILARASAN M， KHANDELWAL R， YADAV C B，et al. Identification and molecular characterization of MYB transcription factor superfamily in C4 model plant foxtail millet（Setaria italica L.）[J]. PLoS One，2014，9（10）：e109920.

[23] FENG Z J， HE G H， ZHENG W J，et al. Foxtail millet NF-Y families：genome-wide survey and evolution analyses identified two functional genes important in abiotic stresses[J]. Frontiers in Plant Science，2015，6：1142.

[24] CHAI W B， SI W N， JI W，et al. Genome-wide investigation and expression profiling of HD-zip transcription factors in foxtail millet（Setaria italica L.）[J]. BioMed Research International，2018，2018：8457614.

[25] ZHANG L， LIU B L， ZHENG G W，et al. Genome-wide characterization of the SiDof gene family in foxtail millet（Setaria italica）[J]. Biosystems，2017，151：27-33.

[26] ZHAO W， ZHANG L L， XU Z S，et al. Genome-wide analysis of MADS-box genes in foxtail millet（Setaria italica L.） and functional assessment of the role of SiMADS51 in the drought stress response[J]. Frontiers in Plant Science，2021，12：659474.

[27] LATA C R， MISHRA A K， MUTHAMILARASAN M，et al. Genome-wide investigation and expression profiling of AP2/ERF transcription factor superfamily in foxtail millet（Setaria italica L.）[J]. PLoS One，2014，9（11）：e113092.

[28] LI Y T， YU S M， ZHANG Q Y，et al. Genome-wide identification and characterization of the CCT gene family in foxtail millet（Setaria italica） response to diurnal rhythm and abiotic stress[J]. Genes，2022，13（10）：1829.

[29] GE L H， DOU Y N， LI M M，et al. SiMYB3 in foxtail millet（Setaria italica） confers tolerance to low-nitrogen stress by regulating root growth in transgenic plants[J]. International Journal of Molecular Sciences，2019，20（22）：5741.

[30] YI F， HUO M Y， LI J R，et al. Time-series transcriptomics reveals a drought-responsive temporal network and crosstalk between drought stress and the circadian clock in foxtail millet[J]. The Plant Journal：for Cell and Molecular Biology，2022，110（4）：1213-1228.

[31] WANG T， SONG H， LI P T，et al. Transcriptome analysis provides insights into grain filling in foxtail millet（Setaria italica L.）[J]. International Journal of Molecular Sciences，2020，21（14）：5031.

[32] LI X R， WANG L L， LI W D，et al. Genome-wide identification and expression profiling of cytochrome P450 monooxygenase superfamily in foxtail millet[J]. International Journal of Molecular Sciences，2023，24（13）：11053.

[33] WANG L L， FU H B， ZHAO J，et al. Genome-wide identification and expression profiling of glutathione S-transferase gene family in foxtail millet（Setaria italica L.）[J]. Plants，2023，12（5）：1138.

[34] CHANG O C， LIN W Y. Variation of growth and transcriptome responses to arbuscular mycorrhizal symbiosis in different foxtail millet lines[J]. Botanical Studies，2023，64（1）：16.

[35] LI S Y， CUI Y J， LIU D，et al. Transcriptome analysis and characterization of genes associated to leaf tannin content in foxtail millet（Setaria italica（L.） P.Beauv.）[J]. BMC Genomics，2022，23（1）：512.

[36] LIU D， CUI Y J， ZHAO Z L，et al. Transcriptome analysis and mining of genes related to shade tolerance in foxtail millet（Setaria italica（L.） P.Beauv.）[J]. Royal Society Open Science，2022，9（10）：220953.

[37] ZHAO Y F， LI Y Q， ZHEN X X，et al. Uncovering the mechanism of anthocyanin accumulation in a purple-leaved variety of foxtail millet（Setaria italica） by transcriptome analysis[J]. PeerJ，2022，10：e14099.

[38] LI H， HAN S L， HUO Y Q，et al. Comparative metabolomic and transcriptomic analysis reveals a coexpression network of the carotenoid metabolism pathway in the panicle of Setaria italica [J]. BMC Plant Biology，2022，22（1）：105.

[39] ZHANG H Y， WANG J Y， ZHAO J，et al. Integrated lipidomic and transcriptomic analysis reveals lipid metabolism in foxtail millet（Setaria italica）[J]. Frontiers in Genetics，2021，12：758003.

[40] YUAN Y H， LIU C J， ZHAO G，et al. Transcriptome analysis reveals the mechanism associated with dynamic changes in fatty acid and phytosterol content in foxtail millet（Setaria italica） during seed development[J]. Food Research International，2021，145：110429.

[41] LI R J， HAN Y Q， ZHANG Q，et al. Transcriptome profiling analysis reveals co-regulation of hormone pathways in foxtail millet during Sclerospora graminicola infection[J]. International Journal of Molecular Sciences，2020，21（4）：1226.

[42] LI P， LI B Y， SENEWEERA S，et al. Photosynthesis and yield response to elevated CO₂，C₄ plant foxtail millet behaves similarly to C₃ species[J]. Plant Science，2019，285：239-247.

[43] CAO X Y， HU L Q， CHEN X Y，et al. Genome-wide analysis and identification of the low potassium stress responsive gene SiMYB3 in foxtail millet（Setariaitalica L.）[J]. BMC Genomics，2019，20（1）：136.

[44] YI F， XIE S J， LIU Y W，et al. Genome-wide characterization of microRNA in foxtail millet（Setaria italica）[J]. BMC Plant Biology，2013，13：212.

[45] QI X， XIE S J， LIU Y W，et al. Genome-wide annotation of genes and noncoding RNAs of foxtail millet in response to simulated drought stress by deep sequencing[J]. Plant Molecular Biology，2013，83（4）：459-473.

[46] YADAV A， KHAN Y， PRASAD M. Dehydration-responsive miRNAs in foxtail millet：genome-wide identification，characterization and expression profiling[J]. Planta，2016，243（3）：749-766.

[47] DENG Y， ZHANG H L， WANG H L，et al. The construction and exploration of a comprehensive microRNA centered regulatory network in foxtail millet（Setaria italica L.）[J]. Frontiers in Plant Science，2022，13：848474.

[48] HOU S Y， MEN Y H， ZHANG Y J，et al. Role of miRNAs in regulation of SA-mediated upregulation of genes involved in folate and methionine metabolism in foxtail millet[J]. Frontiers in Plant Science，2022，13：1023764.

[49] WANG T， SONG H， WEI Y Y，et al. High throughput deep sequencing elucidates the important role of lncRNAs in foxtail millet response to herbicides[J]. Genomics，2020，112（6）：4463-4473.

[50] LIU Z， FAN H， MA Z Y. Comparison of SWEET gene family between maize and foxtail millet through genomic，transcriptomic，and proteomic analyses[J]. The Plant Genome，2022，15（3）：e20226.

[51] YANG Y， CHAI Y M， LIU J Y，et al. Amino acid transporter（AAT） gene family in foxtail millet（Setaria italica L.）：widespread family expansion，functional differentiation，roles in quality formation and response to abiotic stresses[J]. BMC Genomics，2021，22（1）：519.

[52] YU T F， ZHAO W Y， FU J D，et al. Genome-wide analysis of CDPK family in foxtail millet and determination of SiCDPK24 functions in drought stress[J]. Frontiers in Plant Science，2018，9：651.

[53] YADAV C B， MUTHAMILARASAN M， DANGI A，et al. Comprehensive analysis of SET domain gene family in foxtail millet identifies the putative role of SiSET14 in abiotic stress tolerance[J]. Scientific Reports，2016，6：32621.

[54] XING G F， JIN M S， QU R F，et al. Genome-wide investigation of histone acetyltransferase gene family and its responses to biotic and abiotic stress in foxtail millet（Setaria italica（L.）P.Beauv）[J]. BMC Plant Biology，2022，22（1）：292.

[55] JIA G Q， HUANG X H， ZHI H，et al. A haplotype map of genomic variations and genome-wide association studies of agronomic traits in foxtail millet（Setaria italica）[J]. Nature Genetics，2013，45：957-961.

[56] LI C C， WANG G P， LI H Q，et al. High-depth resequencing of 312 accessions reveals the local adaptation of foxtail millet[J]. Theoretical and Applied Genetics，2021，134（5）：1303-1317.

[57] JAISWAL V， GUPTA S， GAHLAUT V，et al. Genome-wide association study of major agronomic traits in foxtail millet（Setaria italica L.） using ddRAD sequencin g[J].Scientific Reports，2019，9（1）：5020.

[58] LIU X D， YANG Y， HOU S Y，et al. The integration of genome-wide association study and homology analysis to explore the genomic regions and candidate genes for panicle-related traits in foxtail millet[J]. International Journal of Molecular Sciences，2022，23（23）：14735.

[59] LI X K， GAO J H， SONG J Y，et al. Multi-omics analyses of 398 foxtail millet accessions reveal genomic regions associated with domestication，metabolite traits，and anti-inflammatory effects[J]. Molecular Plant，2022，15（8）：1367-1383.

[60] DOUST A N， DEVOS K M， GADBERRY M D，et al. The genetic basis for inflorescence variation between foxtail and green millet（Poaceae）[J]. Genetics，2005，169（3）：1659-1672.

[61] FANG X M， DONG K J， WANG X Q，et al. A high density genetic map and QTL for agronomic and yield traits in foxtail millet（Setaria italica（L.） P.Beauv）[J]. BMC Genomics，2016，17：336.

[62] WANG Z L， WANG J， PENG J X，et al. QTL mapping for 11 agronomic traits based on a genome-wide Bin-map in a large F₂ population of foxtail millet（Setaria italica（L.） P.Beauv）[J]. Molecular Breeding，2019，39（2）：18.

[63] ZHI H， HE Q， TANG S，et al. Genetic control and phenotypic characterization of panicle architecture and grain yield-related traits in foxtail millet（Setaria italica）[J]. Theoretical and Applied Genetics，2021，134（9）：3023-3036.

[64] 王育选，张梦超，常丽君，等. 谷子MADS-box基因家族的鉴定和表达分析[J]. 生物信息学，2023，21（1）：51-59.

[65] WANG Y X， ZHANG M C， CHANG L J，et al. Genome-wide identification and expression analysis of MADS-box genes in foxtail millet[J]. Chinese Journal of Bioinformatics，2023，21（1）：51-59.

[66] LAI D L， FAN Y， XUE G X，et al. Genome-wide identification and characterization of the SPL gene family and its expression in the various developmental stages and stress conditions in foxtail millet（Setaria italica）[J]. BMC Genomics，2022，23（1）：389.

[67] HEISLER M G， ATKINSON A， BYLSTRA Y H，et al. SPATULA，a gene that controls development of carpel margin tissues in Arabidopsis，encodes a bHLH protein[J]. Development，2001，128（7）：1089-1098.

[68] ALVAREZ J， SMYTH D R. CRABS CLAW and SPATULA，two Arabidopsis genes that control carpel development in parallel with AGAMOUS [J]．Development，1999，126（11）：2377-2386.

[69] FAN Y， WEI X B， LAI D L，et al. Genome-wide investigation of the GRAS transcription factor family in foxtail millet（Setaria italica L.）[J]. BMC Plant Biology，2021，21（1）：508.

[70] CHOI D， KIM J H， KENDE H. Whole genome analysis of the OsGRF gene family encoding plant-specific putative transcription activators in rice（Oryza sativa L.）[J]. Plant and Cell Physiology，2004，45（7）：897-904.

[71] ZHANG D F， LI B， JIA G Q，et al. Isolation and characterization of genes encoding GRF transcription factors and GIF transcriptional coactivators in maize（Zea mays L.）[J]. Plant Science，2008，175（6）：809-817.

[72] WANG F D， QIU N W， DING Q，et al. Genome-wide identification and analysis of the growth-regulating factor family in Chinese cabbage（Brassica rapa L.ssp.pekinensis）[J]. BMC Genomics，2014，15（1）：807.

[73] ZHANG J F， LI Z F， JIN J J，et al. Genome-wide identification and analysis of the growth-regulating factor family in tobacco（Nicotiana tabacum）[J]. Gene，2018，639：117-127.

[74] 雷彪，曲瑞芳，任超，等. 谷子LBD基因家族的鉴定及其对逆境胁迫的响应[J]. 植物生理学报，2023，59（3）：527-542.

[75] LEI B， QU R F， REN C，et al. Identification of SiLBDs gene family in foxtail millet and its participation in stress response[J]. Plant Physiology Journal，2023，59（3）：527-542.

[76] LI C C， MA J， WANG G P，et al. Exploring the SiCCT gene family and its role in heading date in foxtail millet[J]. Frontiers in Plant Science，2022，13：863298.

[77] 朱垠豪，李晨，葛星辰，等. 谷子KCS基因家族鉴定及其对干旱胁迫的响应[J]. 山西农业科学，2023，51（6）：599-609.

[78] ZHU Y H， LI C， GE X C，et al. Identification and response to drought stress of KCS gene family in foxtail millet[J]. Journal of Shanxi Agricultural Sciences，2023，51（6）：599-609.

[79] 晋敏姗，曲瑞芳，李红英，等. 谷子糖转运蛋白基因SiSTPs的鉴定及其参与谷子抗逆胁迫响应的研究[J]. 作物学报，2022，48（4）：825-839.

[80] JIN M S， QU R F， LI H Y，et al. Identification of sugar transporter gene family SiSTPs in foxtail millet and its participation in stress response[J]. Acta Agronomica Sinica，2022，48（4）：825-839.

[81] WANG H， WANG H Y. The miR156/SPL module，a regulatory hub and versatile toolbox，gears up crops for enhanced agronomic traits[J]. Molecular Plant，2015，8（5）：677-688.

[82] 邢恒荣. 谷子WRKY基因家族鉴定及SiWRKY89基因功能鉴定[D]. 太谷：山西农业大学，2019.

[83] XING H R. Identification of WRKY gene family and functional identification of SiWRKY89 gene in foxtail millet[D]. Taigu：Shanxi Agricultural University，2019.

[84] 尹泽群，耿菁，刘莹，等. 谷子CCCH全基因家族鉴定及非生物胁迫下的表达分析[J/OL]. 分子植物育种：1-13[2023-07-26].

[85] YIN Z Q， GENG J， LIU Y，et al. Genome-wide identification of CCCH and expression analysis under abiotic stress in foxtail millet[J/OL]. Molecular Plant Breeding：1-13[2023-07-26].

[86] 张浩林. 谷子生长调节因子SiGRFs基因家族鉴定及功能初步分析[D]. 上海：上海师范大学，2022[2023-07-26].

[87] ZHANG H L. Identification and functional analysis of millet growth regulator SiGRFs gene family[D]. Shanghai：Shanghai Normal University，2022.

[88] 赵娟莹. 谷子、大豆ANK家族基因鉴定及SiANK37和GmANK114的抗旱耐盐功能解析[D]. 杨凌：西北农林科技大学，2022.

[89] ZHAO J Y. Identification of ANK family genes in millet and soybean and analysis of drought resistance and salt tolerance of SiANK37 and GmANK114[D]. Yangling：Northwest A & F University，2022.

[90] 魏萌涵，解慧芳，赵二源，等. 谷子APX基因家族的全基因组鉴定与表达分析[J/OL]. 分子植物育种:1-11[2021-12-13].

[91] WEI M H， XIE H F， ZHAO E Y，et al. Genome-wide identification and expression analysis of APX gene family in foxtail millet[J/OL]. Molecular Plant Breeding：1-11[2021-12-13].

[92] 石尧. 谷子SULTR基因家族鉴定及SiSULTRC3硒转运功能分析[D]. 太谷：山西农业大学，2022.

[93] SHI Y. Identification of SULTR gene family and analysis of transport function of SiSULTRC3 [D]. Taigu：Shanxi Agricultural University，2022.

[94] 宋敏. 高叶酸谷子种质筛选与SiFPGS基因特征、表达及遗传转化研究[D]. 太谷：山西农业大学，2018.

[95] SONG M. Screening of millet germplasm with high folic acid and study on the characteristics，expression and genetic transformation of SiFPGS gene[D]. Taigu：Shanxi Agricultural University，2018.

[96] ZHANG B， GUO Y， WANG H，et al. Identification and characterization of shaker K⁺ channel gene family in foxtail millet（Setaria italica） and their role in stress response[J]. Frontiers in Plant Science，2022，13：907635.

[97] ALAGARASAN G， DUBEY M， ASWATHY K S，et al. Genome wide identification of orthologous ZIP genes associated with zinc and iron translocation in Setaria italica [J]. Frontiers in Plant Science，2017，8：775.

[98] ROCH G V， MAHARAJAN T， AJEESH KRISHNA T P，et al. Expression of PHT1 family transporter genes contributes for low phosphate stress tolerance in foxtail millet（Setaria italica） genotypes[J]. Planta，2020，252（6）：98.

[99] 李明明，晋敏姗，胡海斌，等. 谷子MYB-CC基因家族的鉴定与表达分析[J]. 激光生物学报，2021，30（3）：259-269.

[100] LI M M， JIN M S， HU H B，et al. Identification and expression analysis of MYB-CC gene family in foxtail millet[J]. Acta Laser Biology Sinica，2021，30（3）：259-269.

[101] 张敏. 谷子苗期响应低氮胁迫生理及转录组特征分析[D]. 太谷：山西农业大学，2021.

[102] ZHANG M. Physiological and transcriptome characteristics of foxtail millet seedlings in response to low nitrogen stress[D]. Taigu：Shanxi Agricultural University，2021.

[103] DECHORGNAT J， NGUYEN C T， ARMENGAUD P，et al. From the soil to the seeds：the long journey of nitrate in plants[J]. Journal of Experimental Botany，2011，62（4）：1349-1359.

基金资助

山西省基础研究计划项目(20210302124367)