《Discover Plants》：In silico characterization of C3H transcription factors reveals key regulators of fungal disease responses in apple (Malus domestica)
编辑推荐：
苹果（Malus domestica）是一种广泛栽培的果树作物，具有遗传和表型多样性，并含有因促进健康价值而受到重视的生物活性化合物。尽管其是重要果树作物，但与生物胁迫（biotic stress）耐受相关信号通路的转录因子（transcription fac
苹果（Malus domestica）是一种广泛栽培的果树作物，具有遗传和表型多样性，并含有因促进健康价值而受到重视的生物活性化合物。尽管其是重要果树作物，但与生物胁迫（biotic stress）耐受相关信号通路的转录因子（transcription factors, TFs）调控尚未被充分表征。C3H转录因子家族是积极参与植物发育和胁迫响应的重要转录因子家族之一。该家族已在多种植物物种中研究，但其在苹果响应生物胁迫中的作用仍表征不足。本研究共鉴定60个C3H成员，并依据其在全部17条染色体上的分布命名为MdC3H01至MdC3H60。随后，根据其与Arabidopsis thaliana、Prunus persica和Oryza sativa的同源性，将这些基因划分为不同组。研究人员与另外两个物种进行共线性分析（synteny analysis），以探讨其进化关系。在共线性分析中鉴定到23对片段重复（segmental duplication）基因对和2对串联重复（tandem duplication）基因对。基因表达分析（gene expression analysis）显示，MdC3H16、MdC3H37在响应Venturia inaequalis、Gymnosporangium yamadae和Alternaria alternata时高度诱导。这些基因作为协调对多种病原防御反应的核心调控因子。本研究广泛研究了C3H在植物发育及抗真菌病原中的作用。这些发现为通过分子育种培育具有改良生物胁迫抗性的抗逆品种提供了宝贵遗传资源。
苹果（Malus domestica）是蔷薇科重要果树，具有较高经济价值并含多糖、多酚、三萜和植物甾醇等生物活性化合物。真菌病害是苹果生产的重要威胁，Venturia inaequalis、Alternaria alternata和Gymnosporangium yamadae分别引起苹果黑星病、链格孢叶斑病和日本锈病，造成产量和品质损失。C3H型锌指蛋白属于转录因子（transcription factors, TFs；调控基因表达的蛋白），参与植物发育、信号转导及生物胁迫（biotic stress；病原等生物因素引起的胁迫）响应，但苹果中C3H家族响应生物胁迫的作用表征不足。因此，研究人员开展全基因组鉴定和转录响应分析。研究鉴定出60个非冗余MdC3H基因，命名为MdC3H01—MdC3H60，并与A. thaliana、P. persica和O. sativa同类基因进行同源性分组和共线性（synteny；保守基因排列关系）分析；表达分析显示MdC3H16、MdC3H37在三种病原诱导下高度上调，提示其作为协调多种病原防御反应的核心调控因子。该研究为抗病分子育种提供了遗传资源。
主要技术方法：研究人员从PlantTFDB和NCBI获取序列，BLAST确认C3H保守结构域；用TBtools-II、GSDS v2.0和MEME进行染色体定位、内含子–外显子及保守基序分析；MEGA11以邻接法构建系统发育树，iTOL可视化，NCBI CDD鉴定保守结构域；MCScanX分析共线性，PlantCARE鉴定启动子顺式元件；ShinyGO进行基因本体（Gene Ontology, GO；基因功能分类体系）富集，STRING构建蛋白–蛋白互作（protein–protein interaction, PPI；蛋白质间功能关联）网络；转录组数据来自NCBI SRA登录号SRP478269、SRP334696、SRP201829（对应Bioprojects PRJNA1053344、PRJNA758843、PRJNA549565），在Galaxy服务器以Bowtie2比对M. domestica参考基因组GCF_042453785.1，featureCounts计算FPKM（fragments per kilobase of exon model per million mapped reads；基因表达定量指标）。
3.1 C3H家族基因在Malus domestica中的鉴定与分布（Identification and distribution of C3H family genes in Malus domestica）：全基因组数据鉴定113个C3H基因，去除53个冗余后得60个非冗余MdC3H，命名MdC3H01—MdC3H60；MdC3H42最短（218 aa），MdC3H29最长（2,583 aa），平均597.6 aa；pI为4.073—11.16，均值7.5；分子量24.04—288 kDa，均值68.8 kDa。通过全基因组鉴定和理化性质分析得出结论。
3.2 MdC3H基因的基因组定位（Genomic localization of MdC3H genes）：60个MdC3H分布于全部17条染色体但不均匀；chr15最多（8个，13.3%），chr3为6个（10.0%）；chr5、chr8、chr10各1个（1.66%）；chr9、chr12、chr16各2个（3.33%）；chr2、chr4、chr7、chr11、chr13、chr14共24个（每条4个）。通过染色体定位分析得出结论。
3.3 比较系统发育与保守结构域分析（Comparative phylogeny of C3H gene family in M. domestica, A. thaliana, Prunus persica and O. sativa and conserved domain analysis of MdC3H genes）：224个来自四个物种的C3H蛋白构建无根邻接树，分为9组；组6最大（43个），组9最小（13个）；组4为37个；组7、8、3分别为24、15、17个；组1、2、5各25个。鉴定34个不同保守结构域；MdC3H14和MdC3H39无保守结构域；zf-CCCH（24个）和ZnF_C3H1（20个）最常见。通过系统发育和NCBI CDD分析得出结论。
3.4 MdC3H基因启动子顺式调控元件预测（Cis-regulatory elements prediction in MdC3H genes promoters）：检测到激素信号、光调控、非生物胁迫、核心启动子和防御相关元件；TATA/CAAT常见；ABRE在40个、Box 4在38个基因中；ARE和G-box各32个；CGTCA/TGACG MeJA响应元件在27个；G-Box、TCT-motif、GT1-motif分别为25、21、20个。提示MdC3H可能响应多种环境胁迫。通过PlantCARE分析得出结论。
3.5 基因结构与基序分析（Gene structure and motif analyses）：MdC3H21和MdC3H29分别含20和21个内含子、21和22个外显子；14个MdC3H基因无内含子，其余含1—11个内含子。MEME鉴定10个保守基序；motif1在49个、motif3在41个蛋白中最普遍，主要位于组6；motif2最长（50 aa），motif4最短（11 aa）。通过GSDS和MEME分析得出结论。
3.6 共线性与基因重复分析（Synteny and gene duplication analyses）：与O. sativa共线性区域少，与A. thaliana中等，与P. persica最强（同科）。25对重复基因中，2对为串联重复（MdC3H22/MdC3H23、MdC3H57/MdC3H58），23对为片段重复；chr15重复最多（7个），chr5和chr10无重复片段。通过MCScanX和Circos分析得出结论。
3.7 基因富集分析与蛋白–蛋白互作网络构建（Gene enrichment analysis and protein–protein interaction network construction）：GO富集显示参与mRNA结合、snoRNA剪接、锌指基序（CCCH型、Nab2型）及DNA/RNA结合；还富集木质素生物合成、春化响应和金属离子响应。PPI网络涉及12个MdC3H蛋白；MdC3H31与MdC3H18、MdC3H53等高度连接；MdC3H40与MdC3H58存在潜在连接。通过ShinyGO和STRING分析得出结论。
3.8 MdC3H基因对Venturia inaequalis、Alternaria alternata和Gymnosporangium yamadae的响应表达谱（Expression profiling of MdC3H genes in response to Venturia inaequalis, Alternaria alternata, and Gymnosporangium yamadae）：V. inaequalis感染下MdC3H06、MdC3H16、MdC3H37、MdC3H42（组1和4）高度上调，MdC3H03、MdC3H26、MdC3H39受抑制；MdC3H49、MdC3H59、MdC3H60早期稳定。A. alternata在6、18、24、48 hpi（hours post-inoculation，接种后小时）下MdC3H25最强上调，MdC3H16、MdC3H37、MdC3H46较低上调；MdC3H07、MdC3H10、MdC3H33、MdC3H56中性；MdC3H03、MdC3H26、MdC3H30、MdC3H59持续下调。G. yamadae在10 dpi（days post-inoculation，接种后天数）下MdC3H16、MdC3H25、MdC3H37、MdC3H42、MdC3H46上调，MdC3H03、MdC3H11、MdC3H26、MdC3H39、MdC3H60下调；30 dpi下MdC3H16、MdC3H25、MdC3H37、MdC3H42高度诱导，MdC3H20、MdC3H26、MdC3H39、MdC3H57、MdC3H60强抑制。通过转录组表达谱分析得出结论。
讨论部分总结：C3H蛋白参与植物生长、信号转导及非生物和生物胁迫响应。研究人员比较了不同物种C3H数量，并对MdC3H25、MdC3H23、MdC3H60、MdC3H03和MdC3H18等功能注释进行讨论；chr15基因数最多与Chen et al.报道一致；系统发育分组、顺式元件、基因结构、保守基序、基因重复、共线性、GO富集和PPI网络均得到讨论。无内含子基因可能转录更快，有利于胁迫下快速诱导；片段重复和串联重复共同促进C3H家族扩张；与P. persica共线性最强符合同科关系。未来需验证胁迫诱导的MdC3H基因是关键防御调控因子还是下游标记，以用于抗病苹果精准育种。
研究结论（翻译）：本研究针对苹果C3H基因及其对真菌侵染转录响应表征有限的问题。研究鉴定60个推定C3H基因，并基于系统发育分析分为9组。研究人员开展多种生物信息学分析，以研究基因结构、基因本体、蛋白互作、保守基序和表达谱。观察到片段重复和串联重复，表明其对MdC3H基因家族进化扩张的贡献。表达谱显示，Venturia inaequalis感染后MdC3H06、MdC3H16、MdC3H37和MdC3H42强上调；Alternaria alternata胁迫下MdC3H25上调最高；Gymnosporangium yamadae胁迫下MdC3H16、MdC3H25、MdC3H37和MdC3H42在两个时间点表达增加，提示其参与真菌感染响应。这些发现为未来MdC3H基因功能表征及其在分子育种改良苹果抗病性中的潜在应用提供基础。
