SERpredict: detection of tissue- or tumor-specific isoforms generated through exonization of transposable elements

Tissue-specific usage of transposable element-derived promoters in mouse development

BACKGROUND

Transposable elements (TEs) are a significant component of eukaryotic genomes and play essential roles in genome evolution. Mounting evidence indicates that TEs are highly transcribed in early embryo development and contribute to distinct biological functions and tissue morphology.

RESULTS

We examine the epigenetic dynamics of mouse TEs during the development of five tissues: intestine, liver, lung, stomach, and kidney. We found that TEs are associated with over 20% of open chromatin regions during development. Close to half of these accessible TEs are only activated in a single tissue and a specific developmental stage. Most accessible TEs are rodent-specific. Across these five tissues, 453 accessible TEs are found to create the transcription start sites of downstream genes in mouse, including 117 protein-coding genes and 144 lincRNA genes, 93.7% of which are mouse-specific. Species-specific TE-derived transcription start sites are found to drive the expression of tissue-specific genes and change their tissue-specific expression patterns during evolution.

CONCLUSION

Our results suggest that TE insertions increase the regulatory potential of the genome, and some TEs have been domesticated to become a crucial component of gene and regulate tissue-specific expression during mouse tissue development.

A single mutation in the ACTR8 gene associated with lineage-specific expression in primates

Background

Alternative splicing (AS) generates various transcripts from a single gene and thus plays a significant role in transcriptomic diversity and proteomic complexity. Alu elements are primate-specific transposable elements (TEs) and can provide a donor or acceptor site for AS. In a study on TE-mediated AS, we recently identified a novel AluSz6-exonized ACTR8 transcript of the crab-eating monkey (Macaca fascicularis). In the present study, we sought to determine the molecular mechanism of AluSz6 exonization of the ACTR8 gene and investigate its evolutionary and functional consequences in the crab-eating monkey.

Results

We performed RT-PCR and genomic PCR to analyze AluSz6 exonization in the ACTR8 gene and the expression of the AluSz6-exonized transcript in nine primate samples, including prosimians, New world monkeys, Old world monkeys, and hominoids. AluSz6 integration was estimated to have occurred before the divergence of simians and prosimians. The Alu-exonized transcript obtained by AS was lineage-specific and expressed only in Old world monkeys and apes, and humans. This lineage-specific expression was caused by a single G duplication in AluSz6, which provides a new canonical 5′ splicing site. We further identified other alternative transcripts that were unaffected by the AluSz6 insertion. Finally, we observed that the alternative transcripts were transcribed into new isoforms with C-terminus deletion, and in silico analysis showed that these isoforms do not have a destructive function.

Conclusions

The single G duplication in the TE sequence is the source of TE exonization and AS, and this mutation may suffer a different fate of ACTR8 gene expression during primate evolution.

An atlas of transposable element-derived alternative splicing in cancer

Transposable element (TE)-derived sequences comprise more than half of the human genome, and their presence has been documented to alter gene expression in a number of different ways, including the generation of alternatively spliced transcript isoforms. Alternative splicing has been associated with tumorigenesis for a number of different cancers. The objective of this study was to broadly characterize the role of human TEs in generating alternatively spliced transcript isoforms in cancer. To do so, we screened for the presence of TE-derived sequences co-located with alternative splice sites that are differentially used in normal versus cancer tissues. We analysed a comprehensive set of alternative splice variants characterized for 614 matched normal-tumour tissue pairs across 13 cancer types, resulting in the discovery of 4820 TE-generated alternative splice events distributed among 723 cancer-associated genes. Short interspersed nuclear elements (Alu) and long interspersed nuclear elements (L1) were found to contribute the majority of TE-generated alternative splice sites in cancer genes. A number of cancer-associated genes, including MYH11, WHSC1 and CANT1, were shown to have overexpressed TE-derived isoforms across a range of cancer types. TE-derived isoforms were also linked to cancer-specific fusion transcripts, suggesting a novel mechanism for the generation of transcriptome diversity via trans-splicing mediated by dispersed TE repeats. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

The Role of Transposable Elements in Emergence of Metazoa

