Evidence for widespread subfunctionalization of splice forms in vertebrate genomes

Translational control of gene function through optically regulated nucleic acids

Translation of mRNA into protein is one of the most fundamental processes within biological systems. Gene expression is tightly regulated both in space and time, often involving complex signaling or gene regulatory networks, as most prominently observed in embryo development. Thus, studies of gene function require tools with a matching level of external control. Light is an excellent conditional trigger as it is minimally invasive, can be easily tuned in wavelength and amplitude, and can be applied with excellent spatial and temporal resolution. To this end, modification of established oligonucleotide-based technologies with optical control elements, in the form of photocaging groups and photoswitches, has rendered these tools capable of navigating the dynamic regulatory pathways of mRNA translation in cellular and in vivo models. In this review, we discuss the different optochemical approaches used to generate photoresponsive nucleic acids that activate and deactivate gene expression and function at the translational level.

Interactome Analysis of KIN (Kin17) Shows New Functions of This Protein

KIN (Kin17) protein is overexpressed in a number of cancerous cell lines, and is therefore considered a possible cancer biomarker. It is a well-conserved protein across eukaryotes and is ubiquitously expressed in all cell types studied, suggesting an important role in the maintenance of basic cellular function which is yet to be well determined. Early studies on KIN suggested that this nuclear protein plays a role in cellular mechanisms such as DNA replication and/or repair; however, its association with chromatin depends on its methylation state. In order to provide a better understanding of the cellular role of this protein, we investigated its interactome by proximity-dependent biotin identification coupled to mass spectrometry (BioID-MS), used for identification of protein-protein interactions. Our analyses detected interaction with a novel set of proteins and reinforced previous observations linking KIN to factors involved in RNA processing, notably pre-mRNA splicing and ribosome biogenesis. However, little evidence supports that this protein is directly coupled to DNA replication and/or repair processes, as previously suggested. Furthermore, a novel interaction was observed with PRMT7 (protein arginine methyltransferase 7) and we demonstrated that KIN is modified by this enzyme. This interactome analysis indicates that KIN is associated with several cell metabolism functions, and shows for the first time an association with ribosome biogenesis, suggesting that KIN is likely a moonlight protein.

AMPA receptor auxiliary subunits emerged during early vertebrate evolution by neo/subfunctionalization of unrelated proteins

In mammalian synapses, the function of ionotropic glutamate receptors is critically modulated by auxiliary subunits. Most of these specifically regulate the synaptic localization and electrophysiological properties of AMPA-type glutamate receptors (AMPARs). Here, we comprehensively investigated the animal evolution of the protein families that contain AMPAR auxiliary subunits (ARASs). We observed that, on average, vertebrates have four times more ARASs than other animal species. We also demonstrated that ARASs belong to four unrelated protein families: CACNG-GSG1, cornichon, shisa and Dispanin C. Our study demonstrates that, despite the ancient origin of these four protein families, the majority of ARASs emerged during vertebrate evolution by independent but convergent processes of neo/subfunctionalization that resulted in the multiple ARASs found in present vertebrate genomes. Importantly, although AMPARs appeared and diversified in the ancestor of bilateral animals, the ARAS expansion did not occur until much later, in early vertebrate evolution. We propose that the surge in ARASs and consequent increase in AMPAR functionalities, contributed to the increased complexity of vertebrate brains and cognitive functions.

Similarities in biological processes can be used to bridge ecology and molecular biology

Much of the research in biology aims to understand the origin of diversity. Naturally, ecological diversity was the first object of study, but we now have the necessary tools to probe diversity at molecular scales. The inherent differences in how we study diversity at different scales caused the disciplines of biology to be organized around these levels, from molecular biology to ecology. Here, we illustrate that there are key properties of each scale that emerge from the interactions of simpler components and that these properties are often shared across different levels of organization. This means that ideas from one level of organization can be an inspiration for novel hypotheses to study phenomena at another level. We illustrate this concept with examples of events at the molecular level that have analogs at the organismal or ecological level and vice versa. Through these examples, we illustrate that biological processes at different organization levels are governed by general rules. The study of the same phenomena at different scales could enrich our work through a multidisciplinary approach, which should be a staple in the training of future scientists.

