Regional transcriptional architecture of Parkinson's disease pathogenesis and network spread

Distinguish risk genes functioning at presynaptic or postsynaptic regions and key connectomes associated with pathological α-synuclein spreading

Previous studies have suggested that pathological α-synuclein (α-Syn) mainly transmits along the neuronal network, but several key questions remain unanswered: (1) How many and which connections in the connectome are necessary for predicting the progression of pathological α-Syn? (2) How to identify risk gene that affects pathology spreading functioning at presynaptic or postsynaptic regions, and are these genes enriched in different cell types? Here, we addressed these key questions with novel mathematical models. Strikingly, the spreading of pathological α-Syn is predominantly determined by the key subnetworks composed of only 2% of the strongest connections in the connectome. We further explored the genes that are responsible for the selective vulnerability of different brain regions to transmission to distinguish the genes that play roles in presynaptic from those in postsynaptic regions. Those risk genes were significantly enriched in microglial cells of presynaptic regions and neurons of postsynaptic regions. Gene regulatory network analyses were then conducted to identify 'key drivers' of genes responsible for selective vulnerability and overlapping with Parkinson's disease risk genes. By identifying and discriminating between key gene mediators of transmission operating at presynaptic and postsynaptic regions, our study has demonstrated for the first time that these are functionally distinct processes.

Connectome-based biophysical models of pathological protein spreading in neurodegenerative diseases

Neurodegenerative diseases are a group of disorders characterized by progressive degeneration or death of neurons. The complexity of clinical symptoms and irreversibility of disease progression significantly affects individual lives, leading to premature mortality. The prevalence of neurodegenerative diseases keeps increasing, yet the specific pathogenic mechanisms remain incompletely understood and effective treatment strategies are lacking. In recent years, convergent experimental evidence supports the "prion-like transmission" assumption that abnormal proteins induce misfolding of normal proteins, and these misfolded proteins propagate throughout the neural networks to cause neuronal death. To elucidate this dynamic process in vivo from a computational perspective, researchers have proposed three connectome-based biophysical models to simulate the spread of pathological proteins: the Network Diffusion Model, the Epidemic Spreading Model, and the agent-based Susceptible-Infectious-Removed model. These models have demonstrated promising predictive capabilities. This review focuses on the explanations of their fundamental principles and applications. Then, we compare the strengths and weaknesses of the models. Building upon this foundation, we introduce new directions for model optimization and propose a unified framework for the evaluation of connectome-based biophysical models. We expect that this review could lower the entry barrier for researchers in this field, accelerate model optimization, and thereby advance the clinical translation of connectome-based biophysical models.

Predictive spread of obsessive-compulsive disorder pathology using the network diffusion model

An increasing body of studies propose that structural abnormalities begin with focal brain regions then propagate to other regions following the architecture of healthy brain network in neuropsychiatric disorders. However, these findings are untested in obsessive-compulsive disorder (OCD), also showing widespread structural brain abnormalities. In this study, we aimed to investigate whether healthy functional brain network contributed to structural brain abnormalities in OCD. The gray matter morphological abnormalities were obtained in 98 patients with OCD in relative to matched healthy controls (n = 130, HCs). The network diffusion model (NDM) was conducted to identify putative seed regions and patterns of disease propagation from seed regions to other brain regions along the functional network in OCD. The NDM has been proved to succeeded in capturing the trans-neuronal propagation of pathology and even in predicting future longitudinal progression of pathology in neurodegenerative diseases. In this study, when seeding at the right anterior cingulate cortex, the NDM best recapitulated the patterns of gray matter morphological abnormalities, suggesting this region was the most likely seed region. Further analyses revealed that pathology preferentially propagated to higher order brain systems from seed region. For non-seed regions, the arrival time of pathology was negatively correlated with their shortest functional paths to the seed (r = -0.46, p < 0.001). These results suggest that gray matter morphological abnormalities are constrained by healthy brain network and reveal temporal sequencing of pathology progression in OCD.

Spatial transcriptomic patterns underlying amyloid-β and tau pathology are associated with cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We study the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies by leveraging two large independent AD cohorts. We identify AD susceptibility genes and gene modules in a gene co-expression network with expression profiles specifically related to regional vulnerability to Aβ and tau pathologies in AD. In addition, we identify distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. These findings may explain the discordance between regional Aβ and tau pathologies. Finally, we propose an analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations.

Association between gene expression and functional-metabolic architecture in Parkinson's disease

