How did alternative splicing evolve?

Exonization of transposed elements: A challenge and opportunity for evolution

Protein-coding genes are composed of exons and introns flanked by untranslated regions. Before the mRNA of a gene can be translated into protein, the splicing machinery removes all the intronic regions and joins the protein-coding exons together. Exonization is a process, whereby genes acquire new exons from non-protein-coding, primarily intronic, DNA sequences. Genomic insertions or point mutations within DNA sequences often generate alternative splice sites, causing the splicing system to include new sequences as exons or to elongate existing exons. Because the alternative splice sites are not as efficient as the originals the new variants usually constitute a minor fraction of mature mRNAs. While the prevailing original splice variant maintains functionality, the additional sequence, free from selection pressure, evolves a new function or eventually vanishes. If the new splice variant is advantageous, selection might operate to optimize the new splice sites and consequently increase the proportion of the alternative splice variant. In some instances, the original splice variant is completely replaced by constitutive splicing of the new form. Because of the fortuitous presence of internal splice site-like structures within their sequences, portions of transposed elements frequently serve as modules of exonization. Their recruitment requires a long and versatile optimization process involving multiple changes over a time span of millions, even hundreds of millions, of years. Comparisons of corresponding genes and mRNAs in phylogenetically related species enables one to chronologically reconstruct such changes, from ancient ancestors to living species, in a stepwise manner. We will review this process using three different exemplary cases: (1) the evolution of a constitutively spliced mammalian-wide repeat (MIR), (2) the evolution of an alternative exon 1 from an alternative 5'-extended primary transcript containing an Alu element, and (3) a rare case of the stepwise exoniztion of an Alu element-derived sequence mediated by A-to-I RNA editing.

SUS1 introns are required for efficient mRNA nuclear export in yeast

Efficient coupling between mRNA synthesis and export is essential for gene expression. Sus1/ENY2, a component of the SAGA and TREX-2 complexes, is involved in both transcription and mRNA export. While most yeast genes lack introns, we previously reported that yeast SUS1 bears two. Here we show that this feature is evolutionarily conserved and critical for Sus1 function. We determine that while SUS1 splicing is inefficient, it responds to cellular conditions, and intronic mutations either promoting or blocking splicing lead to defects in mRNA export and cell growth. Consistent with this, we find that an intron-less SUS1 only partially rescues sus1Δ phenotypes. Remarkably, splicing of each SUS1 intron is also affected by the presence of the other and by SUS1 exonic sequences. Moreover, by following SUS1 RNA and protein levels we establish that nonsense-mediated decay (NMD) pathway and the splicing factor Mud2 both play a role in SUS1 expression. Our data (and those of the accompanying work by Hossain et al.) provide evidence of the involvement of splicing, translation, and decay in the regulation of early events in mRNP biogenesis; and imply the additional requirement for a balance in splicing isoforms from a single gene.

Spliceosome structure and function

Pre-mRNA splicing is catalyzed by the spliceosome, a multimegadalton ribonucleoprotein (RNP) complex comprised of five snRNPs and numerous proteins. Intricate RNA-RNA and RNP networks, which serve to align the reactive groups of the pre-mRNA for catalysis, are formed and repeatedly rearranged during spliceosome assembly and catalysis. Both the conformation and composition of the spliceosome are highly dynamic, affording the splicing machinery its accuracy and flexibility, and these remarkable dynamics are largely conserved between yeast and metazoans. Because of its dynamic and complex nature, obtaining structural information about the spliceosome represents a major challenge. Electron microscopy has revealed the general morphology of several spliceosomal complexes and their snRNP subunits, and also the spatial arrangement of some of their components. X-ray and NMR studies have provided high resolution structure information about spliceosomal proteins alone or complexed with one or more binding partners. The extensive interplay of RNA and proteins in aligning the pre-mRNA's reactive groups, and the presence of both RNA and protein at the core of the splicing machinery, suggest that the spliceosome is an RNP enzyme. However, elucidation of the precise nature of the spliceosome's active site, awaits the generation of a high-resolution structure of its RNP core.

Compensatory signals associated with the activation of human GC 5' splice sites

GC 5′ splice sites (5′ss) are present in ∼1% of human introns, but factors promoting their efficient selection are poorly understood. Here, we describe a case of X-linked agammaglobulinemia resulting from a GC 5′ss activated by a mutation in BTK intron 3. This GC 5′ss was intrinsically weak, yet it was selected in >90% primary transcripts in the presence of a strong and intact natural GT counterpart. We show that efficient selection of this GC 5′ss required a high density of GAA/CAA-containing splicing enhancers in the exonized segment and was promoted by SR proteins 9G8, Tra2β and SC35. The GC 5′ss was efficiently inhibited by splice-switching oligonucleotides targeting either the GC 5′ss itself or the enhancer. Comprehensive analysis of natural GC-AG introns and previously reported pathogenic GC 5′ss showed that their efficient activation was facilitated by higher densities of splicing enhancers and lower densities of silencers than their GT 5′ss equivalents. Removal of the GC-AG introns was promoted to a minor extent by the splice-site strength of adjacent exons and inhibited by flanking Alu repeats, with the first downstream Alus located on average at a longer distance from the GC 5′ss than other transposable elements. These results provide new insights into the splicing code that governs selection of noncanonical splice sites.

Blind sterile 2 (bs2), a hypomorphic mutation in Agps, results in cataracts and male sterility in mice

Blind sterile 2 (bs2) is a spontaneous autosomal recessive mouse mutation exhibiting cataracts and male sterility. Detailed clinical and histological evaluation revealed that bs2 mice have cataracts resulting from severely disrupted lens fiber cells. Analysis of bs2 testes revealed the absence of mature sperm and the presence of large multinucleate cells within the lumens of seminiferous tubules. Linkage analysis mapped the bs2 locus to mouse chromosome 2, approximately 45cM distal from the centromere. Fine mapping established a 3.1Mb bs2 critical region containing 19 candidate genes. Sequence analysis of alkylglycerone-phosphate synthase (Agps), a gene within the bs2 critical region, revealed a G to A substitution at the +5 position of intron 14. This mutation results in two abundantly expressed aberrantly spliced Agps transcripts: Agps(∆exon14) lacking exon 14 or Agps(exon∆13-14) lacking both exons 13 and 14 as well as full-length Agps transcript. Agps is a peroxisomal enzyme which catalyzes the formation of the ether bond during the synthesis of ether lipids. Both aberrantly spliced Agps(∆exon14) and Agps(exon∆13-14) transcripts led to a frame shift, premature stop and putative proteins lacking the enzymatic FAD domain. We present evidence that bs2 mice have significantly decreased levels of ether lipids. Human mutations in Agps result in rhizomelic chondrodysplasia punctata type 3 (RCDP3), a disease for which bs2 is the only genetic model. Thus, bs2 is a hypomorphic mutation in Agps, and represents a useful model for investigation of the tissue specificity of ether lipid requirements which will be particularly valuable for elucidating the mechanism of disease phenotypes resulting from ether lipid depletion.

