Population-dependent Intron Retention and DNA Methylation in Breast Cancer

Landmark Series: The Cancer Genome Atlas and the Study of Breast Cancer Disparities

Race-related variation in breast cancer incidence and mortality are well-documented in the United States. The effect of genetic ancestry on disparities in tumor genomics, risk factors, treatment, and outcomes of breast cancer is less understood. The Cancer Genome Atlas (TCGA) is a publicly available resource that has allowed for the recent emergence of genome analysis research seeking to characterize tumor DNA and protein expression by ancestry as well as the social construction of race and ethnicity. Results from TCGA based studies support previous clinical evidence that demonstrates that American women with African ancestry are more likely to be afflicted with breast cancers featuring aggressive biology and poorer outcomes compared with women with other backgrounds. Data from TCGA based studies suggest that Asian women have tumors with favorable immune microenvironments and may experience better disease-free survival compared with white Americans. TCGA contains limited data on Hispanic/Latinx patients due to small sample size. Overall, TCGA provides important opportunities to define the molecular, biologic, and germline genetic factors that contribute to breast cancer disparities.

Aberrant DNA Methylation, Expression, and Occurrence of Transcript Variants of the ABC Transporter ABCA7 in Breast Cancer

The ABC transporter ABCA7 has been found to be aberrantly expressed in a variety of cancer types, including breast cancer. We searched for specific epigenetic and genetic alterations and alternative splicing variants of ABCA7 in breast cancer and investigated whether these alterations are associated with ABCA7 expression. By analyzing tumor tissues from breast cancer patients, we found CpGs at the exon 5-intron 5 boundary aberrantly methylated in a molecular subtype-specific manner. The detection of altered DNA methylation in tumor-adjacent tissues suggests epigenetic field cancerization. In breast cancer cell lines, DNA methylation levels of CpGs in promoter-exon 1, intron 1, and at the exon 5-intron 5 boundary were not correlated with ABCA7 mRNA levels. By qPCR involving intron-specific and intron-flanking primers, we identified intron-containing ABCA7 mRNA transcripts. The occurrence of intron-containing transcripts was neither molecular subtype-specific nor directly correlated with DNA methylation at the respective exon-intron boundaries. Treatment of breast cancer cell lines MCF-7, BT-474, SK-BR3, and MDA-MB-231 with doxorubicin or paclitaxel for 72 h resulted in altered ABCA7 intron levels. Shotgun proteomics revealed that an increase in intron-containing transcripts was associated with significant dysregulation of splicing factors linked to alternative splicing.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Towards resolution of the intron retention paradox in breast cancer

BACKGROUND

After many years of neglect in the field of alternative splicing, the importance of intron retention (IR) in cancer has come into focus following landmark discoveries of aberrant IR patterns in cancer. Many solid and liquid tumours are associated with drastic increases in IR, and such patterns have been pursued as both biomarkers and therapeutic targets. Paradoxically, breast cancer (BrCa) is the only tumour type in which IR is reduced compared to adjacent normal breast tissue.

METHODS

In this study, we have conducted a pan-cancer analysis of IR with emphasis on BrCa and its subtypes. We explored mechanisms that could cause aberrant and pathological IR and clarified why normal breast tissue has unusually high IR.

RESULTS

Strikingly, we found that aberrantly decreasing IR in BrCa can be largely attributed to normal breast tissue having the highest occurrence of IR events compared to other healthy tissues. Our analyses suggest that low numbers of IR events in breast tumours are associated with poor prognosis, particularly in the luminal B subtype. Interestingly, we found that IR frequencies negatively correlate with cell proliferation in BrCa cells, i.e. rapidly dividing tumour cells have the lowest number of IR events. Aberrant RNA-binding protein expression and changes in tissue composition are among the causes of aberrantly decreasing IR in BrCa.

CONCLUSIONS

