Multi-omic Directed Networks Describe Features of Gene Regulation in Aged Brains and Expand the Set of Genes Driving Cognitive Decline

The promise of multi-omics approaches to discover biological alterations with clinical relevance in Alzheimer's disease

Beyond the core features of Alzheimer's disease (AD) pathology, i.e. amyloid pathology, tau-related neurodegeneration and microglia response, multiple other molecular alterations and pathway dysregulations have been observed in AD. Their inter-individual variations, complex interactions and relevance for clinical manifestation and disease progression remain poorly understood, however. Heterogeneity at both pathophysiological and clinical levels complicates diagnosis, prognosis, treatment and drug design and testing. High-throughput "omics" comprise unbiased and untargeted data-driven methods which allow the exploration of a wide spectrum of disease-related changes at different endophenotype levels without focussing a priori on specific molecular pathways or molecules. Crucially, new methodological and statistical advances now allow for the integrative analysis of data resulting from multiple and different omics methods. These multi-omics approaches offer the unique advantage of providing a more comprehensive characterisation of the AD endophenotype and to capture molecular signatures and interactions spanning various biological levels. These new insights can then help decipher disease mechanisms more deeply. In this review, we describe the different multi-omics tools and approaches currently available and how they have been applied in AD research so far. We discuss how multi-omics can be used to explore molecular alterations related to core features of the AD pathologies and how they interact with comorbid pathological alterations. We further discuss whether the identified pathophysiological changes are relevant for the clinical manifestation of AD, in terms of both cognitive impairment and neuropsychiatric symptoms, and for clinical disease progression over time. Finally, we address the opportunities for multi-omics approaches to help discover novel biomarkers for diagnosis and monitoring of relevant pathophysiological processes, along with personalised intervention strategies in AD.

Cell type-specific changes identified by single-cell transcriptomics in Alzheimer's disease

The rapid advancement of single-cell transcriptomics in neurology has allowed for profiling of post-mortem human brain tissue across multiple diseases. Over the past 3 years, several studies have examined tissue from donors with and without diagnoses of Alzheimer's disease, highlighting key changes in cell type composition and molecular signatures associated with pathology and, in some cases, cognitive decline. Although all of these studies have generated single-cell/nucleus RNA-seq or ATAC-seq data from the full array of major cell classes in the brain, they have each focused on changes in specific cell types. Here, we synthesize the main findings from these studies and contextualize them in the overall space of large-scale omics studies of Alzheimer's disease. Finally, we touch upon new horizons in the field, in particular advancements in high-resolution spatial interrogation of tissue and multi-modal efforts-and how they are likely to further advance mechanistic and target-selection studies on Alzheimer's disease.

Spatiotemporal specificity of correlated DNA methylation and gene expression pairs across different human tissues and stages of brain development

DNA methylation (DNAm) that occurs on promoter regions is primarily considered to repress gene expression. Previous studies indicated that DNAm could also show positive correlations with gene expression. Both DNAm and gene expression profiles are known to be tissue- and development-specific. This study aims to investigate how DNAm and gene expression are coordinated across different human tissues and developmental stages, as well as the biological significance of such correlations. By analyzing 2,239 samples with both DNAm and gene expression data in the same human subjects obtained from six published datasets, we evaluated the correlations between gene and CpG pairs (GCPs) at cis-regions and compared significantly correlated GCPs (cGCPs) across different tissues and brains at different age groups. A total of 37,363 cGCPs was identified in the six datasets; approximately 38% of the cGCPs were positively correlated. The majority (>90%) of cGCPs was tissue- or development-specific. We also observed that the correlation direction can be opposite in different tissues and ages. Further analysis highlights the importance of cGCPs for their cellular functions and potential roles in complex traits and human diseases. For instance, the early developmental brain possessed a highly unique set of cGCPs that were associated with neurogenesis and psychiatric disorders. By assessing the epigenetic factors involved in cGCPs, we discovered novel regulatory mechanisms of positive cGCPs distinct from negative cGCPs, which were related to multiple factors, such as H3K27me3, CTCF, and JARD2. The catalogue of cGCPs compiled can be used to guide functional interpretation of genetic and epigenetic studies.