Systems initially emerged for protecting genomes against insertions of transposable elements and represented by mechanisms of splicing regulation, RNA-interference, and epigenetic factors have played a key role in the evolution of animals. Many studies have shown inherited transpositions of mobile elements in embryogenesis and preservation of their activities in certain tissues of adult organisms. It was supposed that on the emergence of Metazoa the self-regulation mechanisms of transposons related with the gene networks controlling their activity could be involved in intercellular cell coordination in the cascade of successive divisions with differentiated gene expression for generation of tissues and organs. It was supposed that during evolution species-specific features of transposons in the genomes of eukaryotes could form the basis for creation of dynamically related complexes of systems for epigenetic regulation of gene expression. These complexes could be produced due to the influence of noncoding transposon-derived RNAs on DNA methylation, histone modifications, and processing of alternative splicing variants, whereas the mobile elements themselves could be directly involved in the regulation of gene expression in cis and in trans. Transposons are widely distributed in the genomes of eukaryotes; therefore, their activation can change the expression of specific genes. In turn, this can play an important role in cell differentiation during ontogenesis. It is supposed that transposons can form a species-specific pattern for control of gene expression, and that some variants of this pattern can be favorable for adaptation. The presented data indicate the possible influence of transposons in karyotype formation. It is supposed that transposon localization relative to one another and to protein-coding genes can influence the species-specific epigenetic regulation of ontogenesis.

Alu-Derived Alternative Splicing Events Specific to Macaca Lineages in CTSF Gene

Cathepsin F, which is encoded by CTSF, is a cysteine proteinase ubiquitously expressed in several tissues. In a previous study, novel transcripts of the CTSF gene were identified in the crab-eating monkey deriving from the integration of an Alu element-AluYRa1. The occurrence of AluYRa1-derived alternative transcripts and the mechanism of exonization events in the CTSF gene of human, rhesus monkey, and crab-eating monkey were investigated using PCR and reverse transcription PCR on the genomic DNA and cDNA isolated from several tissues. Results demonstrated that AluYRa1 was only integrated into the genome of Macaca species and this lineage-specific integration led to exonization events by producing a conserved 3' splice site. Six transcript variants (V1-V6) were generated by alternative splicing (AS) events, including intron retention and alternative 5' splice sites in the 5' and 3' flanking regions of CTSF_AluYRa1. Among them, V3-V5 transcripts were ubiquitously expressed in all tissues of rhesus monkey and crab-eating monkey, whereas AluYRa1-exonized V1 was dominantly expressed in the testis of the crab-eating monkey, and V2 was only expressed in the testis of the two monkeys. These five transcript variants also had different amino acid sequences in the C-terminal region of CTSF, as compared to reference sequences. Thus, species-specific Alu-derived exonization by lineage-specific integration of Alu elements and AS events seems to have played an important role during primate evolution by producing transcript variants and gene diversification.

Dual function of C/D box small nucleolar RNAs in rRNA modification and alternative pre-mRNA splicing