Gene expression plays a critical role in the pathogenesis of Parkinson's disease (PD). How gene expression profiles are correlated with functional-metabolic architecture remains obscure. We enrolled 34 PD patients and 25 age-and-sex-matched healthy controls for simultaneous 18 F-FDG-PET/functional MRI scanning during resting state. We investigated the functional gradients and the ratio of standard uptake value. Principal component analysis was used to further combine the functional gradients and glucose metabolism into functional-metabolic architecture. Using partial least squares (PLS) regression, we introduced the transcriptomic data from the Allen Institute of Brain Sciences to identify gene expression patterns underlying the affected functional-metabolic architecture in PD. Between-group comparisons revealed significantly higher gradient variation in the visual, somatomotor, dorsal attention, frontoparietal, default mode, and subcortical network (pFDR  < .048) in PD. Increased FDG-uptake was found in the somatomotor and ventral attention network while decreased FDG-uptake was found in the visual network (pFDR  < .008). Spatial correlation analysis showed consistently affected patterns of functional gradients and metabolism (p = 2.47 × 10-8 ). PLS analysis and gene ontological analyses further revealed that genes were mainly enriched for metabolic, catabolic, cellular response to ions, and regulation of DNA transcription and RNA biosynthesis. In conclusion, our study provided genetic pathological mechanism to explain imaging-defined brain functional-metabolic architecture of PD.

Connectome-based modelling of neurodegenerative diseases: towards precision medicine and mechanistic insight

Neurodegenerative diseases are the most common cause of dementia. Although their underlying molecular pathologies have been identified, there is substantial heterogeneity in the patterns of progressive brain alterations across and within these diseases. Recent advances in neuroimaging methods have revealed that pathological proteins accumulate along specific macroscale brain networks, implicating the network architecture of the brain in the system-level pathophysiology of neurodegenerative diseases. However, the extent to which 'network-based neurodegeneration' applies across the wide range of neurodegenerative disorders remains unclear. Here, we discuss the state-of-the-art of neuroimaging-based connectomics for the mapping and prediction of neurodegenerative processes. We review findings supporting brain networks as passive conduits through which pathological proteins spread. As an alternative view, we also discuss complementary work suggesting that network alterations actively modulate the spreading of pathological proteins between connected brain regions. We conclude this Perspective by proposing an integrative framework in which connectome-based models can be advanced along three dimensions of innovation: incorporating parameters that modulate propagation behaviour on the basis of measurable biological features; building patient-tailored models that use individual-level information and allowing model parameters to interact dynamically over time. We discuss promises and pitfalls of these strategies for improving disease insights and moving towards precision medicine.

Spatial transcriptomic patterns underlying regional vulnerability to amyloid-β and tau pathologies and their relationships to cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We studied the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies leveraging two large independent cohorts (n = 715) of participants along the AD continuum. We identified several AD susceptibility genes and gene modules in a gene co-expression network with expression profiles related to regional vulnerability to Aβ and tau pathologies in AD. In particular, we found that the positive APOE -to-tau association was only seen in the AD cohort, whereas patients with AD and frontotemporal dementia shared similar positive MAPT -to-tau association. Some AD candidate genes showed sex-dependent negative gene-to-Aβ and gene-to-tau associations. In addition, we identified distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. Finally, we proposed a novel analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations. Taken together, our study identified distinct gene expression profiles and biochemical pathways that may explain the discordance between regional Aβ and tau pathologies, and filled the gap between gene-to-pathology associations and cognitive dysfunction in individual AD patients that may ultimately help identify novel personalized pathogenetic biomarkers and therapeutic targets.

One Sentence Summary

We identified replicable cognition-related associations between regional gene expression profiles and selectively regional vulnerability to amyloid-β and tau pathologies in AD.

Mitochondrial function-associated genes underlie cortical atrophy in prodromal synucleinopathies

