Alternative splicing and evolution: diversification, exon definition and function

Diversity of teleost leukocyte molecules: role of alternative splicing

Alternative splicing is an important mechanism of gene expression control that also produces a large proteome from a limited number of genes. In the immune system of mammals, numerous relevant genes have been found to undergo alternative splicing that contributes to the complexity of immune response. An increasing number of reports have recently indicated that alternative splicing also occurs in other vertebrates, such as fish. In this review we summarize the general features of such molecular events in cytokines and leukocyte co-receptors and their contribution to diversity and regulation of fish leukocytes.

Complex functions of Mef2 splice variants in the differentiation of endoderm and of a neuronal cell type in a sea anemone

In triploblastic animals, mesoderm gives rise to many tissues and organs, including muscle. By contrast, the representatives of the diploblastic phylum Cnidaria (corals, sea anemones, jellyfish and hydroids) lack mesoderm but possess muscle. In vertebrates and insects, the transcription factor Mef2 plays a pivotal role in muscle differentiation; however, it is also an important regulator of neuron differentiation and survival. In the sea anemone Nematostella vectensis, an organism that lacks mesoderm but has muscles and neurons, Mef2 (Nvmef2) has been reported in single ectodermal cells of likely neural origin. To our surprise, we found that Nvmef2 is alternatively spliced, forming differentially expressed variants. Using morpholino-mediated knockdown and mRNA injection, we demonstrate that specific splice variants of Nvmef2 are required for the proliferation and differentiation of endodermal cells and for the development of ectodermal nematocytes, a neuronal cell type. Moreover, we identified a small conserved motif in the transactivation domain that is crucially involved in the endodermal function of Nvmef2. The identification of a crucial and conserved motif in the transactivation domain predicts a similarly important role in vertebrate Mef2 function. This is the first functional study of a determinant of several mesodermal derivatives in a diploblastic animal. Our data suggest that the involvement of alternative splice variants of Mef2 in endomesoderm and neuron differentiation predates the cnidarian-bilaterian split.

InTRONs in biotech

Eukaryotic gene expression relies on several complex molecular machineries that act in a highly coordinated fashion. These machineries govern all the different steps of mRNA maturation, from gene transcription and pre-mRNA processing in the nucleus to the export of the mRNA to the cytoplasm and its translation. In particular, the pre-mRNA splicing process consists in the joining together of sequences (known as “exons”) that have to be differentiated from their intervening sequences commonly referred to as “introns.” The complex required to perform this process is a very dynamic macromolecular ribonucleoprotein assembly that functions as an enzyme, and is called the “spliceosome.” Because of its flexibility, the splicing process represents one of the main mechanisms of qualitative and quantitative regulation of gene expression in eukaryotic genomes. This flexibility is mainly due to the possibility of alternatively recognizing the various exons that are present in a pre-mRNA molecule and therefore enabling the possibility of obtaining multiple transcripts from the same gene. However, regulation of gene expression by the spliceosome is also achieved through its ability to influence many other gene expression steps that include transcription, mRNA export, mRNA stability, and even protein translation. Therefore, from a biotechnological point of view the splicing process can be exploited to improve production strategies and processes of molecules of interest. In this work, we have aimed to provide an overview on how biotechnology applications may benefit from the introduction of introns within a sequence of interest.

Changes in exon-intron structure during vertebrate evolution affect the splicing pattern of exons

Exon-intron architecture is one of the major features directing the splicing machinery to the short exons that are located within long flanking introns. However, the evolutionary dynamics of exon-intron architecture and its impact on splicing is largely unknown. Using a comparative genomic approach, we analyzed 17 vertebrate genomes and reconstructed the ancestral motifs of both 3' and 5' splice sites, as also the ancestral length of exons and introns. Our analyses suggest that vertebrate introns increased in length from the shortest ancestral introns to the longest primate introns. An evolutionary analysis of splice sites revealed that weak splice sites act as a restrictive force keeping introns short. In contrast, strong splice sites allow recognition of exons flanked by long introns. Reconstruction of the ancestral state suggests these phenomena were not prevalent in the vertebrate ancestor, but appeared during vertebrate evolution. By calculating evolutionary rate shifts in exons, we identified cis-acting regulatory sequences that became fixed during the transition from early vertebrates to mammals. Experimental validations performed on a selection of these hexamers confirmed their regulatory function. We additionally revealed many features of exons that can discriminate alternative from constitutive exons. These features were integrated into a machine-learning approach to predict whether an exon is alternative. Our algorithm obtains very high predictive power (AUC of 0.91), and using these predictions we have identified and successfully validated novel alternatively spliced exons. Overall, we provide novel insights regarding the evolutionary constraints acting upon exons and their recognition by the splicing machinery.

Approaches to Fungal Genome Annotation

Fungal genome annotation is the starting point for analysis of genome content. This generally involves the application of diverse methods to identify features on a genome assembly such as protein-coding and non-coding genes, repeats and transposable elements, and pseudogenes. Here we describe tools and methods leveraged for eukaryotic genome annotation with a focus on the annotation of fungal nuclear and mitochondrial genomes. We highlight the application of the latest technologies and tools to improve the quality of predicted gene sets. The Broad Institute eukaryotic genome annotation pipeline is described as one example of how such methods and tools are integrated into a sequencing center's production genome annotation environment.

Yeast pheromone receptor genes STE2 and STE3 are differently regulated at the transcription and polyadenylation level

The orderly expression of specific genes is the basis for cell differentiation. Saccharomyces cerevisiae has two haploid mating types, a and α cells, in which the mating-specific genes are differentially expressed. When a and α cells are committed to mate, their growth is arrested. Here we show that a cryptic polyadenylation site is present inside the coding region of the a-specific STE2 gene, encoding the receptor for the α-factor. The two cell types produce an incomplete STE2 transcript, but only a cells generate full-length STE2 mRNA. We eliminated the cryptic poly(A) signal, thereby allowing the production of a complete STE2 mRNA in α cells. We mutagenized α cells and isolated a mutant producing full-length STE2 mRNA. The mutation occurred in the ITC1 gene, whose product, together with the product of ISW2, is known to repress STE2 transcriptional initiation. We propose that the regulation of the yeast mating genes is achieved through a concerted mechanism involving transcriptional and posttranscriptional events. In particular, the early poly(A) site in STE2 could contribute to a complete shutoff of its expression in α cells, avoiding autocrine activation and growth arrest. Remarkably, no cryptic poly(A) sites are present in the a-factor receptor STE3 gene, indicating that S. cerevisiae has devised different strategies to regulate the two receptor genes. It is predictable that a correlation between the repression of a gene and the presence of a cryptic poly(A) site could also be found in other organisms, especially when expression of that gene may be harmful.