Elevated ganglioside GM2 activator (GM2A) in human brain tissue reduces neurite integrity and spontaneous neuronal activity

BACKGROUND

Alzheimer's Disease (AD) affects millions globally, but therapy development is lagging. New experimental systems that monitor neuronal functions in conditions approximating the AD brain may be beneficial for identifying new therapeutic strategies.

METHODS

We expose cultured neurons to aqueous-soluble human brain extract from 43 individuals across a spectrum of AD pathology. Multi-electrode arrays (MEAs) and live-cell imaging were used to assess neuronal firing and neurite integrity (NI), respectively, following treatments of rat cortical neurons (MEA) and human iPSC-derived neurons (iN) with human brain extracts.

RESULTS

We observe associations between spontaneous activity and Aβ42:40 levels, between neurite integrity and oligomeric Aβ, and between neurite integrity and tau levels present in the brain extracts. However, these associations with Aβ and tau do not fully account for the effects observed. Proteomic profiling of the brain extracts revealed additional candidates correlated with neuronal structure and activity. Neurotoxicity in MEA and NI assays was associated with proteins implicated in lysosomal storage disorders, while neuroprotection was associated with proteins of the WAVE regulatory complex controlling actin cytoskeleton dynamics. Elevated ganglioside GM2 activator (GM2A) associates with reductions in both NI and MEA activity, and cell-derived GM2A alone is sufficient to induce a loss of neurite integrity and a reduction in neuronal firing.

CONCLUSIONS

The techniques and data herein introduce a system for modeling neuronal vulnerability in response to factors in the human brain and provide insights into proteins potentially contributing to AD pathogenesis.

Thirteen Independent Genetic Loci Associated with Preserved Processing Speed in a Study of Cognitive Resilience in 330,097 Individuals in the UK Biobank

Cognitive resilience is the ability to withstand the negative effects of stress on cognitive functioning and is important for maintaining quality of life while aging. The UK Biobank does not have measurements of the same cognitive phenotype at distal time points. Therefore, we used education years (EY) as a proxy phenotype for past cognitive performance and current cognitive performance was based on processing speed. This represented an average time span of 40 years between past and current cognitive performance in 330,097 individuals. A confounding factor was that EY is highly polygenic and masked the genetics of resilience. To overcome this, we employed Genomics Structural Equation Modelling (GenomicSEM) to perform a genome-wide association study (GWAS)-by-subtraction using two GWAS, one GWAS of EY and resilience and a second GWAS of EY but not resilience, to generate a GWAS of Resilience. Using independent discovery and replication samples, we found 13 independent genetic loci for Resilience. Functional analyses showed enrichment in several brain regions and specific cell types. Gene-set analyses implicated the biological process "neuron differentiation", the cellular component "synaptic part" and the "WNT signalosome". Mendelian randomisation analysis showed a causative effect of white matter volume on cognitive resilience. These results may contribute to the neurobiological understanding of resilience.

Stem cell-derived neurons reflect features of protein networks, neuropathology, and cognitive outcome of their aged human donors

We have generated a controlled and manipulable resource that captures genetic risk for Alzheimer's disease: iPSC lines from 53 individuals coupled with RNA and proteomic profiling of both iPSC-derived neurons and brain tissue of the same individuals. Data collected for each person include genome sequencing, longitudinal cognitive scores, and quantitative neuropathology. The utility of this resource is exemplified here by analyses of neurons derived from these lines, revealing significant associations between specific Aβ and tau species and the levels of plaque and tangle deposition in the brain and, more importantly, with the trajectory of cognitive decline. Proteins and networks are identified that are associated with AD phenotypes in iPSC neurons, and relevant associations are validated in brain. The data presented establish this iPSC collection as a resource for investigating person-specific processes in the brain that can aid in identifying and validating molecular pathways underlying AD.

A multi-omic study for uncovering molecular mechanisms associated with hyperammonemia-induced cerebellar function impairment in rats

