Identifying robust communities and multi-community nodes by combining top-down and bottom-up approaches to clustering

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

Dissecting the human leptomeninges at single-cell resolution

Emerging evidence shows that the meninges conduct essential immune surveillance and immune defense at the brain border, and the dysfunction of meningeal immunity contributes to aging and neurodegeneration. However, no study exists on the molecular properties of cell types within human leptomeninges. Here, we provide single nuclei profiling of dissected postmortem leptomeninges from aged individuals. We detect diverse cell types, including unique meningeal endothelial, mural, and fibroblast subtypes. For immune cells, we show that most T cells express CD8 and bear characteristics of tissue-resident memory T cells. We also identify distinct subtypes of border-associated macrophages (BAMs) that display differential gene expressions from microglia and express risk genes for Alzheimer's Disease (AD), as nominated by genome-wide association studies (GWAS). We discover cell-type-specific differentially expressed genes in individuals with Alzheimer's dementia, particularly in fibroblasts and BAMs. Indeed, when cultured, leptomeningeal cells display the signature of ex vivo AD fibroblasts upon amyloid-β treatment. We further explore ligand-receptor interactions within the leptomeningeal niche and computationally infer intercellular communications in AD. Thus, our study establishes a molecular map of human leptomeningeal cell types, providing significant insight into the border immune and fibrotic responses in AD.

Robust, scalable, and informative clustering for diverse biological networks

Clustering molecular data into informative groups is a primary step in extracting robust conclusions from big data. However, due to foundational issues in how they are defined and detected, such clusters are not always reliable, leading to unstable conclusions. We compare popular clustering algorithms across thousands of synthetic and real biological datasets, including a new consensus clustering algorithm-SpeakEasy2: Champagne. These tests identify trends in performance, show no single method is universally optimal, and allow us to examine factors behind variation in performance. Multiple metrics indicate SpeakEasy2 generally provides robust, scalable, and informative clusters for a range of applications.

A Molecular Basis of Human Brain Connectivity

Neuroimaging is commonly used to infer human brain connectivity, but those measurements are far-removed from the molecular underpinnings at synapses. To uncover the molecular basis of human brain connectivity, we analyzed a unique cohort of 98 individuals who provided neuroimaging and genetic data contemporaneous with dendritic spine morphometric, proteomic, and gene expression data from the superior frontal and inferior temporal gyri. Through cellular contextualization of the molecular data with dendritic spine morphology, we identified hundreds of proteins related to synapses, energy metabolism, and RNA processing that explain between-individual differences in functional connectivity and structural covariation. By integrating data at the genetic, molecular, subcellular, and tissue levels, we bridged the divergent fields of molecular biology and neuroimaging to identify a molecular basis of brain connectivity.

One-Sentence Summary

Dendritic spine morphometry and synaptic proteins unite the divergent fields of molecular biology and neuroimaging.

Inferring protein expression changes from mRNA in Alzheimer's dementia using deep neural networks

Identifying the molecular systems and proteins that modify the progression of Alzheimer's disease and related dementias (ADRD) is central to drug target selection. However, discordance between mRNA and protein abundance, and the scarcity of proteomic data, has limited our ability to advance candidate targets that are mainly based on gene expression. Therefore, by using a deep neural network that predicts protein abundance from mRNA expression, here we attempt to track the early protein drivers of ADRD. Specifically, by applying the clei2block deep learning model to 1192 brain RNA-seq samples, we identify protein modules and disease-associated expression changes that were not directly observed at the mRNA level. Moreover, pseudo-temporal trajectory inference based on the predicted proteome became more closely correlated with cognitive decline and hippocampal atrophy compared to RNA-based trajectories. This suggests that the predicted changes in protein expression could provide a better molecular representation of ADRD progression. Furthermore, overlaying clinical traits on protein pseudotime trajectory identifies protein modules altered before cognitive impairment. These results demonstrate how our method can be used to identify potential early protein drivers and possible drug targets for treating and/or preventing ADRD.

A network generator for covert network structures

We focus on organizational structures in covert networks, such as criminal or terrorist networks. Their members engage in illegal activities and attempt to hide their association and interactions with these networks. Hence, data about such networks are incomplete. We introduce a novel method of rewiring covert networks parameterized by the edge connectivity standard deviation. The generated networks are statistically similar to themselves and to the original network. The higher-level organizational structures are modeled as a multi-layer network while the lowest level uses the Stochastic Block Model. Such synthetic networks provide alternative structures for data about the original network. Using them, analysts can find structures that are frequent, therefore stable under perturbations. Another application is to anonymize generated networks and use them for testing new software developed in open research facilities. The results indicate that modeling edge structure and the hierarchy together is essential for generating networks that are statistically similar but not identical to each other or the original network. In experiments, we generate many synthetic networks from two covert networks. Only a few structures of synthetics networks repeat, with the most stable ones shared by 18% of all synthetic networks making them strong candidates for the ground truth structure.