Identification of constitutively active interleukin 33 (IL-33) splice variant

IL-33/IL-1F11 is a new member of the IL-1 family ligand and provokes T helper-type immune responses. IL-33 is the ligand of ST2 and IL-1 receptor accessory protein (IL-1RAcP) that triggers nuclear factor-κ light chain enhancer of activated B cells (NF-κB) and MAPK signaling. We discovered a novel short splice variant of IL-33 that was termed spIL-33. The new spIL-33 lacks exon 3 containing a proposed caspase-1 cleavage site. We isolated spIL-33 cDNA from the Huh7 human hepatocarcinoma cell line and expressed the recombinant spIL-33 protein in Escherichia coli. The recombinant spIL-33 and pro-IL-33 were not cleaved by caspase-1, unlike IL-18 (IL-1F4). The recombinant spIL-33 was constitutively active, and spIL-33-induced inflammatory cytokine production was caspase-1-independent in HMC-1 and Raw 264.7 cells. The recombinant spIL-33 induced the phosphorylation of IL-1 receptor-associated kinase (IRAK1), NF-κB, p38 MAPK, p44/42 MAPK, and JNK in a time- and dose-dependent manner. Anti-ST2 monoclonal antibody specifically blocked the spIL-33-induced cytokine production. In this study, we identified and characterized a new IL-33 splice variant, which was a constitutively active IL-33 isoform. The existence of constitutively active spIL-33 suggests that the biological activity of IL-33 could be triggered by diverse stimulations during immune responses. Further investigation of the spIL-33 expression pattern may contribute to understanding the involvement of IL-33 in inflammatory disorders.

Gene therapeutic approach using mutation-adapted U1 snRNA to correct a RPGR splice defect in patient-derived cells

Retinitis pigmentosa (RP) is a disease that primarily affects the peripheral retina and ultimately causes visual impairment. X-chromosomal forms of RP are frequently caused by mutations in the retinitis pigmentosa GTPase regulator (RPGR) gene. We show that the novel splice donor site (SDS) mutation c.1245+3A>T in intron 10 of RPGR cosegregates with RP in a five-generation Caucasian family. The mutation causes in-frame skipping of exon 10 from RPGR transcripts in patient-derived primary fibroblasts. To correct the splice defect, we developed a gene therapeutic approach using mutation-adapted U1 small nuclear RNA (U1). U1 is required for SDS recognition of pre-mRNAs and initiates the splice process. The mutation described herein interferes with the recognition of the SDS by U1. To overcome the deleterious effects of the mutation, we generated four U1 isoforms with increasing complementarity to the SDS. Lentiviral particles were used to transduce patient-derived fibroblasts with these U1 variants. Full complementarity of U1 corrects the splice defect partially and increases recognition of the mutant SDS. The therapeutic effect is U1-concentration dependent as we show for endogenously expressed RPGR transcripts in patient-derived cells. U1-based gene therapeutic approaches constitute promising technologies to treat SDS mutations in inherited diseases including X-linked RP.

Comparative analysis of information contents relevant to recognition of introns in many species

Background

The basic process of RNA splicing is conserved among eukaryotic species. Three signals (5' and 3' splice sites and branch site) are commonly used to directly conduct splicing, while other features are also related to the recognition of an intron. Although there is experimental evidence pointing to the significant species specificities in the features of intron recognition, a quantitative evaluation of the divergence of these features among a wide variety of eukaryotes has yet to be conducted.

Results

To better understand the splicing process from the viewpoints of evolution and information theory, we collected introns from 61 diverse species of eukaryotes and analyzed the properties of the nucleotide sequences relevant to splicing. We found that trees individually constructed from the five features (the three signals, intron length, and nucleotide composition within an intron) roughly reflect the phylogenetic relationships among the species but sometimes extensively deviate from the species classification. The degree of topological deviation of each feature tree from the reference trees indicates the lowest discordance for the 5' splicing signal, followed by that for the 3' splicing signal, and a considerably greater discordance for the other three features. We also estimated the relative contributions of the five features to short intron recognition in each species. Again, moderate correlation was observed between the similarities in pattern of short intron recognition and the genealogical relationships among the species. When mammalian introns were categorized into three subtypes according to their terminal dinucleotide sequences, each subtype segregated into a nearly monophyletic group, regardless of the host species, with respect to the 5' and 3' splicing signals. It was also found that GC-AG introns are extraordinarily abundant in some species with high genomic G + C contents, and that the U12-type spliceosome might make a greater contribution than currently estimated in most species.

Conclusions

Overall, the present study indicates that both splicing signals themselves and their relative contributions to short intron recognition are rather susceptible to evolutionary changes, while some poorly characterized properties seem to be preserved within the mammalian intron subtypes. Our findings may afford additional clues to understanding of evolution of splicing mechanisms.

Analysis of mutations that influence pre-mRNA splicing

A rapidly increasing number of human diseases are now recognized as being caused by the selection of wrong splice sites. In most cases, these changes in alternative splice site selection are due to single nucleotide exchanges in splicing regulatory elements. This chapter describes the use of bioinformatics tools to predict the influence of a mutation on alternative pre-mRNA splicing and the experimental testing of these predictions. The bioinformatic analysis determines the influence of a mutation on splicing enhancers and silencers, splice sites and RNA secondary structures. This approach generates hypotheses that are tested using splicing reporter constructs, which are then analyzed in transfection assays. We describe a recombination-based system that allows for the generation of splicing reporter constructs in the first week and their subsequent analysis in the second week.

Genome-wide survey of alternative splicing in the grass Brachypodium distachyon: a emerging model biosystem for plant functional genomics

A draft sequence of the genome of Brachypodium distachyon, the emerging grass model, was recently released. This represents a unique opportunity to determine its functional diversity compared to the genomes of other model species. Using homology mapping of assembled expressed sequence tags with chromosome scale pseudomolecules, we identified 128 alternative splicing events in B. distachyon. Our study identified that retention of introns is the major type of alternative splicing events (53%) in this plant and highlights the prevalence of splicing site recognition for definition of introns in plants. We have analyzed the compositional profiles of exon-intron junctions by base-pairing nucleotides with U1 snRNA which serves as a model for describing the possibility of sequence conservation. The alternative splicing isoforms identified in this study are novel and represent one of the potentially biologically significant means by which B. distachyon controls the function of its genes. Our observations serve as a basis to understand alternative splicing events of cereal crops with more complex genomes, like wheat or barley.

