Exon-intron boundary inhibits m6A deposition, enabling m6A distribution hallmark, longer mRNA half-life and flexible protein coding.

Regional bias of N6-methyladenosine (m6A) mRNA modification avoiding splice site region, calls for an open hypothesis whether exon-intron boundary could affect m6A deposition. By deep learning modeling, we find that exon-intron boundary represses a proportion (12% to 34%) of m6A deposition at adjacent exons (~100 nt to splice site). Experiments validate that m6A signal increases once the host gene does not undergo pre-mRNA splicing to produce the same mRNA. Inhibited m6A sites have higher m6A enhancers and lower m6A silencers locally and show high heterogeneity at different exons genome-widely, with only a small proportion (12% to 15%) of exons showing strong inhibition, enabling more stable mRNAs and flexible protein coding. m6A is majorly responsible for why mRNAs with more exons be more stable. Exon junction complex (EJC) only partially contributes to this exon-intron boundary m6A inhibition in some short internal exons, highlighting additional factors yet to be identified.

Deep learning modeling m6A deposition reveals the importance of downstream cis-element sequences.

The N6-methyladenosine (m6A) modification is deposited to nascent transcripts on chromatin, but its site-specificity mechanism is mostly unknown. Here we model the m6A deposition to pre-mRNA by iM6A (intelligent m6A), a deep learning method, demonstrating that the site-specific m6A methylation is primarily determined by the flanking nucleotide sequences. iM6A accurately models the m6A deposition (AUROC = 0.99) and uncovers surprisingly that the cis-elements regulating the m6A deposition preferentially reside within the 50 nt downstream of the m6A sites. The m6A enhancers mostly include part of the RRACH motif and the m6A silencers generally contain CG/GT/CT motifs. Our finding is supported by both independent experimental validations and evolutionary conservation. Moreover, our work provides evidences that mutations resulting in synonymous codons can affect the m6A deposition and the TGA stop codon favors m6A deposition nearby. Our iM6A deep learning modeling enables fast paced biological discovery which would be cost-prohibitive and unpractical with traditional experimental approaches, and uncovers a key cis-regulatory mechanism for m6A site-specific deposition.

Saturation mutagenesis reveals manifold determinants of exon definition.

To illuminate the extent and roles of exonic sequences in the splicing of human RNA transcripts, we conducted saturation mutagenesis of a 51-nt internal exon in a three-exon minigene. All possible single and tandem dinucleotide substitutions were surveyed. Using high-throughput genetics, 5560 minigene molecules were assayed for splicing in human HEK293 cells. Up to 70% of mutations produced substantial (greater than twofold) phenotypes of either increased or decreased splicing. Of all predicted secondary structural elements, only a single 15-nt stem-loop showed a strong correlation with splicing, acting negatively. The in vitro formation of exon-protein complexes between the mutant molecules and proteins associated with spliceosome formation (U2AF35, U2AF65, U1A, and U1-70K) correlated with splicing efficiencies, suggesting exon definition as the step affected by most mutations. The measured relative binding affinities of dozens of human RNA binding protein domains as reported in the CISBP-RNA database were found to correlate either positively or negatively with splicing efficiency, more than could fit on the 51-nt test exon simultaneously. The large number of these functional protein binding correlations point to a dynamic and heterogeneous population of pre-mRNA molecules, each responding to a particular collection of binding proteins.

Pre-mRNA processing includes N6 methylation of adenosine residues that are retained in mRNA exons and the fallacy of "RNA epigenetics".

By using a cell fraction technique that separates chromatin-associated nascent RNA, newly completed nucleoplasmic mRNA and cytoplasmic mRNA, we have shown in a previous study that residues in exons are methylated (m6A) in nascent pre-mRNA and remain methylated in the same exonic residues in nucleoplasmic and cytoplasmic mRNA. Thus, there is no evidence of a substantial degree of demethylation in mRNA exons that would correspond to so-called "epigenetic" demethylation. The turnover rate of mRNA molecules is faster, depending on m6A content in HeLa cell mRNA, suggesting that specification of mRNA stability may be the major role of m6A exon modification. In mouse embryonic stem cells (mESCs) lacking Mettl3, the major mRNA methylase, the cells continue to grow, making the same mRNAs with unchanged splicing profiles in the absence (>90%) of m6A in mRNA, suggesting no common obligatory role of m6A in splicing. All these data argue strongly against a commonly used "reversible dynamic methylation/demethylation" of mRNA, calling into question the concept of "RNA epigenetics" that parallels the well-established role of dynamic DNA epigenetics.

m6A mRNA modifications are deposited in nascent pre-mRNA and are not required for splicing but do specify cytoplasmic turnover.

