咖啡种植海拔是影响风味品质的关键生态因子。本研究采集云南保山同一农场内1000m、1200m、1400m、1600m四个海拔梯度的小粒咖啡（Coffea arabica cv. Cartim）红果，整合16S rRNA和ITS扩增子测序与UPLC-MS/MS非靶向代谢组学，结合HPLC定量分析，系统探究海拔对咖啡樱桃微生物群落和内源代谢物的影响。结果发现：细菌门水平以假单胞菌门（71%-94%）为主，优势属鞘氨醇单胞菌和Pleomorphomonas在1400m处相对丰度最高；真菌以子囊菌门（78%-89%）为主，枝孢霉属和Strelitziana为优势属。海拔升高显著降低咖啡因、葫芦巴碱和绿原酸含量，而杯测评分在1600m最高（80.56分）。差异代谢物数量随海拔差增大而增加（AT4 vs AT3达116种），且与特定微生物存在强相关性。本研究为利用海拔调控微生物提升咖啡品质提供了理论依据。
今天给大家解读一篇5月发表在《Frontiers in Nutrition》上的题目为“Cherry microbiota and metabolites with planting altitude of Coffea arabica in Baoshan of China.”的文章。该研究旨在探究不同种植海拔（1000米至1600米）如何影响阿拉比卡咖啡果实的微生物群落（细菌和真菌）及其代谢物。研究采集了四个海拔梯度的咖啡果实，分别利用高通量测序分析微生物多样性，利用UPLC-MS/MS分析非挥发性代谢物，并通过HPLC法分析了咖啡因、葫芦巴碱和绿原酸三种关键风味前体物质。结果显示，海拔显著影响了微生物群的α多样性和β多样性，以及微生物在门和属水平上的相对丰度。同时，不同海拔间的差异代谢物数量从22到116个不等。咖啡因、葫芦巴碱和绿原酸含量随海拔升高而显著降低。感官分析表明，海拔1600米的咖啡获得最高杯测分数（80.56），达到精品级，且具有更强的果香和烘焙谷物坚果香气。（请持续关注我们，每天为您解读最新见刊的文献!）想薅生信资料羊毛？直接在对话框回复 “资料”，免费领取干货大礼包！包括数据集、绘图代码、图表复现、思路总结、参考文献……0代码！鼠标点点点即可轻松完成5-10分生信SCI全文复现！
不想做实验，没数据，还想要快速发表文章，没问题的！公共数据库就是我们的数据宝藏！没思路不用担心，作为专业的生信团队，我们很乐意为你们效劳，提供研究路线设计和数据挖掘分析，扫码联系我们吧！
团队成员合影（位于上海陆家嘴中心，可随时预约参观）
题目：《中国保山地区种植海拔对阿拉比卡咖啡（Coffea arabica）果实微生物群及代谢物的影响》Cherry microbiota and metabolites with planting altitude of Coffea arabica in Baoshan of China
发表期刊：Frontiers in Nutrition
影响因子：5.1
研究背景：
- 海拔的重要性
海拔是影响农作物品质的关键生态因子，它通过整合温度、湿度、光照和土壤性质等变量，显著影响农产品（如咖啡）的化学成分、营养元素和药理学效应。
- 咖啡风味的关键影响因素
咖啡的风味由多种因素决定，包括咖啡品种、地理起源、气候条件、种植方式、加工方法等。其中，种植海拔对咖啡风味特征的贡献度高达74.2%。
- 微生物的作用
微生物在咖啡产业中扮演核心角色。咖啡果实、土壤、空气和水中的微生物通过产生酶来降解果胶并产生代谢物，从而改变化学组成，影响最终风味。这些微生物群落受到海拔等环境因素的显著影响，形成了独特的“咖啡微生物风土”。
CNSknowall 平台 Pubmed+AI 快速提炼全文要点
研究思路：
- 样品采集
在同一农场的四个海拔点（1000, 1200, 1400, 1600米）手工采摘成熟度一致的完全成熟红色咖啡果实。
- 微生物分析
利用Illumina高通量测序技术，对咖啡果实的细菌（16S rRNA）和真菌（ITS）进行测序，分析其α多样性（Chao, Shannon等指数）和β多样性（PCA, nMDS），并鉴定不同海拔下的优势菌属。
- 代谢物分析
用UPLC-MS/MS技术对咖啡鲜果中的非挥发性代谢物进行全谱检测和鉴定，并通过OPLS-DA等统计方法筛选不同海拔间的差异代谢物（DCMs）。
- 关键前体物质分析
使用HPLC法定量分析生豆中咖啡因、葫芦巴碱和绿原酸的含量。
- 感官评价
将生豆进行统一烘焙，由3位Q-grader咖啡品鉴师按照SCA标准进行双盲杯测，评估香气、风味、酸度、醇厚度等各项指标。
- 关联分析
将所有数据进行整合，分析优势微生物属与筛选出的28种核心差异代谢物之间的相关性。
研究亮点：
- 微生物与代谢物的关联分析
研究不仅描述了不同海拔下的微生物群落和代谢物差异，还通过相关性分析，揭示了特定优势菌属（如 Sphingomonas, Pleomorphomonas, Cladosporium）与特定差异代谢物（如苯丙氨酸-甜菜红碱、果糖-6-磷酸等）之间的强相关性，为理解“微生物-代谢物-海拔”之间的互作关系提供了线索。
- 关键风味物质的量化
精确测量了咖啡因、葫芦巴碱和绿原酸这三种对咖啡风味（苦味、涩味、香气）起关键作用的前体物质随海拔的变化趋势，发现其含量随海拔升高而降低，与高海拔咖啡更好的风味品质（更高的杯测分数）形成对比。
- 高海拔微生物作为发酵菌种的潜力
研究指出，不同海拔的微生物群落存在显著差异，特别是高海拔咖啡中的独特微生物，这为未来开发具有特定风味的咖啡发酵菌种提供了潜在机会。
研究结果：
- 风味前体物质
1000米海拔的咖啡豆中，咖啡因（12.97%）、葫芦巴碱（8.35%）和绿原酸（37.24%）含量最高。总体趋势是，这三种物质的含量随着海拔升高而显著降低，但在1400米和1600米样品间差异不显著。
- 感官分析
1600米（AT4）的咖啡杯测总得分最高（80.56），被认定为精品级，其香气/风味与其它海拔样品差异显著，展现出强烈的果味、烘焙谷物和坚果香气。
- 微生物多样性
- 细菌
群落以变形菌门（Pseudomonadota）为主。鞘氨醇单胞菌属（Sphingomonas）和嗜临单胞菌属（Pleomorphomonas）是优势菌属，均在1400米（AT3）海拔达到最大相对丰度。
- 真菌
群落以子囊菌门（Ascomycota）为主。枝孢菌属（Cladosporium）和Strelitziana属是优势菌属。此外，Vishniacozyma和Setophoma在1000米（AT1）海拔更重要，而Setophoma和亚隔孢壳属（Didymella）在1400米（AT3）海拔更重要。
- 代谢物变化
在咖啡鲜果中鉴定出1040种代谢物。不同海拔组间存在大量差异代谢物（DCMs）。例如，在最高（1600m）和最低（1000m）海拔比较中，发现了77种DCMs。通过进一步筛选，获得了28种在OPLS-DA模型中VIP值>1且p值≤0.05的核心差异代谢物。
- 相关性分析
-
细菌Sphingomonas与(S)-溴内酯等呈强负相关。
-
细菌Methylobacterium与苯丙氨酸-甜菜红碱呈强正相关，与多种糖苷和糖类（如果糖-6-磷酸）呈强负相关。
-
真菌Vishniacozyma与多种代谢物（如D-2-羟基戊二酸、γ-丁内酯等）呈强正相关。
研究总结：
结论： 种植海拔是影响咖啡果实微生物群落和内源性代谢物特征的关键因素，这些变化共同塑造了咖啡的风味品质。高海拔地区的咖啡具有较低的咖啡因、葫芦巴碱和绿原酸含量，但同时拥有更佳的风味和更高的感官评分。
讨论与启示：
- 海拔的作用机制
海拔对咖啡风味的影响主要通过两条途径实现：一是直接改变咖啡植物的代谢，影响生豆中风味前体物质（如咖啡因、绿原酸）的积累；二是通过塑造独特的微生物生态位（即“风土”），影响果实的微生物附着和发酵过程，从而产生或修饰风味化合物。
- 微生物与风味的潜在联系
研究中发现的不同海拔下优势菌属（如Sphingomonas 具备促生和增强抗逆性能力，Cladosporium 常见于植物体），以及它们与代谢物的显著相关性，暗示了微生物在咖啡风味形成中扮演着重要角色，不仅是发酵阶段，也贯穿于果实生长过程中。
- 实践应用与局限性
- 应用
高海拔咖啡中特有的微生物群系有潜力作为特色发酵的候选菌种，用于开发具有独特风味的高品质咖啡。通过调控微生物群落，或可实现对特定风味的定向改善。
- 局限性
本研究仅采集了保山地区的样品，未覆盖云南其他主要产区（如普洱、临沧），结论的区域代表性有限。同时，与海拔相关的具体环境变量（如温度、湿度、土壤性质）数据未能同步获取，因此无法完全归因于海拔单一变量。文章最后提出了可验证假说，例如通过微生物互换实验和共培养实验，来进一步证明微生物与代谢物之间的因果关系。
结果译文：
校准曲线、线性范围、检测限和定量限、验证参数以及化合物鉴定结果如表1所示。基于校准曲线，计算了咖啡因、葫芦巴碱和绿原酸的含量（图1A）。低海拔咖啡中的咖啡因、葫芦巴碱和绿原酸含量与高海拔咖啡有显著差异，而在1400米和1600米之间则无显著差异。1000米处的咖啡豆具有最高的咖啡因、葫芦巴碱和绿原酸含量，分别为（12.97±0.51）%、（8.35±0.50）%和（37.24±0.28）%。
详细属性的评分如图1B所示。均一性、甜度和干净度的评分为10分，表明这些咖啡样品无缺陷。杯测质量结果显示，不同海拔咖啡样品中AT4的杯测得分最高，为80.56±2.89分，这意味着AT4的咖啡被视为精品级。AT4的香气/香味和风味与其他海拔咖啡有显著差异。此外，AT4具有强烈的水果特征，明显的烤谷物和坚果香气。
结果显示，细菌共优化获得802,385条序列，真菌962,447条，细菌总碱基数为302,431,929，真菌总碱基数为294,154,598。平均序列长度分别为377 bp和306 bp。使用α多样性（图2）评估微生物生物多样性。对于细菌，AT2的Chao和ACE指数最高，达到725.77±255.10和729.58±259.63，而最低值出现在AT3，分别为316.67±78.55和316.75±78.43。Simpson指数范围为0.047±0.035（AT4）至0.21±0.0086（AT3），而Shannon指数范围为2.71±0.30（AT3）至4.67±0.60（AT4）。对于真菌，AT4表现出最高的α多样性（Chao、Shannon和ACE指数），分别达到519.81±93.22、4.44±0.31和519.62±82.97。Simpson指数范围为0.047±0.0020（AT2）至0.071±0.017（AT3）。
使用Mothur（版本1.30.2）进行β多样性分析，获得主成分分析（PCA）、非度量多维尺度分析（nMDS）和置换多元方差分析（PERMANOVA）（补充图S1）。在细菌中，PC1（37.25%）和PC2（19.67%）的PCA显示了四种不同海拔咖啡样品的区分，解释了56.92%的方差。在真菌中，PC1（29.84%）和PC2（16.17%）解释了46.01%的方差。nMDS分析显示了四个不同海拔样品的分布，突出了微生物群落谱的差异。在真菌群落中，四个种植海拔之间观察到了最显著的不相似性。在细菌群落中，AT3和AT4表现出最大的不相似性。此外，PERMANOVA显示细菌的R²为0.47（p=0.009），真菌的R²为0.49（p=0.001）。
从四个咖啡种植海拔中鉴定出六个细菌门，包括Pseudomonas、Actinomycetota、Bacteroidota、Bacillota、Chloroflexota等（图3A）。其中，Pseudomonas是优势细菌门，其相对丰度范围从71.17%（AT2）到94.21%（AT3），在1400米处达到最大相对丰度。在Pseudomonas之后，Actinomycetota也是一个重要的细菌门，范围从3.97%（AT3）到23.02%（AT4）。在AT2中，Bacteroidota占6.76%。鉴定出21个属，包括Sphingomonas、Pleomorphomonas、Pseudomonas、Tardiphaga、Methylobacterium、Bradyrhizobium、Actinomycetospora、Rhizobium、Comamonas、Microbacterium、Bosea、Chryseobacterium、Amnibacterium、Cupriavidus、Novosphingobium、Nakamurella、Microlunatus、Aureimonas等（图3B）。根据属水平的相对丰度，Sphingomonas（12.52%~30.74%）和Pleomorphomonas（7.00%~31.82%）是不同种植海拔咖啡中的优势细菌属。同时，Sphingomonas和Pleomorphomonas在AT3中显示出最大相对丰度。
对于真菌，在门水平上鉴定出六个门（Ascomycota、Basidiomycota、Chytridiomycota、Mucoromycota、未分类真菌等）（图3C）。Ascomycota是优势真菌门，范围从77.84%到88.52%。最大相对丰度为AT2中的88.52%。在Ascomycota之后，Basidiomycota也是一个重要的真菌门，范围从10.30%（AT2）到21.54%（AT1）。鉴定出21个属，包括Cladosporium、Strelitziana、Aureobasidium、Didymella、Vishniacozyma、Hannaella、Fusarium、Phaeosphaeria、Colletotrichum、Setophoma、Cyphellophora、Papiliotrema、Alternaria、Paracamarosporium、Acrocalymma等（图3D）。在属水平上，Cladosporium和Strelitziana真菌显示出最高的相对丰度，范围分别为19.11%（AT1）至27.49%（AT2）和9.57%（AT2）至16.94%（AT4）。此外，Vishniacozyma（8.69%）和Setophoma（8.27%）在AT1中也很重要。Setophoma和Didymella在AT2中分别占9.57%和9.24%，在AT3中分别占11.46%和4.67%。然而，Colletotrichum在AT4中占5.04%。
在不同种植海拔的咖啡中，共鉴定出16个超类别的1040种代谢物（补充图S2）。此外，不同种植海拔咖啡的差异变化代谢物（DCMs）以VIP>1.0、p≤0.05且FC<0.67或FC>1.5为标准，如图4所示。在AT2与AT1比较中检测到的22种DCMs中，17种显示上调（例如环黄芪醇、鳄梨烯、白三烯D5、tsugarioside B、依替米星、泛酸、四环素、三咖啡酰亚精胺、zerano、去乙酰长春碱酰胺、二异丙醇胺、棕榈酰阿糖胞苷、甲硫氨酰-酪氨酸等），5种显示下调（例如甜菊苷、7α-羟基睾酮、异海松酸、镰刀菌素C、毒胡萝卜素；图4A）。
同样，在AT3与AT1比较中鉴定出的68种DCMs中，33种显示上调（例如环黄芪醇、开环异落叶松脂素、D-山梨醇、Lac Dye、sambacin、蔷薇酸、维甲酸、毒胡萝卜素、vignatic acid A、螺环内酯C、乔松酮5-葡萄糖苷和licoflavanone、开环异落叶松脂素等），35种显示下调（例如头孢甘氨酸、岩藻黄质、Ac-(5H-二苄基(A,D)环庚烯-10,11-二氢甘氨酸-Leu-Asp-Ile-Ile-Trp)、槲皮素3-槐糖苷-7-鼠李糖苷、多柔比星、fevicordin B 2-[鼠李糖基-(1->4)-葡萄糖基-(1->6)-葡萄糖苷]、2,3-丁二醇芹菜糖基葡萄糖苷、硝基苄基-6-硫代肌苷等；图4B）。
在AT3与AT2比较中检测到的38种DCMs中，20种显示上调（例如甜菊苷、(3B,20R,22R)-3,20,27-三羟基-1-氧代睡茄-5,24-二烯醇ide 3-葡萄糖苷、(5Alpha,6Beta,14Alpha,20R,22R)-5,6,14,20,27-五羟基-1-氧代睡茄-24-烯醇ide、藤黄酸等），18种显示下调（例如gambafoultalin、(S)-琥珀酰二氢哌啶酰胺、alliofuroside A、谷甾醇β-D-葡萄糖苷、Ile-Ile-Ile-Pro、精氨酰-精氨酰-芳酰-酪氨酰-异亮氨酰-亮氨酸、dyrrhorin A 1-8、微囊藻毒素-Lr、三七皂苷J、三咖啡酰亚精胺、Ac-(5H-二苄基(A,D)环庚烯-10,11-二氢甘氨酸-Leu-Asp-Ile-Ile-Trp)等；图4C）。
在AT4与AT1比较中鉴定出的77种DCMs中，44种上调（例如5-丁基-3,4-二甲基-2-呋喃十一烷酸、葡萄糖基半乳糖基羟赖氨酸、三咖啡酰亚精胺、埃博霉素A、山茶苷B、5-脱氧kievitone水合物、红霉素、davercin、des-Arg(9)-缓激肽、ebratide、亮氨酸-甜菜黄素、酪氨酰-苯丙氨酸等），33种下调（例如豆甾醇、28-吡喃葡萄糖基-3-甲基齐墩果酸、5-过氧羟基二十碳四烯酸、甲硫氨酰-酪氨酸、异青霉素N、lyciumin B、环丙氨嗪、3-羧基-1-羟丙基硫胺素二磷酸、4-吡哆酸、5-庚基四氢-2-氧代-3-呋喃甲酸、三乙二醇双(3-叔丁基-4-羟基-5-甲基苯基)丙酸酯等；图4D）。
在AT4与AT2比较中检测到的57种DCMs中，19种上调（例如四降丙烷二酸、5-过氧羟基二十碳四烯酸、3-羧基-1-羟丙基硫胺素二磷酸、琥珀酸、甲硫氨酰-酪氨酸、D-2-羟基戊二酸、lyciumin B、D-赤藓糖、替马西泮葡萄糖醛酸苷、(4-(2-(3-(环戊氧基)-4-甲氧基苯基)-2-苯基乙基)吡啶等），38种下调（例如10-硝基油酸、氯贝特葡萄糖醛酸苷、藤黄酸等；图4E）。
在AT4与AT3比较中检测到的116种DCMs中，72种显示上调（例如萼花醇A、楤木皂苷A、奇任醇、cyclosquamosin G、cyclosquamosin C、酪氨酰-苯丙氨酸、氯贝特葡萄糖醛酸苷、隐绿原酸、松三糖、蔗糖、α-茄碱、4-氧代-依诺沙星、去乙酰长春碱酰胺、3-(3,5-二氯苯基)-4-羟基-1,3-噻唑-2-酮、硝基苄基-6-硫代肌苷等），44种显示下调（例如(9Z)-(13S)-12,13-环氧十八碳-9,11-二烯酸酯、27-羟基异芒果酸、谷氨酰亮氨酰精氨酸、琥珀酸、lyciumin B、D-2-羟基戊二酸、果糖-6-磷酸、4-[(2-呋喃基甲基)硫基]-2-戊酮、4-吡哆酸、tricin 7-[鼠李糖基-(1->2)-半乳糖醛酸苷]、鸟苷5'-单磷酸等；图4F）。
最后，使用OPLS-DA模型（R²=0.904，Q²=0.826）中的VIP和单因素ANOVA的p值来鉴定不同海拔样品之间的DCMs，阈值VIP>1且p≤0.05，共发现28种DCMs（表2）。此外，分析了这些DCMs与细菌、真菌中相对丰度前五的微生物之间的相关性（图5）。根据相关系数（r）的值，0.8~1.0表示极强相关，0.6~0.8表示强相关（27）。Sphingomonas与(S)-溴烯醇内酯、M4和山茶苷A呈极强负相关。Pleomorphomonas与苯丙氨酸-甜菜黄素呈极强负相关，与棘囊酸-3-O-葡萄糖苷、果糖-6-磷酸和苏马树脂酸呈极强正相关。Pseudomonas与(S)-溴烯醇内酯和fistulosin呈极强负相关。Tardiphaga与苯丙氨酸-甜菜黄素呈极强负相关，与(S)-溴烯醇内酯1、果糖-6-磷酸和苏马树脂酸呈极强正相关。Methylobacterium与27-羟基异芒果酸、棘囊酸-3-O-葡萄糖苷、果糖-6-磷酸、苏马树脂酸呈极强负相关，与苯丙氨酸-甜菜黄素呈极强正相关。Bradyrhizobium与呋喃酮呈极强正相关。Cladosporium与α-茄碱和果糖-6-磷酸呈极强负相关，与glucopiericidin A呈极强正相关。Strelitziana与D-2-羟基戊二酸、苯丙氨酸-甜菜黄素和4-氧代-依诺沙星呈极强负相关，与去乙酰长春碱酰胺和环黄芪醇呈极强正相关。Aureobasidium与α-茄碱呈极强负相关，与glucopiericidin A呈极强正相关。Didymella与glucopiericidin A和赤藓糖醇呈极强正相关，与茶黄素和山茶苷A呈极强负相关。Vishniacozyma与赤藓糖醇呈极强正相关，与姜酮醇和fistulosin呈极强负相关。
更多结果和补充图表：doi： 10.3389/fnut.2026.1817512
长按二维码关注我们，用最短的时间和最高的效率学习更多生信思路!
扫描上方二维码或登录平台官网后添加CNSknowall客服微信咨询！官网地址：https://cnsknowall.com
CNSknowall：24年最新问世的遥遥领先的颠覆性科研数据（0代码生信+统计学）分析平台，同时含有机制图模块（原创3000多素材和机制图模板）+AI一键生成高质量比国自然标书初稿+汉化版Pubmed融合Deepseek高效筛选目标文献同时一键提炼全文核心创新点+SCI文献例句/语料检索模块+全文翻译+文献求助+图片查重+期刊查询+OPenAI官方GPT接口，＞500款CNS级别图表皆可一秒内一键出图，登录即秒变数据分析大神，体验前所未有的便捷数据分析之旅，开启科研天骄之路！
可向下滑动发掘更多科研秘籍！
