Alternative splicing and evolution: diversification, exon definition and function

A transcript perspective on evolution

Alternative splicing is now recognized as a major mechanism for transcriptome and proteome diversity in higher eukaryotes, yet its evolution is poorly understood. Most studies focus on the evolution of exons and introns at the gene level, while only few consider the evolution of transcripts. In this paper, we present a framework for transcript phylogenies where ancestral transcripts evolve along the gene tree by gains, losses, and mutation. We demonstrate the usefulness of our method on a set of 805 genes and two different topics. First, we improve a method for transcriptome reconstruction from ESTs (ASPic), then we study the evolution of function in transcripts. The use of transcript phylogenies allows us to double the precision of ASPic, whereas results on the functional study reveal that conserved transcripts are more likely to share protein domains than functional sites. These studies validate our framework for the study of evolution in large collections of organisms from the perspective of transcripts; for this purpose, we developed and provide a new tool, TrEvoR.

Regulation of vascular function on posttranscriptional level

Posttranscriptional control of gene expression is crucial for regulating plurality of proteins and functional plasticity of the proteome under (patho)physiologic conditions. Alternative splicing as well as micro (mi)RNA-mediated mechanisms play an important role for the regulation of protein expression on posttranscriptional level. Both alternative splicing and miRNAs were shown to influence cardiovascular functions, such as endothelial thrombogenicity and the vascular tone, by regulating the expression of several vascular proteins and their isoforms, such as Tissue Factor (TF) or the endothelial nitric oxide synthase (eNOS). This review will summarize and discuss the latest findings on the (patho)physiologic role of alternative splicing processes as well as of miRNAs on modulation of vascular functions, such as coagulation, thrombosis, and regulation of the vascular tone.

Reactivation of fetal splicing programs in diabetic hearts is mediated by protein kinase C signaling

Diabetic cardiomyopathy is one of the complications of diabetes that eventually leads to heart failure and death. Aberrant activation of PKC signaling contributes to diabetic cardiomyopathy by mechanisms that are poorly understood. Previous reports indicate that PKC is implicated in alternative splicing regulation. Therefore, we wanted to test whether PKC activation in diabetic hearts induces alternative splicing abnormalities. Here, using RNA sequencing we identified a set of 22 alternative splicing events that undergo a developmental switch in splicing, and we confirmed that splicing reverts to an embryonic pattern in adult diabetic hearts. This network of genes has important functions in RNA metabolism and in developmental processes such as differentiation. Importantly, PKC isozymes α/β control alternative splicing of these genes via phosphorylation and up-regulation of the RNA-binding proteins CELF1 and Rbfox2. Using a mutant of CELF1, we show that phosphorylation of CELF1 by PKC is necessary for regulation of splicing events altered in diabetes. In summary, our studies indicate that activation of PKCα/β in diabetic hearts contributes to the genome-wide splicing changes through phosphorylation and up-regulation of CELF1/Rbfox2 proteins. These findings provide a basis for PKC-mediated cardiac pathogenesis under diabetic conditions.

Released selective pressure on a structural domain gives new insights on the functional relaxation of mitochondrial aspartyl-tRNA synthetase

Mammalian mitochondrial aminoacyl-tRNA synthetases are nuclear-encoded enzymes that are essential for mitochondrial protein synthesis. Due to an endosymbiotic origin of the mitochondria, many of them share structural domains with homologous bacterial enzymes of same specificity. This is also the case for human mitochondrial aspartyl-tRNA synthetase (AspRS) that shares the so-called bacterial insertion domain with bacterial homologs. The function of this domain in the mitochondrial proteins is unclear. Here, we show by bioinformatic analyses that the sequences coding for the bacterial insertion domain are less conserved in opisthokont and protist than in bacteria and viridiplantae. The divergence suggests a loss of evolutionary pressure on this domain for non-plant mitochondrial AspRSs. This discovery is further connected with the herein described occurrence of alternatively spliced transcripts of the mRNAs coding for some mammalian mitochondrial AspRSs. Interestingly, the spliced transcripts alternately lack one of the four exons that code for the bacterial insertion domain. Although we showed that the human alternative transcript is present in all tested tissues; co-exists with the full-length form, possesses 5'- and 3'-UTRs, a poly-A tail and is bound to polysomes, we were unable to detect the corresponding protein. The relaxed selective pressure combined with the occurrence of alternative splicing, involving a single structural sub-domain, favors the hypothesis of the loss of function of this domain for AspRSs of mitochondrial location. This evolutionary divergence is in line with other characteristics, established for the human mt-AspRS, that indicate a functional relaxation of non-viridiplantae mt-AspRSs when compared to bacterial and plant ones, despite their common ancestry.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

New gene evolution: little did we know

Genes are perpetually added to and deleted from genomes during evolution. Thus, it is important to understand how new genes are formed and how they evolve to be critical components of the genetic systems that determine the biological diversity of life. Two decades of effort have shed light on the process of new gene origination and have contributed to an emerging comprehensive picture of how new genes are added to genomes, ranging from the mechanisms that generate new gene structures to the presence of new genes in different organisms to the rates and patterns of new gene origination and the roles of new genes in phenotypic evolution. We review each of these aspects of new gene evolution, summarizing the main evidence for the origination and importance of new genes in evolution. We highlight findings showing that new genes rapidly change existing genetic systems that govern various molecular, cellular, and phenotypic functions.

The splicing landscape is globally reprogrammed during male meiosis

Meiosis requires conserved transcriptional changes, but it is not known whether there is a corresponding set of RNA splicing switches. Here, we used RNAseq of mouse testis to identify changes associated with the progression from mitotic spermatogonia to meiotic spermatocytes. We identified ∼150 splicing switches, most of which affect conserved protein-coding exons. The expression of many key splicing regulators changed in the course of meiosis, including downregulation of polypyrimidine tract binding protein (PTBP1) and heterogeneous nuclear RNP A1, and upregulation of nPTB, Tra2β, muscleblind, CELF proteins, Sam68 and T-STAR. The sequences near the regulated exons were significantly enriched in target sites for PTB, Tra2β and STAR proteins. Reporter minigene experiments investigating representative exons in transfected cells showed that PTB binding sites were critical for splicing of a cassette exon in the Ralgps2 mRNA and a shift in alternative 5′ splice site usage in the Bptf mRNA. We speculate that nPTB might functionally replace PTBP1 during meiosis for some target exons, with changes in the expression of other splicing factors helping to establish meiotic splicing patterns. Our data suggest that there are substantial changes in the determinants and patterns of alternative splicing in the mitotic-to-meiotic transition of the germ cell cycle.