C/D box small nucleolar RNAs (SNORDs) are small noncoding RNAs, and their best-understood function is to target the methyltransferase fibrillarin to rRNA (for example, SNORD27 performs 2'-O-methylation of A27 in 18S rRNA). Unexpectedly, we found a subset of SNORDs, including SNORD27, in soluble nuclear extract made under native conditions, where fibrillarin was not detected, indicating that a fraction of the SNORD27 RNA likely forms a protein complex different from canonical snoRNAs found in the insoluble nuclear fraction. As part of this previously unidentified complex,SNORD27 regulates the alternative splicing of the transcription factor E2F7p re-mRNA through direct RNA-RNA interaction without methylating the RNA, likely by competing with U1 small nuclear ribonucleoprotein (snRNP). Furthermore, knockdown of SNORD27 activates previously "silent" exons in several other genes through base complementarity across the entire SNORD27 sequence, not just the antisense boxes. Thus, some SNORDs likely function in both rRNA and pre-mRNA processing, which increases the repertoire of splicing regulators and links both processes.

Transcription of LINE-derived sequences in exercise-induced stress in horses

A large proportion of mammalian genomes is represented by transposable elements (TE), most of them being long interspersed nuclear elements 1 (LINE-1 or L1). An increased expression of LINE-1 elements may play an important role in cellular stress-related conditions exerting drastic effects on the mammalian transcriptome. To understand the impact of TE on the known horse transcriptome, we masked the horse EST database, pointing out that the amount is consistent with other major vertebrates. A previously developed transcript-derived fragments (TDFs) dataset, deriving from exercise-stimulated horse peripheral blood mononuclear cells (PBMCs), was found to be enriched with L1 (26.8% in terms of bp). We investigated the involvement of TDFs in exercise-induced stress through bioinformatics and gene expression analysis. Results indicate that LINE-derived sequences are not only highly but also differentially expressed during physical effort, hinting at interesting scenarios in the regulation of gene expression in relation to exercise.

Gene length and expression level shape genomic novelties

Gene duplication and alternative splicing are important mechanisms in the production of genomic novelties. Previous work has shown that a gene’s family size and the number of splice variants it produces are inversely related, although the underlying reason is not well understood. Here, we report that gene length and expression level together explain this relationship. We found that gene lengths correlate with both gene duplication and alternative splicing: Longer genes are less likely to produce duplicates and more likely to exhibit alternative splicing. We show that gene length is a dynamic property, increasing with evolutionary time—due in part to the insertions of transposable elements—and decreasing following partial gene duplications. However, gene length alone does not account for the relationship between alternative splicing and gene duplication. A gene’s expression level appears both to impose a strong constraint on its length and to restrict gene duplications. Furthermore, high gene expression promotes alternative splicing, in particular for long genes, and alternatively, short genes with low expression levels have large gene families. Our analysis of the human and mouse genomes shows that gene length and expression level are primary genic properties that together account for the relationship between gene duplication and alternative splicing and bias the origin of novelties in the genome.

The extent of Ds1 transposon to enrich transcriptomes and proteomes by exonization

BACKGROUND

Exonization is an event which an intronic transposed element (TE) provides splice sites and leads to alternatively spliced cassette exons. Without disrupting of the inserted gene's function, TEs can expand the proteome diversity by adding the splice variant that encodes a different, yet functional protein. Previously, we found that the main contribution of Ds exonization for gene divergence is not providing genetic messages but incorporating the intron sequences with different reading frame patterns to enrich the plant proteome. Ds1, another member of Ac/Ds transposon system, differs from Ds by providing 3 splice donor sites and 2 acceptor sites for alternative splicing, which may greatly increase the extent for proteome expansion.

RESULTS

In this study, we performed a genome-wide survey of Ds1 exonization events to assess its extent to enrich proteomes in plants. Each Ds1 insertion yielded 11 transcript isoforms by integrating the splice donor and/or acceptor sites, which composed a bulk of all exonized transcript orthologs from the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice). The exonized transcripts were analyzed by the locations of the termination codon (PTC) and the putative targets for the nonsense-mediated decay (NMD) pathway were then excluded. Compared with the Ds element, Ds1 harbors more contents of non-NMD transcripts for protein isoforms.

CONCLUSIONS

The contribution of Ds1 exonization for gene divergence is incorporating the intron sequences with different reading frame patterns to enrich the plant proteome. All these simulation results direct new experimental analysis at the molecular level.

Analysis of new functional profiles of protein isoforms yielded by ds exonization in rice

Insertion of transposable elements (TEs) into introns can lead to their activation as alternatively spliced cassette exons, an event called exonization. Exonization can enrich the complexity of transcriptomes and proteomes. Previously, we performed a genome-wide computational analysis of Ds exonization events in the monocot Oryza sativa (rice). The insertion patterns of Ds increased the number of transcripts and subsequent protein isoforms, which were determined as interior and C-terminal variants. In this study, these variants were scanned with the PROSITE database in order to identify new functional profiles (domains) that were referred to their reference proteins. The new profiles of the variants were expected to be beneficial for a selective advantage and more than 70% variants achieved this. The new functional profiles could be contributed by an exon–intron junction, an intron alone, an intron–TE junction, or a TE alone. A Ds-inserted intron may yield 167 new profiles on average, while some cases can yield thousands of new profiles, of which C-terminal variants were in major. Additionally, more than 90% of the TE-inserted genes were found to gain novel functional profiles in each intron via exonization. Therefore, new functional profiles yielded by the exonization may occur in many local regions of the reference protein.

Structure and Expression Analyses of SVA Elements in Relation to Functional Genes

