Brain high-throughput multi-omics data reveal molecular heterogeneity in Alzheimer's disease.

Unbiased data-driven omic approaches are revealing the molecular heterogeneity of Alzheimer disease. Here, we used machine learning approaches to integrate high-throughput transcriptomic, proteomic, metabolomic, and lipidomic profiles with clinical and neuropathological data from multiple human AD cohorts. We discovered 4 unique multimodal molecular profiles, one of them showing signs of poor cognitive function, a faster pace of disease progression, shorter survival with the disease, severe neurodegeneration and astrogliosis, and reduced levels of metabolomic profiles. We found this molecular profile to be present in multiple affected cortical regions associated with higher Braak tau scores and significant dysregulation of synapse-related genes, endocytosis, phagosome, and mTOR signaling pathways altered in AD early and late stages. AD cross-omics data integration with transcriptomic data from an SNCA mouse model revealed an overlapping signature. Furthermore, we leveraged single-nuclei RNA-seq data to identify distinct cell-types that most likely mediate molecular profiles. Lastly, we identified that the multimodal clusters uncovered cerebrospinal fluid biomarkers poised to monitor AD progression and possibly cognition. Our cross-omics analyses provide novel critical molecular insights into AD.

Multi-ancestry GWAS reveals excitotoxicity associated with outcome after ischaemic stroke.

During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6 h of stroke onset and NIHSS at 24 h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24 h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke.

Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders.

Understanding the tissue-specific genetic controls of protein levels is essential to uncover mechanisms of post-transcriptional gene regulation. In this study, we generated a genomic atlas of protein levels in three tissues relevant to neurological disorders (brain, cerebrospinal fluid and plasma) by profiling thousands of proteins from participants with and without Alzheimer's disease. We identified 274, 127 and 32 protein quantitative trait loci (pQTLs) for cerebrospinal fluid, plasma and brain, respectively. cis-pQTLs were more likely to be tissue shared, but trans-pQTLs tended to be tissue specific. Between 48.0% and 76.6% of pQTLs did not co-localize with expression, splicing, DNA methylation or histone acetylation QTLs. Using Mendelian randomization, we nominated proteins implicated in neurological diseases, including Alzheimer's disease, Parkinson's disease and stroke. This first multi-tissue study will be instrumental to map signals from genome-wide association studies onto functional genes, to discover pathways and to identify drug targets for neurological diseases.

Meningeal lymphatics affect microglia responses and anti-Aß immunotherapy.

Alzheimer's disease (AD) is the most prevalent cause of dementia1. Although there is no effective treatment for AD, passive immunotherapy with monoclonal antibodies against amyloid beta (Aß) is a promising therapeutic strategy2,3. Meningeal lymphatic drainage has an important role in the accumulation of Aß in the brain4, but it is not known whether modulation of meningeal lymphatic function can influence the outcome of immunotherapy in AD. Here we show that ablation of meningeal lymphatic vessels in 5xFAD mice (a mouse model of amyloid deposition that expresses five mutations found in familial AD) worsened the outcome of mice treated with anti-Aß passive immunotherapy by exacerbating the deposition of Aß, microgliosis, neurovascular dysfunction, and behavioural deficits. By contrast, therapeutic delivery of vascular endothelial growth factor C improved clearance of Aß by monoclonal antibodies. Notably, there was a substantial overlap between the gene signature of microglia from 5xFAD mice with impaired meningeal lymphatic function and the transcriptional profile of activated microglia from the brains of individuals with AD. Overall, our data demonstrate that impaired meningeal lymphatic drainage exacerbates the microglial inflammatory response in AD and that enhancement of meningeal lymphatic function combined with immunotherapies could lead to better clinical outcomes.

African Americans Have Differences in CSF Soluble TREM2 and Associated Genetic Variants.

