BHLHE40/41 regulate microglia and peripheral macrophage responses associated with Alzheimer's disease and other disorders of lipid-rich tissues

Genetic and experimental evidence suggests that Alzheimer's disease (AD) risk alleles and genes may influence disease susceptibility by altering the transcriptional and cellular responses of macrophages, including microglia, to damage of lipid-rich tissues like the brain. Recently, sc/nRNA sequencing studies identified similar transcriptional activation states in subpopulations of macrophages in aging and degenerating brains and in other diseased lipid-rich tissues. We collectively refer to these subpopulations of microglia and peripheral macrophages as DLAMs. Using macrophage sc/nRNA-seq data from healthy and diseased human and mouse lipid-rich tissues, we reconstructed gene regulatory networks and identified 11 strong candidate transcriptional regulators of the DLAM response across species. Loss or reduction of two of these transcription factors, BHLHE40/41, in iPSC-derived microglia and human THP-1 macrophages as well as loss of Bhlhe40/41 in mouse microglia, resulted in increased expression of DLAM genes involved in cholesterol clearance and lysosomal processing, increased cholesterol efflux and storage, and increased lysosomal mass and degradative capacity. These findings provide targets for therapeutic modulation of macrophage/microglial function in AD and other disorders affecting lipid-rich tissues.

Autosomal Dominant Alzheimer's Disease Mutations in Human Microglia Are Not Sufficient to Trigger Amyloid Pathology in WT Mice but Might Affect Pathology in 5XFAD Mice

Several genetic variants that affect microglia function have been identified as risk factors for Alzheimer's Disease (AD), supporting the importance of this cell type in disease progression. However, the effect of autosomal dominant mutations in the amyloid precursor protein (APP) or the presenilin (PSEN1/2) genes has not been addressed in microglia in vivo. We xenotransplanted human microglia derived from non-carriers and carriers of autosomal dominant AD (ADAD)-causing mutations in the brain of hCSF1 WT or 5XFAD mice. We observed that ADAD mutations in microglia are not sufficient to trigger amyloid pathology in WT mice. In 5XFAD mice, we observed a non-statistically significant increase in amyloid plaque volume and number of dystrophic neurites, coupled with a reduction in plaque-associated microglia in the brain of mice xenotransplanted with ADAD human microglia compared to mice xenotransplanted with non-ADAD microglia. In addition, we observed a non-statistically significant impairment in working and contextual memory in 5XFAD mice xenotransplanted with ADAD microglia compared to those xenotransplanted with non-ADAD-carrier microglia. We conclude that, although not sufficient to initiate amyloid pathology in the healthy brain, mutations in APP and PSEN1 in human microglia might cause mild changes in pathological and cognitive outcomes in 5XFAD mice in a manner consistent with increased AD risk.

A novel molecular class that recruits HDAC/MECP2 complexes to PU.1 motifs reduces neuroinflammation

Pervasive neuroinflammation occurs in many neurodegenerative diseases, including Alzheimer's disease (AD). SPI1/PU.1 is a transcription factor located at a genome-wide significant AD-risk locus and its reduced expression is associated with delayed onset of AD. We analyzed single-cell transcriptomic datasets from microglia of human AD patients and found an enrichment of PU.1-binding motifs in the differentially expressed genes. In hippocampal tissues from transgenic mice with neurodegeneration, we found vastly increased genomic PU.1 binding. We then screened for PU.1 inhibitors using a PU.1 reporter cell line and discovered A11, a molecule with anti-inflammatory efficacy and nanomolar potency. A11 regulated genes putatively by recruiting a repressive complex containing MECP2, HDAC1, SIN3A, and DNMT3A to PU.1 motifs, thus representing a novel mechanism and class of molecules. In mouse models of AD, A11 ameliorated neuroinflammation, loss of neuronal integrity, AD pathology, and improved cognitive performance. This study uncovers a novel class of anti-inflammatory molecules with therapeutic potential for neurodegenerative disorders.

APOE and immunity: Research highlights

INTRODUCTION

At the Alzheimer's Association's APOE and Immunity virtual conference, held in October 2021, leading neuroscience experts shared recent research advances on and inspiring insights into the various roles that both the apolipoprotein E gene (APOE) and facets of immunity play in neurodegenerative diseases, including Alzheimer's disease and other dementias.

METHODS

The meeting brought together more than 1200 registered attendees from 62 different countries, representing the realms of academia and industry.

RESULTS

During the 4-day meeting, presenters illuminated aspects of the cross-talk between APOE and immunity, with a focus on the roles of microglia, triggering receptor expressed on myeloid cells 2 (TREM2), and components of inflammation (e.g., tumor necrosis factor α [TNFα]).

