A systems approach identifies networks and genes linking sleep and stress: implications for neuropsychiatric disorders

Multi-omic molecular profiling and network biology for precision anaesthesiology: a narrative review

Technological advancement, data democratisation, and decreasing costs have led to a revolution in molecular biology in which the entire set of DNA, RNA, proteins, and various other molecules - the 'multi-omic' profile - can be measured in humans. Sequencing 1 million bases of human DNA now costs US$0.01, and emerging technologies soon promise to reduce the cost of sequencing the whole genome to US$100. These trends have made it feasible to sample the multi-omic profile of millions of people, much of which is publicly available for medical research. Can anaesthesiologists use these data to improve patient care? This narrative review brings together a rapidly growing literature in multi-omic profiling across numerous fields that points to the future of precision anaesthesiology. Here, we discuss how DNA, RNA, proteins, and other molecules interact in molecular networks that can be used for preoperative risk stratification, intraoperative optimisation, and postoperative monitoring. This literature provides evidence for four fundamental insights: (1) Clinically similar patients have different molecular profiles and, as a consequence, different outcomes. (2) Vast, publicly available, and rapidly growing molecular datasets have been generated in chronic disease patients and can be repurposed to estimate perioperative risk. (3) Multi-omic networks are altered in the perioperative period and influence postoperative outcomes. (4) Multi-omic networks can serve as empirical, molecular measurements of a successful postoperative course. With this burgeoning universe of molecular data, the anaesthesiologist-of-the-future will tailor their clinical management to an individual's multi-omic profile to optimise postoperative outcomes and long-term health.

Transcriptomic profiles of stress susceptibility and resilience in the amygdala and hippocampus

A single, severe episode of stress can bring about myriad responses amongst individuals, ranging from cognitive enhancement to debilitating and persistent anxiety; however, the biological mechanisms that contribute to resilience versus susceptibility to stress are poorly understood. The dentate gyrus (DG) of the hippocampus and the basolateral nucleus of the amygdala (BLA) are key limbic regions that are susceptible to the neural and hormonal effects of stress. Previous work has also shown that these regions contribute to individual variability in stress responses; however, the molecular mechanisms underlying the role of these regions in susceptibility and resilience are unknown. In this study, we profiled the transcriptomic signatures of the DG and BLA of rats with divergent behavioral outcomes after a single, severe stressor. We subjected rats to three hours of immobilization with exposure to fox urine and conducted a behavioral battery one week after stress to identify animals that showed persistent, high anxiety-like behavior. We then conducted bulk RNA sequencing of the DG and BLA from susceptible, resilient, and unexposed control rats. Differential gene expression analyses revealed that the molecular signatures separating each of the three groups were distinct and non-overlapping between the DG and BLA. In the amygdala, key genes associated with insulin and hormonal signaling corresponded with vulnerability. Specifically, Inhbb, Rab31 , and Ncoa3 were upregulated in the amygdala of stress-susceptible animals compared to resilient animals. In the hippocampus, increased expression of Cartpt - which encodes a key neuropeptide involved in reward, reinforcement, and stress responses - was strongly correlated with vulnerability to anxiety-like behavior. However, few other genes distinguished stress-susceptible animals from control animals, while a larger number of genes separated stress-resilient animals from control and stress-susceptible animals. Of these, Rnf112, Tbx19 , and UBALD1 distinguished resilient animals from both control and susceptible animals and were downregulated in resilience, suggesting that an active molecular response in the hippocampus facilitates protection from the long-term consequences of severe stress. These results provide novel insight into the mechanisms that bring about individual variability in the behavioral responses to stress and provide new targets for the advancement of therapies for stress-induced neuropsychiatric disorders.

A Nile rat transcriptomic landscape across 22 organs by ultra-deep sequencing and comparative RNA-seq pipeline (CRSP)

RNA sequencing (RNA-seq) has been a widely used high-throughput method to characterize transcriptomic dynamics spatiotemporally. However, RNA-seq data analysis pipelines typically depend on either a sequenced genome and/or corresponding reference transcripts. This limitation is a challenge for species lacking sequenced genomes and corresponding reference transcripts. The Nile rat (Arvicanthis niloticus) has two key features - it is daytime active, and it is prone to diet-induced diabetes, which makes it more similar to humans than regular laboratory rodents. However, at the time of this study, neither a Nile rat genome nor a reference transcript set were available, making it technically challenging to perform large-scale RNA-seq based transcriptomic studies. This genome-independent work progressed concurrently with our generation of a Nile rat genome. A well-annotated genome requires several iterations of manually reviewing curated transcripts and takes years to achieve. Here, we developed a Comparative RNA-Seq Pipeline (CRSP), integrating a comparative species strategy independent of a specific sequenced genome or species-matched reference transcripts. We performed benchmarking to validate that our CRSP tool can accurately quantify gene expression levels. In this study, we generated the first ultra-deep (2.3 billion × 2 paired-end) Nile rat RNA-seq data from 59 biopsy samples representing 22 major organs, providing a unique resource and spatial gene expression reference for Nile rat researchers. Importantly, CRSP is not limited to the Nile rat species and can be applied to any species without prior genomic knowledge. To facilitate a general use of CRSP, we also characterized the number of RNA-seq reads required for accurate estimation via simulation studies. CRSP and documents are available at: https://github.com/pjiang1105/CRSP.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Comparative analysis of hippocampal transcriptional features between major depressive disorder patients and animal models