OBJECTIVE: To evaluate for racial differences in triggering receptor expressed on myeloid cells 2 (TREM2), a key immune mediator in Alzheimer disease, the levels of CSF soluble TREM2 (sTREM2), and the frequency of associated genetic variants were compared in groups of individuals who self-reported their race as African American (AA) or non-Hispanic White (NHW). METHODS: Community-dwelling older research participants underwent measurement of CSF sTREM2 concentrations and genetic analyses. RESULTS: The primary cohort included 91 AAs and 868 NHWs. CSF sTREM2 levels were lower in the AA compared with the NHW group (1,336 ± 470 vs 1,856 ± 624 pg/mL, p < 0.0001). AAs were more likely to carry TREM2 coding variants (15% vs 3%, p < 0.0001), which were associated with lower CSF sTREM2. AAs were less likely to carry the rs1582763 minor allele (8% vs 37%, p < 0.0001), located near MS4A4A, which was associated with higher CSF sTREM2. These findings were replicated in an independent cohort of 23 AAs and 917 NHWs: CSF sTREM2 levels were lower in the AA group (p = 0.03), AAs were more likely to carry coding TREM2 variants (22% vs 4%, p = 0.002), and AAs were less likely to carry the rs1582763 minor allele (16% vs 37%, p = 0.003). CONCLUSIONS: On average, AAs had lower CSF sTREM2 levels compared with NHWs, potentially because AAs are more likely to carry genetic variants associated with lower CSF sTREM2 levels. Importantly, CSF sTREM2 reflects TREM2-mediated microglial activity, a critical step in the immune response to amyloid plaques. These findings suggest that race may be associated with risk for genetic variants that influence Alzheimer disease-related inflammation.

Long runs of homozygosity are associated with Alzheimer's disease.

Long runs of homozygosity (ROH) are contiguous stretches of homozygous genotypes, which are a footprint of inbreeding and recessive inheritance. The presence of recessive loci is suggested for Alzheimer's disease (AD); however, their search has been poorly assessed to date. To investigate homozygosity in AD, here we performed a fine-scale ROH analysis using 10 independent cohorts of European ancestry (11,919 AD cases and 9181 controls.) We detected an increase of homozygosity in AD cases compared to controls [ßAVROH (CI 95%) = 0.070 (0.037-0.104); P = 3.91 × 10-5; ßFROH (CI95%) = 0.043 (0.009-0.076); P = 0.013]. ROHs increasing the risk of AD (OR > 1) were significantly overrepresented compared to ROHs increasing protection (p < 2.20 × 10-16). A significant ROH association with AD risk was detected upstream the HS3ST1 locus (chr4:11,189,482‒11,305,456), (ß (CI 95%) = 1.09 (0.48 ‒ 1.48), p value = 9.03 × 10-4), previously related to AD. Next, to search for recessive candidate variants in ROHs, we constructed a homozygosity map of inbred AD cases extracted from an outbred population and explored ROH regions in whole-exome sequencing data (N = 1449). We detected a candidate marker, rs117458494, mapped in the SPON1 locus, which has been previously associated with amyloid metabolism. Here, we provide a research framework to look for recessive variants in AD using outbred populations. Our results showed that AD cases have enriched homozygosity, suggesting that recessive effects may explain a proportion of AD heritability.

Quantitative endophenotypes as an alternative approach to understanding genetic risk in neurodegenerative diseases.

Endophenotypes, as measurable intermediate features of human diseases, reflect underlying molecular mechanisms. The use of quantitative endophenotypes in genetic studies has improved our understanding of pathophysiological changes associated with diseases. The main advantage of the quantitative endophenotypes approach to study human diseases over a classic case-control study design is the inferred biological context that can enable the development of effective disease-modifying treatments. Here, we summarize recent progress on biomarkers for neurodegenerative diseases, including cerebrospinal fluid and blood-based, neuroimaging, neuropathological, and clinical studies. This review focuses on how endophenotypic studies have successfully linked genetic modifiers to disease risk, disease onset, or progression rate and provided biological context to genes identified in genome-wide association studies. Finally, we review critical methodological considerations for implementing this approach and future directions.

Molecular pathways in patients with systemic lupus erythematosus revealed by gene-centred DNA sequencing.

OBJECTIVES: Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE. METHODS: We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS). RESULTS: We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes. CONCLUSIONS: Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.

Chi3l1/YKL-40 is controlled by the astrocyte circadian clock and regulates neuroinflammation and Alzheimer's disease pathogenesis.

