Publications by authors named "Curtis S Younkin"

Background: Rare coding variants ABI3_rs616338-T and PLCG2_rs72824905-G were identified as risk or protective factors, respectively, for Alzheimer's disease (AD).

Methods: We tested the association of these variants with five neurodegenerative diseases in Caucasian case-control cohorts: 2742 AD, 231 progressive supranuclear palsy (PSP), 838 Parkinson's disease (PD), 306 dementia with Lewy bodies (DLB) and 150 multiple system atrophy (MSA) vs. 3351 controls; and in an African-American AD case-control cohort (181 AD, 331 controls). 1479 AD and 1491 controls were non-overlapping with a prior report.

Results: Using Fisher's exact test, there was significant association of both ABI3_rs616338-T (OR = 1.41, p = 0.044) and PLCG2_rs72824905-G (OR = 0.56, p = 0.008) with AD. These OR estimates were maintained in the non-overlapping replication AD-control analysis, albeit at reduced significance (ABI3_rs616338-T OR = 1.44, p = 0.12; PLCG2_rs72824905-G OR = 0.66, p = 0.19). None of the other cohorts showed significant associations that were concordant with those for AD, although the DLB cohort had suggestive findings (Fisher's test: ABI3_rs616338-T OR = 1.79, p = 0.097; PLCG2_rs72824905-G OR = 0.32, p = 0.124). PLCG2_rs72824905-G showed suggestive association with pathologically-confirmed MSA (OR = 2.39, p = 0.050) and PSP (OR = 1.97, p = 0.061), although in the opposite direction of that for AD. We assessed RNA sequencing data from 238 temporal cortex (TCX) and 224 cerebellum (CER) samples from AD, PSP and control patients and identified co-expression networks, enriched in microglial genes and immune response GO terms, and which harbor PLCG2 and/or ABI3. These networks had higher expression in AD, but not in PSP TCX, compared to controls. This expression association did not survive adjustment for brain cell type population changes.

Conclusions: We validated the associations previously reported with ABI3_rs616338-T and PLCG2_rs72824905-G in a Caucasian AD case-control cohort, and observed a similar direction of effect in DLB. Conversely, PLCG2_rs72824905-G showed suggestive associations with PSP and MSA in the opposite direction. We identified microglial gene-enriched co-expression networks with significantly higher levels in AD TCX, but not in PSP, a primary tauopathy. This co-expression network association appears to be driven by microglial cell population changes in a brain region affected by AD pathology. Although these findings require replication in larger cohorts, they suggest distinct effects of the microglial genes, ABI3 and PLCG2 in neurodegenerative diseases that harbor significant vs. low/no amyloid ß pathology.

Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by tau pathology in neurons and glial cells. Transcriptional regulation has been implicated as a potential mechanism in conferring disease risk and neuropathology for some PSP genetic risk variants. However, the role of transcriptional changes as potential drivers of distinct cell-specific tau lesions has not been explored. In this study, we integrated brain gene expression measurements, quantitative neuropathology traits and genome-wide genotypes from 268 autopsy-confirmed PSP patients to identify transcriptional associations with unique cell-specific tau pathologies. We provide individual transcript and transcriptional network associations for quantitative oligodendroglial (coiled bodies = CB), neuronal (neurofibrillary tangles = NFT), astrocytic (tufted astrocytes = TA) tau pathology, and tau threads and genomic annotations of these findings. We identified divergent patterns of transcriptional associations for the distinct tau lesions, with the neuronal and astrocytic neuropathologies being the most different. We determined that NFT are positively associated with a brain co-expression network enriched for synaptic and PSP candidate risk genes, whereas TA are positively associated with a microglial gene-enriched immune network. In contrast, TA is negatively associated with synaptic and NFT with immune system transcripts. Our findings have implications for the diverse molecular mechanisms that underlie cell-specific vulnerability and disease risk in PSP.

Introduction: Comparative transcriptome analyses in Alzheimer's disease (AD) and other neurodegenerative proteinopathies can uncover both shared and distinct disease pathways.

Methods: We analyzed 940 brain transcriptomes including patients with AD, progressive supranuclear palsy (PSP; a primary tauopathy), and control subjects.

Results: We identified transcriptional coexpression networks implicated in myelination, which were lower in PSP temporal cortex (TCX) compared with AD. Some of these associations were retained even after adjustments for brain cell population changes. These TCX myelination network structures were preserved in cerebellum but they were not differentially expressed in cerebellum between AD and PSP. Myelination networks were downregulated in both AD and PSP, when compared with control TCX samples.

