The Meaning of NMD: Translate or Perish

Sample Size Estimation for Detection of Splicing Events in Transcriptome Sequencing Data

Merging data from multiple samples is required to detect low expressed transcripts or splicing events that might be present only in a subset of samples. However, the exact number of required replicates enabling the detection of such rare events often remains a mystery but can be approached through probability theory. Here, we describe a probabilistic model, relating the number of observed events in a batch of samples with observation probabilities. Therein, samples appear as a heterogeneous collection of events, which are observed with some probability. The model is evaluated in a batch of 54 transcriptomes of human dermal fibroblast samples. The majority of putative splice-sites (alignment gap-sites) are detected in (almost) all samples or only sporadically, resulting in an U-shaped pattern for observation probabilities. The probabilistic model systematically underestimates event numbers due to a bias resulting from finite sampling. However, using an additional assumption, the probabilistic model can predict observed event numbers within a <10% deviation from the median. Single samples contain a considerable amount of uniquely observed putative splicing events (mean 7122 in alignments from TopHat alignments and 86,215 in alignments from STAR). We conclude that the probabilistic model provides an adequate description for observation of gap-sites in transcriptome data. Thus, the calculation of required sample sizes can be done by application of a simple binomial model to sporadically observed random events. Due to the large number of uniquely observed putative splice-sites and the known stochastic noise in the splicing machinery, it appears advisable to include observation of rare splicing events into analysis objectives. Therefore, it is beneficial to take scores for the validation of gap-sites into account.

Regulation of gene expression by translation factor eIF5A: Hypusine-modified eIF5A enhances nonsense-mediated mRNA decay in human cells

Nonsense-mediated mRNA decay (NMD) couples protein synthesis to mRNA turnover. It eliminates defective transcripts and controls the abundance of certain normal mRNAs. Our study establishes a connection between NMD and the translation factor eIF5A (eukaryotic initiation factor 5A) in human cells. eIF5A modulates the synthesis of groups of proteins (the eIF5A regulon), and undergoes a distinctive two-step post-translational modification (hypusination) catalyzed by deoxyhypusine synthase and deoxyhypusine hydroxylase. We show that expression of NMD-susceptible constructs was increased by depletion of the major eIF5A isoform, eIF5A1. NMD was also attenuated when hypusination was inhibited by RNA interference with either of the two eIF5A modifying enzymes, or by treatment with the drugs ciclopirox or deferiprone which inhibit deoxyhypusine hydroxylase. Transcriptome analysis by RNA-Seq identified human genes whose expression is coordinately regulated by eIF5A1, its modifying enzymes, and the pivotal NMD factor, Upf1. Transcripts encoding components of the translation system were highly represented, including some encoding ribosomal proteins controlled by alternative splicing coupled to NMD (AS-NMD). Our findings extend and strengthen the association of eIF5A with NMD, previously inferred in yeast, and show that hypusination is important for this function of human eIF5A. In addition, they advance drug-mediated NMD suppression as a therapeutic opportunity for nonsense-associated diseases. We propose that regulation of mRNA stability contributes to eIF5A's role in selective gene expression.

Human nonsense-mediated RNA decay regulates EMT by targeting the TGF-ß signaling pathway in lung adenocarcinoma

Nonsense-mediated mRNA decay (NMD) is a highly conserved pathway that selectively degrades aberrant RNA transcripts. In this study, we proved that NMD regulates the epithelial-mesenchymal transition (EMT) of lung adenocarcinoma (ADC). Moreover, we found that NMD core factor UP-frameshift 1 tends to be expressed at lower levels in human ADC tissues than in normal lung tissues, thereby raising the possibility that NMD may be downregulated to permit ADC oncogenesis. Our experiments in human ADC cell lines showed that downregulating NMD can promote EMT. Moreover, EMT can be inhibited by upregulating NMD. We tested the role of TGF-ß signaling and found that NMD influences EMT by targeting the TGF-ß signaling pathway. Our findings reveal that NMD is a potential tumor regulatory mechanism and may be a potential therapeutic target for ADC.

Nonsense-mediated mRNA decay in Tetrahymena is EJC independent and requires a protozoa-specific nuclease

