复垦矿地芒属（Miscanthus）种间F1群体农艺与生物量相关性状的全基因组标记–性状分析：LG9重复提示信号与LG3年份依赖关联

《Plants》：Genome-Wide Marker–Trait Analysis of Agronomic and Biomass-Related Traits in an Interspecific Miscanthus F1 Population Grown on Reclaimed Mine Land

　　

芒属（Miscanthus）是适合边际和退化土地的多年生生物能源作物，然而对种间杂种跨年度生物量相关变异的遗传基础仍知之甚少。研究人员在复垦矿地地点对来自Miscanthus sacchariflorus × M. lutarioriparius的100个F₁基因型进行了两个生长季评价，并使用已有包含19个连锁群（linkage groups, LGs）上461个SNP（single nucleotide polymorphism，单核苷酸多态性）标记位置的连锁图谱分析了16个农艺、生物量相关和生理性状。没有单一性状峰值超过其性状特异性5%置换阈值。最清晰的重复提示性信号出现在LG9近端区域：marker2342（11.571 cM）在株高和分蘖长度上显示最强跨年度信号，而相邻marker2343（10.561 cM）在分蘖数上给出最强跨年度信号。四个相关构型性状也通过主成分分析（principal component analysis, PCA）概括。PC1解释其变异的73.2–76.6%，2019年LG3上的marker0580（12.516 cM）超过全基因组阈值（LOD 3.67 vs. 3.47）。该关联在去除单个稀有基因型类别观测后仍显著，但2018年未重现；正式标记×年份互作检验证实四个组成性状均存在时间异质性。生物量和生理信号较弱且重复性较低。因此，该研究识别了一个重复但中等的LG9模式，以及在该复垦地点表达的更强、时间可变的LG3信号；它并未建立特异性赋予复垦矿地适应性的位点。这些区域为在评估其标记辅助育种价值之前，使用更大群体、更密集标记覆盖、更多年份和对比环境进行验证提供了优先事项。

芒属（Miscanthus）是多年生C4草，因高生物量生产和相对低投入而受重视，是边际土地候选生物量作物。复垦矿地是这一特性相关的场景之一，常规作物生产常受限制，而田间研究表明芒属可在该条件下建植并生产生物量。东亚是芒属多样性中心。M. sacchariflorus贡献广泛适应性，M. lutarioriparius以旺盛生长和高生物量潜力著称。M. lutarioriparius已有基因组和群体转录组资源，使这两个物种间杂交可用于研究生长和生物量相关性状的分离变异。生物量产量由多个相关性状决定，包括株高、分蘖数、茎尺寸、叶片发育、物候和单分蘖生物量。M. sinensis和M. sacchariflorus的定位研究已鉴定许多这些组分的位点，但效应不总在年份或环境间重复。对于多年生田间性状，跨季节重复性可与单年效应大小同样重要。因此，多年度数据对解释生长和生物量的基因组信号很重要。基于SNP的连锁图谱已使高度杂合芒属群体的数量变异研究成为可能。RNA-seq衍生标记已用于构建覆盖19个连锁群的图谱，F₁群体已用于定位株高、开花、叶性状和生物量生产力。对M. sacchariflorus × M. lutarioriparius群体知之甚少，尤其在半干旱复垦矿地条件下。该种间杂交此前已在复垦矿地评价生长、生物量生产和生理表现，并在该田间环境进行了转录组分析。然而，其定位遗传变异尚未针对多年度测定的农艺和生理表型进行系统评价。研究人员在此评价了黄土高原复垦矿地地点两个生长季中100个M. sacchariflorus × M. lutarioriparius F₁基因型，并使用已有461个标记连锁图谱作为标记–性状扫描的位置框架。研究人员询问哪些基因组信号跨年份和相关性状重复，哪些具有时间变异性。分析鉴定出LG9近端一个重复但低于阈值的株型信号，以及LG3上marker0580一个更强的2019年特异多变量株型信号。由于所有基因型仅在单一复垦地点评价，这些关联被解释为在该环境表达的遗传效应，而非特异性赋予复垦矿地适应性的位点。

