Sequence conservation, relative isoform frequencies, and nonsense-mediated decay in evolutionarily conserved alternative splicing

Tau-induced deficits in nonsense-mediated mRNA decay contribute to neurodegeneration

INTRODUCTION

While brains of patients with Alzheimer's disease and related tauopathies have evidence of altered RNA processing, we lack a mechanistic understanding of how altered RNA processing arises in these disorders and if such changes are causally linked to neurodegeneration.

METHODS

Using Drosophila melanogaster models of tauopathy, we find that overall activity of nonsense-mediated mRNA decay (NMD), a key RNA quality-control mechanism, is reduced. Genetic manipulation of NMD machinery significantly modifies tau-induced neurotoxicity, suggesting that deficits in NMD are causally linked to neurodegeneration. Mechanistically, we find that deficits in NMD are a consequence of aberrant RNA export and RNA accumulation within nuclear envelope invaginations in tauopathy. We identify a pharmacological activator of NMD that suppresses neurodegeneration in tau transgenic Drosophila, indicating that tau-induced deficits in RNA quality control are druggable.

DISCUSSION

Our studies suggest that NMD activators should be explored for their potential therapeutic value to patients with tauopathies.

Inefficient splicing of long non-coding RNAs is associated with higher transcript complexity in human and mouse

Recent reports show that long non-coding RNAs (lncRNAs) have inefficient splicing and fewer alternative splice variants than mRNAs. Here, we have explored the efficiency of lncRNAs and mRNAs in producing various splice variants, given the number of exons in humans and mice. Intriguingly, lncRNAs produce more splice variants per exon, referred to as Transcript Complexity, than mRNAs. Most lncRNA splice variants are the product of the alternative last exon and exon skipping. LncRNAs and mRNAs with higher transcript complexity have shorter intron lengths. Longer exon length and GC/AG at 5'/3' splice sites are associated with higher transcript complexity in lncRNAs. Lastly, our results indicate that inefficient splicing of lncRNAs may facilitate multiple introns splicing and, thus, more spliced products per exon.

Ubiquitous conservative interaction patterns between post-spliced introns and their mRNAs revealed by genome-wide interspecies comparison

Introns, as important vectors of biological functions, can influence many stages of mRNA metabolism. However, in recent research, post-spliced introns are rarely considered. In this study, the optimal matched regions between introns and their mRNAs in nine model organism genomes were investigated with improved Smith-Waterman local alignment software. Our results showed that the distributions of mRNA optimal matched frequencies were highly consistent or universal. There are optimal matched frequency peaks in the UTR regions, which are obvious, especially in the 3'-UTR. The matched frequencies are relatively low in the CDS regions of the mRNA. The distributions of the optimal matched frequencies around the functional sites are also remarkably changed. The centers of the GC content distributions for different sequences are different. The matched rate distributions are highly consistent and are located mainly between 60% and 80%. The most probable value of the optimal matched segments is about 20 bp for lower eukaryotes and 30 bp for higher eukaryotes. These results show that there are abundant functional units in the introns, and these functional units are correlated structurally with all kinds of sequences of mRNA. The interaction between the post-spliced introns and their corresponding mRNAs may play a key role in gene expression.

FOXA1 regulates alternative splicing in prostate cancer

Dysregulation of alternative splicing in prostate cancer is linked to transcriptional programs activated by AR, ERG, FOXA1, and MYC. Here, we show that FOXA1 functions as the primary orchestrator of alternative splicing dysregulation across 500 primary and metastatic prostate cancer transcriptomes. We demonstrate that FOXA1 binds to the regulatory regions of splicing-related genes, including HNRNPK and SRSF1. By controlling trans-acting factor expression, FOXA1 exploits an "exon definition" mechanism calibrating alternative splicing toward dominant isoform production. This regulation especially impacts splicing factors themselves and leads to a reduction of nonsense-mediated decay (NMD)-targeted isoforms. Inclusion of the NMD-determinant FLNA exon 30 by FOXA1-controlled oncogene SRSF1 promotes cell growth in vitro and predicts disease recurrence. Overall, we report a role for FOXA1 in rewiring the alternative splicing landscape in prostate cancer through a cascade of events from chromatin access, to splicing factor regulation, and, finally, to alternative splicing of exons influencing patient survival.

The endogenous HBZ interactome in ATL leukemic cells reveals an unprecedented complexity of host interacting partners involved in RNA splicing

Adult T-cell leukemia/lymphoma (ATL) is a T-cell lymphoproliferative neoplasm caused by the human T-cell leukemia virus type 1 (HTLV-1). Two viral proteins, Tax-1 and HBZ play important roles in HTLV-1 infectivity and in HTLV-1-associated pathologies by altering key pathways of cell homeostasis. However, the molecular mechanisms through which the two viral proteins, particularly HBZ, induce and/or sustain the oncogenic process are still largely elusive. Previous results suggested that HBZ interaction with nuclear factors may alter cell cycle and cell proliferation. To have a more complete picture of the HBZ interactions, we investigated in detail the endogenous HBZ interactome in leukemic cells by immunoprecipitating the HBZ-interacting complexes of ATL-2 leukemic cells, followed by tandem mass spectrometry analyses. RNA seq analysis was performed to decipher the differential gene expression and splicing modifications related to HTLV-1. Here we compared ATL-2 with MOLT-4, a non HTLV-1 derived leukemic T cell line and further compared with HBZ-induced modifications in an isogenic system composed by Jurkat T cells and stably HBZ transfected Jurkat derivatives. The endogenous HBZ interactome of ATL-2 cells identified 249 interactors covering three main clusters corresponding to protein families mainly involved in mRNA splicing, nonsense-mediated RNA decay (NMD) and JAK-STAT signaling pathway. Here we analyzed in detail the cluster involved in RNA splicing. RNAseq analysis showed that HBZ specifically altered the transcription of many genes, including crucial oncogenes, by affecting different splicing events. Consistently, the two RNA helicases, members of the RNA splicing family, DDX5 and its paralog DDX17, recently shown to be involved in alternative splicing of cellular genes after NF-κB activation by HTLV-1 Tax-1, interacted and partially co-localized with HBZ. For the first time, a complete picture of the endogenous HBZ interactome was elucidated. The wide interaction of HBZ with molecules involved in RNA splicing and the subsequent transcriptome alteration strongly suggests an unprecedented complex role of the viral oncogene in the establishment of the leukemic state.

Alternative splicing patterns of doublesex reveal a missing link between Nix and doublesex in the sex determination cascade of Aedes albopictus