DISCUSSION

This manuscript emphasizes the importance of diversity in current and future research and presents an integrated view of innate immune functions in Alzheimer's disease as well as related promising directions in drug development.

The complex genetic architecture of Alzheimer's disease: novel insights and future directions

BACKGROUND

Alzheimer's disease (AD) is a complex multifactorial neurodegenerative disorder and the most common form of dementia. AD is highly heritable, with heritability estimates of ∼70% from twin studies. Progressively larger genome-wide association studies (GWAS) have continued to expand our knowledge of AD/dementia genetic architecture. Until recently these efforts had identified 39 disease susceptibility loci in European ancestry populations.

RECENT DEVELOPMENTS

Two new AD/dementia GWAS have dramatically expanded the sample sizes and the number of disease susceptibility loci. The first increased total sample size to 1,126,563-with an effective sample size of 332,376-by predominantly including new biobank and population-based dementia datasets. The second, expands on an earlier GWAS from the International Genomics of Alzheimer's Project (IGAP) by increasing the number of clinically-defined AD cases/controls in addition to incorporating biobank dementia datasets, resulting in a total sample size to 788,989 and an effective sample size of 382,472. Collectively both GWAS identified 90 independent variants across 75 AD/dementia susceptibility loci, including 42 novel loci. Pathway analyses indicate the susceptibility loci are enriched for genes involved in amyloid plaque and neurofibrillary tangle formation, cholesterol metabolism, endocytosis/phagocytosis, and the innate immune system. Gene prioritization efforts for the novel loci identified 62 candidate causal genes. Many of the candidate genes from known and newly discovered loci play key roles in macrophages and highlight phagocytic clearance of cholesterol-rich brain tissue debris by microglia (efferocytosis) as a core pathogenetic hub and putative therapeutic target for AD. WHERE NEXT?: While GWAS in European ancestry populations have substantially enhanced our understanding of AD genetic architecture, heritability estimates from population based GWAS cohorts are markedly smaller than those from twin studies. While this missing heritability is likely due to a combination of factors, it highlights that our understanding of AD genetic architecture and genetic risk mechanisms remains incomplete. These knowledge gaps result from several underexplored areas in AD research. First, rare variants remain understudied due to methodological issues in identifying them and the cost of generating sufficiently powered whole exome/genome sequencing datasets. Second, sample sizes of non-European ancestry populations in AD GWAS remain small. Third, GWAS of AD neuroimaging and cerebrospinal fluid endophenotypes remains limited due to low compliance and high costs associated with measuring amyloid-β and tau levels and other disease-relevant biomarkers. Studies generating sequencing data, including diverse populations, and incorporating blood-based AD biomarkers are set to substantially improve our knowledge of AD genetic architecture.

BHLHE40/41 regulate macrophage/microglia responses associated with Alzheimer's disease and other disorders of lipid-rich tissues

Background

Genetic and experimental evidence strongly implicates myeloid cells in the etiology of AD and suggests that AD-associated alleles and genes may modulate disease risk by altering the transcriptional and cellular responses of macrophages (like microglia) to damage of lipid-rich tissues (like the brain). Specifically, recent single-cell/nucleus RNA sequencing (sc/nRNA-seq) studies identified a transcriptionally distinct state of subsets of macrophages in aging or degenerating brains (usually referred to as disease-associated microglia or DAM) and in other diseased lipid-rich tissues (e.g., obese adipose tissue, fatty liver, and atherosclerotic plaques). We collectively refer to these subpopulations as lipid-associated macrophages or LAMs. Importantly, this particular activation state is characterized by increased expression of genes involved in the phagocytic clearance of lipid-rich cellular debris (efferocytosis), including several AD risk genes.

Methods

We used sc/nRNA-seq data from human and mouse microglia from healthy and diseased brains and macrophages from other lipid-rich tissues to reconstruct gene regulatory networks and identify transcriptional regulators whose regulons are enriched for LAM response genes (LAM TFs) across species. We then used gene knock-down/knock-out strategies to validate some of these LAM TFs in human THP-1 macrophages and iPSC-derived microglia in vitro, as well as mouse microglia in vivo.

Results

We nominate 11 strong candidate LAM TFs shared across human and mouse networks (BHLHE41, HIF1A, ID2, JUNB, MAF, MAFB, MEF2A, MEF2C, NACA, POU2F2 and SPI1). We also demonstrate a strong enrichment of AD risk alleles in the cistrome of BHLHE41 (and its close homolog BHLHE40), thus implicating its regulon in the modulation of disease susceptibility. Loss or reduction of BHLHE40/41 expression in human THP-1 macrophages and iPSC-derived microglia, as well as loss of Bhlhe40/41 in mouse microglia led to increased expression of LAM response genes, specifically those involved in cholesterol clearance and lysosomal processing, with a concomitant increase in cholesterol efflux and storage, as well as lysosomal mass and degradative capacity.

