Search

BackgroundMyeloid cells, including monocytes, macrophages, microglia, dendritic cells and neutrophils are a part of innate immunity, playing a major role in orchestrating innate and adaptive immune responses. Microglia are the resident myeloid cells of the central nervous system, and many Alzheimer’s disease (AD) risk loci are found in or near genes that are highly or sometimes uniquely expressed in myeloid cells. Similarly, inflammatory bowel disease (IBD) loci are also enriched for genes expressed by myeloid cells. However, the extent to which there is overlap between the effects of AD and IBD susceptibility loci in myeloid cells remains poorly described, and the substantial IBD genetic maps may help to accelerate AD research.MethodsHere, we leveraged summary statistics from large-scale genome-wide association studies (GWAS) to investigate the causal effect of IBD (including ulcerative colitis and Crohn’s disease) variants on AD and AD endophenotypes. Microglia and monocyte expression Quantitative Trait Locus (eQTLs) were used to examine the functional consequences of IBD and AD risk variants enrichment in two different myeloid cell subtypes.ResultsOur results showed that, whilePTK2Bis implicated in both diseases and both sets of risk loci are enriched for myeloid genes, AD and IBD susceptibility loci largely implicate distinct sets of genes and pathways. AD loci are significantly more enriched for microglial eQTLs than IBD. We also found that genetically determined IBD is associated with a lower risk of AD, which may driven by a negative effect on the accumulation of neurofibrillary tangles (beta=−1.04, p=0.013). In addition, IBD displayed a significant positive genetic correlation with psychiatric disorders and multiple sclerosis, while AD showed a significant positive genetic correlation with amyotrophic lateral sclerosis.ConclusionTo our knowledge, this is the first study to systematically contrast the genetic association between IBD and AD, our findings highlight a possible genetically protective effect of IBD on AD even if the majority of effects on myeloid cell gene expression by the two sets of disease variants are distinct. Thus, IBD myeloid studies may not help to accelerate AD functional studies, but our observation reinforces the role of myeloid cells in the accumulation of tau proteinopathy and provides a new avenue for discovering a protective factor.

Purpose: Cancer cell phenotypes are partially determined by epigenetic specifications, such as DNA methylation. Metastasis development is a late event in cancerogenesis and might be associated with epigenetic alterations.Experimental Design: An in vivo selection approach was used to generate highly aggressive non–small cell lung cancer (NSCLC) cell lines (A549 and HTB56) followed by genome-wide DNA methylation analysis. Furthermore, the therapeutic effects of the epigenetic agent azacytidine on DNA methylation patterns and the in vivo phenotypes were explored.Results: Widespread changes of DNA methylation were observed during development of highly aggressive cell lines. Up to 2.5% of the CpG-rich region was differentially methylated as identified by reduced representation bisulfite sequencing compared with the less aggressive parental cell lines. DNA methyltransferase inhibition by azacytidine reversed the prometastatic phenotype; this was highly associated with the preferential loss of DNA methylation at sites that were hypermethylated during the in vivo selection. Of note, polycomb (PRC2) binding sites were particularly affected by DNA methylation changes after azacytidine exposure that persisted over time.Conclusions: We could show that metastatic capability of NSCLC is closely associated with DNA methylome alterations. Because inhibition of DNA methyltransferase reversed metastasis-prone phenotype, epigenetic modulation seems to be a potential therapeutic approach to prevent metastasis formation. Clin Cancer Res; 20(4); 814–26. ©2013 AACR.

PDF file - 52K

XLSX file - 655K, Supplementary Table S4 - annotated DNA methylation in confirmed promoter DMRs (RRBS) in HTB56.