Sexual development in insects is regulated by a complicated hierarchical cascade of sex determination. The primary signals are diverse, whereas the central nexus doublesex (dsx) gene is relatively conserved within the pathway. Aedes (Stegomyia) albopictus is an important vector with an extensive worldwide distribution. We previously reported that Ae. albopictus dsx (Aalbdsx) yields one male- (Aalbdsx^M ) and three female-specific isoforms (Aalbdsx^F1-3 ); however, the spatiotemporal expression profiles and mechanisms regulating sex-specific alternative splicing require further investigation. In this study, we demonstrated that the Aalbdsx^M messenger RNA (mRNA) represents the default pattern when analyzed in human foreskin fibroblasts and HeLa cells. We combined reverse transcription polymerase chain reaction with RNA immunoprecipitation using specific antibodies against tagged Ae. albopictus male-determining factor AalNix and confirmed that AalNix indirectly regulates dsx pre-mRNA and regulates its alternative splicing. During the early embryo stage (0-2 and 4-8 h), maternal dsx^F and default splicing dsx^M were detected in both sexes; the expression of dsx^M then decreased until sufficient AalNix transcripts accumulated in male embryos at 20-24 h. These findings suggest that one or more potential dsx splicing enhancers can shift dsx^M to dsx^F in both sexes; however, the presence of Nix influences the function of this unknown splicing enhancer and ultimately leads to the formation of dsx^M in males. Finally, our results provide important insight into the regulatory mechanism of dsx alternative splicing in the mosquito.

Evolutionary dynamics of circular RNAs in primates

Many primate genes produce circular RNAs (circRNAs). However, the extent of circRNA conservation between closely related species remains unclear. By comparing tissue-specific transcriptomes across over 70 million years of primate evolution, we identify that within 3 million years circRNA expression profiles diverged such that they are more related to species identity than organ type. However, our analysis also revealed a subset of circRNAs with conserved neural expression across tens of millions of years of evolution. By comparing to species-specific circRNAs, we identified that the downstream intron of the conserved circRNAs display a dramatic lengthening during evolution due to the insertion of novel retrotransposons. Our work provides comparative analyses of the mechanisms promoting circRNAs to generate increased transcriptomic complexity in primates.

Splice variants of metabolic nuclear receptors: Relevance for metabolic disease and therapeutic targeting

Metabolic nuclear receptors are ligand-activated transcription factors which control a wide range of metabolic processes and signaling pathways in response to nutrients and xenobiotics. Targeting these NRs is at the forefront of our endeavours to generate novel treatment options for diabetes, metabolic syndrome and fatty liver disease. Numerous splice variants have been described for these metabolic receptors. Structural changes, as a result of alternative splicing, lead to functional differences among NR isoforms, resulting in the regulation of different metabolic pathways by these NR splice variants. In this review, we describe known splice variants of FXR, LXRs, PXR, RXR, LRH-1, CAR and PPARs. We discuss their structure and functions, and elaborate on the regulation of splice variant abundance by nutritional signals. We conclude that NR splice variants pose an intriguing new layer of complexity in metabolic signaling, which needs to be taken into account in the development of treatment strategies for metabolic diseases.

Expression and clinical importance of a newly discovered alternative splice variant of the gene for acrosin binding protein found in human brain tumors

Background

Acrosin binding protein (ACRBP) is a member of the cancer-testis antigen (CTA) family. Normally, ACRBP mRNA is expressed only in seminiferous tubules, while abnormally it is expressed in various types of cancers in tumor tissues, such as brain tumor.

Objectives

To determine the expression and clinical impact of a newly discovered splice variant of ACRBP in brain tumor.

Methods

Total RNA was extracted and reverse transcribed from 92 tumor specimens and 3 cell lines. Primers were designed to determine the expression of the new splice variant in all the samples. Quantitative real-time PCR (qPCR) was conducted for samples positive in reverse transcriptase-PCR. Association of the expression of ACRBP with the clinicopathological features of the various brain tumors was assessed statistically.

Results

The primers identified a newly discovered splice variant of ACRBP named ACRBP-V5a. The proportions of samples of the various brain tumor types positive for the ACRBP-V5a splicing variant were as follows: astrocytoma 10/33 (30%), glioblastoma 10/30 (33%), medulloblastoma 14/29 (48%), all tumors 34/92 (37%). Although we did not find a significant difference in the proportions of samples of various types of brain tumor tissues positive for the new splice variant (P > 0.05), levels of expression of the ACRBP-V5a splice variant were significantly different for tumor grade (P = 0.01) and tumor type (P = 0.02).

Conclusions

A newly discovered splice variant, ACRBP-V5a, is present in brain tumor. The new splicing variant may have discriminative value and potential importance in molecular-targeted therapy for brain tumors.

Dynamic changes in intron retention are tightly associated with regulation of splicing factors and proliferative activity during B-cell development

Intron retention (IR) has been proposed to modulate the delay between transcription and translation. Here, we provide an exhaustive characterization of IR in differentiated white blood cells from both the myeloid and lymphoid lineage where we observed highest levels of IR in monocytes and B-cells, in addition to previously reported granulocytes. During B-cell differentiation, we found an increase in IR from the bone marrow precursors to cells residing in secondary lymphoid organs. B-cells that undergo affinity maturation to become antibody producing plasma cells steadily decrease retention. In general, we found an inverse relationship between global IR levels and both the proliferative state of cells, and the global levels of expression of splicing factors. IR dynamics during B-cell differentiation appear to be conserved between human and mouse, suggesting that IR plays an important biological role, evolutionary conserved, during blood cell differentiation. By correlating the expression of non-core splicing factors with global IR levels, and analyzing RNA binding protein knockdown and eCLIP data, we identify a few splicing factors likely playing an evolutionary conserved role in IR regulation. Our work provides new insights into the role of IR during hematopoiesis, and on the main factors involved in regulating IR.

Effect of expression alteration in flanking genes on phenotypes of St8sia2-deficient mice

ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 2 (ST8SIA2) synthesizes polysialic acid (PSA), which is essential for brain development. Although previous studies reported that St8sia2-deficient mice that have a mixed 129 and C57BL/6 (B6) genetic background showed mild and variable phenotypes, the reasons for this remain unknown. We hypothesized that this phenotypic difference is caused by diversity in the expression or function of flanking genes of St8sia2. A genomic polymorphism and gene expression analysis in the flanking region revealed reduced expression of insulin-like growth factor 1 receptor (Igf1r) on the B6 background than on that of the 129 strain. This observation, along with the finding that administration of an IGF1R agonist during pregnancy increased litter size, suggests that the decreased expression of Igf1r associated with ST8SIA2 deficiency caused lethality. This study demonstrates the importance of gene expression level in the flanking regions of a targeted null allele having an effect on phenotype.

