Relevance of Extending FGFR3 Gene Analysis in Osteochondrodysplasia to Non-Coding Sequences: A Case Report

Skeletal dysplasia, also called osteochondrodysplasia, is a category of disorders affecting bone development and children's growth. Up to 552 genes, including fibroblast growth factor receptor 3 (FGFR3), have been implicated by pathogenic variations in its genesis. Frequently identified causal mutations in osteochondrodysplasia arise in the coding sequences of the FGFR3 gene: c.1138G>A and c.1138G>C in achondroplasia and c.1620C>A and c.1620C>G in hypochondroplasia. However, in some cases, the diagnostic investigations undertaken thus far have failed to identify the causal anomaly, which strengthens the relevance of the diagnostic strategies being further refined. We observed a Caucasian adult with clinical and radiographic features of achondroplasia, with no common pathogenic variant. Exome sequencing detected an FGFR3(NM_000142.4):c.1075+95C>G heterozygous intronic variation. In vitro studies showed that this variant results in the aberrant exonization of a 90-nucleotide 5' segment of intron 8, resulting in the substitution of the alanine (Ala359) for a glycine (Gly) and the in-frame insertion of 30 amino acids. This change may alter FGFR3's function. Our report provides the first clinical description of an adult carrying this variant, which completes the phenotype description previously provided in children and confirms the recurrence, the autosomal-dominant pathogenicity, and the diagnostic relevance of this FGFR3 intronic variant. We support its inclusion in routinely used diagnostic tests for osteochondrodysplasia. This may increase the detection rate of causal variants and therefore could have a positive impact on patient management. Finally, FGFR3 alteration via non-coding sequence exonization should be considered a recurrent disease mechanism to be taken into account for new drug design and clinical trial strategies.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

Evolutionary origin of the type 2 corticotropin-releasing hormone receptor γ splice variant

Many G protein-coupled receptors have splice variants, with potentially different pharmaceutical properties, expression patterns and roles. The human brain expresses three functional splice variants of the type 2 corticotropin-releasing hormone: CRHR2α, -β and -γ. CRHR2γ has only been reported in humans, but its phylogenetic distribution, and how and when during mammalian evolution it arose, is unknown. Based on genomic sequence analyses, we predict that a functional CRHR2γ is present in all Old World monkeys and apes, and is unique to these species. CRHR2γ arose by exaptation of an intronic sequence-already present in the common ancestor of primates and rodents-after retrotransposition of a short interspersed nuclear element (SINE) and mutations that created a 5' donor splice site and in-frame start codon, 32-43 million years ago. The SINE is not part of the coding sequence, only of the 5' untranslated region and may therefore play a role in translational regulation. Putative regulatory elements and an alternative transcriptional start site were added earlier to this genomic locus by a DNA transposon. The evolutionary history of CRHR2γ confirms some of the earlier reported principles behind the "birth" of alternative exons. The functional significance of CRHR2γ, particularly in the brain, remains to be showed.

Long Interleukin-22 Binding Protein Isoform-1 Is an Intracellular Activator of the Unfolded Protein Response

The human IL22RA2 gene co-produces three protein isoforms in dendritic cells [IL-22 binding protein isoform-1 (IL-22BPi1), IL-22BPi2, and IL-22BPi3]. Two of these, IL-22BPi2 and IL-22BPi3, are capable of neutralizing the biological activity of IL-22. The function of IL-22BPi1, which differs from IL-22BPi2 through an in-frame 32-amino acid insertion provided by an alternatively spliced exon, remains unknown. Using transfected human cell lines, we demonstrate that IL-22BPi1 is secreted detectably, but at much lower levels than IL-22BPi2, and unlike IL-22BPi2 and IL-22BPi3, is largely retained in the endoplasmic reticulum (ER). As opposed to IL-22BPi2 and IL-22BPi3, IL-22BPi1 is incapable of neutralizing or binding to IL-22 measured in bioassay or assembly-induced IL-22 co-folding assay. We performed interactome analysis to disclose the mechanism underlying the poor secretion of IL-22BPi1 and identified GRP78, GRP94, GRP170, and calnexin as main interactors. Structure-function analysis revealed that, like IL-22BPi2, IL-22BPi1 binds to the substrate-binding domain of GRP78 as well as to the middle domain of GRP94. Ectopic expression of wild-type GRP78 enhanced, and ATPase-defective GRP94 mutant decreased, secretion of both IL-22BPi1 and IL-22BPi2, while neither of both affected IL-22BPi3 secretion. Thus, IL-22BPi1 and IL-22BPi2 are bona fide clients of the ER chaperones GRP78 and GRP94. However, only IL-22BPi1 activates an unfolded protein response (UPR) resulting in increased protein levels of GRP78 and GRP94. Cloning of the IL22RA2 alternatively spliced exon into an unrelated cytokine, IL-2, bestowed similar characteristics on the resulting protein. We also found that CD14++/CD16+ intermediate monocytes produced a higher level of IL22RA2 mRNA than classical and non-classical monocytes, but this difference disappeared in immature dendritic cells (moDC) derived thereof. Upon silencing of IL22RA2 expression in moDC, GRP78 levels were significantly reduced, suggesting that native IL22RA2 expression naturally contributes to upregulating GRP78 levels in these cells. The IL22RA2 alternatively spliced exon was reported to be recruited through a single mutation in the proto-splice site of a Long Terminal Repeat retrotransposon sequence in the ape lineage. Our work suggests that positive selection of IL-22BPi1 was not driven by IL-22 antagonism as in the case of IL-22BPi2 and IL-22BPi3, but by capacity for induction of an UPR response.