Anti-apoptosis and cell survival: a review

Type I programmed cell death (PCD) or apoptosis is critical for cellular self-destruction for a variety of processes such as development or the prevention of oncogenic transformation. Alternative forms, including type II (autophagy) and type III (necrotic) represent the other major types of PCD that also serve to trigger cell death. PCD must be tightly controlled since disregulated cell death is involved in the development of a large number of different pathologies. To counter the multitude of processes that are capable of triggering death, cells have devised a large number of cellular processes that serve to prevent inappropriate or premature PCD. These cell survival strategies involve a myriad of coordinated and systematic physiological and genetic changes that serve to ward off death. Here we will discuss the different strategies that are used to prevent cell death and focus on illustrating that although anti-apoptosis and cellular survival serve to counteract PCD, they are nevertheless mechanistically distinct from the processes that regulate cell death.

Functional implications of the emergence of alternative splicing in hnRNP A/B transcripts

The heterogeneous nuclear ribonucleoproteins (hnRNPs) A/B are a family of RNA-binding proteins that participate in various aspects of nucleic acid metabolism, including mRNA trafficking, telomere maintenance, and splicing. They are both regulators and targets of alternative splicing, and the patterns of alternative splicing of their transcripts have diverged between paralogs and between orthologs in different species. Surprisingly, the extent of this splicing variation and its implications for post-transcriptional regulation have remained largely unexplored. Here, we conducted a detailed analysis of hnRNP A/B sequences and expression patterns across six vertebrates. Alternative exons emerged via the introduction of new splice sites, changes in the strengths of existing splice sites, and the accumulation of auxiliary splicing regulatory motifs. Observed isoform expression patterns could be attributed to the frequency and strength of cis-elements. We found a trend toward increased splicing variation in mammals and identified novel alternatively spliced isoforms in human and chicken. Pulldown and translational assays demonstrated that the inclusion of alternative exons altered the affinity of hnRNP A/B proteins for their cognate nucleic acids and modified protein expression levels. As the hnRNPs A/B regulate several key steps in mRNA processing, the involvement of diverse hnRNP isoforms in multiple cellular contexts and species implies concomitant differences in the transcriptional output of these systems. We conclude that the emergence of alternative splicing in the hnRNPs A/B has contributed to the diversification of their roles in the regulation of alternative splicing and has thus added an unexpected layer of regulatory complexity to transcription in vertebrates.

Alu-derived old world monkeys exonization event and experimental validation of the LEPR gene

The leptin receptor (LEPR) is a crucial regulatory protein that interacts with Leptin. In our analysis of LEPR, novel AluJb-derived alternative transcripts were identified in the genome of the rhesus monkey. In order to investigate the occurrence of AluJb-derived alternative transcripts and the mechanism underlying exonization events, we conducted analyses using a number of primate genomic DNAs and adipose RNAs of tissue and primary cells derived from the crab-eating monkey. Our results demonstrate that the AluJb element has been integrated into our common ancestor genome prior to the divergence of simians and prosimians. The lineage-specific exonization event of the LEPR gene in chimpanzees, orangutans, and Old World monkeys appear to have been accomplished via transition mutations of the 5' splicing site (second position of C to T). However, in New World monkeys and prosimians, the AluJb-related LEPR transcript should be silenced by the additional transversion mutation (fourth position of T to G). The AluJb-related transcript of human LEPR should also be silenced by a mutation of the 5' splicing site (first position of G to A) and the insertion of one nucleotide sequence (minus fourth position of A). Our data suggests that lineage-specific exonization events should be determined by the combination event of the formation of splicing sites and protection against site-specific mutation pressures. These evolutionary mechanisms could be major sources for primate diversification.

A humanized Smn gene containing the SMN2 nucleotide alteration in exon 7 mimics SMN2 splicing and the SMA disease phenotype

Proximal spinal muscular atrophy (SMA) is a neurodegenerative disease caused by low levels of the survival motor neuron (SMN) protein. In humans, SMN1 and SMN2 encode the SMN protein. In SMA patients, the SMN1 gene is lost and the remaining SMN2 gene only partially compensates. Mediated by a C>T nucleotide transition in SMN2, the inefficient recognition of exon 7 by the splicing machinery results in low levels of SMN. Because the SMN2 gene is capable of expressing SMN protein, correction of SMN2 splicing is an attractive therapeutic option. Although current mouse models of SMA characterized by Smn knock-out alleles in combination with SMN2 transgenes adequately model the disease phenotype, their complex genetics and short lifespan have hindered the development and testing of therapies aimed at SMN2 splicing correction. Here we show that the mouse and human minigenes are regulated similarly by conserved elements within in exon 7 and its downstream intron. Importantly, the C>T mutation is sufficient to induce exon 7 skipping in the mouse minigene as in the human SMN2. When the mouse Smn gene was humanized to carry the C>T mutation, keeping it under the control of the endogenous promoter, and in the natural genomic context, the resulting mice exhibit exon 7 skipping and mild adult onset SMA characterized by muscle weakness, decreased activity and an alteration of the muscle fibers size. This Smn C>T mouse represents a new model for an adult onset form of SMA (type III/IV) also know as the Kugelberg-Welander disease.

The nature of protein domain evolution: shaping the interaction network

The proteomes that make up the collection of proteins in contemporary organisms evolved through recombination and duplication of a limited set of domains. These protein domains are essentially the main components of globular proteins and are the most principal level at which protein function and protein interactions can be understood. An important aspect of domain evolution is their atomic structure and biochemical function, which are both specified by the information in the amino acid sequence. Changes in this information may bring about new folds, functions and protein architectures. With the present and still increasing wealth of sequences and annotation data brought about by genomics, new evolutionary relationships are constantly being revealed, unknown structures modeled and phylogenies inferred. Such investigations not only help predict the function of newly discovered proteins, but also assist in mapping unforeseen pathways of evolution and reveal crucial, co-evolving inter- and intra-molecular interactions. In turn this will help us describe how protein domains shaped cellular interaction networks and the dynamics with which they are regulated in the cell. Additionally, these studies can be used for the design of new and optimized protein domains for therapy. In this review, we aim to describe the basic concepts of protein domain evolution and illustrate recent developments in molecular evolution that have provided valuable new insights in the field of comparative genomics and protein interaction networks.