Identification of genes underlying phenotypic plasticity of wing size via insulin signaling pathway by network-based analysis in Sogatella furcifera

Background

Phenotypic plasticity is a common and highly adaptive phenomenon where the same genotype produces different phenotypes in response to environmental cues. Sogatella furcifera, a migratory pest of rice exhibits wing dimorphism, is a model insect for studying phenotypic plasticity of wing size. The Insullin-PI3K-Akt-FOXO signaling pathway plays a crucial role in the manipulation of wing size in the migratory insects. However, the regulatory mechanism via the pathway involved in wing dimorphism are still unexplored.

Results

Accompanied by special alternative splicing, genes involved in muscle contraction and energy metabolism were highly expressed in the wing hinges of macropters, demonstrating their adaptation for energy-demanding long-distance flights. Based on FOXO ChIP-Seq analysis, a total of 1259 putative target genes were observed in the wing hinges, including wing morph development, flight muscle and energy metabolism genes. An integrated gene interaction network was built by combining four heterogeneous datasets, and the IIS-PI3K-Akt-FOXO pathway was clustered in a divided functional module. In total, 45 genes in the module directly interacting with the IIS-PI3K-Akt-FOXO pathway showed differential expression levels between the two wing hinges, thus are regarded as potential Insulin pathway mediated wing dimorphism related genes (IWDRGs). Of the 45 IWDRGs, 5 were selected for verification by gene knockdown experiments, and played significant roles in the insect wing size regulation.

Conclusions

We provided valuable insights on the genetic basis of wing dimorphism, and also demonstrated that network analysis is a powerful approach to identify new genes regulating wing dimorphic development via insulin signaling pathway in the migratory insect.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5793-z) contains supplementary material, which is available to authorized users.

Potential relations between post-spliced introns and mature mRNAs in the Caenorhabditis elegans genome

There are potential interactions between introns and their corresponding coding sequences (CDSs) in ribosomal protein genes that have been proposed by our group and the interactions are achieved by sequence matches between the two kinds of sequences. Here, the optimal matching relations between mature mRNAs and their corresponding introns in Caenorhabditis elegans (C.elegans) were investigated by improved Smith-Waterman local alignment software. Our results showed that the remarkably matched regions appear in the untranslated regions (UTRs) of mRNAs, especially in the 3' UTR. The optimal matched segments (OMSs) are highly organized segments. In addition, the optimal matching relations were analysed between mature mRNAs and other introns. The matching strengths in the UTRs are clearly lower than those in their corresponding introns. Our studies indicate that there are potential interactions between mature mRNAs and their corresponding introns and the post-spliced introns should have other novel functions in the gene expression process.

Origin of exon skipping-rich transcriptomes in animals driven by evolution of gene architecture

Background

Alternative splicing, particularly through intron retention and exon skipping, is a major layer of pre-translational regulation in eukaryotes. While intron retention is believed to be the most prevalent mode across non-animal eukaryotes, animals have unusually high rates of exon skipping. However, when and how this high prevalence of exon skipping evolved is unknown. Since exon skipping can greatly expand proteomes, answering these questions sheds light on the evolution of higher organismal complexity in metazoans.

Results

We used RNA-seq data to quantify exon skipping and intron retention frequencies across 65 eukaryotic species, with particular focus on early branching animals and unicellular holozoans. We found that only bilaterians have significantly increased their exon skipping frequencies compared to all other eukaryotic groups. Unlike in other eukaryotes, however, exon skipping in nearly all animals, including non-bilaterians, is strongly enriched for frame-preserving sequences, suggesting that exon skipping involvement in proteome expansion predated the increase in frequency. We also identified architectural features consistently associated with higher exon skipping rates within all studied eukaryotic genomes. Remarkably, these architectures became more prevalent during animal evolution, indicating co-evolution between genome architectures and exon skipping frequencies.

Conclusion

We suggest that the increase of exon skipping rates in animals followed a two-step process. First, exon skipping in early animals became enriched for frame-preserving events. Second, bilaterian ancestors dramatically increased their exon skipping frequencies, likely driven by the interplay between a shift in their genome architectures towards more exon definition and recruitment of frame-preserving exon skipping events to functionally diversify their cell-specific proteomes.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1499-9) contains supplementary material, which is available to authorized users.

Proteomic Evidence for In-Frame and Out-of-Frame Alternatively Spliced Isoforms in Human and Mouse

In order to find evidence for translation of alternatively spliced transcripts, especially those that result in a change in reading frame, we collected exon-skipping cases previously found by RNA-Seq and applied a computational approach to screen millions of mass spectra. These spectra came from seven human and six mouse tissues, five of which are the same between the two organisms: liver, kidney, lung, heart, and brain. Overall, we detected 4 percent of all exon-skipping events found in RNA-seq data, regardless of their effect on reading frame. The fraction of alternative isoforms detected did not differ between out-of-frame and in-frame events. Moreover, the fraction of identified alternative exon-exon junctions and constitutive junctions were similar. Together, our results suggest that both in-frame and out-of-frame translation may be actively used to regulate protein activity or localization.

Unexpected conservation of the RNA splicing apparatus in the highly streamlined genome of Galdieria sulphuraria

Background

Genome reduction in intracellular pathogens and endosymbionts is usually compensated by reliance on the host for energy and nutrients. Free-living taxa with reduced genomes must however evolve strategies for generating functional diversity to support their independent lifestyles. An emerging model for the latter case is the Rhodophyta (red algae) that comprises an ecologically widely distributed, species-rich phylum. Red algae have undergone multiple phases of significant genome reduction, including extremophilic unicellular taxa with limited nuclear gene inventories that must cope with hot, highly acidic environments.

Results

Using genomic data from eight red algal lineages, we identified 155 spliceosomal machinery (SM)-associated genes that were putatively present in the red algal common ancestor. This core SM gene set is most highly conserved in Galdieria species (150 SM genes) and underwent differing levels of gene loss in other examined red algae (53-145 SM genes). Surprisingly, the high SM conservation in Galdieria sulphuraria coincides with the enrichment of spliceosomal introns in this species (2 introns/gene) in comparison to other red algae (< 0.34 introns/gene). Spliceosomal introns in G. sulphuraria undergo alternatively splicing, including many that are differentially spliced upon changes in culture temperature.

Conclusions