Transcriptomic characterization of cold acclimation in larval zebrafish

Background

Temperature is one of key environmental parameters that affect the whole life of fishes and an increasing number of studies have been directed towards understanding the mechanisms of cold acclimation in fish. However, the adaptation of larvae to cold stress and the cold-specific transcriptional alterations in fish larvae remain largely unknown. In this study, we characterized the development of cold-tolerance in zebrafish larvae and investigated the transcriptional profiles under cold stress using RNA-seq.

Results

Pre-exposure of 96 hpf zebrafish larvae to cold stress (16°C) for 24 h significantly increased their survival rates under severe cold stress (12°C). RNA-seq generated 272 million raw reads from six sequencing libraries and about 92% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 1,431 up- and 399 down-regulated genes. Gene ontology enrichment analysis of cold-induced genes revealed that RNA splicing, ribosome biogenesis and protein catabolic process were the most highly overrepresented biological processes. Spliceosome, proteasome, eukaryotic ribosome biogenesis and RNA transport were the most highly enriched pathways for genes up-regulated by cold stress. Moreover, alternative splicing of 197 genes and promoter switching of 64 genes were found to be regulated by cold stress. A shorter isoform of stk16 that lacks 67 amino acids at the N-terminus was specifically generated by skipping the second exon in cold-treated larvae. Alternative promoter usage was detected for per3 gene under cold stress, which leading to a highly up-regulated transcript encoding a truncated protein lacking the C-terminal domains.

Conclusions

These findings indicate that zebrafish larvae possess the ability to build cold-tolerance under mild low temperature and transcriptional and post-transcriptional regulations are extensively involved in this acclimation process.

Alternative RNA splicing and cancer

Alternative splicing of pre-messenger RNA (mRNA) is a fundamental mechanism by which a gene can give rise to multiple distinct mRNA transcripts, yielding protein isoforms with different, even opposing, functions. With the recognition that alternative splicing occurs in nearly all human genes, its relationship with cancer-associated pathways has emerged as a rapidly growing field. In this review, we summarize recent findings that have implicated the critical role of alternative splicing in cancer and discuss current understandings of the mechanisms underlying dysregulated alternative splicing in cancer cells.

Identification of a novel alternative splicing variant of hemocyanin from shrimp Litopenaeus vannamei

Recent evidences suggest that invertebrates express families of immune molecules with high levels of sequence diversity. Hemocyanin is an important non-specific immune molecule present in the hemolymph of both mollusks and arthropods. In the present study, we characterized a novel alternative splicing variant of hemocyanin (cHE1) from Litopenaeus vannamei that produced mRNA transcript of 2579 bp in length. The isoform contained two additional sequences of 296 and 267 bp in the 5'- and 3'-terminus respectively, in comparison to that of wild type hemocyanin (cHE). Sequence of cHE1 shows 100% identity to that of hemocyanin genomic DNA (HE, which does not form an open reading frame), suggesting that cHE1 might be an alternative splicing variant due to intron retention. Moreover, cHE1 could be detected by RT-PCR from five tissues (heart, gill, stomach, intestine and brain), and from shrimps at stages from nauplius to mysis larva. Further, cHE1 mRNA transcripts were significantly increased in hearts after 12h of infection with Vibrio parahemolyticus or poly I: C, while no significant difference in the transcript levels of hepatopancreas cHE was detected in the pathogen-treated shrimps during the period. In summary, these studies suggested a novel splicing variant of hemocyanin in shrimp, which might be involved in shrimp resistance to pathogenic infection.

Microsatellite-encoded domain in rodent Sry functions as a genetic capacitor to enable the rapid evolution of biological novelty

The male program of therian mammals is determined by Sry, a transcription factor encoded by the Y chromosome. Specific DNA binding is mediated by a high mobility group (HMG) box. Expression of Sry in the gonadal ridge activates a Sox9-dependent gene regulatory network leading to testis formation. A subset of Sry alleles in superfamily Muroidea (order Rodentia) is remarkable for insertion of an unstable DNA microsatellite, most commonly encoding (as in mice) a CAG repeat-associated glutamine-rich domain. We provide evidence, based on an embryonic pre-Sertoli cell line, that this domain functions at a threshold length as a genetic capacitor to facilitate accumulation of variation elsewhere in the protein, including the HMG box. The glutamine-rich domain compensates for otherwise deleterious substitutions in the box and absence of nonbox phosphorylation sites to ensure occupancy of DNA target sites. Such compensation enables activation of a male transcriptional program despite perturbations to the box. Whereas human SRY requires nucleocytoplasmic shuttling and coupled phosphorylation, mouse Sry contains a defective nuclear export signal analogous to a variant human SRY associated with inherited sex reversal. We propose that the rodent glutamine-rich domain has (i) fostered accumulation of cryptic intragenic variation and (ii) enabled unmasking of such variation due to DNA replicative slippage. This model highlights genomic contingency as a source of protein novelty at the edge of developmental ambiguity and may underlie emergence of non-Sry-dependent sex determination in the radiation of Muroidea.

Shark IgW C region diversification through RNA processing and isotype switching

Sharks and skates represent the earliest vertebrates with an adaptive immune system based on lymphocyte Ag receptors generated by V(D)J recombination. Shark B cells express two classical Igs, IgM and IgW, encoded by an early, alternative gene organization consisting of numerous autonomous miniloci, where the individual gene cluster carries a few rearranging gene segments and one C region, μ or ω. We have characterized eight distinct Ig miniloci encoding the nurse shark ω H chain. Each cluster consists of VH, D, and JH segments and six to eight C domain exons. Two interspersed secretory exons, in addition to the 3'-most C exon with tailpiece, provide the gene cluster with the ability to generate at least six secreted isoforms that differ as to polypeptide length and C domain combination. All clusters appear to be functional, as judged by the capability for rearrangement and absence of defects in the deduced amino acid sequence. We previously showed that IgW VDJ can perform isotype switching to μ C regions; in this study, we found that switching also occurs between ω clusters. Thus, C region diversification for any IgW VDJ can take place at the DNA level by switching to other ω or μ C regions, as well as by RNA processing to generate different C isoforms. The wide array of pathogens recognized by Abs requires different disposal pathways, and our findings demonstrate complex and unique pathways for C effector function diversity that evolved independently in cartilaginous fishes.

Alternative pre-mRNA splicing in neurons: growing up and extending its reach

