The human microglia responsome: a resource to better understand microglia states in health and disease.

Microglia, the immune cells of the brain, are increasingly implicated in neurodegenerative disorders through genetic studies. However, how genetic risk factors for these diseases are related to microglial gene expression, microglial function, and ultimately disease, is still largely unknown. Microglia change rapidly in response to alterations in their cellular environment, which is regulated through changes in transcriptional programs, which are as yet poorly understood. Here, we compared the effects of a set of inflammatory and restorative stimuli (lipopolysaccharide, interferon-gamma, resiquimod, tumor necrosis factor-alpha, adenosine triphosphate, dexamethasone, and interleukin-4) on human microglial cells from 67 different donors (N = 398 samples) at the gene and transcript level. We show that microglia from different anatomical brain regions show distinct responses to inflammatory stimuli. We observed a greater overlap between human stimulated microglia and human monocytes than with mouse microglia. We define specific microglial signatures across conditions which are highly relevant for a wide range of biological functions and complex human diseases. Finally, we used our stimulation signatures to interpret associations from Alzheimer's disease (AD) genetic studies and microglia by integrating our inflammatory gene expression profiles with common genetic variants to map cis -expression QTLs (eQTLs). Together, we provide the most comprehensive transcriptomic database of the human microglia responsome. Highlights: RNA-sequencing of 398 human microglial samples exposed to six different triggers.Microglia from different anatomical regions show distinct stimulation responses.Responses in human microglia show a greater overlap with human monocytes than murine microglia.Mapping of response Quantitative Trait Loci identifies interactions between genotype and effect of stimulation on gene expression.Our atlas provides a reference map for interpreting microglia signatures in health and disease.

HnRNP K mislocalisation is a novel protein pathology of frontotemporal lobar degeneration and ageing and leads to cryptic splicing.

Heterogeneous nuclear ribonucleoproteins (HnRNPs) are a group of ubiquitously expressed RNA-binding proteins implicated in the regulation of all aspects of nucleic acid metabolism. HnRNP K is a member of this highly versatile hnRNP family. Pathological redistribution of hnRNP K to the cytoplasm has been linked to the pathogenesis of several malignancies but, until now, has been underexplored in the context of neurodegenerative disease. Here we show hnRNP K mislocalisation in pyramidal neurons of the frontal cortex to be a novel neuropathological feature that is associated with both frontotemporal lobar degeneration and ageing. HnRNP K mislocalisation is mutually exclusive to TDP-43 and tau pathological inclusions in neurons and was not observed to colocalise with mitochondrial, autophagosomal or stress granule markers. De-repression of cryptic exons in RNA targets following TDP-43 nuclear depletion is an emerging mechanism of potential neurotoxicity in frontotemporal lobar degeneration and the mechanistically overlapping disorder amyotrophic lateral sclerosis. We silenced hnRNP K in neuronal cells to identify the transcriptomic consequences of hnRNP K nuclear depletion. Intriguingly, by performing RNA-seq analysis we find that depletion of hnRNP K induces 101 novel cryptic exon events. We validated cryptic exon inclusion in an SH-SY5Y hnRNP K knockdown and in FTLD brain exhibiting hnRNP K nuclear depletion. We, therefore, present evidence for hnRNP K mislocalisation to be associated with FTLD and for this to induce widespread changes in splicing.

Truncated stathmin-2 is a marker of TDP-43 pathology in frontotemporal dementia.

