Longitudinal imaging biomarkers are invaluable for understanding the course of neurodegeneration, promising the ability to track disease progression and to detect disease earlier than cross-sectional biomarkers. To properly realize their potential, biomarker trajectory models must be robust to both under-sampling and measurement errors and should be able to integrate multi-modal information to improve trajectory inference and prediction. Here we present a parametric Bayesian multi-task learning based approach to modeling univariate trajectories across subjects that addresses these criteria. Our approach learns multiple subjects' trajectories within a single model that allows for different types of information sharing, that is, coupling, across subjects. It optimizes a combination of uncoupled, fully coupled and kernel coupled models. Kernel-based coupling allows linking subjects' trajectories based on one or more biomarker measures. We demonstrate this using Alzheimer's Disease Neuroimaging Initiative (ADNI) data, where we model longitudinal trajectories of MRI-derived cortical volumes in neurodegeneration, with coupling based on APOE genotype, cerebrospinal fluid (CSF) and amyloid PET-based biomarkers. In addition to detecting established disease effects, we detect disease related changes within the insula that have not received much attention within the literature. Due to its sensitivity in detecting disease effects, its competitive predictive performance and its ability to learn the optimal parameter covariance from data rather than choosing a specific set of random and fixed effects a priori, we propose that our model can be used in place of or in addition to linear mixed effects models when modeling biomarker trajectories. A software implementation of the method is publicly available.

Prior work in late-onset Alzheimer's disease (LOAD) has resulted in discrepant findings as to whether recent consanguinity and outbred autozygosity are associated with LOAD risk. In the current study, we tested the association between consanguinity and outbred autozygosity with LOAD in the largest such analysis to date, in which 20 LOAD GWAS datasets were retrieved through public databases. Our analyses were restricted to eight distinct ethnic groups: African-Caribbean, Ashkenazi-Jewish European, European-Caribbean, French-Canadian, Finnish European, North-Western European, South-Eastern European, and Yoruba African for a total of 21,492 unrelated subjects (11,196 LOAD and 10,296 controls). Recent consanguinity determination was performed using FSuite v1.0.3, according to subjects' ancestral background. The level of autozygosity in the outbred population was assessed by calculating inbreeding estimates based on the proportion (FROH) and the number (NROH) of runs of homozygosity (ROHs). We analyzed all eight ethnic groups using a fixed-effect meta-analysis, which showed a significant association of recent consanguinity with LOAD (N = 21,481; OR = 1.262, P = 3.6 × 10-4), independently of APOE ∗4 (N = 21,468, OR = 1.237, P = 0.002), and years of education (N = 9,257; OR = 1.274, P = 0.020). Autozygosity in the outbred population was also associated with an increased risk of LOAD, both for F ROH (N = 20,237; OR = 1.204, P = 0.030) and N ROH metrics (N = 20,237; OR = 1.019, P = 0.006), independently of APOE ∗4 [(F ROH, N = 20,225; OR = 1.222, P = 0.029) (N ROH, N = 20,225; OR = 1.019, P = 0.007)]. By leveraging the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data, we determined that LOAD subjects do not show an enrichment of rare, risk-enhancing minor homozygote variants compared to the control population. A two-stage recessive GWAS using ADSP data from 201 consanguineous subjects in the discovery phase followed by validation in 10,469 subjects led to the identification of RPH3AL p.A303V (rs117190076) as a rare minor homozygote variant increasing the risk of LOAD [discovery: Genotype Relative Risk (GRR) = 46, P = 2.16 × 10-6; validation: GRR = 1.9, P = 8.0 × 10-4]. These results confirm that recent consanguinity and autozygosity in the outbred population increase risk for LOAD. Subsequent work, with increased samples sizes of consanguineous subjects, should accelerate the discovery of non-additive genetic effects in LOAD.