Our work reveals the unique nature of G. sulphuraria among red algae with respect to the conservation of the spliceosomal machinery and introns. We discuss the possible implications of these findings in the highly streamlined genome of this free-living eukaryote.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1161-x) contains supplementary material, which is available to authorized users.

The PRMT5/WDR77 complex regulates alternative splicing through ZNF326 in breast cancer

We observed overexpression and increased intra-nuclear accumulation of the PRMT5/WDR77 in breast cancer cell lines relative to immortalized breast epithelial cells. Utilizing mass spectrometry and biochemistry approaches we identified the Zn-finger protein ZNF326, as a novel interaction partner and substrate of the nuclear PRMT5/WDR77 complex. ZNF326 is symmetrically dimethylated at arginine 175 (R175) and this modification is lost in a PRMT5 and WDR77-dependent manner. Loss of PRMT5 or WDR77 in MDA-MB-231 cells leads to defects in alternative splicing, including inclusion of A-T rich exons in target genes, a phenomenon that has previously been observed upon loss of ZNF326. We observed that the alternatively spliced transcripts of a subset of these genes, involved in proliferation and tumor cell migration like REPIN1/AP4, ST3GAL6, TRNAU1AP and PFKM are degraded upon loss of PRMT5. In summary, we have identified a novel mechanism through which the PRMT5/WDR77 complex maintains the balance between splicing and mRNA stability through methylation of ZNF326.

Transcriptome analysis of hypoxic cancer cells uncovers intron retention in EIF2B5 as a mechanism to inhibit translation

Cells adjust to hypoxic stress within the tumor microenvironment by downregulating energy-consuming processes including translation. To delineate mechanisms of cellular adaptation to hypoxia, we performed RNA-Seq of normoxic and hypoxic head and neck cancer cells. These data revealed a significant down regulation of genes known to regulate RNA processing and splicing. Exon-level analyses classified > 1,000 mRNAs as alternatively spliced under hypoxia and uncovered a unique retained intron (RI) in the master regulator of translation initiation, EIF2B5. Notably, this intron was expressed in solid tumors in a stage-dependent manner. We investigated the biological consequence of this RI and demonstrate that its inclusion creates a premature termination codon (PTC), that leads to a 65kDa truncated protein isoform that opposes full-length eIF2Bε to inhibit global translation. Furthermore, expression of 65kDa eIF2Bε led to increased survival of head and neck cancer cells under hypoxia, providing evidence that this isoform enables cells to adapt to conditions of low oxygen. Additional work to uncover -cis and -trans regulators of EIF2B5 splicing identified several factors that influence intron retention in EIF2B5: a weak splicing potential at the RI, hypoxia-induced expression and binding of the splicing factor SRSF3, and increased binding of total and phospho-Ser2 RNA polymerase II specifically at the intron retained under hypoxia. Altogether, these data reveal differential splicing as a previously uncharacterized mode of translational control under hypoxia and are supported by a model in which hypoxia-induced changes to cotranscriptional processing lead to selective retention of a PTC-containing intron in EIF2B5.

An Ultraconserved Element (UCE) controls homeostatic splicing of ARGLU1 mRNA

Arginine and Glutamate-Rich protein 1 (ARGLU1) is a protein whose function is poorly understood, but may act in both transcription and pre-mRNA splicing. We demonstrate that the ARGLU1 gene expresses at least three distinct RNA splice isoforms – a fully spliced isoform coding for the protein, an isoform containing a retained intron that is detained in the nucleus, and an isoform containing an alternative exon that targets the transcript for nonsense mediated decay. Furthermore, ARGLU1 contains a long, highly evolutionarily conserved sequence known as an Ultraconserved Element (UCE) that is within the retained intron and overlaps the alternative exon. Manipulation of the UCE, in a reporter minigene or via random mutations in the genomic context using CRISPR/Cas9, changed the splicing pattern. Further, overexpression of the ARGLU1 protein shifted the splicing of endogenous ARGLU1 mRNA, resulting in an increase in the retained intron isoform and nonsense mediated decay susceptible isoform and a decrease in the fully spliced isoform. Taken together with data showing that functional protein knockout shifts splicing toward the fully spliced isoform, our data are consistent with a model in which unproductive splicing complexes assembled at the alternative exon lead to inefficient splicing and intron retention.

Correspondence of D. melanogaster and C. elegans developmental stages revealed by alternative splicing characteristics of conserved exons

Background

We report a statistical study to find correspondence of D. melanogaster and C. elegans developmental stages based on alternative splicing (AS) characteristics of conserved cassette exons using modENCODE RNA-seq data. We identify “stage-associated exons” to capture the AS characteristics of each stage and use these exons to map pairwise stages within and between the two species by an overlap test.

Results

Within fly and worm, adjacent developmental stages are mapped to each other, i.e., a strong diagonal pattern is observed as expected, supporting the validity of our approach. Between fly and worm, two parallel mapping patterns are observed between fly early embryos to early larvae and worm life cycle, and between fly late larvae to adults and worm late embryos to adults. We also apply this approach to compare tissues and cells from fly and worm. Findings include the high similarity between fly/worm adults and fly/worm embryos, groupings of fly cell lines, and strong mappings of fly head tissues to worm late embryos and male adults. Gene ontology and KEGG enrichment analyses provide a detailed functional annotation of the identified stage-associated exons, as well as a functional explanation of the observed correspondence map between fly and worm developmental stages.

Conclusions

Our results suggest that AS dynamics of the exon pairs that share similar DNA sequences are informative for finding transcriptomic similarity of biological samples. Our study is innovative in two aspects. First, to our knowledge, our study is the first comprehensive study of AS events in fly and worm developmental stages, tissues, and cells. AS events provide an alternative perspective of transcriptome dynamics, compared to gene expression events. Second, our results do not entirely rely on the information of orthologous genes. Interesting results are also observed for fly and worm cassette exon pairs with DNA sequence similarity but not in orthologous gene pairs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3600-2) contains supplementary material, which is available to authorized users.

IRFinder: assessing the impact of intron retention on mammalian gene expression

Intron retention (IR) occurs when an intron is transcribed into pre-mRNA and remains in the final mRNA. We have developed a program and database called IRFinder to accurately detect IR from mRNA sequencing data. Analysis of 2573 samples showed that IR occurs in all tissues analyzed, affects over 80% of all coding genes and is associated with cell differentiation and the cell cycle. Frequently retained introns are enriched for specific RNA binding protein sites and are often retained in clusters in the same gene. IR is associated with lower protein levels and intron-retaining transcripts that escape nonsense-mediated decay are not actively translated.

Application of Monoclonal Antibodies Against Mouse Dermokine