No treatment for frontotemporal dementia (FTD), the second most common type of early-onset dementia, is available, but therapeutics are being investigated to target the 2 main proteins associated with FTD pathological subtypes: TDP-43 (FTLD-TDP) and tau (FTLD-tau). Testing potential therapies in clinical trials is hampered by our inability to distinguish between patients with FTLD-TDP and FTLD-tau. Therefore, we evaluated truncated stathmin-2 (STMN2) as a proxy of TDP-43 pathology, given the reports that TDP-43 dysfunction causes truncated STMN2 accumulation. Truncated STMN2 accumulated in human induced pluripotent stem cell-derived neurons depleted of TDP-43, but not in those with pathogenic TARDBP mutations in the absence of TDP-43 aggregation or loss of nuclear protein. In RNA-Seq analyses of human brain samples from the NYGC ALS cohort, truncated STMN2 RNA was confined to tissues and disease subtypes marked by TDP-43 inclusions. Last, we validated that truncated STMN2 RNA was elevated in the frontal cortex of a cohort of patients with FTLD-TDP but not in controls or patients with progressive supranuclear palsy, a type of FTLD-tau. Further, in patients with FTLD-TDP, we observed significant associations of truncated STMN2 RNA with phosphorylated TDP-43 levels and an earlier age of disease onset. Overall, our data uncovered truncated STMN2 as a marker for TDP-43 dysfunction in FTD.

Collagenous Colitis Is Associated With HLA Signature and Shares Genetic Risks With Other Immune-Mediated Diseases.

BACKGROUND & AIMS: Collagenous colitis (CC) is an inflammatory bowel disorder with unknown etiopathogenesis involving HLA-related immune-mediated responses and environmental and genetic risk factors. We carried out an array-based genetic association study in a cohort of patients with CC and investigated the common genetic basis between CC and Crohn's disease (CD), ulcerative colitis (UC), and celiac disease. METHODS: DNA from 804 CC formalin-fixed, paraffin-embedded tissue samples was genotyped with Illumina Immunochip. Matching genotype data on control samples and CD, UC, and celiac disease cases were provided by the respective consortia. A discovery association study followed by meta-analysis with an independent cohort, polygenic risk score calculation, and cross-phenotype analyses were performed. Enrichment of regulatory expression quantitative trait loci among the CC variants was assessed in hemopoietic and intestinal cells. RESULTS: Three HLA alleles (HLA-B∗08:01, HLA-DRB1∗03:01, and HLA-DQB1∗02:01), related to the ancestral haplotype 8.1, were significantly associated with increased CC risk. We also identified an independent protective effect of HLA-DRB1∗04:01 on CC risk. Polygenic risk score quantifying the risk across multiple susceptibility loci was strongly associated with CC risk. An enrichment of expression quantitative trait loci was detected among the CC-susceptibility variants in various cell types. The cross-phenotype analysis identified a complex pattern of polygenic pleiotropy between CC and other immune-mediated diseases. CONCLUSIONS: In this largest genetic study of CC to date with histologically confirmed diagnosis, we strongly implicated the HLA locus and proposed potential non-HLA mechanisms in disease pathogenesis. We also detected a shared genetic risk between CC, celiac disease, CD, and UC, which supports clinical observations of comorbidity.

Tensor decomposition of stimulated monocyte and macrophage gene expression profiles identifies neurodegenerative disease-specific trans-eQTLs.

Recent human genetic studies suggest that cells of the innate immune system have a primary role in the pathogenesis of neurodegenerative diseases. However, the results from these studies often do not elucidate how the genetic variants affect the biology of these cells to modulate disease risk. Here, we applied a tensor decomposition method to uncover disease associated gene networks linked to distal genetic variation in stimulated human monocyte and macrophage gene expression profiles. We report robust evidence that some disease associated genetic variants affect the expression of multiple genes in trans. These include a Parkinson's disease locus influencing the expression of genes mediated by a protease that controls lysosomal function, and Alzheimer's disease loci influencing the expression of genes involved in type 1 interferon signaling, myeloid phagocytosis, and complement cascade pathways. Overall, we uncover gene networks in induced innate immune cells linked to disease associated genetic variants, which may help elucidate the underlying biology of disease.

A common haplotype lowers PU.1 expression in myeloid cells and delays onset of Alzheimer's disease.