Systems biology approaches to unravel the molecular and genetic architecture of Alzheimer's disease and related tauopathies

Over the years, genetic studies have identified multiple genetic risk variants associated with neurodegenerative disorders and helped reveal new biological pathways and genes of interest. However, genetic risk variants commonly reside in non-coding regions and may regulate distant genes rather than the nearest gene, as well as a gene's interaction partners in biological networks. Systems biology and functional genomics approaches provide the framework to unravel the functional significance of genetic risk variants in disease. In this review, we summarize the genetic and transcriptomic studies of Alzheimer's disease and related tauopathies and focus on the advantages of performing systems-level analyses to interrogate the biological pathways underlying neurodegeneration. Finally, we highlight new avenues of multi-omics analysis with single-cell approaches, which provide unparalleled opportunities to systematically explore cellular heterogeneity, and present an example of how to integrate publicly available single-cell datasets. Systems-level analysis has illuminated the function of many disease risk genes, but much work remains to study tauopathies and to understand spatiotemporal gene expression changes of specific cell types.

Do not waste your electrodes-principles of optimal electrode geometry for spike sorting

Objective. This study examines how the geometrical arrangement of electrodes influences spike sorting efficiency, and attempts to formalise principles for the design of electrode systems enabling optimal spike sorting performance.Approach. The clustering performance of KlustaKwik, a popular toolbox, was evaluated using semi-artificial multi-channel data, generated from a library of real spike waveforms recorded in the CA1 region of mouse Hippocampusin vivo.Main results. Based on spike sorting results under various channel configurations and signal levels, a simple model was established to describe the efficiency of different electrode geometries. Model parameters can be inferred from existing spike waveform recordings, which allowed quantifying both the cooperative effect between channels and the noise dependence of clustering performance.Significance. Based on the model, analytical and numerical results can be derived for the optimal spacing and arrangement of electrodes for one- and two-dimensional electrode systems, targeting specific brain areas.

Community Detection in Large-Scale Bipartite Biological Networks

Networks are useful tools to represent and analyze interactions on a large, or genome-wide scale and have therefore been widely used in biology. Many biological networks-such as those that represent regulatory interactions, drug-gene, or gene-disease associations-are of a bipartite nature, meaning they consist of two different types of nodes, with connections only forming between the different node sets. Analysis of such networks requires methodologies that are specifically designed to handle their bipartite nature. Community structure detection is a method used to identify clusters of nodes in a network. This approach is especially helpful in large-scale biological network analysis, as it can find structure in networks that often resemble a "hairball" of interactions in visualizations. Often, the communities identified in biological networks are enriched for specific biological processes and thus allow one to assign drugs, regulatory molecules, or diseases to such processes. In addition, comparison of community structures between different biological conditions can help to identify how network rewiring may lead to tissue development or disease, for example. In this mini review, we give a theoretical basis of different methods that can be applied to detect communities in bipartite biological networks. We introduce and discuss different scores that can be used to assess the quality of these community structures. We then apply a wide range of methods to a drug-gene interaction network to highlight the strengths and weaknesses of these methods in their application to large-scale, bipartite biological networks.

Neuroticism alters the transcriptome of the frontal cortex to contribute to the cognitive decline and onset of Alzheimer's disease

Accumulating evidence has suggested that the molecular transcriptional mechanism contributes to Alzheimer's disease (AD) and its endophenotypes of cognitive decline and neuropathological traits, β-amyloid (Aβ) and phosphorylated tangles (TAU). However, it is unknown what is the impact of the AD risk factors, personality characteristics assessed by the NEO Five-Factor Inventory, on the human brain's transcriptome. Using postmortem human brain samples from 466 subjects, we found that neuroticism has a significant overall impact on the brain transcriptome (omnibus P = 0.005) but not the other four personality characteristics. Focused on those cognitive decline related gene co-expressed modules, neuroticism has nominally significant associations (P < 0.05) with four neuronal modules, which are more related to PHFtau than Aβ across all eight brain regions. Furthermore, the effect of neuroticism on cognitive decline and AD might be mediated through the expression of module 7 and TAU pathology (P = 0.008). To conclude, neuroticism has a broad impact on the transcriptome of human brains, and its effect on cognitive decline and AD may be mediated through gene transcription programs related to TAU pathology.

A cortical immune network map identifies distinct microglial transcriptional programs associated with β-amyloid and Tau pathologies