Our results suggest that IR should be considered for therapeutic manipulation in BrCa patients with aberrantly low IR levels and that further work is needed to understand the cause and impact of high IR in other tumour types.

Epigenetic interplay between methylation and miRNA in bladder cancer: focus on isoform expression

BACKGROUND

Various epigenetic factors are responsible for the non-genetic regulation on gene expression. The epigenetically dysregulated oncogenes or tumor suppressors by miRNA and/or DNA methylation are often observed in cancer cells. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear. In this study, we propose an integrative framework to systematically investigate epigenetic interactions between miRNA and methylation at the alternatively spliced mRNA level in bladder cancer. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear.

RESULTS

The integrative analyses yielded 136 significant combinations (methylation, miRNA and isoform). Further, overall survival analysis on the 136 combinations based on methylation and miRNA, high and low expression groups resulted in 13 combinations associated with survival. Additionally, different interaction patterns were examined.

CONCLUSIONS

Our study provides a higher resolution of molecular insight into the crosstalk between two epigenetic factors, DNA methylation and miRNA. Given the importance of epigenetic interactions and alternative splicing in cancer, it is timely to identify and understand the underlying mechanisms based on epigenetic markers and their interactions in cancer, leading to alternative splicing with primary functional impact.

Echinacea purpurea (L.) Moench treatment of monocytes promotes tonic interferon signaling, increased innate immunity gene expression and DNA repeat hypermethylated silencing of endogenous retroviral sequences

Background

Herbal remedies of Echinacea purpurea tinctures are widely used today to reduce common cold respiratory tract infections.

Methods

Transcriptome, epigenome and kinome profiling allowed a systems biology level characterisation of genomewide immunomodulatory effects of a standardized Echinacea purpurea (L.) Moench extract in THP1 monocytes.

Results

Gene expression and DNA methylation analysis revealed that Echinaforce® treatment triggers antiviral innate immunity pathways, involving tonic IFN signaling, activation of pattern recognition receptors, chemotaxis and immunometabolism. Furthermore, phosphopeptide based kinome activity profiling and pharmacological inhibitor experiments with filgotinib confirm a key role for Janus Kinase (JAK)-1 dependent gene expression changes in innate immune signaling. Finally, Echinaforce® treatment induces DNA hypermethylation at intergenic CpG, long/short interspersed nuclear DNA repeat elements (LINE, SINE) or long termininal DNA repeats (LTR). This changes transcription of flanking endogenous retroviral sequences (HERVs), involved in an evolutionary conserved (epi) genomic protective response against viral infections.

Conclusions

Altogether, our results suggest that Echinaforce® phytochemicals strengthen antiviral innate immunity through tonic IFN regulation of pattern recognition and chemokine gene expression and DNA repeat hypermethylated silencing of HERVs in monocytes. These results suggest that immunomodulation by Echinaforce® treatment holds promise to reduce symptoms and duration of infection episodes of common cold corona viruses (CoV), Severe Acute Respiratory Syndrome (SARS)-CoV, and new occurring strains such as SARS-CoV-2, with strongly impaired interferon (IFN) response and weak innate antiviral defense.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12906-021-03310-5.

Polymorphisms in eicosanoid-related biosynthesis enzymes associated with acute urticaria/angioedema induced by nonsteroidal anti-inflammatory drug hypersensitivity

BACKGROUND

Nonsteroidal anti-inflammatory drugs (NSAIDs) are the main triggers of drug hypersensitivity, with NSAID-induced acute urticaria/angioedema (NIUA) the most frequent phenotype. NSAID hypersensitivity is caused by cyclooxygenase 1 inhibition, which leads to an imbalance in prostaglandin (PG) and cysteinyl leukotriene (CysLT) synthesis. As only susceptible individuals develop NSAID hypersensitivity, genetic factors are believed to be involved; however, no study has assessed the overall genetic variability of key enzymes in PG and CysLT synthesis in NSAID hypersensitivity.

