SpliceAPP: an interactive web server to predict splicing errors arising from human mutations.

BACKGROUND: Splicing variants are a major class of pathogenic mutations, with their severity equivalent to nonsense mutations. However, redundant and degenerate splicing signals hinder functional assessments of sequence variations within introns, particularly at branch sites. We have established a massively parallel splicing assay to assess the impact on splicing of 11,191 disease-relevant variants. Based on the experimental results, we then applied regression-based methods to identify factors determining splicing decisions and their respective weights. RESULTS: Our statistical modeling is highly sensitive, accurately annotating the splicing defects of near-exon intronic variants, outperforming state-of-the-art predictive tools. We have incorporated the algorithm and branchpoint information into a web-based tool, SpliceAPP, to provide an interactive application. This user-friendly website allows users to upload any genetic variants with genome coordinates (e.g., chr15 74,687,208 A G), and the tool will output predictions for splicing error scores and evaluate the impact on nearby splice sites. Additionally, users can query branch site information within the region of interest. CONCLUSIONS: In summary, SpliceAPP represents a pioneering approach to screening pathogenic intronic variants, contributing to the development of precision medicine. It also facilitates the annotation of splicing motifs. SpliceAPP is freely accessible using the link https://bc.imb.sinica.edu.tw/SpliceAPP . Source code can be downloaded at https://github.com/hsinnan75/SpliceAPP .

De novo mutation and skewed X-inactivation in girl with BCAP31-related syndrome.

Full genome analysis of a young girl with deafness, dystonia, central hypomyelination, refractory seizure, and fluctuating liver function impairment revealed a heterozygous, de novo variant in the BCAP31 gene on chromosome Xq28 (NM_001256447.2:c.92G>A), mutations of which caused the X-linked recessive severe neurologic disorder deafness, dystonia, and cerebral hypomyelination. Reverse transcription-polymerase chain reaction of the patient's white blood cells showed the absence of wild-type BCAP31 messenger RNA (mRNA) but the presence of two novel BCAP31 mRNAs. The major alternatively spliced mRNA is due to Exon 2 skipping and the utilization of a new initiation site in Exon 3 that leads to a frameshift and truncated transcript while the minor novel mRNA has a 110 nucleotide insertion to Exon 2. Phasing studies showed that the de novo variant arose in the paternal X chromosome. X chromosome inactivation assay was done and confirmed that the patient's maternal X chromosome was preferentially inactivated, providing evidence that the mutated BCAP31 gene was the one predominantly expressed. According to the American College of Medical Genetics and Genomics guideline, this variant is deemed "pathogenic" (PS2, PS3, PM2, PP3, and PP4) and deleterious. This is the first reported female patient in BCAP31-related syndrome resulted from skewed X-inactivation and a de novo mutation in the active X chromosome.

Large-scale analysis of branchpoint usage across species and cell lines.

The coding sequence of each human pre-mRNA is interrupted, on average, by 11 introns that must be spliced out for proper gene expression. Each intron contains three obligate signals: a 5' splice site, a branch site, and a 3' splice site. Splice site usage has been mapped exhaustively across different species, cell types, and cellular states. In contrast, only a small fraction of branch sites have been identified even once. The few reported annotations of branch site are imprecise as reverse transcriptase skips several nucleotides while traversing a 2-5 linkage. Here, we report large-scale mapping of the branchpoints from deep sequencing data in three different species and in the SF3B1 K700E oncogenic mutant background. We have developed a novel method whereby raw lariat reads are refined by U2snRNP/pre-mRNA base-pairing models to return the largest current data set of branchpoint sequences with quality metrics. This analysis discovers novel modes of U2snRNA:pre-mRNA base-pairing conserved in yeast and provides insight into the biogenesis of intron circles. Finally, matching branch site usage with isoform selection across the extensive panel of ENCODE RNA-seq data sets offers insight into the mechanisms by which branchpoint usage drives alternative splicing.

RNA structure replaces the need for U2AF2 in splicing.