Modeling Neuronal Diseases in Zebrafish in the Era of CRISPR

BACKGROUND

Danio rerio is a powerful experimental model for studies in genetics and development. Recently, CRISPR technology has been applied in this species to mimic various human diseases, including those affecting the nervous system. Zebrafish offer multiple experimental advantages: external embryogenesis, rapid development, transparent embryos, short life cycle, and basic neurobiological processes shared with humans. This animal model, together with the CRISPR system, emerging imaging technologies, and novel behavioral approaches, lay the basis for a prominent future in neuropathology and will undoubtedly accelerate our understanding of brain function and its disorders.

OBJECTIVE

Gather relevant findings from studies that have used CRISPR technologies in zebrafish to explore basic neuronal function and model human diseases.

METHODS

We systematically reviewed the most recent literature about CRISPR technology applications for understanding brain function and neurological disorders in D. rerio. We highlighted the key role of CRISPR in driving forward our understanding of particular topics in neuroscience.

RESULTS

We show specific advances in neurobiology when the CRISPR system has been applied in zebrafish and describe how CRISPR is accelerating our understanding of brain organization.

CONCLUSION

Today, CRISPR is the preferred method to modify genomes of practically any living organism. Despite the rapid development of CRISPR technologies to generate disease models in zebrafish, more efforts are needed to efficiently combine different disciplines to find the etiology and treatments for many brain diseases.

Evolutionary Dynamics of Structural Variation at a Key Locus for Color Pattern Diversification in Cichlid Fishes

Color patterns in African cichlid fishes vary spectacularly. Although phylogenetic analysis showed already 30 years ago that many color patterns evolved repeatedly in these adaptive radiations, only recently have we begun to understand the genomic basis of color variation. Horizontal stripe patterns evolved and were lost several times independently across the adaptive radiations of Lake Victoria, Malawi, and Tanganyika and regulatory evolution of agouti-related peptide 2 (agrp2/asip2b) has been linked to this phenotypically labile trait. Here, we asked whether the agrp2 locus exhibits particular characteristics that facilitate divergence in color patterns. Based on comparative genomic analyses, we discovered several recent duplications, insertions, and deletions. Interestingly, one of these events resulted in a tandem duplication of the last exon of agrp2. The duplication likely precedes the East African radiations that started 8-12 Ma, is not fixed within any of the radiations, and is found to vary even within some species. Moreover, we also observed variation in copy number (two to five copies) and secondary loss of the duplication, illustrating a surprising dynamic at this locus that possibly promoted functional divergence of agrp2. Our work suggests that such instances of exon duplications are a neglected mechanism potentially involved in the repeated evolution and diversification that deserves more attention.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Conserved functions of RNA-binding proteins in muscle

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

A Two-Headed Monster to Avert Disaster: HBS1/SKI7 Is Alternatively Spliced to Build Eukaryotic RNA Surveillance Complexes

The cytosolic RNA exosome, a 3'→5' exoribonuclease complex, contributes to mRNA degradation in eukaryotes, limiting the accumulation of poorly-translated, improperly translated, or aberrant mRNA species. Disruption of cytosolic RNA exosome activity allows aberrant RNA species to accumulate, which can then be detected by host antiviral immune systems as a signature of pathogen infection, activating antiviral defenses. SKI7 is a critical component of the cytosolic RNA exosome in yeast, bridging the catalytic exoribonuclease core with the SKI2/SKI3/SKI8 adaptor complex that guides aberrant RNA substrates into the exosome. The ortholog of SKI7 was only recently identified in humans as an alternative splice form of the HBS1 gene, which encodes a decoding factor translational GTPase that rescues stalled ribosomes. Here, we identify the plant orthologs of HBS1/SKI7. We found that HBS1 and SKI7 are typically encoded by alternative splice forms of a single locus, although some plant lineages have evolved subfunctionalized genes that apparently encode only HBS1 or only SKI7. In all plant lineages examined, the SKI7 gene is subject to regulation by alternative splicing that can yield unproductive transcripts, either by removing deeply conserved SKI7 coding sequences, or by introducing premature stop codons that render SKI7 susceptible to nonsense-mediated decay. Taking a comparative, evolutionary approach, we define crucial features of the SKI7 protein shared by all eukaryotes, and use these deeply conserved features to identify SKI7 proteins in invertebrate lineages. We conclude that SKI7 is a conserved cytosolic RNA exosome subunit across eukaryotic lineages, and that SKI7 is consistently regulated by alternative splicing, suggesting broad coordination of nuclear and cytosolic RNA metabolism.

