Epigenetic control of adaptive or homeostatic splicing during interval-training activities.

Interval-training activities induce adaptive cellular changes without altering their fundamental identity, but the precise underlying molecular mechanisms are not fully understood. In this study, we demonstrate that interval-training depolarization (ITD) of pituitary cells triggers distinct adaptive or homeostatic splicing responses of alternative exons. This occurs while preserving the steady-state expression of the Prolactin and other hormone genes. The nature of these splicing responses depends on the exon's DNA methylation status, the methyl-C-binding protein MeCP2 and its associated CA-rich motif-binding hnRNP L. Interestingly, the steady expression of the Prolactin gene is also reliant on MeCP2, whose disruption leads to exacerbated multi-exon aberrant splicing and overexpression of the hormone gene transcripts upon ITD, similar to the observed hyperprolactinemia or activity-dependent aberrant splicing in Rett Syndrome. Therefore, epigenetic control is crucial for both adaptive and homeostatic splicing and particularly the steady expression of the Prolactin hormone gene during ITD. Disruption in this regulation may have significant implications for the development of progressive diseases.

Variations of intronic branchpoint motif: identification and functional implications in splicing and disease.

The branchpoint (BP) motif is an essential intronic element for spliceosomal pre-mRNA splicing. In mammals, its sequence composition, distance to the downstream exon, and number of BPs per 3´ splice site are highly variable, unlike the GT/AG dinucleotides at the intron ends. These variations appear to provide evolutionary advantages for fostering alternative splicing, satisfying more diverse cellular contexts, and promoting resilience to genetic changes, thus contributing to an extra layer of complexity for gene regulation. Importantly, variants in the BP motif itself or in genes encoding BP-interacting factors cause human genetic diseases or cancers, highlighting the critical function of BP motif and the need to precisely identify functional BPs for faithful interpretation of their roles in splicing. In this perspective, we will succinctly summarize the major findings related to BP motif variations, discuss the relevant issues/challenges, and provide our insights.

Transcriptome-Wide Detection of Intron/Exon Definition in the Endogenous Pre-mRNA Transcripts of Mammalian Cells and Its Regulation by Depolarization.

Pairing of splice sites across an intron or exon is the central point of intron or exon definition in pre-mRNA splicing with the latter mode proposed for most mammalian exons. However, transcriptome-wide pairing within endogenous transcripts has not been examined for the prevalence of each mode in mammalian cells. Here we report such pairings in rat GH3 pituitary cells by measuring the relative abundance of nuclear RNA-Seq reads at the intron start or end (RISE). Interestingly, RISE indexes are positively correlated between 5' and 3' splice sites specifically across introns or exons but inversely correlated with the usage of adjacent exons. Moreover, the ratios between the paired indexes were globally modulated by depolarization, which was disruptible by 5-aza-Cytidine. The nucleotide matrices of the RISE-positive splice sites deviate significantly from the rat consensus, and short introns or exons are enriched with the cross-intron or -exon RISE pairs, respectively. Functionally, the RISE-positive genes cluster for basic cellular processes including RNA binding/splicing, or more specifically, hormone production if regulated by depolarization. Together, the RISE analysis identified the transcriptome-wide regulation of either intron or exon definition between weak splice sites of short introns/exons in mammalian cells. The analysis also provides a way to further track the splicing intermediates and intron/exon definition during the dynamic regulation of alternative splicing by extracellular factors.

Forming cytoplasmic stress granules PURα suppresses mRNA translation initiation of IGFBP3 to promote esophageal squamous cell carcinoma progression.

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal malignancies worldwide. Recently, our group identified purine-rich element binding protein alpha (PURα), a single-stranded DNA/RNA-binding protein, to be significantly associated with the progression of ESCC. Additional immunofluorescence staining demonstrated that PURα forms cytoplasmic stress granules to suppress mRNA translation initiation. The expression level of cytoplasmic PURα in ESCC tumor tissues was significantly higher than that in adjacent epithelia and correlated with a worse patient survival rate by immunohistochemistry. Functionally, PURα strongly preferred to bind to UG-/U-rich motifs and mRNA 3´UTR by CLIP-seq analysis. Moreover, PURα knockout significantly increased the protein level of insulin-like growth factor binding protein 3 (IGFBP3). In addition, it was further demonstrated that PURα-interacting proteins are remarkably associated with translation initiation factors and ribosome-related proteins and that PURα regulates protein expression by interacting with translation initiation factors, such as PABPC1, eIF3B and eIF3F, in an RNA-independent manner, while the interaction with ribosome-related proteins is significantly dependent on RNA. Specifically, PURα was shown to interact with the mRNA 3´UTR of IGFBP3 and inhibit its expression by suppressing mRNA translation initiation. Together, this study identifies cytoplasmic PURα as a modulator of IGFBP3, which could be a promising therapeutic target for ESCC treatment.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity.

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R1C2T3R4A5Y6, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G4/A4 ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G4 increase is consistent with A4 to G4 transitions in evolution. Of the G4/A4-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G3/A3 co-vary with the G4 to some extent. However, corresponding increase of the G4-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G4-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G4/A4 ratio correlates well with the abundance of alternative splicing in the two phyla, and G4 enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

