At least 90% of human genes are alternatively spliced. Alternative splicing has an important function regulating gene expression and miss-splicing can contribute to risk for human diseases, including Alzheimer’s disease (AD). We developed a splicing decision model as a molecular mechanism to identify functional exon skipping events and genetic variation affecting alternative splicing on a genome-wide scale by integrating genomics, transcriptomics, and neuroimaging data in a systems biology approach. In this study, we analyzed RNA-Seq data of hippocampus brain tissue from Alzheimer’s disease (AD; n = 24) and cognitively normal elderly controls (CN; n = 50) and identified three exon skipping events in two genes (RELN and NOS1) as significantly associated with AD (corrected p-value < 0.05 and fold change > 1.5). Next, we identified single-nucleotide polymorphisms (SNPs) affecting exon skipping events using the splicing decision model and then performed an association analysis of SNPs potentially affecting three exon skipping events with a global cortical measure of amyloid-β deposition measured by [18F] Florbetapir position emission tomography (PET) scan as an AD-related quantitative phenotype. A whole-brain voxel-based analysis was also performed. Two exons in RELN and one exon in NOS1 showed significantly lower expression levels in the AD participants compared to CN participants, suggesting that the exons tend to be skipped more in AD. We also showed the loss of the core protein structure due to the skipped exons using the protein 3D structure analysis. The targeted SNP-based association analysis identified one intronic SNP (rs362771) adjacent to the skipped exon 24 in RELN as significantly associated with cortical amyloid-β levels (corrected p-value < 0.05). This SNP is within the splicing regulatory element, i.e., intronic splicing enhancer. The minor allele of rs362771 conferred decreases in cortical amyloid-β levels in the right temporal and bilateral parietal lobes. Our results suggest that exon skipping events and splicing-affecting SNPs in the human hippocampus may contribute to AD pathogenesis. Integration of multiple omics and neuroimaging data provides insights into possible mechanisms underlying AD pathophysiology through exon skipping and may help identify novel therapeutic targets.

TREM2 is among the most well-known Alzheimer’s disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.

Clusterin (CLU) is one of the risk genes most associated with late onset Alzheimer’s disease (AD), and several genetic variants in CLU are associated with AD risk. However, the functional role of known AD risk genetic variants in CLU has been little explored. We investigated the effect of an AD risk variant (rs7982) in the 5th exon of CLU on alternative splicing by using an integrative approach of brain-tissue-based RNA-Seq and whole genome sequencing data from Accelerating Medicines Partnership—Alzheimer’s Disease (AMP-AD). RNA-Seq data were generated from three regions in the temporal lobe of the brain—the temporal cortex, superior temporal gyrus, and parahippocampal gyrus. The rs7982 was significantly associated with intron retention (IR) of the 5th exon of CLU; as the number of alternative alleles (G) increased, the IR rates decreased more significantly in females than in males. Our results suggest a sex-dependent role of rs7982 in AD pathogenesis via splicing regulation.

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e−06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

The Cancer Genome Atlas (TCGA) project is a public resource that provides transcriptomic, DNA sequence, methylation, and clinical data for 33 cancer types. Transforming the large size and high complexity of TCGA cancer genome data into integrated knowledge can be useful to promote cancer research. Alternative splicing (AS) is a key regulatory mechanism of genes in human cancer development and in the interaction with epigenetic factors. Therefore, AS-guided integration of existing TCGA data sets will make it easier to gain insight into the genetic architecture of cancer risk and related outcomes. There are already existing tools analyzing and visualizing alternative mRNA splicing patterns for large-scale RNA-seq experiments. However, these existing web-based tools are limited to the analysis of individual TCGA data sets at a time, such as only transcriptomic information. We implemented CAS-viewer (integrative analysis of Cancer genome data based on Alternative Splicing), a web-based tool leveraging multi-cancer omics data from TCGA. It illustrates alternative mRNA splicing patterns along with methylation, miRNAs, and SNPs, and then provides an analysis tool to link differential transcript expression ratio to methylation, miRNA, and splicing regulatory elements for 33 cancer types. Moreover, one can analyze AS patterns with clinical data to identify potential transcripts associated with different survival outcome for each cancer. CAS-viewer is a web-based application for transcript isoform-driven integration of multi-omics data in multiple cancer types and will aid in the visualization and possible discovery of biomarkers for cancer by integrating multi-omics data from TCGA.

The molecular basis of increased hemoglobin in Andean Aymara highlanders is unknown. We conducted an integrative analysis of whole-genome-sequencing and granulocytes transcriptomics from Aymara and Europeans in Bolivia to explore genetic basis of the Aymara high hemoglobin. Differentially expressed and spliced genes in Aymaras were associated with inflammatory and hypoxia-related pathways. We identified transcripts with 4^th or 5^th exon skipping of NFKB1 (AS-NFKB1), key part of NF-kB complex, and their splicing quantitative trait loci; these were increased in Aymaras. AS-NFKB1 transcripts correlated with both transcripts and protein levels of inflammatory and HIF-regulated genes, including hemoglobin. While overexpression of the AS-NFKB1 variant led to increased expression of inflammatory and HIF-targeted genes; under inflammatory stress, NF-kB protein translocation to the nucleus was attenuated, resulting in reduced expression of these genes. Our study reveals AS-NFKB1 splicing events correlating with increased hemoglobin in Aymara and their possible protective mechanisms against excessive inflammation.