Isolated rapid eye movement sleep behaviour disorder (iRBD) is a sleep disorder characterized by the loss of rapid eye movement sleep muscle atonia and the appearance of abnormal movements and vocalizations during rapid eye movement sleep. It is a strong marker of incipient synucleinopathy such as dementia with Lewy bodies and Parkinson's disease. Patients with iRBD already show brain changes that are reminiscent of manifest synucleinopathies including brain atrophy. However, the mechanisms underlying the development of this atrophy remain poorly understood. In this study, we performed cutting-edge imaging transcriptomics and comprehensive spatial mapping analyses in a multicentric cohort of 171 polysomnography-confirmed iRBD patients [67.7 ± 6.6 (49-87) years; 83% men] and 238 healthy controls [66.6 ± 7.9 (41-88) years; 77% men] with T1-weighted MRI to investigate the gene expression and connectivity patterns associated with changes in cortical thickness and surface area in iRBD. Partial least squares regression was performed to identify the gene expression patterns underlying cortical changes in iRBD. Gene set enrichment analysis and virtual histology were then done to assess the biological processes, cellular components, human disease gene terms, and cell types enriched in these gene expression patterns. We then used structural and functional neighbourhood analyses to assess whether the atrophy patterns in iRBD were constrained by the brain's structural and functional connectome. Moreover, we used comprehensive spatial mapping analyses to assess the specific neurotransmitter systems, functional networks, cytoarchitectonic classes, and cognitive brain systems associated with cortical changes in iRBD. All comparisons were tested against null models that preserved spatial autocorrelation between brain regions and compared to Alzheimer's disease to assess the specificity of findings to synucleinopathies. We found that genes involved in mitochondrial function and macroautophagy were the strongest contributors to the cortical thinning occurring in iRBD. Moreover, we demonstrated that cortical thinning was constrained by the brain's structural and functional connectome and that it mapped onto specific networks involved in motor and planning functions. In contrast with cortical thickness, changes in cortical surface area were related to distinct genes, namely genes involved in the inflammatory response, and to different spatial mapping patterns. The gene expression and connectivity patterns associated with iRBD were all distinct from those observed in Alzheimer's disease. In summary, this study demonstrates that the development of brain atrophy in synucleinopathies is constrained by specific genes and networks.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

Modeling seeding and neuroanatomic spread of pathology in amyotrophic lateral sclerosis

The neurodegenerative disorder amyotrophic lateral sclerosis (ALS) is characterized by the progressive loss of upper and lower motor neurons, with pathological involvement of cerebral motor and extra-motor areas in a clinicopathological spectrum with frontotemporal dementia (FTD). A key unresolved issue is how the non-random distribution of pathology in ALS reflects differential network vulnerability, including molecular factors such as regional gene expression, or preferential spread of pathology via anatomical connections. A system of histopathological staging of ALS based on the regional burden of TDP-43 pathology observed in postmortem brains has been supported to some extent by analysis of distribution of in vivo structural MRI changes. In this paper, computational modeling using a Network Diffusion Model (NDM) was used to investigate whether a process of focal pathological 'seeding' followed by structural network-based spread recapitulated postmortem histopathological staging and, secondly, whether this had any correlation to the pattern of expression of a panel of genes implicated in ALS across the healthy brain. Regionally parcellated T1-weighted MRI data from ALS patients (baseline n=79) was studied in relation to a healthy control structural connectome and a database of associated regional cerebral gene expression. The NDM provided strong support for a structural network-based basis for regional pathological spread in ALS, but no simple relationship to the spatial distribution of ALS-related genes in the healthy brain. Interestingly, OPTN gene was identified as a significant but a weaker non-NDM contributor within the network-gene interaction model (LASSO). Intriguingly, the critical seed regions for spread within the model were not within the primary motor cortex but basal ganglia, thalamus and insula, where NDM recapitulated aspects of the postmortem histopathological staging system. Within the ALS-FTD clinicopathological spectrum, non-primary motor structures may be among the earliest sites of cerebral pathology.

Initiation and progression of α-synuclein pathology in Parkinson’s disease

α-Synuclein aggregation is a critical molecular process that underpins the pathogenesis of Parkinson’s disease. Aggregates may originate at synaptic terminals as a consequence of aberrant interactions between α-synuclein and lipids or evasion of proteostatic defences. The nature of these interactions is likely to influence the emergence of conformers or strains that in turn could explain the clinical heterogeneity of Parkinson’s disease and related α-synucleinopathies. For neurodegeneration to occur, α-synuclein assemblies need to exhibit seeding competency, i.e. ability to template further aggregation, and toxicity which is at least partly mediated by interference with synaptic vesicle or organelle homeostasis. Given the dynamic and reversible conformational plasticity of α-synuclein, it is possible that seeding competency and cellular toxicity are mediated by assemblies of different structure or size along this continuum. It is currently unknown which α-synuclein assemblies are the most relevant to the human condition but recent advances in the cryo-electron microscopic characterisation of brain-derived fibrils and their assessment in stem cell derived and animal models are likely to facilitate the development of precision therapies or biomarkers. This review summarises the main principles of α-synuclein aggregate initiation and propagation in model systems, and their relevance to clinical translation.