Transcriptome protection by the expanded family of hnRNPs.

The family of heterogeneous ribonucleoproteins (hnRNPs) have multiple functions in RNA metabolism. In recent years, several hnRNPs have also been shown to be essential for the maintenance of transcriptome integrity, by preventing intronic cryptic splicing signals from mis-splicing of many endogeneous pre-mRNA transcripts. Here we discuss the possibility for a general role of this family of proteins and their expansion in transcriptome protection.

Multilevel Differential Control of Hormone Gene Expression Programs by hnRNP L and LL in Pituitary Cells.

The pituitary-derived somatolactotrophe GH3 cells secrete both growth hormone (GH) and prolactin (PRL). We have found that the hnRNP L and L-like (LL) paralogs differentially regulate alternative splicing of genes in these cells. Here, we show that hnRNP L is essential for PRL only, but LL is essential for both PRL and GH production. Transcriptome-wide RNA sequencing (RNA-Seq) analysis indicates that they differentially control groups of hormone or hormone-related genes involved in hormone production/regulation at total transcript and alternative exon levels. Interestingly, hnRNP L also specifically binds and prevents the aberrant usage of a nonconserved CA-rich intron piece of Prl pre-mRNA transcripts, and many others involved in endocrine functions, to prevent mostly cryptic last exons and mRNA truncation. Essential for the full hnRNP L effect on specific exons is a proline-rich region that emerged during evolution in vertebrate hnRNP L only but not LL. Together, our data demonstrate that the hnRNP L and its paralog, LL, differentially control hormone gene expression programs at multiple levels, and hnRNP L in particular is critical for protecting the transcriptome from aberrant usage of intronic sequences. The multilevel differential control by hnRNPs likely tailors the transcriptome to help refine and safeguard the different gene expression programs for different hormones.

The matrices and constraints of GT/AG splice sites of more than 1000 species/lineages.

To provide a resource for the splice sites (SS) of different species, we calculated the matrices of nucleotide compositions of about 38 million splice sites from >1000 species/lineages. The matrices are enriched of aGGTAAGT (5'SS) or (Y)6N(C/t)AG(g/a)t (3'SS) overall; however, they are quite diverse among hundreds of species. The diverse matrices remain prominent even under sequence selection pressures, suggesting the existence of diverse constraints as well as U snRNAs and other spliceosomal factors and/or their interactions with the splice sites. Using an algorithm to measure and compare the splice site constraints across all species, we demonstrate their distinct differences quantitatively. As an example of the resource's application to answering specific questions, we confirm that high constraints of particular positions are significantly associated with transcriptome-wide, increased occurrences of alternative splicing when uncommon nucleotides are present. More interestingly, the abundance of alternative splicing in 16 species correlates with the average constraint index of splice sites in a bell curve. This resource will allow users to assess specific sequences/splice sites against the consensus of every Ensembl-annotated species, and to explore the evolutionary changes or relationship to alternative splicing and transcriptome diversity. Web-search or update features are also included.

Widespread Separation of the Polypyrimidine Tract From 3' AG by G Tracts in Association With Alternative Exons in Metazoa and Plants.