A genome-wide survival analysis of 14,406 Alzheimer's disease (AD) cases and 25,849 controls identified eight previously reported AD risk loci and 14 novel loci associated with age at onset. Linkage disequilibrium score regression of 220 cell types implicated the regulation of myeloid gene expression in AD risk. The minor allele of rs1057233 (G), within the previously reported CELF1 AD risk locus, showed association with delayed AD onset and lower expression of SPI1 in monocytes and macrophages. SPI1 encodes PU.1, a transcription factor critical for myeloid cell development and function. AD heritability was enriched within the PU.1 cistrome, implicating a myeloid PU.1 target gene network in AD. Finally, experimentally altered PU.1 levels affected the expression of mouse orthologs of many AD risk genes and the phagocytic activity of mouse microglial cells. Our results suggest that lower SPI1 expression reduces AD risk by regulating myeloid gene expression and cell function.

Alzheimer's disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci.

We used a collection of 708 prospectively collected autopsied brains to assess the methylation state of the brain's DNA in relation to Alzheimer's disease (AD). We found that the level of methylation at 71 of the 415,848 interrogated CpGs was significantly associated with the burden of AD pathology, including CpGs in the ABCA7 and BIN1 regions, which harbor known AD susceptibility variants. We validated 11 of the differentially methylated regions in an independent set of 117 subjects. Furthermore, we functionally validated these CpG associations and identified the nearby genes whose RNA expression was altered in AD: ANK1, CDH23, DIP2A, RHBDF2, RPL13, SERPINF1 and SERPINF2. Our analyses suggest that these DNA methylation changes may have a role in the onset of AD given that we observed them in presymptomatic subjects and that six of the validated genes connect to a known AD susceptibility gene network.

Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Fine-mapping the genetic association of the major histocompatibility complex in multiple sclerosis: HLA and non-HLA effects.

The major histocompatibility complex (MHC) region is strongly associated with multiple sclerosis (MS) susceptibility. HLA-DRB1*15:01 has the strongest effect, and several other alleles have been reported at different levels of validation. Using SNP data from genome-wide studies, we imputed and tested classical alleles and amino acid polymorphisms in 8 classical human leukocyte antigen (HLA) genes in 5,091 cases and 9,595 controls. We identified 11 statistically independent effects overall: 6 HLA-DRB1 and one DPB1 alleles in class II, one HLA-A and two B alleles in class I, and one signal in a region spanning from MICB to LST1. This genomic segment does not contain any HLA class I or II genes and provides robust evidence for the involvement of a non-HLA risk allele within the MHC. Interestingly, this region contains the TNF gene, the cognate ligand of the well-validated TNFRSF1A MS susceptibility gene. The classical HLA effects can be explained to some extent by polymorphic amino acid positions in the peptide-binding grooves. This study dissects the independent effects in the MHC, a critical region for MS susceptibility that harbors multiple risk alleles.

Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis.

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8 × 10(-8)), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3' UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1 × 10(-11) in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry.

Alzheimer disease susceptibility loci: evidence for a protein network under natural selection.

Recent genome-wide association studies have identified a number of susceptibility loci for Alzheimer disease (AD). To understand the functional consequences and potential interactions of the associated loci, we explored large-scale data sets interrogating the human genome for evidence of positive natural selection. Our findings provide significant evidence for signatures of recent positive selection acting on several haplotypes carrying AD susceptibility alleles; interestingly, the genes found in these selected haplotypes can be assembled, independently, into a molecular complex via a protein-protein interaction (PPI) network approach. These results suggest a possible coevolution of genes encoding physically-interacting proteins that underlie AD susceptibility and are coexpressed in different tissues. In particular, PICALM, BIN1, CD2AP, and EPHA1 are interconnected through multiple interacting proteins and appear to have coordinated evidence of selection in the same human population, suggesting that they may be involved in the execution of a shared molecular function. This observation may be AD-specific, as the 12 loci associated with Parkinson disease do not demonstrate excess evidence of natural selection. The context for selection is probably unrelated to AD itself; it is likely that these genes interact in another context, such as in immune cells, where we observe cis-regulatory effects at several of the selected AD loci.