Discussion: Downregulation of myelination networks may underlie both PSP and AD pathophysiology, but may be more pronounced in PSP. These data also highlight conservation of transcriptional networks across brain regions and the influence of cell type changes on these networks.

Introduction: We hypothesized that common Alzheimer's disease (AD)-associated variants within the triggering receptor expressed on myeloid (TREM) gene cluster influence disease through gene expression.

Methods: Expression microarrays on temporal cortex and cerebellum from ∼400 neuropathologically diagnosed subjects and two independent RNAseq replication cohorts were used for expression quantitative trait locus analysis.

Results: A variant within a DNase hypersensitive site 5' of TREM2, rs9357347-C, associates with reduced AD risk and increased TREML1 and TREM2 levels (uncorrected P = 6.3 × 10 and 4.6 × 10, respectively). Meta-analysis on expression quantitative trait locus results from three independent data sets (n = 1006) confirmed these associations (uncorrected P = 3.4 × 10 and 3.5 × 10, Bonferroni-corrected P = 6.7 × 10 and 7.1 × 10, respectively).

Discussion: Our findings point to rs9357347 as a functional regulatory variant that contributes to a protective effect observed at the TREM locus in the International Genomics of Alzheimer's Project genome-wide association study meta-analysis and suggest concomitant increase in TREML1 and TREM2 brain levels as a potential mechanism for protection from AD.

Previous genome-wide association studies (GWAS), conducted by our group and others, have identified loci that harbor risk variants for neurodegenerative diseases, including Alzheimer's disease (AD). Human disease variants are enriched for polymorphisms that affect gene expression, including some that are known to associate with expression changes in the brain. Postulating that many variants confer risk to neurodegenerative disease via transcriptional regulatory mechanisms, we have analyzed gene expression levels in the brain tissue of subjects with AD and related diseases. Herein, we describe our collective datasets comprised of GWAS data from 2,099 subjects; microarray gene expression data from 773 brain samples, 186 of which also have RNAseq; and an independent cohort of 556 brain samples with RNAseq. We expect that these datasets, which are available to all qualified researchers, will enable investigators to explore and identify transcriptional mechanisms contributing to neurodegenerative diseases.

To determine the effects of single nucleotide polymorphisms (SNPs) identified in a genome-wide association study of progressive supranuclear palsy (PSP), we tested their association with brain gene expression, CpG methylation and neuropathology. In 175 autopsied PSP subjects, we performed associations between seven PSP risk variants and temporal cortex levels of 20 genes in-cis, within ±100 kb. Methylation measures were collected using reduced representation bisulfite sequencing in 43 PSP brains. To determine whether SNP/expression associations are due to epigenetic modifications, CpG methylation levels of associated genes were tested against relevant variants. Quantitative neuropathology endophenotypes were tested for SNP associations in 422 PSP subjects. Brain levels of LRRC37A4 and ARL17B were associated with rs8070723; MOBP with rs1768208 and both ARL17A and ARL17B with rs242557. Expression associations for LRRC37A4 and MOBP were available in an additional 100 PSP subjects. Meta-analysis revealed highly significant associations for PSP risk alleles of rs8070723 and rs1768208 with higher LRRC37A4 and MOBP brain levels, respectively. Methylation levels of one CpG in the 3' region of ARL17B associated with rs242557 and rs8070723. Additionally, methylation levels of an intronic ARL17A CpG associated with rs242557 and that of an intronic MOBP CpG with rs1768208. MAPT and MOBP region risk alleles also associated with higher levels of neuropathology. Strongest associations were observed for rs242557/coiled bodies and tufted astrocytes; and for rs1768208/coiled bodies and tau threads. These findings suggest that PSP variants at MAPT and MOBP loci may confer PSP risk via influencing gene expression and tau neuropathology. MOBP, LRRC37A4, ARL17A and ARL17B warrant further assessment as candidate PSP risk genes. Our findings have implications for the mechanism of action of variants at some of the top PSP risk loci.

