紫花苜蓿MAG2基因家族的全基因组鉴定与特征分析

《Genes》：Genome-Wide Identification and Characterization of the MAG2 Gene Family in Medicago sativa

　　

背景/目的：囊泡运输介导蛋白质及其他细胞组分在细胞内细胞器之间的转运。MAG2复合体作为介导内质网–高尔基体（ER–Golgi）逆向囊泡运输的关键拴系因子，已被证实参与植物生长、发育及胁迫应答。方法：本研究利用中牧4号参考基因组，通过隐马尔可夫模型搜索（HMM）和BLASTP分析，在紫花苜蓿（*Medicago sativa*）中鉴定到7个MAG2家族基因。结果：系统发育分析将4个基因归入MAG2亚家族，3个基因归入MAG2L亚家族。染色体分布和共线性分析表明，片段重复可能促进了苜蓿中MAG2家族的扩张。基因结构与保守基序分析显示，所鉴定成员之间具有高度保守性，其中基序10（Motif 10）特异存在于MAG2亚家族中。亚细胞定位预测分析预测4个MsMAG2蛋白定位于叶绿体，2个定位于细胞核，1个定位于过氧化物酶体。启动子序列分析鉴定出大量与激素、光照及环境胁迫应答相关的顺式调控元件，提示*MsMAG2*基因可能参与植物发育与环境适应。此外，全基因组关联分析鉴定出与分蘖可能相关的*MsMAG2L2*。结论：这些结果为苜蓿MAG2基因家族提供了全面的全基因组特征分析，并为进一步研究*MsMAG2L2*在分枝发育中的功能关联奠定了基础。

论文解读：紫花苜蓿MAG2基因家族的全基因组鉴定与功能关联分析

**1. 研究背景与问题**
真核细胞含有众多由磷脂双分子层包裹的细胞器，它们为特定生化反应提供了相对独立的反应空间。蛋白质、脂质等细胞组分的交换主要依赖囊泡运输，这一过程维持了细胞器身份并协调细胞内物质分配。囊泡运输包括出芽、定向移动、拴系以及膜融合等多个步骤，其中拴系因子（tethering factor）介导囊泡与靶膜的初始接触。MAG2复合体是定位于内质网（ER）的多亚基拴系复合体，调控内质网–高尔基体（ER–Golgi）之间的逆向囊泡运输。在拟南芥中，MAG2功能缺失会导致储藏蛋白前体在内质网中积累，并形成MAG2小体。已有研究表明，MAG2及相关蛋白（如MAL、MIP1/2/3）参与植物生长发育、激素信号响应及非生物胁迫应答。然而，在重要豆科牧草紫花苜蓿（*Medicago sativa*）中，MAG2基因家族的系统鉴定和功能研究尚未开展。紫花苜蓿被称为“牧草之王”，具有高蛋白和高营养含量，是全球广泛栽培的多年生豆科植物。因此，鉴定苜蓿MAG2家族并解析其进化特征与潜在功能，对理解该家族在植物发育及环境适应中的作用具有重要意义。

**2. 主要技术方法**
研究人员利用已发表的中牧4号参考基因组及FigShare平台下载的基因组注释文件。通过PFAM数据库获得MAG2保守结构域（PF04437）的隐马尔可夫模型，使用TBtools中的HMM Search（E-value阈值1×10^-10）和BLASTP比对筛选候选基因，并手动过滤保留覆盖度>50%且具有完整开放阅读框的序列。系统发育分析采用MEGA7.0进行多序列比对，基于最大似然法（ML）构建进化树（bootstrap=1000）。保守基序分析使用MEME工具，基因结构分析结合GFF文件在TBtools中可视化。利用MCScanX进行共线性分析。启动子顺式元件分析提取转录起始位点上游2000bp序列。全基因组关联分析（GWAS）使用100份苜蓿种质资源组成的关联群体，经VCFtools质控后获得5,653,034个高质量SNP，采用一般线性模型（GLM）和混合线性模型（MLM）进行关联分析，显著阈值设为LOD≥5。蛋白互作网络预测使用STRING数据库，以蒺藜苜蓿（*Medicago truncatula*）为参考物种。

**3. 研究结果**
**3.1 MAG2基因家族鉴定**
研究人员在苜蓿中鉴定到7个MAG2基因，命名为*MsMAG2-1*至*MsMAG2-4*、*MsMAG2L1*至*MsMAG2L3*。亚细胞定位预测显示：4个蛋白（MsMAG2-1/2/3/4）定位于叶绿体，2个蛋白（MsMAG2L1、MsMAG2L3）定位于细胞核，1个蛋白（MsMAG2L2）定位于过氧化物酶体。蛋白长度从505个氨基酸（MsMAG2L2）到839个氨基酸（MsMAG2L3）不等；分子量范围为59,022.97 Da至96,895.75 Da；等电点范围为5.52至6.66，均为酸性蛋白。染色体定位显示，7个基因分布在第3号和第7号染色体的不同单倍型上，其中MAG2亚家族全部位于第7号染色体（chr7_1、chr7_2、chr7_3、chr7_4），MAG2L亚家族全部位于第3号染色体（chr3_1、chr3_2、chr3_3），这种分布与苜蓿同源四倍体基因组特征一致。