Alternative pre-mRNA splicing determines the protein output of most neuronally expressed genes. Many examples have been described of protein function being modulated by coding changes in different mRNA isoforms. Several recent studies demonstrate that, through the coupling of splicing to other processes of mRNA metabolism, alternative splicing can also act as an on/off switch for gene expression. Other regulated splicing events may determine how an mRNA is utilized in its later cytoplasmic life by changing its localization or translation. These studies make clear that the multiple steps of post-transcriptional gene regulation are strongly linked. Together, these regulatory process play key roles in all aspects of the cell biology of neurons, from their initial differentiation, to their choice of connections, and finally to their function with mature circuits.

Lariat sequencing in a unicellular yeast identifies regulated alternative splicing of exons that are evolutionarily conserved with humans

Alternative splicing is a potent regulator of gene expression that vastly increases proteomic diversity in multicellular eukaryotes and is associated with organismal complexity. Although alternative splicing is widespread in vertebrates, little is known about the evolutionary origins of this process, in part because of the absence of phylogenetically conserved events that cross major eukaryotic clades. Here we describe a lariat-sequencing approach, which offers high sensitivity for detecting splicing events, and its application to the unicellular fungus, Schizosaccharomyces pombe, an organism that shares many of the hallmarks of alternative splicing in mammalian systems but for which no previous examples of exon-skipping had been demonstrated. Over 200 previously unannotated splicing events were identified, including examples of regulated alternative splicing. Remarkably, an evolutionary analysis of four of the exons identified here as subject to skipping in S. pombe reveals high sequence conservation and perfect length conservation with their homologs in scores of plants, animals, and fungi. Moreover, alternative splicing of two of these exons have been documented in multiple vertebrate organisms, making these the first demonstrations of identical alternative-splicing patterns in species that are separated by over 1 billion y of evolution.

Evolution of the plasma and tissue kallikreins, and their alternative splicing isoforms

Kallikreins are secreted serine proteases with important roles in human physiology. Human plasma kallikrein, encoded by the KLKB1 gene on locus 4q34-35, functions in the blood coagulation pathway, and in regulating blood pressure. The human tissue kallikrein and kallikrein-related peptidases (KLKs) have diverse expression patterns and physiological roles, including cancer-related processes such as cell growth regulation, angiogenesis, invasion, and metastasis. Prostate-specific antigen (PSA), the product of the KLK3 gene, is the most widely used biomarker in clinical practice today. A total of 15 KLKs are encoded by the largest contiguous cluster of protease genes in the human genome (19q13.3-13.4), which makes them ideal for evolutionary analysis of gene duplication events. Previous studies on the evolution of KLKs have traced mammalian homologs as well as a probable early origin of the family in aves, amphibia and reptilia. The aim of this study was to address the evolutionary and functional relationships between tissue KLKs and plasma kallikrein, and to examine the evolution of alternative splicing isoforms. Sequences of plasma and tissue kallikreins and their alternative transcripts were collected from the NCBI and Ensembl databases, and comprehensive phylogenetic analysis was performed by Bayesian as well as maximum likelihood methods. Plasma and tissue kallikreins exhibit high sequence similarity in the trypsin domain (>50%). Phylogenetic analysis indicates an early divergence of KLKB1, which groups closely with plasminogen, chymotrypsin, and complement factor D (CFD), in a monophyletic group distinct from trypsin and the tissue KLKs. Reconstruction of the earliest events leading to the diversification of the tissue KLKs is not well resolved, indicating rapid expansion in mammals. Alternative transcripts of each KLK gene show species-specific divergence, while examination of sequence conservation indicates that many annotated human KLK isoforms are missing the catalytic triad that is crucial for protease activity.

Molecular characterization, alternative splicing and expression analysis of bovine DBC1

Deleted in breast cancer 1 (DBC1, KIAA1967, p30 DBC) is a novel protein that has been recently shown to bind and regulate SIRT1. Loss of function of DBC1 increased SIRT1 deacetylase activity, which promotes "browning" of WAT by deacetylating peroxisome proliferator-activated receptor (PPARγ) on Lys268 and Lys293. In the present study, we have cloned and characterized the bovine DBC1 gene. Two transcript variants of bovine DBC1 were identified, designated DBC1-A and DBC1-B, respectively, which were both located in nucleus. Protein sequence analysis showed that DBC1-A was well conserved across species. Expression analysis of DBC1 in seven different tissues of calves and bulls by RT-PCR indicated that the two transcripts were ubiquitously expressed. However, the relatively level of DBC1-A was higher when compared to DBC1-B in all examined tissues. Surprisingly, the expression of DBC1-A was extraordinary high in calves adipose tissue, which implicated its potential key role in regulating calve adipocyte development. These findings provide new insight into our understanding of the biochemical characteristics and physiological role of bovine DBC1.

A compensatory effect upon splicing results in normal function of the CYP2A6*14 allele