SINE-VNTR-Alu (SVA) elements are present in hominoid primates and are divided into 6 subfamilies (SVA-A to SVA-F) and active in the human population. Using a bioinformatic tool, 22 SVA element-associated genes are identified in the human genome. In an analysis of genomic structure, SVA elements are detected in the 5' untranslated region (UTR) of HGSNAT (SVA-B), MRGPRX3 (SVA-D), HYAL1 (SVA-F), TCHH (SVA-F), and ATXN2L (SVA-F) genes, while some elements are observed in the 3'UTR of SPICE1 (SVA-B), TDRKH (SVA-C), GOSR1 (SVA-D), BBS5 (SVA-D), NEK5 (SVA-D), ABHD2 (SVA-F), C1QTNF7 (SVA-F), ORC6L (SVA-F), TMEM69 (SVA-F), and CCDC137 (SVA-F) genes. They could contribute to exon extension or supplying poly A signals. LEPR (SVA-C), ALOX5 (SVA-D), PDS5B (SVA-D), and ABCA10 (SVA-F) genes also showed alternative transcripts by SVA exonization events. Dominant expression of HYAL1_SVA appeared in lung tissues, while HYAL1_noSVA showed ubiquitous expression in various human tissues. Expression of both transcripts (TDRKH_SVA and TDRKH_noSVA) of the TDRKH gene appeared to be ubiquitous. Taken together, these data suggest that SVA elements cause transcript isoforms that contribute to modulation of gene regulation in various human tissues.

Genome-wide survey of ds exonization to enrich transcriptomes and proteomes in plants

Insertion of transposable elements (TEs) into introns can lead to their activation as alternatively spliced cassette exons, an event called exonization which can enrich the complexity of transcriptomes and proteomes. Previously, we performed the first experimental assessment of TE exonization by inserting a Ds element into each intron of the rice epsps gene. Exonization of Ds in plants was biased toward providing splice donor sites from the beginning of the inserted Ds sequence. Additionally, Ds inserted in the reverse direction resulted in a continuous splice donor consensus region by offering 4 donor sites in the same intron. The current study involved genome-wide computational analysis of Ds exonization events in the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice). Up to 71% of the exonized transcripts were putative targets for the nonsense-mediated decay (NMD) pathway. The insertion patterns of Ds and the polymorphic splice donor sites increased the transcripts and subsequent protein isoforms. Protein isoforms contain protein sequence due to unspliced intron-TE region and/or a shift of the reading frame. The number of interior protein isoforms would be twice that of C-terminal isoforms, on average. TE exonization provides a promising way for functional expansion of the plant proteome.

Ds transposon is biased towards providing splice donor sites for exonization in transgenic tobacco

Insertion of transposed elements into introns can lead to their activation as alternatively spliced cassette exons, an event called exonization, which can enrich the complexity of transcriptomes and proteomes. In this study, the first exonization event was detected when the modified rice EPSPS marker gene was inserted with the Ac transposon 5' end, which provided a splice donor site to yield abundant novel transcripts. To assess the contribution of splice donor and acceptor sites of transposon sequences, we inserted a Ds element into each intron of the EPSPS marker gene. This process yielded 14 constructs, with the Ds transposon inserted in the forward and reverse direction in each of the 7 introns of the EPSPS marker gene. The constructs were transformed into tobacco plants, and novel transcripts were identified by RT-PCR with specific primers. Exonization of Ds in EPSPS was biased towards providing splice donor sites of the inserted Ds sequence. Additionally, when the Ds inserted in reverse direction, a continuous splice donor consensus region was determined by offering 4 donor sites in the same intron. Information on these exonization events may help enhance gene divergence and functional genomic studies.

Quantification of Hepatic UDP glucuronosyltransferase 1A splice variant expression and correlation of UDP glucuronosyltransferase 1A1 variant expression with glucuronidation activity