Microglial dysfunction has been proposed as one of the many cellular mechanisms that can contribute to the development of Alzheimer's disease (AD). Here, using a transcriptional network map of the human frontal cortex, we identify five modules of co-expressed genes related to microglia and assess their role in the neuropathologic features of AD in 540 subjects from two cohort studies of brain aging. Two of these transcriptional programs-modules 113 and 114-relate to the accumulation of β-amyloid, while module 5 relates to tau pathology. We replicate these associations in brain epigenomic data and in two independent datasets. In terms of tau, we propose that module 5, a marker of activated microglia, may lead to tau accumulation and subsequent cognitive decline. We validate our model further by showing that three representative module 5 genes (ACADVL, TRABD, and VASP) encode proteins that are upregulated in activated microglia in AD.

Overlapping Community Detection Based on Membership Degree Propagation

A community in a complex network refers to a group of nodes that are densely connected internally but with only sparse connections to the outside. Overlapping community structures are ubiquitous in real-world networks, where each node belongs to at least one community. Therefore, overlapping community detection is an important topic in complex network research. This paper proposes an overlapping community detection algorithm based on membership degree propagation that is driven by both global and local information of the node community. In the method, we introduce a concept of membership degree, which not only stores the label information, but also the degrees of the node belonging to the labels. Then the conventional label propagation process could be extended to membership degree propagation, with the results mapped directly to the overlapping community division. Therefore, it obtains the partition result and overlapping node identification simultaneously and greatly reduces the computational time. The proposed algorithm was applied to a synthetic Lancichinetti–Fortunato–Radicchi (LFR) dataset and nine real-world datasets and compared with other up-to-date algorithms. The experimental results show that our proposed algorithm is effective and outperforms the comparison methods on most datasets. Our proposed method significantly improved the accuracy and speed of the overlapping node prediction. It can also substantially alleviate the computational complexity of community structure detection in general.

The genetics of circulating BDNF: towards understanding the role of BDNF in brain structure and function in middle and old ages

Brain-derived neurotrophic factor (BDNF) plays an important role in brain development and function. Substantial amounts of BDNF are present in peripheral blood, and may serve as biomarkers for Alzheimer’s disease incidence as well as targets for intervention to reduce Alzheimer’s disease risk. With the exception of the genetic polymorphism in the BDNF gene, Val66Met, which has been extensively studied with regard to neurodegenerative diseases, the genetic variation that influences circulating BDNF levels is unknown. We aimed to explore the genetic determinants of circulating BDNF levels in order to clarify its mechanistic involvement in brain structure and function and Alzheimer’s disease pathophysiology in middle-aged and old adults. Thus, we conducted a meta-analysis of genome-wide association study of circulating BDNF in 11 785 middle- and old-aged individuals of European ancestry from the Age, Gene/Environment Susceptibility-Reykjavik Study (AGES), the Framingham Heart Study (FHS), the Rotterdam Study and the Study of Health in Pomerania (SHIP-Trend). Furthermore, we performed functional annotation analysis and related the genetic polymorphism influencing circulating BDNF to common Alzheimer’s disease pathologies from brain autopsies. Mendelian randomization was conducted to examine the possible causal role of circulating BDNF levels with various phenotypes including cognitive function, stroke, diabetes, cardiovascular disease, physical activity and diet patterns. Gene interaction networks analysis was also performed. The estimated heritability of BDNF levels was 30% (standard error = 0.0246, P-value = 4 × 10⁻⁴⁸). We identified seven novel independent loci mapped near the BDNF gene and in BRD3, CSRNP1, KDELC2, RUNX1 (two single-nucleotide polymorphisms) and BDNF-AS. The expression of BDNF was associated with neurofibrillary tangles in brain tissues from the Religious Orders Study and Rush Memory and Aging Project (ROSMAP). Seven additional genes (ACAT1, ATM, NPAT, WDR48, TTC21A, SCN114 and COX7B) were identified through expression and protein quantitative trait loci analyses. Mendelian randomization analyses indicated a potential causal role of BDNF in cardioembolism. Lastly, Ingenuity Pathway Analysis placed circulating BDNF levels in four major networks. Our study provides novel insights into genes and molecular pathways associated with circulating BDNF levels and highlights the possible involvement of plaque instability as an underlying mechanism linking BDNF with brain neurodegeneration. These findings provide a foundation for a better understanding of BDNF regulation and function in the context of brain aging and neurodegenerative pathophysiology.

A multiresolution framework to characterize single-cell state landscapes