Nonsense-mediated mRNA decay (NMD) is essential for removing premature termination codon-containing transcripts from cells. Studying the NMD pathway in model organisms can help to elucidate the NMD mechanism in humans and improve our understanding of how this biologically important process has evolved. Ciliates are among the earliest branching eukaryotes; their NMD mechanism is poorly understood and may be primordial. We demonstrate that highly conserved Upf proteins (Upf1a, Upf2 and Upf3) are involved in the NMD pathway of the ciliate, Tetrahymena thermophila. We further show that a novel protozoa-specific nuclease, Smg6L, is responsible for destroying many NMD-targeted transcripts. Transcriptome-wide identification and characterization of NMD-targeted transcripts in vegetative Tetrahymena cells showed that many have exon-exon junctions downstream of the termination codon. However, Tetrahymena may lack a functional exon junction complex (EJC), and the Tetrahymena ortholog of an EJC core component, Mago nashi (Mag1), is dispensable for NMD. Therefore, NMD is EJC independent in this early branching eukaryote.

Estimating the prevalence of functional exonic splice regulatory information

In addition to coding information, human exons contain sequences necessary for correct splicing. These elements are known to be under purifying selection and their disruption can cause disease. However, the density of functional exonic splicing information remains profoundly uncertain. Several groups have experimentally investigated how mutations at different exonic positions affect splicing. They have found splice information to be distributed widely in exons, with one estimate putting the proportion of splicing-relevant nucleotides at >90%. These results suggest that splicing could place a major pressure on exon evolution. However, analyses of sequence conservation have concluded that the need to preserve splice regulatory signals only slightly constrains exon evolution, with a resulting decrease in the average human rate of synonymous evolution of only 1–4%. Why do these two lines of research come to such different conclusions? Among other reasons, we suggest that the methods are measuring different things: one assays the density of sites that affect splicing, the other the density of sites whose effects on splicing are visible to selection. In addition, the experimental methods typically consider short exons, thereby enriching for nucleotides close to the splice junction, such sites being enriched for splice-control elements. By contrast, in part owing to correction for nucleotide composition biases and to the assumption that constraint only operates on exon ends, the conservation-based methods can be overly conservative.

Frequent GU wobble pairings reduce translation efficiency in Plasmodium falciparum

Plasmodium falciparum genome has 81% A+T content. This nucleotide bias leads to extreme codon usage bias and culminates in frequent insertion of asparagine homorepeats in the proteome. Using recodonized GFP sequences, we show that codons decoded via G:U wobble pairing are suboptimal codons that are negatively associated to protein translation efficiency. Despite this, one third of all codons in the genome are GU wobble codons, suggesting that codon usage in P. falciparum has not been driven to maximize translation efficiency, but may have evolved as translational regulatory mechanism. Particularly, asparagine homorepeats are generally encoded by locally clustered GU wobble AAT codons, we demonstrated that this GU wobble-rich codon context is the determining factor that causes reduction of protein level. Moreover, insertion of clustered AAT codons also causes destabilization of the transcripts. Interestingly, more frequent asparagine homorepeats insertion is seen in single-exon genes, suggesting transcripts of these genes may have been programmed for rapid mRNA decay to compensate for the inefficiency of mRNA surveillance regulation on intronless genes. To our knowledge, this is the first study that addresses P. falciparum codon usage in vitro and provides new insights on translational regulation and genome evolution of this parasite.

Pharmacological inhibition of the spliceosome subunit SF3b triggers exon junction complex-independent nonsense-mediated decay

Spliceostatin A, meayamycin, and pladienolide B are small molecules that target the SF3b subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP). These compounds are attracting much attention as tools to manipulate splicing and for use as potential anti-cancer drugs. We investigated the effects of these inhibitors on mRNA transport and stability in human cells. Upon splicing inhibition, unspliced pre-mRNAs accumulated in the nucleus, particularly within enlarged nuclear speckles. However, a small fraction of the pre-mRNA molecules were exported to the cytoplasm. We identified the export adaptor ALYREF as being associated with intron-containing transcripts and show its requirement for the nucleo-cytoplasmic transport of unspliced pre-mRNA. In contrast, the exon junction complex (EJC) core protein eIF4AIII failed to form a stable complex with intron-containing transcripts. Despite the absence of EJC, unspliced transcripts in the cytoplasm were degraded by nonsense-mediated decay (NMD), suggesting that unspliced transcripts are degraded by an EJC-independent NMD pathway. Collectively, our results indicate that although blocking the function of SF3b elicits a massive accumulation of unspliced pre-mRNAs in the nucleus, intron-containing transcripts can still bind the ALYREF export factor and be transported to the cytoplasm, where they trigger an alternative NMD pathway.