At the end of introns, the polypyrimidine tract (Py) is often close to the 3' AG in a consensus (Y)20NCAGgt in humans. Interestingly, we have found that they could also be separated by purine-rich elements including G tracts in thousands of human genes. These regulatory elements between the Py and 3' AG (REPA) mainly regulate alternative 3' splice sites (3' SS) and intron retention. Here we show their widespread distribution and special properties across kingdoms. The purine-rich 3' SS are found in up to about 60% of the introns among more than 1,000 species/lineages by whole genome analysis, and up to 18% of these introns contain the REPA G-tracts (REPAG) in about 0.6 million of 3' SS in total. In particular, they are significantly enriched over their 3' SS and genome backgrounds in metazoa and plants, and highly associated with alternative splicing of genes in diverse functional clusters. Cryptic splice sites harboring such G- and the other purine-triplets tend to be enriched (2-9 folds over the disrupted canonical 3' SS) and aberrantly used in cancer patients carrying mutations of the SF3B1 or U2AF35, factors critical for branch point (BP) or 3' AG recognition, respectively. Moreover, the REPAGs are significantly associated with reduced occurrences of BP motifs between the -24 and -4 positions, in particular absent between the -7 and -5 positions in several model organisms examined. The more distant BPs are associated with increased occurrences of alternative splicing in humans and zebrafish. The REPAGs appear to have evolved in a species- or phylum-specific way. Thus, there is widespread separation of the Py and 3' AG by REPAGs that have evolved differentially. This special 3' SS arrangement likely contributes to the generation of diverse transcript or protein isoforms in biological functions or diseases through alternative or aberrant splicing.

Identification of Claudin 1 Transcript Variants in Human Invasive Breast Cancer.

BACKGROUND: The claudin 1 tight junction protein, solely responsible for the barrier function of epithelial cells, is frequently down regulated in invasive human breast cancer. The underlying mechanism is largely unknown, and no obvious mutations in the claudin 1 gene (CLDN1) have been identified to date in breast cancer. Since many genes have been shown to undergo deregulation through splicing and mis-splicing events in cancer, the current study was undertaken to investigate the occurrence of transcript variants for CLDN1 in human invasive breast cancer. METHODS: RT-PCR analysis of CLDN1 transcripts was conducted on RNA isolated from 12 human invasive breast tumors. The PCR products from each tumor were resolved by agarose gel electrophoresis, cloned and sequenced. Genomic DNA was also isolated from each of the 12 tumors and amplified using PCR CLDN1 specific primers. Sanger sequencing and single nucleotide polymorphism (SNP) analyses were conducted. RESULTS: A number of CLDN1 transcript variants were identified in these breast tumors. All variants were shorter than the classical CLDN1 transcript. Sequence analysis of the PCR products revealed several splice variants, primarily in exon 1 of CLDN1; resulting in truncated proteins. One variant, V1, resulted in a premature stop codon and thus likely led to nonsense mediated decay. Interestingly, another transcript variant, V2, was not detected in normal breast tissue samples. Further, sequence analysis of the tumor genomic DNA revealed SNPs in 3 of the 4 coding exons, including a rare missense SNP (rs140846629) in exon 2 which represents an Ala124Thr substitution. To our knowledge this is the first report of CLDN1 transcript variants in human invasive breast cancer. These studies suggest that alternate splicing may also be a mechanism by which claudin 1 is down regulated at both the mRNA and protein levels in invasive breast cancer and may provide novel insights into how CLDN1 is reduced or silenced in human breast cancer.

Multiple effects of curcumin on promoting expression of the exon 7-containing SMN2 transcript.

Survival of motor neuron 2 (SMN2) is a modifier gene for spinal muscular atrophy (SMA), a neurodegenerative disease caused by insufficient SMN protein mostly due to SMN1 defect. SMN2 is nearly identical to SMN1 but unfortunately only able to produce a small amount of SMN protein due to exon 7 skipping. The exon 7-containing SMN2 transcript (SMN2_E7+) can be increased by a dietary compound, curcumin, but the involved molecular changes are not clear. Here we have found that in fibroblast cells of a SMA type II patient, curcumin enhanced the inclusion of SMN2 exon 7. Examination of the potential splicing factors showed that curcumin specifically increased the protein and transcript levels of SRSF1. The increased SRSF1 protein was mainly nuclear and hyperphosphorylated. Interestingly, the curcumin effects on the SMN2 and SRSF1 transcripts were inhibited by a protein deacetylase inhibitor, trichostatin A. Moreover, in support of its role in the SMN2 splicing, knocking down SRSF1 reduced the inclusion of SMN2 exon 7. Thus, curcumin appears to have multiple effects on the SMN2 transcript and its splicing regulators, including the change of alternative splicing and transcript/protein level as well as phosphorylation. Protein deacetylases and phosphatases are likely involved in these effects. Interestingly, the effects all seem to favor production of the SMN2_E7+ transcript in SMA patient cells.

Diverse regulation of 3' splice site usage.