Dissecting the cellular heterogeneity embedded in single-cell transcriptomic data is challenging. Although many methods and approaches exist, identifying cell states and their underlying topology is still a major challenge. Here, we introduce the concept of multiresolution cell-state decomposition as a practical approach to simultaneously capture both fine- and coarse-grain patterns of variability. We implement this concept in ACTIONet, a comprehensive framework that combines archetypal analysis and manifold learning to provide a ready-to-use analytical approach for multiresolution single-cell state characterization. ACTIONet provides a robust, reproducible, and highly interpretable single-cell analysis platform that couples dominant pattern discovery with a corresponding structural representation of the cell state landscape. Using multiple synthetic and real data sets, we demonstrate ACTIONet's superior performance relative to existing alternatives. We use ACTIONet to integrate and annotate cells across three human cortex data sets. Through integrative comparative analysis, we define a consensus vocabulary and a consistent set of gene signatures discriminating against the transcriptomic cell types and subtypes of the human prefrontal cortex.

Identifying the molecular systems that influence cognitive resilience to Alzheimer's disease in genetically diverse mice

Individual differences in cognitive decline during normal aging and Alzheimer's disease (AD) are common, but the molecular mechanisms underlying these distinct outcomes are not fully understood. We utilized a combination of genetic, molecular, and behavioral data from a mouse population designed to model human variation in cognitive outcomes to search for the molecular mechanisms behind this population-wide variation. Specifically, we used a systems genetics approach to relate gene expression to cognitive outcomes during AD and normal aging. Statistical causal-inference Bayesian modeling was used to model systematic genetic perturbations matched with cognitive data that identified astrocyte and microglia molecular networks as drivers of cognitive resilience to AD. Using genetic mapping, we identified Fgf2 as a potential regulator of the astrocyte network associated with individual differences in short-term memory. We also identified several immune genes as regulators of a microglia network associated with individual differences in long-term memory, which was partly mediated by amyloid burden. Finally, significant overlap between mouse and two different human coexpression networks provided strong evidence of translational relevance for the genetically diverse AD-BXD panel as a model of late-onset AD. Together, this work identified two candidate molecular pathways enriched for microglia and astrocyte genes that serve as causal AD cognitive biomarkers, and provided a greater understanding of processes that modulate individual and population-wide differences in cognitive outcomes during AD.

Genetics of Gene Expression in the Aging Human Brain Reveal TDP-43 Proteinopathy Pathophysiology

Here, we perform a genome-wide screen for variants that regulate the expression of gene co-expression modules in the aging human brain; we discover and replicate such variants in the TMEM106B and RBFOX1 loci. The TMEM106B haplotype is known to influence the accumulation of TAR DNA-binding protein 43 kDa (TDP-43) proteinopathy, and the haplotype's large-scale transcriptomic effects include the dysregulation of lysosomal genes and alterations in synaptic gene splicing that are also seen in the pathophysiology of TDP-43 proteinopathy. Further, a variant near GRN, another TDP-43 proteinopathy susceptibility gene, shows concordant effects with the TMEM106B haplotype. Leveraging neuropathology data from the same participants, we also show that TMEM106B and APOE-amyloid-β effects converge to alter myelination and lysosomal gene expression, which then contributes to TDP-43 accumulation. These results advance our mechanistic understanding of the TMEM106B TDP-43 risk haplotype and uncover a transcriptional program that mediates the converging effects of APOE-amyloid-β and TMEM106B on TDP-43 aggregation in older adults.

Deconvolving the contributions of cell-type heterogeneity on cortical gene expression

Complexity of cell-type composition has created much skepticism surrounding the interpretation of bulk tissue transcriptomic studies. Recent studies have shown that deconvolution algorithms can be applied to computationally estimate cell-type proportions from gene expression data of bulk blood samples, but their performance when applied to brain tissue is unclear. Here, we have generated an immunohistochemistry (IHC) dataset for five major cell-types from brain tissue of 70 individuals, who also have bulk cortical gene expression data. With the IHC data as the benchmark, this resource enables quantitative assessment of deconvolution algorithms for brain tissue. We apply existing deconvolution algorithms to brain tissue by using marker sets derived from human brain single cell and cell-sorted RNA-seq data. We show that these algorithms can indeed produce informative estimates of constituent cell-type proportions. In fact, neuronal subpopulations can also be estimated from bulk brain tissue samples. Further, we show that including the cell-type proportion estimates as confounding factors is important for reducing false associations between Alzheimer's disease phenotypes and gene expression. Lastly, we demonstrate that using more accurate marker sets can substantially improve statistical power in detecting cell-type specific expression quantitative trait loci (eQTLs).

Meta-Analysis of the Alzheimer's Disease Human Brain Transcriptome and Functional Dissection in Mouse Models