Conclusions

Taken together, this study nominates transcriptional regulators of the LAM response, experimentally validates BHLHE40/41 in human and mouse macrophages/microglia, and provides novel targets for therapeutic modulation of macrophage/microglia function in AD and other disorders of lipid-rich tissues.

Microglial efferocytosis: Diving into the Alzheimer's disease gene pool

Genome-wide association studies and functional genomics studies have linked specific cell types, genes, and pathways to Alzheimer's disease (AD) risk. In particular, AD risk alleles primarily affect the abundance or structure, and thus the activity, of genes expressed in macrophages, strongly implicating microglia (the brain-resident macrophages) in the etiology of AD. These genes converge on pathways (endocytosis/phagocytosis, cholesterol metabolism, and immune response) with critical roles in core macrophage functions such as efferocytosis. Here, we review these pathways, highlighting relevant genes identified in the latest AD genetics and genomics studies, and describe how they may contribute to AD pathogenesis. Investigating the functional impact of AD-associated variants and genes in microglia is essential for elucidating disease risk mechanisms and developing effective therapeutic approaches.

APOE4 confers transcriptomic and functional alterations to primary mouse microglia

Common genetic variants in more than forty loci modulate risk for Alzheimer's disease (AD). AD risk alleles are enriched within enhancers active in myeloid cells, suggesting that microglia, the brain-resident macrophages, may play a key role in the etiology of AD. A major genetic risk factor for AD is Apolipoprotein E (APOE) genotype, with the ε4/ε4 (E4) genotype increasing risk for AD by approximately 15 fold compared to the most common ε3/ε3 (E3) genotype. However, the impact of APOE genotype on microglial function has not been thoroughly investigated. To address this, we cultured primary microglia from mice in which both alleles of the mouse Apoe gene have been humanized to encode either human APOE ε3 or APOE ε4. Relative to E3 microglia, E4 microglia exhibit altered morphology, increased endolysosomal mass, increased cytokine/chemokine production, and increased lipid and lipid droplet accumulation at baseline. These changes were accompanied by decreased translation and increased phosphorylation of eIF2ɑ and eIF2ɑ-kinases that participate in the integrated stress response, suggesting that E4 genotype leads to elevated levels of cellular stress in microglia relative to E3 genotype. Using live-cell imaging and flow cytometry, we also show that E4 microglia exhibited increased phagocytic uptake of myelin and other substrates compared to E3 microglia. While transcriptomic profiling of myelin-challenged microglia revealed a largely overlapping response profile across genotypes, differential enrichment of genes in interferon signaling, extracellular matrix and translation-related pathways was identified in E4 versus E3 microglia both at baseline and following myelin challenge. Together, our results suggest E4 genotype confers several important functional alterations to microglia even prior to myelin challenge, providing insight into the molecular and cellular mechanisms by which APOE4 may increase risk for AD.

Genome-wide association study and functional validation implicates JADE1 in tauopathy

Primary age-related tauopathy (PART) is a neurodegenerative pathology with features distinct from but also overlapping with Alzheimer disease (AD). While both exhibit Alzheimer-type temporal lobe neurofibrillary degeneration alongside amnestic cognitive impairment, PART develops independently of amyloid-β (Aβ) plaques. The pathogenesis of PART is not known, but evidence suggests an association with genes that promote tau pathology and others that protect from Aβ toxicity. Here, we performed a genetic association study in an autopsy cohort of individuals with PART (n = 647) using Braak neurofibrillary tangle stage as a quantitative trait. We found some significant associations with candidate loci associated with AD (SLC24A4, MS4A6A, HS3ST1) and progressive supranuclear palsy (MAPT and EIF2AK3). Genome-wide association analysis revealed a novel significant association with a single nucleotide polymorphism on chromosome 4 (rs56405341) in a locus containing three genes, including JADE1 which was significantly upregulated in tangle-bearing neurons by single-soma RNA-seq. Immunohistochemical studies using antisera targeting JADE1 protein revealed localization within tau aggregates in autopsy brains with four microtubule-binding domain repeats (4R) isoforms and mixed 3R/4R, but not with 3R exclusively. Co-immunoprecipitation in post-mortem human PART brain tissue revealed a specific binding of JADE1 protein to four repeat tau lacking N-terminal inserts (0N4R). Finally, knockdown of the Drosophila JADE1 homolog rhinoceros (rno) enhanced tau-induced toxicity and apoptosis in vivo in a humanized 0N4R mutant tau knock-in model, as quantified by rough eye phenotype and terminal deoxynucleotidyl transferase dUTP nick end-labeling (TUNEL) in the fly brain. Together, these findings indicate that PART has a genetic architecture that partially overlaps with AD and other tauopathies and suggests a novel role for JADE1 as a modifier of neurofibrillary degeneration.