Regulation of glial activation and neuroinflammation are critical factors in the pathogenesis of Alzheimer's disease (AD). YKL-40, a primarily astrocytic protein encoded by the gene Chi3l1, is a widely studied cerebrospinal fluid biomarker that increases with aging and early in AD. However, the function of Chi3l1/YKL-40 in AD is unknown. In a cohort of patients with AD, we observed that a variant in the human CHI3L1 gene, which results in decreased CSF YKL-40 expression, was associated with slower AD progression. At baseline, Chi3l1 deletion in mice had no effect on astrocyte activation while modestly promoting microglial activation. In a mouse APP/PS1 model of AD, Chi3l1 deletion decreased amyloid plaque burden and increased periplaque expression of the microglial lysosomal marker CD68, suggesting that Chi3l1 may suppress glial phagocytic activation and promote amyloid accumulation. Accordingly, Chi3l1 knockdown increased phagocytosis of zymosan particles and of ß-amyloid peptide in both astrocytes and microglia in vitro. We further observed that expression of Chi3l1 is regulated by the circadian clock, as deletion of the core clock proteins BMAL1 or CLOCK/NPAS2 strongly suppresses basal Chi3l1 expression, whereas deletion of the negative clock regulators PER1/PER2 increased Chi3l1 expression. Basal Chi3l1 mRNA was nonrhythmic because of a long mRNA half-life in astrocytes. However, inflammatory induction of Chi3l1 was gated by the clock. Our findings reveal Chi3l1/YKL-40 as a modulator of glial phagocytic activation and AD pathogenesis in both mice and humans and suggest that the astrocyte circadian clock regulates inflammatory Chi3l1 induction.

Impact of Permeable Membrane on the Hydrocyclone Separation Performance for Oily Water Treatment.

In the oil industry and academy, the treatment of water contaminated with oil using conventional hydrocyclones and membranes has been an alternative to meet the requirements established by environmental control agencies. However, such equipment is not fully efficient in the treatment of much diluted oily water, with both presenting restrictions in their performance. In this sense, the present work proposes to study the separation process of oily water using a new configuration of hydrocyclone, equipped with a porous ceramic membrane in the conical part's wall (filtering hydrocyclone). For the theoretical study, a Eulerian-Eulerian approach was applied to solve the mass and momentum conservation equations, and the turbulence model, using the computational fluid dynamics technique. The results of the velocity, pressure and volumetric fraction of the involved phases, and the separation performance of the hydrocyclone, are presented, analyzed, and compared with those obtained with a conventional hydrocyclone. The results confirmed the high potential of the proposed equipment to be used in the separation of the water and oil mixture.

Multi-ancestry genetic study in 5,876 patients identifies an association between excitotoxic genes and early outcomes after acute ischemic stroke.