Cortical cell loss is a core feature of Huntington's disease (HD), beginning many years before clinical motor diagnosis, during the premanifest stage. However, it is unclear how genetic topography relates to cortical cell loss. Here, we explore the biological processes and cell types underlying this relationship and validate these using cell-specific post-mortem data. Eighty premanifest participants on average 15 years from disease onset and 71 controls were included. Using volumetric and diffusion MRI we extracted HD-specific whole brain maps where lower grey matter volume and higher grey matter mean diffusivity, relative to controls, were used as proxies of cortical cell loss. These maps were combined with gene expression data from the Allen Human Brain Atlas (AHBA) to investigate the biological processes relating genetic topography and cortical cell loss. Cortical cell loss was positively correlated with the expression of developmental genes (i.e. higher expression correlated with greater atrophy and increased diffusivity) and negatively correlated with the expression of synaptic and metabolic genes that have been implicated in neurodegeneration. These findings were consistent for diffusion MRI and volumetric HD-specific brain maps. As wild-type huntingtin is known to play a role in neurodevelopment, we explored the association between wild-type huntingtin (HTT) expression and developmental gene expression across the AHBA. Co-expression network analyses in 134 human brains free of neurodegenerative disorders were also performed. HTT expression was correlated with the expression of genes involved in neurodevelopment while co-expression network analyses also revealed that HTT expression was associated with developmental biological processes. Expression weighted cell-type enrichment (EWCE) analyses were used to explore which specific cell types were associated with HD cortical cell loss and these associations were validated using cell specific single nucleus RNAseq (snRNAseq) data from post-mortem HD brains. The developmental transcriptomic profile of cortical cell loss in preHD was enriched in astrocytes and endothelial cells, while the neurodegenerative transcriptomic profile was enriched for neuronal and microglial cells. Astrocyte-specific genes differentially expressed in HD post-mortem brains relative to controls using snRNAseq were enriched in the developmental transcriptomic profile, while neuronal and microglial-specific genes were enriched in the neurodegenerative transcriptomic profile. Our findings suggest that cortical cell loss in preHD may arise from dual pathological processes, emerging as a consequence of neurodevelopmental changes, at the beginning of life, followed by neurodegeneration in adulthood, targeting areas with reduced expression of synaptic and metabolic genes. These events result in age-related cell death across multiple brain cell types.

Neuroanatomical normative modelling can capture individual variability in Alzheimer's Disease (AD). We used neuroanatomical normative modelling to track individuals' disease progression in people with mild cognitive impairment (MCI) and patients with AD.

Identifying genetic risk factors underpinning different aspects of Alzheimer's disease has the potential to provide important insights into pathogenesis. Moving away from simple case-control definitions, there is considerable interest in using quantitative endophenotypes, such as those derived from imaging as outcome measures. Previous genome-wide association studies of imaging-derived biomarkers in sporadic late-onset Alzheimer's disease focused only on phenotypes derived from single imaging modalities. In contrast, we computed a novel multi-modal neuroimaging phenotype comprising cortical amyloid burden and bilateral hippocampal volume. Both imaging biomarkers were used as input to a disease progression modelling algorithm, which estimates the biomarkers' long-term evolution curves from population-based longitudinal data. Among other parameters, the algorithm computes the shift in time required to optimally align a subjects' biomarker trajectories with these population curves. This time shift serves as a disease progression score and it was used as a quantitative trait in a discovery genome-wide association study with n = 944 subjects from the Alzheimer's Disease Neuroimaging Initiative database diagnosed as Alzheimer's disease, mild cognitive impairment or healthy at the time of imaging. We identified a genome-wide significant locus implicating LCORL (rs6850306, chromosome 4; P = 1.03 × 10-8). The top variant rs6850306 was found to act as an expression quantitative trait locus for LCORL in brain tissue. The clinical role of rs6850306 in conversion from healthy ageing to mild cognitive impairment or Alzheimer's disease was further validated in an independent cohort comprising healthy, older subjects from the National Alzheimer's Coordinating Center database. Specifically, possession of a minor allele at rs6850306 was protective against conversion from mild cognitive impairment to Alzheimer's disease in the National Alzheimer's Coordinating Center cohort (hazard ratio = 0.593, 95% confidence interval = 0.387-0.907, n = 911, PBonf = 0.032), in keeping with the negative direction of effect reported in the genome-wide association study (βdisease progression score = -0.07 ± 0.01). The implicated locus is linked to genes with known connections to Alzheimer's disease pathophysiology and other neurodegenerative diseases. Using multimodal imaging phenotypes in association studies may assist in unveiling the genetic drivers of the onset and progression of complex diseases.

Optimizing HIV therapies is crucial since the virus rapidly develops mutations to evade drug pressure. Recent studies have shown that genotypic information might not be sufficient for the design of therapies and that other clinical and demographical factors may play a role in therapy failure. This study is designed to assess the improvement in prediction achieved when such information is taken into account. We use these factors to generate a prediction engine using a variety of machine learning methods and to determine which clinical conditions are most misleading in terms of predicting the outcome of a therapy.