**3.2 MAG2蛋白的系统发育结构**
基于苜蓿和拟南芥MAG2蛋白构建的系统发育树显示，7个MsMAG2蛋白分为两大类：CLASS I（MAG2亚家族）和CLASS II（MAG2L亚家族）。与拟南芥仅含有1个AtMAG2和1个AtMAG2L相比，苜蓿MAG2亚家族有4个成员（57.1%），MAG2L亚家族有3个成员（42.9%）。每个单倍型仅含1个成员，提示片段重复可能是家族扩张的主要机制。

**3.3 保守基序与基因结构分析**
MEME预测共鉴定出10个保守基序。所有MAG2蛋白均含有Motif1、2、3、4、6、7、8和9，表明这些基序可能在MAG2蛋白的结构和功能中发挥关键作用。Motif10仅存在于MAG2亚家族的所有成员中，而MAG2L亚家族均不含该基序，提示Motif10与亚家族功能分化相关。Motif5在MAG2L亚家族中出现不规律，存在跨亚家族变异。

**3.4 基因共线性分析**
种内共线性分析检测到7对共线性基因关系，包括*MsMAG2-1*与*MsMAG2-4*、*MsMAG2-1*与*MsMAG2-2*、*MsMAG2-4*与*MsMAG2-2*，以及*MsMAG2L2*与*MsMAG2L3*、*MsMAG2L2*与*MsMAG2L1*等。值得注意的是，共线性只发生在同一亚家族内部，不同亚家族之间未检测到共线性，且两个亚家族内所有基因间均存在共线性，表明MAG2家族在进化过程中高度保守，与基序分析结果一致。

**3.5 顺式作用元件分析**
对7个*MsMAG2*基因启动子区（起始密码子上游2000 bp）的分析发现，它们含有F-box元件以及大量MYC和MYB转录因子结合位点。各基因所含功能元件的数量从29个（*MsMAG2L1*）到53个（*MsMAG2-4*）不等。总共鉴定出144个与光响应、激素响应及胁迫响应相关的元件，其中激素响应元件最多（86个），其次为光响应元件（61个），胁迫响应元件最少（17个）。在激素响应元件中，茉莉酸响应元件（45个）最多，脱落酸响应元件（28个）次之，水杨酸响应元件（7个）和生长素响应元件（5个）较少，赤霉素响应元件最少（1个）。胁迫相关元件中包含8个冷胁迫响应元件和6个干旱响应元件。这些结果提示*MsMAG2*基因可能主要受茉莉酸和脱落酸调控，部分基因可能参与抗寒和抗旱调节。

**3.6 全基因组关联分析**
利用100份苜蓿种质资源的分蘖性状数据进行GWAS，在染色体chr3_3:34372876位置检测到一个显著SNP。基因注释和比对分析表明该SNP位于*MsMAG2L2*基因（Msa0449710）内。该结果提示*MsMAG2L2*与分蘖相关，但并未证明其直接调控作用。

**3.7 MsMAG2L2蛋白互作网络预测**
基于蒺藜苜蓿同源蛋白A0A072V5P0的STRING预测表明，MsMAG2L2与10个蛋白存在潜在互作（p值<1.0×10^-16）。其中A0A072TMC2被注释为膜融合蛋白Use1，G7JT86被注释为推定的囊泡运输蛋白。这些互作蛋白的功能注释提示MsMAG2L2可能参与膜融合和囊泡运输相关过程。

**4. 讨论与结论**
研究结果表明，苜蓿MAG2家族相对拟南芥出现扩张，这与其同源四倍体基因组背景和进化过程中重复基因的保留有关。基因结构和基序的保守性反映了MAG2蛋白在囊泡运输中的功能保守性。然而，预测的亚细胞定位（叶绿体、细胞核、过氧化物酶体）与拟南芥MAG2复合体的ER–Golgi定位不一致，这可能反映了基于序列预测算法的局限性，需要荧光蛋白融合瞬时表达等实验验证。启动子顺式元件分析提示*MsMAG2*基因可能响应多种发育和环境信号，但需注意顺式元件的存在仅表示潜在转录响应能力，并非实际调控证据。GWAS将*MsMAG2L2*定位为分蘖相关候选基因，结合其蛋白互作网络中含膜融合和囊泡运输相关蛋白，推测其可能通过影响发育调控因子的运输或定位而参与分枝发育。但上述推断均为初步结果，需通过表达分析、亚细胞定位、基因编辑或过表达等实验加以验证。总体而言，苜蓿MAG2家族在进化上保守且经历了扩张，基因结构、启动子组成和GWAS结果共同表明*MsMAG2*基因不仅可能参与保守的囊泡运输过程，还可能参与分蘖等农艺性状的调控。本研究为深入解析苜蓿MAG2家族基因的生物学功能提供了有价值的框架和遗传资源，相关成果发表在《Genes》期刊。

**结论翻译**：本研究在紫花苜蓿中鉴定了MAG2基因家族，并将7个成员划分为MAG2和MAG2L两个亚家族。系统发育和染色体分析显示这些基因在基因组中分布不均，片段重复可能驱动了家族扩张。保守基序和基因结构分析表明MAG2家族在紫花苜蓿中进化保守，Motif 10是MAG2亚家族的特征基序。值得注意的是，GWAS分析鉴定出*MsMAG2L2*为与分蘖性状显著相关的候选基因，提示其可能参与紫花苜蓿分枝发育的调控。然而，本研究中的基因分析仍属理论性工作，尚需进一步的实验验证。