Zebrafish Models of Rare Hereditary Pediatric Diseases

Recent advances in sequencing technologies have made it significantly easier to find the genetic roots of rare hereditary pediatric diseases. These novel methods are not panaceas, however, and they often give ambiguous results, highlighting multiple possible causative mutations in affected patients. Furthermore, even when the mapping results are unambiguous, the affected gene might be of unknown function. In these cases, understanding how a particular genotype can result in a phenotype also needs carefully designed experimental work. Model organism genetics can offer a straightforward experimental setup for hypothesis testing. Containing orthologs for over 80% of the genes involved in human diseases, zebrafish (Danio rerio) has emerged as one of the top disease models over the past decade. A plethora of genetic tools makes it easy to create mutations in almost any gene of the zebrafish genome and these mutant strains can be used in high-throughput preclinical screens for active molecules. As this small vertebrate species offers several other advantages as well, its popularity in biomedical research is bound to increase, with “aquarium to bedside” drug development pipelines taking a more prevalent role in the near future.

Alternative splicing and the evolution of phenotypic novelty

Alternative splicing, a mechanism of post-transcriptional RNA processing whereby a single gene can encode multiple distinct transcripts, has been proposed to underlie morphological innovations in multicellular organisms. Genes with developmental functions are enriched for alternative splicing events, suggestive of a contribution of alternative splicing to developmental programmes. The role of alternative splicing as a source of transcript diversification has previously been compared to that of gene duplication, with the relationship between the two extensively explored. Alternative splicing is reduced following gene duplication with the retention of duplicate copies higher for genes which were alternatively spliced prior to duplication. Furthermore, and unlike the case for overall gene number, the proportion of alternatively spliced genes has also increased in line with the evolutionary diversification of cell types, suggesting alternative splicing may contribute to the complexity of developmental programmes. Together these observations suggest a prominent role for alternative splicing as a source of functional innovation. However, it is unknown whether the proliferation of alternative splicing events indeed reflects a functional expansion of the transcriptome or instead results from weaker selection acting on larger species, which tend to have a higher number of cell types and lower population sizes.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

A targeted gene expression platform allows for rapid analysis of chemical-induced antioxidant mRNA expression in zebrafish larvae

Chemical-induced oxidative stress and the biochemical pathways that protect against oxidative damage are of particular interest in the field of toxicology. To rapidly identify oxidative stress-responsive gene expression changes in zebrafish, we developed a targeted panel of antioxidant genes using the Affymetrix QuantiGene Plex (QGP) platform. The genes contained in our panel include eight putative Nrf2 (Nfe2l2a)-dependent antioxidant genes (hmox1a, gstp1, gclc, nqo1, prdx1, gpx1a, sod1, sod2), a stress response gene (hsp70), an inducible DNA damage repair gene (gadd45bb), and three reference genes (actb1, gapdh, hprt1). We tested this platform on larval zebrafish exposed to tert-butyl hydroperoxide (tBHP) and cadmium (Cd), two model oxidative stressors with different modes of action, and compared our results with those obtained using the more common quantitative PCR (qPCR) method. Both methods showed that exposure to tBHP and Cd induced expression of prdx1, gstp1, and hmox1a (2- to 12-fold increase via QGP), indicative of an activated Nrf2 response in larval zebrafish. Both compounds also elicited a general stress response as reflected by elevation of hsp70 and gadd45bb, with Cd being the more potent inducer. Transient changes were observed in sod2 and gpx1a expression, whereas nqo1, an Nrf2-responsive gene in mammalian cells, was minimally affected by either tBHP or Cd chemical exposures. Developmental expression analysis of the target genes by QGP revealed marked upregulation of sod2 between 0-96hpf, and to a lesser extent, of sod1 and gstp1. Once optimized, QGP analysis of these experiments was accomplished more rapidly, using far less tissue, and at lower total costs than qPCR analysis. In summary, the QGP platform as applied to higher-throughput zebrafish studies provides a reasonable cost-effective alternative to qPCR or more comprehensive transcriptomics approaches to rapidly assess the potential for chemicals to elicit oxidative stress as a mechanism of chemical toxicity.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