Dermokine is one of the most highly expressed proteins in differentiating keratinocytes. Mouse dermokine has been reported to be encoded by 22 exons, and its expression leads to three transcripts, β, γ, and α, which are transcribed from two different transcriptional start sites. The α isoform represents the carboxyl-terminal domain of the β isoform, whereas the γ isoform lacks this domain. To reveal the distributions and expression levels of each isoform in mice, we generated rat monoclonal antibodies against dermokine-β/γ and dermokine-β/α. In immunofluorescence studies, the expression levels of dermokine in the cytosol of the cultured mouse keratinocytes were significantly elevated by high levels of extracellular calcium. In Western blot analyses, the expression levels of dermokine-β and dermokine-α were increased in the presence of high calcium. Finally, we developed a monoclonal antibody-based sensitive sandwich enzyme-linked immunosorbent assay (ELISA) and showed that the secreted dermokine-β into the culture medium from mouse keratinocytes was significantly increased in a manner dependent on the extracellular calcium concentration. These dermokine-specific antibodies have allowed us to gain new insights into the role of each dermokine isoform in cutaneous homeostasis.

Evolutionarily conserved autoregulation of alternative pre-mRNA splicing by ribosomal protein L10a

Alternative splicing of pre-mRNAs can regulate expression of protein-coding genes by generating unproductive mRNAs rapidly degraded by nonsense-mediated mRNA decay (NMD). Many of the genes directly regulated by alternative splicing coupled with NMD (AS-NMD) are related to RNA metabolism, but the repertoire of genes regulated by AS-NMD in vivo is to be determined. Here, we analyzed transcriptome data of wild-type and NMD-defective mutant strains of the nematode worm Caenorhabditis elegans and demonstrate that eight of the 82 cytoplasmic ribosomal protein (rp) genes generate unproductively spliced mRNAs. Knockdown of any of the eight rp genes exerted a dynamic and compensatory effect on alternative splicing of its own transcript and inverse effects on that of the other rp genes. A large subunit protein L10a, termed RPL-1 in nematodes, directly and specifically binds to an evolutionarily conserved 39-nt stretch termed L10ARE between the two alternative 5′ splice sites in its own pre-mRNA to switch the splice site choice. Furthermore, L10ARE-mediated splicing autoregulation of the L10a-coding gene is conserved in vertebrates. These results indicate that L10a is an evolutionarily conserved splicing regulator and that homeostasis of a subset of the rp genes are regulated at the level of pre-mRNA splicing in vivo.

A proteogenomic approach to understand splice isoform functions through sequence and expression-based computational modeling

The products of multi-exon genes are a mixture of alternatively spliced isoforms, from which the translated proteins can have similar, different or even opposing functions. It is therefore essential to differentiate and annotate functions for individual isoforms. Computational approaches provide an efficient complement to expensive and time-consuming experimental studies. The input data of these methods range from DNA sequence, to RNA selection pressure, to expressed sequence tags, to full-length complementary DNA, to exon array, to RNA-seq expression, to proteomic data. Notably, RNA-seq technology generates quantitative profiling of transcript expression at the genome scale, with an unprecedented amount of expression data available for developing isoform function prediction methods. Integrative analysis of these data at different molecular levels enables a proteogenomic approach to systematically interrogate isoform functions. Here, we briefly review the state-of-the-art methods according to their input data sources, discuss their advantages and limitations and point out potential ways to improve prediction accuracies.

Correlated Evolution of Nucleotide Positions within Splice Sites in Mammals

Splice sites (SSs)--short nucleotide sequences flanking introns--are under selection for spliceosome binding, and adhere to consensus sequences. However, non-consensus nucleotides, many of which probably reduce SS performance, are frequent. Little is known about the mechanisms maintaining such apparently suboptimal SSs. Here, we study the correlations between strengths of nucleotides occupying different positions of the same SS. Such correlations may arise due to epistatic interactions between positions (i.e., a situation when the fitness effect of a nucleotide in one position depends on the nucleotide in another position), their evolutionary history, or to other reasons. Within both the intronic and the exonic parts of donor SSs, nucleotides that increase (decrease) SS strength tend to co-occur with other nucleotides increasing (respectively, decreasing) it, consistent with positive epistasis. Between the intronic and exonic parts of donor SSs, the correlations of nucleotide strengths tend to be negative, consistent with negative epistasis. In the course of evolution, substitutions at a donor SS tend to decrease the strength of its exonic part, and either increase or do not change the strength of its intronic part. In acceptor SSs, the situation is more complicated; the correlations between adjacent positions appear to be driven mainly by avoidance of the AG dinucleotide which may cause aberrant splicing. In summary, both the content and the evolution of SSs is shaped by a complex network of interdependences between adjacent nucleotides that respond to a range of sometimes conflicting selective constraints.

On the path to genetic novelties: insights from programmed DNA elimination and RNA splicing

Understanding how genetic novelties arise is a central goal of evolutionary biology. To this end, programmed DNA elimination and RNA splicing deserve special consideration. While programmed DNA elimination reshapes genomes by eliminating chromatin during organismal development, RNA splicing rearranges genetic messages by removing intronic regions during transcription. Small RNAs help to mediate this class of sequence reorganization, which is not error-free. It is this imperfection that makes programmed DNA elimination and RNA splicing excellent candidates for generating evolutionary novelties. Leveraging a number of these two processes' mechanistic and evolutionary properties, which have been uncovered over the past years, we present recently proposed models and empirical evidence for how splicing can shape the structure of protein-coding genes in eukaryotes. We also chronicle a number of intriguing similarities between the processes of programmed DNA elimination and RNA splicing, and highlight the role that the variation in the population-genetic environment may play in shaping their target sequences.

Regulation of gene expression through production of unstable mRNA isoforms

Alternative splicing is universally accredited for expanding the information encoded within the transcriptome. In recent years, several tightly regulated alternative splicing events have been reported which do not lead to generation of protein products, but lead to unstable mRNA isoforms. Instead these transcripts are targets for NMD (nonsense-mediated decay) or retained in the nucleus and degraded. In the present review I discuss the regulation of these events, and how many have been implicated in control of gene expression that is instrumental to a number of developmental paradigms. I further discuss their relevance to disease settings and conclude by highlighting technologies that will aid identification of more candidate events in future.

SNPs in APOBEC3 cytosine deaminases and their association with Visna/Maedi disease progression