Depression risk is exacerbated by genetic factors and stress exposure; however, the biological mechanisms through which these factors interact to confer depression risk are poorly understood. One putative biological mechanism implicates variability in the ability of cortisol, released in response to stress, to trigger a cascade of adaptive genomic and non-genomic processes through glucocorticoid receptor (GR) activation. Here, we demonstrate that common genetic variants in long-range enhancer elements modulate the immediate transcriptional response to GR activation in human blood cells. These functional genetic variants increase risk for depression and co-heritable psychiatric disorders. Moreover, these risk variants are associated with inappropriate amygdala reactivity, a transdiagnostic psychiatric endophenotype and an important stress hormone response trigger. Network modeling and animal experiments suggest that these genetic differences in GR-induced transcriptional activation may mediate the risk for depression and other psychiatric disorders by altering a network of functionally related stress-sensitive genes in blood and brain.

Incomplete hippocampal inversion (IHI), also called hippocampal malrotation, is an atypical presentation of the hippocampus present in about 20% of healthy individuals. Here we conducted the first genome-wide association study (GWAS) in IHI to elucidate the genetic underpinnings that may contribute to the incomplete inversion during brain development. A total of 1381 subjects contributed to the discovery cohort obtained from the IMAGEN database. The incidence rate of IHI was 26.1%. Loci with P<1e-5 were followed up in a validation cohort comprising 161 subjects from the PING study. Summary statistics from the discovery cohort were used to compute IHI heritability as well as genetic correlations with other traits. A locus on 18q11.2 (rs9952569; OR = 1.999; Z = 5.502; P = 3.755e-8) showed a significant association with the presence of IHI. A functional annotation of the locus implicated genes AQP4 and KCTD1. However, neither this locus nor the other 16 suggestive loci reached a significant p-value in the validation cohort. The h2 estimate was 0.54 (sd: 0.30) and was significant (Z = 1.8; P = 0.036). The top three genetic correlations of IHI were with traits representing either intelligence or education attainment and reached nominal P< = 0.013.

KLOTHO∗VS heterozygosity (KL∗VSHET+) was recently shown to be associated with reduced risk of Alzheimer's disease (AD) in APOE∗4 carriers. Additional studies suggest that KL∗VSHET+ protects against amyloid burden in cognitively normal older subjects, but sample sizes were too small to draw definitive conclusions. We performed a well-powered meta-analysis across 5 independent studies, comprising 3581 pre-clinical participants ages 60-80, to investigate whether KL∗VSHET+ reduces the risk of having an amyloid-positive positron emission tomography scan. Analyses were stratified by APOE∗4 status. KL∗VSHET+ reduced the risk of amyloid positivity in APOE∗4 carriers (odds ratio = 0.67 [0.52-0.88]; p = 3.5 × 10-3), but not in APOE∗4 non-carriers (odds ratio = 0.94 [0.73-1.21]; p = 0.63). The combination of APOE∗4 and KL∗VS genotypes should help enrich AD clinical trials for pre-symptomatic subjects at increased risk of developing amyloid aggregation and AD. KL-related pathways may help elucidate protective mechanisms against amyloid accumulation and merit exploration for novel AD drug targets. Future investigation of the biological mechanisms by which KL interacts with APOE∗4 and AD are warranted.

Hearing loss (HL) is one of the most common age-related diseases. Here, we investigate the central auditory correlates of HL in people with normal cognition and mild cognitive impairment (MCI) and test their association with genetic markers with the aim of revealing pathogenic mechanisms.

Microarray technologies are established approaches for high throughput gene expression, methylation and genotyping analysis. An accurate mapping of the array probes is essential to generate reliable biological findings. However, manufacturers of the microarray platforms typically provide incomplete and outdated annotation tables, which often rely on older genome and transcriptome versions that differ substantially from up-to-date sequence databases. Here, we present the Re-Annotator, a re-annotation pipeline for microarray probe sequences. It is primarily designed for gene expression microarrays but can also be adapted to other types of microarrays. The Re-Annotator uses a custom-built mRNA reference database to identify the positions of gene expression array probe sequences. We applied Re-Annotator to the Illumina Human-HT12 v4 microarray platform and found that about one quarter (25%) of the probes differed from the manufacturer's annotation. In further computational experiments on experimental gene expression data, we compared Re-Annotator to another probe re-annotation tool, ReMOAT, and found that Re-Annotator provided an improved re-annotation of microarray probes. A thorough re-annotation of probe information is crucial to any microarray analysis. The Re-Annotator pipeline is freely available at http://sourceforge.net/projects/reannotator along with re-annotated files for Illumina microarrays HumanHT-12 v3/v4 and MouseRef-8 v2.