Supplementary Figures S1-6. Supplementary Fig S1: Combined patient characteristics and subjected methods Supplementary Fig. S2 A: ROC curves were constructed by using the somatic scores from the sequenced genotypes. Supplementary Fig. S3: Circos plots of representative Ewing sarcomas. Supplementary Fig. S4: Sanger sequencing results confirmed (A) somatic inactivating mutations of the STAG2 gene (positions according to RefSeq NG 033796.2) and (B) the activating N546 mutation of FGFR1. Supplementary Fig. S5: FGFR-1 expression in Ewing Sarcomas by means of real time RT PCR (A, normalized to FGFR-1 expression in VH-64). (B) Array-based analysis of FGFR1 expression in ES compared to normal tissue (NBA; SI Materials and Methods). (C) Knock-down by means of lentiviral transduced FGFR1.shRNA effectively reduced FGFR1 mRNA expression by 70 % or more in Ewing sarcoma cell lines as compared to scrambled shRNA. (D) Knockdown was confirmed on the protein level by means of western blotting. Supplementary Fig. S6: Amplification of the FGFR-1 gene (A) as well as gene expression (B) are associated with a trend towards inferior survival in patients with ES. (C) Treatment with Ponatinib inhibited growth of ES cell lines significantly. (D) FGFR1 tyrosine kinase inhibitor therapy in a patient with relapsed Ewing sarcoma.

Supplementary Materials. Additional Material and Methods section and clinical data of the Ponatib treated patient.

Purpose: A low mutation rate seems to be a general feature of pediatric cancers, in particular in oncofusion gene-driven tumors. Genetically, Ewing sarcoma is defined by balanced chromosomal EWS/ETS translocations, which give rise to oncogenic chimeric proteins (EWS-ETS). Other contributing somatic mutations involved in disease development have only been observed at low frequency.Experimental Design: Tumor samples of 116 Ewing sarcoma patients were analyzed here. Whole-genome sequencing was performed on two patients with normal, primary, and relapsed tissue. Whole-exome sequencing was performed on 50 Ewing sarcoma and 22 matched normal tissues. A discovery dataset of 14 of these tumor/normal pairs identified 232 somatic mutations. Recurrent nonsynonymous mutations were validated in the 36 remaining exomes. Transcriptome analysis was performed in a subset of 14 of 50 Ewing sarcomas and DNA copy number gain and expression of FGFR1 in 63 of 116 Ewing sarcomas.Results: Relapsed tumors consistently showed a 2- to 3-fold increased number of mutations. We identified several recurrently mutated genes at low frequency (ANKRD30A, CCDC19, KIAA0319, KIAA1522, LAMB4, SLFN11, STAG2, TP53, UNC80, ZNF98). An oncogenic fibroblast growth factor receptor 1 (FGFR1) mutation (N546K) was detected, and the FGFR1 locus frequently showed copy number gain (31.7%) in primary tumors. Furthermore, high-level FGFR1 expression was noted as a characteristic feature of Ewing sarcoma. RNA interference of FGFR1 expression in Ewing sarcoma lines blocked proliferation and completely suppressed xenograft tumor growth. FGFR1 tyrosine kinase inhibitor (TKI) therapy in a patient with Ewing sarcoma relapse significantly reduced 18-FDG-PET activity.Conclusions: FGFR1 may constitute a promising target for novel therapeutic approaches in Ewing sarcoma. Clin Cancer Res; 21(21); 4935–46. ©2015 AACR.

PDF file - 118K

PDF file - 1572K, Figure S4 4A+B): Column charts depict hypo- and hypermethylated regions detected by RRBS in high metastatic A549 (A) and HTB56 (B) cells after 6 days of 5-Azacytidine-treatment and after 7 days of 5-Azacytidine-release. 4C+D) Smoothed scatter plot of 5-Azacytidine-treated (250nM) versus untreated high metastatic HTB56 cells. Shown are methylation levels for CpG-sites analyzed by RRBS. C: X-axis: Methylation levels in HTB56 high metastatic cells (d0=untreated) Y-axis: Methylation levels in HTB56 high metastatic cells (d6_250= after six days of 5-Azacytidinetreatment, 5-Azacytidine-concentration of 250 nM) D: X-axis: Methylation levels in HTB56 high metastatic cells (d0=untreated) Y-axis: Methylation levels in HTB56 high metastatic cells (d13_250= after thirteen days of 5- Azacytidine-treatment, 5-Azacytidine-concentration of 250 nM) Figure S5 Unsupervised hierarchical clustering shows that 5-Azacytidine exposed cells cluster more closely together than individual cell lines. Figure S6 Chromosomal distribution of clusters and DMRs. The curves show the density distribution of hypomethylated DMRs and CpG clusters in relation to their chromosomal position: 0 indicates chromosomal end and 1 the centromeric region. Hypomethylated DMRs (green), hypermethylated DMRs (red) and CpG clusters (black) are shown. DNA demethylation upon 5-Azacytidine occurs primarily on chromosome ends.