Big defensins, a diverse family of antimicrobial peptides that follows different patterns of expression in hemocytes of the oyster Crassostrea gigas

Background

Big defensin is an antimicrobial peptide composed of a highly hydrophobic N-terminal region and a cationic C-terminal region containing six cysteine residues involved in three internal disulfide bridges. While big defensin sequences have been reported in various mollusk species, few studies have been devoted to their sequence diversity, gene organization and their expression in response to microbial infections.

Findings

Using the high-throughput Digital Gene Expression approach, we have identified in Crassostrea gigas oysters several sequences coding for big defensins induced in response to a Vibrio infection. We showed that the oyster big defensin family is composed of three members (named Cg-BigDef1, Cg-BigDef2 and Cg-BigDef3) that are encoded by distinct genomic sequences. All Cg-BigDefs contain a hydrophobic N-terminal domain and a cationic C-terminal domain that resembles vertebrate β-defensins. Both domains are encoded by separate exons. We found that big defensins form a group predominantly present in mollusks and closer to vertebrate defensins than to invertebrate and fungi CSαβ-containing defensins. Moreover, we showed that Cg-BigDefs are expressed in oyster hemocytes only and follow different patterns of gene expression. While Cg-BigDef3 is non-regulated, both Cg-BigDef1 and Cg-BigDef2 transcripts are strongly induced in response to bacterial challenge. Induction was dependent on pathogen associated molecular patterns but not damage-dependent. The inducibility of Cg-BigDef1 was confirmed by HPLC and mass spectrometry, since ions with a molecular mass compatible with mature Cg-BigDef1 (10.7 kDa) were present in immune-challenged oysters only. From our biochemical data, native Cg-BigDef1 would result from the elimination of a prepropeptide sequence and the cyclization of the resulting N-terminal glutamine residue into a pyroglutamic acid.

Conclusions

We provide here the first report showing that big defensins form a family of antimicrobial peptides diverse not only in terms of sequences but also in terms of genomic organization and regulation of gene expression.

Expanding the action of duplex RNAs into the nucleus: redirecting alternative splicing

Double-stranded RNAs are powerful agents for silencing gene expression in the cytoplasm of mammalian cells. The potential for duplex RNAs to control expression in the nucleus has received less attention. Here, we investigate the ability of small RNAs to redirect splicing. We identify RNAs targeting an aberrant splice site that restore splicing and production of functional protein. RNAs can target sequences within exons or introns and affect the inclusion of exons within SMN2 and dystrophin, genes responsible for spinal muscular atrophy and Duchenne muscular dystrophy, respectively. Duplex RNAs recruit argonaute 2 (AGO2) to pre-mRNA transcripts and altered splicing requires AGO2 expression. AGO2 promotes transcript cleavage in the cytoplasm, but recruitment of AGO2 to pre-mRNAs does not reduce transcript levels, exposing a difference between cytoplasmic and nuclear pathways. Involvement of AGO2 in splicing, a classical nuclear process, reinforces the conclusion from studies of RNA-mediated transcriptional silencing that RNAi pathways can be adapted to function in the mammalian nucleus. These data provide a new strategy for controlling splicing and expand the reach of small RNAs within the nucleus of mammalian cells.

Functional consequences of developmentally regulated alternative splicing

Genome-wide analyses of metazoan transcriptomes have revealed an unexpected level of mRNA diversity that is generated by alternative splicing. Recently, regulatory networks have been identified through which splicing promotes dynamic remodelling of the transcriptome to promote physiological changes, which involve robust and coordinated alternative splicing transitions. The regulation of splicing in yeast, worms, flies and vertebrates affects a variety of biological processes. The functional classes of genes that are regulated by alternative splicing include both those with widespread homeostatic activities and those with cell-type-specific functions. Alternative splicing can drive determinative physiological change or can have a permissive role by providing mRNA variability that is used by other regulatory mechanisms.

A detailed history of intron-rich eukaryotic ancestors inferred from a global survey of 100 complete genomes

Protein-coding genes in eukaryotes are interrupted by introns, but intron densities widely differ between eukaryotic lineages. Vertebrates, some invertebrates and green plants have intron-rich genes, with 6–7 introns per kilobase of coding sequence, whereas most of the other eukaryotes have intron-poor genes. We reconstructed the history of intron gain and loss using a probabilistic Markov model (Markov Chain Monte Carlo, MCMC) on 245 orthologous genes from 99 genomes representing the three of the five supergroups of eukaryotes for which multiple genome sequences are available. Intron-rich ancestors are confidently reconstructed for each major group, with 53 to 74% of the human intron density inferred with 95% confidence for the Last Eukaryotic Common Ancestor (LECA). The results of the MCMC reconstruction are compared with the reconstructions obtained using Maximum Likelihood (ML) and Dollo parsimony methods. An excellent agreement between the MCMC and ML inferences is demonstrated whereas Dollo parsimony introduces a noticeable bias in the estimations, typically yielding lower ancestral intron densities than MCMC and ML. Evolution of eukaryotic genes was dominated by intron loss, with substantial gain only at the bases of several major branches including plants and animals. The highest intron density, 120 to 130% of the human value, is inferred for the last common ancestor of animals. The reconstruction shows that the entire line of descent from LECA to mammals was intron-rich, a state conducive to the evolution of alternative splicing.

In eukaryotes, protein-coding genes are interrupted by non-coding introns. The intron densities widely differ, from 6–7 introns per kilobase of coding sequence in vertebrates, some invertebrates and plants, to only a few introns across the entire genome in many unicellular forms. We applied a robust statistical methodology, Markov Chain Monte Carlo, to reconstruct the history of intron gain and loss throughout the evolution of eukaryotes using a set of 245 homologous genes from 99 genomes that represent the diversity of eukaryotes. Intron-rich ancestors were confidently inferred for each major eukaryotic group including 53% to 74% of the human intron density for the last eukaryotic common ancestor, and 120% to 130% of the human value for the last common ancestor of animals. Evolution of eukaryotic genes involved primarily intron loss, with substantial gain only at the bases of several major branches including plants and animals. Thus, the common ancestor of all extant eukaryotes was a complex organism with a gene architecture resembling those in multicellular organisms. The line of descent from the last common ancestor to mammals was an uninterrupted intron-rich state that, given the error-prone splicing in intron-rich organisms, was conducive to the elaboration of functional alternative splicing.

