《International Journal of Molecular Sciences》：Genome-Wide Analysis of FAR-RELATED SEQUENCES (FRS) Genes Related to Light Response in Soybean (Glycine max)
Sujun Ye,
Lixin Ma,
Yinhua Lv,
Wenmin Lin,
Yang Tang and
Xiaoya Lin
编辑推荐：
本文通过全基因组分析，在大豆中鉴定出49个GmFRS基因，并系统研究了其在光信号、胁迫适应和发育调控中的潜在功能。该研究揭示了FRS家族在植物感知和响应光环境（尤其是远红光）中的核心作用，为理解大豆等作物在复杂环境下的适应机制及分子育种提供了关键靶点。
FAR1-RELATED SEQUENCES (FRS)基因家族是植物特异性转录因子家族，在光信号转导、胁迫适应和发育调控中发挥关键作用，直接影响作物生长和产量。本研究旨在对大豆(Glycine max)中的FRS基因家族进行全基因组综合分析，以填补相关研究空白，并为大豆育种和产量提升提供理论依据。
鉴定与物理化学分析
研究利用BLAST比对和结构域验证，从大豆基因组中准确鉴定出49个高置信度的GmFRS基因。这些基因编码的蛋白质在氨基酸数量、分子量和理论等电点(pI)等方面表现出多样性。亚细胞定位预测表明，除GmFRS45被预测定位于质膜外，其余成员均主要定位于细胞核，符合其作为转录因子的功能。
系统发育与进化分析
通过构建大豆、拟南芥(Arabidopsis thaliana)、番茄(Solanum lycopersicum)和玉米(Zea mays) FRS蛋白的系统发育树，108个FRS基因被划分为七个亚组(I–VII)。值得注意的是，亚组VII缺乏拟南芥同源物，但包含两个大豆成员(GmFRS8, GmFRS22)，它们与番茄和玉米的蛋白聚类，提示可能存在功能保守性。系统发育分析还表明，拟南芥中的关键光信号转导因子AtFHY3和AtFAR1分别与大豆中的GmFRS7/GmFRS12和GmFRS23/GmFRS49亲缘关系最近，暗示它们可能具有相似功能。
染色体分布与复制事件
49个GmFRS基因不均匀地分布在17条大豆染色体上，其中15号染色体上成员最多(9个)。共线性分析识别出15对同线基因对，所有配对均为片段复制/全基因组复制(WGD)事件，表明GmFRS家族的扩张主要由此驱动。与非同线基因复制的结合，共同促进了该家族的扩张。Ka/Ks比值分析显示，所有基因对的比值均显著小于1，表明GmFRS家族在进化过程中受到强烈的纯化选择，这是功能重要性基因的典型特征。
保守结构域与基因结构
对GmFRS蛋白的保守基序和结构域分析揭示了其结构多样性。大多数成员包含三个核心结构域：N端的FAR1 DNA结合域、中央的MULE转座酶结构域和C端的SWIM锌指结构域。然而，也发现了多种结构变异类型，例如亚组VI的成员是仅包含FAR1 DNA结合域的截短蛋白，而GmFRS5和GmFRS46则完全缺失该结构域。基因结构分析显示，同一亚组内的基因具有相对保守的外显子-内含子结构，其复杂度与保守蛋白结构域的组成显著相关。例如，包含三个完整结构域的亚组I、II和IV成员通常具有一个较长的外显子。
蛋白结构特征
二级结构预测表明，GmFRS蛋白主要由α-螺旋和无规则卷曲构成。通过AlphaFold 3预测了所有49个GmFRS蛋白的三维结构，发现同一亚组内的蛋白具有相似的整体折叠构象。对关键成员(GmFRS7, GmFRS12, GmFRS23, GmFRS49)的结构解析进一步确认了保守的FAR1 DNA结合域（黄色）和MULE转座酶结构域（红色）的存在。
启动子顺式作用元件分析
对GmFRS基因起始密码子上游2000 bp的启动子序列进行分析，共鉴定出1134个顺式作用元件。这些元件主要与光响应(44.8%)、植物激素(29.4%)、胁迫响应(18.0%)以及植物生长发育(7.8%)相关。其中最丰富的元件是光响应元件Box-4。此外，与脱落酸(ABA)响应相关的ABRE元件、与茉莉酸(JA)响应相关的CGTCA/TGACG motif、与厌氧诱导相关的ARE元件以及与干旱响应相关的MBS元件等也广泛存在。这复杂的顺式调控网络表明GmFRS基因能够整合多种信号来调控其表达。
亚细胞定位验证
通过农杆菌(Agrobacterium tumefaciens)介导的瞬时转化实验，在本氏烟草(Nicotiana benthamiana)叶片表皮细胞中观察了GmFRS7、GmFRS12、GmFRS23和GmFRS49与绿色荧光蛋白(GFP)融合蛋白的亚细胞定位。荧光显微镜观察显示，这四个融合蛋白均主要定位于细胞核，与生物信息学预测及其拟南芥同源物AtFHY3/AtFAR1的已知定位一致，支持了它们作为转录因子的功能预测。
不同组织中的表达模式
基于大豆表达图谱(Soybean Expression Atlas)的RNA-seq数据分析显示，GmFRS基因在不同器官中呈现差异表达。总体而言，GmFRS转录本在花中最为丰富(占总量的32.7%)，其次是在茎中(26.5%)。有8个基因(GmFRS7, GmFRS12, GmFRS13, GmFRS20, GmFRS26, GmFRS28, GmFRS32, GmFRS48)在所有检测器官中均表现出相对较高的转录丰度，提示它们可能在多种生理过程中发挥功能。而另有7个基因(GmFRS10, GmFRS11, GmFRS21, GmFRS27, GmFRS29, GmFRS43, GmFRS45)的表达量几乎检测不到。
光响应下的转录表达
鉴于启动子中含有丰富的光响应元件，研究进一步通过qRT-PCR分析了GmFRS基因在不同红/远红光比例(R:FR)光照下的表达响应。在长日照条件下，将大豆植株暴露于高R:FR (R:FR = 4.5)和低R:FR (R:FR = 0.45)的光照中，并在光照处理前4小时(-4 h)、光照开始时(0 h)和光照4小时后(4 h)采集叶片样品。结果表明，所有检测的42个GmFRS基因均对富远红光条件有响应，其表达水平在远红光照射后上升。大多数基因(35/42)在持续远红光照射下也能被激活。另一方面，富红光条件也能激活几乎所有的FRS基因（有5个例外）。然而，对于大多数基因而言，红光的激活效应不如远红光显著（GmFRS20除外）。其中，GmFRS7, GmFRS15, GmFRS24, GmFRS28, GmFRS36和GmFRS49对远红光表现出明显的偏好性敏感性，暗示这6个基因可能在避荫反应中起主要作用。
结论与展望
本研究首次对大豆FRS基因家族进行了全基因组范围内的系统鉴定和综合分析，成功鉴定出49个GmFRS基因。研究揭示了该家族在进化、结构和表达上的特征，特别是其与光信号（尤其是远红光）响应的紧密关联。启动子分析提示GmFRS基因可能整合光、激素和胁迫等多种信号。组织特异性表达和光诱导表达谱为进一步阐明GmFRS基因在大豆生长发育和环境适应中的功能提供了重要线索。该研究为深入解析大豆FRS家族在光信号转导、发育调控及胁迫适应中的分子机制奠定了坚实基础，也为利用关键GmFRS基因进行大豆分子育种改良提供了潜在靶点。未来研究应通过CRISPR/Cas9基因敲除、DNA结合特异性分析及蛋白互作研究等手段，对核心候选基因（如GmFRS7, GmFRS12, GmFRS23, GmFRS49）的功能进行验证。