The UDP glucuronosyltransferase (UGT) 1A gene cluster encodes nine UGT1A family members via splicing of individual first exons to common exons 2 through 5. Each of these nine UGT1As can also undergo alternative splicing at their 3' ends by using an alternate exon 5, resulting in 27 different UGT1A mRNA species with each UGT1A gene encoding three different combinations of 5A and 5B UGT1A exons. To examine the importance of UGT1A exon 5 splice variants on overall UGT1A activity, a nested quantitative polymerase chain reaction assay was developed to accurately assess the combined expression of exon 5 splice variants (termed v2/v3) versus the expression of wild-type (termed v1) for each specific UGT1A. v1 expression was 16-, 17-, 57- and 29-fold higher than that observed for the levels of v2/v3 for UGTs 1A1, 1A4, 1A6, and 1A9, respectively, in normal human liver specimens. In a series of 58 normal human liver specimens, the expression of both UGT1A1 v1 and v2/v3 mRNAs was positively correlated with raloxifene glucuronidation activity in corresponding microsomes prepared from the same specimens (p < 0.0001, r² = 0.720; p = 0.0002, r² = 0.241, respectively), with expression of both variants lower in individuals homozygous for the UGT1A1*28 allele (42% for v1, p = 0.041; 53% for v2/v3, p = 0.0075). The expression of UGT1A1 v2/v3 was 1.6-fold higher than v1 (p = 0.03) in HepG2 cells, and short interfering RNA knockdown of HepG2 v2/v3 increased raloxifene glucuronidation activity by 83%. Together, these data suggest that hepatic UGT1A v2/v3 mRNA species are minor form variants in human livers from most individuals.

Transposable elements as tools for reshaping the genome: it is a huge world after all!

Transposable elements (TEs) are discrete pieces of DNA that can move from one site to another within genomes and sometime between genomes. They are found in all major branches of life. Because of their wide distribution and considerable diversity, they are a considerable source of genomic variation and as such, they constitute powerful drivers of genome evolution. Moreover, it is becoming clear that the epigenetic regulation of certain genes is derived from defense mechanisms against the activity of ancestral transposable elements. TEs now tend to be viewed as natural molecular tools that can reshape the genome, which challenges the idea that TEs are natural tools used to answer biological questions. In the first part of this chapter, we review the classification and distribution of TEs, and look at how they have contributed to the structural and transcriptional reshaping of genomes. In the second part, we describe methodological innovations that have modified their contribution as molecular tools.

Expressing genes do not forget their LINEs: transposable elements and gene expression

Historically the accumulated mass of mammalian transposable elements (TEs), particularly those located within gene boundaries, was viewed as a genetic burden potentially detrimental to the genomic landscape. This notion has been strengthened by the discovery that transposable sequences can alter the architecture of the transcriptome, not only through insertion, but also long after the integration process is completed. Insertions previously considered harmless are now known to impact the expression of host genes via modification of the transcript quality or quantity, transcriptional interference, or by the control of pathways that affect the mRNA life-cycle. Conversely, several examples of the evolutionary advantageous impact of TEs on the host gene structure that diversified the cellular transcriptome are reported. TE-induced changes in gene expression can be tissue- or disease-specific, raising the possibility that the impact of TE sequences may vary during development, among normal cell types, and between normal and disease-affected tissues. The understanding of the rules and abundance of TE-interference with gene expression is in its infancy, and its contribution to human disease and/or evolution remains largely unexplored.

Molecular characterization of alternative transcripts of the horse BMAL1 gene

The horse BMAL1 gene encodes the brain and muscle Arnt-like protein 1, which is a key regulator of circadian rhythmic systems in most organs and cells. The first exon of the horse-specific BMAL1 gene is produced by an exonization event of LINE3 (CR1) and SINE (MIR) was detected by bioinformatic analysis. Alternative variants generated by cassette exon event in various horse tissues were also detected by RT-PCR amplification and sequencing. The cDNA sequences of the horse transcripts (BMAL1a, BMAL1b) contain additional 21 bp and 71 bp fragments relative to horse BMAL1. Quantitative real-time RT-PCR was performed to compare the expression patterns between transcript variants in various horse tissues. The results of these experiments showed splice variants that were widely expressed in most tissues. Furthermore, they were highly expressed in cerebellum, heart, and kidney.

Widespread establishment and regulatory impact of Alu exons in human genes

