基因组结构预测印度谷螟（Plodia interpunctella）对颗粒体病毒（PiGV）的抗性且未检测到生活史权衡

《Molecular Ecology》：Genomic Structure Predicts Granulovirus Resistance Without a Detectable Life-History Trade-Off in Plodia interpunctella

　　

理解病原体抗性的遗传结构及其与宿主生活史性状的关系是进化生物学中的一个基本问题，然而抗性成本在何种条件下表现为可检测的权衡仍知之甚少。在此，研究人员检验了印度谷螟（Plodia interpunctella，Indian meal moth）对其颗粒体病毒（granulovirus，PiGV）的抗性，涉及12个在无有意病原体暴露下维持的近交系（inbred lines），并以一个独立组装的远交原种（outbred stock population）作为基因组参考。研究人员使用来自剂量–反应测定（dose–response assays）的LD₅₀值量化抗性，使用由蛹重、发育时间和生长速率推导的复合主成分（principal component，PC）轴总结发育生活史变异，并评估基因组结构以及与抗性的窗口水平关联。抗性在近交系间约变化两倍。复合生活史轴以及任何单个生活史性状均与LD₅₀无显著相关。相反，基因组PC1，即编码序列（coding sequence，CDS）限制性基因组分化的主要轴，与抗性显著相关（r = ?0.65，p = 0.022），而总体成对距离则不然。没有单个PCA（principal component analysis）–LD₅₀窗口通过多重检验校正；相反，抗性相关信号广泛分布于整个基因组，并且跨五条染色体的趋同信号未显示对具有已知免疫功能基因的富集。该模式与抗性的弥散、潜在多基因基础一致，尽管有限数量的系限制了位点水平推断。结合先前证明该系统中抗性–生长权衡的工作，研究人员的结果表明，抗性可以在基因组上结构化，而不与本文测量的发育生活史变异发生可检测耦合，并且抗性成本的表达取决于遗传和环境背景。

病原体抗性的遗传结构及其与宿主生活史性状的关系，是进化生物学的基本问题之一。抗性通常伴随适合度成本，如生长速率降低、发育时间延长或繁殖受损，从而在病原压力下影响适应。尽管在植物和昆虫宿主—病原系统中已反复记录这类权衡，但抗性成本在何种条件下表现为可检测权衡，以及其遗传结构如何，仍了解不足。资源分配理论预测免疫成本在资源真正受限时更明显；印度谷螟（Plodia interpunctella）—颗粒体病毒（granulovirus，PiGV）系统中的既往研究显示，抗性—生长速率相关依赖营养环境，资源质量、历史偶然性和选择时间尺度均影响权衡是否可检测。Bartlett et al.还发现直接选择发育时间未产生相应抗性变化，提示关系依赖选择方向。该宿主—病原对因世代短、可近交、感染终点明确而适合解析抗性遗传基础；实验进化和近交系工作已证明抗性—发育时间/生长权衡具有遗传性和组成性，可在无有意病原暴露种群中表达。然而，该物种尚无全基因组水平的抗性变异表征。近期染色体水平参考基因组使研究人员能够在表型和进化背景已知的系统中检查抗性遗传结构。为此，研究人员在12个无有意病原暴露维持的近交系中检验PiGV抗性，并以独立组装的远交原种（outbred stock population）作为基因组参考。核心结论是：基因组PC1（principal component 1，主成分1）与抗性显著相关，而复合生活史轴不相关；抗性相关信号弥散分布于基因组，跨五条染色体的趋同信号未富集已知免疫基因，提示弥散、潜在多基因基础且未检测到与所测发育生活史变异的耦合。

关键技术方法：12个近交系源自共同远交实验室原种，曾经历5代严格全同胞交配，研究人员再实施3代严格全同胞交配（每系5个平行重复并传代最可育重复），表型测定在G5–G13；独立远交stock来自加州Oakland和Berkeley两地家庭侵染并合并维持，作为基因组参考。感染测定采用PiGV系列稀释（1×至10^?3）与无病毒对照，每种群325头幼虫，拟合logistic回归估计LD₅₀。生活史测蛹重、发育时间和生长速率并做PCA（principal component analysis，主成分分析）。基因组用NovaSeq X Plus测序，每系/stock为2–3个体池DNA，比对参考基因组GCF_027563975.1；开展CDS（coding sequence，编码序列）限制PCA、ROH（runs of homozygosity，纯合片段）量化、250 kb窗口局部PCA和WZA（weighted Z-score analysis，加权Z分数分析），并行Mantel检验。

3.1 近交系间基因组结构。通过邻居连接树、Mantel检验、CDS限制PCA、全基因组变体稳健性检验、合并stock PCA和ROH量化，结果表明近交系间存在明显分化；总体基因组距离与LD₅₀成对差异不显著相关（r=0.16，p=0.14），但基因组PC1与LD₅₀显著负相关（r=?0.65，p=0.022）。该结果在基因组范围变体中方向一致，stock与近交系沿PC1清晰分离并形成两个与双来源一致的簇。近交系ROH覆盖均值26.1%（14.4%–32.1%），stock均值7.3%（1.3%–26.6%），组间差异显著；ROH覆盖与LD₅₀正相关（r=0.68，p=0.016），但因与基因组PC1混杂，未解释为纯合本身增加抗性。