The regulation of splice site (SS) usage is important for alternative pre-mRNA splicing and thus proper expression of protein isoforms in cells; its disruption causes diseases. In recent years, an increasing number of novel regulatory elements have been found within or nearby the 3'SS in mammalian genes. The diverse elements recruit a repertoire of trans-acting factors or form secondary structures to regulate 3'SS usage, mostly at the early steps of spliceosome assembly. Their mechanisms of action mainly include: (1) competition between the factors for RNA elements, (2) steric hindrance between the factors, (3) direct interaction between the factors, (4) competition between two splice sites, or (5) local RNA secondary structures or longer range loops, according to the mode of protein/RNA interactions. Beyond the 3'SS, chromatin remodeling/transcription, posttranslational modifications of trans-acting factors and upstream signaling provide further layers of regulation. Evolutionarily, some of the 3'SS elements seem to have emerged in mammalian ancestors. Moreover, other possibilities of regulation such as that by non-coding RNA remain to be explored. It is thus likely that there are more diverse elements/factors and mechanisms that influence the choice of an intron end. The diverse regulation likely contributes to a more complex but refined transcriptome and proteome in mammals.

Differential expression, distinct localization and opposite effect on Golgi structure and cell differentiation by a novel splice variant of human PRMT5.

Alternative splicing contributes greatly to the proteomic diversity of metazoans. Protein arginine methyltransferase 5 (PRMT5) methylates arginines of Golgi components and other factors exerting diverse effects on cell growth/differentiation, but the underlying molecular basis for its subcellular distribution and diverse roles has not been fully understood. Here we show the detailed properties of an evolutionarily emerged splice variant of human PRMT5 (PRMT5S) that is distinct from the original isoform (PRMT5L). The isoforms are differentially expressed among mammalian cells and tissues. The PRMT5S is distributed all over the cell but PRMT5L mainly colocalizes with Giantin, a Golgi marker. PRMT5 knockdown led to an enlarged Giantin pattern, which was prevented by the expression of either isoform. Rescuing PRMT5S also increased the percentage of cells with an interphase Giantin pattern compacted at one end of the nucleus, consistent with its cell cycle-arresting effect, while rescuing PRMT5L increased that of the mitotic Giantin patterns of dynamically fragmented structures. Moreover, the isoforms are differentially expressed during neuronal or dendritic cell differentiation, and their ectopic expression showed an opposite effect on dendritic cell differentiation. Furthermore, besides their differential regulation of gene expression, both isoforms also similarly regulate over a thousand genes particularly those involved in apoptosis and differentiation. Taking these properties together, we propose that their differential expression and subcellular localization contribute to spatial and temporal regulation of arginine methylation and gene expression to exert different effects. The novel PRMT5S likely contributes to the observed diverse effects of PRMT5 in cells.

An endogenous 'non-specific' protein detected by a His-tag antibody is human transcription regulator YY1.