We present a consensus atlas of the human brain transcriptome in Alzheimer's disease (AD), based on meta-analysis of differential gene expression in 2,114 postmortem samples. We discover 30 brain coexpression modules from seven regions as the major source of AD transcriptional perturbations. We next examine overlap with 251 brain differentially expressed gene sets from mouse models of AD and other neurodegenerative disorders. Human-mouse overlaps highlight responses to amyloid versus tau pathology and reveal age- and sex-dependent expression signatures for disease progression. Human coexpression modules enriched for neuronal and/or microglial genes broadly overlap with mouse models of AD, Huntington's disease, amyotrophic lateral sclerosis, and aging. Other human coexpression modules, including those implicated in proteostasis, are not activated in AD models but rather following other, unexpected genetic manipulations. Our results comprise a cross-species resource, highlighting transcriptional networks altered by human brain pathophysiology and identifying correspondences with mouse models for AD preclinical studies.

Gene expression and DNA methylation are extensively coordinated with MRI-based brain microstructural characteristics

Cognitive function relies on both molecular levels and cellular structures. However, systematic relationships between these two components of cognitive function, and their joint contribution to disease, are largely unknown. We utilize postmortem neuroimaging in tandem with gene expression and DNA methylation, from 222 deeply-phenotyped persons in a longitudinal aging cohort. Expression of hundreds of genes and methylation at thousands of loci are related to the microstructure of extensive regions of this same set of brains, as assessed by MRI. The genes linked to brain microstructure perform functions related to cell motility, transcriptional regulation and nuclear processes, and are selectively associated with Alzheimer’s phenotypes. Similar methodology can be applied to other diseases to identify their joint molecular and structural basis, or to infer molecular levels in the brain on the basis of neuroimaging for precision medicine applications.

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.

Classifying tumors by supervised network propagation

Motivation

Network propagation has been widely used to aggregate and amplify the effects of tumor mutations using knowledge of molecular interaction networks. However, propagating mutations through interactions irrelevant to cancer leads to erosion of pathway signals and complicates the identification of cancer subtypes.

Results

To address this problem we introduce a propagation algorithm, Network-Based Supervised Stratification (NBS²), which learns the mutated subnetworks underlying tumor subtypes using a supervised approach. Given an annotated molecular network and reference tumor mutation profiles for which subtypes have been predefined, NBS² is trained by adjusting the weights on interaction features such that network propagation best recovers the provided subtypes. After training, weights are fixed such that mutation profiles of new tumors can be accurately classified. We evaluate NBS² on breast and glioblastoma tumors, demonstrating that it outperforms the best network-based approaches in classifying tumors to known subtypes for these diseases. By interpreting the interaction weights, we highlight characteristic molecular pathways driving selected subtypes.

Availability and implementation

The NBS² package is freely available at: https://github.com/wzhang1984/NBSS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Exploring the psychology of suicidal ideation: A theory driven network analysis

Two leading theories within the field of suicide prevention are the interpersonal psychological theory of suicidal behaviour (IPT) and the integrated motivational-volitional (IMV) model. The IPT posits that suicidal thoughts emerge from high levels of perceived burdensomeness and thwarted belongingness. The IMV model is a multivariate framework that conceptualizes defeat and entrapment as key drivers of suicide ideation. We applied network analysis to cross-sectional data collected as part of the Scottish Wellbeing Study, in which a nationally representative sample of 3508 young adults (18-34 years) completed a battery of psychological measures. Network analysis can help us to understand how the different theoretical components interact and how they relate to suicide ideation. Within a network that included only the core factors from both models, internal entrapment and perceived burdensomeness were most strongly related to suicide ideation. The core constructs defeat, external entrapment and thwarted belonginess were mainly related to other factors than suicide ideation. Within the network of all available psychological factors, 12 of the 20 factors were uniquely related to suicide ideation, with perceived burdensomeness, internal entrapment, depressive symptoms and history of suicide ideation explaining the most variance. None of the factors was isolated, and we identified four larger clusters: mental wellbeing, interpersonal needs, personality, and suicide-related factors. Overall, the results suggest that relationships between suicide ideation and psychological risk factors are complex, with some factors contributing direct risk, and others having indirect impact.

Neuropathological correlates and genetic architecture of microglial activation in elderly human brain

Microglia, the resident immune cells of the brain, have important roles in brain health. However, little is known about the regulation and consequences of microglial activation in the aging human brain. Here we report that the proportion of morphologically activated microglia (PAM) in postmortem cortical tissue is strongly associated with β-amyloid, tau-related neuropathology, and the rate of cognitive decline. Effect sizes for PAM measures are substantial, comparable to that of APOE ε4, the strongest genetic risk factor for Alzheimer's disease, and mediation models support an upstream role for microglial activation in Alzheimer's disease via accumulation of tau. Further, we identify a common variant (rs2997325) influencing PAM that also affects in vivo microglial activation measured by [11C]-PBR28 PET in an independent cohort. Thus, our analyses begin to uncover pathways regulating resident neuroinflammation and identify overlaps of PAM's genetic architecture with those of Alzheimer's disease and several other traits.