Pol III binding in six mammals shows conservation among amino acid isotypes despite divergence among tRNA genes

RNA polymerase III (pol III) transcription of transfer RNA (tRNA) genes is essential for generating the tRNA adapter molecules that link genetic sequence and protein translation. By mapping pol III occupancy genome-wide in the livers of mouse, rat, human, macaque, dog and opossum, we found that pol III binding to individual tRNA genes varies substantially in strength and location. However, taking into account tRNA redundancies by grouping pol III occupancy into 46 anticodon isoacceptor families or 21 amino acid-based isotype classes shows strong conservation. Similarly, pol III occupancy of amino-acid isotypes is almost invariant among transcriptionally and evolutionarily diverse tissues in mouse. Thus, synthesis of functional tRNA isotypes has been highly constrained, though the usage of individual tRNA genes has evolved rapidly.

A G-to-T transversion at the splice acceptor site of dystrophin exon 14 shows multiple splicing outcomes that are not exemplified by transition mutations

Mutations at splicing consensus sequences have been shown to induce splicing errors such as exon skipping or cryptic splice site activation. Here, we identified eight splicing products caused by a G-to-T transversion mutation at the splice acceptor site of exon 14 of the dystrophin gene (c.1603-1G>T). Unexpectedly, the most abundant product showed skipping of the two consecutive exons 14 and 15, and exon 14 skipping was observed as the second most abundant product. To examine the cause of this splicing multiplicity, minigenes containing dystrophin exons 14 and 15 with their flanking introns were constructed and subjected to in vitro splicing. Minigenes with the wild-type sequence or a G>A transition at position c.1603-1 produced only the mature mRNA. On the other hand, the minigenes with a G>T or G>C transversion mutation produced multiple splicing products. A time-course analysis of the in vitro splicing revealed that splicing of the middle intron, intron 14, was the first step in transcript maturation for all four minigene constructs. The identity of the mutant nucleotide, but not its position, is a factor leading to multiple splicing outcomes. Our results suggest that exon skipping therapy for Duchenne's muscular dystrophy should be carefully monitored for their splicing outcomes.

Alternative splicing of T-box transcription factor genes

T-box (TBX) transcription factors are an ancient gene family with critical roles in embryogenesis. Currently, TBX3, TBX5, and TBX20 are TBX genes defined to have multiple protein isoforms created by alternative splicing and characterized by expression and functional studies. These proteins are important for development as mutations lead to severe developmental disorders in humans and mice. Cumulative studies suggest that alternative splicing of these genes can regulate TBX activities during multiple biological processes including cardiogenesis, limb development, and cancer mechanisms. This mini-review focuses on how alternative splicing adds complexity to transcriptional regulation of target genes controlled by TBX transcription factors.

A limited role for gene duplications in the evolution of platypus venom

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the "venome" of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

Independent effects of alternative splicing and structural constraint on the evolution of mammalian coding exons

Alternative splicing (AS) is known to significantly affect exon-level protein evolutionary rates in mammals. Particularly, alternatively spliced exons (ASEs) have a higher nonsynonymous-to-synonymous substitution rate (dN/dS) ratio than constitutively spliced exons (CSEs), possibly because the former are required only occasionally for normal biological functions. Meanwhile, intrinsically disordered regions (IDRs), the protein regions lacking fixed 3D structures, are also reported to have an increased evolutionary rate due to lack of structural constraint. Interestingly, IDRs tend to be located in alternative protein regions. Yet which of these two factors is the major determinant of the increased dN/dS in mammalian ASEs remains unclear. By comparing human-macaque and human-mouse one-to-one orthologous genes, we demonstrate that AS and protein structural disorder have independent effects on mammalian exon evolution. We performed analyses of covariance to demonstrate that the slopes of the (dN/dS-percentage of IDR) regression lines differ significantly between CSEs and ASEs. In other words, the dN/dS ratios of both ASEs and CSEs increase with the proportion of IDR (PIDR), whereas ASEs have higher dN/dS ratios than CSEs when they have similar PIDRs. Since ASEs and IDRs may less frequently overlap with protein domains (which also affect dN/dS), we also examined the correlations between dN/dS ratio and exon type/PIDR by controlling for the density of protein domain. We found that the effects of exon type and PIDR on dN/dS are both independent of domain density. Our results imply that nature can select for different biological features with regard to ASEs and IDRs, even though the two biological features tend to be localized in the same protein regions.

Key features of the two-intron Saccharomyces cerevisiae gene SUS1 contribute to its alternative splicing

Alternative pre-mRNA splicing allows dramatic expansion of the eukaryotic proteome and facilitates cellular response to changes in environmental conditions. The Saccharomyces cerevisiae gene SUS1, which encodes a protein involved in mRNA export and histone H2B deubiquitination, contains two introns; non-canonical sequences in the first intron contribute to its retention, a common form of alternative splicing in plants and fungi. Here we show that the pattern of SUS1 splicing changes in response to environmental change such as temperature elevation, and the retained intron product is subject to nonsense-mediated decay. The activities of different splicing factors determine the pattern of SUS1 splicing, including intron retention and exon skipping. Unexpectedly, removal of the 3′ intron is affected by splicing of the upstream intron, suggesting that cross-exon interactions influence intron removal. Production of different SUS1 isoforms is important for cellular function, as we find that the temperature sensitivity and histone H2B deubiquitination defects observed in sus1Δ cells are only partially suppressed by SUS1 cDNA, but SUS1 that is able to undergo splicing complements these phenotypes. These data illustrate a role for S. cerevisiae alternative splicing in histone modification and cellular function and reveal important mechanisms for splicing of yeast genes containing multiple introns.

Loss of exon identity is a common mechanism of human inherited disease