Patients with liver cirrhosis may develop covert or minimal hepatic encephalopathy (MHE). Hyperammonemia (HA) and peripheral inflammation play synergistic roles in inducing the cognitive and motor alterations in MHE. The cerebellum is one of the main cerebral regions affected in MHE. Rats with chronic HA show some motor and cognitive alterations reproducing neurological impairment in cirrhotic patients with MHE. Neuroinflammation and altered neurotransmission and signal transduction in the cerebellum from hyperammonemic (HA) rats are associated with motor and cognitive dysfunction, but underlying mechanisms are not completely known. The aim of this work was to use a multi-omic approach to study molecular alterations in the cerebellum from hyperammonemic rats to uncover new molecular mechanisms associated with hyperammonemia-induced cerebellar function impairment. We analyzed metabolomic, transcriptomic, and proteomic data from the same cerebellums from control and HA rats and performed a multi-omic integrative analysis of signaling pathway enrichment with the PaintOmics tool. The histaminergic system, corticotropin-releasing hormone, cyclic GMP-protein kinase G pathway, and intercellular communication in the cerebellar immune system were some of the most relevant enriched pathways in HA rats. In summary, this is a good approach to find altered pathways, which helps to describe the molecular mechanisms involved in the alteration of brain function in rats with chronic HA and to propose possible therapeutic targets to improve MHE symptoms.

Identifying Mechanisms of Normal Cognitive Aging Using a Novel Mouse Genetic Reference Panel

Developing strategies to maintain cognitive health is critical to quality of life during aging. The basis of healthy cognitive aging is poorly understood; thus, it is difficult to predict who will have normal cognition later in life. Individuals may have higher baseline functioning (cognitive reserve) and others may maintain or even improve with age (cognitive resilience). Understanding the mechanisms underlying cognitive reserve and resilience may hold the key to new therapeutic strategies for maintaining cognitive health. However, reserve and resilience have been inconsistently defined in human studies. Additionally, our understanding of the molecular and cellular bases of these phenomena is poor, compounded by a lack of longitudinal molecular and cognitive data that fully capture the dynamic trajectories of cognitive aging. Here, we used a genetically diverse mouse population (B6-BXDs) to characterize individual differences in cognitive abilities in adulthood and investigate evidence of cognitive reserve and/or resilience in middle-aged mice. We tested cognitive function at two ages (6 months and 14 months) using y-maze and contextual fear conditioning. We observed heritable variation in performance on these traits (h 2 RIx̄ = 0.51-0.74), suggesting moderate to strong genetic control depending on the cognitive domain. Due to the polygenetic nature of cognitive function, we did not find QTLs significantly associated with y-maze, contextual fear acquisition (CFA) or memory, or decline in cognitive function at the genome-wide level. To more precisely interrogate the molecular regulation of variation in these traits, we employed RNA-seq and identified gene networks related to transcription/translation, cellular metabolism, and neuronal function that were associated with working memory, contextual fear memory, and cognitive decline. Using this method, we nominate the Trio gene as a modulator of working memory ability. Finally, we propose a conceptual framework for identifying strains exhibiting cognitive reserve and/or resilience to assess whether these traits can be observed in middle-aged B6-BXDs. Though we found that earlier cognitive reserve evident early in life protects against cognitive impairment later in life, cognitive performance and age-related decline fell along a continuum, with no clear genotypes emerging as exemplars of exceptional reserve or resilience - leading to recommendations for future use of aging mouse populations to understand the nature of cognitive reserve and resilience.

Considerations for integrative multi-omic approaches to explore Alzheimer's disease mechanisms

The past decade has seen the maturation of multiple different forms of high-dimensional molecular profiling to the point that these methods could be deployed in initially hundreds and more recently thousands of human samples. In the field of Alzheimer's disease (AD), these profiles have been applied to the target organ: the aging brain. In a growing number of cases, the same samples were profiled with multiple different approaches, yielding genetic, transcriptomic, epigenomic and proteomic data. Here, we review lessons learned so far as we move beyond quantitative trait locus (QTL) analyses which map the effect of genetic variation on molecular features to integrate multiple levels of "omic" data in an effort to identify the molecular drivers of AD. One thing is clear: no single layer of molecular or "omic" data is sufficient to capture the variance of AD or aging-related cognitive decline. Nonetheless, reproducible findings are emerging from current efforts, and there is evidence of convergence using different approaches. Thus, we are on the cusp of an acceleration of truly integrative studies as the availability of large numbers of well-characterized brain samples profiled in three or more dimensions enables the testing, comparison and refinement of analytic methods with which to dissect the molecular architecture of the aging brain.