Statistical testing in transcriptomic-neuroimaging studies: A how-to and evaluation of methods assessing spatial and gene specificity

Multiscale integration of gene transcriptomic and neuroimaging data is becoming a widely used approach for exploring the molecular underpinnings of large‐scale brain organization in health and disease. Proper statistical evaluation of determined associations between imaging‐based phenotypic and transcriptomic data is key in these explorations, in particular to establish whether observed associations exceed “chance level” of random, nonspecific effects. Recent approaches have shown the importance of statistical models that can correct for spatial autocorrelation effects in the data to avoid inflation of reported statistics. Here, we discuss the need for examination of a second category of statistical models in transcriptomic‐neuroimaging analyses, namely those that can provide “gene specificity.” By means of a couple of simple examples of commonly performed transcriptomic‐neuroimaging analyses, we illustrate some of the potentials and challenges of transcriptomic‐imaging analyses, showing that providing gene specificity on observed transcriptomic‐neuroimaging effects is of high importance to avoid reports of nonspecific effects. Through means of simulations we show that the rate of reported nonspecific effects (i.e., effects that cannot be specifically linked to a specific gene or gene‐set) can run as high as 60%, with only less than 5% of transcriptomic‐neuroimaging associations observed through ordinary linear regression analyses showing both spatial and gene specificity. We provide a discussion, a tutorial, and an easy‐to‐use toolbox for the different options of null models in transcriptomic‐neuroimaging analyses.

Imaging Transcriptomics of Brain Disorders

Noninvasive neuroimaging is a powerful tool for quantifying diverse aspects of brain structure and function in vivo, and it has been used extensively to map the neural changes associated with various brain disorders. However, most neuroimaging techniques offer only indirect measures of underlying pathological mechanisms. The recent development of anatomically comprehensive gene expression atlases has opened new opportunities for studying the transcriptional correlates of noninvasively measured neural phenotypes, offering a rich framework for evaluating pathophysiological hypotheses and putative mechanisms. Here, we provide an overview of some fundamental methods in imaging transcriptomics and outline their application to understanding brain disorders of neurodevelopment, adulthood, and neurodegeneration. Converging evidence indicates that spatial variations in gene expression are linked to normative changes in brain structure during age-related maturation and neurodegeneration that are in part associated with cell-specific gene expression markers of gene expression. Transcriptional correlates of disorder-related neuroimaging phenotypes are also linked to transcriptionally dysregulated genes identified in ex vivo analyses of patient brains. Modeling studies demonstrate that spatial patterns of gene expression are involved in regional vulnerability to neurodegeneration and the spread of disease across the brain. This growing body of work supports the utility of transcriptional atlases in testing hypotheses about the molecular mechanism driving disease-related changes in macroscopic neuroimaging phenotypes.

Convergent structural network and gene signatures for MRgFUS thalamotomy in patients with Parkinson's disease

MRgFUS has just been made available for the 1.7 million Parkinson's disease patients in China. Despite its non-invasive and rapid therapeutic advantages for involuntary tremor, some concerns have emerged about outcomes variability, non-specificity, and side-effects, as little is known about its impact on the long-term plasticity of brain structure. We sought to dissect the characteristics of long-term changes in brain structure caused by MRgFUS lesion and explored potential biological mechanisms. One-year multimodal imaging follow-ups were conducted for nine tremor-dominant Parkinson's disease patients undergoing unilateral MRgFUS thalamotomy. A structural connectivity map was generated for each patient to analyze dynamic changes in brain structure. The human brain transcriptome was extracted and spatially registered for connectivity vulnerability. Genetic functional enrichment analysis was performed and further clarified using in vivo emission computed tomography data. MRgFUS not only abolished tremors but also significantly disrupted the brain network topology. Network-based statistics identified a U-shape MRgFUS-sensitive subnetwork reflective of hand tremor recovery and surgical process, accompanied by relevant cerebral blood flow and gray matter alteration. Using human brain gene expression data, we observed that dopaminergic signatures were responsible for the preferential vulnerability associated with these architectural alterations. Additional PET/SPECT data not only validated these gene signatures, but also suggested that structural alteration was significantly correlated with D1 and D2 receptors, DAT, and F-DOPA measures. There was a long-term dynamic loop between structural alteration and dopaminergic signature for MRgFUS thalamotomy, which may be closely related to the long-term improvements in clinical tremor.

