《Theoretical and Applied Genetics》：The stability of fatty acid composition in sunflower oil is dependent on environment and affected by structural variation
编辑推荐：
在向日葵（Helianthus annuus L.）中，种子中脂肪酸组成（主要为油酸、亚油酸、硬脂酸和棕榈酸）对油品质至关重要。尽管如此，该性状的遗传基础及其与环境的互作仍知之甚少。理解这种互作对于气候变化背景下向日葵的改良至关重要。本研究中，研究人员将来自向
在向日葵（Helianthus annuus L.）中，种子中脂肪酸组成（主要为油酸、亚油酸、硬脂酸和棕榈酸）对油品质至关重要。尽管如此，该性状的遗传基础及其与环境的互作仍知之甚少。理解这种互作对于气候变化背景下向日葵的改良至关重要。本研究中，研究人员将来自向日葵SAM群体（向日葵关联作图群体）的脂肪酸组成测定数据（涵盖8个环境，跨越广泛的地理梯度）整合到全基因组关联分析（GWAS）中。SAM群体由287个品种组成，代表向日葵约90%的多样性，拥有220万个高质量SNP（最小等位基因频率MAF>5%）。为增强统计效力，研究人员采用四种不同输入进行多变量GWAS：（i）每个环境内的平均脂肪酸组成；（ii）剔除高油酸品种后每个环境内的平均脂肪酸组成；（iii）通过重复样本间标准误量化的环境内性状稳定性（α稳定性）；（iv）Eberhart和Russell的β，该指标量化跨环境的性状稳定性（β稳定性）。所有四种分析均产生高度显著关联的SNP。研究人员发现，高油酸品种表现出高β性状稳定性，导致分析（i）和（iv）之间的标记有大量重叠，且分析（i）中信号在不同环境间相当一致。对于分析（ii）和（iii），显著标记在不同试验间趋于变化。在所有分析的显著SNP中，鉴定出147个候选基因，包括有潜力的候选基因如15个脂肪酸代谢基因、6个热休克蛋白和22个转录因子。最后，研究人员发现5号染色体上由两个侧翼倒位序列构成的大片段渐渗，与佐治亚州试验中的α稳定性重合，提示其在高温条件下脂肪酸组成稳定性中发挥作用。
**论文解读：向日葵油脂肪酸组成稳定性的遗传基础与环境互作**
**研究背景与问题**
向日葵（Helianthus annuus L.）是全球重要的油料作物，其种子油中脂肪酸组成（主要为油酸Oleic acid、亚油酸Linoleic acid、硬脂酸Stearic acid和棕榈酸Palmitic acid）直接影响油的工业用途和营养价值。然而，这一复杂性状的遗传基础及其与环境（如温度、水分等）的互作机制尚不清楚。现有研究多集中于单一脂肪酸的基因组关联分析（GWAS）或数量性状位点（QTL）定位，但缺乏多变量同时分析多个脂肪酸的GWAS，也未曾系统研究脂肪酸表型稳定性在气候梯度下的遗传控制。此外，某些高油酸（HO）品种的HO表型在世代间不稳定，成因未知，成为育种中的难题。因此，在气候变化背景下，解析脂肪酸组成及其稳定性的遗传机制，对向日葵定向改良至关重要。
**研究内容与结论**
研究人员利用向日葵关联作图群体（SAM群体，287个品种，覆盖约90%向日葵遗传多样性），在北美8个不同环境（包括加拿大不列颠哥伦比亚省、美国明尼苏达州、爱荷华州和佐治亚州）田间种植，测定种子中四种主要脂肪酸含量，并基于220万个高质量SNP（MAF>5%）开展多变量GWAS。分析包括四种输入：环境内平均脂肪酸组成（含或不含高油酸品种）、环境内α稳定性（重复标准误）和跨环境β稳定性（Eberhart和Russell的β）。主要结论：高油酸品种表现出较高的β稳定性，且与平均脂肪酸组成的GWAS信号高度重叠；剔除高油酸品种后，多环境支持的QTL减少；在α稳定性分析中，不同环境支持的位点差异较大，但5号染色体上一个已知的大片段渐渗（含两个倒位序列）与佐治亚州高温试验中的α稳定性显著关联。共鉴定147个候选基因，包括15个脂肪酸代谢基因、6个热休克蛋白和22个转录因子。该研究为通过分子育种维持向日葵脂肪酸组成稳定性提供了重要靶点。论文发表在《Theoretical and Applied Genetics》。
**关键技术方法**
1. **群体与样本来源**：SAM群体（287个品种，包括近交系、开放授粉品种的自交后代、地方品种等），来源于北美8个环境（不列颠哥伦比亚省2010年；明尼苏达州2015年、2016年早播和晚播；爱荷华州2010年、2013年、2014年；佐治亚州2010年）。
2. **基因分型**：从HA 412HOv2.0参考基因组调用SNP，经GATK质控、过滤（minQ=100，MAF≥0.05，最大缺失率0.9），BEAGLE v5.3填补缺失，最终获得2.2M高质量SNP。
3. **多变量GWAS**：采用GEMMA v0.98.5进行多变量分析，分别以四种输入（平均脂肪酸组成、剔除高油酸品种的平均脂肪酸组成、α稳定性、β稳定性）进行关联分析，同时使用vcf2gwas v0.8.7进行单变量GWAS。
4. **稳定性量化**：α稳定性定义为同一环境内重复样本间脂肪酸组成的标准误；β稳定性采用Eberhart和Russell（1966）的方法，即各环境均值排序后的回归斜率。
5. **候选基因鉴定**：基于显著SNP上下游1 Mb区域，筛选在发育种子中高表达（log2倍变化>2，与10种其他组织比较）且具有MapMan功能注释的基因。
**研究结果**
* **表型分析**：通过相关性分析发现，油酸与亚油酸呈显著负相关（因FAD2-1基因不同等位基因导致），且油酸和亚油酸组成与纬度、水汽压亏缺、日最低温（夜间温度）和日最高温等气候变量显著相关，尤其在播种后第三个月（种子发育期）。硬脂酸与其他脂肪酸的关联模式不同，其与水分相关气候变量无显著相关。GGE双标图显示，当包含高油酸品种时，油酸与棕榈酸之间存在权衡；剔除高油酸品种后，则呈现饱和脂肪酸（棕榈酸和硬脂酸）与不饱和脂肪酸（油酸和亚油酸）之间的直接权衡。
* **GWAS on mean fatty acid composition**：在所有环境中均检测到高度显著关联，最高峰位于14号染色体（FAD2-1突变位点所在单倍块），此外1、11、13、16号染色体上的峰获得多个环境支持。多变量GWAS结果与β稳定性GWAS结果在峰位置上有大量重叠。
* **GWAS on trait stability across environments (β stability)**：β稳定性的GWAS在1、11、13、14、16号染色体上检测到显著关联，与平均脂肪酸组成的QTL重合，且13号染色体上的稳定性关联SNP分布区域更广。
* **Omission of high oleic varieties in GWAS**：剔除含FAD2-1重复的28个高油酸品种后，显著SNP的p值大幅降低，且多环境支持的QTL减少。β稳定性GWAS未达到严格Bonferroni阈值，但4个位点满足宽松阈值。
* **GWAS on trait α stability within environments**：在6个有重复试验的环境中进行α稳定性GWAS，除4号染色体外，所有染色体均出现显著峰，但仅5、7、13号染色体上的峰获得多个环境支持。5号染色体上145 Mb和185 Mb处的两个显著峰位于已知的148–177 Mb倒位/渐渗区两侧，提示该渐渗与种子内脂肪酸组成均匀性相关。
* **Candidate gene analysis**：从四种GWAS分析中分别鉴定出候选基因：平均FA分析（含FAD突变）64个、平均FA分析（剔除FAD突变）68个、β稳定性分析22个、α稳定性分析58个。多个基因在不同分析中重叠，如脂肪酸代谢基因（油质蛋白oleosin、ω-3/ω-6脂肪酸去饱和酶）、热休克蛋白、转录因子，甚至一个siRNA生物合成调控复合体基因。α稳定性与β稳定性的候选基因重叠极少，提示两者遗传基础不同。
**讨论与结论**
讨论部分指出，高油酸品种的β稳定性较高，且与平均脂肪酸组成的GWAS信号重叠，说明高油酸育种可能同时选择了稳定性，或稳定性在育种过程中被间接选择。FAD2-1位点对脂肪酸组成起主导作用，但其他位点也贡献显著效应。α稳定性与β稳定性的遗传基础差异明显，前者可能更多受高油酸品系向野生型脂肪酸组成回复的影响。5号染色体上的大片段渐渗（已知影响舌状花宽度、种子大小和果皮强度）在佐治亚州高温试验中与α稳定性关联，提示该结构变异在高温条件下维持脂肪酸组成均匀性中的作用。研究结论：脂肪酸组成和稳定性受FAD2-1位点控制，并受基因组多个位点影响；跨环境β稳定性遗传较为简单，与平均脂肪酸组成位点重叠；环境内α稳定性则在不同环境间变异较大，但少数位点（如5号染色体结构变异）重要。该工作为向日葵脂肪酸谱稳定性的育种提供了靶点，对全球重要油料作物——向日葵的改良具有关键意义。