The protein diversity that exists today has resulted from various evolutionary processes. It is well known that gene duplication (GD) along with the accumulation of mutations are responsible, among other factors, for an increase in the number of different proteins. The gene structure in eukaryotes requires the removal of non-coding sequences, introns, to produce mature mRNAs. This process, known as cis-splicing, referred to here as splicing, is regulated by several factors which can lead to numerous splicing arrangements, commonly designated as alternative splicing (AS). AS, producing several transcripts isoforms form a single gene, also increases the protein diversity. However, the evolution and manner for increasing protein variation differs between AS and GD. An important question is how are patterns of AS affected after a GD event. Here, we review the current knowledge of AS and GD, focusing on their evolutionary relationship. These two processes are now considered the main contributors to the increasing protein diversity and therefore their relationship is a relevant, yet understudied, area of evolutionary study.

Evolutionary and Functional Analysis of Membrane-Bound NAC Transcription Factor Genes in Soybean

Functional divergence is thought to be an important evolutionary driving force for the retention of duplicate genes. We reconstructed the evolutionary history of soybean (Glycine max) membrane-bound NAC transcription factor (NTL) genes. NTLs are thought to be components of stress signaling and unique in their requirement for proteolytic cleavage to free them from the membrane. Most of the 15 GmNTL genes appear to have evolved under strong purifying selection. By analyzing the phylogenetic tree and gene synteny, we identified seven duplicate gene pairs generated by the latest whole-genome duplication. The members of each pair were shown to have variously diverged at the transcriptional (organ specificity and responsiveness to stress), posttranscriptional (alternative splicing), and protein (proteolysis-mediated membrane release and transactivation activity) levels. The dormant (full-length protein) and active (protein without a transmembrane motif) forms of one pair of duplicated gene products (GmNTL1/GmNLT11) were each separately constitutively expressed in Arabidopsis (Arabidopsis thaliana). The heteroexpression of active but not dormant forms of these proteins caused improved tolerance to abiotic stresses, suggesting that membrane release was required for their functionality. Arabidopsis carrying the dormant form of GmNTL1 was more tolerant to hydrogen peroxide, which induces its membrane release. Tolerance was not increased in the line carrying dormant GmNTL11, which was not released by hydrogen peroxide treatment. Thus, NTL-release pattern changes may cause phenotypic divergence. It was concluded that a variety of functional divergences contributed to the retention of these GmNTL duplicates.

Coregulation of tandem duplicate genes slows evolution of subfunctionalization in mammals

Gene duplication is a fundamental process in genome evolution. However, most young duplicates are degraded by loss-of-function mutations, and the factors that allow some duplicate pairs to survive long-term remain controversial. One class of models to explain duplicate retention invokes sub- or neofunctionalization, whereas others focus on sharing of gene dosage. RNA-sequencing data from 46 human and 26 mouse tissues indicate that subfunctionalization of expression evolves slowly and is rare among duplicates that arose within the placental mammals, possibly because tandem duplicates are coregulated by shared genomic elements. Instead, consistent with the dosage-sharing hypothesis, most young duplicates are down-regulated to match expression levels of single-copy genes. Thus, dosage sharing of expression allows for the initial survival of mammalian duplicates, followed by slower functional adaptation enabling long-term preservation.

Functional divergence and convergence between the transcript network and gene network in lung adenocarcinoma

INTRODUCTION

Alternative RNA splicing is a critical regulatory mechanism during tumorigenesis. However, previous oncological studies mainly focused on the splicing of individual genes. Whether and how transcript isoforms are coordinated to affect cellular functions remain underexplored. Also of great interest is how the splicing regulome cooperates with the transcription regulome to facilitate tumorigenesis. The answers to these questions are of fundamental importance to cancer biology.