OBJECTIVES

To evaluate simultaneously variants in the main genes involved in PG and CysLT biosynthesis in NIUA.

METHODS

Two independent cohorts of patients were recruited in Spain, alongside NSAID-tolerant controls. The discovery cohort included only patients with NIUA; the replication cohort included patients with NSAID-exacerbated respiratory disease (NERD). A set of tagging single-nucleotide polymorphisms (tagSNPs) in PTGS1, PTGS2, ALOX5 and LTC4S was genotyped using mass spectrometry coupled with endpoint polymerase chain reaction.

RESULTS

The study included 1272 individuals. Thirty-five tagSNPs were successfully genotyped in the discovery cohort, with three being significantly associated after Bonferroni correction (rs10306194 and rs1330344 in PTGS1; rs28395868 in ALOX5). These polymorphisms were genotyped in the replication cohort: rs10306194 and rs28395868 remained associated with NIUA, and rs28395868 was marginally associated with NERD. Odds ratios (ORs) in the combined analysis (discovery and replication NIUA populations) were 1·7 for rs10306194 [95% confidence interval (CI) 1·34-2·14; P_corrected = 2·83 × 10^-4 ) and 2·19 for rs28395868 (95% CI 1·43-3·36; P_corrected = 0·002).

CONCLUSIONS

Variants of PTGS1 and ALOX5 may play a role in NIUA and NERD, supporting the proposed mechanisms of NSAID-hypersensitivity and shedding light on their genetic basis.

Comprehensive analysis based on DNA methylation and RNA-seq reveals hypermethylation of the up-regulated WT1 gene with potential mechanisms in PAM50 subtypes of breast cancer

Background

Breast cancer (BC), one of the most widespread cancers worldwide, caused the deaths of more than 600,000 women in 2018, accounting for about 15% of all cancer-associated deaths in women that year. In this study, we aimed to discover potential prognostic biomarkers and explore their molecular mechanisms in different BC subtypes using DNA methylation and RNA-seq.

Methods

We downloaded the DNA methylation datasets and the RNA expression profiles of primary tissues of the four BC molecular subtypes (luminal A, luminal B, basal-like, and HER2-enriched), as well as the survival information from The Cancer Genome Atlas (TCGA). The highly expressed and hypermethylated genes across all the four subtypes were screened. We examined the methylation sites and the downstream co-expressed genes of the selected genes and validated their prognostic value using a different dataset (GSE20685). For selected transcription factors, the downstream genes were predicted based on the Gene Transcription Regulation Database (GTRD). The tumor microenvironment was also evaluated based on the TCGA dataset.

Results

We found that Wilms tumor gene 1 (WT1), a transcription factor, was highly expressed and hypermethylated in all the four BC subtypes. All the WT1 methylation sites exhibited hypermethylation. The methylation levels of the TSS200 and 1stExon regions were negatively correlated with WT1 expression in two BC subtypes, while that of the gene body region was positively associated with WT1 expression in three BC subtypes. Patients with low WT1 expression had better overall survival (OS). Five genes including COL11A1, GFAP, FGF5, CD300LG, and IGFL2 were predicted as the downstream genes of WT1. Those five genes were dysregulated in the four BC subtypes. Patients with a favorable 6-gene signature (low expression of WT1 and its five predicted downstream genes) exhibited better OS than that with an unfavorable 6-gene signature. We also found a correlation between WT1 and tamoxifen using STITCH. Higher infiltration rates of CD8 T cells, plasma cells, and monocytes were found in the lower quartile WT1 group and the favorable 6-gene signature group. In conclusion, we demonstrated that WT1 is hypermethylated and up-regulated in the four BC molecular subtypes and a 6-gene signature may predict BC prognosis.

Holding on to Junk Bonds: Intron Retention in Cancer and Therapy