Mapping the human brain's cortical-subcortical functional network organization

Understanding complex systems such as the human brain requires characterization of the system's architecture across multiple levels of organization - from neurons, to local circuits, to brain regions, and ultimately large-scale brain networks. Here we focus on characterizing the human brain's large-scale network organization, as it provides an overall framework for the organization of all other levels. We developed a highly principled approach to identify cortical network communities at the level of functional systems, calibrating our community detection algorithm using extremely well-established sensory and motor systems as guides. Building on previous network partitions, we replicated and expanded upon well-known and recently-identified networks, including several higher-order cognitive networks such as a left-lateralized language network. We expanded these cortical networks to subcortex, revealing 358 highly-organized subcortical parcels that take part in forming whole-brain functional networks. Notably, the identified subcortical parcels are similar in number to a recent estimate of the number of cortical parcels (360). This whole-brain network atlas - released as an open resource for the neuroscience community - places all brain structures across both cortex and subcortex into a single large-scale functional framework, with the potential to facilitate a variety of studies investigating large-scale functional networks in health and disease.

Genome-wide RNAseq study of the molecular mechanisms underlying microglia activation in response to pathological tau perturbation in the rTg4510 tau transgenic animal model

BACKGROUND

Activation of microglia, the resident immune cells of the central nervous system, is a prominent pathological hallmark of Alzheimer's disease (AD). However, the gene expression changes underlying microglia activation in response to tau pathology remain elusive. Furthermore, it is not clear how murine gene expression changes relate to human gene expression networks.

METHODS

Microglia cells were isolated from rTg4510 tau transgenic mice and gene expression was profiled using RNA sequencing. Four age groups of mice (2-, 4-, 6-, and 8-months) were analyzed to capture longitudinal gene expression changes that correspond to varying levels of pathology, from minimal tau accumulation to massive neuronal loss. Statistical and system biology approaches were used to analyze the genes and pathways that underlie microglia activation. Differentially expressed genes were compared to human brain co-expression networks.

RESULTS

Statistical analysis of RNAseq data indicated that more than 4000 genes were differentially expressed in rTg4510 microglia compared to wild type microglia, with the majority of gene expression changes occurring between 2- and 4-months of age. These genes belong to four major clusters based on their temporal expression pattern. Genes involved in innate immunity were continuously up-regulated, whereas genes involved in the glutamatergic synapse were down-regulated. Up-regulated innate inflammatory pathways included NF-κB signaling, cytokine-cytokine receptor interaction, lysosome, oxidative phosphorylation, and phagosome. NF-κB and cytokine signaling were among the earliest pathways activated, likely driven by the RELA, STAT1 and STAT6 transcription factors. The expression of many AD associated genes such as APOE and TREM2 was also altered in rTg4510 microglia cells. Differentially expressed genes in rTg4510 microglia were enriched in human neurodegenerative disease associated pathways, including Alzheimer's, Parkinson's, and Huntington's diseases, and highly overlapped with the microglia and endothelial modules of human brain transcriptional co-expression networks.

CONCLUSION

This study revealed temporal transcriptome alterations in microglia cells in response to pathological tau perturbation and provides insight into the molecular changes underlying microglia activation during tau mediated neurodegeneration.

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.

Seasonal plasticity of cognition and related biological measures in adults with and without Alzheimer disease: Analysis of multiple cohorts

BACKGROUND

There are few data concerning the association between season and cognition and its neurobiological correlates in older persons-effects with important translational and therapeutic implications for the diagnosis and treatment of Alzheimer disease (AD). We aimed to measure these effects.

METHODS AND FINDINGS