Whole genome sequencing identifies loci specifically associated with thoracic aortic wall defects and abdominal aortic aneurysms in patients with European ancestry

Objective

The objective of this study was to better understand the pathophysiology and underlying genetic mechanisms behind two abdominal aortic aneurysm (AAA) subtypes using computed tomographic imaging in combination with whole genome sequencing.

Methods

Patients with a known AAA and European ancestry were included in this investigation and underwent genetic and image analysis. Patients with AAAs and indications of descending thoracic aortic pathology (aortic dissection, penetrating aortic ulcers, intramural hematoma, atheromas, ulcerative plaque, and intramural ulceration, and intimal flaps/tears) were classified as having thoracic aortic disease, grouped together, and compared with patients with an AAA and a normal descending thoracic aorta. Whole genome sequencing was then performed on the 93 patients who had imaging features consistent with thoracic aortic disease and the 126 patients with a normal descending thoracic aorta.

Results

The results of this study suggest one variant-level, four gene-level, and one gene set-level associations in patients with thoracic aortic disease who also had an AAA. The variant rs79508780 located in TSEN54 achieved study-wide significance (P = 1.71E-06). BATF3 and SMLR1 were significantly associated and EFCAB3 and TAF4 were reached suggestive assocation with a diseased descending thoracic aorta (P = 5.23E-26, P = 1.86E-25, P = 1.54E-05, and P = 8.31E-05, respectively). Gene sets were also compiled using MSigDB and trait-based index single nucleotide variation from major genome-wide association studies. GO_DNA_DOUBLE_STRAND_BREAK_PROCESSING, a gene set related to double-stranded DNA break repair, was significantly associated with thoracic aortic disease in AAA patients (P = 1.80E-06).

Conclusions

This pilot study provides further evidence that an AAA may be the end result of multiple degenerative pathways. Genetic variations in vitamin D signaling, cholesterol metabolism, extracellular matrix breakdown, and double-stranded DNA break repair pathways were associated with European patients who had an AAA and thoracic aortic disease. Additionally, this study provides support for the application of a radiogenomic approach for the investigation of other potential pathologies that could lead to the development of an AAA or influence future management decisions. (JVS–Vascular Science.)

Clinical Relevance

In this study, we provide evidence that abdominal aortic aneurysms (AAAs) may be a result of multiple pathophysiologies rather than a single disease. We have identified genetic variants involved in vitamin D signaling, cholesterol metabolism, extracellular matrix breakdown, and double-stranded DNA break repair associated with structural defects in the aortic wall in patients with AAAs who are of European descent. Patients with AAAs and structural defects in the thoracic aorta have been previously linked to differential behavior after endovascular aneurysm repair. These patients with wall defects exhibited greater sac regression, a marker of surgical success, after endovascular aneurysm repair. Our study demonstrates the usefulness of a radiogenomic approach for elucidating mechanisms behind the formation and future behavior of AAAs that could aid surgeons in making future procedural and management decisions.

Multi-omics integration analysis identifies novel genes for alcoholism with potential overlap with neurodegenerative diseases

Identification of causal variants and genes underlying genome-wide association study (GWAS) loci is essential to understand the biology of alcohol use disorder (AUD) and drinks per week (DPW). Multi-omics integration approaches have shown potential for fine mapping complex loci to obtain biological insights to disease mechanisms. In this study, we use multi-omics approaches, to fine-map AUD and DPW associations at single SNP resolution to demonstrate that rs56030824 on chromosome 11 significantly reduces SPI1 mRNA expression in myeloid cells and lowers risk for AUD and DPW. Our analysis also identifies MAPT as a candidate causal gene specifically associated with DPW. Genes prioritized in this study show overlap with causal genes associated with neurodegenerative disorders. Multi-omics integration analyses highlight, genetic similarities and differences between alcohol intake and disordered drinking, suggesting molecular heterogeneity that might inform future targeted functional and cross-species studies.

Beyond association: successes and challenges in linking non-coding genetic variation to functional consequences that modulate Alzheimer's disease risk