Identification and characterization of DSPIa, a novel isoform of human desmoplakin

Desmoplakin is a ubiquitous component of desmosomes and desmosome-like structures, such as the cardiomyocyte area composita. Two major isoforms, desmoplakin I (DSPI) and desmoplakin II (DSPII) are encoded by alternative mRNA transcripts differentially spliced from the same gene. The resulting proteins are identical in amino acid sequence with the exception that DSPII contains only one third of the central alpha-helical rod domain present in DSPI. Here we describe a novel minor isoform of desmoplakin that is also produced by alternative splicing of the desmoplakin gene and that we name desmoplakin Ia (DSPIa). DSPIa is an alternatively spliced DSPI mRNA with a unique splice donor site that is 90% homologous to and downstream of the DSPII specific donor. The resulting DSPIa mRNA is in-frame and encodes a protein that has a central alpha-helical rod domain of intermediate size and that is 156 amino acids larger than DSPII and 443 amino acids smaller than DSPI. We demonstrate, through recombinant expression and short interfering RNA knockdown, that the DSPIa protein is readily detectable, albeit at substantially lower levels than the dominant isoforms, DSPI and DSPII. DSPIa mRNA has a similar tissue distribution to that of DSPI and of DSPII.

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.

Understanding protein non-folding

This review describes the family of intrinsically disordered proteins, members of which fail to form rigid 3-D structures under physiological conditions, either along their entire lengths or only in localized regions. Instead, these intriguing proteins/regions exist as dynamic ensembles within which atom positions and backbone Ramachandran angles exhibit extreme temporal fluctuations without specific equilibrium values. Many of these intrinsically disordered proteins are known to carry out important biological functions which, in fact, depend on the absence of a specific 3-D structure. The existence of such proteins does not fit the prevailing structure-function paradigm, which states that a unique 3-D structure is a prerequisite to function. Thus, the protein structure-function paradigm has to be expanded to include intrinsically disordered proteins and alternative relationships among protein sequence, structure, and function. This shift in the paradigm represents a major breakthrough for biochemistry, biophysics and molecular biology, as it opens new levels of understanding with regard to the complex life of proteins. This review will try to answer the following questions: how were intrinsically disordered proteins discovered? Why don't these proteins fold? What is so special about intrinsic disorder? What are the functional advantages of disordered proteins/regions? What is the functional repertoire of these proteins? What are the relationships between intrinsically disordered proteins and human diseases?

An autoinhibitory effect of the homothorax domain of Meis2

Myeloid ecotropic insertion site (Meis)2 is a homeodomain protein containing a conserved homothorax (Hth) domain that is present in all Meis and Prep family proteins and in the Drosophila Hth protein. The Hth domain mediates interaction with Pbx homeodomain proteins, allowing for efficient DNA binding. Here we show that, like Meis1, Meis2 has a strong C-terminal transcriptional activation domain, which is required for full activation of transcription by homeodomain protein complexes composed of Meis2 and Pbx1. We also show that the activity of the activation domain is inhibited by the Hth domain, and that this autoinhibition can be partially relieved by the interaction of Pbx1 with the Hth domain of Meis2. Targeting of the Hth domain to DNA suggests that it is not a portable trans-acting repression domain. However, the Hth domain can inhibit a linked activation domain, and this inhibition is not limited to the Meis2 activation domain. Database searching reveals that the Meis3.2 splice variant, which is found in several vertebrate species, disrupts the Hth domain by removing 17 codons from the 5'-end of exon 6. We show that the equivalent deletion in Meis2 derepresses the C-terminal activation domain and weakens interaction with Pbx1. This work suggests that the transcriptional activity of all members of the Meis/Prep Hth protein family is subject to autoinhibition by their Hth domains, and that the Meis3.2 splice variant encodes a protein that bypasses this autoinhibitory effect.

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.

Evolution of alternative splicing in primate brain transcriptomes

Alternative splicing is a predominant form of gene regulation in higher eukaryotes. The evolution of alternative splicing provides an important mechanism for the acquisition of novel gene functions. In this work, we carried out a genome-wide phylogenetic survey of lineage-specific splicing patterns in the primate brain, via high-density exon junction array profiling of brain transcriptomes of humans, chimpanzees and rhesus macaques. We identified 509 genes showing splicing differences among these species. RT-PCR analysis of 40 exons confirmed the predicted splicing evolution of 33 exons. Of these 33 exons, outgroup analysis using rhesus macaques confirmed 13 exons with human-specific increase or decrease in transcript inclusion levels after humans diverged from chimpanzees. Some of the human-specific brain splicing patterns disrupt domains critical for protein-protein interactions, and some modulate translational efficiency of their host genes. Strikingly, for exons showing splicing differences across species, we observed a significant increase in the rate of silent substitutions within exons, coupled with accelerated sequence divergence in flanking introns. This indicates that evolution of cis-regulatory signals is a major contributor to the emergence of human-specific splicing patterns. In one gene (MAGOH), using minigene reporter assays, we demonstrated that the combination of two human-specific cis-sequence changes created its human-specific splicing pattern. Together, our data reveal widespread human-specific changes of alternative splicing in the brain and suggest an important role of splicing in the evolution of neuronal gene regulation and functions.

Exon definition complexes contain the tri-snRNP and can be directly converted into B-like precatalytic splicing complexes

The first step in splicing of pre-mRNAs with long introns is exon definition, where U1 and U2 snRNPs bind at opposite ends of an exon. After exon definition, these snRNPs must form a complex across the upstream intron to allow splicing catalysis. Exon definition and conversion of cross-exon to cross-intron spliceosomal complexes are poorly understood. Here we demonstrate that, in addition to U1 and U2 snRNPs, cross-exon complexes contain U4, U5, and U6 (which form the tri-snRNP). Tri-snRNP docking involves the formation of U2/U6 helix II. This interaction is stabilized by a 5' splice site (SS)-containing oligonucleotide, which can bind the tri-snRNP and convert the cross-exon complex into a cross-intron, B-like complex. Our data suggest that the switch from cross-exon to cross-intron complexes can occur directly when an exon-bound tri-snRNP interacts with an upstream 5'SS, without prior formation of a cross-intron A complex, revealing an alternative spliceosome assembly pathway.

Genome-wide analysis of alternative splicing evolution among Mus subspecies