We analyzed data from 3,353 participants from 3 observational community-based cohort studies of older persons (the Rush Memory and Aging Project [MAP], the Religious Orders Study [ROS], and the Minority Aging Research Study [MARS]) and 2 observational memory-clinic-based cohort studies (Centre de Neurologie Cognitive [CNC] study at Lariboisière Hospital and the Sunnybrook Dementia Study [SDS]). We performed neuropsychological testing and, in subsets of participants, evaluated cerebrospinal fluid AD biomarkers, standardized structured autopsy measures, and/or prefrontal cortex gene expression by RNA sequencing. We examined the association between season and these variables using nested multiple linear and logistic regression models. There was a robust association between season and cognition that was replicated in multiple cohorts (amplitude = 0.14 SD [a measure of the magnitude of seasonal variation relative to overall variability; 95% CI 0.07-0.23], p = 0.007, in the combined MAP, ROS, and MARS cohorts; amplitude = 0.50 SD [95% CI 0.07-0.66], p = 0.017, in the SDS cohort). Average composite global cognitive function was higher in the summer and fall compared to winter and spring, with the difference equivalent in cognitive effect to 4.8 years' difference in age (95% CI 2.1-8.4, p = 0.002). Further, the odds of meeting criteria for mild cognitive impairment or dementia were higher in the winter and spring (odds ratio 1.31 [95% CI 1.10-1.57], p = 0.003). These results were robust against multiple potential confounders including depressive symptoms, sleep, physical activity, and thyroid status and persisted in cases with AD pathology. Moreover, season had a marked effect on cerebrospinal fluid Aβ 42 level (amplitude 0.30 SD [95% CI 0.10-0.64], p = 0.003), which peaked in the summer, and on the brain expression of 4 cognition-associated modules of co-expressed genes (m6: amplitude = 0.44 SD [95% CI 0.21-0.65], p = 0.0021; m13: amplitude = 0.46 SD [95% CI 0.27-0.76], p = 0.0009; m109: amplitude = 0.43 SD [95% CI 0.24-0.67], p = 0.0021; and m122: amplitude 0.46 SD [95% CI 0.20-0.71], p = 0.0012), which were in phase or anti-phase to the rhythms of cognition and which were in turn associated with binding sites for several seasonally rhythmic transcription factors including BCL11A, CTCF, EGR1, MEF2C, and THAP1. Limitations include the evaluation of each participant or sample once per annual cycle, reliance on self-report for measurement of environmental and behavioral factors, and potentially limited generalizability to individuals in equatorial regions or in the southern hemisphere.

CONCLUSIONS

Season has a clinically significant association with cognition and its neurobiological correlates in older adults with and without AD pathology. There may be value in increasing dementia-related clinical resources in the winter and early spring, when symptoms are likely to be most pronounced. Moreover, the persistence of robust seasonal plasticity in cognition and its neurobiological correlates, even in the context of concomitant AD pathology, suggests that targeting environmental or behavioral drivers of seasonal cognitive plasticity, or the key transcription factors and genes identified in this study as potentially mediating these effects, may allow us to substantially improve cognition in adults with and without AD.

Polygenic analysis of inflammatory disease variants and effects on microglia in the aging brain

BACKGROUND

The role of the innate immune system in Alzheimer's disease (AD) and neurodegenerative disease susceptibility has recently been highlighted in genetic studies. However, we do not know whether risk for inflammatory disease predisposes unaffected individuals to late-life cognitive deficits or AD-related neuropathology. We investigated whether genetic risk scores for seven immune diseases and central nervous system traits were related to cognitive decline (nmax = 1601), classical AD neuropathology (nmax = 985), or microglial density (nmax = 184).

METHODS

Longitudinal cognitive decline, postmortem amyloid and tau neuropathology, microglial density, and gene module expression from bulk brain tissue were all measured in participants from two large cohorts (the Rush Religious Orders Study and Memory and Aging Project; ROS/MAP) of elderly subjects (mean age at entry 78 +/- 8.7 years). We analyzed data primarily using robust regression methods. Neuropathologists were blind to clinical data.

RESULTS

The AD genetic risk scores, including and excluding APOE effects, were strongly associated with cognitive decline in all domains (min Puncor = 3.2 × 10- 29). Multiple sclerosis (MS), Parkinson's disease, and schizophrenia risk did not influence cognitive decline in older age, but the rheumatoid arthritis (RA) risk score alone was significantly associated with microglial density after correction (t146 = - 3.88, Puncor = 1.6 × 10- 4). Post-hoc tests found significant effects of the RA genetic risk score in multiple regions and stages of microglial activation (min Puncor = 1.5 × 10- 6). However, these associations were driven by only one or two variants, rather than cumulative polygenicity. Further, individual MS (Pone-sided < 8.4 × 10- 4) and RA (Pone-sided = 3 × 10- 4) variants associated with higher microglial density were also associated with increased expression of brain immune gene modules.

CONCLUSIONS

Our results demonstrate that global risk of inflammatory disease does not strongly influence aging-related cognitive decline but that susceptibility variants that influence peripheral immune function also alter microglial density and immune gene expression in the aging brain, opening a new perspective on the control of microglial and immune responses within the central nervous system. Further study on the molecular mechanisms of peripheral immune disease risk influencing glial cell activation will be required to identify key regulators of these pathways.

A molecular network of the aging human brain provides insights into the pathology and cognitive decline of Alzheimer's disease

There is a need for new therapeutic targets with which to prevent Alzheimer’s disease (AD), a major contributor to aging-related cognitive decline. Here, we report the construction and validation of a molecular network of the aging human frontal cortex. Using RNA sequence data from 478 individuals, we first build a molecular network using modules of coexpressed genes and then relate these modules to AD and its neuropathologic and cognitive endophenotypes. We confirm these associations in two independent AD datasets as well as in epigenomic data. We also illustrate the use of the network in prioritizing amyloid-associated genes for in vitro validation in human neurons and astrocytes. These analyses based on unique cohorts enable us to resolve the role of distinct cortical modules that have a direct effect on the accumulation of AD pathology from those that have a direct effect on cognitive decline, exemplifying a network approach to complex diseases.