Activated microglia facilitate the transmission of α-synuclein in Parkinson's disease

Parkinson's disease (PD) is characterized by the loss of dopaminergic neurons in the substantia nigra pars compacta and abnormal aggregates of α-synuclein protein called Lewy bodies. To date, there is no drug that can definitely slow down or stop the progression of this disease. The discovery of the cell-to-cell transmission of pathologic α-synuclein seeds offers the possibility to explore novel treatment strategies to prevent the spread of α-synuclein, with the purpose of slowing down the progression of PD in its tracks. Although recent studies have made tremendous progress in understanding how α-synuclein spreads throughout the brain, neuroinflammation seems to play a crucial role in the development of α-synuclein pathology in PD. The activation of microglia, one of the hallmarks of the neuroinflammatory process, is suggested to influence the neuron-to-neuron transmission of α-synuclein. This review summarizes how activated microglia facilitate this process, and focuses on the following mechanisms including the activation of microglia in PD, the reduced ability of activated microglia to clear α-synuclein and increased migratory capacity of microglia in PD, as well as the cooperation between microglia and exosomes in mediating α-synuclein release and propagation. In conclusion, this article help collate information on microglia in-relation to PD.

Graph Models of Pathology Spread in Alzheimer's Disease: An Alternative to Conventional Graph Theoretic Analysis

Background: Graph theory and connectomics are new techniques for uncovering disease-induced changes in the brain's structural network. Most prior studied have focused on network statistics as biomarkers of disease. However, an emerging body of work involves exploring how the network serves as a conduit for the propagation of disease factors in the brain and has successfully mapped the functional and pathological consequences of disease propagation. In Alzheimer's disease (AD), progressive deposition of misfolded proteins amyloid and tau is well-known to follow fiber projections, under a "prion-like" trans-neuronal transmission mechanism, through which misfolded proteins cascade along neuronal pathways, giving rise to network spread. Methods: In this review, we survey the state of the art in mathematical modeling of connectome-mediated pathology spread in AD. Then we address several open questions that are amenable to mathematically precise parsimonious modeling of pathophysiological processes, extrapolated to the whole brain. We specifically identify current formal models of how misfolded proteins are produced, aggregate, and disseminate in brain circuits, and attempt to understand how this process leads to stereotyped progression in Alzheimer's and other related diseases. Conclusion: This review serves to unify current efforts in modeling of AD progression that together have the potential to explain observed phenomena and serve as a test-bed for future hypothesis generation and testing in silico. Impact statement Graph theory is a powerful new approach that is transforming the study of brain processes. There do not exist many focused reviews of the subfield of graph modeling of how Alzheimer's and other dementias propagate within the brain network, and how these processes can be mapped mathematically. By providing timely and topical review of this subfield, we fill a critical gap in the community and present a unified view that can serve as an in silico test-bed for future hypothesis generation and testing. We also point to several open and unaddressed questions and controversies that future practitioners can tackle.

Integrated transcriptomic and neuroimaging brain model decodes biological mechanisms in aging and Alzheimer's disease

Both healthy aging and Alzheimer's disease (AD) are characterized by concurrent alterations in several biological factors. However, generative brain models of aging and AD are limited in incorporating the measures of these biological factors at different spatial resolutions. Here, we propose a personalized bottom-up spatiotemporal brain model that accounts for the direct interplay between hundreds of RNA transcripts and multiple macroscopic neuroimaging modalities (PET, MRI). In normal elderly and AD participants, the model identifies top genes modulating tau and amyloid-β burdens, vascular flow, glucose metabolism, functional activity, and atrophy to drive cognitive decline. The results also revealed that AD and healthy aging share specific biological mechanisms, even though AD is a separate entity with considerably more altered pathways. Overall, this personalized model offers novel insights into the multiscale alterations in the elderly brain, with important implications for identifying effective genetic targets for extending healthy aging and treating AD progression.

Differential patterns of gray matter volumes and associated gene expression profiles in cognitively-defined Alzheimer's disease subgroups