BACKGROUND

Major depressive disorder (MDD) is a psychiatric disorder caused by various etiologies. Chronic stress models are used to simulate the heterogeneous pathogenic processes of depression. However, few studies have compared transcriptional features between stress models and MDD patients.

METHODS

We generated hippocampal transcriptional profiles of the chronic social defeat model by RNA sequencing and downloaded raw data of the same brain region from public databases of the chronic unpredictable mild stress model, the learned helplessness model, and MDD patients. Differential expression and gene co-expression analyses were integrated to compare transcriptional features between stress models and MDD patients.

RESULTS

Each stress model shared 11.4% to 16.3% of differentially expressed genes with MDD patients. Functional analysis at the gene expression level identified altered ensheathment of neurons in both stress models and MDD patients. At the gene network level, each stress model shared 20.9% to 41.6% of co-expressed genes with MDD patients. Functional analysis based on these genes found that axon guidance signaling is the most significantly enriched pathway that was shared by all stress models and MDD patients.

LIMITATIONS

This study was limited by considering only a single brain region and a single sex of stress model animals.

CONCLUSIONS

Our results show that hippocampal transcriptional features of stress models partially overlap with those of MDD patients. The canonical pathways of MDD patients, including ensheathment of neurons, PTEN signaling, and axonal guidance signaling, were shared with all stress models. Our findings provide further clues to understand the molecular mechanisms of depression.

NREM delta power and AD-relevant tauopathy are associated with shared cortical gene networks

Reduced NREM sleep in humans is associated with AD neuropathology. Recent work has demonstrated a reduction in NREM sleep in preclinical AD, pointing to its potential utility as an early marker of dementia. We test the hypothesis that reduced NREM delta power and increased tauopathy are associated with shared underlying cortical molecular networks in preclinical AD. We integrate multi-omics data from two extensive public resources, a human Alzheimer's disease cohort from the Mount Sinai Brain Bank (N = 125) reflecting AD progression and a (C57BL/6J × 129S1/SvImJ) F2 mouse population in which NREM delta power was measured (N = 98). Two cortical gene networks, including a CLOCK-dependent circadian network, are associated with NREM delta power and AD tauopathy progression. These networks were validated in independent mouse and human cohorts. Identifying gene networks related to preclinical AD elucidate possible mechanisms associated with the early disease phase and potential targets to alter the disease course.

Sleep deficits and cannabis use behaviors: an analysis of shared genetics using linkage disequilibrium score regression and polygenic risk prediction

STUDY OBJECTIVES

Estimate the genetic relationship of cannabis use with sleep deficits and an eveningness chronotype.

METHODS

We used linkage disequilibrium score regression (LDSC) to analyze genetic correlations between sleep deficits and cannabis use behaviors. Secondly, we generated sleep deficit polygenic risk score (PRS) and estimated their ability to predict cannabis use behaviors using linear and logistic regression. Summary statistics came from existing genome-wide association studies of European ancestry that were focused on sleep duration, insomnia, chronotype, lifetime cannabis use, and cannabis use disorder (CUD). A target sample for PRS prediction consisted of high-risk participants and participants from twin/family community-based studies (European ancestry; n = 760, male = 64%; mean age = 26.78 years). Target data consisted of self-reported sleep (sleep duration, feeling tired, and taking naps) and cannabis use behaviors (lifetime ever use, number of lifetime uses, past 180-day use, age of first use, and lifetime CUD symptoms).

RESULTS

Significant genetic correlation between lifetime cannabis use and an eveningness chronotype (rG = 0.24, p < 0.001), as well as between CUD and both short sleep duration (<7 h; rG = 0.23, p = 0.017) and insomnia (rG = 0.20, p = 0.020). Insomnia PRS predicted earlier age of first cannabis use (OR = 0.92, p = 0.036) and increased lifetime CUD symptom count (OR = 1.09, p = 0.012).

CONCLUSION

Cannabis use is genetically associated with both sleep deficits and an eveningness chronotype, suggesting that there are genes that predispose individuals to both cannabis use and sleep deficits.

Parallel metabolomics and lipidomics enables the comprehensive study of mouse brain regional metabolite and lipid patterns

