果蝇Clock基因PolyQ长度自然变异：地理梯度与功能特性

《Molecular Ecology》：Natural Variation in Drosophila melanogaster ClockPolyQ Length: Geographic Gradients and Functional Properties

　　

摘要专业翻译 昼夜节律钟基因的自然变异为理解生物体如何响应环境异质性提供了强有力的框架。Clock（Clk）基因编码昼夜节律的核心转录调控因子，其含有一个多态性多聚谷氨酰胺（polyQ）区段，该区段的进化意义尚不清楚。本研究整合群体基因组学、行为学及分子分析，探讨黑腹果蝇（Drosophila melanogaster）中Clk polyQ变异的功能与地理分布模式。利用来自127个欧洲种群的数据，研究人员鉴定出11个Clk polyQ等位基因，并发现其频率呈现地理结构化分布，最一致的模式为沿东西向梯度：Q25和Q27等位基因在经度及大陆度（continentality）生物气候轴上呈现稳健的渐变群（cline），而纬度与海拔趋势较弱。近等基因系（near-isogenic lines）的行为学检测显示，polyQ长度在热胁迫条件下调节昼夜节律功能：多数等位基因在不同温度下维持稳定的自由运行周期（free-running period，FRP），而中等长度的Q25等位基因表现出最强（尽管幅度有限）的温度敏感性。在实验室条件下，昼夜节律相位（circadian phase）对升高的温度表现出显著的等位基因特异性敏感性，但相位变异与地理变量之间未呈现一致关系。在分子层面，荧光素酶报告基因检测显示，较长的polyQ等位基因具有更高的转录活性，将polyQ长度与CLK介导的基因表达联系起来。综上所述，这些结果表明Clk polyQ长度的自然变异对昼夜节律调控具有可测量的功能后果，且其地理分布模式与潜在的气候变异相一致，凸显了低复杂度区域（low-complexity regions）在环境梯度下以情境依赖方式调节时钟功能的潜力。

论文解读

一、研究背景与科学问题

昼夜节律是生物体内在的计时系统，通过将生理和行为与24小时环境周期同步，赋予生物体重要的适应性优势。在果蝇中，Clock（Clk）基因编码的CLK蛋白是昼夜节律转录调控网络的核心，与CYCLE（CYC）形成异源二聚体，通过E-box增强子元件激活period（per）和timeless（tim）等下游时钟基因的表达。CLK位于这一转录层级的最顶端，其结构或活性的变异可能影响整个昼夜节律网络。

尽管昼夜节律分子网络已被广泛研究，Clk基因多态性对时钟功能适应性调节的贡献仍不清楚。Clk蛋白C端谷氨酰胺丰富结构域中含有一段由CAG/CAA串联密码子编码的多聚谷氨酰胺（polyQ）区段，该区段在不同物种中均存在长度变异，且常与环境梯度相关。已有研究在鲑鱼、雀鲷科鱼类及鸟类中发现了Clk polyQ长度变异与纬度、迁徙时间或发育时程的关联。然而，在果蝇中，Clk polyQ长度变异是否影响昼夜节律的温度补偿功能、是否呈现与气候环境一致的地理分布格局，此前尚缺乏系统研究。

本研究旨在回答两个核心问题：第一，Clk polyQ长度变异是否影响不同热条件下的昼夜节律表现；第二，该变异在欧洲自然种群中的地理分布是否与潜在的气候变异一致，从而评估其在局部适应中的潜在作用。

二、研究设计与主要技术方法

研究人员利用DrosEU联盟的群体基因组数据，对169个欧洲野生果蝇种群进行Pool-seq测序分析，经质控后保留127个种群，覆盖纬度35.1°至62.6°N、经度?8.4°至38.7°E、海拔11至870米的范围。通过IGV可视化比对读取Clk polyQ区域的插入缺失多态性，鉴定出11个长度变异等位基因（Q21至Q33），并计算各等位基因频率。为分离空间与季节变异，地理和气候关联分析仅使用101个春季采集的种群。利用WorldClim数据库的19个生物气候变量进行主成分分析（PCA），并通过拟二项（quasibinomial）广义线性模型（GLM）检验等位基因频率与地理及气候变量的关联。

在功能验证方面，研究人员从DrosEU同系雌性品系中筛选出六个纯合polyQ等位基因（Q21、Q25、Q29、Q30、Q31、Q33），通过四轮回交将Clk位点渗入到缺陷背景中，构建近等基因系（NILs）。采用果蝇活动监测（DAM）系统记录昼夜节律行为，在24°C和18°C下检测自由运行周期（FRP）和光照-黑暗（LD）条件下的活动相位。分子层面，利用S2细胞进行双荧光素酶报告基因检测，评估不同polyQ长度的CLK变体对per E-box报告基因的转录激活能力。

三、主要研究结果

3.1 Clk polyQ等位基因频率的欧洲空间分析