PDF file - 247K, Supplementary Figure 3 Unsupervised hierarchical clustering shows that independent samples at each step of selection cluster more closely together than individual cell lines A549 (A) and HTB56 (B). A high correlation between biological replicates from normal and highly aggressive cell lines from both cell types were observed. The clustering is based on a distance matrix, which has been derived from pearson product-moment correlation coefficients (r) by calculating 1-r.

PDF file - 219K, Supplementary Figure 2 Chart depicts viability of high and low metastatic cells before 5-Azacytidine-treatment and after 7 days of 5-Azacytidine-release.

PDF file - 233K, Supplementary Figure 1 Genomic changes in low and high metastatic A549 cells were analyzed by 500K SNP arrays.

Supplementary Tables S1-6. Table S1: Validation of somatic mutations by Sanger sequencing for construction of ROC Table S2: Characteristics of the discovery data set for exome sequencing Table S3: Somatic mutations in the discovery data set (filtered with MutSigCV1.6) Table S4: Recurrent non-synonymous mutations (filtered with MutSigCV1.6) Table S5: Characteristics of the expression data set. Table S6 : Shared mutations in consecutive sarcomas

PDF file - 96K, Supplementary Table 1: Mutations detected in parental and highly metastatic NSCLC cell lines A549 (A) and HTB56 (B) by SNP analysis Supplementary Table 2: Mutations detected in parental and highly metastatic NSCLC cell lines A549 and HTB56 by Exome sequencing Supplementary Table 3: Conversion rate uniquely mapped reads

BackgroundDepression is a common psychiatric illness and global public health problem. However, our limited understanding of the biological basis of depression has hindered the development of novel treatments and interventions.MethodsTo identify new candidate genes for therapeutic development, we examined single-nucleus RNA sequencing (snucRNAseq) data from the dorsolateral prefrontal cortex (N=424) in relation to ante-mortem depressive symptoms. To complement these direct analyses, we also used genome- wide association study (GWAS) results for depression (N=500,199) along with genetic tools for inferring the expression of 22,159 genes in 7 cell types and 55 cell subtypes to perform transcriptome-wide association studies (TWAS) of depression followed by Mendelian randomization (MR).ResultsOur single-nucleus TWAS analysis identified 71 causal genes in depression that have a role in specific neocortical cell subtypes; 59 of 71 genes were novel compared to previous studies. Depression TWAS genes showed a cell type specific pattern, with the greatest enrichment being in both excitatory and inhibitory neurons as well as astrocytes. Gene expression in different neuron subtypes have different directions of effect on depression risk. Compared to lower genetically correlated traits (e.g. body mass index) with depression, higher correlated traits (e.g., neuroticism) have more common TWAS genes with depression. In parallel, we performed differential gene expression analysis in relation to depression in 55 cortical cell subtypes, and we found that genes such asANKRD36,MADD,TAOK3,SCAIandCHUKare associated with depression in specific cell subtypes.ConclusionsThese two sets of analyses illustrate the utility of large snucRNAseq data to uncover both genes whose expression is altered in specific cell subtypes in the context of depression and to enhance the interpretation of well-powered GWAS so that we can prioritize specific susceptibility genes for further analysis and therapeutic development.