关键技术方法（不超过250字）：样本队列来源于2010年构建的种间F₁定位群体，母本为M. sacchariflorus（accession 02487，中国甘肃），父本为M. lutarioriparius（accession 04081，中国湖南长沙），杂交真实性此前用诊断SSR（simple sequence repeat，简单序列重复）标记确认。田间试验2018年建于中国山西孝义铝矿复垦地、黄土高原（37.12° N，111.45° E），本分析纳入100个F₁基因型。主要方法包括：使用已有461个SNP标记、19个连锁群的连锁图谱作为位置框架；对16个农艺、生物量相关和生理性状进行两年表型评价；以标记观察基因型类别为分类因子进行全基因组标记–性状扫描，计算LOD（logarithm of odds，优势对数）和PVE（phenotypic variance explained，表型方差解释率），用1000次表型置换确定5%全基因组阈值；用模拟评估检测力；用PCA整合四个株型性状；用Welch t检验、Cohen’s d和Benjamini–Hochberg FDR（false discovery rate，错误发现率）进行基因型类别比较；用重复测量GEE（generalized estimating equations，广义估计方程）检验标记×年份互作。

研究结果：
3.1 表型变异与性状分布：研究人员在种间芒属F₁群体中评价16个农艺、生物量相关和生理性状，年度表型为三株测量的基因型级算术均值。株高、分蘖长、分蘖径、分蘖数、叶长、叶宽、叶面积和生理性状对100个定位F₁基因型有完整年度观测；单分蘖鲜重和全株鲜重在两年中各有96个基因型可用，因4个基因型缺少两项生物量测量，叶角仅2019年测定（N=100）。群体在所有测量性状上表现广泛连续变异。全株鲜重表型范围广，部分后代值超出亲本点估计，但因未正式检验亲本变异和超亲分离，故描述为广泛表型分离而非超亲遗传证据。单分蘖鲜重也呈连续数量变异。株型性状和生理性状也广泛变异。跨年相关显示四个主要株型性状中等表型稳定性：株高r=0.491、分蘖长r=0.527、分蘖径r=0.485、分蘖数r=0.684，均p<0.001；而光合速率r=0.171、WUE_e r=0.014、WUE_i r=0.048，时间可塑性较大。多数性状连续，但部分生物量、形态和生理性状显著偏离正态，研究人员对部分性状采用秩基逆正态变换评估主要关联稳健性。
3.2 用于全基因组标记–性状分析的遗传连锁图谱：所用图谱包含461个定位SNP位置，分布于19个连锁群，与芒属单倍体补体19条染色体（2n=38）一致。本研究未重建图谱，而使用此前构建图谱保留的标记位置。由于亲本高度杂合、后代为全同胞F₁家系，标记分离在CP（cross-pollinator，异交授粉者）框架下解释。保留的100基因型矩阵中，160个标记显示AA/AB类别，124个显示AB/BB类别，173个显示AA/AB/BB三类，4个为AB不变。这些是归档分析数据集中观察到的基因型类别构型，而非重建的正式CP分离码，原始亲本相位分配不可用。图谱总遗传长度3215.56 cM，各连锁群从LG3的95.03 cM到LG9的231.95 cM；每连锁群标记数12至48，平均24.26，相邻标记平均间距4.54至15.61 cM。标记覆盖不均匀，LG5、LG14和LG19最大相邻标记间隙分别为41.05 cM、91.63 cM和38.83 cM，而LG9含48个标记、平均间距4.94 cM。因上游标记过滤和分离记录不可用，大间隙原因无法确定。不均匀密度限制定位分辨率及区分紧密连锁效应的能力。
3.3 全基因组标记–性状关联分析：研究人员对16个性状在19个连锁群保留标记集上进行了全基因组扫描，跨年调整基因型均值（原文称为BLUE值）作为两年性状主要表型，并检查年度特异扫描评估时间一致性。性状特异5%全基因组阈值由1000次表型置换确定。没有单性状关联峰超过相应置换阈值，因此在当前群体大小和标记密度下未检测到全基因组显著单性状标记关联。尽管无显著单性状关联，若干提示信号由生物学相关性状共享。株型最清晰模式在LG9近端：marker2342（11.571 cM）对株高（LOD=2.26；PVE=9.9%）和分蘖长（LOD=2.51；PVE=10.9%）显示最强BLUE关联，相邻marker2343（10.561 cM）为分蘖数最强信号（LOD=2.13；PVE=9.4%）。生物量相关性状中，LG12上marker0907（57.272 cM）为单分蘖鲜重和全株鲜重最强BLUE标记，分别解释约10.8%和8.5%表型变异。生理性状中，LG3上marker0277（5.114 cM）为胞内CO₂浓度和WUE_i最强BLUE关联，PVE分别为10.6%和10.7%。蒸腾速率另在LG9的marker3163有提示关联。因均未过阈值，这些被视为提示性关联而非显著QTL（quantitative trait locus，数量性状位点）。LG9近端因跨相关株型性状重复出现而被进一步检查，但该重复不视为独立重复，因组分性状相关。
3.4 LG9重复提示株型簇内的基因型效应：在LG9近端簇内，marker2342对株高和分蘖长最强，相邻marker2343对分蘖数最强。研究人员按marker2342基因型类别分别检查2018、2019和跨年BLUE表型差异。AA类别个体通常比AB类别更高且分蘖更长。2018年AA平均株高约123.5 cm，AB约112.6 cm，差约10.9 cm；2019年为154.8 vs. 143.0 cm，跨年BLUE差约11.3 cm。分蘖长同向：AA在2018比AB长8.6 cm，2019长11.4 cm，跨年BLUE差约10.0 cm。FDR校正后，株高在2018和BLUE数据集显著，分蘖长在2018、2019和BLUE数据集显著。正式重复测量标记×年份检验对株高（p=0.825）和分蘖长（p=0.417）不显著，对分蘖径不显著（p=0.788），但对分蘖数显著（p=0.00567；FDR校正p=0.0227）。因此marker2342不应被描述为在所有株型性状上一致稳定。相邻marker2343捕获最强BLUE分蘖数信号，但观察类别高度不平衡（AA=1，AB=99），未进一步解释其基因型类别效应。LG9近端被解释为重复但中等的提示关联区域，跨年一致性最强证据涉及株高和分蘖长，而非整个株型。
3.5 多变量分析识别LG3上显著但年份依赖的株型关联：因株高、分蘖长、分蘖径和分蘖数强相关并共同贡献整体株型，研究人员进行PCA整合为营养生长复合指标。四性状标准化后PCA。PC1解释2018、2019和BLUE数据集总变异的73.20%、74.17%和76.57%。四性状PC1载荷均为正，前三者贡献较强，高PC1代表植株更高、分蘖更长更粗且分蘖数更多。以PC1为复合表型的全基因组扫描显示比多数单株型性状更强信号。2018年最强PC1关联在LG8 marker0399（38.237 cM），LOD=3.06，低于5%阈值3.47。2019年PC1扫描在LG3 marker0580（12.516 cM）LOD=3.67，超过阈值3.47，解释PC1约15.5%，为本分析最强全基因组关联。BLUE PC1最高峰也在marker0580（LOD=2.95；约12.7%方差解释），但未超过BLUE特异阈值3.39。marker0580保留基因型类别为BB=80、AB=19、AA=1。排除少于5个体的类别后，marker0580与2019 PC1关联仍显著，LOD=3.24，超过相应5%置换阈值3.16，故2019关联并非仅由单个稀有观测驱动。2019年BB比AB个体株高150.35 vs. 131.75 cm（差18.60 cm，p=0.00023）；分蘖长89.83 vs. 73.28 cm（差16.55 cm，p=0.000012）；分蘖径6.90 vs. 5.83 mm（p=0.000039）；分蘖数4.27 vs. 2.84（p=0.0041）。四组分性状一致。marker0580信号在2018年显著较弱，当年PC1峰在LG8。排除单例AA后正式重复测量标记×年份检验确认四株型性状时间异质性：株高p=0.00647（FDR校正p=0.00862），分蘖长p=0.00116（FDR校正p=0.00232），分蘖径p=0.000228（FDR校正p=0.000912），分蘖数p=0.0191（FDR校正p=0.0191）。LG3信号在2019年强且全基因组显著，但跨两季时间不稳定。
3.6 生物量和生理性状的额外提示关联：生物量和生理性状检测到若干较弱关联，均未超过对应5%全基因组置换阈值，视为次要提示信号。LG12上marker0907（57.272 cM）为单分蘖鲜重和全株鲜重最强BLUE标记（LOD=2.39和1.84；PVE=10.8%和8.5%），也是2019年两性状最强标记，但2018年未一致重现，且部分对秩基逆正态变换敏感，年份依赖明显。因测量鲜重而非干重，收获时短期植物水分状态差异也可能贡献变异，LG12信号被视为初步。生理性状中，LG3上marker0277（5.114 cM）为胞内CO₂浓度和WUE_i最强BLUE标记，LOD约2.4，PVE为10.6%和10.7%。两性状逆正态变换后仍保留最强BLUE峰，但年度特异支持差，因此解释为共享但提示性的生理信号而非独立重复。远端LG9 marker3163有蒸腾速率提示峰，但观察类别高度不平衡（AB=97，BB=3），未进一步解释。总体生物量和生理发现比株型信号更弱且时间一致性更差，需在更大群体和更多环境验证。

