Deep Sequence Analysis of Non-Small Cell Lung Cancer: Integrated Analysis of Gene Expression, Alternative Splicing, and Single Nucleotide Variations in Lung Adenocarcinomas with and without Oncogenic KRAS Mutations

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

A candidate regulatory variant at the TREM gene cluster associates with decreased Alzheimer's disease risk and increased TREML1 and TREM2 brain gene expression

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−3 and 4.6 × 10−2, respectively). Meta-analysis on expression quantitative trait locus results from three independent data sets (n = 1006) confirmed these associations (uncorrected P = 3.4 × 10−2 and 3.5 × 10−3, Bonferroni-corrected P = 6.7 × 10−2 and 7.1 × 10−3, 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

Concordant association of insulin degrading enzyme gene (IDE) variants with IDE mRNA, abeta, and alzheimer's disease.

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 Ass40 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

Genome-wide association study of corticobasal degeneration identifies risk variants shared with progressive supranuclear palsy

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 x 10(-12)),8p12 at lnc-KIF13B-1, a long non-coding RNA (rs643472;P = 3.41 x 10(-8)),and 2p22 at SOS1 (rs963731;P = 1.76 x 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 x 10(-7)) and MAPT H1c (17q21;rs242557;P = 7.91 x 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)

Association of MAPT haplotypes with Alzheimer’s disease risk and MAPT brain gene expression levels

Introduction: MAPT encodes for tau, the predominant component of neurofibrillary tangles that are neuropathological hallmarks of Alzheimer’s disease (AD). Genetic association of MAPT variants with late-onset AD (LOAD) risk has been inconsistent, although insufficient power and incomplete assessment of MAPT haplotypes may account for this. Methods: We examined the association of MAPT haplotypes with LOAD risk in more than 20,000 subjects (n-cases = 9,814, n-controls = 11,550) from Mayo Clinic (n-cases = 2,052, n-controls = 3,406) and the Alzheimer’s Disease Genetics Consortium (ADGC, n-cases = 7,762, n-controls = 8,144). We also assessed associations with brain MAPT gene expression levels measured in the cerebellum (n = 197) and temporal cortex (n = 202) of LOAD subjects. Six single nucleotide polymorphisms (SNPs) which tag MAPT haplotypes with frequencies greater than 1% were evaluated. Results: H2-haplotype tagging rs8070723-G allele associated with reduced risk of LOAD (odds ratio, OR = 0.90, 95% confidence interval, CI = 0.85-0.95, p = 5.2E-05) with consistent results in the Mayo (OR = 0.81, p = 7.0E-04) and ADGC (OR = 0.89, p = 1.26E-04) cohorts. rs3785883-A allele was also nominally significantly associated with LOAD risk (OR = 1.06, 95% CI = 1.01-1.13, p = 0.034). Haplotype analysis revealed significant global association with LOAD risk in the combined cohort (p = 0.033), with significant association of the H2 haplotype with reduced risk of LOAD as expected (p = 1.53E-04) and suggestive association with additional haplotypes. MAPT SNPs and haplotypes also associated with brain MAPT levels in the cerebellum and temporal cortex of AD subjects with the strongest associations observed for the H2 haplotype and reduced brain MAPT levels (β = -0.16 to -0.20, p = 1.0E-03 to 3.0E-03). Conclusions: These results confirm the previously reported MAPT H2 associations with LOAD risk in two large series, that this haplotype has the strongest effect on brain MAPT expression amongst those tested and identify additional haplotypes with suggestive associations, which require replication in independent series. These biologically congruent results provide compelling evidence to screen the MAPT region for regulatory variants which confer LOAD risk by influencing its brain gene expression

Divergent brain gene expression patterns associate with distinct cell-specific tau neuropathology traits in progressive supranuclear palsy.

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

TLR5 decoy receptor as a novel anti-amyloid therapeutic for Alzheimer\u27s disease.

There is considerable interest in harnessing innate immunity to treat Alzheimer\u27s disease (AD). Here, we explore whether a decoy receptor strategy using the ectodomain of select TLRs has therapeutic potential in AD. AAV-mediated expression of human TLR5 ectodomain (sTLR5) alone or fused to human IgG4 Fc (sTLR5Fc) results in robust attenuation of amyloid β (Aβ) accumulation in a mouse model of Alzheimer-type Aβ pathology. sTLR5Fc binds to oligomeric and fibrillar Aβ with high affinity, forms complexes with Aβ, and blocks Aβ toxicity. Oligomeric and fibrillar Aβ modulates flagellin-mediated activation of human TLR5 but does not, by itself, activate TLR5 signaling. Genetic analysis shows that rare protein coding variants in human TLR5 may be associated with a reduced risk of AD. Further, transcriptome analysis shows altered TLR gene expression in human AD. Collectively, our data suggest that TLR5 decoy receptor-based biologics represent a novel and safe Aβ-selective class of biotherapy in AD