The Apolipoprotein B mRNA-editing catalytic polypeptide-like 3 (APOBEC3) genes are able to inhibit the replication of a wide range of exogenous retroviruses, as well as endogenous retroviruses and retrotransposons. Three APOBEC3 genes, named APOBEC3Z1, APOBEC3Z2 and APOBEC3Z3, have been described in sheep. In this work the three genes have been screened in order to identify polymorphisms. No polymorphism was detected for the A3Z2 and A3Z3 genes but 16 SNPs and a 3-bp deletion were found in the A3Z1 gene. A thermoestability prediction analysis was applied to the detected amino acidic SNPs by three different programs. This analysis revealed a number of polymorphisms that could affect the protein stability. The SNPs of the 3'UTR were tested to detect alterations on the predicted microRNA target sites. Two new microRNA target sites were discovered for one of the alleles. Two SNPs were selected for association studies in relation with the retroviral disease Visna/Maedi in Latxa and Assaf sheep breeds. Although association analyses resulted unconclusive, probably due to the unsuitability of the SNP allele frequency distribution of the selected polymorphisms in the analyzed breeds, these genes remain good candidates for association studies.

Building robust transcriptomes with master splicing factors

Coherent splicing networks arise from many discrete splicing decisions regulated in unison. Here, we examine the properties of robust, context-specific splicing networks. We propose that a subset of key splicing regulators, or "master splicing factors," respond to environmental cues to establish and maintain tissue transcriptomes during development.

A comprehensive survey of non-canonical splice sites in the human transcriptome

We uncovered the diversity of non-canonical splice sites at the human transcriptome using deep transcriptome profiling. We mapped a total of 3.7 billion human RNA-seq reads and developed a set of stringent filters to avoid false non-canonical splice site detections. We identified 184 splice sites with non-canonical dinucleotides and U2/U12-like consensus sequences. We selected 10 of the herein identified U2/U12-like non-canonical splice site events and successfully validated 9 of them via reverse transcriptase-polymerase chain reaction and Sanger sequencing. Analyses of the 184 U2/U12-like non-canonical splice sites indicate that 51% of them are not annotated in GENCODE. In addition, 28% of them are conserved in mouse and 76% are involved in alternative splicing events, some of them with tissue-specific alternative splicing patterns. Interestingly, our analysis identified some U2/U12-like non-canonical splice sites that are converted into canonical splice sites by RNA A-to-I editing. Moreover, the U2/U12-like non-canonical splice sites have a differential distribution of splicing regulatory sequences, which may contribute to their recognition and regulation. Our analysis provides a high-confidence group of U2/U12-like non-canonical splice sites, which exhibit distinctive features among the total human splice sites.

Complex tissue-specific patterns and distribution of multiple RAGE splice variants in different mammals

The receptor for advanced glycosylation end products (RAGE) is a multiligand receptor involved in diverse cell signaling pathways. Previous studies show that this gene expresses several splice variants in human, mouse, and dog. Alternative splicing (AS) plays an important role in expanding transcriptomic and proteomic diversity, and it has been related to disease. AS is also one of the main evolutionary mechanisms in mammalian genomes. However, limited information is available regarding the AS of RAGE in a wide context of mammalian tissues. In this study, we examined in detail the different RAGE mRNAs generated by AS from six mammals, including two primates (human and monkey), two artiodactyla (cow and pig), and two rodentia (mouse and rat) in 6–18 different tissues including fetal, adult, and tumor. By nested reverse transcription-polymerase chain reaction (RT-PCR) we identified a high number of splice variants including noncoding transcripts and predicted coding ones with different potential protein modifications affecting mainly the transmembrane and ligand-binding domains that could influence their biological function. However, analysis of RNA-seq data enabled detecting only the most abundant splice variants. More than 80% of the detected RT-PCR variants (87 of 101 transcripts) are novel (different exon/intron structure to the previously described ones), and interestingly, 20–60% of the total transcripts (depending on the species) are noncoding ones that present tissue specificity. Our results suggest that RAGE undergoes extensive AS in mammals, with different expression patterns among adult, fetal, and tumor tissues. Moreover, most splice variants seem to be species specific, especially the noncoding variants, with only two (canonical human Tv1-RAGE, and human N-truncated or Tv10-RAGE) conserved among the six different species. This could indicate a special evolution pattern of this gene at mRNA level.

Origin of spliceosomal introns and alternative splicing

In this work we review the current knowledge on the prehistory, origins, and evolution of spliceosomal introns. First, we briefly outline the major features of the different types of introns, with particular emphasis on the nonspliceosomal self-splicing group II introns, which are widely thought to be the ancestors of spliceosomal introns. Next, we discuss the main scenarios proposed for the origin and proliferation of spliceosomal introns, an event intimately linked to eukaryogenesis. We then summarize the evidence that suggests that the last eukaryotic common ancestor (LECA) had remarkably high intron densities and many associated characteristics resembling modern intron-rich genomes. From this intron-rich LECA, the different eukaryotic lineages have taken very distinct evolutionary paths leading to profoundly diverged modern genome structures. Finally, we discuss the origins of alternative splicing and the qualitative differences in alternative splicing forms and functions across lineages.

Evolution of splicing regulatory networks in Drosophila

The proteome expanding effects of alternative pre-mRNA splicing have had a profound impact on eukaryotic evolution. The events that create this diversity can be placed into four major classes: exon skipping, intron retention, alternative 5′ splice sites, and alternative 3′ splice sites. Although the regulatory mechanisms and evolutionary pressures among alternative splicing classes clearly differ, how these differences affect the evolution of splicing regulation remains poorly characterized. We used RNA-seq to investigate splicing differences in D. simulans, D. sechellia, and three strains of D. melanogaster. Regulation of exon skipping and tandem alternative 3′ splice sites (NAGNAGs) were more divergent than other splicing classes. Splicing regulation was most divergent in frame-preserving events and events in noncoding regions. We further determined the contributions of cis- and trans-acting changes in splicing regulatory networks by comparing allele-specific splicing in F1 interspecific hybrids, because differences in allele-specific splicing reflect changes in cis-regulatory element activity. We find that species-specific differences in intron retention and alternative splice site usage are primarily attributable to changes in cis-regulatory elements (median ∼80% cis), whereas species-specific exon skipping differences are driven by both cis- and trans-regulatory divergence (median ∼50% cis). These results help define the mechanisms and constraints that influence splicing regulatory evolution and show that networks regulating the four major classes of alternative splicing diverge through different genetic mechanisms. We propose a model in which differences in regulatory network architecture among classes of alternative splicing affect the evolution of splicing regulation.

Alternative pre-mRNA splicing