The Alu element has been a major source of new exons during primate evolution. Thousands of human genes contain spliced exons derived from Alu elements. However, identifying Alu exons that have acquired genuine biological functions remains a major challenge. We investigated the creation and establishment of Alu exons in human genes, using transcriptome profiles of human tissues generated by high-throughput RNA sequencing (RNA-Seq) combined with extensive RT-PCR analysis. More than 25% of Alu exons analyzed by RNA-Seq have estimated transcript inclusion levels of at least 50% in the human cerebellum, indicating widespread establishment of Alu exons in human genes. Genes encoding zinc finger transcription factors have significantly higher levels of Alu exonization. Importantly, Alu exons with high splicing activities are strongly enriched in the 5'-UTR, and two-thirds (10/15) of 5'-UTR Alu exons tested by luciferase reporter assays significantly alter mRNA translational efficiency. Mutational analysis reveals the specific molecular mechanisms by which newly created 5'-UTR Alu exons modulate translational efficiency, such as the creation or elongation of upstream ORFs that repress the translation of the primary ORFs. This study presents genomic evidence that a major functional consequence of Alu exonization is the lineage-specific evolution of translational regulation. Moreover, the preferential creation and establishment of Alu exons in zinc finger genes suggest that Alu exonization may have globally affected the evolution of primate and human transcriptomes by regulating the protein production of master transcriptional regulators in specific lineages.

All y'all need to know 'bout retroelements in cancer

Genetic instability is one of the principal hallmarks and causative factors in cancer. Human transposable elements (TE) have been reported to cause human diseases, including several types of cancer through insertional mutagenesis of genes critical for preventing or driving malignant transformation. In addition to retrotransposition-associated mutagenesis, TEs have been found to contribute even more genomic rearrangements through non-allelic homologous recombination. TEs also have the potential to generate a wide range of mutations derivation of which is difficult to directly trace to mobile elements, including double strand breaks that may trigger mutagenic genomic rearrangements. Genome-wide hypomethylation of TE promoters and significantly elevated TE expression in almost all human cancers often accompanied by the loss of critical DNA sensing and repair pathways suggests that the negative impact of mobile elements on genome stability should increase as human tumors evolve. The biological consequences of elevated retroelement expression, such as the rate of their amplification, in human cancers remain obscure, particularly, how this increase translates into disease-relevant mutations. This review is focused on the cellular mechanisms that control human TE-associated mutagenesis in cancer and summarizes the current understanding of TE contribution to genetic instability in human malignancies.

The role of transposable elements in the evolution of non-mammalian vertebrates and invertebrates

Background

Transposable elements (TEs) have played an important role in the diversification and enrichment of mammalian transcriptomes through various mechanisms such as exonization and intronization (the birth of new exons/introns from previously intronic/exonic sequences, respectively), and insertion into first and last exons. However, no extensive analysis has compared the effects of TEs on the transcriptomes of mammals, non-mammalian vertebrates and invertebrates.

Results

We analyzed the influence of TEs on the transcriptomes of five species, three invertebrates and two non-mammalian vertebrates. Compared to previously analyzed mammals, there were lower levels of TE introduction into introns, significantly lower numbers of exonizations originating from TEs and a lower percentage of TE insertion within the first and last exons. Although the transcriptomes of vertebrates exhibit significant levels of exonization of TEs, only anecdotal cases were found in invertebrates. In vertebrates, as in mammals, the exonized TEs are mostly alternatively spliced, indicating that selective pressure maintains the original mRNA product generated from such genes.

Conclusions

Exonization of TEs is widespread in mammals, less so in non-mammalian vertebrates, and very low in invertebrates. We assume that the exonization process depends on the length of introns. Vertebrates, unlike invertebrates, are characterized by long introns and short internal exons. Our results suggest that there is a direct link between the length of introns and exonization of TEs and that this process became more prevalent following the appearance of mammals.

Characteristics of transposable element exonization within human and mouse

Insertion of transposed elements within mammalian genes is thought to be an important contributor to mammalian evolution and speciation. Insertion of transposed elements into introns can lead to their activation as alternatively spliced cassette exons, an event called exonization. Elucidation of the evolutionary constraints that have shaped fixation of transposed elements within human and mouse protein coding genes and subsequent exonization is important for understanding of how the exonization process has affected transcriptome and proteome complexities. Here we show that exonization of transposed elements is biased towards the beginning of the coding sequence in both human and mouse genes. Analysis of single nucleotide polymorphisms (SNPs) revealed that exonization of transposed elements can be population-specific, implying that exonizations may enhance divergence and lead to speciation. SNP density analysis revealed differences between Alu and other transposed elements. Finally, we identified cases of primate-specific Alu elements that depend on RNA editing for their exonization. These results shed light on TE fixation and the exonization process within human and mouse genes.