RNA secondary structure plays an integral role in catalytic, ribosomal, small nuclear, micro, and transfer RNAs. Discovering a prevalent role for secondary structure in pre-mRNAs has proven more elusive. By utilizing a variety of computational and biochemical approaches, we present evidence for a class of nuclear introns that relies upon secondary structure for correct splicing. These introns are defined by simple repeat expansions of complementary AC and GT dimers that co-occur at opposite boundaries of an intron to form a bridging structure that enforces correct splice site pairing. Remarkably, this class of introns does not require U2AF2, a core component of the spliceosome, for its processing. Phylogenetic analysis suggests that this mechanism was present in the ancestral vertebrate lineage prior to the divergence of tetrapods from teleosts. While largely lost from land dwelling vertebrates, this class of introns is found in 10% of all zebrafish genes.

The effects of structure on pre-mRNA processing and stability.

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

Integrative analysis of micro-RNA, gene expression, and survival of glioblastoma multiforme.

Glioblastoma multiforme (GBM), the most common type of malignant brain tumor, is highly fatal. Limited understanding of its rapid progression necessitates additional approaches that integrate what is known about the genomics of this cancer. Using a discovery set (n = 348) and a validation set (n = 174) of GBM patients, we performed genome-wide analyses that integrated mRNA and micro-RNA expression data from GBM as well as associated survival information, assessing coordinated variability in each as this reflects their known mechanistic functions. Cox proportional hazards models were used for the survival analyses, and nonparametric permutation tests were performed for the micro-RNAs to investigate the association between the number of associated genes and its prognostication. We also utilized mediation analyses for micro-RNA-gene pairs to identify their mediation effects. Genome-wide analyses revealed a novel pattern: micro-RNAs related to more gene expressions are more likely to be associated with GBM survival (P = 4.8 × 10(-5)). Genome-wide mediation analyses for the 32,660 micro-RNA-gene pairs with strong association (false discovery rate [FDR] < 0.01%) identified 51 validated pairs with significant mediation effect. Of the 51 pairs, miR-223 had 16 mediation genes. These 16 mediation genes of miR-223 were also highly associated with various other micro-RNAs and mediated their prognostic effects as well. We further constructed a gene signature using the 16 genes, which was highly associated with GBM survival in both the discovery and validation sets (P = 9.8 × 10(-6)). This comprehensive study discovered mediation effects of micro-RNA to gene expression and GBM survival and provided a new analytic framework for integrative genomics.

RNA structure in splicing: An evolutionary perspective.

Pre-mRNA splicing is a key post-transcriptional regulation process in which introns are excised and exons are ligated together. A novel class of structured intron was recently discovered in fish. Simple expansions of complementary AC and GT dimers at opposite boundaries of an intron were found to form a bridging structure, thereby enforcing correct splice site pairing across the intron. In some fish introns, the RNA structures are strong enough to bypass the need of regulatory protein factors for splicing. Here, we discuss the prevalence and potential functions of highly structured introns. In humans, structured introns usually arise through the co-occurrence of C and G-rich repeats at intron boundaries. We explore the potentially instructive example of the HLA receptor genes. In HLA pre-mRNA, structured introns flank the exons that encode the highly polymorphic ß sheet cleft, making the processing of the transcript robust to variants that disrupt splicing factor binding. While selective forces that have shaped HLA receptor are fairly atypical, numerous other highly polymorphic genes that encode receptors contain structured introns. Finally, we discuss how the elevated mutation rate associated with the simple repeats that often compose structured intron can make structured introns themselves rapidly evolving elements.

Differences in oral habit and lymphocyte subpopulation affect malignant transformation of patients with oral precancer.