It is widely accepted that at least 10% of all mutations causing human inherited disease disrupt splice-site consensus sequences. In contrast to splice-site mutations, the role of auxiliary cis-acting elements such as exonic splicing enhancers (ESE) and exonic splicing silencers (ESS) in human inherited disease is still poorly understood. Here we use a top-down approach to determine rates of loss or gain of known human exonic splicing regulatory (ESR) sequences associated with either disease-causing mutations or putatively neutral single nucleotide polymorphisms (SNPs). We observe significant enrichment toward loss of ESEs and gain of ESSs among inherited disease-causing variants relative to neutral polymorphisms, indicating that exon skipping may play a prominent role in aberrant gene regulation. Both computational and biochemical approaches underscore the relevance of exonic splicing enhancer loss and silencer gain in inherited disease. Additionally, we provide direct evidence that both SRp20 (SRSF3) and possibly PTB (PTBP1) are involved in the function of a splicing silencer that is created de novo by a total of 83 different inherited disease mutations in 67 different disease genes. Taken together, we find that ~25% (7154/27,681) of known mis-sense and nonsense disease-causing mutations alter functional splicing signals within exons, suggesting a much more widespread role for aberrant mRNA processing in causing human inherited disease than has hitherto been appreciated.

Alternative splicing in spinal muscular atrophy underscores the role of an intron definition model

Humans have two nearly identical copies of the Survival Motor Neuron (SMN) gene: SMN1 and SMN2. The two SMN genes code for identical proteins; however, SMN2 predominantly generates a shorter transcript due to skipping of exon 7, the last coding exon. Skipping of SMN2 exon 7 leads to production of a truncated SMN protein that is highly unstable. The inability of SMN2 to compensate for the loss of SMN1 results in spinal muscular atrophy (SMA), the second most prevalent genetic cause of infant mortality. Since SMN2 is almost universally present in SMA patients, correction of SMN2 exon 7 splicing holds the promise for cure. Consistently, SMN2 exon 7 splicing has emerged as one of the best studied splicing systems in humans. The vast amount of recent literature provides a clue that SMN2 exon 7 splicing is regulated by an intron definition mechanism, which does not require cross-exon communication as prerequisite for exon inclusion. Our conclusion is based on the prominent role of intronic cis-elements, some of them have emerged as the frontrunners among potential therapeutic targets of SMA. Further, the widely expressed T-cell-restricted intracellular antigen-1 (TIA1), a member of the Q-rich domain containing RNA-binding proteins, has recently been found to regulate SMN exon 7 splicing by binding to intron 7 sequences away from the 5′ ss. These findings make a strong argument for an "intron definition model", according to which regulatory sequences within a downstream intron are capable of enforcing exon inclusion even in the absence of a defined upstream 3′ ss of an alternatively spliced exon.

Pervasive haplotypic variation in the spliceo-transcriptome of the human major histocompatibility complex

The human major histocompatibility complex (MHC) on chromosome 6p21 is a paradigm for genomics, showing remarkable polymorphism and striking association with immune and non-immune diseases. The complex genomic landscape of the MHC, notably strong linkage disequilibrium, has made resolving causal variants very challenging. A promising approach is to investigate gene expression levels considered as tractable intermediate phenotypes in mapping complex diseases. However, how transcription varies across the MHC, notably relative to specific haplotypes, remains unknown. Here, using an original hybrid tiling and splice junction microarray that includes alternate allele probes, we draw the first high-resolution strand-specific transcription map for three common MHC haplotypes (HLA-A1-B8-Cw7-DR3, HLA-A3-B7-Cw7-DR15, and HLA-A26-B18-Cw5-DR3-DQ2) strongly associated with autoimmune diseases including type 1 diabetes, systemic lupus erythematosus, and multiple sclerosis. We find that haplotype-specific differences in gene expression are common across the MHC, affecting 96 genes (46.4%), most significantly the zing finger protein gene ZFP57. Differentially expressed probes are correlated with polymorphisms between haplotypes, consistent with cis effects that we directly demonstrate for ZFP57 in a cohort of healthy volunteers (P = 1.2 × 10(-14)). We establish that alternative splicing is significantly more frequent in the MHC than genome-wide (72.5% vs. 62.1% of genes, P ≤ 1 × 10(-4)) and shows marked haplotypic differences. We also unmask novel and abundant intergenic transcription involving 31% of transcribed blocks identified. Our study reveals that the renowned MHC polymorphism also manifests as transcript diversity, and our novel haplotype-based approach marks a new step toward identification of regulatory variants involved in the control of MHC-associated phenotypes and diseases.

Novel α1 and γ2 GABAA receptor subunit mutations in families with idiopathic generalized epilepsy

Epilepsy is a heterogeneous neurological disease affecting approximately 50 million people worldwide. Genetic factors play an important role in both the onset and severity of the condition, with mutations in several ion-channel genes being implicated, including those encoding the GABA(A) receptor. Here, we evaluated the frequency of additional mutations in the GABA(A) receptor by direct sequencing of the complete open reading frame of the GABRA1 and GABRG2 genes from a cohort of French Canadian families with idiopathic generalized epilepsy (IGE). Using this approach, we have identified three novel mutations that were absent in over 400 control chromosomes. In GABRA1, two mutations were found, with the first being a 25-bp insertion that was associated with intron retention (i.e. K353delins18X) and the second corresponding to a single point mutation that replaced the aspartate 219 residue with an asparagine (i.e. D219N). Electrophysiological analysis revealed that K353delins18X and D219N altered GABA(A) receptor function by reducing the total surface expression of mature protein and/or by curtailing neurotransmitter effectiveness. Both defects would be expected to have a detrimental effect on inhibitory control of neuronal circuits. In contrast, the single point mutation identified in the GABRG2 gene, namely P83S, was indistinguishable from the wildtype subunit in terms of surface expression and functionality. This finding was all the more intriguing as the mutation exhibited a high degree of penetrance in three generations of one French Canadian family. Further experimentation will be required to understand how this mutation contributes to the occurrence of IGE in these individuals.

Predicting mutually exclusive spliced exons based on exon length, splice site and reading frame conservation, and exon sequence homology

BACKGROUND

Alternative splicing of pre-mature RNA is an important process eukaryotes utilize to increase their repertoire of different protein products. Several types of different alternative splice forms exist including exon skipping, differential splicing of exons at their 3'- or 5'-end, intron retention, and mutually exclusive splicing. The latter term is used for clusters of internal exons that are spliced in a mutually exclusive manner.

RESULTS