Human exonization through differential nucleosome occupancy

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, whereas the roles of nucleosome binding in shaping changes through evolution remain to be addressed. Here we performed a comparative study to clarify the roles of nucleosome occupancy in exon origination. By profiling a high-resolution, cross-species mononucleosome landscape for mammalian tissues, we found nucleosome occupancy profiles are conserved across tissues and species. Further, through a phylogenetic approach, we found that the feature of differential nucleosome occupancy appears prior to the origination of new exons and, presumably, facilitates the origin of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content, which suggests the function of nucleosome binding in exonization.

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, but the roles of nucleosome occupancy in shaping changes through evolution remain to be addressed. Here we present high-resolution nucleosome occupancy profiles for multiple tissues derived from human, macaque, tree shrew, mouse, and pig. Genome-wide comparison reveals conserved nucleosome occupancy profiles across both different species and tissue types. Notably, we found significantly higher levels of nucleosome occupancy in exons than in introns, a pattern correlated with the different exon–intron GC content. We then determined whether this biased occupancy may play roles in the origination of new exons through evolution, rather than being a downstream effect of exonization, through a comparative approach to sequentially trace the order of the exonization and biased nucleosome binding. By identifying recently evolved exons in human but not in macaque using matched RNA sequencing, we found that higher exonic nucleosome occupancy also existed in macaque regions orthologous to these exons. Presumably, such biased nucleosome occupancy facilitates the origination of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content. These data thus support a model that sites bound by nucleosomes are more likely to evolve into exons, which we term the “nucleosome-first” model.

A Multilayered Control of the Human Survival Motor Neuron Gene Expression by Alu Elements

Humans carry two nearly identical copies of Survival Motor Neuron gene: SMN1 and SMN2. Mutations or deletions of SMN1, which codes for SMN, cause spinal muscular atrophy (SMA), a leading genetic disease associated with infant mortality. Aberrant expression or localization of SMN has been also implicated in other pathological conditions, including male infertility, inclusion body myositis, amyotrophic lateral sclerosis and osteoarthritis. SMN2 fails to compensate for the loss of SMN1 due to skipping of exon 7, leading to the production of SMNΔ7, an unstable protein. In addition, SMNΔ7 is less functional due to the lack of a critical C-terminus of the full-length SMN, a multifunctional protein. Alu elements are specific to primates and are generally found within protein coding genes. About 41% of the human SMN gene including promoter region is occupied by more than 60 Alu-like sequences. Here we discuss how such an abundance of Alu-like sequences may contribute toward SMA pathogenesis. We describe the likely impact of Alu elements on expression of SMN. We have recently identified a novel exon 6B, created by exonization of an Alu-element located within SMN intron 6. Irrespective of the exon 7 inclusion or skipping, transcripts harboring exon 6B code for the same SMN6B protein that has altered C-terminus compared to the full-length SMN. We have demonstrated that SMN6B is more stable than SMNΔ7 and likely functions similarly to the full-length SMN. We discuss the possible mechanism(s) of regulation of SMN exon 6B splicing and potential consequences of the generation of exon 6B-containing transcripts.

The Ds1 Transposon Provides Messages That Yield Unique Profiles of Protein Isoforms and Acts Synergistically With Ds to Enrich Proteome Complexity via Exonization

In exonization events, Ds1 may provide donor and/or acceptor sites for splicing after inserting into genes and be incorporated into new transcripts with new exon(s). In this study, the protein variants of Ds1 exonization yielding additional functional profile(s) were studied. Unlike Ds exonization, which creates new profiles mostly by incorporating flanking intron sequences with the Ds message, Ds1 exonization additionally creates new profiles through the presence or absence of Ds1 messages. The number of unique functional profiles harboring Ds1 messages is 1.3-fold more than that of functional profiles without Ds1 messages. The highly similar 11 protein isoforms at a single insertion site also contribute to proteome complexity enrichment by exclusively creating new profiles. Particularly, Ds1 exonization produces 459 unique profiles, of which 129 cannot be built by Ds. We thus conclude that Ds and Ds1 are independent but synergistic in their capacity to enrich proteome complexity through exonization.

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.

Current Bacterial Gene Encoding Capsule Biosynthesis Protein CapI Contains Nucleotides Derived from Exonization