Alzheimer's disease (AD) is the most common type of dementia, affecting millions of people worldwide; however, no disease-modifying treatments are currently available. Genome-wide association studies (GWASs) have identified more than 40 loci associated with AD risk. However, most of the disease-associated variants reside in non-coding regions of the genome, making it difficult to elucidate how they affect disease susceptibility. Nonetheless, identification of the regulatory elements, genes, pathways and cell type/tissue(s) impacted by these variants to modulate AD risk is critical to our understanding of disease pathogenesis and ability to develop effective therapeutics. In this review, we provide an overview of the methods and approaches used in the field to identify the functional effects of AD risk variants in the causal path to disease risk modification as well as describe the most recent findings. We first discuss efforts in cell type/tissue prioritization followed by recent progress in candidate causal variant and gene nomination. We discuss statistical methods for fine-mapping as well as approaches that integrate multiple levels of evidence, such as epigenomic and transcriptomic data, to identify causal variants and risk mechanisms of AD-associated loci. Additionally, we discuss experimental approaches and data resources that will be needed to validate and further elucidate the effects of these variants and genes on biological pathways, cellular phenotypes and disease risk. Finally, we discuss future steps that need to be taken to ensure that AD GWAS functional mapping efforts lead to novel findings and bring us closer to finding effective treatments for this devastating disease.

Integration of Alzheimer's disease genetics and myeloid genomics identifies disease risk regulatory elements and genes

Genome-wide association studies (GWAS) have identified more than 40 loci associated with Alzheimer's disease (AD), but the causal variants, regulatory elements, genes and pathways remain largely unknown, impeding a mechanistic understanding of AD pathogenesis. Previously, we showed that AD risk alleles are enriched in myeloid-specific epigenomic annotations. Here, we show that they are specifically enriched in active enhancers of monocytes, macrophages and microglia. We integrated AD GWAS with myeloid epigenomic and transcriptomic datasets using analytical approaches to link myeloid enhancer activity to target gene expression regulation and AD risk modification. We identify AD risk enhancers and nominate candidate causal genes among their likely targets (including AP4E1, AP4M1, APBB3, BIN1, MS4A4A, MS4A6A, PILRA, RABEP1, SPI1, TP53INP1, and ZYX) in twenty loci. Fine-mapping of these enhancers nominates candidate functional variants that likely modify AD risk by regulating gene expression in myeloid cells. In the MS4A locus we identified a single candidate functional variant and validated it in human induced pluripotent stem cell (hiPSC)-derived microglia and brain. Taken together, this study integrates AD GWAS with multiple myeloid genomic datasets to investigate the mechanisms of AD risk alleles and nominates candidate functional variants, regulatory elements and genes that likely modulate disease susceptibility.

The variant call format provides efficient and robust storage of GWAS summary statistics