Recent work has shown that people with common epilepsies have characteristic patterns of cortical thinning, and that these changes may be progressive over time. Leveraging a large multicenter cross-sectional cohort, we investigated whether regional morphometric changes occur in a sequential manner, and whether these changes in people with mesial temporal lobe epilepsy and hippocampal sclerosis (MTLE-HS) correlate with clinical features.

Amyloid-β is thought to facilitate the spread of tau throughout the neocortex in Alzheimer's disease, though how this occurs is not well understood. This is because of the spatial discordance between amyloid-β, which accumulates in the neocortex, and tau, which accumulates in the medial temporal lobe during ageing. There is evidence that in some cases amyloid-β-independent tau spreads beyond the medial temporal lobe where it may interact with neocortical amyloid-β. This suggests that there may be multiple distinct spatiotemporal subtypes of Alzheimer's-related protein aggregation, with potentially different demographic and genetic risk profiles. We investigated this hypothesis, applying data-driven disease progression subtyping models to post-mortem neuropathology and in vivo PET-based measures from two large observational studies: the Alzheimer's Disease Neuroimaging Initiative (ADNI) and the Religious Orders Study and Rush Memory and Aging Project (ROSMAP). We consistently identified 'amyloid-first' and 'tau-first' subtypes using cross-sectional information from both studies. In the amyloid-first subtype, extensive neocortical amyloid-β precedes the spread of tau beyond the medial temporal lobe, while in the tau-first subtype, mild tau accumulates in medial temporal and neocortical areas prior to interacting with amyloid-β. As expected, we found a higher prevalence of the amyloid-first subtype among apolipoprotein E (APOE) ε4 allele carriers while the tau-first subtype was more common among APOE ε4 non-carriers. Within tau-first APOE ε4 carriers, we found an increased rate of amyloid-β accumulation (via longitudinal amyloid PET), suggesting that this rare group may belong within the Alzheimer's disease continuum. We also found that tau-first APOE ε4 carriers had several fewer years of education than other groups, suggesting a role for modifiable risk factors in facilitating amyloid-β-independent tau. Tau-first APOE ε4 non-carriers, in contrast, recapitulated many of the features of primary age-related tauopathy. The rate of longitudinal amyloid-β and tau accumulation (both measured via PET) within this group did not differ from normal ageing, supporting the distinction of primary age-related tauopathy from Alzheimer's disease. We also found reduced longitudinal subtype consistency within tau-first APOE ε4 non-carriers, suggesting additional heterogeneity within this group. Our findings support the idea that amyloid-β and tau may begin as independent processes in spatially disconnected regions, with widespread neocortical tau resulting from the local interaction of amyloid-β and tau. The site of this interaction may be subtype-dependent: medial temporal lobe in amyloid-first, neocortex in tau-first. These insights into the dynamics of amyloid-β and tau may inform research and clinical trials that target these pathologies.