The clinical presentation of Alzheimer's disease (AD) varies widely across individuals but the neurobiological mechanisms underlying this heterogeneity are largely unknown. Here, we compared regional gray matter (GM) volumes and associated gene expression profiles between cognitively-defined subgroups of amyloid-β positive individuals clinically diagnosed with AD dementia (age: 66 ± 7, 47% male, MMSE: 21 ± 5). All participants underwent neuropsychological assessment with tests covering memory, executive-functioning, language and visuospatial-functioning domains. Subgroup classification was achieved using a psychometric framework that assesses which cognitive domain shows substantial relative impairment compared to the intra-individual average across domains, which yielded the following subgroups in our sample; AD-Memory (n = 41), AD-Executive (n = 117), AD-Language (n = 33), AD-Visuospatial (n = 171). We performed voxel-wise contrasts of GM volumes derived from 3Tesla structural MRI between subgroups and controls (n = 127, age 58 ± 9, 42% male, MMSE 29 ± 1), and observed that differences in regional GM volumes compared to controls closely matched the respective cognitive profiles. Specifically, we detected lower medial temporal lobe GM volumes in AD-Memory, lower fronto-parietal GM volumes in AD-Executive, asymmetric GM volumes in the temporal lobe (left < right) in AD-Language, and lower GM volumes in posterior areas in AD-Visuospatial. In order to examine possible biological drivers of these differences in regional GM volumes, we correlated subgroup-specific regional GM volumes to brain-wide gene expression profiles based on a stereotactic characterization of the transcriptional architecture of the human brain as provided by the Allen human brain atlas. Gene-set enrichment analyses revealed that variations in regional expression of genes involved in processes like mitochondrial respiration and metabolism of proteins were associated with patterns of regional GM volume across multiple subgroups. Other gene expression vs GM volume-associations were only detected in particular subgroups, e.g., genes involved in the cell cycle for AD-Memory, specific sets of genes related to protein metabolism in AD-Language, and genes associated with modification of gene expression in AD-Visuospatial. We conclude that cognitively-defined AD subgroups show neurobiological differences, and distinct biological pathways may be involved in the emergence of these differences.

Network model of pathology spread recapitulates neurodegeneration and selective vulnerability in Huntington's Disease

Huntington's Disease (HD), an autosomal dominant genetic disorder caused by a mutation in the Huntingtin gene (HTT), displays a stereotyped topography in the human brain and a stereotyped progression, initially appearing in the striatum. Like other degenerative diseases, spatial topography of HD is divorced from where implicated genes are expressed, a dissociation whose mechanistic underpinning is not currently understood. Cell autonomous molecular factors characterized by gene expression signatures, including proteolytic and post translational modifications, play a role in vulnerability to disease. Non-autonomous mechanisms, likely involving the brain's anatomic or functional connectivity patterns, might also be responsible for selective vulnerability in HD. Leveraging a large dataset of 635 subjects from a multinational study, this paper tests various cell-autonomous and non-autonomous models that can explain HD topography. We test whether the expression patterns of implicated genes is sufficient to explain regional HD atrophy, or whether the network transmission of protein products is required to explain them. We find that network models are capable of predicting, to a high degree, observed atrophy in human subjects. Lastly, we propose a model of anterograde network transmission, and show that it is the most parsimonious yet most likely to explain observed atrophy patterns in HD. Collectively, these data indicate that pathology spread in HD may be mediated by the brain's intrinsic structural network organization. This is the first study to systematically and quantitatively test multiple hypotheses of pathology spread in living human subjects with HD.

Disentangling the Molecular Pathways of Parkinson's Disease using Multiscale Network Modeling

Parkinson's disease (PD) is a complex neurodegenerative disorder. The identification of genetic variants has shed light on the molecular pathways for inherited PD, while the disease mechanism for idiopathic PD remains elusive, partly due to a lack of robust tools. The complexity of PD arises from the heterogeneity of clinical symptoms, pathologies, environmental insults contributing to the disease, and disease comorbidities. Molecular networks have been increasingly used to identify molecular pathways and drug targets in complex human diseases. Here, we review recent advances in molecular network approaches and their application to PD. We discuss how network modeling can predict functions of PD genetic risk factors through network context and assist in the discovery of network-based therapeutics for neurodegenerative diseases.

Multimodal phenotypic axes of Parkinson's disease