Bayesian integrative analysis of epigenomic and transcriptomic data identifies Alzheimer's disease candidate genes and networks

Biomedical research studies have generated large multi-omic datasets to study complex diseases like Alzheimer’s disease (AD). An important aim of these studies is the identification of candidate genes that demonstrate congruent disease-related alterations across the different data types measured by the study. We developed a new method to detect such candidate genes in large multi-omic case-control studies that measure multiple data types in the same set of samples. The method is based on a gene-centric integrative coefficient quantifying to what degree consistent differences are observed in the different data types. For statistical inference, a Bayesian hierarchical model is used to study the distribution of the integrative coefficient. The model employs a conditional autoregressive prior to integrate a functional gene network and to share information between genes known to be functionally related. We applied the method to an AD dataset consisting of histone acetylation, DNA methylation, and RNA transcription data from human cortical tissue samples of 233 subjects, and we detected 816 genes with consistent differences between persons with AD and controls. The findings were validated in protein data and in RNA transcription data from two independent AD studies. Finally, we found three subnetworks of jointly dysregulated genes within the functional gene network which capture three distinct biological processes: myeloid cell differentiation, protein phosphorylation and synaptic signaling. Further investigation of the myeloid network indicated an upregulation of this network in early stages of AD prior to accumulation of hyperphosphorylated tau and suggested that increased CSF1 transcription in astrocytes may contribute to microglial activation in AD. Thus, we developed a method that integrates multiple data types and external knowledge of gene function to detect candidate genes, applied the method to an AD dataset, and identified several disease-related genes and processes demonstrating the usefulness of the integrative approach.

Recent technological advances have led to a new generation of studies that interrogate multiple molecular levels in the same target tissue of a set of subjects, generating complex multi-omic datasets with which to study disease mechanism. These datasets of genetic, epigenomic, transcriptomic, and other data have the potential to reveal novel biological insights; however, integrative analyses remain challenging and require new computational methods. We developed an integrative Bayesian approach to detect genes with consistent differences between case and control samples across multiple data types. The method further integrates prior knowledge about gene function in the form of a gene functional similarity network to improve statistical inference by sharing information between related genes. We applied our method to an Alzheimer’s disease dataset of epigenomic and transcriptomic data and detected and then validated several novel and known candidate genes as well as three major disease-related biological processes. One of these processes reflected microglial activation and included the cytokine CSF1. Single-nucleus data revealed that CSF1 was primarily upregulated in astrocytes, implicating the involvement of this cell type in microglial activation. Hence, we demonstrated that integrative analysis approaches to multi-omic datasets can improve candidate gene detection and thereby generate new insights into complex diseases.

Ageotypes: Distinct Biomolecular Trajectories in Human Aging

Many studies have demonstrated that biological age (BA) varies significantly among individuals of similar chronological age. A recent study by Ahadi et al. used longitudinal and deep multi-omic profiling to identify individuals with distinct BA phenotypes or 'ageotypes'. These ageotypes open new avenues to creating diagnostic and treatment strategies that may slow the aging process based on the unique biochemistry of each individual.

Quantitative Systems Pharmacology for Neuroscience Drug Discovery and Development: Current Status, Opportunities, and Challenges