Systems biology analysis of RNA sequencing data from the aging human cortex identifies a molecular network which prioritizes groups of genes that influence cognitive decline or neuropathology in Alzheimer’s disease.

A semi-synchronous label propagation algorithm with constraints for community detection in complex networks

Community structure is an important feature of a complex network, where detection of the community structure can shed some light on the properties of such a complex network. Amongst the proposed community detection methods, the label propagation algorithm (LPA) emerges as an effective detection method due to its time efficiency. Despite this advantage in computational time, the performance of LPA is affected by randomness in the algorithm. A modified LPA, called CLPA-GNR, was proposed recently and it succeeded in handling the randomness issues in the LPA. However, it did not remove the tendency for trivial detection in networks with a weak community structure. In this paper, an improved CLPA-GNR is therefore proposed. In the new algorithm, the unassigned and assigned nodes are updated synchronously while the assigned nodes are updated asynchronously. A similarity score, based on the Sørensen-Dice index, is implemented to detect the initial communities and for breaking ties during the propagation process. Constraints are utilised during the label propagation and community merging processes. The performance of the proposed algorithm is evaluated on various benchmark and real-world networks. We find that it is able to avoid trivial detection while showing substantial improvement in the quality of detection.

Deconvoluting simulated metagenomes: the performance of hard- and soft- clustering algorithms applied to metagenomic chromosome conformation capture (3C)

BACKGROUND

Chromosome conformation capture, coupled with high throughput DNA sequencing in protocols like Hi-C and 3C-seq, has been proposed as a viable means of generating data to resolve the genomes of microorganisms living in naturally occuring environments. Metagenomic Hi-C and 3C-seq datasets have begun to emerge, but the feasibility of resolving genomes when closely related organisms (strain-level diversity) are present in the sample has not yet been systematically characterised.

METHODS

We developed a computational simulation pipeline for metagenomic 3C and Hi-C sequencing to evaluate the accuracy of genomic reconstructions at, above, and below an operationally defined species boundary. We simulated datasets and measured accuracy over a wide range of parameters. Five clustering algorithms were evaluated (2 hard, 3 soft) using an adaptation of the extended B-cubed validation measure.

RESULTS

When all genomes in a sample are below 95% sequence identity, all of the tested clustering algorithms performed well. When sequence data contains genomes above 95% identity (our operational definition of strain-level diversity), a naive soft-clustering extension of the Louvain method achieves the highest performance.

DISCUSSION

Previously, only hard-clustering algorithms have been applied to metagenomic 3C and Hi-C data, yet none of these perform well when strain-level diversity exists in a metagenomic sample. Our simple extension of the Louvain method performed the best in these scenarios, however, accuracy remained well below the levels observed for samples without strain-level diversity. Strain resolution is also highly dependent on the amount of available 3C sequence data, suggesting that depth of sequencing must be carefully considered during experimental design. Finally, there appears to be great scope to improve the accuracy of strain resolution through further algorithm development.

Detecting Community Structure by Using a Constrained Label Propagation Algorithm

Community structure is considered one of the most interesting features in complex networks. Many real-world complex systems exhibit community structure, where individuals with similar properties form a community. The identification of communities in a network is important for understanding the structure of said network, in a specific perspective. Thus, community detection in complex networks gained immense interest over the last decade. A lot of community detection methods were proposed, and one of them is the label propagation algorithm (LPA). The simplicity and time efficiency of the LPA make it a popular community detection method. However, the LPA suffers from instability detection due to randomness that is induced in the algorithm. The focus of this paper is to improve the stability and accuracy of the LPA, while retaining its simplicity. Our proposed algorithm will first detect the main communities in a network by using the number of mutual neighbouring nodes. Subsequently, nodes are added into communities by using a constrained LPA. Those constraints are then gradually relaxed until all nodes are assigned into groups. In order to refine the quality of the detected communities, nodes in communities can be switched to another community or removed from their current communities at various stages of the algorithm. We evaluated our algorithm on three types of benchmark networks, namely the Lancichinetti-Fortunato-Radicchi (LFR), Relaxed Caveman (RC) and Girvan-Newman (GN) benchmarks. We also apply the present algorithm to some real-world networks of various sizes. The current results show some promising potential, of the proposed algorithm, in terms of detecting communities accurately. Furthermore, our constrained LPA has a robustness and stability that are significantly better than the simple LPA as it is able to yield deterministic results.