Alternative splicing, the combination of different exons to produce a variety of transcripts from a single gene, contributes enormously to transcriptome diversity in mammals, and the majority of genes encode alternatively spliced products. Previous research comparing mouse, rat and human has shown that a significant proportion of splice forms are not conserved across species, suggesting that alternative transcripts are an important source of evolutionary novelty. Here, we studied the evolution of alternative splicing in the early stages of species divergence in the house mouse. We sequenced the testis transcriptomes of three Mus musculus subspecies and Mus spretus using Illumina technology. On the basis of a genome-wide analysis of read coverage differences among subspecies, we identified several hundred candidate alternatively spliced regions. We conservatively estimate that 6.5% of testis-expressed genes show alternative splice differences between at least one pair of M. musculus subspecies, a proportion slightly higher than the proportion of genes differentially expressed among subspecies. These results suggest that differences in both the structure and abundance of transcripts contribute to early transcriptome divergence.

Molecular study of ovine glucose 6-phosphate dehydrogenase gene expression in respect to different energy intake

Glucose 6-phosphate dehydrogenase (G6PD) plays an important role in a ruminant's metabolism catalyzing the first committed reaction in the pentose phosphate pathway as it provides necessary compounds of NADPH for the synthesis of fatty acids. The cloning of ovine (Ovis aries) G6PD gene revealed the presence of two cDNA transcripts (oG6PD(A) and oG6PD(B)), with oG6PD(B) being a product of alternative splicing and with no similarity to any other previously reported G6PD transcript. Here, we attempt to study the effect of energy balance in ovine G6PD transcript expression, trying simultaneously to find out any potential physiological role of the oG6PD(B) transcript. Changes of energy balance that lead to synergistic changes in the expression of both transcripts, but in opposite directions and not in a proportional way. Negative energy balance favours the presence of the oG6PD(B) transcript leading to a significant increase of its expression, compared to oG6PD(A) expression (P<0.05). In contrast, positive energy balance leads to a significant increase of oG6PD(A) compared to oG6PD(B) expression (P<0.05). In either condition oG6PD(B) expression is unchanged. Regression analysis showed that there is an energy balance threshold where the expression of both transcripts shows no change.

Chromatin density and splicing destiny: on the cross-talk between chromatin structure and splicing

How are short exonic sequences recognized within the vast intronic oceans in which they reside? Despite decades of research, this remains one of the most fundamental, yet enigmatic, questions in the field of pre-mRNA splicing research. For many years, studies aiming to shed light on this process were focused at the RNA level, characterizing the manner by which splicing factors and auxiliary proteins interact with splicing signals, thereby enabling, facilitating and regulating splicing. However, we increasingly understand that splicing is not an isolated process; rather it occurs co-transcriptionally and is presumably also regulated by transcription-related processes. In fact, studies by our group and others over the past year suggest that DNA structure in terms of nucleosome positioning and specific histone modifications, which have a well established role in transcription, may also have a role in splicing. In this review we discuss evidence for the coupling between transcription and splicing, focusing on recent findings suggesting a link between chromatin structure and splicing, and highlighting challenges this emerging field is facing.

Survey of the transcriptome of Aspergillus oryzae via massively parallel mRNA sequencing

Aspergillus oryzae, an important filamentous fungus used in food fermentation and the enzyme industry, has been shown through genome sequencing and various other tools to have prominent features in its genomic composition. However, the functional complexity of the A. oryzae transcriptome has not yet been fully elucidated. Here, we applied direct high-throughput paired-end RNA-sequencing (RNA-Seq) to the transcriptome of A. oryzae under four different culture conditions. With the high resolution and sensitivity afforded by RNA-Seq, we were able to identify a substantial number of novel transcripts, new exons, untranslated regions, alternative upstream initiation codons and upstream open reading frames, which provide remarkable insight into the A. oryzae transcriptome. We were also able to assess the alternative mRNA isoforms in A. oryzae and found a large number of genes undergoing alternative splicing. Many genes and pathways that might be involved in higher levels of protein production in solid-state culture than in liquid culture were identified by comparing gene expression levels between different cultures. Our analysis indicated that the transcriptome of A. oryzae is much more complex than previously anticipated, and these results may provide a blueprint for further study of the A. oryzae transcriptome.

BAK and BAX deletion using zinc-finger nucleases yields apoptosis-resistant CHO cells

Anoxic and metabolic stresses in large-scale cell culture during biopharmaceutical production can induce apoptosis. Strategies designed to ameliorate the problem of apoptosis in cell culture have focused on mRNA knockdown of pro-apoptotic proteins and over-expression of anti-apoptotic ones. Apoptosis in cell culture involves mitochondrial permeabilization by the pro-apoptotic Bak and Bax proteins; activity of either protein is sufficient to permit apoptosis. We demonstrate here the complete and permanent elimination of both the Bak and Bax proteins in combination in Chinese hamster ovary (CHO) cells using zinc-finger nuclease-mediated gene disruption. Zinc-finger nuclease cleavage of BAX and BAK followed by inaccurate DNA repair resulted in knockout of both genes. Cells lacking Bax and Bak grow normally but fail to activate caspases in response to apoptotic stimuli. When grown using scale-down systems under conditions that mimic growth in large-scale bioreactors they are significantly more resistant to apoptosis induced by starvation, staurosporine, and sodium butyrate. When grown under starvation conditions, BAX- and BAK-deleted cells produce two- to fivefold more IgG than wild-type CHO cells. Under normal growth conditions in suspension culture in shake flasks, double-knockout cultures achieve equal or higher cell densities than unmodified wild-type cultures and reach viable cell densities relevant for large-scale industrial protein production.

Overlapping splicing regulatory motifs--combinatorial effects on splicing

Regulation of splicing in eukaryotes occurs through the coordinated action of multiple splicing factors. Exons and introns contain numerous putative binding sites for splicing regulatory proteins. Regulation of splicing is presumably achieved by the combinatorial output of the binding of splicing factors to the corresponding binding sites. Although putative regulatory sites often overlap, no extensive study has examined whether overlapping regulatory sequences provide yet another dimension to splicing regulation. Here we analyzed experimentally-identified splicing regulatory sequences using a computational method based on the natural distribution of nucleotides and splicing regulatory sequences. We uncovered positive and negative interplay between overlapping regulatory sequences. Examination of these overlapping motifs revealed a unique spatial distribution, especially near splice donor sites of exons with weak splice donor sites. The positively selected overlapping splicing regulatory motifs were highly conserved among different species, implying functionality. Overall, these results suggest that overlap of two splicing regulatory binding sites is an evolutionary conserved widespread mechanism of splicing regulation. Finally, over-abundant motif overlaps were experimentally tested in a reporting minigene revealing that overlaps may facilitate a mode of splicing that did not occur in the presence of only one of the two regulatory sequences that comprise it.