During the first hours after stroke onset neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between NIH stroke scale (NIHSS) within six hours of stroke onset and NIHSS at 24h (ΔNIHSS). A total of 5,876 individuals from seven countries (Spain, Finland, Poland, United States, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of ΔNIHSS variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture than that of stroke risk. Seven loci (2p25.1, 2q31.2, 2q33.3, 4q34.3, 5q33.2, 6q26 and 7p21.1) were genome-wide significant and explained 2.1% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each loci. eQTL mapping and SMR indicate that ADAM23 (log Bayes Factor (LBF)=6.34) was driving the association for 2q33.3. Gene based analyses suggested that GRIA1 (LBF=5.26), which is predominantly expressed in brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated PARK2 (LBF=5.30) and ABCB5 (LBF=5.70) for the 6q26 and 7p21.1 loci. Human brain single nuclei RNA-seq indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23 , a pre-synaptic protein, and GRIA1 , a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provides the first evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischemic stroke. RESEARCH INTO CONTEXT: Evidence before this study: No previous genome-wide association studies have investigated the genetic architecture of early outcomes after ischemic stroke.Added Value of this study: This is the first study that investigated genetic influences on early outcomes after ischemic stroke using a genome-wide approach, revealing seven genome-wide significant loci. A unique aspect of this genetic study is the inclusion of all of the major ethnicities by recruiting from participants throughout the world. Most genetic studies to date have been limited to populations of European ancestry.Implications of all available evidence: The findings provide the first evidence that genes implicating excitotoxicity contribute to human acute ischemic stroke, and demonstrates proof of principle that GWAS of acute ischemic stroke patients can reveal mechanisms involved in ischemic brain injury.

Functional genomic analyses uncover APOE-mediated regulation of brain and cerebrospinal fluid beta-amyloid levels in Parkinson disease.

Alpha-synuclein is the main protein component of Lewy bodies, the pathological hallmark of Parkinson's disease. However, genetic modifiers of cerebrospinal fluid (CSF) alpha-synuclein levels remain unknown. The use of CSF levels of amyloid beta1-42, total tau, and phosphorylated tau181 as quantitative traits in genetic studies have provided novel insights into Alzheimer's disease pathophysiology. A systematic study of the genomic architecture of CSF biomarkers in Parkinson's disease has not yet been conducted. Here, genome-wide association studies of CSF biomarker levels in a cohort of individuals with Parkinson's disease and controls (N = 1960) were performed. PD cases exhibited significantly lower CSF biomarker levels compared to controls. A SNP, proxy for APOE ε4, was associated with CSF amyloid beta1-42 levels (effect = - 0.5, p = 9.2 × 10-19). No genome-wide loci associated with CSF alpha-synuclein, total tau, or phosphorylated tau181 levels were identified in PD cohorts. Polygenic risk score constructed using the latest Parkinson's disease risk meta-analysis were associated with Parkinson's disease status (p = 0.035) and the genomic architecture of CSF amyloid beta1-42 (R2 = 2.29%; p = 2.5 × 10-11). Individuals with higher polygenic risk scores for PD risk presented with lower CSF amyloid beta1-42 levels (p = 7.3 × 10-04). Two-sample Mendelian Randomization revealed that CSF amyloid beta1-42 plays a role in Parkinson's disease (p = 1.4 × 10-05) and age at onset (p = 7.6 × 10-06), an effect mainly mediated by variants in the APOE locus. In a subset of PD samples, the APOE ε4 allele was associated with significantly lower levels of CSF amyloid beta1-42 (p = 3.8 × 10-06), higher mean cortical binding potentials (p = 5.8 × 10-08), and higher Braak amyloid beta score (p = 4.4 × 10-04). Together these results from high-throughput and hypothesis-free approaches converge on a genetic link between Parkinson's disease, CSF amyloid beta1-42, and APOE.

A new domestic cat genome assembly based on long sequence reads empowers feline genomic medicine and identifies a novel gene for dwarfism.

The domestic cat (Felis catus) numbers over 94 million in the USA alone, occupies households as a companion animal, and, like humans, suffers from cancer and common and rare diseases. However, genome-wide sequence variant information is limited for this species. To empower trait analyses, a new cat genome reference assembly was developed from PacBio long sequence reads that significantly improve sequence representation and assembly contiguity. The whole genome sequences of 54 domestic cats were aligned to the reference to identify single nucleotide variants (SNVs) and structural variants (SVs). Across all cats, 16 SNVs predicted to have deleterious impacts and in a singleton state were identified as high priority candidates for causative mutations. One candidate was a stop gain in the tumor suppressor FBXW7. The SNV is found in cats segregating for feline mediastinal lymphoma and is a candidate for inherited cancer susceptibility. SV analysis revealed a complex deletion coupled with a nearby potential duplication event that was shared privately across three unrelated cats with dwarfism and is found within a known dwarfism associated region on cat chromosome B1. This SV interrupted UDP-glucose 6-dehydrogenase (UGDH), a gene involved in the biosynthesis of glycosaminoglycans. Importantly, UGDH has not yet been associated with human dwarfism and should be screened in undiagnosed patients. The new high-quality cat genome reference and the compilation of sequence variation demonstrate the importance of these resources when searching for disease causative alleles in the domestic cat and for identification of feline biomedical models.

Examination of the Effect of Rare Variants in TREM2, ABI3, and PLCG2 in LOAD Through Multiple Phenotypes.

BACKGROUND: Rare variants in PLCG2 (p.P522R), ABI3 (p.S209F), and TREM2 (p.R47H, p.R62H) have been associated with late onset Alzheimer's disease (LOAD) risk in Caucasians. After the initial report, several studies have found positive results in cohorts of different ethnic background and with different phenotype. OBJECTIVE: In this study, we aim to evaluate the association of rare coding variants in PLCG2, ABI3, and TREM2 with LOAD risk and their effect at different time points of the disease. METHODS: We used a European American cohort to assess the association of the variants prior onset (using CSF Aß42, tau, and pTau levels, and amyloid imaging as endophenotypes) and after onset (measured as rate of memory decline). RESULTS: We confirm the association with LOAD risk of TREM2 p.R47H, p.R62H and ABI3 p.S209F variants, and the protective effect of PLCG2 p.P522R. In addition, ABI3 and TREM2 gene-sets showed significant association with LOAD risk. TREM2 p.R47H and PLCG2 p.P522R variants were also statistically associated with increase of amyloid imaging and AD progression, respectively. We did not observe any association of ABI3 p.S209F with any of the other AD endophenotypes. CONCLUSION: The results of this study highlight the importance of including biomarkers and alternative phenotypes to better understand the role of novel candidate genes with the disease.

The TMEM106B FTLD-protective variant, rs1990621, is also associated with increased neuronal proportion.

Apart from amyloid ß deposition and tau neurofibrillary tangles, Alzheimer's disease (AD) is a neurodegenerative disorder characterized by neuronal loss and astrocytosis in the cerebral cortex. The goal of this study is to investigate genetic factors associated with the neuronal proportion in health and disease. To identify cell-autonomous genetic variants associated with neuronal proportion in cortical tissues, we inferred cellular population structure from bulk RNA-Seq derived from 1536 individuals. We identified the variant rs1990621 located in the TMEM106B gene region as significantly associated with neuronal proportion (p value = 6.40 × 10-07) and replicated this finding in an independent dataset (p value = 7.41 × 10-04) surpassing the genome-wide threshold in the meta-analysis (p value = 9.42 × 10-09). This variant is in high LD with the TMEM106B non-synonymous variant p.T185S (rs3173615; r2 = 0.98) which was previously identified as a protective variant for frontotemporal lobar degeneration (FTLD). We stratified the samples by disease status, and discovered that this variant modulates neuronal proportion not only in AD cases, but also several neurodegenerative diseases and in elderly cognitively healthy controls. Furthermore, we did not find a significant association in younger controls or schizophrenia patients, suggesting that this variant might increase neuronal survival or confer resilience to the neurodegenerative process. The single variant and gene-based analyses also identified an overall genetic association between neuronal proportion, AD and FTLD risk. These results suggest that common pathways are implicated in these neurodegenerative diseases, that implicate neuronal survival. In summary, we identified a protective variant in the TMEM106B gene that may have a neuronal protection effect against general aging, independent of disease status, which could help elucidate the relationship between aging and neuronal survival in the presence or absence of neurodegenerative disorders. Our findings suggest that TMEM106B could be a potential target for neuronal protection therapies to ameliorate cognitive and functional deficits.

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations.

Parietal cortex RNA-sequencing (RNA-seq) data were generated from individuals with and without Alzheimer disease (AD; ncontrol = 13; nAD = 83) from the Knight Alzheimer Disease Research Center (Knight ADRC). Using this and an independent (Mount Sinai Brain Bank (MSBB)) AD RNA-seq dataset, cortical circular RNA (circRNA) expression was quantified in the context of AD. Significant associations were identified between circRNA expression and AD diagnosis, clinical dementia severity and neuropathological severity. It was demonstrated that most circRNA-AD associations are independent of changes in cognate linear messenger RNA expression or estimated brain cell-type proportions. Evidence was provided for circRNA expression changes occurring early in presymptomatic AD and in autosomal dominant AD. It was also observed that AD-associated circRNAs co-expressed with known AD genes. Finally, potential microRNA-binding sites were identified in AD-associated circRNAs for miRNAs predicted to target AD genes. Together, these results highlight the importance of analyzing non-linear RNAs and support future studies exploring the potential roles of circRNAs in AD pathogenesis.

The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer's disease risk.

Soluble triggering receptor expressed on myeloid cells 2 (sTREM2) in cerebrospinal fluid (CSF) has been associated with Alzheimer's disease (AD). TREM2 plays a critical role in microglial activation, survival, and phagocytosis; however, the pathophysiological role of sTREM2 in AD is not well understood. Understanding the role of sTREM2 in AD may reveal new pathological mechanisms and lead to the identification of therapeutic targets. We performed a genome-wide association study (GWAS) to identify genetic modifiers of CSF sTREM2 obtained from the Alzheimer's Disease Neuroimaging Initiative. Common variants in the membrane-spanning 4-domains subfamily A (MS4A) gene region were associated with CSF sTREM2 concentrations (rs1582763; P = 1.15 × 10-15); this was replicated in independent datasets. The variants associated with increased CSF sTREM2 concentrations were associated with reduced AD risk and delayed age at onset of disease. The single-nucleotide polymorphism rs1582763 modified expression of the MS4A4A and MS4A6A genes in multiple tissues, suggesting that one or both of these genes are important for modulating sTREM2 production. Using human macrophages as a proxy for microglia, we found that MS4A4A and TREM2 colocalized on lipid rafts at the plasma membrane, that sTREM2 increased with MS4A4A overexpression, and that silencing of MS4A4A reduced sTREM2 production. These genetic, molecular, and cellular findings suggest that MS4A4A modulates sTREM2. These findings also provide a mechanistic explanation for the original GWAS signal in the MS4A locus for AD risk and indicate that TREM2 may be involved in AD pathogenesis not only in TREM2 risk-variant carriers but also in those with sporadic disease.

Overlap in the Genetic Architecture of Stroke Risk, Early Neurological Changes, and Cardiovascular Risk Factors.

Background and Purpose- The genetic relationships between stroke risk, stroke severity, and early neurological changes are complex and not completely understood. Genetic studies have identified 32 all stroke risk loci. Polygenic risk scores can be used to compare the genetic architecture of related traits. In this study, we compare the genetic architecture of stroke risk, stroke severity, and early neurological changes with that of 2 stroke risk factors: type 2 diabetes mellitus (T2DM) and hypertension. Methods- We assessed the degree of overlap in the genetic architecture of stroke risk, T2DM, hypertension, and 2 acute stroke phenotypes based on the National Institutes of Health Stroke Scale (NIHSS), which ranges from 0 for no stroke symptoms to 21 to 42 for a severe stroke: baseline (within 6 hours after onset) and change in NIHSS (ΔNIHSS=NIHSS at baseline-NIHSS at 24 hours). This was done by (1) single-nucleotide polymorphism by single-nucleotide polymorphism comparison, (2) weighted polygenic risk scores with sentinel variants, and (3) whole-genome polygenic risk scores using multiple P thresholds. Results- We found evidence of genetic architecture overlap between stroke risk and T2DM ( P=2.53×10-169), hypertension ( P=3.93×10-04), and baseline NIHSS ( P=0.03). However, there was no evidence of overlap between ΔNIHSS and stroke risk, T2DM, or hypertension. Conclusions- The genetic architecture of stroke risk is correlated with that of T2DM, hypertension, and initial stroke severity (NIHSS within 6 hours of stroke onset). However, the genetic architecture of early neurological change after stroke (ΔNIHSS) is not correlated with that of ischemic stroke risk, T2DM, or hypertension. Thus, stroke risk and early neurological change after stroke have distinct genetic architectures.

TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers.

BACKGROUND: Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript. METHODS: We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million. RESULTS: The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007). CONCLUSION: Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