讨论部分：4.1 种间芒属群体生物量相关性状表型变异与遗传复杂性。F₁群体显示广泛且多连续变异。早期研究亦表明芒属生物量反映株高、分蘖生产、茎尺寸、发育和环境响应的综合效应。无单性状峰超过5%全基因组置换阈值，但若干解释约8–13%表型变异并在相关性状或数据集重现。模拟检测力显示，模拟5%全基因组阈值为LOD 3.52，接近经验阈值；直接基因分型标记解释10% PVE时仅23.8%模拟检测到，15% PVE为53.9%，20% PVE为80.7%，25% PVE为94.9%；在N=100和现标记集下，约需20% PVE才能达80%检测力。观察到的8–13% PVE峰落在检测力有限范围。该有限检测力不能证明低于阈值峰为真位点。由于模拟假设直接基因分型因果标记，实际检测力可能更低；PVE由同一个小群体估计，单标记R²可能向上偏，应视为近似效应指标。单性状结果被视为提示信号，跨性状、年份或替代表型重复可优先后续，但不单独建立真位点，尤其性状相关时。
4.2 LG9重复提示株型关联簇。LG9近端含最可重复单性状模式。marker2342（11.571 cM）对株高和分蘖长有提示信号，基因型类别差异方向在2018、2019和跨年调整值相似。分蘖数峰在相邻marker2343（10.561 cM），但其AA类仅1个体，未进一步解释。LOD低于全基因组阈值，故不描述为显著QTL。正式标记×年份检验显示marker2342对株高和分蘖长无检测到的互作，而分蘖数有显著互作。LG9模式应更狭义解释为纵向营养生长性状的重复提示证据，而非统一稳定株型位点。因株高、分蘖长和分蘖数是相关株型组分，跨性状重复非独立重复。以往定位研究也报道Miscanthus中株高和发育性状效应在年份或环境间一致性变化。在当前标记密度下，无法区分多效性与连锁，也无法确定marker2342和marker2343标记一个位点还是分开的连锁效应。因此LG9近端被视为需更密定位和独立验证的重复提示关联区域。
4.3 多变量分析揭示LG3强但年份依赖株型关联。marker0580行为不同。四株型性状在PC1正载荷，PC1捕获2018、2019和BLUE数据集中73.20–76.57%变异。2019年marker0580（12.516 cM）达LOD 3.67，高于5%置换阈值3.47。去除单稀有类别后仍显著，主要基因型类别在四组分性状上一致差异。2018年未观察同样模式，BLUE峰未过全基因组阈值。正式重复测量检验经FDR校正确认四株型性状显著标记×年份互作。该时间异质性不能简单解释为季节性环境响应。2018年为建植年，2019年为已建植多年生群体；年降水也不同（2018年388 mm vs. 2019年443 mm），且季节内降雨分布不同。marker0580相关时间差异可能反映建植阶段、株龄、气候变异、年龄依赖表达或其组合。以往研究也报道Miscanthus生长和生物量表现跨生长条件和环境变化。marker0580应视为强2019年特异关联，需在更多年份和环境重复，而非时间稳定选择标记。
4.4 生物量和生理性状提示遗传信号。生物量和生理信号较不令人信服。LG12 marker0907对单分蘖鲜重和全株鲜重显示次要、年份依赖提示信号，主要在2019年，但2018年弱且部分对表型变换敏感。因测鲜重而非干重，这些性状也可能反映收获时短期水分状态。以往全基因组关联研究也在多环境评价的多样芒属种质中研究生物量产量遗传基础。LG3 marker0277为胞内CO₂浓度和WUE_i最强BLUE峰，逆正态变换后位置保留，但年度特异支持差。该杂交在半干旱矿地条件下的转录组研究为这些生理响应提供背景，但未验证标记关联。远端LG9 marker3163蒸腾信号因次要观察类别仅3个体，未在基因型类别水平解释。这些区域仅作为需独立验证的次要提示发现。
4.5 对芒属改良的意义和研究局限。从育种角度看，统计强度和重复性对两个主要区域给出不同图景。LG9近端信号中等，marker2342对株高和分蘖长无检测到的标记×年份互作；marker0580在2019年更强且四株型性状均有显著标记×年份互作。两者均未准备用于标记辅助。研究有技术和实验局限：群体仅100个F₁基因型，461个定位SNP位置间距不均，限制定位分辨率和中等效应位点检测；标记密度、基因型 calling质量和标记顺序是异交芒属群体定位精度的重要决定因素。尽管保留基因型矩阵可报告观察类别构型，原始亲本相位分配和完整上游过滤记录不可用，无法回顾重建完全相间的区间定位分析。因此结果报告为定位标记–性状关联，而非传统区间定位QTL。此外，2018年为建植年，2019年为已建植多年生群体，时间差异不能完全归因于季节性气候变异。复垦田内残余空间异质性也不能排除，因无正式种植前空间调查或精细空间协变量。所有基因型在同一复垦矿地条件下评价，未包括对比农业或非复垦地点；检测关联代表在该地点表达的遗传变异，而非证明赋予矿地特异适应的位点。未来验证应使用更大群体、更密标记、更多年份和对比环境，并受益于重复、标准化的高通量动态株型表型。

结论翻译：这项两年田间研究在复垦矿地地点评价的芒属种间F₁群体中，识别了一小组用于后续跟进的基因组区域。没有单性状峰达到全基因组置换阈值。LG9近端区间显示与株型性状重复提示关联，涉及marker2342对株高和分蘖长以及相邻marker2343对分蘖数；这些单性状信号均未超过相应全基因组阈值。相反，LG3上marker0580在2019年对株型PC1达到全基因组显著，并在排除单例AA标记类别后仍显著，尽管正式标记×年份检验证明时间异质性。生物量和生理关联较弱且重复性较低。因为第一季代表建植，所有基因型仅在一个复垦地点测试，且连锁图谱相对稀疏，这些发现应解释为在该实验条件下表达的标记–性状关系，而非已验证的矿地适应或标记辅助选择位点。需要在更大群体、更密集标记覆盖、更多年份和对比环境中验证，然后才能对LG9和LG3区域进行生物学解析或评估其育种用途。