芒果作为广受欢迎的热带水果，其果实的形状千变万化，从圆形、椭圆形、S形到肾形，不一而足。这份多样化的“外貌”，不仅是吸引消费者的颜值担当，更是影响其市场定价、包装效率乃至运输损耗的关键因素。然而，尽管果形如此重要，其背后的遗传与分子调控机制却像一个尚未打开的“黑匣子”。相比之下，西红柿、黄瓜等作物果形研究的“基因地图”已相对清晰，但芒果领域的研究相对匮乏。消费者偏爱外形美观、形状对称的水果，市场也对果实的尺寸和均匀度有严格标准。因此，揭开芒果果形形成之谜，培育出满足市场需求的优良品种，对于整个芒果产业的可持续发展和竞争力提升至关重要。
本研究正是为了破解这一难题而展开。研究者们集成了基因组学、转录组学、表型组学等多种技术手段，对芒果果形进行了系统性的“基因解码”，相关成果发表在园艺学领域期刊《Scientia Horticulturae》上。
为系统探究芒果果形的遗传基础，研究人员主要采用了以下几个关键技术方法：首先，他们利用PacBio长读长与Illumina短读长测序相结合的策略，组装了圆形果品种‘金苹芒’ (‘JPM’) 的高质量染色体级别基因组。其次，针对长果形品种‘红象牙’ (‘HXY’) 与圆形果‘JPM’，在授粉后多个时间点进行采样，通过高通量RNA测序（RNA-seq）获取了全转录组表达数据，并利用加权基因共表达网络分析（WGCNA）筛选与果形相关的基因模块。第三，利用三维（3D）激光扫描技术对75份来自云南丽江芒果种质资源圃、具有不同果形的代表性种质资源进行了9个果形相关性状（如果重、纵径、横径、侧径、体积、果形指数等）的精确测量。最后，对这75份种质进行了全基因组重测序，获得了海量的单核苷酸多态性（SNP）数据，并以此为遗传标记，利用贝叶斯信息和连锁不平衡迭代嵌套关键途径（BLINK）模型，对关键果形性状进行了全基因组关联研究（GWAS），以定位相关的遗传位点。
研究结果
3.1. ‘JPM’品种的基因组组装与注释
研究团队成功组装了圆形果芒果品种‘JPM’的高质量基因组。最终的染色体级别基因组大小为350.66 Mb，包含34,035个蛋白质编码基因，重复序列占比39.28%。BUSCO评估显示基因组完整度达97.70%。通过基因组比较与共线性分析发现，‘JPM’与另一圆形果品种‘Alphonso’基因组间存在强共线性，而与长果形品种‘HXY’的共线性较弱，存在较多倒位事件。这张图展示了‘JPM’基因组的染色体特征及其与‘Alphonso’、‘HXY’基因组的共线性比较，直观呈现了不同果形品种间的基因组差异。
3.2. 基于RNA-Seq表达模式的两种芒果品系比较分析
通过对‘HXY’和‘JPM’果实发育过程（从授粉后1天到110天，共14个时间点）的形态观察和转录组测序，研究人员发现授粉后第3天（DAP3）、第5天（DAP5）和第10天（DAP10）是果形指数（FSI，纵径/横径）趋于稳定的关键发育阶段。在这三个阶段，‘HXY’与‘JPM’之间存在大量差异表达基因。结合WGCNA分析，鉴定出与果形指数最相关的蓝色和黄色基因模块。最终，通过整合差异表达分析与WGCNA，共筛选出467个最可能影响芒果果形发育的关键基因。GO富集分析显示，这些基因主要富集在激素响应调节、细胞骨架依赖性胞质分裂等通路。该图展示了两种果形芒果在不同发育阶段的果实形态、尺寸变化趋势、差异基因韦恩图以及WGCNA模块与性状的关联热图，清晰地揭示了果形发育的关键时期和相关基因网络。
3.3. 种质资源重测序揭示了海量SNP
对75份不同果形的芒果种质进行重测序，获得了超过2600万个SNP。基于高质量SNP构建的系统发育树和主成分分析（PCA）将75份种质划分为4个主要类群，且每个类群在果形上表现出明显的偏好性，例如中国本土选育品种多呈象牙形，而引进品种更偏向圆形，这表明芒果种群进化与果形存在潜在关联。这张图呈现了SNP在染色体上的分布、类型构成以及基于SNP构建的种质系统发育树，并标注了各分支上的代表性果形，直观展示了遗传多样性与果形表型的对应关系。
3.4. 芒果果形多样性分析
利用3D扫描技术获取的9个果形相关性状数据呈正态分布，且性状间存在中等到强的相关性。基于这些性状的PCA分析同样将75份种质分为4个类群，分别代表了大而圆、小而圆以及两种不同类型的长形果，证实了研究材料在果形上具有广泛的代表性。该图通过密度分布图和PCA散点图，展示了所有种质在果形性状上的变异范围和自然聚类情况。
3.5. 6个果形相关性状的GWAS分析
对果重（FW）、果实体积（V）、纵径（LD）、横径（CD）、侧径（SD）和果形指数（FSI）这6个性状进行GWAS分析。以P值小于10-6为阈值，共鉴定出99个与性状显著关联的SNP位点，其中与FSI相关的有75个，与V相关的有24个。在关联区域（SNP上下游各100 kb）内共注释到59个蛋白编码基因。结合这些基因在‘HXY’和‘JPM’果实发育过程中的表达模式，最终筛选出44个关键候选基因。其中，包括LOC123222869(TCP4-like)、LOC123213947(bZIP61-like)、LOC123210819(ABC-G10)、LOC123200572(ARF3) 和LOC123228419(TTL1)在内的8个基因被提议为调控果形变异的关键候选者。该图包含GWAS的曼哈顿图、QQ图以及关键候选基因的表达热图，是定位果形遗传位点和筛选核心基因的直接证据。
结论与讨论
本研究通过整合多组学分析，系统解析了芒果果实形状的遗传基础。主要结论如下：第一，成功组装了圆形果芒果代表品种‘JPM’的高质量基因组，填补了该果形类型基因组资源的空白。第二，明确了DAP3、DAP5和DAP10是芒果果形发育的关键时期，并通过转录组分析筛选出467个果形相关基因。第三，基于75份种质的重测序和3D表型数据，通过GWAS定位到99个与果形性状显著关联的遗传位点。第四，通过整合GWAS与转录组数据，最终鉴定出44个关键候选基因，并重点提出了8个（如TCP4-like、bZIP61-like等）作为调控果形变异的核心候选基因。
讨论部分深入阐述了这些发现的意义。研究指出，芒果果形的多样性与其在印度和东南亚的两个独立驯化中心有关，中国早期引进的品种多源自东南亚，其象牙形果状与此吻合。本研究所用的75份中国栽培品种在遗传和表型上均呈现出丰富的多样性，并可划分为四个明显类群，这为针对不同市场需求进行定向育种提供了材料基础。对于筛选出的关键候选基因，研究者结合其在其他物种中的同源基因功能进行了深入探讨。例如，TCP4-like基因在拟南芥和番茄中已被证实是调控器官大小和形状的关键因子，其在‘JPM’果实发育过程中持续高表达，可能通过抑制细胞增殖和伸长来促进圆形果形成。bZIP61-like作为一个转录因子，可能通过整合生长素和乙烯等激素信号，调控下游细胞壁重塑相关基因的表达来影响果形。ABC-G10属于细胞分裂素转运蛋白家族，可能通过调控激素的长距离运输影响发育。TTL1则可能通过调节微管动力学来协调细胞扩张与增殖，从而影响最终果实大小和形状。
这项研究的重大意义在于，它首次通过如此系统的多组学整合策略，大规模揭示了芒果果形性状的遗传架构和关键调控基因。这不仅深化了我们对芒果这一重要热带果树发育生物学的理解，更重要的是，为未来开展基于基因组的精准育种（genome-based breeding）提供了直接的靶点和丰富的遗传资源。育种者可以利用这些鉴定的遗传位点和候选基因，通过分子标记辅助选择（marker-assisted selection）或基因编辑等技术，更高效地培育出果形优美、均匀一致、更符合市场与加工需求的芒果新品种，从而提升产业效益和竞争力。尽管本研究采用的72份种质群体规模在一定程度上限制了对微效基因的检测，但其所鉴定出的主效基因为后续功能验证和育种应用奠定了坚实的基础。未来，扩大种质群体、结合多环境数据和高密度标记将进一步细化果形的遗传结构，并验证这些关键基因位点的效应。