Synonymous polymorphisms at splicing regulatory sites are associated with CpGs in neurodegenerative disease-related genes

Neuronal plasticity is associated with alternative splicing and epigenetic modulation. Recent evidence reveals the association of cytosine methylation with alternative splicing and splicing regulatory mechanisms. Single nucleotide polymorphisms (SNPs) are generally less frequent in conserved coding regions and probably in splice sites, compared to non-coding regions. CpG polymorphisms in coding regions and splice sites and their association with splicing regulatory elements have not been investigated till presently. We currently analyzed the CpG variability in 28 genes (361 constitutive and 105 alternative exons and the corresponding splice sites) associated with neurodegenerative diseases (ND). CpG polymorphisms in the splice sites of these genes are particularly frequent when compared to those at AG sequences. Moreover, in both constitutive and alternative exons, polymorphisms in CpGs are more frequent than in AG, GT sequences. On the contrary, in the polypyrimidine acceptor sequence C/T conservation is prominent indicating that in this locus the sequence of cytosines and thymines is preserved. Bioinformatic analysis of the splicing-associated regulatory elements in these exons and splice sites reveals that 18 out of a total of 39 SNPs which could strongly affect splicing (>1.5 score difference) contain CpG sequences. Cytosines are considerably more frequent and variable than expected at the position preceding the GT splice donors, while sites of epigenetic modification are absent from acceptors. The high CpG frequency in polymorphic splicing-associated sites implicates the involvement of epigenetic mechanisms in splicing selection decisions regulated by these sites, and indicates the complexity of genetic studies involving these, tentatively critical, polymorphisms in ND.

Genome-wide evidence for selection acting on single amino acid repeats

Low complexity and homopolymer sequences within coding regions are known to evolve rapidly. While their expansion may be deleterious, there is increasing evidence for a functional role associated with these amino acid sequences. Homopolymer sequences are thought to evolve mostly through replication slippage and, therefore, they may be expected to be longer in regions with relaxed selective constraint. Within the coding sequences of eukaryotes, alternatively spliced exons are known to evolve under relaxed constraints in comparison to those exons that are constitutively spliced because they are not included in all of the mature mRNA of a gene. This relaxed exposure to selection leads to faster rates of evolution for alternatively spliced exons in comparison to constitutively spliced exons. Here, we have tested the effect of splicing on the structure (composition, length) of homopolymer sequences in relation to the splicing pattern in which they are found. We observed a significant relationship between alternative splicing and homopolymer sequences with alternatively spliced genes being enriched in number and length of homopolymer sequences. We also observed lower codon diversity and longer homocodons, suggesting a balance between slippage and point mutations linked to the constraints imposed by selection.

RNA processing and its regulation: global insights into biological networks

In recent years views of eukaryotic gene expression have been transformed by the finding that enormous diversity can be generated at the RNA level. Advances in technologies for characterizing RNA populations are revealing increasingly complete descriptions of RNA regulation and complexity; for example, through alternative splicing, alternative polyadenylation and RNA editing. New biochemical strategies to map protein-RNA interactions in vivo are yielding transcriptome-wide insights into mechanisms of RNA processing. These advances, combined with bioinformatics and genetic validation, are leading to the generation of functional RNA maps that reveal the rules underlying RNA regulation and networks of biologically coherent transcripts. Together these are providing new insights into molecular cell biology and disease.

Predicting functional alternative splicing by measuring RNA selection pressure from multigenome alignments

High-throughput methods such as EST sequencing, microarrays and deep sequencing have identified large numbers of alternative splicing (AS) events, but studies have shown that only a subset of these may be functional. Here we report a sensitive bioinformatics approach that identifies exons with evidence of a strong RNA selection pressure ratio (RSPR)--i.e., evolutionary selection against mutations that change only the mRNA sequence while leaving the protein sequence unchanged--measured across an entire evolutionary family, which greatly amplifies its predictive power. Using the UCSC 28 vertebrate genome alignment, this approach correctly predicted half to three-quarters of AS exons that are known binding targets of the NOVA splicing regulatory factor, and predicted 345 strongly selected alternative splicing events in human, and 262 in mouse. These predictions were strongly validated by several experimental criteria of functional AS such as independent detection of the same AS event in other species, reading frame-preservation, and experimental evidence of tissue-specific regulation: 75% (15/20) of a sample of high-RSPR exons displayed tissue specific regulation in a panel of ten tissues, vs. only 20% (4/20) among a sample of low-RSPR exons. These data suggest that RSPR can identify exons with functionally important splicing regulation, and provides biologists with a dataset of over 600 such exons. We present several case studies, including both well-studied examples (GRIN1) and novel examples (EXOC7). These data also show that RSPR strongly outperforms other approaches such as standard sequence conservation (which fails to distinguish amino acid selection pressure from RNA selection pressure), or pairwise genome comparison (which lacks adequate statistical power for predicting individual exons).

Large-scale discovery of insertion hotspots and preferential integration sites of human transposed elements

Throughout evolution, eukaryotic genomes have been invaded by transposable elements (TEs). Little is known about the factors leading to genomic proliferation of TEs, their preferred integration sites and the molecular mechanisms underlying their insertion. We analyzed hundreds of thousands nested TEs in the human genome, i.e. insertions of TEs into existing ones. We first discovered that most TEs insert within specific ‘hotspots’ along the targeted TE. In particular, retrotransposed Alu elements contain a non-canonical single nucleotide hotspot for insertion of other Alu sequences. We next devised a method for identification of integration sequence motifs of inserted TEs that are conserved within the targeted TEs. This method revealed novel sequences motifs characterizing insertions of various important TE families: Alu, hAT, ERV1 and MaLR. Finally, we performed a global assessment to determine the extent to which young TEs tend to nest within older transposed elements and identified a 4-fold higher tendency of TEs to insert into existing TEs than to insert within non-TE intergenic regions. Our analysis demonstrates that TEs are highly biased to insert within certain TEs, in specific orientations and within specific targeted TE positions. TE nesting events also reveal new characteristics of the molecular mechanisms underlying transposition.

Entropic contributions to the splicing process