We have implemented an extension to the WebScipio software to search for mutually exclusive exons. Here, the search is based on the precondition that mutually exclusive exons encode regions of the same structural part of the protein product. This precondition provides restrictions to the search for candidate exons concerning their length, splice site conservation and reading frame preservation, and overall homology. Mutually exclusive exons that are not homologous and not of about the same length will not be found. Using the new algorithm, mutually exclusive exons in several example genes, a dynein heavy chain, a muscle myosin heavy chain, and Dscam were correctly identified. In addition, the algorithm was applied to the whole Drosophila melanogaster X chromosome and the results were compared to the Flybase annotation and an ab initio prediction. Clusters of mutually exclusive exons might be subsequent to each other and might encode dozens of exons.

CONCLUSIONS

This is the first implementation of an automatic search for mutually exclusive exons in eukaryotes. Exons are predicted and reconstructed in the same run providing the complete gene structure for the protein query of interest. WebScipio offers high quality gene structure figures with the clusters of mutually exclusive exons colour-coded, and several analysis tools for further manual inspection. The genome scale analysis of all genes of the Drosophila melanogaster X chromosome showed that WebScipio is able to find all but two of the 28 annotated mutually exclusive spliced exons and predicts 39 new candidate exons. Thus, WebScipio should be able to identify mutually exclusive spliced exons in any query sequence from any species with a very high probability. WebScipio is freely available to academics at http://www.webscipio.org.

SNP discovery in black cottonwood (Populus trichocarpa) by population transcriptome resequencing

The western black cottonwood (Populus trichocarpa) was the first tree to have its genome fully sequenced and has emerged as the model species for the study of secondary growth and wood formation. It is also a good candidate species for the production of lignocellulosic biofuels. Here, we present and make available to the research community the results of the sequencing of the transcriptome of developing xylem in 20 accessions with high-throughput next generation sequencing technology. We found over 0.5 million putative single nucleotide polymorphisms (SNPs) in 26,595 genes that are expressed in developing secondary xylem. More than two-thirds of all SNPs were found in annotated exons, with 18% and 14% in regions of the genome annotated as introns and intergenic, respectively, where only 3% and 4% of sequence reads mapped. This suggests that the current annotation of the poplar genome is remarkably incomplete and that there are many transcripts and novel genes waiting to be annotated. We hope that this resource will stimulate further research in expression profiling, detection of alternative splicing and adaptive evolution in poplar.

Autoregulatory circuit of human rpL3 expression requires hnRNP H1, NPM and KHSRP

Alternative pre-mRNA splicing (AS) is a major mechanism that allows proteomic variability in eukaryotic cells. However, many AS events result in mRNAs containing a premature termination codon, which are degraded by nonsense-mediated mRNA decay (NMD) pathway. We have previously demonstrated that human rpL3 autoregulates its expression through the association of AS with NMD. In fact, overexpression of rpL3 promotes downregulation of canonical splicing and upregulation of alternative splicing that produces an NMD-targeted mRNA isoform. The result of these events is a decreased production of rpL3. We have also identified heterogeneous nuclear ribonucleoprotein (hnRNP) H1 as a splicing factor involved in the regulation of rpL3 alternative splicing and identified its regulatory cis-elements within intron 3 transcript. Here, we report that NPM and KHSRP are two newly identified proteins involved in the regulation of rpL3 gene expression via AS-NMD. We demonstrate that hnRNP H1, KHSRP and NPM can be found associated, and present also in ribonucleoproteins (RNPs) including rpL3 and intron 3 RNA in vivo, and describe protein–protein and RNA–protein interactions. Moreover, our data provide an insight on the crucial role of hnRNP H1 in the regulation of the alternative splicing of the rpL3 gene.

SON controls cell-cycle progression by coordinated regulation of RNA splicing

It has been suspected that cell-cycle progression might be functionally coupled with RNA processing. However, little is known about the role of the precise splicing control in cell-cycle progression. Here, we report that SON, a large Ser/Arg (SR)-related protein, is a splicing cofactor contributing to efficient splicing of cell-cycle regulators. Downregulation of SON leads to severe impairment of spindle pole separation, microtubule dynamics, and genome integrity. These molecular defects result from inadequate RNA splicing of a specific set of cell-cycle-related genes that possess weak splice sites. Furthermore, we show that SON facilitates the interaction of SR proteins with RNA polymerase II and other key spliceosome components, suggesting its function in efficient cotranscriptional RNA processing. These results reveal a mechanism for controlling cell-cycle progression through SON-dependent constitutive splicing at suboptimal splice sites, with strong implications for its role in cancer and other human diseases.

Evolutionary connections between coding and splicing regulatory regions in the fibronectin EDA exon

Research on exonic coding sequences has demonstrated that many substitutions at the amino acid level may also reflect profound changes at the level of splicing regulatory regions. These results have revealed that, for many alternatively spliced exons, there is considerable pressure to strike a balance between two different and sometimes conflicting forces: the drive to improve the quality and production efficiency of proteins and the maintenance of proper exon recognition by the splicing machinery. Up to now, the systems used to investigate these connections have mostly focused on short alternatively spliced exons that contain a high density of splicing regulatory elements. Although this is obviously a desirable feature in order to maximize the chances of spotting connections, it also complicates the process of drawing straightforward evolutionary pathways between different species (because of the numerous alternative pathways through which the same end point can be achieved). The alternatively spliced fibronectin extra domain A exon (also referred to as EDI or EIIIA) does not have these limitations, as its inclusion is already known to depend on a single exonic splicing enhancer element within its sequence. In this study, we have compared the rat and human fibronectin EDA exons with regard to RNA structure, exonic splicing enhancer strengths, and SR protein occupancy. The results gained from these analyses have then been used to perform an accurate evaluation of EDA sequences observed in a wide range of animal species. This comparison strongly suggests the existence of an evolutionary connection between changes at the nucleotide levels and the need to maintain efficient EDA recognition in different species.

Role of the ubiquitin-like protein Hub1 in splice-site usage and alternative splicing