Progressive functional decline in the epilepsies is largely unexplained. We formed the ENIGMA-Epilepsy consortium to understand factors that influence brain measures in epilepsy, pooling data from 24 research centres in 14 countries across Europe, North and South America, Asia, and Australia. Structural brain measures were extracted from MRI brain scans across 2149 individuals with epilepsy, divided into four epilepsy subgroups including idiopathic generalized epilepsies (n =367), mesial temporal lobe epilepsies with hippocampal sclerosis (MTLE; left, n = 415; right, n = 339), and all other epilepsies in aggregate (n = 1026), and compared to 1727 matched healthy controls. We ranked brain structures in order of greatest differences between patients and controls, by meta-analysing effect sizes across 16 subcortical and 68 cortical brain regions. We also tested effects of duration of disease, age at onset, and age-by-diagnosis interactions on structural measures. We observed widespread patterns of altered subcortical volume and reduced cortical grey matter thickness. Compared to controls, all epilepsy groups showed lower volume in the right thalamus (Cohen's d = -0.24 to -0.73; P < 1.49 × 10-4), and lower thickness in the precentral gyri bilaterally (d = -0.34 to -0.52; P < 4.31 × 10-6). Both MTLE subgroups showed profound volume reduction in the ipsilateral hippocampus (d = -1.73 to -1.91, P < 1.4 × 10-19), and lower thickness in extrahippocampal cortical regions, including the precentral and paracentral gyri, compared to controls (d = -0.36 to -0.52; P < 1.49 × 10-4). Thickness differences of the ipsilateral temporopolar, parahippocampal, entorhinal, and fusiform gyri, contralateral pars triangularis, and bilateral precuneus, superior frontal and caudal middle frontal gyri were observed in left, but not right, MTLE (d = -0.29 to -0.54; P < 1.49 × 10-4). Contrastingly, thickness differences of the ipsilateral pars opercularis, and contralateral transverse temporal gyrus, were observed in right, but not left, MTLE (d = -0.27 to -0.51; P < 1.49 × 10-4). Lower subcortical volume and cortical thickness associated with a longer duration of epilepsy in the all-epilepsies, all-other-epilepsies, and right MTLE groups (beta, b < -0.0018; P < 1.49 × 10-4). In the largest neuroimaging study of epilepsy to date, we provide information on the common epilepsies that could not be realistically acquired in any other way. Our study provides a robust ranking of brain measures that can be further targeted for study in genetic and neuropathological studies. This worldwide initiative identifies patterns of shared grey matter reduction across epilepsy syndromes, and distinctive abnormalities between epilepsy syndromes, which inform our understanding of epilepsy as a network disorder, and indicate that certain epilepsy syndromes involve more widespread structural compromise than previously assumed.

Medical artificial intelligence (AI) offers great potential for recognizing signs of health conditions in retinal images and expediting the diagnosis of eye diseases and systemic disorders1. However, the development of AI models requires substantial annotation and models are usually task-specific with limited generalizability to different clinical applications2. Here, we present RETFound, a foundation model for retinal images that learns generalizable representations from unlabelled retinal images and provides a basis for label-efficient model adaptation in several applications. Specifically, RETFound is trained on 1.6 million unlabelled retinal images by means of self-supervised learning and then adapted to disease detection tasks with explicit labels. We show that adapted RETFound consistently outperforms several comparison models in the diagnosis and prognosis of sight-threatening eye diseases, as well as incident prediction of complex systemic disorders such as heart failure and myocardial infarction with fewer labelled data. RETFound provides a generalizable solution to improve model performance and alleviate the annotation workload of experts to enable broad clinical AI applications from retinal imaging.

The brain's functional connectome is dynamic, constantly reconfiguring in an individual-specific manner. However, which characteristics of such reconfigurations are subject to genetic effects, and to what extent, is largely unknown. Here, we identified heritable dynamic features, quantified their heritability, and determined their association with cognitive phenotypes. In resting-state fMRI, we obtained multivariate features, each describing a temporal or spatial characteristic of connectome dynamics jointly over a set of connectome states. We found strong evidence for heritability of temporal features, particularly, Fractional Occupancy (FO) and Transition Probability (TP), representing the duration spent in each connectivity configuration and the frequency of shifting between configurations, respectively. These effects were robust against methodological choices of number of states and global signal regression. Genetic effects explained a substantial proportion of phenotypic variance of these features (h2=0.39, 95% CI= [.24,.54] for FO; h2=0.43, 95% CI=[.29,.57] for TP). Moreover, these temporal phenotypes were associated with cognitive performance. Contrarily, we found no robust evidence for heritability of spatial features of the dynamic states (i.e., states' Modularity and connectivity pattern). Genetic effects may therefore primarily contribute to how the connectome transitions across states, rather than the precise spatial instantiation of the states in individuals. In sum, genetic effects impact the dynamic trajectory of state transitions (captured by FO and TP), and such temporal features may act as endophenotypes for cognitive abilities.

Nomograms are important clinical tools applied widely in both developing and aging populations. They are generally constructed as normative models identifying cases as outliers to a distribution of healthy controls. Currently used normative models do not account for genetic heterogeneity. Hippocampal volume (HV) is a key endophenotype for many brain disorders. Here, we examine the impact of genetic adjustment on HV nomograms and the translational ability to detect dementia patients. Using imaging data from 35,686 healthy subjects aged 44-82 from the UK Biobank (UKB), we built HV nomograms using Gaussian process regression (GPR), which - compared to a previous method - extended the application age by 20 years, including dementia critical age ranges. Using HV polygenic scores (HV-PGS), we built genetically adjusted nomograms from participants stratified into the top and bottom 30% of HV-PGS. This shifted the nomograms in the expected directions by ~100 mm3 (2.3% of the average HV), which equates to 3 years of normal aging for a person aged ~65. Clinical impact of genetically adjusted nomograms was investigated by comparing 818 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database diagnosed as either cognitively normal (CN), having mild cognitive impairment (MCI) or Alzheimer's disease (AD) patients. While no significant change in the survival analysis was found for MCI-to-AD conversion, an average of 68% relative decrease was found in intra-diagnostic-group variance, highlighting the importance of genetic adjustment in untangling phenotypic heterogeneity.