3.2 抗性变异及其遗传基础。通过logistic广义线性模型（generalized linear model，GLM）、每群LD₅₀估计、窗口局部PCA和WZA分析，发现剂量是感染主导预测因子（χ²=1267.24，df=1，p<0.001），种群（χ²=40.28，df=12，p<0.001）及剂量×种群交互（χ²=23.12，df=12，p=0.027）显著。LD₅₀从IL-4的0.439到IL-24的0.826，约两倍范围；基因组PC1负相关，排除IL-4后减弱但未消失（r=?0.60，p=0.052），等位频率标准化PCA相似（r=?0.58，p=0.049）。2164个窗口均未通过Bonferroni校正，最小p=0.000458，阈值2.31×10^?5；PVE（proportion of variance explained，解释方差比例）分布显示top 1%窗口贡献4.7%，top 10%为29.3%，top 25%为55.2%，观察累积PVE未超过95%零分布包络并低于下界，提示信号分散。21个窗口按PVE和WZA联合排名形成10个连续区域，跨五条染色体；这些区域未富集已知免疫基因（严格集观察=1，零均值=1.08，p=0.61；宽泛集观察=8，零均值=5.96，p=0.27），故不提名候选基因。

3.3 种群间生活史变异。通过单因素方差分析（analysis of variance，ANOVA）、Tukey HSD、生活史PCA和羽化率分析，发现蛹重（F_12,770=23.89，p<0.001）、发育时间（F_12,770=26.03，p<0.001）和生长速率（F_12,770=9.63，p<0.001）在种群间差异显著。生活史PC1解释67.8%方差，三性状均负载荷；蛹重与发育时间正相关（r=0.74，p=0.004）。近交系羽化率均值0.77（0.61–0.85），stock为0.87；羽化率在近交系间无显著差异（χ²=11.32，df=11，p=0.42），与LD₅₀无显著相关（r=?0.49，p=0.11）。

3.4 生活史与抗性关系。通过Pearson相关和k均值聚类，发现生活史PC1与LD₅₀不显著相关（r=0.23，p=0.47），单个性状亦然：生长速率r=?0.11，p=0.73；蛹重r=?0.28，p=0.35；发育时间r=?0.32，p=0.29。生活史PC1与基因组PC1（r=?0.37，p=0.235）和PC2（r=?0.08，p=0.80）均不显著相关。k=3聚类得到高抗n=3、低抗n=5、中抗n=4；簇内生长速率与LD₅₀关系不显著，低抗簇接近显著（r=?0.88，p=0.051，n=5），被视为假设生成。总体上未检测到一致抗性—生活史关系。

讨论总结：讨论认为，基因组PC1与LD₅₀的关联反映特定基因组分化轴，而非总体基因组分歧；ROH与LD₅₀正相关使基因组结构和纯合性贡献难以分离。n=12、测定日与种群混杂、池DNA导致的个体内/个体间杂合性模糊、祖先历史不完整、注释不完整和250 kb窗口含多基因均限制位点推断。窗口分析未支持少数主效位点，更符合弥散、潜在多基因结构，但不能排除低功效。若确认，该结构更兼容定量共进化框架而非简单基因对基因模型。抗性与生活史不沿同一主要分化轴，但本研究不解析位点水平多效性。讨论还提到围食膜基质和RNAi（RNA interference，RNA干扰）可能参与抗性，而窗口关联不易检测。当代stock与近交系基因组不同，且非祖先种群，不能推断祖先等位频率；stock的LD₅₀=0.5185位于近交系较低范围。无生活史相关不意味着存在“Darwinian demons”，既往工作表明权衡依赖营养、环境和选择方向，本研究系与Bartlett et al.的系不同但同祖先，提示可检测性依赖遗传背景和繁殖历史。资源丰富条件可能掩盖成本，但Bartlett et al.在相近饮食下检测到权衡，故资源并非唯一原因。抗性变异可通过空间或时间异质选择维持而不必在所有条件下表现固定权衡。未来需个体测序、更大系系和更接近自然环境的条件。

研究结论部分翻译：本研究扩展了该系统中Bartlett et al.、Boots和Begon等既往工作，提供了印度谷螟抗性变异的首个全基因组表征。未来检验这些近交系基因型在可变生态条件（包括资源限制、温度波动和竞争压力）下的表现，将有助于确定抗性生活史权衡是否在更接近自然环境的条件下出现。该系统未来基因组工作也将受益于个体水平测序而非池提取，以消除上述个体内/个体间杂合性模糊，并受益于额外世代有记录的同胞配对交配，以更充分地将各系近交历史与其共享源种群区分开。