Alternative splicing of pre-messenger RNAs diversifies gene products in eukaryotes and is guided by factors that enable spliceosomes to recognize particular splice sites. Here we report that alternative splicing of Saccharomyces cerevisiae SRC1 pre-mRNA is promoted by the conserved ubiquitin-like protein Hub1. Structural and biochemical data show that Hub1 binds non-covalently to a conserved element termed HIND, which is present in the spliceosomal protein Snu66 in yeast and mammals, and Prp38 in plants. Hub1 binding mildly alters spliceosomal protein interactions and barely affects general splicing in S. cerevisiae. However, spliceosomes that lack Hub1, or are defective in Hub1-HIND interaction, cannot use certain non-canonical 5' splice sites and are defective in alternative SRC1 splicing. Hub1 confers alternative splicing not only when bound to HIND, but also when experimentally fused to Snu66, Prp38, or even the core splicing factor Prp8. Our study indicates a novel mechanism for splice site utilization that is guided by non-covalent modification of the spliceosome by an unconventional ubiquitin-like modifier.

Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition

Maternally deposited mRNAs direct early development before the initiation of zygotic transcription during mid-blastula transition (MBT). To study mechanisms regulating this developmental event in zebrafish, we applied mRNA deep sequencing technology and generated comprehensive information and valuable resources on transcriptome dynamics during early embryonic (egg to early gastrulation) stages. Genome-wide transcriptome analysis documented at least 8000 maternal genes and identified the earliest cohort of zygotic transcripts. We determined expression levels of maternal and zygotic transcripts with the highest resolution possible using mRNA-seq and clustered them based on their expression pattern. We unravel delayed polyadenylation in a large cohort of maternal transcripts prior to the MBT for the first time in zebrafish. Blocking polyadenylation of these transcripts confirms their role in regulating development from the MBT onward. Our study also identified a large number of novel transcribed regions in annotated and unannotated regions of the genome, which will facilitate reannotation of the zebrafish genome. We also identified splice variants with an estimated frequency of 50%-60%. Taken together, our data constitute a useful genomic information and valuable transcriptome resource for gene discovery and for understanding the mechanisms of early embryogenesis in zebrafish.

Networking in a global world: establishing functional connections between neural splicing regulators and their target transcripts

Recent genome-wide analyses have indicated that almost all primary transcripts from multi-exon human genes undergo alternative pre-mRNA splicing (AS). Given the prevalence of AS and its importance in expanding proteomic complexity, a major challenge that lies ahead is to determine the functional specificity of isoforms in a cellular context. A significant fraction of alternatively spliced transcripts are regulated in a tissue- or cell-type-specific manner, suggesting that these mRNA variants likely function in the generation of cellular diversity. Complementary to these observations, several tissue-specific splicing factors have been identified, and a number of methodological advances have enabled the identification of large repertoires of target transcripts regulated by these proteins. An emerging theme is that tissue-specific splicing factors regulate coherent sets of splice variants in genes known to function in related biological pathways. This review focuses on the recent progress in our understanding of neural-specific splicing factors and their regulatory networks and outlines existing and emerging strategies for uncovering important biological roles for the isoforms that comprise these networks.

Massive-scale RNA-Seq analysis of non ribosomal transcriptome in human trisomy 21

Hybridization- and tag-based technologies have been successfully used in Down syndrome to identify genes involved in various aspects of the pathogenesis. However, these technologies suffer from several limits and drawbacks and, to date, information about rare, even though relevant, RNA species such as long and small non-coding RNAs, is completely missing. Indeed, none of published works has still described the whole transcriptional landscape of Down syndrome. Although the recent advances in high-throughput RNA sequencing have revealed the complexity of transcriptomes, most of them rely on polyA enrichment protocols, able to detect only a small fraction of total RNA content. On the opposite end, massive-scale RNA sequencing on rRNA-depleted samples allows the survey of the complete set of coding and non-coding RNA species, now emerging as novel contributors to pathogenic mechanisms. Hence, in this work we analysed for the first time the complete transcriptome of human trisomic endothelial progenitor cells to an unprecedented level of resolution and sensitivity by RNA-sequencing. Our analysis allowed us to detect differential expression of even low expressed genes crucial for the pathogenesis, to disclose novel regions of active transcription outside yet annotated loci, and to investigate a plethora of non-polyadenylated long as well as short non coding RNAs. Novel splice isoforms for a large subset of crucial genes, and novel extended untranslated regions for known genes--possibly novel miRNA targets or regulatory sites for gene transcription--were also identified in this study. Coupling the rRNA depletion of samples, followed by high-throughput RNA-sequencing, to the easy availability of these cells renders this approach very feasible for transcriptome studies, offering the possibility of investigating in-depth blood-related pathological features of Down syndrome, as well as other genetic disorders.

Gene structure and splicing in schistosomes

Schistosomes are blood dwelling platyhelminths with a complex life cycle and persist in the definitive host during decades, indicating that they are very successful parasites. The challenge of infecting two hosts from different evolutionary branches, namely an invertebrate snail and a vertebrate, suggests that Schistosomes must display a very sophisticated genetic program to circumvent all the barriers imposed by the hosts' immune systems. Recent large-scale genome and transcriptome data from Schistosomes are facilitating the analysis of gene structure and splicing. Studying the structure of genes coding for secreted proteins is of particular interest since these proteins mediate processes in the host-parasite interface. Recent description of Micro-Exon Genes (MEGs), polymorphic mucin genes (SmPoMucs) and venom allergen-like (SmVALs) proteins with unusual gene structure apparently oriented towards generation of protein variability through alternate splicing, and the presence of multiple copies of these genes, indicate that the parasite developed a sophisticated system to interact with its hosts. This opens up opportunities for further studies with the use of proteomic techniques to better characterize the protein variability created by these systems and their role in parasite survival. Complete description of the functions of these variable proteins will greatly contribute to our understanding of host-parasite interactions.

Alternative splicing and nonsense-mediated RNA decay contribute to the regulation of SHOX expression

The human SHOX gene is composed of seven exons and encodes a paired-related homeodomain transcription factor. SHOX mutations or deletions have been associated with different short stature syndromes implying a role in growth and bone formation. During development, SHOX is expressed in a highly specific spatiotemporal expression pattern, the underlying regulatory mechanisms of which remain largely unknown. We have analysed SHOX expression in diverse embryonic, fetal and adult human tissues and detected expression in many tissues that were not known to express SHOX before, e.g. distinct brain regions. By using RT-PCR and comparing the results with RNA-Seq data, we have identified four novel exons (exon 2a, 7-1, 7-2 and 7-3) contributing to different SHOX isoforms, and also established an expression profile for the emerging new SHOX isoforms. Interestingly, we found the exon 7 variants to be exclusively expressed in fetal neural tissues, which could argue for a specific role of these variants during brain development. A bioinformatical analysis of the three novel 3′UTR exons yielded insights into the putative role of the different 3′UTRs as targets for miRNA binding. Functional analysis revealed that inclusion of exon 2a leads to nonsense-mediated RNA decay altering SHOX expression in a tissue and time specific manner. In conclusion, SHOX expression is regulated by different mechanisms and alternative splicing coupled with nonsense-mediated RNA decay constitutes a further component that can be used to fine tune the SHOX expression level.