Understanding the biologic role of N6-methyladenosine (m6A) RNA modifications in mRNA requires an understanding of when and where in the life of a pre-mRNA transcript the modifications are made. We found that HeLa cell chromatin-associated nascent pre-mRNA (CA-RNA) contains many unspliced introns and m6A in exons but very rarely in introns. The m6A methylation is essentially completed upon the release of mRNA into the nucleoplasm. Furthermore, the content and location of each m6A modification in steady-state cytoplasmic mRNA are largely indistinguishable from those in the newly synthesized CA-RNA or nucleoplasmic mRNA. This result suggests that quantitatively little methylation or demethylation occurs in cytoplasmic mRNA. In addition, only ∼10% of m6As in CA-RNA are within 50 nucleotides of 5' or 3' splice sites, and the vast majority of exons harboring m6A in wild-type mouse stem cells is spliced the same in cells lacking the major m6A methyltransferase Mettl3. Both HeLa and mouse embryonic stem cell mRNAs harboring m6As have shorter half-lives, and thousands of these mRNAs have increased half-lives (twofold or more) in Mettl3 knockout cells compared with wild type. In summary, m6A is added to exons before or soon after exon definition in nascent pre-mRNA, and while m6A is not required for most splicing, its addition in the nascent transcript is a determinant of cytoplasmic mRNA stability.

A majority of m6A residues are in the last exons, allowing the potential for 3' UTR regulation.

We adapted UV CLIP (cross-linking immunoprecipitation) to accurately locate tens of thousands of m(6)A residues in mammalian mRNA with single-nucleotide resolution. More than 70% of these residues are present in the 3'-most (last) exons, with a very sharp rise (sixfold) within 150-400 nucleotides of the start of the last exon. Two-thirds of last exon m(6)A and >40% of all m(6)A in mRNA are present in 3' untranslated regions (UTRs); contrary to earlier suggestions, there is no preference for location of m(6)A sites around stop codons. Moreover, m(6)A is significantly higher in noncoding last exons than in next-to-last exons harboring stop codons. We found that m(6)A density peaks early in the 3' UTR and that, among transcripts with alternative polyA (APA) usage in both the brain and the liver, brain transcripts preferentially use distal polyA sites, as reported, and also show higher proximal m(6)A density in the last exons. Furthermore, when we reduced m6A methylation by knocking down components of the methylase complex and then examined 661 transcripts with proximal m6A peaks in last exons, we identified a set of 111 transcripts with altered (approximately two-thirds increased proximal) APA use. Taken together, these observations suggest a role of m(6)A modification in regulating proximal alternative polyA choice.

Quantitative evaluation of all hexamers as exonic splicing elements.

We describe a comprehensive quantitative measure of the splicing impact of a complete set of RNA 6-mer sequences by deep sequencing successfully spliced transcripts. All 4096 6-mers were substituted at five positions within two different internal exons in a 3-exon minigene, and millions of successfully spliced transcripts were sequenced after transfection of human cells. The results allowed the assignment of a relative splicing strength score to each mutant molecule. The effect of 6-mers on splicing often depended on their location; much of this context effect could be ascribed to the creation of different overlapping sequences at each site. Taking these overlaps into account, the splicing effect of each 6-mer could be quantified, and 6-mers could be designated as enhancers (ESEseqs) and silencers (ESSseqs), with an ESRseq score indicating their strength. Some 6-mers exhibited positional bias relative to the two splice sites. The distribution and conservation of these ESRseqs in and around human exons supported their classification. Predicted RNA secondary structure effects were also seen: Effective enhancers, silencers and 3' splice sites tend to be single stranded, and effective 5' splice sites tend to be double stranded. 6-mers that may form positive or negative synergy with another were also identified. Chromatin structure may also influence the splicing enhancement observed, as a good correspondence was found between splicing performance and the predicted nucleosome occupancy scores of 6-mers. This approach may prove of general use in defining nucleic acid regulatory motifs, substitute for functional SELEX in most cases, and provide insights about splicing mechanisms.

Context-dependent splicing regulation: exon definition, co-occurring motif pairs and tissue specificity.

Splicing is a crucial process in gene expression in higher organisms because: 1) most vertebrate genes contain introns; and 2) alternative splicing is primarily responsible for increasing proteomic complexity and functional diversity. Intron definition, the coordination across an intron, is a mandatory step in the splicing process. However, exon definition, the coordination across an exon, is also thought to be required for the splicing of most vertebrate exons. Recent investigations of exon definition complexes provide insights into splicing dynamics. That splicing regulators act in a context-dependent mode is supported by a large collection of evidence. Splicing contexts generally can be classified as cis-element and trans-element contexts. A widespread cis-element context is defined by co-occurring motif pairs to which splicing regulatory factors bind to direct specific molecular interactions. Splicing regulation is also coordinated by trans-element contexts as exemplified by tissue specific splicing, where alternative exons can be coordinately regulated by a few splicing factors, the expression and/or activity of which are concertedly higher or lower in the corresponding tissues.

Intronic motif pairs cooperate across exons to promote pre-mRNA splicing.

