昼夜节律钟适应的调控架构解析：黑腹果蝇全基因组eQTL分析

《Heredity》：Mapping the regulatory architecture of circadian clock adaptation: A genome-wide eQTL analysis in Drosophila melanogaster

　　

昼夜节律钟（circadian clock）使生物体能够将内部日常节律与环境信号对齐，对生存和适应性具有重要影响。尽管黑腹果蝇（Drosophila melanogaster）中这一系统的分子框架已通过涉及十个核心钟基因的转录-翻译反馈环（transcription-translation feedback loops, TTFLs）得以充分表征，但其表达自然变异的遗传基础仍知之甚少。在此，研究人员利用自然表达变异，通过全基因组关联分析（genome-wide association mapping）鉴定表达数量性状基因座（expression quantitative trait loci, eQTLs）。利用果蝇遗传参考面板（Drosophila Genetic Reference Panel, DGRP），研究人员在光照开始后两小时的单一时间点测量了所有核心钟基因的相对表达量。研究人员鉴定了109个显著SNPs（single nucleotide polymorphisms, 单核苷酸多态性）和28个插入缺失（indels）与时钟网络中的表达变异相关。表达水平变化广泛，其中Pdp1ε变异最大（极端品系间差异86倍），cyc变异最小（11.3倍）。仅有三个显著SNPs位于时钟基因自身内部，全部位于Clk，而大多数关联代表位于多种分子功能基因中的反式eQTL（trans-eQTLs）。候选调控因子包括转录因子如Abd-B、tai和E5；RNA结合蛋白包括Pum、Bru-3和Mbl；以及若干长链非编码和反义RNA。在gbb和BMP通路转录因子Mad中也检测到变异。与此一致，Mad敲低降低了vri表达。综上所述，这些结果揭示了昼夜节律基因表达自然变异背后的复杂调控架构。

昼夜节律钟使生物体能够将内部日常节律与环境信号对齐，是进化适应的重要机制。黑腹果蝇（Drosophila melanogaster）的昼夜节律分子框架已被充分解析，其核心是由十个基因构成的转录-翻译反馈环（transcription-translation feedback loops, TTFLs）。然而，这些核心钟基因表达的自然变异如何受遗传调控，此前尚未得到系统研究。为填补这一空白，研究人员利用果蝇遗传参考面板（Drosophila Genetic Reference Panel, DGRP）开展了全基因组eQTL（expression quantitative trait loci, 表达数量性状基因座）分析，以期绘制时钟基因表达的调控架构，并判断自然变异主要源于顺式还是反式调控。

研究基于来自北卡罗来纳州罗利市的120个DGRP近交品系。研究人员在ZT2（光照开始后2小时）收集混合性别成虫头部，提取RNA并逆转录，通过实时定量PCR检测十个核心钟基因的表达，以RpL32为内参计算相对表达量。使用DGRP2软件进行全基因组关联分析，显著性阈值设为p<10^-5，最小等位基因频率（minor allele frequency, MAF）≥0.05。对标准化表达矩阵进行主成分分析（PCA）；利用ReMap 2022 ChIP-seq数据及JASPAR、TFinder预测转录因子结合位点；并通过GAL4/UAS系统在PDF神经元中特异性敲低Mad，结合ANOVA进行表达与行为表型分析。

**时钟基因表达变异与共表达模块**

十个核心钟基因的相对表达量在不同DGRP品系间变化广泛，其中Pdp1ε变异最大（极端品系间差异86倍），cyc变异最小（11.3倍）。相关性分析揭示两个共表达模块：其一包含vri、cwo、Pdp1ε、sgg和dbt，两两Pearson相关系数介于0.53至0.88；另一个独立模块由per和cry构成（r=0.63）；Clk则与所有其他时钟基因几乎不相关（|r|≤0.10）。根据Li和Ji方法估计，十个基因的有效独立性状数为5.82，说明数据保留了充分的独立变异。

**全基因组eQTL鉴定以反式调控为主**

全基因组关联分析共鉴定出109个显著SNPs和28个indels（p<10^-5）。仅有三个显著SNPs位于时钟基因内部，全部位于Clk，其中两个达到或接近Bonferroni校正阈值。绝大多数关联属于反式eQTL（trans-eQTL），定位在远端基因组区域，往往与靶基因不在同一染色体上。许多反式eQTL位于转录因子、RNA结合蛋白或长链非编码RNA（lncRNA）基因附近，例如Abd-B、tai、E5、pumilio（pum）、bru-3、Muscleblind（Mbl）等，提示存在一个广泛分布的调控网络。

**候选转录因子结合与Mad功能验证**

由于部分显著反式eQTL位于转录因子基因内，研究人员利用ReMap 2022 ChIP-seq数据检测到MAD在vri启动子区域存在多个结合位点，并在per启动子中检测到Abd-B结合位点。JASPAR和TFinder分析进一步在高置信度水平预测了vri启动子中的MAD基序、Clk启动子中的E5基序以及per启动子中的多个Abd-B基序。针对Mad的实验表明，在PDF神经元中敲低Mad会使vri表达显著降低约2.1倍（ANOVA, p<1.13×10^-8），但不改变自由运行周期，验证了Mad对vri的调控作用。

**复合表达性状的遗传关联**

对PC1和PC2进行关联分析后，研究人员分别鉴定了57个和174个显著变异。PC1（代表激酶/激活子模块）的位点富集于转录因子（如kn、cad、acj6）、lncRNA和神经元信号基因；PC2（代表per-cry模块）的位点富集于转录和翻译后调控因子（如foxo、mirr、Gcn5、shep）以及神经调控基因。PC1和PC2的信号互不重叠，均不定位到核心钟基因，进一步表明共享的表达变异轴由外部调控因子主导。

讨论部分指出，自然变异对昼夜节律基因表达的调控以分布式反式调控为主，顺式eQTL极少且仅见于Clk。反式eQTL集中于转录因子、RNA结合蛋白和lncRNA，暗示时钟基因表达受一个超越经典时钟网络的广泛调控景观调制。BMP通路中的Mad作为vri的调控因子得到功能验证，为反式eQTL的生物学意义提供了直接证据。研究还发现，同一基因座内部存在多个关联变异的聚类现象，结合DGRP中连锁不平衡的快速衰减，提示这些关联可能反映等位基因异质性。与以往eQTL研究相比，个体变异很少重复，但反式调控为主的整体架构具有可重复性。局限性包括单次时间点采样、单生物重复、未控制性别等，但基因型-性别互作分析表明大部分时钟基因的品系间变异在性别间共享。

综上所述，研究结论为：这些结果揭示了昼夜节律基因表达自然变异背后的复杂调控架构。