Does constructive neutral evolution play an important role in the origin of cellular complexity? Making sense of the origins and uses of biological complexity

Recently, constructive neutral evolution has been touted as an important concept for the understanding of the emergence of cellular complexity. It has been invoked to help explain the development and retention of, amongst others, RNA splicing, RNA editing and ribosomal and mitochondrial respiratory chain complexity. The theory originated as a welcome explanation of isolated small scale cellular idiosyncrasies and as a reaction to 'overselectionism'. Here I contend, that in its extended form, it has major conceptual problems, can not explain observed patterns of complex processes, is too easily dismissive of alternative selectionist models, underestimates the creative force of complexity as such, and--if seen as a major evolutionary mechanism for all organisms--could stifle further thought regarding the evolution of highly complex biological processes.

A novel soluble immune-type receptor (SITR) in teleost fish: carp SITR is involved in the nitric oxide-mediated response to a protozoan parasite

Background

The innate immune system relies upon a wide range of germ-line encoded receptors including a large number of immunoglobulin superfamily (IgSF) receptors. Different Ig-like immune receptor families have been reported in mammals, birds, amphibians and fish. Most innate immune receptors of the IgSF are type I transmembrane proteins containing one or more extracellular Ig-like domains and their regulation of effector functions is mediated intracellularly by distinct stimulatory or inhibitory pathways.

Methodology/Principal Findings

Carp SITR was found in a substracted cDNA repertoire from carp macrophages, enriched for genes up-regulated in response to the protozoan parasite Trypanoplasma borreli. Carp SITR is a type I protein with two extracellular Ig domains in a unique organisation of a N-proximal V/C2 (or I-) type and a C-proximal V-type Ig domain, devoid of a transmembrane domain or any intracytoplasmic signalling motif. The carp SITR C-proximal V-type Ig domain, in particular, has a close sequence similarity and conserved structural characteristics to the mammalian CD300 molecules. By generating an anti-SITR antibody we could show that SITR protein expression was restricted to cells of the myeloid lineage. Carp SITR is abundantly expressed in macrophages and is secreted upon in vitro stimulation with the protozoan parasite T. borreli. Secretion of SITR protein during in vivo T. borreli infection suggests a role for this IgSF receptor in the host response to this protozoan parasite. Overexpression of carp SITR in mouse macrophages and knock-down of SITR protein expression in carp macrophages, using morpholino antisense technology, provided evidence for the involvement of carp SITR in the parasite-induced NO production.

Conclusion/Significance

We report the structural and functional characterization of a novel soluble immune-type receptor (SITR) in a teleost fish and propose a role for carp SITR in the NO-mediated response to a protozoan parasite.

Developments in RNA splicing and disease

Pre-mRNA processing, including 5'-end capping, splicing, editing, and polyadenylation, consists of a series of orchestrated and primarily cotranscriptional steps that ensure both the high fidelity and extreme diversity characteristic of eukaryotic gene expression. Alternative splicing and editing allow relatively small genomes to encode vast proteomic arrays while alternative 3'-end formation enables variations in mRNA localization, translation, and stability. Of course, this mechanistic complexity comes at a high price. Mutations in the myriad of RNA sequence elements that regulate mRNA biogenesis, as well as the trans-acting factors that act upon these sequences, underlie a number of human diseases. In this review, we focus on one of these key RNA processing steps, splicing, to highlight recent studies that describe both conventional and novel pathogenic mechanisms that underlie muscle and neurological diseases.

Multiple Roles of Alu-Related Noncoding RNAs

Repetitive Alu and Alu-related elements are present in primates, tree shrews (Scandentia), and rodents and have expanded to 1.3 million copies in the human genome by nonautonomous retrotransposition. Pol III transcription from these elements occurs at low levels under normal conditions but increases transiently after stress, indicating a function of Alu RNAs in cellular stress response. Alu RNAs assemble with cellular proteins into ribonucleoprotein complexes and can be processed into the smaller scAlu RNAs. Alu and Alu-related RNAs play a role in regulating transcription and translation. They provide a source for the biogenesis of miRNAs and, embedded into mRNAs, can be targeted by miRNAs. When present as inverted repeats in mRNAs, they become substrates of the editing enzymes, and their modification causes the nuclear retention of these mRNAs. Certain Alu elements evolved into unique transcription units with specific expression profiles producing RNAs with highly specific cellular functions.

Age-dependent gain of alternative splice forms and biased duplication explain the relation between splicing and duplication

We analyze here the relation between alternative splicing and gene duplication in light of recent genomic data, with a focus on the human genome. We show that the previously reported negative correlation between level of alternative splicing and family size no longer holds true. We clarify this pattern and show that it is sufficiently explained by two factors. First, genes progressively gain new splice variants with time. The gain is consistent with a selectively relaxed regime, until purifying selection slows it down as aging genes accumulate a large number of variants. Second, we show that duplication does not lead to a loss of splice forms, but rather that genes with low levels of alternative splicing tend to duplicate more frequently. This leads us to reconsider the role of alternative splicing in duplicate retention.

Expression of ribosomal protein L22e family members in Drosophila melanogaster: rpL22-like is differentially expressed and alternatively spliced