The substantial progress made in the basic sciences of the brain has yet to be adequately translated to successful clinical therapeutics to treat central nervous system (CNS) diseases. Possible explanations include the lack of quantitative and validated biomarkers, the subjective nature of many clinical endpoints, and complex pharmacokinetic/pharmacodynamic relationships, but also the possibility that highly selective drugs in the CNS do not reflect the complex interactions of different brain circuits. Although computational systems pharmacology modeling designed to capture essential components of complex biological systems has been increasingly accepted in pharmaceutical research and development for oncology, inflammation, and metabolic disorders, the uptake in the CNS field has been very modest. In this article, a cross-disciplinary group with representatives from academia, pharma, regulatory, and funding agencies make the case that the identification and exploitation of CNS therapeutic targets for drug discovery and development can benefit greatly from a system and network approach that can span the gap between molecular pathways and the neuronal circuits that ultimately regulate brain activity and behavior. The National Institute of Neurological Disorders and Stroke (NINDS), in collaboration with the National Institute on Aging (NIA), National Institute of Mental Health (NIMH), National Institute on Drug Abuse (NIDA), and National Center for Advancing Translational Sciences (NCATS), convened a workshop to explore and evaluate the potential of a quantitative systems pharmacology (QSP) approach to CNS drug discovery and development. The objective of the workshop was to identify the challenges and opportunities of QSP as an approach to accelerate drug discovery and development in the field of CNS disorders. In particular, the workshop examined the potential for computational neuroscience to perform QSP-based interrogation of the mechanism of action for CNS diseases, along with a more accurate and comprehensive method for evaluating drug effects and optimizing the design of clinical trials. Following up on an earlier white paper on the use of QSP in general disease mechanism of action and drug discovery, this report focuses on new applications, opportunities, and the accompanying limitations of QSP as an approach to drug development in the CNS therapeutic area based on the discussions in the workshop with various stakeholders.

The Many Faces of Gene Regulation in Cancer: A Computational Oncogenomics Outlook

Cancer is a complex disease at many different levels. The molecular phenomenology of cancer is also quite rich. The mutational and genomic origins of cancer and their downstream effects on processes such as the reprogramming of the gene regulatory control and the molecular pathways depending on such control have been recognized as central to the characterization of the disease. More important though is the understanding of their causes, prognosis, and therapeutics. There is a multitude of factors associated with anomalous control of gene expression in cancer. Many of these factors are now amenable to be studied comprehensively by means of experiments based on diverse omic technologies. However, characterizing each dimension of the phenomenon individually has proven to fall short in presenting a clear picture of expression regulation as a whole. In this review article, we discuss some of the more relevant factors affecting gene expression control both, under normal conditions and in tumor settings. We describe the different omic approaches that we can use as well as the computational genomic analysis needed to track down these factors. Then we present theoretical and computational frameworks developed to integrate the amount of diverse information provided by such single-omic analyses. We contextualize this within a systems biology-based multi-omic regulation setting, aimed at better understanding the complex interplay of gene expression deregulation in cancer.

Genetic risk for Alzheimer's dementia predicts motor deficits through multi-omic systems in older adults

Alzheimer's disease manifests with both cognitive and motor deficits. However, the degree to which genetic risk of Alzheimer's dementia contributes to late-life motor impairment, and the specific molecular systems underlying these associations, are uncertain. Here, we adopted an integrative multi-omic approach to assess genetic influence on motor impairment in older adults and identified key molecular pathways that may mediate this risk. We built a polygenic risk score for clinical diagnosis of Alzheimer's dementia (AD-PRS) and examined its relationship to several motor phenotypes in 1885 older individuals from two longitudinal aging cohorts. We found that AD-PRS was associated with a previously validated composite motor scores and their components. The major genetic risk factor for sporadic Alzheimer's dementia, the APOE/TOMM40 locus, was not a major driver of these associations. To identify specific molecular features that potentially medicate the genetic risk into motor dysfunction, we examined brain multi-omics, including transcriptome, DNA methylation, histone acetylation (H3K9AC), and targeted proteomics, as well as diverse neuropathologies. We found that a small number of factors account for the majority of the influence of AD-PRS on motor function, which comprises paired helical filament tau-tangle density, H3K9AC in specific chromosomal regions encoding genes involved in neuromuscular process. These multi-omic factors have the potential to elucidate key molecular mechanisms developing motor impairment in the context of Alzheimer's dementia.

Molecular Signatures of the Aging Brain: Finding the Links Between Genes and Phenotypes