Individuals with Parkinson’s disease present with a complex clinical phenotype, encompassing sleep, motor, cognitive, and affective disturbances. However, characterizations of PD are typically made for the “average” patient, ignoring patient heterogeneity and obscuring important individual differences. Modern large-scale data sharing efforts provide a unique opportunity to precisely investigate individual patient characteristics, but there exists no analytic framework for comprehensively integrating data modalities. Here we apply an unsupervised learning method—similarity network fusion—to objectively integrate MRI morphometry, dopamine active transporter binding, protein assays, and clinical measurements from n = 186 individuals with de novo Parkinson’s disease from the Parkinson’s Progression Markers Initiative. We show that multimodal fusion captures inter-dependencies among data modalities that would otherwise be overlooked by field standard techniques like data concatenation. We then examine how patient subgroups derived from the fused data map onto clinical phenotypes, and how neuroimaging data is critical to this delineation. Finally, we identify a compact set of phenotypic axes that span the patient population, demonstrating that this continuous, low-dimensional projection of individual patients presents a more parsimonious representation of heterogeneity in the sample compared to discrete biotypes. Altogether, these findings showcase the potential of similarity network fusion for combining multimodal data in heterogeneous patient populations.

Differences in network controllability and regional gene expression underlie hallucinations in Parkinson's disease

Visual hallucinations are common in Parkinson's disease and are associated with poorer prognosis. Imaging studies show white matter loss and functional connectivity changes with Parkinson's visual hallucinations, but the biological factors underlying selective vulnerability of affected parts of the brain network are unknown. Recent models for Parkinson's disease hallucinations suggest they arise due to a shift in the relative effects of different networks. Understanding how structural connectivity affects the interplay between networks will provide important mechanistic insights. To address this, we investigated the structural connectivity changes that accompany visual hallucinations in Parkinson's disease and the organizational and gene expression characteristics of the preferentially affected areas of the network. We performed diffusion-weighted imaging in 100 patients with Parkinson's disease (81 without hallucinations, 19 with visual hallucinations) and 34 healthy age-matched controls. We used network-based statistics to identify changes in structural connectivity in Parkinson's disease patients with hallucinations and performed an analysis of controllability, an emerging technique that allows quantification of the influence a brain region has across the rest of the network. Using these techniques, we identified a subnetwork of reduced connectivity in Parkinson's disease hallucinations. We then used the Allen Institute for Brain Sciences human transcriptome atlas to identify regional gene expression patterns associated with affected areas of the network. Within this network, Parkinson's disease patients with hallucinations showed reduced controllability (less influence over other brain regions), than Parkinson's disease patients without hallucinations and controls. This subnetwork appears to be critical for overall brain integration, as even in controls, nodes with high controllability were more likely to be within the subnetwork. Gene expression analysis of gene modules related to the affected subnetwork revealed that down-weighted genes were most significantly enriched in genes related to mRNA and chromosome metabolic processes (with enrichment in oligodendrocytes) and upweighted genes to protein localization (with enrichment in neuronal cells). Our findings provide insights into how hallucinations are generated, with breakdown of a key structural subnetwork that exerts control across distributed brain regions. Expression of genes related to mRNA metabolism and membrane localization may be implicated, providing potential therapeutic targets.

Dementia risk in Parkinson's disease is associated with interhemispheric connectivity loss and determined by regional gene expression

Parkinson's dementia is a common and devastating part of Parkinson's disease. Whilst timing and severity vary, dementia in Parkinson's is often preceded by visual dysfunction. White matter changes, representing axonal loss, occur early in the disease process. Clarifying which white matter connections are affected in Parkinson's with visual dysfunction and why specific connections are vulnerable will provide important mechanistic insights. Here, we use diffusion tractography in 100 Parkinson's patients (33 low visual performers) and 34 controls to identify patterns of connectivity loss in Parkinson's with visual dysfunction. We examine the relationship between regional transcription and connectivity loss, using the Allen Institute for Brain Science transcriptome atlas. We show that interhemispheric connections are preferentially affected in Parkinson's low visual performers. Interhemispheric connection loss was associated with downweighted genes related to the smoothened signalling pathway (enriched in glutamatergic neurons) and upweighted metabolic genes. Risk genes for Parkinson's but not Alzheimer's or Dementia with Lewy bodies were over-represented in upweighted genes associated with interhemispheric connection loss. Our findings suggest selective vulnerability in Parkinson's patients at highest risk of dementia (those with visual dysfunction), where differences in gene expression and metabolic dysfunction, affecting longer connections with higher metabolic burden, drive connectivity loss.

Origins of atrophy in Parkinson linked to early onset and local transcription patterns