RESULTS

Here, we report a comparative study between the transcript-based network (TN) and the gene-based network (GN) derived from the transcriptomes of paired tumor-normal tissues from 77 lung adenocarcinoma patients. We demonstrate that the two networks differ significantly from each other in terms of patient clustering and the number and functions of network modules. Interestingly, the majority (89.5%) of multi-transcript genes have their transcript isoforms distributed in at least two TN modules, suggesting regulatory and functional divergences between transcript isoforms. Furthermore, TN and GN modules share onlŷ50%-60% of their biological functions. TN thus appears to constitute a regulatory layer separate from GN. Nevertheless, our results indicate that functional convergence and divergence both occur between TN and GN, implying complex interactions between the two regulatory layers. Finally, we report that the expression profiles of module members in both TN and GN shift dramatically yet concordantly during tumorigenesis. The mechanisms underlying this coordinated shifting remain unclear yet are worth further explorations.

CONCLUSION

We show that in lung adenocarcinoma, transcript isoforms per se are coordinately regulated to conduct biological functions not conveyed by the network of genes. However, the two networks may interact closely with each other by sharing the same or related biological functions. Unraveling the effects and mechanisms of such interactions will significantly advance our understanding of this deadly disease.

Evolution of selenophosphate synthetases: emergence and relocation of function through independent duplications and recurrent subfunctionalization

Selenoproteins are proteins that incorporate selenocysteine (Sec), a nonstandard amino acid encoded by UGA, normally a stop codon. Sec synthesis requires the enzyme Selenophosphate synthetase (SPS or SelD), conserved in all prokaryotic and eukaryotic genomes encoding selenoproteins. Here, we study the evolutionary history of SPS genes, providing a map of selenoprotein function spanning the whole tree of life. SPS is itself a selenoprotein in many species, although functionally equivalent homologs that replace the Sec site with cysteine (Cys) are common. Many metazoans, however, possess SPS genes with substitutions other than Sec or Cys (collectively referred to as SPS1). Using complementation assays in fly mutants, we show that these genes share a common function, which appears to be distinct from the synthesis of selenophosphate carried out by the Sec- and Cys- SPS genes (termed SPS2), and unrelated to Sec synthesis. We show here that SPS1 genes originated through a number of independent gene duplications from an ancestral metazoan selenoprotein SPS2 gene that most likely already carried the SPS1 function. Thus, in SPS genes, parallel duplications and subsequent convergent subfunctionalization have resulted in the segregation to different loci of functions initially carried by a single gene. This evolutionary history constitutes a remarkable example of emergence and evolution of gene function, which we have been able to trace thanks to the singular features of SPS genes, wherein the amino acid at a single site determines unequivocally protein function and is intertwined to the evolutionary fate of the entire selenoproteome.

The evolutionary fate of alternatively spliced homologous exons after gene duplication

Alternative splicing and gene duplication are the two main processes responsible for expanding protein functional diversity. Although gene duplication can generate new genes and alternative splicing can introduce variation through alternative gene products, the interplay between the two processes is complex and poorly understood. Here, we have carried out a study of the evolution of alternatively spliced exons after gene duplication to better understand the interaction between the two processes. We created a manually curated set of 97 human genes with mutually exclusively spliced homologous exons and analyzed the evolution of these exons across five distantly related vertebrates (lamprey, spotted gar, zebrafish, fugu, and coelacanth). Most of these exons had an ancient origin (more than 400 Ma). We found examples supporting two extreme evolutionary models for the behaviour of homologous axons after gene duplication. We observed 11 events in which gene duplication was accompanied by splice isoform separation, that is, each paralog specifically conserved just one distinct ancestral homologous exon. At other extreme, we identified genes in which the homologous exons were always conserved within paralogs, suggesting that the alternative splicing event cannot easily be separated from the function in these genes. That many homologous exons fall in between these two extremes highlights the diversity of biological systems and suggests that the subtle balance between alternative splicing and gene duplication is adjusted to the specific cellular context of each gene.