State-of-the-art rare variant association testing methods aggregate the contribution of rare variants in biologically relevant genomic regions to boost statistical power. However, testing single genes separately does not consider the complex interaction landscape of genes, nor the downstream effects of non-synonymous variants on protein structure and function. Here we present the NETwork Propagation-based Assessment of Genetic Events (NETPAGE), an integrative approach aimed at investigating the biological pathways through which rare variation results in complex disease phenotypes. We applied NETPAGE to sporadic, late-onset Alzheimer's disease (AD), using whole-genome sequencing from the AD Neuroimaging Initiative (ADNI) cohort, as well as whole-exome sequencing from the AD Sequencing Project (ADSP). NETPAGE is based on network propagation, a framework that models information flow on a graph and simulates the percolation of genetic variation through tissue-specific gene interaction networks. The result of network propagation is a set of smoothed gene scores that can be tested for association with disease status through sparse regression. The application of NETPAGE to AD enabled the identification of a set of connected genes whose smoothed variation profile was robustly associated to case-control status, based on gene interactions in the hippocampus. Additionally, smoothed scores significantly correlated with risk of conversion to AD in Mild Cognitive Impairment (MCI) subjects. Lastly, we investigated tissue-specific transcriptional dysregulation of the core genes in two independent RNA-seq datasets, as well as significant enrichments in terms of gene sets with known connections to AD. We present a framework that enables enhanced genetic association testing for a wide range of traits, diseases, and sample sizes.

Age-related hearing loss was recently established as the largest modifiable risk factor for Alzheimer's disease (AD), however, the reasons for this link remain unclear. We investigate shared underlying genetic associations using results from recent large genome-wide association studies (GWAS) on adult hearing difficulty and AD. Genetic correlation and Mendelian randomization (MR) analysis do not support a genetic correlation between the disorders, but suggest a direct causal link from AD genetic risk to hearing difficulty, driven by APOE. Systematic MR analyses on the effect of other traits revealed shared effects of glutamine, gamma-glutamylglutamine, and citrate levels on reduced risk of both hearing difficulty and AD. In addition, pathway analysis on GWAS risk variants suggests shared function in neuronal signalling pathways as well as etiology of diabetes and cardiovascular disease. However, after multiple testing corrections, neither analysis led to statistically significant associations. Altogether, our genetic-driven analysis suggests hearing difficulty and AD are linked by a shared vulnerability in molecular pathways rather than by a shared genetic architecture.

Frontotemporal dementia is a heterogeneous neurodegenerative disorder characterized by neuronal loss in the frontal and temporal lobes. Despite progress in understanding which genes are associated with the aetiology of frontotemporal dementia, the biological basis of how mutations in these genes lead to cell loss in specific cortical regions remains unclear. In this work we combined gene expression data for 16,772 genes from the Allen Institute for Brain Science atlas with brain maps of gray matter atrophy in symptomatic C9orf72, GRN and MAPT mutation carriers obtained from the Genetic Frontotemporal dementia Initiative study. No significant association was seen between C9orf72, GRN and MAPT expression and the atrophy patterns in the respective genetic groups. After adjusting for spatial autocorrelation, between 1,000 and 5,000 genes showed a negative or positive association with the atrophy pattern within each individual genetic group, with the most significantly associated genes being TREM2, SSBP3 and GPR158 (negative association in C9orf72, GRN and MAPT respectively) and RELN, MXRA8 and LPA (positive association in C9orf72, GRN and MAPT respectively). An overrepresentation analysis identified a negative association with genes involved in mitochondrial function, and a positive association with genes involved in vascular and glial cell function in each of the genetic groups. A set of 423 and 700 genes showed significant positive and negative association, respectively, with atrophy patterns in all three maps. The gene set with increased expression in spared cortical regions was enriched for neuronal and microglial genes, while the gene set with increased expression in atrophied regions was enriched for astrocyte and endothelial cell genes. Our analysis suggests that these cell types may play a more active role in the onset of neurodegeneration in frontotemporal dementia than previously assumed, and in the case of the positively-associated cell marker genes, potentially through emergence of neurotoxic astrocytes and alteration in the blood-brain barrier respectively.