A rare regulatory variant in the MEF2D gene affects gene regulation and splicing and is associated with a SLE sub-phenotype in Swedish cohorts.

Systemic lupus erythematosus (SLE) is an autoimmune disorder with heterogeneous clinical presentation and complex etiology involving the interplay between genetic, epigenetic, environmental and hormonal factors. Many common SNPs identified by genome wide-association studies (GWAS) explain only a small part of the disease heritability suggesting the contribution from rare genetic variants, undetectable in GWAS, and complex epistatic interactions. Using targeted re-sequencing of coding and conserved regulatory regions within and around 215 candidate genes selected on the basis of their known role in autoimmunity and genes associated with canine immune-mediated diseases, we identified a rare regulatory variant rs200395694:G > T located in intron 4 of the MEF2D gene encoding the myocyte-specific enhancer factor 2D transcription factor and associated with SLE in Swedish cohorts (504 SLE patients and 839 healthy controls, p = 0.014, CI = 1.1-10). Fisher's exact test revealed an association between the genetic variant and a triad of disease manifestations including Raynaud, anti-U1-ribonucleoprotein (anti-RNP), and anti-Smith (anti-Sm) antibodies (p = 0.00037) among the patients. The DNA-binding activity of the allele was further studied by EMSA, reporter assays, and minigenes. The region has properties of an active cell-specific enhancer, differentially affected by the alleles of rs200395694:G > T. In addition, the risk allele exerts an inhibitory effect on the splicing of the alternative tissue-specific isoform, and thus may modify the target gene set regulated by this isoform. These findings emphasize the potential of dissecting traits of complex diseases and correlating them with rare risk alleles with strong biological effects.