《Neurobiology of Sleep and Circadian Rhythms》：Molecular correlates of sleep deprivation in the mouse brain identified by meta-analysis of microarray data
编辑推荐：
摘要：研究睡眠剥夺(Sleep Deprivation, SD)后脑内转录水平变化有助于揭示清醒与睡眠状态下不同的分子机制。由于单个研究受限于小样本量，检测差异表达基因(Differentially Expressed Genes, DEGs)的统计效能不足。
摘要：研究睡眠剥夺(Sleep Deprivation, SD)后脑内转录水平变化有助于揭示清醒与睡眠状态下不同的分子机制。由于单个研究受限于小样本量，检测差异表达基因(Differentially Expressed Genes, DEGs)的统计效能不足。本研究通过对已发表的脑表达数据集进行Meta分析，共纳入245只小鼠的173张微阵列(microarray)数据。在FDR校正q<0.01水平下鉴定出498个随睡眠剥夺显著变化的基因，其中96个为原始研究已报道基因，其余402个为本次新发现的候选睡眠基因；在这402个新基因中，14个与人类睡眠性状相关联，3个在基因敲除(knockout, KO)小鼠中表现出睡眠表型。候选基因验证显示Rasd1（亦称Dexras1）在睡眠剥夺后显著上调，且Rasd1 KO小鼠的行为性睡眠时长、分布及睡眠片段(Sleep Bout)结构均发生改变。本研究增进了对睡眠分子关联的理解，并为睡眠研究领域提供了资源数据库。
本文解读对象为 Abdalla OHMH等发表于《Neurobiology of Sleep and Circadian Rhythms》的研究"Molecular correlates of sleep deprivation in the mouse brain identified by meta-analysis of microarray data"。
【研究背景与立题依据】
既往转录组学研究虽已初步揭示睡眠—觉醒转换中能量代谢、应激应答及突触可塑性相关基因的变化，但单个微阵列(microarray)实验因样本量小、平台与处理方法各异，导致鉴定出的睡眠相关基因重现性差，缺乏一致的 pathway层面认知。整合多数据集的转录组Meta分析在基础生物学中应用较少，却有望通过增大统计效能检出低幅值但一致性高的差异表达基因(Differentially Expressed Genes, DEGs)，并衔接小鼠模型与人类睡眠性状全基因组关联分析(Genome-Wide Association Study, GWAS)数据。为此，研究人员开展此项睡眠剥夺(Sleep Deprivation, SD)小鼠脑转录组Meta分析。
【主要关键技术方法】
研究人员从GEO与ArrayExpress筛选并纳入6篇文献、7个数据集、共173张Affymetrix微阵列(小鼠脑皮层、海马、下丘脑、全脑及分选细胞)，涵盖245只C57BL/6J及部分其他品系成年小鼠(8–52周龄)、SD时长3–12 h(仅光照期 gentle handling 或 novel object exposure法)。原始.CEL文件用RMA-sketch归一化，对每对SD vs 睡眠对照计算Hedges' g效应量，采用逆方差加权随机效应模型(Random Effects Model, REM)合并23个配对比较的效应量，FDR(Benjamini–Hochberg)得q值。显著基因(q<0.01)行Gene Ontology(GO)过代表分析(背景限定脑富集基因)；将上/下调基因映射人直系同源基因，用Peak set Enrichment in Gene Sets(PEGS)法与人类睡眠性状GWAS显著SNP区间(±5 kb–2 Mb)交叉富集。独立做3 h SD小鼠(n=8/组)Cerebral Cortex组织qPCR验证8个候选基因(几何均值Arbp/Tbp/Gapdh归一化)。Rasd1 KO小鼠(回交C57BL/6J ≥13代纯合)用被动红外(Passive Infrared, PIR)传感器连续记录基线7天+6 h SD(ZT0–ZT6)后恢复，以连续≥4个10 s无活动bin定义行为睡眠，Split-Plot ANOVA分析基因型×时点交互。
【研究结果】
Transcriptomic meta-analysis identifies new candidate sleep genes
REM Meta分析在q<0.01阈值下鉴定498个SD显著差异基因(219个上调、279个下调)，含96个原研究已报基因(Fos、Cirbp、Homer1、Per2、Usp2等)，另402个为新候选；多数新基因在单个研究中未被检出，因Meta分析提高效能捕获低幅一致变化。所有基因在≥19/23个配对比较中方向一致，86%在≥21/23中一致。
Pathway analysis
对498个显著基因做GO过代表分析(脑富集基因为背景)，富集到5条通路：mRNA加工(mRNA processing)、昼夜节律(circadian rhythm)、蛋白降解调控(regulation of protein catabolic process)、内质网未折叠蛋白反应(endoplasmic reticulum unfolded protein response)及激酶活性负调控(negative regulation of kinase activity)，含11个已知睡眠相关基因。
Identification of new candidate sleep genes
将402个新基因交叉MGI Mammalian Phenotype(MP)、IMPC及人类GWAS Catalogue/Sleep Disorders Knowledge Portal，筛出14个与人类睡眠性状(时型chronotype、失眠insomnia、睡眠时长sleep duration、日间嗜睡daytime sleepiness)相关基因(Ahsa2、Arntl、Bcap29、Eif4b、Lrrk2、Pccb、Pllp、Rasd1、Rbm4、Srsf10、Tmem219、Ubr5、Zbtb16、Zmynd8)；Fisher精确检验显示富集极显著(P=2.15×10??)。依此结合小鼠表型证据选8个深入验证：Arntl(Bmal1)、Fastkd5、Hsf1、Mag、Maoa、Naglu、Rasd1(Dexras1)、Tipin。
Measures of heterogeneity
各基因Hedges' g在23个配对比较中方向高度一致(最低19/23一致)，未见单一研究造成异质性偏移；新候选基因效应量适中但跨研究变异小于经典大效应睡眠基因。
Comparison with Human GWAS Data
PEGS分析示SD上/下调基因与失眠、长/短睡眠时长及一般睡眠时长GWAS显著SNP区间存在8个显著富集窗口(P<0.05)；其中短睡眠时长SNP ±2 Mb窗口基因集GO分析再现了凋亡负调控、未折叠蛋白反应等与本研究通路一致条目。
Validation of new candidate sleep genes
独立3 h SD小鼠皮层qPCR显示8个候选中Rasd1显著上调(约4.6倍，P=0.0003)，其余7个中6个变化方向与Meta分析一致(Binomial P=0.035)，仅Maoa方向微反。证实Meta分析结果可在独立小样本中复现。
Behavioural analysis of Rasd1 KO mice
基线：Rasd1 KO雌鼠ZT1、ZT4睡眠比例升高，ZT20–ZT21降低，伴睡眠片段延长(ZT1)或缩短(ZT20–ZT21)；雄鼠差异不显著。6 h SD后恢复：KO雄鼠SD后首6 h睡眠片段数↓、平均睡眠片段长度↑(更连贯睡眠)，睡眠总时长与野生型无差异；KO雌鼠恢复各指标同野生型。表明Rasd1参与基础睡眠结构维持(雌性为主)及SD后睡眠反弹片段化调节(雄性为主)。
【讨论与结论翻译】
本Meta分析整合245只小鼠173张微阵列，鉴定498个SD差异基因(含402个新候选)及14个人类睡眠性状关联基因，并通过qPCR与Rasd1 KO小鼠行为睡眠验证。Rasd1(又名Dexras1)在SD后上调，其缺失改变基础睡眠时段分布与SD诱导的睡眠反弹片段结构，且具有性别特异性。研究提供迄今最大规模小鼠SD转录组整合资源与新睡眠候选基因目录，支持昼夜与稳态过程在睡眠调控中的交互，为后续EEG/EMG精细睡眠分析与分子机制研究奠定基础。