Histidine-tags have been used for a wide variety of experiments including protein purification, Western blots, immunoprecipitation and immunohistochemistry. In our previous studies, we have repeatedly detected a 'non-specific' endogenous protein of about 60 kD in Western blots of protein lysates from HEK293T or HeLa cells using the anti-His-tag antibody (His-probe (H3), catalogue #, SC-8036, Santa Cruz Biotech. Co.) (Yu et al., J. Biol. Chem. 284 (2009) 1505-1513). Here we have immunoprecipitated the protein from HeLa nuclear extracts using the anti-His-tag antibody, excised the 60 kD band and subjected it to LC-MS/MS (Fig. 1). The deduced sequences of two peptides of the protein match the human transcriptional regulator YY1 (Yin and Yang 1, UniProt ID, P25490, Fig. 2), which contains 11 histidine residues in a stretch (from amino acid 70 to 80) at its NH2-terminal region without known functions (Lee et al., Nucleic Acids Res. 23 (1995) 925-931; Bushmeyer et al., J. Biol. Chem. 270 (1995) 30213-30220). Since genes encoding other Histidine-repeat proteins also exist in the genome (Salichs et al., PLoS Genet. 5 (2009) e1000397), it is possible that YY1 might not be the only endogenous protein that could be expressed and recognized by the antibody in different sources of samples in future experiments. The presence of various endogenous histidine-repeat proteins suggests that data from experiments particularly immunostaining using His-tag antibodies need to be interpreted with caution. This might also be useful to the broader scientific community by providing an example for the interpretation of 'non-specific' bands in Western blots.

Increased stability of heterogeneous ribonucleoproteins by a deacetylase inhibitor.

Splicing factors are often influenced by various signaling pathways, contributing to the dynamic changes of protein isoforms in cells. Heterogeneous ribonucleoproteins (hnRNPs) regulate many steps of RNA metabolism including pre-mRNA splicing but their control by cell signaling particularly through acetylation and ubiquitination pathways remains largely unknown. Here we show that TSA, a deacetylase inhibitor, reduced the ratio of Bcl-x splice variants Bcl-xL/xS in MDA-MB-231 breast cancer cells. This TSA effect was independent of TGFß1; however, only in the presence of TGFß1 was TSA able to change the splicing regulators hnRNP F/H by slightly reducing their mRNA transcripts but strongly preventing protein degradation. The latter was also efficiently prevented by lactacystin, a proteasome inhibitor, suggesting their protein stability control by both acetylation and ubiquitination pathways. Three lysines K87, K98 and K224 of hnRNP F are potential targets of the mutually exclusive acetylation or ubiquitination (K(Ac/Ub)) in the protein modification database PhosphoSitePlus. Mutating each of them but not a control non-K(Ac/Ub) (K68) specifically abolished the TSA enhancement of protein stability. Moreover, mutating K98 (K98R) and K224 (K224R) also abolished the TSA regulation of alternative splicing of a Bcl-x mini-gene. Furthermore, about 86% (30 of 35) of the multi-functional hnRNP proteins in the database contain lysines that are potential sites for acetylation/ubiquitination. We demonstrate that the degradation of three of them (A1, I and L) are also prevented by TSA. Thus, the deacetylase inhibitor TSA enhances hnRNP F stability through the K(Ac/Ub) lysines, with some of them essential for its regulation of alternative splicing. Such a regulation of protein stability is perhaps common for a group of hnRNPs and RNA metabolism.

Increase of a group of PTC(+) transcripts by curcumin through inhibition of the NMD pathway.

Nonsense-mediated mRNA decay (NMD), an mRNA surveillance mechanism, eliminates premature termination codon-containing (PTC⁺) transcripts. For instance, it maintains the homeostasis of splicing factors and degrades aberrant transcripts of human genetic disease genes. Here we examine the inhibitory effect on the NMD pathway and consequent increase of PTC+ transcripts by the dietary compound curcumin. We have found that several PTC⁺ transcripts including that of serine/arginine-rich splicing factor 1 (SRSF1) were specifically increased in cells by curcumin. We also observed a similar curcumin effect on the PTC⁺ mutant transcript from a Tay-Sachs-causing HEXA allele or from a beta-globin reporter gene. The curcumin effect was accompanied by significantly reduced expression of the NMD factors UPF1, 2, 3A and 3B. Consistently, in chromatin immunoprecipitation assays, curcumin specifically reduced the occupancy of acetyl-histone H3 and RNA polymerase II at the promoter region (-376 to -247nt) of human UPF1, in a time- and dosage-dependent way. Importantly, knocking down UPF1 abolished or substantially reduced the difference of PTC(+) transcript levels between control and curcumin-treated cells. The disrupted curcumin effect was efficiently rescued by expression of exogenous Myc-UPF1 in the knockdown cells. Together, our data demonstrate that a group of PTC⁺ transcripts are stabilized by a dietary compound curcumin through the inhibition of UPF factor expression and the NMD pathway.

Evolutionary emergence of a novel splice variant with an opposite effect on the cell cycle.

Alternative splicing contributes greatly to the diversification of mammalian proteomes, but the molecular basis for the evolutionary emergence of splice variants remains poorly understood. We have recently found a novel class of splicing regulatory elements between the polypyrimidine tract (Py) and 3' AG (REPA) at intron ends in many human genes, including the multifunctional PRMT5 (for protein arginine methyltransferase 5) gene. The PRMT5 element is comprised of two G tracts that arise in most mammals and accompany significant exon skipping in human transcripts. The G tracts inhibit splicing by recruiting heterogeneous nuclear ribonucleoprotein (hnRNP) H and F (H/F) to reduce U2AF65 binding to the Py, causing exon skipping. The resulting novel shorter variant PRMT5S exhibits a histone H4R3 methylation effect similar to that seen with the original longer PRMT5L isoform but exhibits a distinct localization and preferential control of critical genes for cell cycle arrest at interphase in comparison to PRMT5L. This report thus provides a molecular mechanism for the evolutionary emergence of a novel splice variant with an opposite function in a fundamental cell process. The presence of REPA elements in a large group of genes implies their wider impact on different cellular processes for increased protein diversity in humans.

Evolutionarily emerged G tracts between the polypyrimidine tract and 3' AG are splicing silencers enriched in genes involved in cancer.

BACKGROUND: The 3' splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position "N" can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes. RESULTS: Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3' splice sites of annotated introns. This identified several thousand such 3'SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3' splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3'SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer. CONCLUSION: We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3'AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.