Intron retention (IR) in cancer was for a long time overlooked by the scientific community, as it was previously considered to be an artifact of a dysfunctional spliceosome. Technological advancements made in the last decade offer unique opportunities to explore the role of IR as a widespread phenomenon that contributes to the transcriptional diversity of many cancers. Numerous studies in cancer have shed light on dysregulation of cellular mechanisms that lead to aberrant and pathologic IR. IR is not merely a mechanism of gene regulation, but rather it can mediate cancer pathogenesis and therapeutic resistance in various human diseases. The burden of IR in cancer is governed by perturbations to mechanisms known to regulate this phenomenon and include epigenetic variation, mutations within the gene body, and splicing factor dysregulation. This review summarizes possible causes for aberrant IR and discusses the role of IR in therapy or as a consequence of disease treatment. As neoepitopes originating from retained introns can be presented on the cancer cell surface, the development of personalized cancer vaccines based on IR-derived neoepitopes should be considered. Ultimately, a deeper comprehension about the origins and consequences of aberrant IR may aid in the development of such personalized cancer vaccines.

Alternative methylation of intron motifs is associated with cancer-related gene expression in both canine mammary tumor and human breast cancer

Background

Canine mammary tumor (CMT) has long been considered as a good animal model for human breast cancer (HBC) due to their pathological and biological similarities. However, only a few aspects of the epigenome have been explored in both HBC and CMT. Moreover, DNA methylation studies have mainly been limited to the promoter regions of genes.

Results

Genome-wide methylation analysis was performed in CMT and adjacent normal tissues and focused on the intron regions as potential targets for epigenetic regulation. As expected, many tumor suppressors and oncogenes were identified. Of note, most cancer-associated biological processes were enriched in differentially methylated genes (DMGs) that included intron DMRs (differentially methylated regions). Interestingly, two PAX motifs, PAX5 (tumor suppressive) and PAX6 (oncogenic), were frequently found in hyper- and hypomethylated intron DMRs, respectively. Hypermethylation at the PAX5 motifs in the intron regions of CDH5 and LRIG1 genes were found to be anti-correlated with gene expression, while CDH2 and ADAM19 genes harboring hypomethylated PAX6 motifs in their intron region were upregulated. These results were validated from the specimens originally MBD-sequenced as well as additional clinical samples. We also comparatively investigated the intron methylation and downstream gene expression of these genes using human breast invasive carcinoma (BRCA) datasets in TCGA (The Cancer Genome Atlas) public database. Regional alteration of methylation was conserved in the corresponding intron regions and, consequently, gene expression was also altered in HBC.

Conclusions

This study provides good evidence for the conservation of epigenetic regulation in CMT and HBC, and suggests that intronic methylation can be an important factor in better understanding gene regulation in both CMT and HBC.

Intron Retention as a Mode for RNA-Seq Data Analysis

Intron retention (IR) is an alternative splicing mode whereby introns, rather than being spliced out as usual, are retained in mature mRNAs. It was previously considered a consequence of mis-splicing and received very limited attention. Only recently has IR become of interest for transcriptomic data analysis owing to its recognized roles in gene expression regulation and associations with complex diseases. In this article, we first review the function of IR in regulating gene expression in a number of biological processes, such as neuron differentiation and activation of CD4⁺ T cells. Next, we briefly review its association with diseases, such as Alzheimer's disease and cancers. Then, we describe state-of-the-art methods for IR detection, including RNA-seq analysis tools IRFinder and iREAD, highlighting their underlying principles and discussing their advantages and limitations. Finally, we discuss the challenges for IR detection and potential ways in which IR detection methods could be improved.

Characterization of Global DNA Methylation in Different Gene Regions Reveals Candidate Biomarkers in Pigs with High and Low Levels of Boar Taint