It has been recently argued that depletion attraction may play an important role in different aspects of cellular organization, ranging from the organization of transcriptional activity in transcription factories to the formation of nuclear bodies. In this paper, we suggest a new application of these ideas in the context of the splicing process, a crucial step of messenger RNA maturation in eukaryotes. We shall show that entropy effects and the resulting depletion attraction may explain the relevance of the aspecific intron length variable in the choice of splice-site recognition modality. On top of that, some qualitative features of the genome architecture of higher eukaryotes can find evolutionary realistic motivation in the light of our model.

The pivotal roles of TIA proteins in 5' splice-site selection of alu exons and across evolution

More than 5% of alternatively spliced internal exons in the human genome are derived from Alu elements in a process termed exonization. Alus are comprised of two homologous arms separated by an internal polypyrimidine tract (PPT). In most exonizations, splice sites are selected from within the same arm. We hypothesized that the internal PPT may prevent selection of a splice site further downstream. Here, we demonstrate that this PPT enhanced the selection of an upstream 5′ splice site (5′ss), even in the presence of a stronger 5′ss downstream. Deletion of this PPT shifted selection to the stronger downstream 5′ss. This enhancing effect depended on the strength of the downstream 5′ss, on the efficiency of base-pairing to U1 snRNA, and on the length of the PPT. This effect of the PPT was mediated by the binding of TIA proteins and was dependent on the distance between the PPT and the upstream 5′ss. A wide-scale evolutionary analysis of introns across 22 eukaryotes revealed an enrichment in PPTs within ∼20 nt downstream of the 5′ss. For most metazoans, the strength of the 5′ss inversely correlated with the presence of a downstream PPT, indicative of the functional role of the PPT. Finally, we found that the proteins that mediate this effect, TIA and U1C, and in particular their functional domains, are highly conserved across evolution. Overall, these findings expand our understanding of the role of TIA1/TIAR proteins in enhancing recognition of exons, in general, and Alu exons, in particular.

Human genes are composed of functional regions, termed exons, separated by non-functional regions, termed introns. Intronic sequences may gradually accumulate mutations and subsequently become recognized by the splicing machinery as exons, a process termed exonization. Alu elements are prone to undergo exonization: more than 5% of alternatively spliced internal exons in the human genome originate from Alu elements. A typical Alu element is ∼300 nucleotides long, consisting of two arms separated by a polypyrimdine tract (PPT). Interestingly, in most cases, exonization occurs almost exclusively within either the right arm or the left, not both. Here we found that the PPT between the two arms serves as a binding site for TIA proteins and prevents the exon selection process from expanding into downstream regions. To obtain a wider overview of TIA function, we performed a cross-evolutionary analysis within 22 eukaryotes of this protein and of U1C, a protein known to interact with it, and found that functional regions of both these proteins were highly conserved. These findings highlight the pivotal role of TIA proteins in 5′ splice-site selection of Alu exons and exon recognition in general.

SnoRNAs from the filamentous fungus Neurospora crassa: structural, functional and evolutionary insights

BACKGROUND

SnoRNAs represent an excellent model for studying the structural and functional evolution of small non-coding RNAs involved in the post-transcriptional modification machinery for rRNAs and snRNAs in eukaryotic cells. Identification of snoRNAs from Neurospora crassa, an important model organism playing key roles in the development of modern genetics, biochemistry and molecular biology will provide insights into the evolution of snoRNA genes in the fungus kingdom.

RESULTS

Fifty five box C/D snoRNAs were identified and predicted to guide 71 2'-O-methylated sites including four sites on snRNAs and three sites on tRNAs. Additionally, twenty box H/ACA snoRNAs, which potentially guide 17 pseudouridylations on rRNAs, were also identified. Although not exhaustive, the study provides the first comprehensive list of two major families of snoRNAs from the filamentous fungus N. crassa. The independently transcribed strategy dominates in the expression of box H/ACA snoRNA genes, whereas most of the box C/D snoRNA genes are intron-encoded. This shows that different genomic organizations and expression modes have been adopted by the two major classes of snoRNA genes in N. crassa . Remarkably, five gene clusters represent an outstanding organization of box C/D snoRNA genes, which are well conserved among yeasts and multicellular fungi, implying their functional importance for the fungus cells. Interestingly, alternative splicing events were found in the expression of two polycistronic snoRNA gene hosts that resemble the UHG-like genes in mammals. Phylogenetic analysis further revealed that the extensive separation and recombination of two functional elements of snoRNA genes has occurred during fungus evolution.

CONCLUSION

This is the first genome-wide analysis of the filamentous fungus N. crassa snoRNAs that aids in understanding the differences between unicellular fungi and multicellular fungi. As compared with two yeasts, a more complex pattern of methylation guided by box C/D snoRNAs in multicellular fungus than in unicellular yeasts was revealed, indicating the high diversity of post-transcriptional modification guided by snoRNAs in the fungus kingdom.

Divergence of exonic splicing elements after gene duplication and the impact on gene structures

An analysis of human exonic splicing elements in duplicated genes reveals their important role in the generation of new gene structures.

Background

The origin of new genes and their contribution to functional novelty has been the subject of considerable interest. There has been much progress in understanding the mechanisms by which new genes originate. Here we examine a novel way that new gene structures could originate, namely through the evolution of new alternative splicing isoforms after gene duplication.

Results

We studied the divergence of exonic splicing enhancers and silencers after gene duplication and the contributions of such divergence to the generation of new splicing isoforms. We found that exonic splicing enhancers and exonic splicing silencers diverge especially fast shortly after gene duplication. About 10% and 5% of paralogous exons undergo significantly asymmetric evolution of exonic splicing enhancers and silencers, respectively. When compared to pre-duplication ancestors, we found that there is a significant overall loss of exonic splicing enhancers and the magnitude increases with duplication age. Detailed examination reveals net gains and losses of exonic splicing enhancers and silencers in different copies and paralog clusters after gene duplication. Furthermore, we found that exonic splicing enhancer and silencer changes are mainly caused by synonymous mutations, though nonsynonymous changes also contribute. Finally, we found that exonic splicing enhancer and silencer divergence results in exon splicing state transitions (from constitutive to alternative or vice versa), and that the proportion of paralogous exon pairs with different splicing states also increases over time, consistent with previous predictions.

Conclusions

Our results suggest that exonic splicing enhancer and silencer changes after gene duplication have important roles in alternative splicing divergence and that these changes contribute to the generation of new gene structures.

Association of IL-15 polymorphisms with bone mineral density in postmenopausal Korean women