Since the proposition of introns-early hypothesis, although many studies have shown that most eukaryotic ancestors possessed intron-rich genomes, evidence of intron existence in genomes of ancestral bacteria has still been absent. While not a single intron has been found in all protein-coding genes of current bacteria, analyses on bacterial genes horizontally transferred into eukaryotes at ancient time may provide evidence of intron existence in bacterial ancestors. In this study, a bacterial gene encoding capsule biosynthesis protein CapI was found in the genome of sea anemone, Nematostella vectensis. This horizontally transferred gene contains a phase 1 intron of 40 base pairs. The nucleotides of this intron have high sequence identity with those encoding amino acids in current bacterial CapI gene, indicating that the intron and the amino acid-coding nucleotides are originated from the same ancestor sequence. Moreover, 5′-splice site of this intron is located in a GT-poor region associated with a closely following AG-rich region, suggesting that deletion mutation at 5′-splice site has been employed to remove this intron and the intron-like amino acid-coding nucleotides in current bacterial CapI gene are derived from exonization. These data suggest that bacterial CapI gene contained intron(s) at ancient time. This is the first report providing the result of sequence analysis to suggest possible existence of spliceosomal introns in ancestral bacterial genes. The methodology employed in this study may be used to identify more such evidence that would aid in settlement of the dispute between introns-early and introns-late theories.

Gain of a New Exon by a Lineage-Specific Alu Element-Integration Event in the BCS1L Gene during Primate Evolution

BCS1L gene encodes mitochondrial protein and is a member of conserved AAA protein family. This gene is involved in the incorporation of Rieske FeS and Qcr10p into complex III of respiratory chain. In our previous study, AluYRa2-derived alternative transcript in rhesus monkey genome was identified. However, this transcript has not been reported in human genome. In present study, we conducted evolutionary analysis of AluYRa2-exonized transcript with various primate genomic DNAs and cDNAs from humans, rhesus monkeys, and crab-eating monkeys. Remarkably, our results show that AluYRa2 element has only been integrated into genomes of Macaca species. This Macaca lineage-specific integration of AluYRa2 element led to exonization event in the first intron region of BCS1L gene by producing a conserved 3' splice site. Intriguingly, in rhesus and crab-eating monkeys, more diverse transcript variants by alternative splicing (AS) events, including exon skipping and different 5' splice sites from humans, were identified. Alignment of amino acid sequences revealed that AluYRa2-exonized transcript has short N-terminal peptides. Therefore, AS events play a major role in the generation of various transcripts and proteins during primate evolution. In particular, lineage-specific integration of Alu elements and species-specific Alu-derived exonization events could be important sources of gene diversification in primates.

Alport syndrome caused by a COL4A5 deletion and exonization of an adjacent AluY

Mutation-induced activation of splice sites in intronic repetitive sequences has contributed significantly to the evolution of exon–intron structure and genetic disease. Such events have been associated with mutations within transposable elements, most frequently in mutation hot-spots of Alus. Here, we report a case of Alu exonization resulting from a 367-nt genomic COL4A5 deletion that did not encompass any recognizable transposed element, leading to the Alport syndrome. The deletion brought to proximity the 5′ splice site of COL4A5 exon 33 and a cryptic 3′ splice site in an antisense AluY copy in intron 32. The fusion exon was depleted of purines and purine-rich splicing enhancers, but had low levels of intramolecular secondary structure, was flanked by short introns and had strong 5′ and Alu-derived 3′ splice sites, apparently compensating poor composition and context of the new exon. This case demonstrates that Alu splice sites can be activated by outlying deletions, highlighting Alu versatility in shaping the exon–intron organization and expanding the spectrum of mutational mechanisms that introduce repetitive sequences in mRNAs.

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.

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.

Intron sliding in tetraspanins

Specific questions about intron evolution are precisely addressed applying a phylogenomic approach to suitable gene families. With this approach we have recently reported that the appearance of most human tetraspanins occurred in the common ancestor of vertebrates and coincides in nearly all cases with the concomitant acquisition of new introns. We observed that indels at the ends of the DNA exonic sequences with no involvement of the corresponding intronic sequence, were the cause of two discordant intron positions between orthologous tetraspanins. Here, we discuss a putative intron sliding occurrence in which a new acquired intron junction (intron 1a) in the ancestor of chordates could have been shifted to new positions (introns 1b and 1c) during the expansion of the tetraspanin family in vertebrates. Such a mechanism could be responsible for generating some of the variation of function in this important family of membrane spanning proteins.

The birth of new exons: mechanisms and evolutionary consequences

A significant amount of literature was dedicated to hypotheses concerning the origin of ancient introns and exons, but accumulating evidence indicates that new exons are also constantly being added to evolving genomes. Several mechanisms contribute to the creation of novel exons in metazoan genomes, including whole gene and single exon duplications, but perhaps the most intriguing are events of exonization, where intronic sequences become exons de novo. Exonizations of intronic sequences, particularly those originating from repetitive elements, are now widely documented in many genomes including human, mouse, dog, and fish. Such de novo appearance of exons is very frequently associated with alternative splicing, with the new exon-containing variant typically being the rare one. This allows the new variant to be evolutionarily tested without compromising the original one, and provides an evolutionary strategy for generation of novel functions with minimum damage to the existing functional repertoire. This review discusses the molecular mechanisms leading to exonization, its extent in vertebrate genomes, and its evolutionary implications.