Introns as relict retrotransposons: implications for the evolutionary origin of eukaryotic mRNA splicing mechanisms

The reference genome and full-length transcriptome of pakchoi provide insights into cuticle formation and heat adaption

Brassica rapa includes various vegetables with high economic value. Among them, green petiole type pakchoi (B. rapa ssp. chinensis) is one of the major vegetables grown in southern China. Compared with other B. rapa varieties, green petiole type pakchoi shows a higher level of heat resistance, which is partially derived from the rich epicuticular wax. Here we sequence a high-quality genome of green petiole type pakchoi, which has been widely used as the parent in breeding. Our results reveal that long terminal repeat retrotransposon insertion plays critical roles in promoting the genome expansion and transcriptional diversity of pakchoi genes through preferential insertions, particularly in cuticle biosynthetic genes. After whole-genome triplication, over-retained pakchoi genes escape stringent selection pressure, and among them a set of cuticle-related genes are retained. Using bulked-segregant analysis of a heat-resistant pakchoi cultivar, we identify a frame-shift deletion across the third exon and the subsequent intron of BrcCER1 in candidate regions. Using Nanopore long-read sequencing, we analyze the full-length transcriptome of two pakchoi cultivars with opposite sensitivity to high temperature. We find that the heat-resistant pakchoi cultivar can mitigate heat-caused leaf damage by activating an unfolded protein response, as well as by inhibiting chloroplast development and energy metabolism, which are presumably mediated by both transcriptional regulation and splicing factors. Our study provides valuable resources for Brassica functional genomics and breeding research, and deepens our understanding of plant stress resistance.

LTR-retrotransposons in plants: Engines of evolution

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.

Deconstructing the dogma: a new view of the evolution and genetic programming of complex organisms

Since the birth of molecular biology it has been generally assumed that most genetic information is transacted by proteins, and that RNA plays an intermediary role. This led to the subsidiary assumption that the vast tracts of noncoding sequences in the genomes of higher organisms are largely nonfunctional, despite the fact that they are transcribed. These assumptions have since become articles of faith, but they are not necessarily correct. I propose an alternative evolutionary history whereby developmental and cognitive complexity has arisen by constructing sophisticated RNA-based regulatory networks that interact with generic effector complexes to control gene expression patterns and the epigenetic trajectories of differentiation and development. Environmental information can also be conveyed into this regulatory system via RNA editing, especially in the brain. Moreover, the observations that RNA-directed epigenetic changes can be inherited raises the intriguing question: has evolution learnt how to learn?

Signatures of domain shuffling in the human genome

To elucidate the role of exon shuffling in shaping the complexity of the human genome/proteome, we have systematically analyzed intron phase distributions in the coding sequence of human protein domains. We found that introns at the boundaries of domains show high excess of symmetrical phase combinations (i.e., 0-0, 1-1, and 2-2), whereas nonboundary introns show no excess symmetry. This suggests that exon shuffling has primarily involved rearrangement of structural and functional domains as a whole. Furthermore, we found that domains flanked by phase 1 introns have dramatically expanded in the human genome due to domain shuffling and that 1-1 symmetrical domains and domain families are nonrandomly distributed with respect to their age. The predominance and extracellular location of 1-1 symmetrical domains among domains specific to metazoans suggests that they are associated with the rise of multicellularity. On the other hand, 0-0 symmetrical domains tend to be over-represented among ancient protein domains that are shared between the eukaryotic and prokaryotic kingdoms, which is compatible with the suggestion of primordial domain shuffling in the progenote. To see whether the human data reflect general genomic patterns of metazoans, similar analyses were done for the nematode Caenorhabditis elegans. Although the C. elegans data generally concur with the human patterns, we identified fewer intron-bounded domains in this organism, consistent with the lower complexity of C. elegans genes. [The following individuals kindly provided reagents, samples, or unpublished information as indicated in the paper: Z. Gu and R. Stevens.]

Features of spliceosome evolution and function inferred from an analysis of the information at human splice sites