DNA methylation of different gene components, including different exons and introns, or different lengths of exons and introns is associated with differences in gene expression. To investigate the methylation of porcine gene components associated with the boar taint (BT) trait, this study used reduced representation bisulfite sequencing (RRBS) data from nine porcine testis samples in three BT groups (low, medium and high BT). The results showed that the methylation levels of the first exons and first introns were lower than those of the other exons and introns. The first exons/introns of CpG island regions had even lower levels of methylation. A total of 123 differentially methylated promoters (DMPs), 194 differentially methylated exons (DMEs) and 402 differentially methylated introns (DMIs) were identified, of which 80 DMPs (DMP-CpGis), 112 DMEs (DME-CpGis) and 166 DMIs (DMI-CpGis) were discovered in CpG islands. Importantly, GPX1 contained one each of DMP, DME, DMI, DMP-CpGi, DME-CpGi and DMI-CpGi. Gene-GO term relationships and pathways analysis showed DMP-CpGi-related genes are mainly involved in methylation-related biological functions. In addition, gene-gene interaction networks consisted of nodes that were hypo-methylated GPX1, hypo-methylated APP, hypo-methylated ATOX1, hyper-methylated ADRB2, hyper-methylated RPS6KA1 and hyper-methylated PNMT. They could be used as candidate biomarkers for reducing boar taint in pigs, after further validation in large cohorts.

The changing paradigm of intron retention: regulation, ramifications and recipes

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

Dissecting splicing decisions and cell-to-cell variability with designed sequence libraries

Most human genes are alternatively spliced, allowing for a large expansion of the proteome. The multitude of regulatory inputs to splicing limits the potential to infer general principles from investigating native sequences. Here, we create a rationally designed library of >32,000 splicing events to dissect the complexity of splicing regulation through systematic sequence alterations. Measuring RNA and protein splice isoforms allows us to investigate both cause and effect of splicing decisions, quantify diverse regulatory inputs and accurately predict (R² = 0.73–0.85) isoform ratios from sequence and secondary structure. By profiling individual cells, we measure the cell-to-cell variability of splicing decisions and show that it can be encoded in the DNA and influenced by regulatory inputs, opening the door for a novel, single-cell perspective on splicing regulation.

Alternative splicing is regulated by multiple mechanisms. Here the authors employed designed splice site libraries and massively parallel reporter assays to dissect the regulatory complexity and cell-to-cell variability of splicing decisions and to build accurate predictive models.

CpG promoter hypo-methylation and up-regulation of microRNA-190b in hormone receptor-positive breast cancer

Estrogen receptor-positive breast cancer is subdivided into subtypes LuminalA and LuminalB, based on different expression patterns. MicroRNA-190b has been reported to be up-regulated in estrogen receptor-positive breast cancers. In this study we aimed to investigate the role of CpG promoter methylation in regulating miR-190b expression and its impact on clinical presentation and prognosis. DNA methylation analysis for the promotor of microRNA-190b was performed by pyrosequencing 549 primary breast tumors, of which 62 were carriers of the BRCA2 ^999del5 founder mutation, 71 proximal normal breast samples and 16 breast derived cell lines. MicroRNA-190b expression was analysed in 67 primary breast tumors, 14 paired normal breast samples and 16 breast derived cell lines. Tissue microarrays (TMAs) were available for ER (n = 436), PR (n = 436), HER-2 (N = 258) and Ki67 (n = 248). MiR-190b had reduced promoter methylation in estrogen receptor-positive breast cancers (P = 1.02e-12, Median values: ER+ 24.3, ER- 38.26) and miR-190b's expression was up-regulated in a correlative manner (P = 1.83e-06, Spearman's rho -0.62). Through breast cancer specific survival analysis, we demonstrated that LuminalA patients exhibiting miR-190b hypo-methylation had better survival than other patients (P = 0.034, HR = 0.29, 95% CI 0.09-0.91). We, furthermore, demonstrated that miR-190b hypo-methylation occurs less frequently in ER+ tumors from BRCA2 ^999del5 mutation carriers than in non-mutated individuals (P = 0.038, Χ ² = 4.32, n = 335). Our results suggest that upregulation of miR-190b may occur through loss of promoter DNA methylation during the development of estrogen-receptor (ER) positive breast cancers, and that miR-190b hypo-methylation leads to increased breast cancer specific survival within the LuminalA- subtype but not LuminalB.