Global profiling of the metabolome and lipidome of specific brain regions is essential to understanding the cellular and molecular mechanisms regulating brain activity. Given the limited amount of starting material, conventional mouse studies comparing brain regions have mainly targeted a set of known metabolites in large brain regions (e.g., cerebrum, cortex). In this work, we developed a multimodal analytical pipeline enabling parallel analyses of metabolomic and lipidomic profiles from anatomically distinct mouse brain regions starting with less than 0.2 mg of protein content. This analytical pipeline is composed of (1) sonication-based tissue homogenization, (2) parallel metabolite and lipid extraction, (3) BCA-based sample normalization, (4) ultrahigh performance liquid chromatography-mass spectrometry-based multimodal metabolome and lipidome profiling, (5) streamlined data processing, and (6) chord plot-based data visualization. We applied this pipeline to the study of four brain regions in males including the amygdala, dorsal hippocampus, nucleus accumbens and ventral tegmental area. With this novel approach, we detected over 5000 metabolic and 6000 lipid features, among which 134 metabolites and 479 lipids were directly confirmed via automated MS² spectral matching. Interestingly, our analysis identified unique metabolic and lipid profiles in each brain regions. Furthermore, we identified functional relationships amongst metabolic and lipid subclasses, potentially underlying cellular and functional differences across all four brain regions. Overall, our novel workflow generates comprehensive region-specific metabolomic and lipidomic profiles using very low amount of brain sub-regional tissue sample, which could be readily integrated with region-specific genomic, transcriptomic, and proteomic data to reveal novel insights into the molecular mechanisms underlying the activity of distinct brain regions.

Shared Transcriptional Signatures in Major Depressive Disorder and Mouse Chronic Stress Models

BACKGROUND

Most of our knowledge of the biological basis of major depressive disorder (MDD) is derived from studies of chronic stress models in rodents. While these models capture certain aspects of the behavioral and neuroendocrine features of MDD, the extent to which they reproduce the molecular pathology of the human syndrome remains unknown.

METHODS

We systematically compared transcriptional signatures in two brain regions implicated in depression-medial prefrontal cortex and nucleus accumbens-of humans with MDD and of 3 chronic stress models in mice: chronic variable stress, adult social isolation, and chronic social defeat stress. We used differential expression analysis combined with weighted gene coexpression network analysis to create interspecies gene networks and assess the capacity of each stress paradigm to recapitulate the transcriptional organization of gene networks in human MDD.

RESULTS

Our results show significant overlap between transcriptional alterations in medial prefrontal cortex and nucleus accumbens in human MDD and the 3 mouse chronic stress models, with each of the chronic stress paradigms capturing distinct aspects of MDD abnormalities. Chronic variable stress and adult social isolation better reproduce differentially expressed genes, while chronic social defeat stress and adult social isolation better reproduce gene networks characteristic of human MDD. We also identified several gene networks and their constituent genes that are most significantly associated with human MDD and mouse stress models.

CONCLUSIONS

This study demonstrates the ability of 3 chronic stress models in mice to recapitulate distinct aspects of the broad molecular pathology of human MDD, with no one mouse model apparently better than another.

Hybrid analysis of gene dynamics predicts context-specific expression and offers regulatory insights

Motivation

To understand the regulatory pathways underlying diseases, studies often investigate the differential gene expression between genetically or chemically differing cell populations. Differential expression analysis identifies global changes in transcription and enables the inference of functional roles of applied perturbations. This approach has transformed the discovery of genetic drivers of disease and possible therapies. However, differential expression analysis does not provide quantitative predictions of gene expression in untested conditions. We present a hybrid approach, termed Differential Expression in Python (DiffExPy), that uniquely combines discrete, differential expression analysis with in silico differential equation simulations to yield accurate, quantitative predictions of gene expression from time-series data.

Results

To demonstrate the distinct insight provided by DiffExpy, we applied it to published, in vitro, time-series RNA-seq data from several genetic PI3K/PTEN variants of MCF10a cells stimulated with epidermal growth factor. DiffExPy proposed ensembles of several minimal differential equation systems for each differentially expressed gene. These systems provide quantitative models of expression for several previously uncharacterized genes and uncover new regulation by the PI3K/PTEN pathways. We validated model predictions on expression data from conditions that were not used for model training. Our discrete, differential expression analysis also identified SUZ12 and FOXA1 as possible regulators of specific groups of genes that exhibit late changes in expression. Our work reveals how DiffExPy generates quantitatively predictive models with testable, biological hypotheses from time-series expression data.

Availability and implementation