Several ribosomal protein families contain paralogues whose roles may be equivalent or specialized to include extra-ribosomal functions. RpL22e family members rpL22 and rpL22-like are differentially expressed in Drosophila melanogaster: rpL22-like mRNA is gonad specific whereas rpL22 is expressed ubiquitously, suggesting distinctive paralogue functions. To determine if RpL22-like has a divergent role in gonads, rpL22-like expression was analysed by qRT-PCR and western blots, respectively, showing enrichment of rpL22-like mRNA and a 34 kDa (predicted) protein in testis, but not in ovary. Immunohistochemistry of the reproductive tract corroborated testis-specific expression. RpL22-like detection in 80S/polysome fractions from males establishes a role for this tissue-specific paralogue as a ribosomal component. Unpredictably, expression profiles revealed a low abundant, alternative mRNA variant (designated 'rpL22-like short') that would encode a novel protein lacking the C-terminal ribosomal protein signature but retaining part of the N-terminal domain. This variant results from splicing of a retained intron (defined by non-canonical splice sites) within rpL22-like mRNA. Polysome association and detection of a low abundant 13.5 kDa (predicted) protein in testis extracts suggests variant mRNA translation. Collectively, our data show that alternative splicing of rpL22-like generates structurally distinct protein products: ribosomal component RpL22-like and a novel protein with a role distinct from RpL22-like.

When one is better than two: RNA with dual functions

The central dogma of biology, until not long ago, held that genetic information stored on DNA molecules was translated into the final protein products through RNA as intermediate molecules. Then, an additional level of complexity in the regulation of genome expression was added, implicating new classes of RNA molecules called non-coding RNA (ncRNA). These ncRNA are also often referred to as functional RNA in that, although they do not contain the capacity to encode proteins, do have a function as RNA molecules. They have been thus far considered as truly non-coding RNA since no ORF long enough to be considered, nor protein, have been associated with them. However, the recent identification and characterization of bifunctional RNA, i.e. RNA for which both coding capacity and activity as functional RNA have been reported, suggests that a definite categorization of some RNA molecules is far from being straightforward. Indeed, several RNA primarily classified as non-protein-coding RNA has been showed to hold coding capacities and associated peptides. Conversely, mRNA, usually regarded as strictly protein-coding, may act as functional RNA molecules. Here, we describe several examples of these bifunctional RNA that have been already characterized from bacteria to mammals. We also extend this concept to fortuitous acquisition of dual function in pathological conditions and to the recently highlighted duality between information carried by a gene and its pseudogenes counterparts.

How do RNA sequence, DNA sequence, and chromatin properties regulate splicing?

Recent genome-wide studies have revealed a remarkable correspondence between nucleosome positions and exon-intron boundaries, and several studies have implicated specific histone modifications in regulating alternative splicing. In addition, recent progress in cracking the ‘splicing code’ shows that sequence motifs carried on the nascent RNA molecule itself are sufficient to accurately predict tissue-specific alternative splicing patterns. Together, these studies shed light on the complex interplay between RNA sequence, DNA sequence, and chromatin properties in regulating splicing.

Genetic causes and treatment of isolated growth hormone deficiency-an update

Isolated growth hormone deficiency is the most common pituitary hormone deficiency and can result from congenital or acquired causes, although the majority of cases are idiopathic with no identifiable etiology. Known genes involved in the genetic etiology of isolated growth hormone deficiency include those that encode growth hormone (GH1), growth-hormone-releasing hormone receptor (GHRHR) and transcription factor SOX3. However, mutations are identified in a relatively small percentage of patients, which suggests that other, yet unidentified, genetic factors are involved. Among the known factors, heterozygous mutations in GH1 remain the most frequent cause of isolated growth hormone deficiency. The identification of mutations has clinical implications for the management of patients with this condition, as individuals with heterozygous GH1 mutations vary in phenotype and can, in some cases, develop additional pituitary hormone deficiencies. Lifelong follow-up of these patients is, therefore, recommended. Further studies in the genetic etiology of isolated growth hormone deficiency will help to elucidate mechanisms implicated in the control of growth and may influence future treatment options. Advances in pharmacogenomics will also optimize the treatment of isolated growth hormone deficiency and other conditions associated with short stature, for which recombinant human growth hormone is a licensed therapy.

Assessment of orthologous splicing isoforms in human and mouse orthologous genes

Background

Recent discoveries have highlighted the fact that alternative splicing and alternative transcripts are the rule, rather than the exception, in metazoan genes. Since multiple transcript and protein variants expressed by the same gene are, by definition, structurally distinct and need not to be functionally equivalent, the concept of gene orthology should be extended to the transcript level in order to describe evolutionary relationships between structurally similar transcript variants. In other words, the identification of true orthology relationships between gene products now should progress beyond primary sequence and "splicing orthology", consisting in ancestrally shared exon-intron structures, is required to define orthologous isoforms at transcript level.

Results

As a starting step in this direction, in this work we performed a large scale human- mouse gene comparison with a twofold goal: first, to assess if and to which extent traditional gene annotations such as RefSeq capture genuine splicing orthology; second, to provide a more detailed annotation and quantification of true human-mouse orthologous transcripts defined as transcripts of orthologous genes exhibiting the same splicing patterns.

Conclusions

We observed an identical exon/intron structure for 32% of human and mouse orthologous genes. This figure increases to 87% using less stringent criteria for gene structure similarity, thus implying that for about 13% of the human RefSeq annotated genes (and about 25% of the corresponding transcripts) we could not identify any mouse transcript showing sufficient similarity to be confidently assigned as a splicing ortholog. Our data suggest that current gene and transcript data may still be rather incomplete - with several splicing variants still unknown. The observation that alternative splicing produces large numbers of alternative transcripts and proteins, some of them conserved across species and others truly species-specific, suggests that, still maintaining the conventional definition of gene orthology, a new concept of "splicing orthology" can be defined at transcript level.

Loss-of-function variants in the genomes of healthy humans

Genetic variants predicted to seriously disrupt the function of human protein-coding genes-so-called loss-of-function (LOF) variants-have traditionally been viewed in the context of severe Mendelian disease. However, recent large-scale sequencing and genotyping projects have revealed a surprisingly large number of these variants in the genomes of apparently healthy individuals--at least 100 per genome, including more than 30 in a homozygous state--suggesting a previously unappreciated level of variation in functional gene content between humans. These variants are mostly found at low frequency, suggesting that they are enriched for mildly deleterious polymorphisms suppressed by negative natural selection, and thus represent an attractive set of candidate variants for complex disease susceptibility. However, they are also enriched for sequencing and annotation artefacts, so overall present serious challenges for clinical sequencing projects seeking to identify severe disease genes amidst the 'noise' of technical error and benign genetic polymorphism. Systematic, high-quality catalogues of LOF variants present in the genomes of healthy individuals, built from the output of large-scale sequencing studies such as the 1000 Genomes Project, will help to distinguish between benign and disease-causing LOF variants, and will provide valuable resources for clinical genomics.