果蝇寿命性别差异中Y染色体效应的重新审视：基于三种果蝇转座元件丰度与异染色质多态性的比较研究

《Aging Cell》：Exploring the Role of Transposable Elements to Sex Gap in Longevity in Three Drosophila Species

　　

在果蝇（Drosophila）中，与其他动物物种类似，雌性往往比雄性寿命更长，这一现象被称为寿命性别差异（sex gap in longevity, SGL）。一种可能的解释涉及转座元件（transposable elements, TE）的活性。TE通常被表观遗传机制（epigenetic mechanisms）抑制，但这种调控随年龄增长而减弱。Y染色体富含TE，因此与年龄相关的TE活性在年长雄性中可能更为显著，可能影响寿命（毒性Y染色体假说，toxic Y hypothesis）。利用三种不同TE含量的果蝇物种，研究人员表明SGL在野生型种群中变异，反映了种内和种间差异。转录组数据显示，在黑腹果蝇（D. melanogaster）和铃木果蝇（D. suzukii）中，TE来源的reads随年龄增长而增加，但在拟暗果蝇（D. simulans）中未观察到这一现象，且雌雄之间无差异。然而，研究人员发现所测试的三种果蝇物种中，年长雄性与年长雌性相比，上调的TE拷贝数量更高，但这不一定反映转座活性增加。此外，寿命检测未显示TE表达增加与寿命缩短之间存在关联。最后，在不同SGL的品系之间替换Y染色体导致黑腹果蝇中雄性寿命逐渐且品系特异性地缩短，尽管Y染色体大小相似。这伴随基因表达变化和TE表达增加。总之，研究人员的工作呼吁重新审视毒性Y染色体假说，表明它在基因组TE含量低的物种中可能不太相关。值得注意的是，Y染色体对寿命的效应似乎更多地由Y连锁的异染色质多态性（Y-linked heterochromatic polymorphisms）和上位性相互作用（epistatic interactions）驱动，而非TE的总体丰度。

在多数动物物种中，不同性别表现出不同的死亡率模式，导致寿命差异，即寿命性别差异（sex gap in longevity, SGL）。针对该现象，研究人员已提出性选择假说（sexual selection hypothesis）和母系遗传诅咒假说（mother's curse hypothesis）等解释。此外，异配子性别（heterogametic sex）通常寿命较短，机制可能涉及单X染色体的隐性有害等位基因无条件表达，即未保护X效应（unguarded X effect）。毒性Y染色体假说（toxic Y hypothesis）进一步提出，因重组抑制而累积大量转座元件（transposable elements, TE）的Y染色体（或W染色体）可能在衰老过程中释放TE，损害异配性别寿命。TE通常受表观遗传机制（epigenetic mechanisms）（包括DNA甲基化、组蛋白修饰及小RNA）抑制，但随年龄增长，沉默机制减弱，TE转录与转座活性增加。已有研究显示，D. melanogaster年长雄性较雌性更快丢失H₃K₉me₂抑制性组蛋白修饰，导致TE表达上升；异常性染色体构型亦提示Y染色体数量影响生存。然而，近期利用CRISPR/Cas9技术构建不同大小Y染色体的研究表明，Y染色体大小及异染色质（heterochromatin）与TE数量并不导致寿命性别差异；D. pseudoobscura中Y连锁TE变异仅引起TE表达轻微上升，未加速衰老。因此，Y染色体是否通过TE丰度及增龄性TE激活影响雄性寿命，仍是争议焦点。

为厘清SGL的自然变异并验证毒性Y染色体假说，研究人员选取三种TE含量迥异的果蝇物种：黑腹果蝇（D. melanogaster，TE含量20%）、拟暗果蝇（D. simulans，12%）和铃木果蝇（D. suzukii，47%），开展了寿命检测、转录组分析及Y染色体替换回交研究。研究证实SGL在野生型种群中呈现显著种内与种间变异；TE表达随年龄上调，但物种间及性别间存在异质性；Y染色体替换实验揭示Y染色体对寿命的影响更可能由异染色质多态性和上位性互作（epistatic interactions）驱动，而非总体TE丰度。这些发现修正了毒性Y染色体假说，强调Y染色体效应具有物种和品系特异性。

方法上，研究人员主要采用三项技术。第一，寿命检测（longevity assays）：对三个果蝇物种的28至132个野生型隔离系（isofemale strains）逐日追踪生存率，以计算中位寿命和SGL值；并对Y染色体替换品系开展长达10代或5代的寿命追踪。第二，RNA测序（RNA-seq）：对每物种两个近交系（inbred strains）的年轻及年长果蝇（去除性腺后的躯体）进行转录组测序。第三，Y染色体替换回交（Y chromosome replacement introgression）：利用果蝇雄性不重组特性，在不同品系间连续回交以替换Y染色体，并结合流式细胞术（flow cytometry）估算基因组大小。

**研究结果**