Corticobasal degeneration (CBD) is a neurodegenerative disorder affecting movement and cognition, definitively diagnosed only at autopsy. Here, we conduct a genome-wide association study (GWAS) in CBD cases (n=152) and 3,311 controls, and 67 CBD cases and 439 controls in a replication stage. Associations with meta-analysis were 17q21 at MAPT (P=1.42 × 10(-12)), 8p12 at lnc-KIF13B-1, a long non-coding RNA (rs643472; P=3.41 × 10(-8)), and 2p22 at SOS1 (rs963731; P=1.76 × 10(-7)). Testing for association of CBD with top progressive supranuclear palsy (PSP) GWAS single-nucleotide polymorphisms (SNPs) identified associations at MOBP (3p22; rs1768208; P=2.07 × 10(-7)) and MAPT H1c (17q21; rs242557; P=7.91 × 10(-6)). We previously reported SNP/transcript level associations with rs8070723/MAPT, rs242557/MAPT, and rs1768208/MOBP and herein identified association with rs963731/SOS1. We identify new CBD susceptibility loci and show that CBD and PSP share a genetic risk factor other than MAPT at 3p22 MOBP (myelin-associated oligodendrocyte basic protein).

Objective: Recent genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) identified 9 novel risk loci. Discovery of functional variants within genes at these loci is required to confirm their role in Alzheimer disease (AD). Single nucleotide polymorphisms that influence gene expression (eSNPs) constitute an important class of functional variants. We therefore investigated the influence of the novel LOAD risk loci on human brain gene expression.

Methods: We measured gene expression levels in the cerebellum and temporal cortex of autopsied AD subjects and those with other brain pathologies (∼400 total subjects). To determine whether any of the novel LOAD risk variants are eSNPs, we tested their cis-association with expression of 6 nearby LOAD candidate genes detectable in human brain (ABCA7, BIN1, CLU, MS4A4A, MS4A6A, PICALM) and an additional 13 genes ±100 kb of these SNPs. To identify additional eSNPs that influence brain gene expression levels of the novel candidate LOAD genes, we identified SNPs ±100 kb of their location and tested for cis-associations.

Results: CLU rs11136000 (p = 7.81 × 10(-4)) and MS4A4A rs2304933/rs2304935 (p = 1.48 × 10(-4)-1.86 × 10(-4)) significantly influence temporal cortex expression levels of these genes. The LOAD-protective CLU and risky MS4A4A locus alleles associate with higher brain levels of these genes. There are other cis-variants that significantly influence brain expression of CLU and ABCA7 (p = 4.01 × 10(-5)-9.09 × 10(-9)), some of which also associate with AD risk (p = 2.64 × 10(-2)-6.25 × 10(-5)).

Conclusions: CLU and MS4A4A eSNPs may at least partly explain the LOAD risk association at these loci. CLU and ABCA7 may harbor additional strong eSNPs. These results have implications in the search for functional variants at the novel LOAD risk loci.

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n=197, temporal cortex n=202) and with other brain pathologies (non-AD, cerebellar n=177, temporal cortex n=197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ± 100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non-ADs (q<0.05, p=7.70 × 10(-5)-1.67 × 10(-82)). Of these, 2,089 were also significant in the temporal cortex (p=1.85 × 10(-5)-1.70 × 10(-141)). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10(-6)). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9-3.3 fold enrichment (p<10(-6)) of significant cisSNPs with suggestive AD-risk association (p<10(-3)) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non-CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

KRAS mutations are highly prevalent in non-small cell lung cancer (NSCLC), and tumors harboring these mutations tend to be aggressive and resistant to chemotherapy. We used next-generation sequencing technology to identify pathways that are specifically altered in lung tumors harboring a KRAS mutation. Paired-end RNA-sequencing of 15 primary lung adenocarcinoma tumors (8 harboring mutant KRAS and 7 with wild-type KRAS) were performed. Sequences were mapped to the human genome, and genomic features, including differentially expressed genes, alternate splicing isoforms and single nucleotide variants, were determined for tumors with and without KRAS mutation using a variety of computational methods. Network analysis was carried out on genes showing differential expression (374 genes), alternate splicing (259 genes), and SNV-related changes (65 genes) in NSCLC tumors harboring a KRAS mutation. Genes exhibiting two or more connections from the lung adenocarcinoma network were used to carry out integrated pathway analysis. The most significant signaling pathways identified through this analysis were the NFκB, ERK1/2, and AKT pathways. A 27 gene mutant KRAS-specific sub network was extracted based on gene-gene connections from the integrated network, and interrogated for druggable targets. Our results confirm previous evidence that mutant KRAS tumors exhibit activated NFκB, ERK1/2, and AKT pathways and may be preferentially sensitive to target therapeutics toward these pathways. In addition, our analysis indicates novel, previously unappreciated links between mutant KRAS and the TNFR and PPARγ signaling pathways, suggesting that targeted PPARγ antagonists and TNFR inhibitors may be useful therapeutic strategies for treatment of mutant KRAS lung tumors. Our study is the first to integrate genomic features from RNA-Seq data from NSCLC and to define a first draft genomic landscape model that is unique to tumors with oncogenic KRAS mutations.