Alternative pre-mRNA splicing is an integral part of gene regulation in eukaryotes. Here we provide a basic overview of the various types of alternative splicing, as well as the functional role, highlighting how alternative splicing varies across phylogeny. Regulated alternative splicing can affect protein function and ultimately impact biological outcomes. We examine the possibility that portions of alternatively spliced transcripts are the result of stochastic processes rather than regulated. We discuss the implications of misregulated alternative splicing and explore of the role of alternative splicing in human disease.

Nonsense-mediated decay of alternative precursor mRNA splicing variants is a major determinant of the Arabidopsis steady state transcriptome

The nonsense-mediated decay (NMD) surveillance pathway can recognize erroneous transcripts and physiological mRNAs, such as precursor mRNA alternative splicing (AS) variants. Currently, information on the global extent of coupled AS and NMD remains scarce and even absent for any plant species. To address this, we conducted transcriptome-wide splicing studies using Arabidopsis thaliana mutants in the NMD factor homologs UP FRAMESHIFT1 (UPF1) and UPF3 as well as wild-type samples treated with the translation inhibitor cycloheximide. Our analyses revealed that at least 17.4% of all multi-exon, protein-coding genes produce splicing variants that are targeted by NMD. Moreover, we provide evidence that UPF1 and UPF3 act in a translation-independent mRNA decay pathway. Importantly, 92.3% of the NMD-responsive mRNAs exhibit classical NMD-eliciting features, supporting their authenticity as direct targets. Genes generating NMD-sensitive AS variants function in diverse biological processes, including signaling and protein modification, for which NaCl stress-modulated AS-NMD was found. Besides mRNAs, numerous noncoding RNAs and transcripts derived from intergenic regions were shown to be NMD responsive. In summary, we provide evidence for a major function of AS-coupled NMD in shaping the Arabidopsis transcriptome, having fundamental implications in gene regulation and quality control of transcript processing.

The relationship between gene isoform multiplicity, number of exons and protein divergence

At present we know that phenotypic differences between organisms arise from a variety of sources, like protein sequence divergence, regulatory sequence divergence, alternative splicing, etc. However, we do not have yet a complete view of how these sources are related. Here we address this problem, studying the relationship between protein divergence and the ability of genes to express multiple isoforms. We used three genome-wide datasets of human-mouse orthologs to study the relationship between isoform multiplicity co-occurrence between orthologs (the fact that two orthologs have more than one isoform) and protein divergence. In all cases our results showed that there was a monotonic dependence between these two properties. We could explain this relationship in terms of a more fundamental one, between exon number of the largest isoform and protein divergence. We found that this last relationship was present, although with variations, in other species (chimpanzee, cow, rat, chicken, zebrafish and fruit fly). In summary, we have identified a relationship between protein divergence and isoform multiplicity co-occurrence and explained its origin in terms of a simple gene-level property. Finally, we discuss the biological implications of these findings for our understanding of inter-species phenotypic differences.

First functional analysis of a novel splicing mutation in the B3GALTL gene by an ex vivo approach in Tunisian patients with typical Peters plus syndrome

Peters plus syndrome is a rare recessive autosomal disorder comprising ocular anterior segment dysgenesis, short stature, hand abnormalities and distinctive facial features. It was related only to mutations in the B3GALTL gene in the 13q12.3 region. In this study, we undertook the first functional analysis of a novel c.597-2 A>G splicing mutation within the B3GALTL gene using an ex-vivo approach. The results showed a complete skipping of exon 8 in the B3GALTL cDNA, which altered the open reading frame of the mutant transcript and generated a PTC within exon 9. This finding potentially elicits the nonsense mRNA to degradation by NMD (nonsense-mediated mRNA decay). The theoretical consequences of splice site mutations, predicted with the bioinformatics tool Human Splice Finder, were investigated and evaluated in relation to ex-vivo results. The findings confirmed the key role played by the B3GALTL gene in typical Peters-plus syndromes and the utility of mRNA analysis to understand the primary impacts of this mutation and the phenotype of the disease.

Regulation of gene expression in mammalian nervous system through alternative pre-mRNA splicing coupled with RNA quality control mechanisms

Eukaryotic gene expression is orchestrated on a genome-wide scale through several post-transcriptional mechanisms. Of these, alternative pre-mRNA splicing expands the proteome diversity and modulates mRNA stability through downstream RNA quality control (QC) pathways including nonsense-mediated decay (NMD) of mRNAs containing premature termination codons and nuclear retention and elimination (NRE) of intron-containing transcripts. Although originally identified as mechanisms for eliminating aberrant transcripts, a growing body of evidence suggests that NMD and NRE coupled with deliberate changes in pre-mRNA splicing patterns are also used in a number of biological contexts for deterministic control of gene expression. Here we review recent studies elucidating molecular mechanisms and biological significance of these gene regulation strategies with a specific focus on their roles in nervous system development and physiology. This article is part of a Special Issue entitled 'RNA and splicing regulation in neurodegeneration'.

Evolutionary dynamics of gene and isoform regulation in Mammalian tissues

Most mammalian genes produce multiple distinct messenger RNAs through alternative splicing, but the extent of splicing conservation is not clear. To assess tissue-specific transcriptome variation across mammals, we sequenced complementary DNA from nine tissues from four mammals and one bird in biological triplicate, at unprecedented depth. We find that while tissue-specific gene expression programs are largely conserved, alternative splicing is well conserved in only a subset of tissues and is frequently lineage-specific. Thousands of previously unknown, lineage-specific, and conserved alternative exons were identified; widely conserved alternative exons had signatures of binding by MBNL, PTB, RBFOX, STAR, and TIA family splicing factors, implicating them as ancestral mammalian splicing regulators. Our data also indicate that alternative splicing often alters protein phosphorylatability, delimiting the scope of kinase signaling.

Introns in UTRs: why we should stop ignoring them

Although introns in 5'- and 3'-untranslated regions (UTRs) are found in many protein coding genes, rarely are they considered distinctive entities with specific functions. Indeed, mammalian transcripts with 3'-UTR introns are often assumed nonfunctional because they are subject to elimination by nonsense-mediated decay (NMD). Nonetheless, recent findings indicate that 5'- and 3'-UTR intron status is of significant functional consequence for the regulation of mammalian genes. Therefore these features should be ignored no longer.

Identification of androgen receptor variants in testis from humans and other vertebrates