ZCCHC17 is a putative master regulator of synaptic gene dysfunction in Alzheimer’s Disease (AD), and ZCCHC17 protein declines early in AD brain tissue, before significant gliosis or neuronal loss. Here, we investigate the function of ZCCHC17 and its role in AD pathogenesis. Co-immunoprecipitation of ZCCHC17 followed by mass spectrometry analysis in human iPSC-derived neurons reveals that ZCCHC17’s binding partners are enriched for RNA splicing proteins. ZCCHC17 knockdown results in widespread RNA splicing changes that significantly overlap with splicing changes found in AD brain tissue, with synaptic genes commonly affected. ZCCHC17 expression correlates with cognitive resilience in AD patients, and we uncover an APOE4 dependent negative correlation of ZCCHC17 expression with tangle burden. Furthermore, a majority of ZCCHC17 interactors also co-IP with known tau interactors, and we find significant overlap between alternatively spliced genes in ZCCHC17 knockdown and tau overexpression neurons. These results demonstrate ZCCHC17’s role in neuronal RNA processing and its interaction with pathology and cognitive resilience in AD, and suggest that maintenance of ZCCHC17 function may be a therapeutic strategy for preserving cognitive function in the setting of AD pathology.SignificanceAbnormal RNA processing is an important component of AD pathophysiology. We show here that ZCCHC17, a previously identified putative master regulator of synaptic dysfunction in AD, plays a role in neuronal RNA processing, and illustrate that ZCCHC17 dysfunction is sufficient to explain some of the splicing abnormalities seen in AD brain tissue, including synaptic gene splicing abnormalities. Using data from human patients, we demonstrate that ZCCHC17 mRNA levels correlate with cognitive resilience in the setting of AD pathology. These results suggest that maintenance of ZCCHC17 function may be a therapeutic strategy for supporting cognitive function in AD patients, and motivate future work examining a possible role of abnormal RNA processing in AD-related cognitive decline.

Background Motor resilience proteins have not been identified. This proteome-wide discovery study sought to identify proteins that may provide motor resilience. Methods We studied the brains of older decedents with annual motor testing, postmortem brain pathologies, and proteome-wide data. Parkinsonism was assessed using 26 items of a modified United Parkinson Disease Rating Scale. We used linear mixed-effect models to isolate motor resilience, defined as the person-specific estimate of progressive parkinsonism after controlling for age, sex, and 10 brain pathologies. A total of 8 356 high-abundance proteins were quantified from dorsal lateral prefrontal cortex using tandem mass tag and liquid chromatography–mass spectrometry. Results There were 391 older adults (70% female), mean age 80 years at baseline and 89 years at death. Five proteins were associated with motor resilience: A higher level of AP1B1 (Estimate −0.504, SE 0.121, p = 3.12 × 10−5) and GNG3 (Estimate −0.276, SE 0.068, p = 4.82 × 10−5) was associated with slower progressive parkinsonism. By contrast, a higher level of TTC38 (Estimate 0.140, SE 0.029, p = 1.87 × 10−6), CARKD (Estimate 0.413, SE 0.100, p = 3.50 × 10−5), and ABHD14B (Estimate 0.175, SE 0.044, p = 6.48 × 10−5) was associated with faster progressive parkinsonism. Together, these 5 proteins accounted for almost 25% of the variance of progressive parkinsonism above the 17% accounted for by 10 indices of brain pathologies. Discussion Cortical proteins may provide more or less motor resilience in older adults. These proteins are high-value therapeutic targets for drug discovery that may lead to interventions that maintain motor function despite the accumulation of as yet untreatable brain pathologies.