DiffExPy is available on GitHub (https://github.com/bagherilab/diffexpy).

Supplementary information

Supplementary data are available at Bioinformatics online.

Onset of regular cannabis use and young adult insomnia: an analysis of shared genetic liability

STUDY OBJECTIVES

Estimate the genetic and environmental influences on the relationship between onset of regular cannabis use and young adult insomnia.

METHODS

In a population-based twin cohort of 1882 twins (56% female, mean age = 22.99, SD = 2.97) we explored the genetic/environmental etiology of the relationship between onset of regular cannabis use and insomnia-related outcomes via multivariate twin models.

RESULTS

Controlling for sex, current depression symptoms, and prior diagnosis of an anxiety or depression disorder, adult twins who reported early onset for regular cannabis use (age 17 or younger) were more likely to have insomnia (β = 0.07, p = 0.024) and insomnia with short sleep on weekdays (β = 0.08, p = 0.003) as young adults. We found significant genetic contributions for the onset of regular cannabis use (a2 = 76%, p < 0.001), insomnia (a2 = 44%, p < 0.001), and insomnia with short sleep on weekdays (a2 = 37%, p < 0.001). We found significant genetic correlations between onset of regular use and both insomnia (rA = 0.20, p = 0.047) and insomnia with short sleep on weekdays (rA = 0.25, p = 0.008) but no significant environmental associations between these traits.

CONCLUSIONS

We found common genetic liabilities for early onset of regular cannabis use and insomnia, implying pleiotropic influences of genes on both traits.

Recent advances in understanding the genetics of sleep

Sleep is a ubiquitous and complex behavior in both its manifestation and regulation. Despite its essential role in maintaining optimal performance, health, and well-being, the genetic mechanisms underlying sleep remain poorly understood. Here, we review the forward genetic approaches undertaken in the last four years to elucidate the genes and gene pathways affecting sleep and its regulation. Despite an increasing number of studies and mining large databases, a coherent picture on “sleep” genes has yet to emerge. We highlight the results achieved by using unbiased genetic screens mainly in humans, mice, and fruit flies with an emphasis on normal sleep and make reference to lessons learned from the circadian field.

Shared Causal Paths underlying Alzheimer's dementia and Type 2 Diabetes

Although Alzheimer's disease (AD) is a central nervous system disease and type 2 diabetes MELLITUS (T2DM) is a metabolic disorder, an increasing number of genetic epidemiological studies show clear link between AD and T2DM. The current approach to uncovering the shared pathways between AD and T2DM involves association analysis; however such analyses lack power to discover the mechanisms of the diseases. As an alternative, we developed novel causal inference methods for genetic studies of AD and T2DM and pipelines for systematic multi-omic casual analysis to infer multilevel omics causal networks for the discovery of common paths from genetic variants to AD and T2DM. The proposed pipelines were applied to 448 individuals from the ROSMAP Project. We identified 13 shared causal genes, 16 shared causal pathways between AD and T2DM, and 754 gene expression and 101 gene methylation nodes that were connected to both AD and T2DM in multi-omics causal networks.

Onset of regular cannabis use and adult sleep duration: Genetic variation and the implications of a predictive relationship

BACKGROUND

Limited evidence suggests that early cannabis use is associated with sleep problems. Research is needed to understand the developmental impact of early regular cannabis use on later adult sleep duration.

METHODS

In a sample of 1656 adult twins (56% female, Mean age = 25.79yrs), linear mixed effects models were used to analyze the influence of retrospectively assessed age of onset of regular cannabis use on adult sleep duration controlling for sex, depression, and current substance use. Twin analyses provided genetic and environmental variance estimates as well as insights into the association and potential casual relationships between these traits.

RESULTS

Earlier age of onset for regular cannabis use was significantly associated with shorter adult sleep duration on both weekdays (β = -0.13, 95% CI = [-0.23, -0.04]) and weekends (β = -0.18, 95% CI = [-0.27, -0.08]). Additive genetics significantly contributed to the onset of regular cannabis use (a² = 76%, 95% CI = [68, 85]) and adult weekend sleep duration (a² = 20%, 95% CI = [11, 32]). We found evidence of a significant genetic correlation (r_A = -0.31, 95% CI = [-0.41, -0.15]) between these two traits and our best fitting model was consistent with early onset of regular cannabis use causing shorter adult weekend sleep duration (β = -0.11, 95% CI = [-0.18, -0.03]).

CONCLUSIONS

Our results are consistent with the hypothesis that early onset of regular cannabis use may have a negative impact on adult sleep duration.

Parkinson's Disease is Associated with Dysregulations of a Dopamine-Modulated Gene Network Relevant to Sleep and Affective Neurobehaviors in the Striatum

In addition to the characteristic motor symptoms, Parkinson’s disease (PD) often involves a constellation of sleep and mood symptoms. However, the mechanisms underlying these comorbidities are largely unknown. We have previously reconstructed gene networks in the striatum of a population of (C57BL/6J x A/J) F2 mice and associated the networks to sleep and affective phenotypes, providing a resource for integrated analyses to investigate perturbed sleep and affective functions at the gene network level. Combining this resource with PD-relevant transcriptomic datasets from humans and mice, we identified four networks that showed elevated gene expression in PD patients, including a circadian clock and mitotic network that was altered similarly in mouse models of PD. We then utilized multiple types of omics data from public databases and linked this gene network to postsynaptic dopamine signaling in the striatum, CDK1-modulated transcriptional regulation, and the genetic susceptibility of PD. These findings suggest that dopamine deficiency, a key aspect of PD pathology, perturbs a circadian/mitotic gene network in striatal neurons. Since the normal functions of this network were relevant to sleep and affective behaviors, these findings implicate that dysregulation of functional gene networks may be involved in the emergence of non-motor symptoms in PD. Our analyses present a framework for integrating multi-omics data from diverse sources in mice and humans to reveal insights into comorbid symptoms of complex diseases.

REM sleep's unique associations with corticosterone regulation, apoptotic pathways, and behavior in chronic stress in mice

Sleep disturbances are common in stress-related disorders but the nature of these sleep disturbances and how they relate to changes in the stress hormone corticosterone and changes in gene expression remained unknown. Here we demonstrate that in response to chronic mild stress, rapid–eye-movement sleep (REMS), a sleep state involved in emotion regulation and fear conditioning, changed first and more so than any other measured sleep characteristic. Transcriptomic profiles related to REMS continuity and theta oscillations overlapped with those for corticosterone, as well as with predictors for anhedonia, and were enriched for apoptotic pathways. These data highlight the central role of REMS in response to stress and warrant further investigation into REMS’s involvement in stress-related mental health disorders.

One of sleep’s putative functions is mediation of adaptation to waking experiences. Chronic stress is a common waking experience; however, which specific aspect of sleep is most responsive, and how sleep changes relate to behavioral disturbances and molecular correlates remain unknown. We quantified sleep, physical, endocrine, and behavioral variables, as well as the brain and blood transcriptome in mice exposed to 9 weeks of unpredictable chronic mild stress (UCMS). Comparing 46 phenotypic variables revealed that rapid–eye-movement sleep (REMS), corticosterone regulation, and coat state were most responsive to UCMS. REMS theta oscillations were enhanced, whereas delta oscillations in non-REMS were unaffected. Transcripts affected by UCMS in the prefrontal cortex, hippocampus, hypothalamus, and blood were associated with inflammatory and immune responses. A machine-learning approach controlling for unspecific UCMS effects identified transcriptomic predictor sets for REMS parameters that were enriched in 193 pathways, including some involved in stem cells, immune response, and apoptosis and survival. Only three pathways were enriched in predictor sets for non-REMS. Transcriptomic predictor sets for variation in REMS continuity and theta activity shared many pathways with corticosterone regulation, in particular pathways implicated in apoptosis and survival, including mitochondrial apoptotic machinery. Predictor sets for REMS and anhedonia shared pathways involved in oxidative stress, cell proliferation, and apoptosis. These data identify REMS as a core and early element of the response to chronic stress, and identify apoptosis and survival pathways as a putative mechanism by which REMS may mediate the response to stressful waking experiences.

Variation among intact tissue samples reveals the core transcriptional features of human CNS cell classes

It is widely assumed that cells must be physically isolated to study their molecular profiles. However, intact tissue samples naturally exhibit variation in cellular composition, which drives covariation of cell-class-specific molecular features. By analyzing transcriptional covariation in 7,221 intact CNS samples from 840 neurotypical individuals, representing billions of cells, we reveal the core transcriptional identities of major CNS cell classes in humans. By modeling intact CNS transcriptomes as a function of variation in cellular composition, we identify cell-class-specific transcriptional differences in Alzheimer's disease, among brain regions, and between species. Among these, we show that PMP2 is expressed by human but not mouse astrocytes and significantly increases mouse astrocyte size upon ectopic expression in vivo, causing them to more closely resemble their human counterparts. Our work is available as an online resource ( http://oldhamlab.ctec.ucsf.edu/ ) and provides a generalizable strategy for determining the core molecular features of cellular identity in intact biological systems.

A systems genetics resource and analysis of sleep regulation in the mouse

Sleep is essential for optimal brain functioning and health, but the biological substrates through which sleep delivers these beneficial effects remain largely unknown. We used a systems genetics approach in the BXD genetic reference population (GRP) of mice and assembled a comprehensive experimental knowledge base comprising a deep "sleep-wake" phenome, central and peripheral transcriptomes, and plasma metabolome data, collected under undisturbed baseline conditions and after sleep deprivation (SD). We present analytical tools to interactively interrogate the database, visualize the molecular networks altered by sleep loss, and prioritize candidate genes. We found that a one-time, short disruption of sleep already extensively reshaped the systems genetics landscape by altering 60%-78% of the transcriptomes and the metabolome, with numerous genetic loci affecting the magnitude and direction of change. Systems genetics integrative analyses drawing on all levels of organization imply α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor trafficking and fatty acid turnover as substrates of the negative effects of insufficient sleep. Our analyses demonstrate that genetic heterogeneity and the effects of insufficient sleep itself on the transcriptome and metabolome are far more widespread than previously reported.

Cross-species systems analysis identifies gene networks differentially altered by sleep loss and depression

To understand the transcriptomic organization underlying sleep and affective function, we studied a population of (C57BL/6J × 129S1/SvImJ) F2 mice by measuring 283 affective and sleep phenotypes and profiling gene expression across four brain regions. We identified converging molecular bases for sleep and affective phenotypes at both the single-gene and gene-network levels. Using publicly available transcriptomic datasets collected from sleep-deprived mice and patients with major depressive disorder (MDD), we identified three cortical gene networks altered by the sleep/wake state and depression. The network-level actions of sleep loss and depression were opposite to each other, providing a mechanistic basis for the sleep disruptions commonly observed in depression, as well as the reported acute antidepressant effects of sleep deprivation. We highlight one particular network composed of circadian rhythm regulators and neuronal activity-dependent immediate-early genes. The key upstream driver of this network, Arc, may act as a nexus linking sleep and depression. Our data provide mechanistic insights into the role of sleep in affective function and MDD.

Quantitative phosphoproteomic analysis of the molecular substrates of sleep need

Sleep and wake have global effects on brain physiology, from molecular changes^1-4 and neuronal activities to synaptic plasticity^3-7. Sleep-wake homeostasis is maintained by the generation of a sleep need that accumulates during waking and dissipates during sleep^8-11. Here we investigate the molecular basis of sleep need using quantitative phosphoproteomic analysis of the sleep-deprived and Sleepy mouse models of increased sleep need. Sleep deprivation induces cumulative phosphorylation of the brain proteome, which dissipates during sleep. Sleepy mice, owing to a gain-of-function mutation in the Sik3 gene ¹² , have a constitutively high sleep need despite increased sleep amount. The brain proteome of these mice exhibits hyperphosphorylation, similar to that seen in the brain of sleep-deprived mice. Comparison of the two models identifies 80 mostly synaptic sleep-need-index phosphoproteins (SNIPPs), in which phosphorylation states closely parallel changes of sleep need. SLEEPY, the mutant SIK3 protein, preferentially associates with and phosphorylates SNIPPs. Inhibition of SIK3 activity reduces phosphorylation of SNIPPs and slow wave activity during non-rapid-eye-movement sleep, the best known measurable index of sleep need, in both Sleepy mice and sleep-deprived wild-type mice. Our results suggest that phosphorylation of SNIPPs accumulates and dissipates in relation to sleep need, and therefore SNIPP phosphorylation is a molecular signature of sleep need. Whereas waking encodes memories by potentiating synapses, sleep consolidates memories and restores synaptic homeostasis by globally downscaling excitatory synapses^4-6. Thus, the phosphorylation-dephosphorylation cycle of SNIPPs may represent a major regulatory mechanism that underlies both synaptic homeostasis and sleep-wake homeostasis.

The MICALs are a Family of F-actin Dismantling Oxidoreductases Conserved from Drosophila to Humans

Cellular form and function – and thus normal development and physiology – are specified via proteins that control the organization and dynamic properties of the actin cytoskeleton. Using the Drosophila model, we have recently identified an unusual actin regulatory enzyme, Mical, which is directly activated by F-actin to selectively post-translationally oxidize and destabilize filaments – regulating numerous cellular behaviors. Mical proteins are also present in mammals, but their actin regulatory properties, including comparisons among different family members, remain poorly defined. We now find that each human MICAL family member, MICAL-1, MICAL-2, and MICAL-3, directly induces F-actin dismantling and controls F-actin-mediated cellular remodeling. Specifically, each human MICAL selectively associates with F-actin, which directly induces MICALs catalytic activity. We also find that each human MICAL uses an NADPH-dependent Redox activity to post-translationally oxidize actin’s methionine (M) M44/M47 residues, directly dismantling filaments and limiting new polymerization. Genetic experiments also demonstrate that each human MICAL drives F-actin disassembly in vivo, reshaping cells and their membranous extensions. Our results go on to reveal that MsrB/SelR reductase enzymes counteract each MICAL’s effect on F-actin in vitro and in vivo. Collectively, our results therefore define the MICALs as an important phylogenetically-conserved family of catalytically-acting F-actin disassembly factors.

Sex-specific transcriptional signatures in human depression

Major depressive disorder (MDD) is a leading cause of disease burden worldwide. While the incidence, symptoms and treatment of MDD all point toward major sex differences, the molecular mechanisms underlying this sexual dimorphism remain largely unknown. Here, combining differential expression and gene coexpression network analyses, we provide a comprehensive characterization of male and female transcriptional profiles associated with MDD across six brain regions. We overlap our human profiles with those from a mouse model, chronic variable stress, and capitalize on converging pathways to define molecular and physiological mechanisms underlying the expression of stress susceptibility in males and females. Our results show a major rearrangement of transcriptional patterns in MDD, with limited overlap between males and females, an effect seen in both depressed humans and stressed mice. We identify key regulators of sex-specific gene networks underlying MDD and confirm their sex-specific impact as mediators of stress susceptibility. For example, downregulation of the female-specific hub gene Dusp6 in mouse prefrontal cortex mimicked stress susceptibility in females, but not males, by increasing ERK signaling and pyramidal neuron excitability. Such Dusp6 downregulation also recapitulated the transcriptional remodeling that occurs in prefrontal cortex of depressed females. Together our findings reveal marked sexual dimorphism at the transcriptional level in MDD and highlight the importance of studying sex-specific treatments for this disorder.

Dissecting the Molecular Mechanisms of Neurodegenerative Diseases through Network Biology

Neurodegenerative diseases are rarely caused by a mutation in a single gene but rather influenced by a combination of genetic, epigenetic and environmental factors. Emerging high-throughput technologies such as RNA sequencing have been instrumental in deciphering the molecular landscape of neurodegenerative diseases, however, the interpretation of such large amounts of data remains a challenge. Network biology has become a powerful platform to integrate multiple omics data to comprehensively explore the molecular networks in the context of health and disease. In this review article, we highlight recent advances in network biology approaches with an emphasis in brain-networks that have provided insights into the molecular mechanisms leading to the most prevalent neurodegenerative diseases including Alzheimer’s (AD), Parkinson’s (PD) and Huntington’s diseases (HD). We discuss how integrative approaches using multi-omics data from different tissues have been valuable for identifying biomarkers and therapeutic targets. In addition, we discuss the challenges the field of network medicine faces toward the translation of network-based findings into clinically actionable tools for personalized medicine applications.

"Gestaltomics": Systems Biology Schemes for the Study of Neuropsychiatric Diseases

The integration of different sources of biological information about what defines a behavioral phenotype is difficult to unify in an entity that reflects the arithmetic sum of its individual parts. In this sense, the challenge of Systems Biology for understanding the “psychiatric phenotype” is to provide an improved vision of the shape of the phenotype as it is visualized by “Gestalt” psychology, whose fundamental axiom is that the observed phenotype (behavior or mental disorder) will be the result of the integrative composition of every part. Therefore, we propose the term “Gestaltomics” as a term from Systems Biology to integrate data coming from different sources of information (such as the genome, transcriptome, proteome, epigenome, metabolome, phenome, and microbiome). In addition to this biological complexity, the mind is integrated through multiple brain functions that receive and process complex information through channels and perception networks (i.e., sight, ear, smell, memory, and attention) that in turn are programmed by genes and influenced by environmental processes (epigenetic). Today, the approach of medical research in human diseases is to isolate one disease for study; however, the presence of an additional disease (co-morbidity) or more than one disease (multimorbidity) adds complexity to the study of these conditions. This review will present the challenge of integrating psychiatric disorders at different levels of information (Gestaltomics). The implications of increasing the level of complexity, for example, studying the co-morbidity with another disease such as cancer, will also be discussed.

Systems Genetic Analyses Highlight a TGFβ-FOXO3 Dependent Striatal Astrocyte Network Conserved across Species and Associated with Stress, Sleep, and Huntington's Disease

Recent systems-based analyses have demonstrated that sleep and stress traits emerge from shared genetic and transcriptional networks, and clinical work has elucidated the emergence of sleep dysfunction and stress susceptibility as early symptoms of Huntington's disease. Understanding the biological bases of these early non-motor symptoms may reveal therapeutic targets that prevent disease onset or slow disease progression, but the molecular mechanisms underlying this complex clinical presentation remain largely unknown. In the present work, we specifically examine the relationship between these psychiatric traits and Huntington's disease (HD) by identifying striatal transcriptional networks shared by HD, stress, and sleep phenotypes. First, we utilize a systems-based approach to examine a large publicly available human transcriptomic dataset for HD (GSE3790 from GEO) in a novel way. We use weighted gene coexpression network analysis and differential connectivity analyses to identify transcriptional networks dysregulated in HD, and we use an unbiased ranking scheme that leverages both gene- and network-level information to identify a novel astrocyte-specific network as most relevant to HD caudate. We validate this result in an independent HD cohort. Next, we computationally predict FOXO3 as a regulator of this network, and use multiple publicly available in vitro and in vivo experimental datasets to validate that this astrocyte HD network is downstream of a signaling pathway important in adult neurogenesis (TGFβ-FOXO3). We also map this HD-relevant caudate subnetwork to striatal transcriptional networks in a large (n = 100) chronically stressed (B6xA/J)F2 mouse population that has been extensively phenotyped (328 stress- and sleep-related measurements), and we show that this striatal astrocyte network is correlated to sleep and stress traits, many of which are known to be altered in HD cohorts. We identify causal regulators of this network through Bayesian network analysis, and we highlight their relevance to motor, mood, and sleep traits through multiple in silico approaches, including an examination of their protein binding partners. Finally, we show that these causal regulators may be therapeutically viable for HD because their downstream network was partially modulated by deep brain stimulation of the subthalamic nucleus, a medical intervention thought to confer some therapeutic benefit to HD patients. In conclusion, we show that an astrocyte transcriptional network is primarily associated to HD in the caudate and provide evidence for its relationship to molecular mechanisms of neural stem cell homeostasis. Furthermore, we present a unified systems-based framework for identifying gene networks that are associated with complex non-motor traits that manifest in the earliest phases of HD. By analyzing and integrating multiple independent datasets, we identify a point of molecular convergence between sleep, stress, and HD that reflects their phenotypic comorbidity and reveals a molecular pathway involved in HD progression.

Sleep as spatiotemporal integration of biological processes that evolved to periodically reinforce neurodynamic and metabolic homeostasis: The 2m3d paradigm of sleep

Sleep continues to perplex scientists and researchers. Despite decades of sleep research, we still lack a clear understanding of the biological functions and evolution of sleep. In this review, we will examine sleep from a functional and phylogenetic perspective and describe some important conceptual gaps in understanding sleep. Classical theories of the biology and evolution of sleep emphasize sensory activation, energy balance, and metabolic homeostasis. Advances in electrophysiology, functional neuroimaging, and neuroplasticity allow us to view sleep within the framework of neural dynamics. With this paradigm shift, we have come to realize the importance of neurodynamic homeostasis in shaping the biology of sleep. Evidently, animals sleep to achieve neurodynamic and metabolic homeostasis. We are not aware of any framework for understanding sleep where neurodynamic, metabolic, homeostatic, chronophasic, and afferent variables are all taken into account. This motivated us to propose the two-mode three-drive (2m3d) paradigm of sleep. In the 2m3d paradigm, local neurodynamic/metabolic (N/M) processes switch between two modes-m0 and m1-in response to three drives-afferent, chronophasic, and homeostatic. The spatiotemporal integration of local m0/m1 operations gives rise to the global states of sleep and wakefulness. As a framework of evolution, the 2m3d paradigm allows us to view sleep as a robust adaptive strategy that evolved so animals can periodically reinforce neurodynamic and metabolic homeostasis while remaining sensitive to their internal and external environment.

Comparing Alzheimer's and Parkinson's diseases networks using graph communities structure

BACKGROUND

Recent advances in large datasets analysis offer new insights to modern biology allowing system-level investigation of pathologies. Here we describe a novel computational method that exploits the ever-growing amount of "omics" data to shed light on Alzheimer's and Parkinson's diseases. Neurological disorders exhibit a huge number of molecular alterations due to a complex interplay between genetic and environmental factors. Classical reductionist approaches are focused on a few elements, providing a narrow overview of the etiopathogenic complexity of multifactorial diseases. On the other hand, high-throughput technologies allow the evaluation of many components of biological systems and their behaviors. Analysis of Parkinson's Disease (PD) and Alzheimer's Disease (AD) from a network perspective can highlight proteins or pathways common but differently represented that can be discriminating between the two pathological conditions, thus highlight similarities and differences.

RESULTS

In this work we propose a strategy that exploits network community structure identified with a state-of-the-art network community discovery algorithm called InfoMap, which takes advantage of information theory principles. We used two similarity measurements to quantify functional and topological similarities between the two pathologies. We built a Similarity Matrix to highlight similar communities and we analyzed statistically significant GO terms found in clustered areas of the matrix and in network communities. Our strategy allowed us to identify common known and unknown processes including DNA repair, RNA metabolism and glucose metabolism not detected with simple GO enrichment analysis. In particular, we were able to capture the connection between mitochondrial dysfunction and metabolism (glucose and glutamate/glutamine).

CONCLUSIONS

This approach allows the identification of communities present in both pathologies which highlight common biological processes. Conversely, the identification of communities without any counterpart can be used to investigate processes that are characteristic of only one of the two pathologies. In general, the same strategy can be applied to compare any pair of biological networks.