Large-scale analysis of exonized mammalian-wide interspersed repeats in primate genomes

Transposable elements (TEs) are major sources of new exons in higher eukaryotes. Almost half of the human genome is derived from TEs, and many types of TEs have the potential to exonize. In this work, we conducted a large-scale analysis of human exons derived from mammalian-wide interspersed repeats (MIRs), a class of old TEs which was active prior to the radiation of placental mammals. Using exon array data of 328 MIR-derived exons and RT-PCR analysis of 39 exons in 10 tissues, we identified 15 constitutively spliced MIR exons, and 15 MIR exons with tissue-specific shift in splicing patterns. Analysis of RNAs from multiple species suggests that the splicing events of many strongly included MIR exons have been established before the divergence of primates and rodents, while a small percentage result from recent exonization during primate evolution. Interestingly, exon array data suggest substantially higher splicing activities of MIR exons when compared with exons derived from Alu elements, a class of primate-specific retrotransposons. This appears to be a universal difference between exons derived from young and old TEs, as it is also observed when comparing Alu exons to exons derived from LINE1 and LINE2, two other groups of old TEs. Together, this study significantly expands current knowledge about exonization of TEs. Our data imply that with sufficient evolutionary time, numerous new exons could evolve beyond the evolutionary intermediate state and contribute functional novelties to modern mammalian genomes.

Diverse splicing patterns of exonized Alu elements in human tissues

Exonization of Alu elements is a major mechanism for birth of new exons in primate genomes. Prior analyses of expressed sequence tags show that almost all Alu-derived exons are alternatively spliced, and the vast majority of these exons have low transcript inclusion levels. In this work, we provide genomic and experimental evidence for diverse splicing patterns of exonized Alu elements in human tissues. Using Exon array data of 330 Alu-derived exons in 11 human tissues and detailed RT-PCR analyses of 38 exons, we show that some Alu-derived exons are constitutively spliced in a broad range of human tissues, and some display strong tissue-specific switch in their transcript inclusion levels. Most of such exons are derived from ancient Alu elements in the genome. In SEPN1, mutations of which are linked to a form of congenital muscular dystrophy, the muscle-specific inclusion of an Alu-derived exon may be important for regulating SEPN1 activity in muscle. Realtime qPCR analysis of this SEPN1 exon in macaque and chimpanzee tissues indicates human-specific increase in its transcript inclusion level and muscle specificity after the divergence of humans and chimpanzees. Our results imply that some Alu exonization events may have acquired adaptive benefits during the evolution of primate transcriptomes.

New exons have been created and added to existing functional genes during eukaryotic genome evolution. Alu elements, a class of primate-specific retrotransposons, are a major source of new exons in primates. However, recent analyses of expressed sequence tags suggest that the vast majority of Alu-derived exons are low-abundance splice forms and represent non-functional evolutionary intermediates. In order to elucidate the evolutionary impact of Alu-derived exons, we investigated the splicing of 330 Alu-derived exons in 11 human tissues using data from high-density exon arrays with multiple oligonucleotide probes for every exon in the human genome. Our exon array analysis and further RT-PCR experiments reveal surprisingly diverse splicing patterns of these exons. Some Alu-derived exons are constitutively spliced, and some are strongly tissue-specific. In SEPN1, a gene implicated in a form of congenital muscular dystrophy, our data suggest that the muscle-specific inclusion of an Alu-derived exon results from a human-specific splicing change after the divergence of humans and chimpanzees. Our study provides novel insight into the evolutionary significance of Alu exonization events. A subset of Alu-derived exons, especially those derived from more ancient Alu elements in the genome, may have contributed to functional novelties during primate evolution.

Great exaptations

Long interspersed nuclear elements (LINEs) are among the most successful parasitic genetic sequences in higher organisms. Recent work has discovered many instances of LINE incorporation into exons, reminding us of the hazards they pose to genes in their vicinity as well as their potential to be co-opted for the host's purposes.