Background and ObjectivesCognitive resilience is a well-recognized concept, but knowledge gaps about its underlying mechanisms have made it difficult to develop instruments that identify older adults with high or low resilience. We tested whether aggregating cortical peptides associated with cognitive resilience into an index can identify adults with higher or lower cognitive resilience.MethodsWe used data from 1,192 older decedents, including annual clinical testing, indices of 10 Alzheimer disease (AD) and related dementia (ADRD) pathologies, and 226 proteotypic peptides measured in the dorsal lateral prefrontal cortex. We used linear mixed-effects models to identify peptides that were related to cognitive resilience (i.e., cognitive decline not explained by ADRD pathologies [false discovery rate ResultsWe constructed a resilience index from 52 of 226 peptides related to cognitive resilience. A higher index was associated with slower cognitive decline (estimate 0.05, SE 0.003, p < 0.001) and slower motor decline (estimate 0.005, SE 0.001, p < 0.001). Most resilience peptides (70%) were specific to cognitive decline, but 30% also provided resilience for motor decline. A higher index was also related to a lower burden of AD pathologies (odds ratio [OR] 0.41, SE 0.01, p < 0.001) and modified the association of AD pathology with cognition in that a higher index modified the negative effects of AD pathology on AD dementia proximate to death (OR 0.70, SE 0.14, p = 0.010). Up to 90% of cognitive resilience peptides were related to AD pathologic phenotypes.DiscussionCortical proteins may provide some degree of cognitive resilience. These multifunctional proteins also seem to provide resilience to other AD clinical phenotypes and have independent associations with ADRD pathologies. Resilience proteins may be high-value therapeutic targets for drug discovery of interventions that maintain brain health in aging adults via multiple pathways.

SummarySORL1 is strongly implicated in the pathogenesis of Alzheimer’s disease (AD) through human genetic studies that point to an association of reduced SORL1 levels with higher risk for AD. To interrogate the role(s) of SORL1 in human brain cells, SORL1 null iPSCs were generated, followed by differentiation to neuron, astrocyte, microglia, and endothelial cell fates. Loss of SORL1 led to alterations in both overlapping and distinct pathways across cell types, with the greatest effects in neurons and astrocytes. Intriguingly, SORL1 loss led to a dramatic neuron-specific reduction in APOE levels. Further, analyses of iPSCs derived from a human aging cohort revealed a neuron-specific linear correlation between SORL1 and APOE RNA and protein levels, a finding validated in human post-mortem brain. Pathway analysis implicated intracellular transport pathways and TGF- β/SMAD signaling in the function of SORL1 in neurons. In accord, enhancement of retromer-mediated trafficking and autophagy rescued elevated phospho-tau observed in SORL1 null neurons but did not rescue APOE levels, suggesting that these phenotypes are separable. Stimulation and inhibition of SMAD signaling modulated APOE RNA levels in a SORL1-dependent manner. These studies provide a mechanistic link between two of the strongest genetic risk factors for AD.

The transfer of mitochondrial DNA into the nuclear genomes of eukaryotes (Numts) has been linked to lifespan in non-human species1–3and recently demonstrated to occur in rare instances from one human generation to the next4. Here we investigated numtogenesis dynamics in humans in two ways. First, we quantified Numts in 1,187 post-mortem brain and blood samples from different individuals. Compared to circulating immune cells (n=389), post-mitotic brain tissue (n=798) contained more Numts, consistent with their potential somatic accumulation. Within brain samples we observed a 5.5-fold enrichment of somatic Numt insertions in the dorsolateral prefrontal cortex compared to cerebellum samples, suggesting that brain Numts arose spontaneously during development or across the lifespan. Moreover, more brain Numts was linked to earlier mortality. The brains of individuals with no cognitive impairment who died at younger ages carried approximately 2 more Numts per decade of life lost than those who lived longer. Second, we tested the dynamic transfer of Numts using a repeated-measures WGS design in a human fibroblast model that recapitulates several molecular hallmarks of aging5. These longitudinal experiments revealed a gradual accumulation of one Numt every ∼13 days. Numtogenesis was independent of large-scale genomic instability and unlikely driven cell clonality. Targeted pharmacological perturbations including chronic glucocorticoid signaling or impairing mitochondrial oxidative phosphorylation (OxPhos) only modestly increased the rate of numtogenesis, whereas patient-derivedSURF1-mutant cells exhibiting mtDNA instability accumulated Numts 4.7- fold faster than healthy donors. Combined, our data document spontaneous numtogenesis in human cells and demonstrate an association between brain cortical somatic Numts and human lifespan. These findings open the possibility that mito-nuclear horizontal gene transfer among human post-mitotic tissues produce functionally-relevant human Numts over timescales shorter than previously assumed.

MotivationDroplet-based single-cell RNA sequencing (scRNA-seq) is widely used in biomedical research to interrogate the transcriptomes of single cells on a large scale. Pooling and processing cells from different samples together can reduce costs and batch effects. In order to pool cells, cells are often first labeled with hashtag oligonucleotides (HTOs). These HTOs are sequenced along with the cells’ RNA in the droplets and are subsequently used to computationally assign each droplet to its sample of origin, which is referred to as demultiplexing. Accurate demultiplexing is crucial and can be challenging due to background HTOs, low-quality cells/cell debris, and multiplets.ResultsA new demultiplexing method, demuxmix, based on negative binomial regression mixture models is introduced. The method implements two significant improvements. First, demuxmix’s probabilistic classification framework provides error probabilities for droplet assignments that can be used to discard uncertain droplets and inform about the quality of the HTO data and the demultiplexing success. Second, demuxmix utilizes the positive association between detected genes in the RNA library and HTO counts to explain parts of the variance in the HTO data resulting in improved droplet assignments. The improved performance of demuxmix compared to existing demultiplexing methods is assessed on real and simulated data. Finally, the feasibility of accurately demultiplexing experimental designs where non-labeled cells are pooled with labeled cells is demonstrated.AvailabilityR/Bioconductor package demuxmix (https://doi.org/doi:10.18129/B9.bioc.demuxmix)

The heterogeneity of the older population suggests the existence of subsets of individuals which share certain brain molecular features and respond differently to risk factors for Alzheimer’s disease, but this population structure remains poorly defined. Here, we performed an unsupervised clustering of individuals with multi-region brain transcriptomes to assess whether a broader approach, simultaneously considering data from multiple regions involved in cognition would uncover such subsets. We implemented a canonical correlation-based analysis in a Discovery cohort of 459 participants from two longitudinal studies of cognitive aging that have RNA sequence profiles in three brain regions. 690 additional participants that have data in only one or two of these regions were used in the Replication effort. These clustering analyses identified two meta-clusters, MC-1 and MC-2. The two sets of participants differ primarily in their trajectories of cognitive decline, with MC-2 having a delay of 3 years to the median age of incident dementia. This is due, in part, to a greater impact of tau pathology on neuronal chromatin architecture and to broader brain changes including greater loss of white matter integrity in MC-1. Further evidence of biological differences includes a significantly larger impact ofAPOEε4risk on cognitive decline in MC-1. These findings suggest that our proposed population structure captures an aspect of the more distributed molecular state of the aging brain that either enhances the effect of risk factors in MC-1 or of protective effects in MC-2. These observations may inform the design of therapeutic development efforts and of trials as both become increasingly more targeted molecularly.One Sentence Summary:There are two types of aging brains, with one being more vulnerable toAPOEε4and subsequent neuronal dysfunction and cognitive loss.

INTRODUCTIONNot all APOE ε4 carriers who survive to advanced age develop Alzheimer’s disease (AD); factors attenuating the risk of ε4 on AD may exist.METHODSGuided by the top ε4-attenuating signals from methylome-wide association analyses (N=572, ε4+ and ε4-) of neurofibrillary tangles and neuritic plaques, we conducted a meta-analysis for pathological AD within the ε4+ subgroups (N=235) across four independent collections of brains. Cortical RNA-seq and microglial morphology measurements were used in functional analyses.RESULTSThree out of the four significant CpG dinucleotides were captured by one principle component (PC1), which interacts with ε4 on AD, and is associated with expression of innate immune genes and activated microglia. In ε4 carriers, reduction in each unit of PC1 attenuated the odds of AD by 58% (OR=2.39, 95%CI=[1.64,3.46], P=7.08×10−6).DISCUSSIONAn epigenomic factor associated with a reduced proportion of activated microglia appears to attenuate the risk of ε4 on AD.

The relationship between genetic variation and gene expression in individual brain cell types and subtypes has remained elusive. Here, we generated single-nucleus RNA sequencing data from the dorsolateral prefrontal cortex of 424 individuals of advanced age; analyzing 1.5 million nuclear transcriptomes, we assessed the effect of genetic variants on RNA expression incis(cis-eQTL) for 7 cell types and 81 cell subtypes. This effort identified 10,004 eGenes at the cell type level and 8,138 eGenes at the cell subtype level. Many eGenes are only detected within cell subtypes. A new variant influencesAPOEexpression only in microglia and is associated with greater cerebral amyloid angiopathy but not Alzheimer pathology, accounting for the effect ofAPOEε4, providing mechanistic insights into both pathologies. While eQTLs are readily detected, only aTMEM106Bvariant robustly affects the proportion of cell subtypes. Integration of these results with GWAS highlighted the targeted cell type and likely causal gene within susceptibility loci for Alzheimer’s, Parkinson’s, schizophrenia, and educational attainment.

The polygenic and multi-cellular nature of multiple sclerosis (MS) immunopathology necessitates cell-type-specific molecular studies in order to improve our understanding of the diverse mechanisms underlying immune cell dysfunction in MS. Here, by generating a dataset of 1,075 transcriptomes from 209 participants (167 MS and 42 healthy), we assessed MS-associated transcriptional changes in six implicated cell-type-states: naïve and memory helper T cells and classical monocytes purified from peripheral blood, each in their primary (ex vivo, unstimulated) andin vitrostimulated states. Our data suggest that primary profiles show larger MS-associated differences than the post-stimulation contexts. We further identified shared and distinct changes in individual genes, biological pathways, and co-expressed gene modules in MS T cells and monocytes, and prioritized genes such asZBTB16as MS-associated regulators in both cell types. Of six identified MS-associated co-expressed gene modules, three (two lymphoid and one myeloid) were replicated in independent data from peripheral blood mononuclear cells (PBMC) and monocyte-derived macrophages. A subsequentin silicodrug screen prioritized small-molecule compounds for reversing the perturbation of the MS-associated modules. The effects of glucocorticoid receptor agonists as the top-identified therapeutic class for the replicated T cell modules were validated using targetedin silicoanalyses andin vitroexperiments, suggesting the coordinated dysregulation of glucocorticoid-responsive genes in MS T cells. In summary, our study identifies and validates individual genes and co-expressed gene modules from T and myeloid cells that are perturbed in MS, offering new targets for therapeutic discovery and biomarker development to guide the management of MS.

Human microglia play a pivotal role in neurological diseases, but few targeted therapies that directly modulate microglial state or function exist due to an incomplete understanding of microglial heterogeneity. We use single-cell RNA sequencing to profile live human microglia from autopsies or surgical resections across diverse neurological diseases and central nervous system regions. We observe a central divide between oxidative and heterocyclic metabolism and identify subsets associated with antigen presentation, motility, and proliferation. Specific subsets are enriched in susceptibility genes for neurodegenerative diseases or the disease-associated microglial signature. We validate subtypes in situ with an RNAscope-immunofluorescence pipeline and leverage our dataset as a classification resource, finding that iPSC model systems recapitulate substantial in vivo heterogeneity. Finally, we identify and validate candidates for chemically inducing subtype-specific states in vitro, showing that Camptothecin downregulates the transcriptional signature of disease-enriched subsets and upregulates a signature previously shown to be depleted in Alzheimer’s.

Despite a growing focus on neuroimmune mechanisms of Alzheimer’s disease (AD), the role of peripheral monocytes remains largely unknown. Circulating monocytes communicate with the brain’s resident myeloid cells, microglia, via chemical signaling and can directly infiltrate the brain parenchyma.1 Thus, molecular signatures of monocytes may serve as indicators of neuropathological events unfolding in the CNS.2–5 However, no studies have yet directly tested the association of monocyte gene expression on longitudinal cognitive decline or postmortem neuropathology and brain gene expression in aging. Here we present a resource of RNA sequencing of purified CD14+ human monocytes - including an eQTL map - from over 200 elderly individuals, most with accompanying bulk brain RNA sequencing profiles, longitudinal cognitive assessments, and detailed postmortem neuropathological examinations. We tested the direct correlation of gene expression between monocytes and bulk brain tissue, finding very few significant signals driven largely by genetic variation. However, we did identify sets of monocyte-expressed genes that were highly predictive of postmortem microglial activation, diffuse amyloid plaque deposition, and cerebrovascular disease. Our findings prioritize potential blood-based molecular biomarkers for AD; they also reveal the previously unknown architecture of shared gene expression between the CNS and peripheral immune system in aging.

Identifying genes underlying memory function will help characterize cognitively resilient and high-risk declining subpopulations contributing to precision medicine strategies. We estimated episodic memory trajectories in 35,245 ethnically diverse older adults representing eight independent cohorts. We conducted apolipoprotein E (APOE)-stratified genome-wide association study (GWAS) analyses and combined individual cohorts' results via meta-analysis. Three independent transcriptomics datasets were used to further interpret GWAS signals. We identified DCDC2 gene significantly associated with episodic memory (Pmeta = 3.3 x 10