An information analysis of the 5' (donor) and 3' (acceptor) sequences spanning the ends of nearly 1800 human introns has provided evidence for structural features of splice sites that bear upon spliceosome evolution and function: (1) 82% of the sequence information (i.e. sequence conservation) at donor junctions and 97% of the sequence information at acceptor junctions is confined to the introns, allowing codon choices throughout exons to be largely unrestricted. The distribution of information at intron-exon junctions is also described in detail and compared with footprints. (2) Acceptor sites are found to possess enough information to be located in the transcribed portion of the human genome, whereas donor sites possess about one bit less than the information needed to locate them independently. This difference suggests that acceptor sites are located first in humans and, having been located, reduce by a factor of two the number of alternative sites available as donors. Direct experimental evidence exists to support this conclusion. (3) The sequences of donor and acceptor splice sites exhibit a striking similarity. This suggests that the two junctions derive from a common ancestor and that during evolution the information of both sites shifted onto the intron. If so, the protein and RNA components that are found in contemporary spliceosomes, and which are responsible for recognizing donor and acceptor sequences, should also be related. This conclusion is supported by the common structures found in different parts of the spliceosome.

Statistical analysis and prediction of the exonic structure of human genes

Nonhomologous fully sequenced human protein-coding genes were studied. Three sets of exon-exon junctions were formed defined by the intron (shadow) position relative to the reading frame. For the analysis of intron shadow signals in exons, information content and discrimination energy approaches were used with the correction allowing one to ignore the influence of a protein-coding message. The corrected formulas allow one to define the consensuses for the three types of intron shadow signals as a AG/guwn, cAG/GUnn, and cAG/gunU, and provide better recognition than the original formulas. The analysis of the codon usage in the signal positions leads to the conclusion that the prevalence of some amino acids in corresponding protein sites is caused by the signal requirements and not vice versa. The distribution of potential intron shadow signals in exons contradicts the hypothesis of intron insertion into suitable preexisting sites. There exists a correlation between the intron types and/or the exon length modulo 3.

Evolutionary dynamics of transposable elements in prokaryotes and eukaryotes

This paper summarizes some recent theories about the evolution of transposable genetic elements in outbreeding, sexual eukaryotic organisms. The evolutionary possibilities available to self-replicating transposable elements are shown to vary depending on the reproductive biology of the host genome. This effect can be used to explain, in part, the differences in abundance of transposable elements between prokaryotes and eukaryotes. It is argued that the pattern of sexual outbreeding seen in mammals and plants is especially favorable to the spread of transposons. Moreover, because transposon spread is facilitated by zygote formation, the evolutionary origin of sexual conjugation may have been due to selection on transposon-encoded genes. Finally, evidence is also presented that introns could have originated as transposable genetic elements.

How big is the universe of exons?

If genes have been assembled from exon subunits, the frequency with which exons are reused leads to an estimate of the size of the underlying exon universe. An exon database was constructed from available protein sequences, and homologous exons were identified on the basis of amino acid identity; statistically significant matches were determined by Monte Carlo methods. It is estimated that only 1000 to 7000 exons were needed to construct all proteins.

Self-splicing introns as a source for transposable genetic elements

Previous theories have suggested that some introns with the ability to self-splice are derived from transposable elements. However, an interpretation is given here that suggests retrotransposons and retroviruses (transposable elements which move via RNA intermediates) have evolved from self-splicing introns. This is based on the involvement of RNA intermediates, the ancestral nature of the self-splicing reaction, and the assumed presence of introns in an RNA world. Conserved sequences within the introns, essential for splicing, and their wide phylogenetic distribution also make it unlikely that they are descended from transposable elements. Mitochondrial plasmids of Neurospora species containing features of both introns and retrotransposons have a central role in the resolution of the problem and are considered here to support the view that introns are, or have been, sources of mobile elements. The possibility of other transposable elements arising from introns is also considered.

A general model for the evolution of nuclear pre-mRNA introns

We present an overview of the evolution of eukaryotic split gene structure and pre-mRNA splicing mechanisms. We have drawn together several seemingly conflicting ideas and we show that they can all be incorporated in a single unified theory of intron evolution. The resulting model is consistent with the notion that introns, as a class, are very ancient, having originated in the "RNA world"; it also supports the concept that introns may have played a crucial role in the construction of many eukaryotic genes and it accommodates the idea that introns are related to mobile insertion elements. Our conclusion is that introns could have a profound effect on the course of eukaryotic gene evolution, but that the origin and maintenance of intron sequences depends, largely, on natural selection acting on the intron sequences themselves.