There is enormous clinical value in inferring the brain regions initially atrophied in Parkinson disease for individual patients and understanding its relationship with clinical and genetic risk factors. The aim of this study is to leverage a new seed-inference algorithm demonstrated for Alzheimer's disease to the Parkinsonian context and to cluster patients in meaningful subgroups based on these incipient atrophy patterns. Instead of testing brain regions separately as the likely initiation site for each patient, we solve an L1-penalized optimization problem that can return a more predictive heterogeneous, multi-locus seed patterns. A cluster analysis of the individual seed patterns reveals two distinct subgroups (S1 versus S2). The S1 subgroup is characterized by the involvement of the brainstem and ventral nuclei, and S2 by cortex and striatum. Post hoc analysis in features not included in the clustering shows significant differences between subgroups regarding age of onset and local transcriptional patterns of Parkinson-related genes. Top genes associated with regional microglial abundance are strongly associated with subgroup S1 but not with S2. Our results suggest two distinct aetiological mechanisms operative in Parkinson disease. The interplay between immune-related genes, lysosomal genes, microglial abundance and atrophy initiation sites may explain why the age of onset for patients in S1 is on average 4.5 years later than for those in S2. We highlight and compare the most prominently affected brain regions for both subgroups. Altogether, our findings may improve current screening strategies for early Parkinson onsetters.

Network diffusion modeling predicts neurodegeneration in traumatic brain injury

Objective

Traumatic brain injury (TBI) is a heterogeneous disease with multiple neurological deficits that evolve over time. It is also associated with an increased incidence of neurodegenerative diseases. Accordingly, clinicians need better tools to predict a patient’s long‐term prognosis.

Methods

Diffusion‐weighted and anatomical MRI data were collected from 17 adolescents (mean age = 15y8mo) with moderate‐to‐severe TBI and 19 healthy controls. Using a network diffusion model (NDM), we examined the effect of progressive deafferentation and gray matter thinning in young TBI patients. Moreover, using a novel automated inference method, we identified several injury epicenters in order to determine the neural degenerative patterns in each TBI patient.

Results

We were able to identify the subject‐specific patterns of degeneration in each patient. In particular, the hippocampus, temporal cortices, and striatum were frequently found to be the epicenters of degeneration across the TBI patients. Orthogonal transformation of the predicted degeneration, using principal component analysis, identified distinct spatial components in the temporal–hippocampal network and the cortico‐striatal network, confirming the vulnerability of these networks to injury. The NDM model, best predictive of the degeneration, was significantly correlated with time since injury, indicating that NDM can potentially capture the pathological progression in the chronic phase of TBI.

Interpretation

These findings suggest that network spread may help explain patterns of distant gray matter thinning, which would be consistent with Wallerian degeneration of the white matter connections (i.e., “diaschisis”) from diffuse axonal injuries and multifocal contusive injuries, and the neurodegenerative patterns of abnormal protein aggregation and transmission, which are hallmarks of brain changes in TBI. NDM approaches could provide highly subject‐specific biomarkers relevant for disease monitoring and personalized therapies in TBI.

Dynamical Role of Pivotal Brain Regions in Parkinson Symptomatology Uncovered with Deep Learning

BACKGROUND

The release of a broad, longitudinal anatomical dataset by the Parkinson's Progression Markers Initiative promoted a surge of machine-learning studies aimed at predicting disease onset and progression. However, the excessive number of features used in these models often conceals their relationship to the Parkinsonian symptomatology.

OBJECTIVES

The aim of this study is two-fold: (i) to predict future motor and cognitive impairments up to four years from brain features acquired at baseline; and (ii) to interpret the role of pivotal brain regions responsible for different symptoms from a neurological viewpoint.

METHODS

We test several deep-learning neural network configurations, and report our best results obtained with an autoencoder deep-learning model, run on a 5-fold cross-validation set. Comparison with Existing Methods: Our approach improves upon results from standard regression and others. It also includes neuroimaging biomarkers as features.

RESULTS

The relative contributions of pivotal brain regions to each impairment change over time, suggesting a dynamical reordering of culprits as the disease progresses. Specifically, the Putamen is initially the most critical region accounting for the overall cognitive state, only being surpassed by the Substantia Nigra in later years. The Pallidum is the first region to influence motor scores, followed by the parahippocampal and ambient gyri, and the anterior orbital gyrus.

CONCLUSIONS

While the causal link between regional brain atrophy and Parkinson symptomatology is poorly understood, our methods demonstrate that the contributions of pivotal regions to cognitive and motor impairments are more dynamical than generally appreciated.