GWAS summary statistics are fundamental for a variety of research applications yet no common storage format has been widely adopted. Existing tabular formats ambiguously or incompletely store information about genetic variants and associations, lack essential metadata and are typically not indexed yielding poor query performance and increasing the possibility of errors in data interpretation and post-GWAS analyses. To address these issues, we adapted the variant call format to store GWAS summary statistics (GWAS-VCF) and developed open-source tools to use this format in downstream analyses. We provide open access to over 10,000 complete GWAS summary datasets converted to this format ( https://gwas.mrcieu.ac.uk ).

Alzheimer's-associated PU.1 expression levels regulate microglial inflammatory response

More than forty loci contribute to genetic risk for Alzheimer's disease (AD). These risk alleles are enriched in myeloid cell enhancers suggesting that microglia, the brain-resident macrophages, contribute to AD risk. We have previously identified SPI1/PU.1, a master regulator of myeloid cell development in the brain and periphery, as a genetic risk factor for AD. Higher expression of SPI1 is associated with increased risk for AD, while lower expression is protective. To investigate the molecular and cellular phenotypes associated with higher and lower expression of PU.1 in microglia, we used stable overexpression and knock-down of PU.1 in BV2, an immortalized mouse microglial cell line. Transcriptome analysis suggests that reduced PU.1 expression suppresses expression of homeostatic genes similar to the disease-associated microglia response to amyloid plaques in mouse models of AD. Moreover, PU.1 knock-down resulted in activation of protein translation, antioxidant action and cholesterol/lipid metabolism pathways with a concomitant decrease of pro-inflammatory gene expression. PU.1 overexpression upregulated and knock-down downregulated phagocytic uptake in BV2 cells independent of the nature of the engulfed material. However, cells with reduced PU.1 expression retained their ability to internalize myelin similar to control albeit with a delay, which aligns with their anti-inflammatory profile. Here we identified several microglial responses that are modulated by PU.1 expression levels and propose that risk association of PU.1 to AD is driven by increased pro-inflammatory response due to increased viability of cells under cytotoxic conditions. In contrast, low expression of PU.1 leads to increased cell death under cytotoxic conditions accompanied by reduced pro-inflammatory signaling that decreased A1 reactive astrocytes signature supporting the protective effect of SPI1 genotype in AD. These findings inform future in vivo validation studies and design of small molecule screens for therapeutic discovery in AD.

Causal Associations Between Modifiable Risk Factors and the Alzheimer's Phenome

OBJECTIVE

The purpose of this study was to infer causal relationships between 22 previously reported risk factors for Alzheimer's disease (AD) and the "AD phenome": AD, AD age of onset (AAOS), hippocampal volume, cortical surface area and thickness, cerebrospinal fluid (CSF) levels of amyloid-β (Aβ42 ), tau, and ptau181 , and the neuropathological burden of neuritic plaques, neurofibrillary tangles (NFTs), and vascular brain injury (VBI).

METHODS

Polygenic risk scores (PRS) for the 22 risk factors were computed in 26,431 AD cases/controls and the association with AD was evaluated using logistic regression. Two-sample Mendelian randomization (MR) was used to infer the causal effect of risk factors on the AD phenome.

RESULTS

PRS for increased education and diastolic blood pressure were associated with reduced risk for AD. MR indicated that only education was causally associated with reduced risk of AD, delayed AAOS, and increased cortical surface area and thickness. Total- and LDL-cholesterol levels were causally associated with increased neuritic plaque burden, although the effects were driven by single nucleotide polymorphisms (SNPs) within the APOE locus. Diastolic blood pressure and pulse pressure are causally associated with increased risk of VBI. Furthermore, total cholesterol was associated with decreased hippocampal volume; smoking initiation with decreased cortical thickness; type 2 diabetes with an earlier AAOS; and sleep duration with increased cortical thickness.

INTERPRETATION

Our comprehensive examination of the genetic evidence for the causal relationships between previously reported risk factors in AD using PRS and MR supports a causal role for education, blood pressure, cholesterol levels, smoking, and diabetes with the AD phenome. ANN NEUROL 2021;89:54-65.

Risk prediction of late-onset Alzheimer's disease implies an oligogenic architecture

Genetic association studies have identified 44 common genome-wide significant risk loci for late-onset Alzheimer's disease (LOAD). However, LOAD genetic architecture and prediction are unclear. Here we estimate the optimal P-threshold (Poptimal) of a genetic risk score (GRS) for prediction of LOAD in three independent datasets comprising 676 cases and 35,675 family history proxy cases. We show that the discriminative ability of GRS in LOAD prediction is maximised when selecting a small number of SNPs. Both simulation results and direct estimation indicate that the number of causal common SNPs for LOAD may be less than 100, suggesting LOAD is more oligogenic than polygenic. The best GRS explains approximately 75% of SNP-heritability, and individuals in the top decile of GRS have ten-fold increased odds when compared to those in the bottom decile. In addition, 14 variants are identified that contribute to both LOAD risk and age at onset of LOAD.

Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation

The Alzheimer's Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10-7), an immunoglobulin gene whose antibodies interact with β-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10-7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10-6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

Correction: Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation

A correction to this paper has been published and can be accessed via a link at the top of the paper.

Functional annotation of genomic variants in studies of late-onset Alzheimer's disease

Motivation

Annotation of genomic variants is an increasingly important and complex part of the analysis of sequence-based genomic analyses. Computational predictions of variant function are routinely incorporated into gene-based analyses of rare-variants, though to date most studies use limited information for assessing variant function that is often agnostic of the disease being studied.

Results

In this work, we outline an annotation process motivated by the Alzheimer's Disease Sequencing Project, illustrate the impact of including tissue-specific transcript sets and sources of gene regulatory information and assess the potential impact of changing genomic builds on the annotation process. While these factors only impact a small proportion of total variant annotations (∼5%), they influence the potential analysis of a large fraction of genes (∼25%).

Availability and implementation

Individual variant annotations are available via the NIAGADS GenomicsDB, at https://www.niagads.org/genomics/ tools-and-software/databases/genomics-database. Annotations are also available for bulk download at https://www.niagads.org/datasets. Annotation processing software is available at http://www.icompbio.net/resources/software-and-downloads/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Molecular basis for the loss-of-function effects of the Alzheimer's disease-associated R47H variant of the immune receptor TREM2

Triggering receptor expressed on myeloid cells 2 (TREM2) is an immune receptor expressed on the surface of microglia, macrophages, dendritic cells, and osteoclasts. The R47H TREM2 variant is a significant risk factor for late-onset Alzheimer's disease (AD), and the molecular basis of R47H TREM2 loss of function is an emerging area of TREM2 biology. Here, we report three high-resolution structures of the extracellular ligand-binding domains (ECDs) of R47H TREM2, apo-WT, and phosphatidylserine (PS)-bound WT TREM2 at 1.8, 2.2, and 2.2 Å, respectively. The structures reveal that Arg⁴⁷ plays a critical role in maintaining the structural features of the complementarity-determining region 2 (CDR2) loop and the putative positive ligand-interacting surface (PLIS), stabilizing conformations capable of ligand interaction. This is exemplified in the PS-bound structure, in which the CDR2 loop and PLIS drive critical interactions with PS via surfaces that are disrupted in the variant. Together with in vitro and in vivo characterization, our structural findings elucidate the molecular mechanism underlying loss of ligand binding, putative oligomerization, and functional activity of R47H TREM2. They also help unravel how decreased in vitro and in vivo stability of TREM2 contribute to loss of function in disease.

A Tale of Two Genes: Microglial Apoe and Trem2

Microglial cell function is implicated in the etiology of Alzheimer's disease by human genetics. In this issue of Immunity, Krasemann et al. (2017) describe a gene expression signature associated with an APOE- and TREM2-dependent response of microglia to brain tissue damage that accumulates in aging and disease, defining an axis that might be amenable to therapeutic targeting.

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.

Comparison of phosphodiesterase 10A, dopamine receptors D1 and D2 and dopamine transporter ligand binding in the striatum of the R6/2 and BACHD mouse models of Huntington's disease

BACKGROUND

Phosphodiesterase 10A (PDE10A) is expressed at high levels in the striatum and has been proposed both as a biomarker for Huntington's disease pathology and as a target for intervention.

OBJECTIVE

PDE10A radiotracers have been successfully used to measure changes in binding density in Huntington's disease patients, but little is known about PDE10A binding in mouse models that are used extensively to model pathology and test therapeutic interventions.

METHODS

Our study investigated changes in PDE10A binding using the selective tracer 3H-7980 at specific ages of two Huntington's disease transgenic mouse models: R6/2, a short-lived model carrying exon-1 of mutant HTT and BACHD, a longer-lived model carrying full-length mutant HTT. PDE10A binding was compared to binding of known markers of striatal atrophy in Huntington's disease, e.g. dopamine transporter (DAT) and dopamine receptors D1 and D2.

RESULTS

We found that in the R6/2 model at 6 weeks of age, mice showed high variability of binding, however binding of all ligands was significantly decreased at 8 and 12 weeks of age. In contrast, no changes were detectable in the BACHD model at 8, 10 or 12 month of age.

CONCLUSIONS

These findings suggest that radiotracer binding of PDE10A, DAT, D1 and D2 receptor in the R6/2 model may be a good indicator of striatal pathological changes that are observed in Huntington's disease patients, and that the first 12 months in the BACHD model may be more reflective of early stages of the disease.

The Huntington's disease mutation impairs Huntingtin's role in the transport of NF-κB from the synapse to the nucleus

Expansion of a polyglutamine (polyQ) tract in the Huntingtin (Htt) protein causes Huntington's disease (HD), a fatal inherited neurodegenerative disorder. Loss of the normal function of Htt is thought to be an important pathogenetic component of HD. However, the function of wild-type Htt is not well defined. Htt is thought to be a multifunctional protein that plays distinct roles in several biological processes, including synaptic transmission, intracellular transport and neuronal transcription. Here, we show with biochemical and live cell imaging studies that wild-type Htt stimulates the transport of nuclear factor κ light-chain-enhancer of activated B cells (NF-κB) out of dendritic spines (where NF-κB is activated by excitatory synaptic input) and supports a high level of active NF-κB in neuronal nuclei (where NF-κB stimulates the transcription of target genes). We show that this novel function of Htt is impaired by the polyQ expansion and thus may contribute to the etiology of HD.

Context-dependent regulation of NeuroD activity and protein accumulation

NeuroD/BETA2 (referred to as NeuroD hereafter) is a basic helix-loop-helix (bHLH) transcription factor that is required for the development and survival of a subset of neurons and pancreatic endocrine cells in mice. Gain-of-function analyses demonstrated that NeuroD can (i) convert epidermal fate into neuronal fate when overexpressed in Xenopus embryos, and (ii) activate the insulin promoter in pancreatic beta cell lines in response to glucose stimulation. In glucose-stimulated INS-1 pancreatic beta cells, mutations of S259, S266, and S274 to alanines inhibited the ability of NeuroD to activate the insulin promoter. Phosphorylation of those serine residues by ERK1/2 was required for NeuroD activity in that assay. To determine whether the same residues are implicated in the neurogenic activity of NeuroD, we mutated the conserved S259, S266, and S274 of Xenopus NeuroD to alanines (S259A, S266A, and S274A), and performed an ectopic neurogenesis assay in Xenopus embryos. In contrast to what has been observed in the pancreatic beta cell line, the S266A and S274A mutant forms of Xenopus NeuroD displayed significantly increased abilities to form ectopic neurons, while S259A had little effect. In addition, S266A and S274A of Xenopus NeuroD resulted in increased accumulation of protein in the injected embryos while the corresponding mutations on mouse NeuroD did not have the same effect in an insulinoma cell line. Our results demonstrate that the consequence of NeuroD protein modification is context-dependent at both the molecular and functional levels.

NeuroD1 in the endocrine pancreas: Localization and dual function as an activator and repressor

The basic helix-loop-helix transcription factor NeuroD1 regulates cell fate in the nervous system but previously has not been considered to function similarly in the endocrine pancreas due to its reported expression in all islet cell types in the newborn mouse. Because we found that NeuroD1 potently represses somatostatin expression in vitro, its pattern of expression was examined in both strains of mice in which lacZ has been introduced into the NeuroD1 locus by homologous recombination. Analysis of adult transgenic mice revealed that NeuroD1 is predominantly expressed in beta-cells and either absent or expressed below the limit of lacZ detection in mature alpha-, delta-, or PP cells. Consistent with a previous report, NeuroD1 colocalizes with glucagon as well as insulin in immature islets of the newborn mouse. However, no colocalization of NeuroD1with somatostatin was detected in the newborn. In vitro, ectopic expression of NeuroD1 in TRM-6/PDX-1, a human pancreatic delta-cell line, resulted in potent repression of somatostatin concomitant with induction of the beta-cell hormones insulin and islet amyloid polypeptide. Additionally, NeuroD1 induced expression of Nkx2.2, a transcription factor expressed in beta- but not delta-cells. Transfection studies using insulin and somatostatin promoters confirm the ability of NeuroD1 to act as both a transcriptional repressor and activator in the same cell, suggesting a more complex role for NeuroD1 in the establishment and/or maintenance of mature endocrine cells than has been recognized previously.

Stimulation of NeuroD activity by huntingtin and huntingtin-associated proteins HAP1 and MLK2

NeuroD (ND) is a basic helix-loop-helix transcription factor important for neuronal development and survival. By using a yeast two-hybrid screen, we identified two proteins that interact with ND, huntingtin-associated protein 1 (HAP1) and mixed-lineage kinase 2 (MLK2), both of which are known to interact with huntingtin (Htt). Htt is a ubiquitous protein important for neuronal transcription, development, and survival, and loss of its function has been implicated in the pathogenesis of Huntington's disease, a neurodegenerative disorder. However, the mechanism by which Htt exerts its neuron-specific function at the molecular level is unknown. Here we report that Htt interacts with ND via HAP1, and that MLK2 phosphorylates and stimulates the activity of ND. Furthermore, we show that Htt and HAP1 facilitate the activation of ND by MLK2. To our knowledge, ND is the first example of a neuron-specific transcription factor involved in neuronal development and survival whose activity is modulated by Htt. We propose that Htt, together with HAP1, may function as a scaffold for the activation of ND by MLK2.

The NeuroD1/BETA2 sequences essential for insulin gene transcription colocalize with those necessary for neurogenesis and p300/CREB binding protein binding

NeuroD1/BETA2 is a key regulator of pancreatic islet morphogenesis and insulin hormone gene transcription in islet beta cells. This factor also appears to be involved in neurogenic differentiation, because NeuroD1/BETA2 is able to induce premature differentiation of neuronal precursors and convert ectoderm into fully differentiated neurons upon ectopic expression in Xenopus embryos. We have identified amino acid sequences in mammalian and Xenopus NeuroD1/BETA2 that are necessary for insulin gene expression and ectopic neurogenesis. Our results indicate that evolutionarily conserved sequences spanning the basic helix-loop-helix (amino acids [aa] 100 to 155) and C-terminal (aa 156 to 355) regions are important for both of these processes. The transactivation domains (AD1, aa 189 to 299; AD2, aa 300 to 355) were within the carboxy-terminal region, as analyzed by using GAL4:NeuroD1/BETA2 chimeras. Selective activation of mammalian insulin gene enhancer-driven expression and ectopic neurogenesis in Xenopus embryos was regulated by two independent and separable domains of NeuroD1/BETA2, located between aa 156 to 251 and aa 252 to 355. GAL4:NeuroD1/BETA2 constructs spanning these sequences demonstrated that only aa 252 to 355 contained activation domain function, although both aa 156 to 251 and 300 to 355 were found to interact with the p300/CREB binding protein (CBP) coactivator. These results implicate p300/CBP in NeuroD1/BETA2 function and further suggest that comparable mechanisms are utilized to direct target gene transcription during differentiation and in adult islet beta cells.