In-utero epigenetic factors are associated with early-onset myopia in young children

Objectives

To assess whether epigenetic mechanisms affecting gene expression may be involved in the pathogenesis of early-onset myopia, we performed genome-wide DNA methylation analyses of umbilical cord tissues, and assessed any associations between CpG site-specific methylation and the development of the disorder when the children were 3 years old.

Methods

Genome-wide DNA methylation profiling of umbilical cord samples from 519 Singaporean infants involved in a prospective birth cohort ‘Growing Up in Singapore Towards healthy Outcomes’ (GUSTO) was performed using the Illumina Infinium HumanMethylation450K chip microarray. Multivariable logistic regression models were used to assess any associations between site-specific CpG methylation of umbilical cord tissue at birth and myopia risk in 3 year old children, adjusting for potential confounders. Gene expression of genes located near CpG sites that demonstrated statistically significant associations were measured in relevant ocular tissues using human and mouse fetal and adult eye samples.

Results

We identified statistically significant associations between DNA methylation levels at five CpG sites and early-onset myopia risk after correcting for multiple comparisons using a false discovery rate of 5%. Two statistically significant CpG sites were identified in intergenic regions: 8p23(p = 1.70×10⁻⁷) and 12q23.2(p = 2.53×10⁻⁷). The remaining 3 statistically significant CpG sites were identified within the following genes: FGB (4q28, p = 3.60×10⁻⁷), PQLC1 (18q23, p = 8.9×10⁻⁷) and KRT12 (17q21.2, p = 1.2×10⁻⁶). Both PQLC1 and KRT12 were found to be significantly expressed in fetal and adult cornea and sclera tissues in both human and mouse.

Conclusions

We identified five CpG methylation sites that demonstrate a statistically significant association with increased risk of developing early-onset myopia. These findings suggest that variability in the neonatal cord epigenome may influence early-onset myopia risk in children. Further studies of the epigenetic influences on myopia risk in larger study populations, and the associations with adulthood myopia risk are warranted.

Gene-Specific Intron Retention Serves as Molecular Signature that Distinguishes Melanoma from Non-Melanoma Cancer Cells in Greek Patients

Background: Skin cancer represents the most common human malignancy, and it includes BCC, SCC, and melanoma. Since melanoma is one of the most aggressive types of cancer, we have herein attempted to develop a gene-specific intron retention signature that can distinguish BCC and SCC from melanoma biopsy tumors. Methods: Intron retention events were examined through RT-sqPCR protocols, using total RNA preparations derived from BCC, SCC, and melanoma Greek biopsy specimens. Intron-hosted miRNA species and their target transcripts were predicted via the miRbase and miRDB bioinformatics platforms, respectively. Ιntronic ORFs were recognized through the ORF Finder application. Generation and visualization of protein interactomes were achieved by the IntAct and Cytoscape softwares, while tertiary protein structures were produced by using the I-TASSER online server. Results: c-MYC and Sestrin-1 genes proved to undergo intron retention specifically in melanoma. Interaction maps of proteins encoded by genes being potentially targeted by retained intron-accommodated miRNAs were generated and SRPX2 was additionally delivered to our melanoma-specific signature. Novel ORFs were identified in MCT4 and Sestrin-1 introns, with potentially critical roles in melanoma development. Conclusions: The property of c-MYC, Sestrin-1, and SRPX2 genes to retain specific introns could be clinically used to molecularly differentiate non-melanoma from melanoma tumors.

Identification of exon skipping events associated with Alzheimer's disease in the human hippocampus

BACKGROUND

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).

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.