在127个欧洲种群中鉴定出11个Clk polyQ长度变异等位基因，形成仅以谷氨酰胺密码子数量差异的嵌套长度系列。Q33和Q27是最常见等位基因，几乎出现在所有种群中；Q25和Q31也较常见；Q21和Q32最为稀有。对五个最常见等位基因频率进行PCA分析，前两个主成分解释了58.1%的总方差，种群沿前两轴连续分布而非离散聚类，PC2主要由Q28频率变异驱动。

3.2 Clk polyQ等位基因频率的地理与环境关联

在101个春季采集种群中，经度关联最为稳健：Q25与经度显著正相关（p=0.0013），向东频率增加；Q27与经度显著负相关（p=0.0044），向东频率降低。Q28与纬度呈较弱关联（p=0.036），随纬度升高频率降低。生物气候PCA显示，PC2代表大陆度梯度——高PC2对应冬季温和湿润、温度季节性低的西部海洋性气候，低PC2对应寒冷、季节性强的东部大陆性气候。PC2与经度呈显著负相关。在气候模型中，Q27与PC2显著正相关（OR=1.228，p=0.0004），在温和的西部环境中频率更高；Q25与PC2显著负相关（OR=0.787，p=0.0040），在寒冷、季节性强的东部地区频率更高。PC1（温度-降水梯度）的关联在过离散校正后均不显著。

3.3 昼夜节律运动活动分析

对264只近等基因系果蝇在18°C和24°C下进行行为检测。两因素ANOVA显示，温度、等位基因及其交互作用对FRP均无显著主效应。贝叶斯模型分析表明，Q25在高温下表现出可信的周期缩短（ΔFRP=?0.68小时，95% CI：?1.27至?0.09，P(Δ<0)=0.99），Q31呈中等程度缩短但不确定性较大，其余等位基因的95%可信区间均包含零。总体而言，FRP受温度影响较弱，Q25表现出最强的温度敏感性，提示其温度补偿能力相对降低。

在LD条件下的相位分析中，温度对所有七个预设相位参考均有强烈且一致的影响：降温使傍晚峰提前、早晨峰延后，晨昏间隔缩短。在24°C下，基因型在傍晚起始、傍晚峰、傍晚结束及晨昏间隔等参考上存在显著差异（如傍晚峰F_5,346=9.48，p=1.7×10^?8）。事后比较将基因型分为两组：Q29、Q33和Q25表现出最晚的傍晚峰和最长的晨昏间隔，而Q21和Q30表现出最早的傍晚峰和最短的间隔，Q31居中。在18°C下，各基因型在所有相位参考上均无显著差异，收敛为共同表型。Freedman-Lane置换检验确认了傍晚参考和晨昏间隔的显著基因型×温度交互作用。这些结果表明，Clk polyQ长度变异在温暖条件下产生可重复的傍晚活动相位差异，而在凉爽条件下影响甚微。

3.4 Clk polyQ变异对转录活性的影响

荧光素酶报告基因检测显示，polyQ长度对报告基因表达有显著总体效应（Kruskal-Wallis检验：χ2=14.58，df=3，p=0.002）。较长等位基因polyQ31和polyQ33的转录活性显著高于polyQ21和polyQ28（p<0.05），而polyQ31与polyQ33之间、polyQ21与polyQ28之间无显著差异。这表明较长polyQ区段增强CLK对per E-box报告基因的转录激活能力。

四、讨论与结论

本研究表明，Clk polyQ长度多态性同时影响CLK的分子活性和果蝇昼夜节律行为的热稳定性。等位基因的地理分布，特别是Q25和Q27的显著经度渐变群，与欧洲东西向的大陆度气候梯度一致，但当前数据无法区分选择、种群历史或奠基者效应等机制。polyQ长度与功能之间存在非线性关系：最短（Q21）和最长（Q33）等位基因均表现出良好的温度补偿，而中等长度的Q25最为温度敏感。这种非线性凸显了低复杂度区域通过调节蛋白动态和环境响应性来发挥适应性潜力的能力，无需改变整体蛋白结构即可实现功能微调。

在分子机制层面，polyQ区段可能通过液-液相分离（liquid-liquid phase separation）形成动态无膜凝聚体来调节CLK活性。较长的polyQ区段可能在更宽的热条件下促进并稳定CLK凝聚体的形成，从而支持更强的转录激活，这与Q31和Q33较高的报告基因活性一致。这一机制可能是昼夜节律时钟架构中广泛保守的特征。

Q25等位基因在东部种群中的意外丰富可能反映了额外或间接的选择压力，例如该等位基因的结构灵活性可能对滞育或冷响应等其他温度相关性状有利，或通过杂合优势维持多态性。此外，polyQ区段中CAG/CAA混合密码子的翻译速度差异及相邻结构域的构象影响也可能共同塑造CLK功能。

总体而言，本研究证实了Clk polyQ长度多态性对果蝇昼夜节律系统的分子活性和行为热稳定性均具有可测量的功能后果，其地理分布模式与气候变异一致，为理解低复杂度区域在关键调控蛋白中如何驱动环境适应的进化提供了新视角。