Selective inhibition of RNA polymerase I transcription as a potential approach to treat African trypanosomiasis

Trypanosoma brucei relies on an essential Variant Surface Glycoprotein (VSG) coat for survival in the mammalian bloodstream. High VSG expression within an expression site body (ESB) is mediated by RNA polymerase I (Pol I), which in other eukaryotes exclusively transcribes ribosomal RNA genes (rDNA). As T. brucei is reliant on Pol I for VSG transcription, we investigated Pol I transcription inhibitors for selective anti-trypanosomal activity. The Pol I inhibitors quarfloxin (CX-3543), CX-5461, and BMH-21 are currently under investigation for treating cancer, as rapidly dividing cancer cells are particularly dependent on high levels of Pol I transcription compared with nontransformed cells. In T. brucei all three Pol I inhibitors have IC50 concentrations for cell proliferation in the nanomolar range: quarfloxin (155 nM), CX-5461 (279 nM) or BMH-21 (134 nM) compared with IC50 concentrations in the MCF10A human breast epithelial cell line (4.44 μM, 6.89 μM or 460 nM, respectively). T. brucei was therefore 29-fold more sensitive to quarfloxin, 25-fold more sensitive to CX-5461 and 3.4-fold more sensitive to BMH-21. Cell death in T. brucei was due to rapid inhibition of Pol I transcription, as within 15 minutes treatment with the inhibitors rRNA precursor transcript was reduced 97-98% and VSG precursor transcript 91-94%. Incubation with Pol I transcription inhibitors also resulted in disintegration of the ESB as well as the nucleolus subnuclear structures, within one hour. Rapid ESB loss following the block in Pol I transcription argues that the ESB is a Pol I transcription nucleated structure, similar to the nucleolus. In addition to providing insight into Pol I transcription and ES control, Pol I transcription inhibitors potentially also provide new approaches to treat trypanosomiasis.

Trypanosoma brucei is protected by an essential Variant Surface Glycoprotein (VSG) coat in the mammalian bloodstream. The active VSG gene is transcribed by RNA polymerase I (Pol I), which typically only transcribes rDNA. Pol I transcription inhibitors are under clinical trials for cancer chemotherapy. As T. brucei relies on Pol I for VSG transcription, we investigated its susceptibility to these drugs. We show that quarfloxin (CX-3543), CX-5461, and BMH-21 are effective against T. brucei at nanomolar concentrations. T. brucei death was due to rapid and specific inhibition of Pol I transcription. Incubation with Pol I transcription inhibitors also resulted in disappearance of Pol I subnuclear structures like the nucleolus and the VSG expression site body (ESB). Rapid ESB loss followed the Pol I transcription block, arguing that the ESB is nucleated by Pol I transcription. Pol I transcription inhibitors could therefore potentially function as novel drugs against trypanosomiasis.

Exon junction complex proteins bind nascent transcripts independently of pre-mRNA splicing in Drosophila melanogaster

Although it is currently understood that the exon junction complex (EJC) is recruited on spliced mRNA by a specific interaction between its central protein, eIF4AIII, and splicing factor CWC22, we found that eIF4AIII and the other EJC core proteins Y14 and MAGO bind the nascent transcripts of not only intron-containing but also intronless genes on Drosophila polytene chromosomes. Additionally, Y14 ChIP-seq demonstrates that association with transcribed genes is also splicing-independent in Drosophila S2 cells. The association of the EJC proteins with nascent transcripts does not require CWC22 and that of Y14 and MAGO is independent of eIF4AIII. We also show that eIF4AIII associates with both polysomal and monosomal RNA in S2 cell extracts, whereas Y14 and MAGO fractionate separately. Cumulatively, our data indicate a global role of eIF4AIII in gene expression, which would be independent of Y14 and MAGO, splicing, and of the EJC, as currently understood.