BACKGROUND/PURPOSE: In Taiwan, the combination of betel quid chewing, alcohol consumption, and smoking habits increases oral cancer risk by 123-fold compared to persons without these habits. Lymphocyte populations in patients may potentially affect the malignant transformation of oral precancer. METHODS: A total of 28 patients with oral precancer from our previous cohort were enrolled in this study, and their personal information and oral habits were documented. Their lymphocyte populations (CD4+, CD8+, CD19+, and CD56+) and activation markers (CD25 and CD69) were determined by flow cytometry from 1999 to 2004. After follow up till December 2014, data of patients with/without malignant transformation were recorded, and the relation between oral habits and percentage of initial lymphocyte markers was evaluated using the Student t test and Fisher's exact test. RESULTS: Ten precancer patients developed oral squamous cell carcinoma with a mean period of malignant transformation of 6.8 ± 2.1 years. Patients with malignant transformation had a mean age of 48.4 ± 5.0 years (n = 10), relatively more than that of patients without malignant transformation (41.6 ± 6.3 years, n = 18) (p < 0.05). An increase was noted in the population of peripheral blood mononuclear cells expressing CD4+CD69+, CD19+CD69+, and CD56+CD69+ (p < 0.05) in precancer patients with malignant transformation. Alcohol consumption showed an association with the malignant transformation of patients with precancer (p = 0.030), whereas betel quid and smoking showed little effect. CONCLUSION: These results suggest that age, alcohol consumption, and early activation of T cells, B cells, and natural killer cells are crucial in the malignant transformation of oral precancer. Analysis of patient's lymphocyte populations may help predict the malignant transformation of oral precancer.

Transient CPEB dimerization and translational control.

During oocyte development, the cytoplasmic polyadenylation element-binding protein (CPEB) nucleates a set of factors on mRNA that controls cytoplasmic polyadenylation and translation. The regulation of polyadenylation is mediated in part through serial phosphorylations of CPEB, which control both the dynamic integrity of the cytoplasmic polyadenylation apparatus and CPEB stability, events necessary for meiotic progression. Because the precise stoichiometry between CPEB and CPE-containing RNA is responsible for the temporal order of mRNA polyadenylation during meiosis, we hypothesized that, if CPEB production exceeded the amount required to bind mRNA, the excess would be sequestered in an inactive form. One attractive possibility for the sequestration is protein dimerization. We demonstrate that not only does CPEB form a dimer, but dimerization requires its RNA-binding domains. Dimer formation prevents CPEB from being UV cross-linked to RNA, which establishes a second pool of CPEB that is inert for polyadenylation and translational control. During oocyte maturation, the dimers are degraded much more rapidly than the CPEB monomers, due to their greater affinity for polo-like kinase 1 (plx1) and the ubiquitin E3 ligase ß-TrCP. Because dimeric CPEB also binds cytoplasmic polyadenylation factors with greater affinity than monomeric CPEB, it may act as a hub or reservoir for the polyadenylation machinery. We propose that the balance between CPEB and its target mRNAs is maintained by CPEB dimerization, which inactivates spare proteins and prevents them from inducing polyadenylation of RNAs with low affinity binding sites. In addition, the dimers might serve as molecular hubs that release polyadenylation factors for translational activation upon CPEB dimer destruction.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of a viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectrometry. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of the first viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectroscopy. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

The nuclear experience of CPEB: implications for RNA processing and translational control.

CPEB is a sequence-specific RNA binding protein that promotes polyadenylation-induced translation in early development, during cell cycle progression and cellular senescence, and following neuronal synapse stimulation. It controls polyadenylation and translation through other interacting molecules, most notably the poly(A) polymerase Gld2, the deadenylating enzyme PARN, and the eIF4E-binding protein Maskin. Here, we report that CPEB shuttles between the nucleus and cytoplasm and that its export occurs via the CRM1-dependent pathway. In the nucleus of Xenopus oocytes, CPEB associates with lampbrush chromosomes and several proteins involved in nuclear RNA processing. CPEB also interacts with Maskin in the nucleus as well as with CPE-containing mRNAs. Although the CPE does not regulate mRNA export, it influences the degree to which mRNAs are translationally repressed in the cytoplasm. Moreover, CPEB directly or indirectly mediates the alternative splicing of at least one pre-mRNA in mouse embryo fibroblasts as well as certain mouse tissues. We propose that CPEB, together with Maskin, binds mRNA in the nucleus to ensure tight translational repression upon export to the cytoplasm. In addition, we propose that nuclear CPEB regulates specific pre-mRNA alternative splicing.