*2.1 种群和物种间SGL存在差异*
通过277个野生型隔离系的寿命检测，研究人员发现59%的品系呈显著负SGL，仅5%呈显著正值。D. melanogaster、D. simulans和D. suzukii中分别有58%、65%和35%的品系表现雄性寿命缩短，平均分别缩短10.5天、12.5天和9天。贝叶斯层次模型分析表明，种群来源解释总变异的38.2%，物种效应解释8.6%，证实SGL存在显著种内与种间变异。

*2.2 D. melanogaster和D. suzukii随年龄增长高丰度表达TE家族*
利用TEtools分析发现，年老果蝇TE来源reads百分比随年龄显著增加。在家族水平上，D. melanogaster和D. suzukii多数TE家族随年龄上调，而D. simulans因TE含量较低未现统计显著差异。不同性别或品系间差异表达TE家族重叠极少，提示衰老相关TE激活并非家族特异的。

*2.3 年龄相关的TE来源reads和表达TE家族增加不受性别影响*
广义线性模型（GLM）分析显示，年龄对TE来源reads百分比影响显著，但性别及年龄与性别交互作用均不显著。TE家族层面，除少数品系外，年长雄性与雌性在差异表达家族数量上无显著差异，表明衰老驱动的TE表达上升在雌雄间整体相似。

*2.4 D. melanogaster和D. suzukii随年龄增长高丰度表达TE拷贝*
采用Telescope进行拷贝特异性分析，研究人员发现D. melanogaster和D. suzukii中随年龄上调的TE拷贝显著多于下调拷贝；D. simulans中未观察到此趋势。部分在雄性中上调的TE拷贝在雌性中也存在，提示去抑制模式不仅限于Y染色体。

*2.5 年长雄性与年长雌性相比有更多TE拷贝表达*
拷贝特异性分析显示，所测三个物种中，年长雄性较年长雌性均检测到更多上调TE拷贝，该差异在年轻个体中已部分显现。全长度TE拷贝（full-length TE copies）在年长雄性上调TE中显著富集。置换检验未发现DETEs在差异表达基因（DEGs）附近（±1 kb）存在显著富集或耗竭。

*2.6 TE表达增加与雄性寿命缩短无稳定关联*
尽管部分品系中观察到雄性寿命缩短及年长雄性DETEs增加，但D. suzukii未出现显著性别寿命差异。各物种内SGL值与基因组TE含量均无显著相关性，说明整体TE丰度与寿命变异无直接联系。

*2.7 D. melanogaster中Y染色体替换与SGL增大相关*
在D. melanogaster的dmgoth101和dmsj7品系间进行Y染色体替换回交，无论A组还是B组杂交，均随代进程出现雄性寿命逐渐下降，SGL绝对值增大；而dmgoth101与dmgoth63回交未现此现象。在D. simulans中，SGL值随世代波动，未出现一致下降，表明Y染色体替换效应具品系和物种特异性。

*2.8 雄性寿命与Y染色体大小无关*
流式细胞术分析显示，D. melanogaster回交后代与亲本间基因组大小无显著差异，提示替换前后Y染色体大小相似；D. simulans中部分组合虽检测到基因组大小差异，但未出现渐进性雄性寿命缩短，进一步说明Y染色体大小并非决定雄性寿命的关键因素。

*2.9 Y染色体替换与基因表达全局变化和TE表达增加相关*
RNA-seq显示，F10代相对于F1代出现大量DEGs，尤以雄性为甚。基因本体（GO）富集分析表明这些DEGs涉及转录调控、染色质组织等。TE表达分析表明，F10代中TE来源reads百分比显著高于F1代，且雄性始终高于雌性，提示Y染色体替换伴随全局转录组变化及TE表达上升。

**讨论与结论**
研究人员在讨论中指出，本研究证实SGL是野生型果蝇种群的变异性状。拷贝特异性分析揭示TE拷贝随年龄上调，符合异染色质去抑制预期；但年长雄性和雌性间TE表达的整体差异在家族层面被稀释，仅在拷贝层面显现。研究未能证实TE表达增加与雄性寿命缩短之间存在普遍关联，也未发现SGL与基因组TE含量的关系。鉴于Y染色体替换实验中Y染色体大小相似但雄性寿命仍下降，且效应具品系特异性，研究人员推断Y连锁TE丰度不足以解释观察结果；Y染色体与基因组的互作及Y连锁异染色质多态性可能是更重要的驱动因素。这提示毒性Y染色体假说在基因组TE含量较低的物种中相关性可能降低。

研究结论：SGL在野生型果蝇中变异显著，表现为种群和物种水平的显著差异；Y染色体替换可导致D. melanogaster部分品系中雄性寿命渐进性下降，且该效应与Y染色体大小和总体TE丰度无关；Y染色体对寿命的影响更可能由Y连锁异染色质多态性和上位性互作驱动。因此，毒性Y染色体假说亟需重新审视。