The androgen receptor (AR) is a ligand-activated transcription factor member of the nuclear receptor superfamily. The existence of alternatively spliced variants is well recognised for several members of this superfamily, most of them having functional importance. For example, several testicular oestrogen receptor variants have been suggested to play a role in the regulation of spermatogenesis. However, information on AR variants is mostly related to cancer and androgen insensitivity syndrome (AIS) cases. The objective of this study was to investigate the expression of AR variants in the testis from humans and other vertebrates. Four AR variants [ARΔ2(Stop) , ARΔ2(23Stop) , ARΔ3 and ARΔ4(120)] were identified in human testis. ARΔ2(Stop) and ARΔ3, with exon 2 or 3 deleted, respectively, were also expressed in human liver, lung, kidney and heart. In addition, ARΔ2(Stop) was expressed in rat and gilthead seabream testis, while an ARΔ3 was detected in African clawed frog testis. This is the first report revealing the existence of AR variants in the testis of evolutionarily distant vertebrate species and in nonpathological tissues. These data suggest the functional importance of these novel AR forms and demonstrate a complexity in AR signalling that is not exclusive of pathological conditions.

RNA secondary structure mediates alternative 3'ss selection in Saccharomyces cerevisiae

Alternative splicing is the mechanism by which different combinations of exons in the pre-mRNA give rise to distinct mature mRNAs. This process is mediated by splicing factors that bind the pre-mRNA and affect the recognition of its splicing signals. Saccharomyces species lack many of the regulatory factors present in metazoans. Accordingly, it is generally assumed that the amount of alternative splicing is limited. However, there is recent compelling evidence that yeast have functional alternative splicing, mainly in response to environmental conditions. We have previously shown that sequence and structure properties of the pre-mRNA could explain the selection of 3' splice sites (ss) in Saccharomyces cerevisiae. In this work, we extend our previous observations to build a computational classifier that explains most of the annotated 3'ss in the CDS and 5' UTR of this organism. Moreover, we show that the same rules can explain the selection of alternative 3'ss. Experimental validation of a number of predicted alternative 3'ss shows that their usage is low compared to annotated 3'ss. The majority of these alternative 3'ss introduce premature termination codons (PTCs), suggesting a role in expression regulation. Furthermore, a genome-wide analysis of the effect of temperature, followed by experimental validation, yields only a small number of changes, indicating that this type of regulation is not widespread. Our results are consistent with the presence of alternative 3'ss selection in yeast mediated by the pre-mRNA structure, which can be responsive to external cues, like temperature, and is possibly related to the control of gene expression.

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.

Altered expression of dermokine in skin disorders

BACKGROUND

Although dermokine-β, a glycoprotein expressed in epithelial cells, does not have significant homology to other proteins, its carboxyl-terminal domain shares a high pI value with many cytokines, suggesting similar functions.

OBJECTIVE

To better understand the biology of dermokine, we here determined its localization under pathological conditions and examined factors that regulate its expression.

METHODS

We generated an anti-human dermokine-β/γ monoclonal antibody cross-reacting with the mouse protein. Using this antibody, immunohistological staining and Western blotting of dermokine-β/γ were performed with various tissue samples.

RESULTS

Although human dermokine-β/γ was expressed in almost all granular layers, upper spinous layers of the skin were also stained with anti-dermokine-β/γ antibody in inflammatory skin disorders. Dermokine-β/γ was expressed in keratoacanthoma and a part of well-differentiated squamous cell carcinoma (SCC). However, dermokine-β/γ was not detected in poorly differentiated SCC or tumours derived from non-keratinocytes. In mice, dermokine-β/γ-expressed keratinocytes were increased in models of contact hypersensitivity, ultraviolet-irradiated skin injury and wound healing. Consistent with expanded distribution in inflammatory skin diseases, proinflammatory cytokines such as interleukin-1β, interleukin-12, and tumour necrosis factor-α augmented dermokine-β/γ expression in cultured human keratinocytes. In contrast, growth factors including epidermal growth factor, insulin-like growth factor-I, keratinocyte growth factor and transforming growth factor-α significantly reduced dermokine expression.

CONCLUSION

These results provide novel insights into the physiological and pathological significance of dermokine in the epidermis.

Alternative splicing and nonsense-mediated decay modulate expression of important regulatory genes in Arabidopsis

Alternative splicing (AS) coupled to nonsense-mediated decay (NMD) is a post-transcriptional mechanism for regulating gene expression. We have used a high-resolution AS RT-PCR panel to identify endogenous AS isoforms which increase in abundance when NMD is impaired in the Arabidopsis NMD factor mutants, upf1-5 and upf3-1. Of 270 AS genes (950 transcripts) on the panel, 102 transcripts from 97 genes (32%) were identified as NMD targets. Extrapolating from these data around 13% of intron-containing genes in the Arabidopsis genome are potentially regulated by AS/NMD. This cohort of naturally occurring NMD-sensitive AS transcripts also allowed the analysis of the signals for NMD in plants. We show the importance of AS in introns in 5' or 3'UTRs in modulating NMD-sensitivity of mRNA transcripts. In particular, we identified upstream open reading frames overlapping the main start codon as a new trigger for NMD in plants and determined that NMD is induced if 3'-UTRs were >350 nt. Unexpectedly, although many intron retention transcripts possess NMD features, they are not sensitive to NMD. Finally, we have shown that AS/NMD regulates the abundance of transcripts of many genes important for plant development and adaptation including transcription factors, RNA processing factors and stress response genes.

Arabidopsis plants having defects in nonsense-mediated mRNA decay factors UPF1, UPF2, and UPF3 show photoperiod-dependent phenotypes in development and stress responses

Nonsense-mediated mRNA decay (NMD) is an important mRNA quality surveillance pathway in all eukaryotes that eliminates aberrant mRNAs derived from various sources. Three NMD factor proteins, UPF1, UPF2, and UPF3 are required for the NMD process and were found to be also involved in certain stress responses in mammalian and yeast cells. Using Arabidopsis thaliana mutants of UPF1 and UPF3 and UPF2-silenced lines (irUPF2), we examined the involvement of UPF1, UPF2, and UPF3 in development and in response to stresses, wounding and infection by Pseudomonas syringae pv. tomato strain DC3000. Under the long (16 h) photoperiod condition, Arabidopsis having a defect in NMD factors exhibited altered morphologies of various organs, disturbed homeostasis of wounding-induced jasmonic acid and pathogen-elicited salicylic acid, and abnormal wounding- and methyl jasmonate-induced changes in the transcript levels of two defense-related genes, LOX2 and VSP2. Importantly, when plants were cultivated under the short (10 h) photoperiod condition, mutants of UPF1 and UPF3 and irUPF2 showed smaller differences from the wild-type plants in growth and stress-induced responses. These data suggest a complex regulatory network, likely composed of light signaling and NMD factor-mediated pathways, in influencing plant development and adaption to environmental stresses.

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ∼15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.