BACKGROUND: A very early step in splice site recognition is exon definition, a process that is as yet poorly understood. Communication between the two ends of an exon is thought to be required for this step. We report genome-wide evidence for exons being defined through the combinatorial activity of motifs located in flanking intronic regions. RESULTS: Strongly co-occurring motifs were found to specifically reside in four intronic regions surrounding a large number of human exons. These paired motifs occur around constitutive and alternative exons but not pseudo exons. Most co-occurring motifs are limited to intronic regions within 100 nucleotides of the exon. They are preferentially associated with weaker exons. Their pairing is conserved in evolution and they exhibit a lower frequency of single nucleotide polymorphism when paired. Paired motifs display specificity with respect to distance from the exon borders and in constitutive versus alternative splicing. Many resemble binding sites for heterogeneous nuclear ribonucleoproteins. Specific pairs are associated with tissue-specific genes, the higher expression of which coincides with that of the pertinent RNA binding proteins. Tested pairs acted synergistically to enhance exon inclusion, and this enhancement was found to be exon-specific. CONCLUSIONS: The exon-flanking sequence pairs identified here by genomic analysis promote exon inclusion and may play a role in the exon definition step in pre-mRNA splicing. We propose a model in which multiple concerted interactions are required between exonic sequences and flanking intronic sequences to effect exon definition.

Splicing of designer exons reveals unexpected complexity in pre-mRNA splicing.

Pre-messengerRNA (mRNA) splicing requires the accurate recognition of splice sites by the cellular RNA processing machinery. In addition to sequences that comprise the branchpoint and the 3' and 5' splice sites, the cellular splicing machinery relies on additional information in the form of exonic and intronic splicing enhancer and silencer sequences. The high abundance of these motifs makes it difficult to investigate their effects using standard genetic perturbations, since their disruption often leads to the formation of yet new elements. To lessen this problem, we have designed synthetic exons comprised of multiple copies of a single prototypical exonic enhancer and a single prototypical exonic silencer sequence separated by neutral spacer sequences. The spacer sequences buffer the exon against the formation of new elements as the number and order of the original elements are varied. Over 100 such designer exons were constructed by random ligation of enhancer, silencer, and neutral elements. Each exon was positioned as the central exon in a 3-exon minigene and tested for exon inclusion after transient transfection. The level of inclusion of the test exons was seen to be dependent on the provision of enhancers and could be decreased by the provision of silencers. In general, there was a good quantitative correlation between the proportion of enhancers and splicing. However, widely varying inclusion levels could be produced by different permutations of the enhancer and silencer elements, indicating that even in this simplified system splicing decisions rest on complex interplays of yet to be determined parameters.

Activated expression of an Arabidopsis HD-START protein confers drought tolerance with improved root system and reduced stomatal density.

Drought is one of the most important environmental constraints limiting plant growth and agricultural productivity. To understand the underlying mechanism of drought tolerance and to identify genes for improving this important trait, we conducted a gain-of-function genetic screen for improved drought tolerance in Arabidopsis thaliana. One mutant with improved drought tolerance was isolated and designated as enhanced drought tolerance1. The mutant has a more extensive root system than the wild type, with deeper roots and more lateral roots, and shows a reduced leaf stomatal density. The mutant had higher levels of abscisic acid and Pro than the wild type and demonstrated an increased resistance to oxidative stress and high levels of superoxide dismutase. Molecular genetic analysis and recapitulation experiments showed that the enhanced drought tolerance is caused by the activated expression of a T-DNA tagged gene that encodes a putative homeodomain-START transcription factor. Moreover, overexpressing the cDNA of the transcription factor in transgenic tobacco also conferred drought tolerance associated with improved root architecture and reduced leaf stomatal density. Therefore, we have revealed functions of the homeodomain-START factor that were gained upon altering its expression pattern by activation tagging and provide a key regulator that may be used to improve drought tolerance in plants.

Positive selection acting on splicing motifs reflects compensatory evolution.

We have used comparative genomics to characterize the evolutionary behavior of predicted splicing regulatory motifs. Using base substitution rates in intronic regions as a calibrator for neutral change, we found a strong avoidance of synonymous substitutions that disrupt predicted exonic splicing enhancers or create predicted exonic splicing silencers. These results attest to the functionality of the hexameric motif set used and suggest that they are subject to purifying selection. We also found that synonymous substitutions in constitutive exons tend to create exonic splicing enhancers and to disrupt exonic splicing silencers, implying positive selection for these splicing promoting events. We present evidence that this positive selection is the result of splicing-positive events compensating for splicing-negative events as well as for mutations that weaken splice-site sequences. Such compensatory events include nonsynonymous mutations, synonymous mutations, and mutations at splice sites. Compensation was also seen from the fact that orthologous exons tend to maintain the same number of predicted splicing motifs. Our data fit a splicing compensation model of exon evolution, in which selection for splicing-positive mutations takes place to counter the effect of an ongoing splicing-negative mutational process, with the exon as a whole being conserved as a unit of splicing. In the course of this analysis, we observed that synonymous positions in general are conserved relative to intronic sequences, suggesting that messenger RNA molecules are rich in sequence information for functions beyond protein coding and splicing.