A synonymous variant in the first exon of CYP2A6, rs1137115 (51G>A), defines the common reference allele CYP2A6*1A, and is associated with lower mRNA expression and slower in-vivo nicotine metabolism. Another common allele, CYP2A6*14, differs from CYP2A6*1A by a single variant, rs28399435 (86G>A, S29N). However, CYP2A6*14 shows in-vivo activity comparable with that of full-function alleles, and significantly higher than CYP2A6*1A. rs1137115A is predicted to create an exonic splicing suppressor site overlapping an exonic splicing enhancer (ESE) site in the first exon of CYP2A6, whereas rs28399435A is predicted to strengthen another adjacent ESE, potentially compensating for rs1137115A. Using an allelic expression assay to assess cDNAs produced from rs1137115 heterozygous liver biopsy samples, lower expression of the CYP2A6*1A allele is confirmed while CYP2A6*14 expression is found to be indistinguishable from that of rs1137115G alleles. Quantitative PCR assays to determine the relative abundance of spliced and unspliced or partially spliced CYP2A6 mRNAs in liver biopsy samples show that *1A/*1A homozygotes have a significantly lower ratio, due to both a reduction in spliced forms and an increase in unspliced or partially spliced CYP2A6. These results show the importance of common genetic variants that effect exonic splicing suppressor and ESEs to explain human variation regarding clinically-relevant phenotypes.

Splice variants and promoter methylation status of the Bovine Vasa Homology (Bvh) gene may be involved in bull spermatogenesis

Background

Vasa is a member of the DEAD-box protein family that plays an indispensable role in mammalian spermatogenesis, particularly during meiosis. Bovine vasa homology (Bvh) of Bos taurus has been reported, however, its function in bovine testicular tissue remains obscure. This study aimed to reveal the functions of Bvh and to determine whether Bvh is a candidate gene in the regulation of spermatogenesis in bovine, and to illustrate whether its transcription is regulated by alternative splicing and DNA methylation.

Results

Here we report the molecular characterization, alternative splicing pattern, expression and promoter methylation status of Bvh. The full-length coding region of Bvh was 2190 bp, which encodes a 729 amino acid (aa) protein containing nine consensus regions of the DEAD box protein family. Bvh is expressed only in the ovary and testis of adult cattle. Two splice variants were identified and termed Bvh-V4 (2112 bp and 703 aa) and Bvh-V45 (2040 bp and 679 aa). In male cattle, full-length Bvh (Bvh-FL), Bvh-V4 and Bvh-V45 are exclusively expressed in the testes in the ratio of 2.2:1.6:1, respectively. Real-time PCR revealed significantly reduced mRNA expression of Bvh-FL, Bvh-V4 and Bvh-V45 in testes of cattle-yak hybrids, with meiotic arrest compared with cattle and yaks with normal spermatogenesis (P < 0.01). The promoter methylation level of Bvh in the testes of cattle-yak hybrids was significantly greater than in cattle and yaks (P < 0.01).

Conclusion

In the present study, Bvh was isolated and characterized. These data suggest that Bvh functions in bovine spermatogenesis, and that transcription of the gene in testes were regulated by alternative splice and promoter methylation.

Features of evolutionarily conserved alternative splicing events between Brassica and Arabidopsis

Alternative splicing (AS) generates multiple types of mRNA from a single type of pre-mRNA by differential intron splicing. It can result in new protein isoforms or down-regulation of gene expression by transcript decay. The evolutionary conservation of AS events in plants is largely unexplored and only a small number of AS events have been identified as conserved between divergent species. We performed a large-scale analysis of cDNA data from Brassica and Arabidopsis to identify and further characterize conserved AS events. We identified 537 conserved AS events in 485 genes. Alternative donor and acceptor events are significantly overrepresented among conserved events, whereas intron retention and exon skipping events are underrepresented. Conserved AS events are significantly shorter, less likely to be in the 3'UTR, and they are enriched for genes whose products function in the chloroplast. AS modified a functional domain for about half of the genes with conserved events. We further characterized three genes with conserved AS events. This study identifies many AS events that are conserved between Brassica and Arabidopsis, revealing features of conserved AS events. Many of the conserved AS events may have important, but uncharacterized, functions.

Hypoxia-dependent mRNA expression pattern of splicing factor YT521 and its impact on oncological important target gene expression

The ubiquitously expressed splicing factor YT521 (YTHDC1) is characterized by alternatively spliced isoforms with regulatory impact on cancer-associated gene expression. Our recent findings account for the prognostic significance of YT521 in endometrial cancer. In this study, we investigated the hypoxia-dependency of YT521 expression as well as its differential isoform activities on oncological important target genes. YT521's potential regulatory influence on splicing was investigated by a minigene assay for the specific target gene CD44. Functional splicing analysis was performed by YT521 knock-down or overexpression, respectively. In addition, YT521 expression was determined under hypoxia. The two protein-generating YT521 mRNA isoforms 1 and 2 caused a comparable, specific induction of CD44v alternative splicing (P < 0.01). In a number of oncological target genes, YT521 upregulation significantly altered BRCA2 expression pattern, while YT521 knock-down created a significant regulatory impact on PGR expression, respectively. Hypoxia induced a specific switch towards the processing of two non-protein-coding mRNA variants, of which one is described for the first time in this study. The presented study underlines the comparable regulatory potential of both YT521 isoforms 1 and 2, on the investigated target genes in vivo and in vitro. Hypoxia induces a specific switch in YT521 expression pattern towards the two non-protein coding mRNA variants, the already characterized isoform 3 and the newly discovered exon 8-skipping isoform. The altered YT521 alternative splicing is functionally coupled with nonsense-mediated decay and can be interpreted as regulated unproductive splicing and transcription with consecutive impact on the processing of specific cancer-associated genes, such as BRCA2 and PGR.

TWEAK-independent Fn14 self-association and NF-κB activation is mediated by the C-terminal region of the Fn14 cytoplasmic domain

The tumor necrosis factor (TNF) superfamily member TNF-like weak inducer of apoptosis (TWEAK) is a pro-inflammatory and pro-angiogenic cytokine implicated in physiological tissue regeneration and wound repair. TWEAK binds to a 102-amino acid type I transmembrane cell surface receptor named fibroblast growth factor-inducible 14 (Fn14). TWEAK:Fn14 engagement activates several intracellular signaling cascades, including the NF-κB pathway, and sustained Fn14 signaling has been implicated in the pathogenesis of chronic inflammatory diseases and cancer. Although several groups are developing TWEAK- or Fn14-targeted agents for therapeutic use, much more basic science research is required before we fully understand the TWEAK/Fn14 signaling axis. For example, we and others have proposed that TWEAK-independent Fn14 signaling may occur in cells when Fn14 levels are highly elevated, but this idea has never been tested directly. In this report, we first demonstrate TWEAK-independent Fn14 signaling by showing that an Fn14 deletion mutant that is unable to bind TWEAK can activate the NF-κB pathway in transfected cells. We then show that ectopically-expressed, cell surface-localized Fn14 can self-associate into Fn14 dimers, and we show that Fn14 self-association is mediated by an 18-aa region within the Fn14 cytoplasmic domain. Endogenously-expressed Fn14 as well as ectopically-overexpressed Fn14 could also be detected in dimeric form when cell lysates were subjected to SDS-PAGE under non-reducing conditions. Additional experiments revealed that Fn14 dimerization occurs during cell lysis via formation of an intermolecular disulfide bond at cysteine residue 122. These findings provide insight into the Fn14 signaling mechanism and may aid current studies to develop therapeutic agents targeting this small cell surface receptor.

From the ultrasonic to the infrared: molecular evolution and the sensory biology of bats

Great advances have been made recently in understanding the genetic basis of the sensory biology of bats. Research has focused on the molecular evolution of candidate sensory genes, genes with known functions [e.g., olfactory receptor (OR) genes] and genes identified from mutations associated with sensory deficits (e.g., blindness and deafness). For example, the FoxP2 gene, underpinning vocal behavior and sensorimotor coordination, has undergone diversification in bats, while several genes associated with audition show parallel amino acid substitutions in unrelated lineages of echolocating bats and, in some cases, in echolocating dolphins, representing a classic case of convergent molecular evolution. Vision genes encoding the photopigments rhodopsin and the long-wave sensitive opsin are functional in bats, while that encoding the short-wave sensitive opsin has lost functionality in rhinolophoid bats using high-duty cycle laryngeal echolocation, suggesting a sensory trade-off between investment in vision and echolocation. In terms of olfaction, bats appear to have a distinctive OR repertoire compared with other mammals, and a gene involved in signal transduction in the vomeronasal system has become non-functional in most bat species. Bitter taste receptors appear to have undergone a "birth-and death" evolution involving extensive gene duplication and loss, unlike genes coding for sweet and umami tastes that show conservation across most lineages but loss in vampire bats. Common vampire bats have also undergone adaptations for thermoperception, via alternative splicing resulting in the evolution of a novel heat-sensitive channel. The future for understanding the molecular basis of sensory biology is promising, with great potential for comparative genomic analyses, studies on gene regulation and expression, exploration of the role of alternative splicing in the generation of proteomic diversity, and linking genetic mechanisms to behavioral consequences.

Evolution of the vertebrate bone matrix: an expression analysis of the network forming collagen paralogues in amphibian osteoblasts

The emergence of vertebrates is closely associated to the evolution of mineralized bone tissue. However, the molecular basis underlying the origin and subsequent diversification of the skeletal mineralized matrix is still poorly understood. One efficient way to tackle this issue is to compare the expression, between vertebrate species, of osteoblastic genes coding for bone matrix proteins. In this work, we have focused on the evolution of the network forming collagen family which contains the Col8a1, Col8a2, and Col10a1 genes. Both phylogeny and synteny reveal that these three paralogues are vertebrate-specific and derive from two independent duplications in the vertebrate lineage. To shed light on the evolution of this family, we have analyzed the osteoblastic expression of the network forming collagens in endochondral and intramembraneous skeletal elements of the amphibian Xenopus tropicalis. Remarkably, we find that amphibian osteoblasts express Col10a1, a gene strongly expressed in osteoblasts in actinopterygians but not in amniotes. In addition, while Col8a1 is known to be robustly expressed in mammalian osteoblasts, the expression levels of its amphibian orthologue are dramatically reduced. Our work reveals that while a skeletal expression of network forming collagen members is widespread throughout vertebrates, osteoblasts from divergent vertebrate lineages express different combinations of network forming collagen paralogues.

Global regulation of alternative splicing by adenosine deaminase acting on RNA (ADAR)

Alternative mRNA splicing is a major mechanism for gene regulation and transcriptome diversity. Despite the extent of the phenomenon, the regulation and specificity of the splicing machinery are only partially understood. Adenosine-to-inosine (A-to-I) RNA editing of pre-mRNA by ADAR enzymes has been linked to splicing regulation in several cases. Here we used bioinformatics approaches, RNA-seq and exon-specific microarray of ADAR knockdown cells to globally examine how ADAR and its A-to-I RNA editing activity influence alternative mRNA splicing. Although A-to-I RNA editing only rarely targets canonical splicing acceptor, donor, and branch sites, it was found to affect splicing regulatory elements (SREs) within exons. Cassette exons were found to be significantly enriched with A-to-I RNA editing sites compared with constitutive exons. RNA-seq and exon-specific microarray revealed that ADAR knockdown in hepatocarcinoma and myelogenous leukemia cell lines leads to global changes in gene expression, with hundreds of genes changing their splicing patterns in both cell lines. This global change in splicing pattern cannot be explained by putative editing sites alone. Genes showing significant changes in their splicing pattern are frequently involved in RNA processing and splicing activity. Analysis of recently published RNA-seq data from glioblastoma cell lines showed similar results. Our global analysis reveals that ADAR plays a major role in splicing regulation. Although direct editing of the splicing motifs does occur, we suggest it is not likely to be the primary mechanism for ADAR-mediated regulation of alternative splicing. Rather, this regulation is achieved by modulating trans-acting factors involved in the splicing machinery.

The differential transcription network between embryo and endosperm in the early developing maize seed

Transcriptome analysis of early-developing maize (Zea mays) seed was conducted using Illumina sequencing. We mapped 11,074,508 and 11,495,788 paired-end reads from endosperm and embryo, respectively, at 9 d after pollination to define gene structure and alternative splicing events as well as transcriptional regulators of gene expression to quantify transcript abundance in both embryo and endosperm. We identified a large number of novel transcribed regions that did not fall within maize annotated regions, and many of the novel transcribed regions were tissue-specifically expressed. We found that 50.7% (8,556 of 16,878) of multiexonic genes were alternatively spliced, and some transcript isoforms were specifically expressed either in endosperm or in embryo. In addition, a total of 46 trans-splicing events, with nine intrachromosomal events and 37 interchromosomal events, were found in our data set. Many metabolic activities were specifically assigned to endosperm and embryo, such as starch biosynthesis in endosperm and lipid biosynthesis in embryo. Finally, a number of transcription factors and imprinting genes were found to be specifically expressed in embryo or endosperm. This data set will aid in understanding how embryo/endosperm development in maize is differentially regulated.

Gene set enrichment analysis of RNA-Seq data: integrating differential expression and splicing

BACKGROUND

RNA-Seq has become a key technology in transcriptome studies because it can quantify overall expression levels and the degree of alternative splicing for each gene simultaneously. To interpret high-throughout transcriptome profiling data, functional enrichment analysis is critical. However, existing functional analysis methods can only account for differential expression, leaving differential splicing out altogether.

RESULTS

In this work, we present a novel approach to derive biological insight by integrating differential expression and splicing from RNA-Seq data with functional gene set analysis. This approach designated SeqGSEA, uses count data modelling with negative binomial distributions to first score differential expression and splicing in each gene, respectively, followed by two strategies to combine the two scores for integrated gene set enrichment analysis. Method comparison results and biological insight analysis on an artificial data set and three real RNA-Seq data sets indicate that our approach outperforms alternative analysis pipelines and can detect biological meaningful gene sets with high confidence, and that it has the ability to determine if transcription or splicing is their predominant regulatory mechanism.

CONCLUSIONS

By integrating differential expression and splicing, the proposed method SeqGSEA is particularly useful for efficiently translating RNA-Seq data to biological discoveries.

The impact of splicing on protein domain architecture

Many proteins are composed of protein domains, functional units of common descent. Multidomain forms are common in all eukaryotes making up more than half of the proteome and the evolution of novel domain architecture has been accelerated in metazoans. It is also becoming increasingly clear that alternative splicing is prevalent among vertebrates. Given that protein domains are defined as structurally, functionally and evolutionarily distinct units, one may speculate that some alternative splicing events may lead to clean excisions of protein domains, thus generating a number of different domain architectures from one gene template. However, recent findings indicate that smaller alternative splicing events, in particular in disordered regions, might be more prominent than domain architectural changes. The problem of identifying protein isoforms is, however, still not resolved. Clearly, many splice forms identified through detection of mRNA sequences appear to produce 'nonfunctional' proteins, such as proteins with missing internal secondary structure elements. Here, we review the state of the art methods for identification of functional isoforms and present a summary of what is known, thus far, about alternative splicing with regard to protein domain architectures.

Aberrant Alternative Splicing Events in Parkinson's Disease

Alternative splicing (AS) using a sole gene to express multiple transcripts with diverse protein coding sequences and/or RNA regulatory elements raises genomic complexities. In the nervous system, several thousand AS events play important roles in ion transportation, receptor recognition, neurotransmission, memory, and learning. Not surprisingly, AS influences human physiology, development, and disease. Many research studies have focused on aberrant AS in nervous system diseases, including Parkinson's disease (PD), the second most common progressive neurodegenerative disorder of the central nervous system. PD affects the lives of several million people globally. It is caused by protein aggregation, such as in Lewy bodies, and the loss of dopamine-containing neurons in the substantia nigra of the midbrain. To our knowledge, six genes, including PARK2, SNCAIP, LRRK2, SNCA, SRRM2, and MAPT, are involved in aberrant AS events in PD patients. In this review, we highlight the relevance of aberrant AS in PD and discuss the use of an aberrant AS profile as a potential diagnostic or prognostic marker for PD and as a possible means of applying therapy.

Chronic myelomonocytic leukemia: a review of the molecular biology, prognostic models and treatment

SUMMARY Chronic myelomonocytic leukemia (CMML) is a genetically heterogeneous hematologic neoplasm that manifests with features of both a myelodysplastic syndrome and a myeloproliferative neoplasm. Recent advances in the characterization of recurrent genetic markers have resulted in a better understanding of the leukemia-initiating events and have distinguished CMML from other clonal hematopoietic malignancies. Although these mutations may lead to CMML-specific therapies in the relatively near future, the current state of therapy for CMML is based on treatments designed for the myelodysplastic syndromes. Here we review the recurrent genetic mutations and, if known, their clinical significance. We also review the treatment and available CMML-specific prognostic models and novel therapies moving forward.

Integrator complex plays an essential role in adipose differentiation

The dynamic process of adipose differentiation involves stepwise expressions of transcription factors and proteins specific to the mature fat cell phenotype. In this study, it was revealed that expression levels of IntS6 and IntS11, subunits of the Integrator complex, were increased in 3T3-L1 cells in the period when the cells reached confluence and differentiated into adipocytes, while being reduced to basal levels after the completion of differentiation. Suppression of IntS6 or IntS11 expression using siRNAs in 3T3-L1 preadipocytes markedly inhibited differentiation into mature adipocytes, based on morphological findings as well as mRNA analysis of adipocyte-specific genes such as Glut4, perilipin and Fabp4. Although Pparγ2 protein expression was suppressed in IntS6 or IntS11-siRNA treated cells, adenoviral forced expression of Pparγ2 failed to restore the capacity for differentiation into mature adipocytes. Taken together, these findings demonstrate that increased expression of Integrator complex subunits is an indispensable event in adipose differentiation. Although further study is necessary to elucidate the underlying mechanism, the processing of U1, U2 small nuclear RNAs may be involved in cell differentiation steps.

Molecular and bioinformatical characterization of a novel superfamily of cysteine-rich peptides from arthropods

The full-length cDNA sequences of two novel cysteine-rich peptides (referred to as HsVx1 and MmKTx1) were obtained from scorpions. The two peptides represent a novel class of cysteine-rich peptides with a unique cysteine pattern. The genomic sequence of HsVx1 is composed of three exons interrupted by two introns that are localized in the mature peptide encoding region and inserted in phase 1 and phase 2, respectively. Such a genomic organization markedly differs from those of other peptides from scorpions described previously. Genome-wide search for the orthologs of HsVx1 identified 59 novel cysteine-rich peptides from arthropods. These peptides share a consistent cysteine pattern with HsVx1. Genomic comparison revealed extensive intron length differences and intronic number and position polymorphisms among the genes of these peptides. Further analysis identified 30 cases of intron sliding, 1 case of intron gain and 22 cases of intron loss occurred with the genes of the HsVx1 and HsVx1-like peptides. It is interesting to see that three HsVx1-like peptides XP_001658928, XP_001658929 and XP_001658930 were derived from a single gene (XP gene): the former two were generated from alternative splicing; the third one was encoded by a DNA region in the reverse complementary strand of the third intron of the XP gene. These findings strongly suggest that the genes of these cysteine-rich peptides were evolved by intron sliding, intron gain/loss, gene recombination and alternative splicing events in response to selective forces without changing their cysteine pattern. The evolution of these genes is dominated by intron sliding and intron loss.

The orthologue of the fruitfly sex behaviour gene fruitless in the mosquito Aedes aegypti: evolution of genomic organisation and alternative splicing

In Drosophila melanogaster the doublesex (dsx) and fruitless (fru) regulatory genes act at the bottom of the somatic sex determination pathway. Both are regulated via alternative splicing by an upstream female-specific TRA/TRA-2 complex, recognizing a common cis element. dsx controls somatic sexual differentiation of non-neural as well as of neural tissues. fru, on the other hand, expresses male-specific functions only in neural system where it is required to built the neural circuits underlying proper courtship behaviour. In the mosquito Aedes aegypti sex determination is different from Drosophila. The key male determiner M, which is located on one of a pair of homomorphic sex chromosomes, controls sex-specific splicing of the mosquito dsx orthologue. In this study we report the genomic organization and expression of the fru homologue in Ae. aegypti (Aeafru). We found that it is sex-specifically spliced suggesting that it is also under the control of the sex determination pathway. Comparative analyses between the Aeafru and Anopheles gambiae fru (Angfru) genomic loci revealed partial conservation of exon organization and extensive divergence of intron lengths. We find that Aeadsx and Aeafru share novel cis splicing regulatory elements conserved in the alternatively spliced regions. We propose that in Aedes aegypti sex-specific splicing of dsx and fru is most likely under the control of splicing regulatory factors which are different from TRA and TRA-2 found in other dipteran insects and discuss the potential use of fru and dsx for developing new genetic strategies in vector control.

Alternative splicing: a pivotal step between eukaryotic transcription and translation

Alternative splicing was discovered simultaneously with splicing over three decades ago. Since then, an enormous body of evidence has demonstrated the prevalence of alternative splicing in multicellular eukaryotes, its key roles in determining tissue- and species-specific differentiation patterns, the multiple post- and co-transcriptional regulatory mechanisms that control it, and its causal role in hereditary disease and cancer. The emerging evidence places alternative splicing in a central position in the flow of eukaryotic genetic information, between transcription and translation, in that it can respond not only to various signalling pathways that target the splicing machinery but also to transcription factors and chromatin structure.

Function of alternative splicing

Almost all polymerase II transcripts undergo alternative pre-mRNA splicing. Here, we review the functions of alternative splicing events that have been experimentally determined. The overall function of alternative splicing is to increase the diversity of mRNAs expressed from the genome. Alternative splicing changes proteins encoded by mRNAs, which has profound functional effects. Experimental analysis of these protein isoforms showed that alternative splicing regulates binding between proteins, between proteins and nucleic acids as well as between proteins and membranes. Alternative splicing regulates the localization of proteins, their enzymatic properties and their interaction with ligands. In most cases, changes caused by individual splicing isoforms are small. However, cells typically coordinate numerous changes in 'splicing programs', which can have strong effects on cell proliferation, cell survival and properties of the nervous system. Due to its widespread usage and molecular versatility, alternative splicing emerges as a central element in gene regulation that interferes with almost every biological function analyzed.

Direct competition between hnRNP C and U2AF65 protects the transcriptome from the exonization of Alu elements

There are ~650,000 Alu elements in transcribed regions of the human genome. These elements contain cryptic splice sites, so they are in constant danger of aberrant incorporation into mature transcripts. Despite posing a major threat to transcriptome integrity, little is known about the molecular mechanisms preventing their inclusion. Here, we present a mechanism for protecting the human transcriptome from the aberrant exonization of transposable elements. Quantitative iCLIP data show that the RNA-binding protein hnRNP C competes with the splicing factor U2AF65 at many genuine and cryptic splice sites. Loss of hnRNP C leads to formation of previously suppressed Alu exons, which severely disrupt transcript function. Minigene experiments explain disease-associated mutations in Alu elements that hamper hnRNP C binding. Thus, by preventing U2AF65 binding to Alu elements, hnRNP C plays a critical role as a genome-wide sentinel protecting the transcriptome. The findings have important implications for human evolution and disease.

Regulation of gene expression in mammalian nervous system through alternative pre-mRNA splicing coupled with RNA quality control mechanisms

Eukaryotic gene expression is orchestrated on a genome-wide scale through several post-transcriptional mechanisms. Of these, alternative pre-mRNA splicing expands the proteome diversity and modulates mRNA stability through downstream RNA quality control (QC) pathways including nonsense-mediated decay (NMD) of mRNAs containing premature termination codons and nuclear retention and elimination (NRE) of intron-containing transcripts. Although originally identified as mechanisms for eliminating aberrant transcripts, a growing body of evidence suggests that NMD and NRE coupled with deliberate changes in pre-mRNA splicing patterns are also used in a number of biological contexts for deterministic control of gene expression. Here we review recent studies elucidating molecular mechanisms and biological significance of these gene regulation strategies with a specific focus on their roles in nervous system development and physiology. This article is part of a Special Issue entitled 'RNA and splicing regulation in neurodegeneration'.

Pre-mRNA splicing is a determinant of nucleosome organization

Chromatin organization affects alternative splicing and previous studies have shown that exons have increased nucleosome occupancy compared with their flanking introns. To determine whether alternative splicing affects chromatin organization we developed a system in which the alternative splicing pattern switched from inclusion to skipping as a function of time. Changes in nucleosome occupancy were correlated with the change in the splicing pattern. Surprisingly, strengthening of the 5' splice site or strengthening the base pairing of U1 snRNA with an internal exon abrogated the skipping of the internal exons and also affected chromatin organization. Over-expression of splicing regulatory proteins also affected the splicing pattern and changed nucleosome occupancy. A specific splicing inhibitor was used to show that splicing impacts nucleosome organization endogenously. The effect of splicing on the chromatin required a functional U1 snRNA base pairing with the 5' splice site, but U1 pairing was not essential for U1 snRNA enhancement of transcription. Overall, these results suggest that splicing can affect chromatin organization.

The evolutionary landscape of alternative splicing in vertebrate species

How species with similar repertoires of protein-coding genes differ so markedly at the phenotypic level is poorly understood. By comparing organ transcriptomes from vertebrate species spanning ~350 million years of evolution, we observed significant differences in alternative splicing complexity between vertebrate lineages, with the highest complexity in primates. Within 6 million years, the splicing profiles of physiologically equivalent organs diverged such that they are more strongly related to the identity of a species than they are to organ type. Most vertebrate species-specific splicing patterns are cis-directed. However, a subset of pronounced splicing changes are predicted to remodel protein interactions involving trans-acting regulators. These events likely further contributed to the diversification of splicing and other transcriptomic changes that underlie phenotypic differences among vertebrate species.

Prominin-2 and Other Relatives of CD133

Several molecules related to prominin-1/CD133, which was first characterized as a marker of mouse neuroepithelial stem cells and human hematopoietic stem cells, have been identified in various species. In mammals, a second prominin gene, prominin-2, has been identified and characterized, whereas in nonmammalian species, up to three prominin genes are potentially expressed. The structural similarities between prominin-1 and prominin-2 are, to some extent, reflected by their biochemical properties; both proteins are selectively concentrated in specific plasma membrane subdomains that protrude into the extracellular space and are released in small extracellular membrane vesicles. In contrast to the apically confined prominin-1, prominin-2 is distributed in a nonpolarized apico-basolateral fashion in polarized epithelial cells and appears to be expressed in separate epithelial cells. Their distinctive localization in plasma membrane protrusions is a hallmark of prominins, validating the naming of the family after its first identified member. Insights into the distinctive and/or complementary roles of the two prominins may be obtained by analyzing the evolutionary history of these proteins and the characteristics of orthologs and paralogs in more distantly related species. In addition, the characterization of prominins may shed light on the still elusive function of CD133.

Quantification of stochastic noise of splicing and polyadenylation in Entamoeba histolytica

Alternative splicing and polyadenylation were observed pervasively in eukaryotic messenger RNAs. These alternative isoforms could either be consequences of physiological regulation or stochastic noise of RNA processing. To quantify the extent of stochastic noise in splicing and polyadenylation, we analyzed the alternative usage of splicing and polyadenylation sites in Entamoeba histolytica using RNA-Seq. First, we identified a large number of rarely spliced alternative junctions and then showed that the occurrence of these alternative splicing events is correlated with splicing site sequence, occurrence of constitutive splicing events and messenger RNA abundance. Our results implied the majority of these alternative splicing events are likely to be stochastic error of splicing machineries, and we estimated the corresponding error rates. Second, we observed extensive microheterogeneity of polyadenylation cleavage sites, and the extent of such microheterogeneity is correlated with the occurrence of constitutive cleavage events, suggesting most of such microheterogeneity is likely to be stochastic. Overall, we only observed a small fraction of alternative splicing and polyadenylation isoforms that are unlikely to be solely stochastic, implying the functional relevance of alternative splicing and polyadenylation in E. histolytica is limited. Lastly, we revised the gene models and annotated their 3′UTR in AmoebaDB, providing valuable resources to the community.

Mechanisms of Drosophila Dscam mutually exclusive splicing regulation

Alternative splicing of pre-mRNA is a major mechanism to increase protein diversity in higher eukaryotes. Dscam, the Drosophila homologue of human DSCAM (Down's syndrome cell adhesion molecule), generates up to 38016 isoforms through mutually exclusive splicing in four variable exon clusters. This enormous molecular diversity is functionally important for wiring of the nervous system and phagocytosis of invading pathogens. Current models explaining this complex splicing regulation include a default repressed state of the variable exon clusters to prevent the splicing together of adjacent exons, the presence of RNA secondary structures important for the release of one specific variable exon from the repressed state and combinatorial interaction of RNA-binding proteins for choosing a specific exon.

Caspase-1 and IL-1β processing in a teleost fish

Interleukine-1β (IL-1β) is the most studied pro-inflammatory cytokine, playing a central role in the generation of systemic and local responses to infection, injury, and immunological challenges. In mammals, IL-1β is synthesized as an inactive 31 kDa precursor that is cleaved by caspase-1 generating a 17.5 kDa secreted active mature form. The caspase-1 cleavage site strictly conserved in all mammalian IL-1β sequences is absent in IL-1β sequences reported for non-mammalian vertebrates. Recently, fish caspase-1 orthologues have been identified in sea bass (Dicentrarchus labrax) and sea bream (Sparus aurata) but very little is known regarding their processing and activity. In this work it is shown that sea bass caspase-1 auto-processing is similar to that of the human enzyme, resulting in active p24/p10 and p20/p10 heterodimers. Moreover, the presence of alternatively spliced variants of caspase-1 in sea bass is reported. The existence of caspase-1 isoforms in fish and in mammals suggests that they have been evolutionarily maintained and therefore are likely to play a regulatory role in the inflammatory response, as shown for other caspases. Finally, it is shown that sea bass and avian IL-1β are specifically cleaved by caspase-1 at different but phylogenetically conserved aspartates, distinct from the cleavage site of mammalian IL-1β.

Complex evolutionary relationships among four classes of modular RNA-binding splicing regulators in eukaryotes: the hnRNP, SR, ELAV-like and CELF proteins

Alternative RNA splicing in multicellular organisms is regulated by a large group of proteins of mainly unknown origin. To predict the functions of these proteins, classification of their domains at the sequence and structural level is necessary. We have focused on four groups of splicing regulators, the heterogeneous nuclear ribonucleoprotein (hnRNP), serine-arginine (SR), embryonic lethal, abnormal vision (ELAV)-like, and CUG-BP and ETR-like factor (CELF) proteins, that show increasing diversity among metazoa. Sequence and phylogenetic analyses were used to obtain a broader understanding of their evolutionary relationships. Surprisingly, when we characterised sequence similarities across full-length sequences and conserved domains of ten metazoan species, we found some hnRNPs were more closely related to SR, ELAV-like and CELF proteins than to other hnRNPs. Phylogenetic analyses and the distribution of the RRM domains suggest that these proteins diversified before the last common ancestor of the metazoans studied here through domain acquisition and duplication to create genes of mixed evolutionary origin. We propose that these proteins were derived independently rather than through the expansion of a single protein family. Our results highlight inconsistencies in the current classification system for these regulators, which does not adequately reflect their evolutionary relationships, and suggests that a domain-based classification scheme may have more utility.

Next-generation sequencing in the clinic: promises and challenges

The advent of next generation sequencing (NGS) technologies has revolutionized the field of genomics, enabling fast and cost-effective generation of genome-scale sequence data with exquisite resolution and accuracy. Over the past years, rapid technological advances led by academic institutions and companies have continued to broaden NGS applications from research to the clinic. A recent crop of discoveries have highlighted the medical impact of NGS technologies on Mendelian and complex diseases, particularly cancer. However, the ever-increasing pace of NGS adoption presents enormous challenges in terms of data processing, storage, management and interpretation as well as sequencing quality control, which hinder the translation from sequence data into clinical practice. In this review, we first summarize the technical characteristics and performance of current NGS platforms. We further highlight advances in the applications of NGS technologies towards the development of clinical diagnostics and therapeutics. Common issues in NGS workflows are also discussed to guide the selection of NGS platforms and pipelines for specific research purposes.

Whence genes in pieces: reconstruction of the exon-intron gene structures of the last eukaryotic common ancestor and other ancestral eukaryotes

In eukaryotes, protein-coding sequences are interrupted by non-coding sequences known as introns. During mRNA maturation, introns are excised by the spliceosome and the coding regions, exons, are spliced to form the mature coding region. The intron densities widely differ between eukaryotic lineages, from 6 to 7 introns per kb of coding sequence in vertebrates, some invertebrates and green plants, to only a few introns across the entire genome in many unicellular eukaryotes. Evolutionary reconstructions using maximum likelihood methods suggest intron-rich ancestors for each major group of eukaryotes. For the last common ancestor of animals, the highest intron density of all extant and extinct eukaryotes was inferred, at 120-130% of the human intron density. Furthermore, an intron density within 53-74% of the human values was inferred for the last eukaryotic common ancestor. Accordingly, evolution of eukaryotic genes in all lines of descent involved primarily intron loss, with substantial gain only at the bases of several branches including plants and animals. These conclusions have substantial biological implications indicating that the common ancestor of all modern eukaryotes was a complex organism with a gene architecture resembling those in multicellular organisms. Alternative splicing most likely initially appeared as an inevitable result of splicing errors and only later was employed to generate structural and functional diversification of proteins.