Interleukin-15 (IL-15) has been suggested to participate in bone metabolism by stimulating osteoclast differentiation and mediating inflammatory bone loss. This study investigated the effect of IL-15 gene polymorphisms on the bone mineral density (BMD) and bone fracture rates of postmenopausal women. Sequencing of the IL-15 gene in 24 Koreans revealed 16 single-nucleotide polymorphisms (SNPs), of which five were selected for further study. Postmenopausal Korean women (n = 844) were genotyped for these SNPs, and their BMDs and risk of fractures were assessed. It was found that the +20A > G, +13467C > A, +13653A > T, and +13815A > T IL-15 gene polymorphisms were significantly associated with the BMD of the lumbar spine and femoral neck and that their effects were gene-dose dependent. BMD was reduced when the minor allele of +13467A and +13653T or the common allele of +20A and +13815A was present. Haplotype (ht) analyses revealed that ht1 (GCAT) and ht2 (AATA) were associated with BMD of the lumbar spine and femoral neck. However, there was no association between the risk of fracture and IL-15 SNPs or hts. These results suggest that the +20A > G, +13467C > A, +13653A > T, and +13815A > T SNPs in the IL-15 gene affect BMD and, thus, could be genetic markers of osteoporosis.

Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3

Background

Alternative splicing (AS) of maturing mRNA can generate structurally and functionally distinct transcripts from the same gene. Recent bioinformatic analyses of available genome databases inferred a positive correlation between intron length and AS. To study the interplay between intron length and AS empirically and in more detail, we analyzed the diversity of alternatively spliced transcripts (ASTs) in the Drosophila RNA-binding Bruno-3 (Bru-3) gene. This gene was known to encode thirteen exons separated by introns of diverse sizes, ranging from 71 to 41,973 nucleotides in D. melanogaster. Although Bru-3's structure is expected to be conducive to AS, only two ASTs of this gene were previously described.

Results

Cloning of RT-PCR products of the entire ORF from four species representing three diverged Drosophila lineages provided an evolutionary perspective, high sensitivity, and long-range contiguity of splice choices currently unattainable by high-throughput methods. Consequently, we identified three new exons, a new exon fragment and thirty-three previously unknown ASTs of Bru-3. All exon-skipping events in the gene were mapped to the exons surrounded by introns of at least 800 nucleotides, whereas exons split by introns of less than 250 nucleotides were always spliced contiguously in mRNA. Cases of exon loss and creation during Bru-3 evolution in Drosophila were also localized within large introns. Notably, we identified a true de novo exon gain: exon 8 was created along the lineage of the obscura group from intronic sequence between cryptic splice sites conserved among all Drosophila species surveyed. Exon 8 was included in mature mRNA by the species representing all the major branches of the obscura group. To our knowledge, the origin of exon 8 is the first documented case of exonization of intronic sequence outside vertebrates.

Conclusion

We found that large introns can promote AS via exon-skipping and exon turnover during evolution likely due to frequent errors in their removal from maturing mRNA. Large introns could be a reservoir of genetic diversity, because they have a greater number of mutable sites than short introns. Taken together, gene structure can constrain and/or promote gene evolution.

A universal nonmonotonic relationship between gene compactness and expression levels in multicellular eukaryotes

Analysis of gene architecture and expression levels of four organisms, Homo sapiens, Caenorhabditis elegans, Drosophila melanogaster, and Arabidopsis thaliana, reveals a surprising, nonmonotonic, universal relationship between expression level and gene compactness. With increasing expression level, the genes tend at first to become longer but, from a certain level of expression, they become more and more compact, resulting in an approximate bell-shaped dependence. There are two leading hypotheses to explain the compactness of highly expressed genes. The selection hypothesis predicts that gene compactness is predominantly driven by the level of expression, whereas the genomic design hypothesis predicts that expression breadth across tissues is the driving force. We observed the connection between gene expression breadth in humans and gene compactness to be significantly weaker than the connection between expression level and compactness, a result that is compatible with the selection hypothesis but not the genome design hypothesis. The initial gene elongation with increasing expression level could be explained, at least in part, by accumulation of regulatory elements enhancing expression, in particular, in introns. This explanation is compatible with the observed positive correlation between intron density and expression level of a gene. Conversely, the trend toward increasing compactness for highly expressed genes could be caused by selection for minimization of energy and time expenditure during transcription and splicing and for increased fidelity of transcription, splicing, and/or translation that is likely to be particularly critical for highly expressed genes. Regardless of the exact nature of the forces that shape the gene architecture, we present evidence that, at least, in animals, coding and noncoding parts of genes show similar architectonic trends.

Characterization of a novel human HMBOX1 splicing variant lacking the homeodomain and with attenuated transcription repressor activity

We previously identified the human HMBOX1 (homeobox containing 1) gene, which represents a distinct group of HNF transcription factors, and is a potential transcription repressor with highly expression in pancreas. In our present work, we found that there exists a splicing variant of HMBOX1, designated HMBOX1b. Compared to HMBOX1, HMBOX1b encodes a 304 amino acids protein that shares the N-terminal region but has no homeodomain and the C-terminal region because of an alternative exon 7 which results in reading frame shifting. Unlike the highly pancreatic expression of HMBOX1, HMBOX1b was ubiquitous expressed in all human tissues detected by RT-PCR. Immunofluorescence staining showed that HMBOX1b accumulated in both cytoplasm and nucleus, and transcriptional reporter assays indicated that HMBOX1b only retained faint transcriptional repressive activity. Taken together, our findings suggest a distinct role of HMBOX1b, and the control of mRNA splicing might be involved in homeobox genes regulation.

Effects of splice sites on the intron retention in histamine H3 receptors from rats and mice

In the alternative splicing, intron retention, of histamine H(3) receptors in rats and mice, the short transcript isoforms that are excised alternatively spliced introns are easily detected in a very low level in rats and are undetectable in mice using the regular PCR protocol. The retained introns have common 5' splice site and different 3' splice sites. The detailed mechanism for the special alternative splicing remains largely unclear. In this study, we developed a minigene splicing system to recapitulate natural alternative splicing of the receptors and investigated the effects of 5' and 3' splice sites on intron retention in HeLa cells. Mutating weak 5' and 3' splice sites of the alternatively spliced introns toward the canonical consensus sequences promoted the splicing of the corresponding introns in rat and mouse minigenes. The effect of splice site strength was context-dependent and much more significant for the 3' splice site of the longer alternative intron than for the 3' splice site of the shorter alternative intron and the common 5' splice sites; it was also more significant in the rat minigene than in the mouse minigene. Mutating the 3' splice site of the longer alternative intron resulted in almost complete splicing of the intron and made the corresponding isoform to become the nearly exclusive transcript in the rat minigene.