Mechanism and modeling of human disease-associated near-exon intronic variants that perturb RNA splicing.

It is estimated that 10%-30% of disease-associated genetic variants affect splicing. Splicing variants may generate deleteriously altered gene product and are potential therapeutic targets. However, systematic diagnosis or prediction of splicing variants is yet to be established, especially for the near-exon intronic splice region. The major challenge lies in the redundant and ill-defined branch sites and other splicing motifs therein. Here, we carried out unbiased massively parallel splicing assays on 5,307 disease-associated variants that overlapped with branch sites and collected 5,884 variants across the 5' splice region. We found that strong splice sites and exonic features preserve splicing from intronic sequence variation. Whereas the splice-altering mechanism of the 3' intronic variants is complex, that of the 5' is mainly splice-site destruction. Statistical learning combined with these molecular features allows precise prediction of altered splicing from an intronic variant. This statistical model provides the identity and ranking of biological features that determine splicing, which serves as transferable knowledge and out-performs the benchmarking predictive tool. Moreover, we demonstrated that intronic splicing variants may associate with disease risks in the human population. Our study elucidates the mechanism of splicing response of intronic variants, which classify disease-associated splicing variants for the promise of precision medicine.

Serum- and glucocorticoid-inducible kinase 1 enhances zif268 expression through the mediation of SRF and CREB1 associated with spatial memory formation.

Serum- and glucocorticoid-inducible kinase 1 (SGK1) has been shown to play an important role in spatial memory formation, but the molecular mechanism underlying this effect of SGK1 was not known. zif268 is an immediate early gene that is induced by water maze learning. To investigate the role of SGK1 in the regulation of zif268 expression, the dominant negative mutant of SGK1, SGK1 S422A, was infused to the hippocampal CA1 area of rats, and was found to decrease significantly the mRNA level of zif268 in both naïve animals and trained animals. SGK1 was also found to phosphorylate serum response factor (SRF) at Ser73, Ser75, and Ser99, and phosphorylate CREB1 at Ser133. Inhibition of SGK1 phosphorylation sites on SRF and CREB1 with alanine substitution significantly diminished SGK1-enhanced zif268 expression in the promoter-luciferase assay. SGK1 also phosphorylates Elk-1 and SGK1 phosphorylation of Elk-1 decreased the transcriptional activity of Elk-1. But SGK1 phosphorylation of Elk-1 did not affect SGK1-enhanced zif268 expression. Moreover, the phosphorylation of SGK1 was increased in rat CA1 area after water maze learning, accompanied by increased phosphorylation of SRF at Ser99 and increased phosphorylation of CREB1 at Ser133. All these effects were antagonized by SGK1 S422A transfection. These results together suggest that SGK1 enhances zif268 expression through the mediation of SRF and CREB1, and these signaling pathways are associated with spatial memory formation in rats.

Epigenome-wide profiling of DNA methylation in paired samples of adipose tissue and blood.

Many epigenetic association studies have attempted to identify DNA methylation markers in blood that are able to mirror those in target tissues. Although some have suggested potential utility of surrogate epigenetic markers in blood, few studies have collected data to directly compare DNA methylation across tissues from the same individuals. Here, epigenomic data were collected from adipose tissue and blood in 143 subjects using Illumina HumanMethylation450 BeadChip array. The top axis of epigenome-wide variation differentiates adipose tissue from blood, which is confirmed internally using cross-validation and externally with independent data from the two tissues. We identified 1,285 discordant genes and 1,961 concordant genes between blood and adipose tissue. RNA expression data of the two classes of genes show consistent patterns with those observed in DNA methylation. The discordant genes are enriched in biological functions related to immune response, leukocyte activation or differentiation, and blood coagulation. We distinguish the CpG-specific correlation from the within-subject correlation and emphasize that the magnitude of within-subject correlation does not guarantee the utility of surrogate epigenetic markers. The study reinforces the critical role of DNA methylation in regulating gene expression and cellular phenotypes across tissues, and highlights the caveats of using methylation markers in blood to mirror the corresponding profile in the target tissue.