Background: Glutathione S-transferase omega-1 and 2 genes (GSTO1, GSTO2), residing within an Alzheimer and Parkinson disease (AD and PD) linkage region, have diverse functions including mitigation of oxidative stress and may underlie the pathophysiology of both diseases. GSTO polymorphisms were previously reported to associate with risk and age-at-onset of these diseases, although inconsistent follow-up study designs make interpretation of results difficult. We assessed two previously reported SNPs, GSTO1 rs4925 and GSTO2 rs156697, in AD (3,493 ADs vs. 4,617 controls) and PD (678 PDs vs. 712 controls) for association with disease risk (case-controls), age-at-diagnosis (cases) and brain gene expression levels (autopsied subjects).

Results: We found that rs156697 minor allele associates with significantly increased risk (odds ratio = 1.14, p = 0.038) in the older ADs with age-at-diagnosis > 80 years. The minor allele of GSTO1 rs4925 associates with decreased risk in familial PD (odds ratio = 0.78, p = 0.034). There was no other association with disease risk or age-at-diagnosis. The minor alleles of both GSTO SNPs associate with lower brain levels of GSTO2 (p = 4.7 × 10-11-1.9 × 10-27), but not GSTO1. Pathway analysis of significant genes in our brain expression GWAS, identified significant enrichment for glutathione metabolism genes (p = 0.003).

Conclusion: These results suggest that GSTO locus variants may lower brain GSTO2 levels and consequently confer AD risk in older age. Other glutathione metabolism genes should be assessed for their effects on AD and other chronic, neurologic diseases.

Background: The insulin-degrading enzyme gene (IDE) is a strong functional and positional candidate for late onset Alzheimer's disease (LOAD).

Methodology/principal Findings: We examined conserved regions of IDE and its 10 kb flanks in 269 AD cases and 252 controls thereby identifying 17 putative functional polymorphisms. These variants formed eleven haplotypes that were tagged with ten variants. Four of these showed significant association with IDE transcript levels in samples from 194 LOAD cerebella. The strongest, rs6583817, which has not previously been reported, showed unequivocal association (p = 1.5x10(-8), fold-increase = 2.12,); the eleven haplotypes were also significantly associated with transcript levels (global p = 0.003). Using an in vitro dual luciferase reporter assay, we found that rs6583817 increases reporter gene expression in Be(2)-C (p = 0.006) and HepG2 (p = 0.02) cell lines. Furthermore, using data from a recent genome-wide association study of two Croatian isolated populations (n = 1,879), we identified a proxy for rs6583817 that associated significantly with decreased plasma Abeta40 levels (ss = -0.124, p = 0.011) and total measured plasma Abeta levels (b = -0.130, p = 0.009). Finally, rs6583817 was associated with decreased risk of LOAD in 3,891 AD cases and 3,605 controls. (OR = 0.87, p = 0.03), and the eleven IDE haplotypes (global p = 0.02) also showed significant association.

Conclusions: Thus, a previously unreported variant unequivocally associated with increased IDE expression was also associated with reduced plasma Abeta40 and decreased LOAD susceptibility. Genetic association between LOAD and IDE has been difficult to replicate. Our findings suggest that targeted testing of expression SNPs (eSNPs) strongly associated with altered transcript levels in autopsy brain samples may be a powerful way to identify genetic associations with LOAD that would otherwise be difficult to detect.

By analyzing late-onset Alzheimer's disease (LOAD) in a genome-wide association study (313,504 SNPs, three series, 844 cases and 1,255 controls) and evaluating the 25 SNPs with the most significant allelic association in four additional series (1,547 cases and 1,209 controls), we identified a SNP (rs5984894) on Xq21.3 in PCDH11X that is strongly associated with LOAD in individuals of European descent from the United States. Analysis of rs5984894 by multivariable logistic regression adjusted for sex gave global P values of 5.7 x 10(-5) in stage 1, 4.8 x 10(-6) in stage 2 and 3.9 x 10(-12) in the combined data. Odds ratios were 1.75 (95% CI = 1.42-2.16) for female homozygotes (P = 2.0 x 10(-7)) and 1.26 (95% CI = 1.05-1.51) for female heterozygotes (P = 0.01) compared to female noncarriers. For male hemizygotes (P = 0.07) compared to male noncarriers, the odds ratio was 1.18 (95% CI = 0.99-1.41).