Aging is associated with cognitive decline and increased vulnerability to neurodegenerative diseases. The progressive extension of the average human lifespan is bound to lead to a corresponding increase in the fraction of cognitively impaired elderly individuals among the human population, with an enormous societal and economic burden. At the cellular and tissue levels, cognitive decline is linked to a reduction in specific neuronal subpopulations, a widespread decrease in synaptic plasticity and an increase in neuroinflammation due to an enhanced activation of astrocytes and microglia, but the molecular mechanisms underlying these functional changes during normal aging and in neuropathological conditions remain poorly understood. In this review, we summarize very recent and outstanding progress in elucidating the molecular changes associated with cognitive decline through the genome-wide profiling of aging brain cells at different molecular levels (genomic, epigenomic, transcriptomic, proteomic). We discuss how the correlation of different molecular and phenotypic traits driven by mathematical and computational analyses of large datasets has led to the prediction of key molecular nodes of neurodegenerative pathways, and provide a few examples of candidate regulators of cognitive decline identified with these approaches. Furthermore, we highlight the dysregulation of the synaptic transcriptome in neuronal cells and of the inflammatory transcriptome in glial cells as some of the key events during normal and neuropathological human brain aging.

Aging of the frontal lobe

Healthy aging is associated with numerous deficits in cognitive function, which have been attributed to changes within the prefrontal cortex (PFC). This chapter summarizes some of the most prominent cognitive changes associated with age-related alterations in the anatomy and physiology of the PFC. Specifically, aging of the PFC results in deficient aspects of cognitive control, including sustained attention, selective attention, inhibitory control, working memory, and multitasking abilities. Yet, not all cognitive functions associated with the PFC exhibit age-related declines, such as arithmetic, comprehension, emotion perception, and emotional control. Moreover, not all older adults exhibit declines in cognition. Multiple life-course and lifestyle factors, as well as genetics, play a role in the trajectory of cognitive performance across the life span. Thus many adults retain cognitive function well into advanced age. Moreover, the brain remains plastic throughout life and there is increasing evidence that most age-related declines in cognition can be remediated by various methods such as physical exercise, cognitive training, or noninvasive brain stimulation. Overall, because cognitive aging is associated with numerous life-course and lifestyle factors, successful aging likely begins in early life, while maintaining cognition or remediating declines is a life-long process.

The Molecular and Neuropathological Consequences of Genetic Risk for Alzheimer's Dementia

Alzheimer's dementia commonly impacts the health of older adults and lacks any preventative therapy. While Alzheimer's dementia risk has a substantial genetic component, the specific molecular mechanisms and neuropathologies triggered by most of the known genetic variants are unclear. Resultantly, they have shown limited influence on drug development portfolios to date. To facilitate our understanding of the consequences of Alzheimer's dementia susceptibility variants, we examined their relationship to a wide range of clinical, molecular and neuropathological features. Because the effect size of individual variants is typically small, we utilized a polygenic (overall) risk approach to identify the global impact of Alzheimer's dementia susceptibility variants. Under this approach, each individual has a polygenic risk score (PRS) that we related to clinical, molecular and neuropathological phenotypes. Applying this approach to 1,272 individuals who came to autopsy from one of two longitudinal aging cohorts, we observed that an individual's PRS was associated with cognitive decline and brain pathologies including beta-amyloid, tau-tangles, hippocampal sclerosis, and TDP-43, MIR132, four proteins including VGF, IGFBP5, and STX1A, and many chromosomal regions decorated with acetylation on histone H3 lysine 9 (H3K9Ac). While excluding the APOE/TOMM40 region (containing the single largest genetic risk factor for late-onset Alzheimer's dementia) in the calculation of the PRS resulted in a slightly weaker association with the molecular signatures, results remained significant. These PRS-associated brain pathologies and molecular signatures appear to mediate genetic risk, as they attenuated the association of the PRS with cognitive decline. Notably, the PRS induced changes in H3K9Ac throughout the genome, implicating it in large-scale chromatin changes. Thus, the PRS for Alzheimer's dementia (AD-PRS) showed effects on diverse clinical, molecular, and pathological systems, ranging from the epigenome to specific proteins. These convergent targets of a large number of genetic risk factors for Alzheimer's dementia will help define the experimental systems and models needed to test therapeutic targets, which are expected to be broadly effective in the aging population that carries diverse genetic risks for Alzheimer's dementia.