An unusual two-step control of CPEB destruction by Pin1.

Cytoplasmic polyadenylation is a conserved mechanism that controls mRNA translation and stability. A key protein that promotes polyadenylation-induced translation of mRNAs in maturing Xenopus oocytes is the cytoplasmic polyadenylation element binding protein (CPEB). During this meiotic transition, CPEB is subjected to phosphorylation-dependent ubiquitination and partial destruction, which is necessary for successive waves of polyadenylation of distinct mRNAs. Here we identify the peptidyl-prolyl cis-trans isomerase Pin1 as an important factor mediating CPEB destruction. Pin1 interacts with CPEB in an unusual manner in which it occurs prior to CPEB phosphorylation and prior to Pin1 activation by serine 71 dephosphorylation. Upon induction of maturation, CPEB becomes phosphorylated, which occurs simultaneously with Pin1 dephosphorylation. At this time, the CPEB-Pin1 interaction requires cdk1-catalyzed CPEB phosphorylation on S/T-P motifs. Subsequent CPEB ubiquitination and destruction are mediated by a conformational change induced by Pin1 isomerization of CPEB. Similar to M phase progression in maturing Xenopus oocytes, the destruction of CPEB during the mammalian cell cycle requires Pin1 as well. These data identify Pin1 as a new and essential factor regulating CPEB degradation.

Cytoplasmic polyadenylation and cytoplasmic polyadenylation element-dependent mRNA regulation are involved in Xenopus retinal axon development.

BACKGROUND: Translation in axons is required for growth cone chemotropic responses to many guidance cues. Although locally synthesized proteins are beginning to be identified, how specific mRNAs are selected for translation remains unclear. Control of poly(A) tail length by cytoplasmic polyadenylation element (CPE) binding protein 1 (CPEB1) is a conserved mechanism for mRNA-specific translational regulation that could be involved in regulating translation in axons. RESULTS: We show that cytoplasmic polyadenylation is required in Xenopus retinal ganglion cell (RGC) growth cones for translation-dependent, but not translation-independent, chemotropic responses in vitro, and that inhibition of CPE binding through dominant-negative interference severely reduces axon outgrowth in vivo. CPEB1 mRNA transcripts are present at low levels in RGCs but, surprisingly, CPEB1 protein was not detected in eye or brain tissue, and CPEB1 loss-of-function does not affect chemotropic responses or pathfinding in vivo. UV cross-linking experiments suggest that CPE-binding proteins other than CPEB1 in the retina regulate retinal axon development. CONCLUSION: These results indicate that cytoplasmic polyadenylation and CPE-mediated translational regulation are involved in retinal axon development, but that CPEB1 may not be the key regulator of polyadenylation in the developing retina.

CPEB3 and CPEB4 in neurons: analysis of RNA-binding specificity and translational control of AMPA receptor GluR2 mRNA.

CPEB is a sequence-specific RNA-binding protein that promotes polyadenylation-induced translation in oocytes and neurons. Vertebrates contain three additional genes that encode CPEB-like proteins, all of which are expressed in the brain. Here, we use SELEX, RNA structure probing, and RNA footprinting to show that CPEB and the CPEB-like proteins interact with different RNA sequences and thus constitute different classes of RNA-binding proteins. In transfected neurons, CPEB3 represses the translation of a reporter RNA in tethered function assays; in response to NMDA receptor activation, translation is stimulated. In contrast to CPEB, CPEB3-mediated translation is unlikely to involve cytoplasmic polyadenylation, as it requires neither the cis-acting AAUAAA nor the trans-acting cleavage and polyadenylation specificity factor, both of which are necessary for CPEB-induced polyadenylation. One target of CPEB3-mediated translation is GluR2 mRNA; not only does CPEB3 bind this RNA in vitro and in vivo, but an RNAi knockdown of CPEB3 in neurons results in elevated levels of GluR2 protein. These results indicate that CPEB3 is a sequence-specific translational regulatory protein.