Evolution of tissue-specific expression of ancestral genes across vertebrates and insects.

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

Genome-wide transcriptomic changes reveal the genetic pathways involved in insect migration.

Insects are capable of extraordinary feats of long-distance movement that have profound impacts on the function of terrestrial ecosystems. The ability to undertake these movements arose multiple times through the evolution of a suite of traits that make up the migratory syndrome, however the underlying genetic pathways involved remain poorly understood. Migratory hoverflies (Diptera: Syrphidae) are an emerging model group for studies of migration. They undertake seasonal movements in huge numbers across large parts of the globe and are important pollinators, biological control agents and decomposers. Here, we assembled a high-quality draft genome of the marmalade hoverfly (Episyrphus balteatus). We leveraged this genomic resource to undertake a genome-wide transcriptomic comparison of actively migrating Episyrphus, captured from a high mountain pass as they flew south to overwinter, with the transcriptomes of summer forms which were non-migratory. We identified 1543 genes with very strong evidence for differential expression. Interrogation of this gene set reveals a remarkable range of roles in metabolism, muscle structure and function, hormonal regulation, immunity, stress resistance, flight and feeding behaviour, longevity, reproductive diapause and sensory perception. These features of the migrant phenotype have arisen by the integration and modification of pathways such as insulin signalling for diapause and longevity, JAK/SAT for immunity, and those leading to octopamine production and fuelling to boost flight capabilities. Our results provide a powerful genomic resource for future research, and paint a comprehensive picture of global expression changes in an actively migrating insect, identifying key genomic components involved in this important life-history strategy.

Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision.

Retinal photoreceptors have a distinct transcriptomic profile compared to other neuronal subtypes, likely reflecting their unique cellular morphology and function in the detection of light stimuli by way of the ciliary outer segment. We discovered a layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite the fact that both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors. Its deletion in zebrafish results in widespread down-regulation of microexon inclusion from early developmental stages, followed by other transcriptomic alterations, severe photoreceptor defects, and blindness. These results shed light on the transcriptomic specialization and functionality of photoreceptors, uncovering unique cell type-specific roles for Srrm3 and microexons with implications for retinal diseases.

The X-linked splicing regulator MBNL3 has been co-opted to restrict placental growth in eutherians.

Understanding the regulatory interactions that control gene expression during the development of novel tissues is a key goal of evolutionary developmental biology. Here, we show that Mbnl3 has undergone a striking process of evolutionary specialization in eutherian mammals resulting in the emergence of a novel placental function for the gene. Mbnl3 belongs to a family of RNA-binding proteins whose members regulate multiple aspects of RNA metabolism. We find that, in eutherians, while both Mbnl3 and its paralog Mbnl2 are strongly expressed in placenta, Mbnl3 expression has been lost from nonplacental tissues in association with the evolution of a novel promoter. Moreover, Mbnl3 has undergone accelerated protein sequence evolution leading to changes in its RNA-binding specificities and cellular localization. While Mbnl2 and Mbnl3 share partially redundant roles in regulating alternative splicing, polyadenylation site usage and, in turn, placenta maturation, Mbnl3 has also acquired novel biological functions. Specifically, Mbnl3 knockout (M3KO) alone results in increased placental growth associated with higher Myc expression. Furthermore, Mbnl3 loss increases fetal resource allocation during limiting conditions, suggesting that location of Mbnl3 on the X chromosome has led to its role in limiting placental growth, favoring the maternal side of the parental genetic conflict.

Environmental Enrichment Induces Epigenomic and Genome Organization Changes Relevant for Cognition.

In early development, the environment triggers mnemonic epigenomic programs resulting in memory and learning experiences to confer cognitive phenotypes into adulthood. To uncover how environmental stimulation impacts the epigenome and genome organization, we used the paradigm of environmental enrichment (EE) in young mice constantly receiving novel stimulation. We profiled epigenome and chromatin architecture in whole cortex and sorted neurons by deep-sequencing techniques. Specifically, we studied chromatin accessibility, gene and protein regulation, and 3D genome conformation, combined with predicted enhancer and chromatin interactions. We identified increased chromatin accessibility, transcription factor binding including CTCF-mediated insulation, differential occupancy of H3K36me3 and H3K79me2, and changes in transcriptional programs required for neuronal development. EE stimuli led to local genome re-organization by inducing increased contacts between chromosomes 7 and 17 (inter-chromosomal). Our findings support the notion that EE-induced learning and memory processes are directly associated with the epigenome and genome organization.

ACME dissociation: a versatile cell fixation-dissociation method for single-cell transcriptomics.

Single-cell sequencing technologies are revolutionizing biology, but they are limited by the need to dissociate live samples. Here, we present ACME (ACetic-MEthanol), a dissociation approach for single-cell transcriptomics that simultaneously fixes cells. ACME-dissociated cells have high RNA integrity, can be cryopreserved multiple times, and are sortable and permeable. As a proof of principle, we provide single-cell transcriptomic data of different species, using both droplet-based and combinatorial barcoding single-cell methods. ACME uses affordable reagents, can be done in most laboratories and even in the field, and thus will accelerate our knowledge of cell types across the tree of life.

Autism-Misregulated eIF4G Microexons Control Synaptic Translation and Higher Order Cognitive Functions.

Microexons represent the most highly conserved class of alternative splicing, yet their functions are poorly understood. Here, we focus on closely related neuronal microexons overlapping prion-like domains in the translation initiation factors, eIF4G1 and eIF4G3, the splicing of which is activity dependent and frequently disrupted in autism. CRISPR-Cas9 deletion of these microexons selectively upregulates synaptic proteins that control neuronal activity and plasticity and further triggers a gene expression program mirroring that of activated neurons. Mice lacking the Eif4g1 microexon display social behavior, learning, and memory deficits, accompanied by altered hippocampal synaptic plasticity. We provide evidence that the eIF4G microexons function as a translational brake by causing ribosome stalling, through their propensity to promote the coalescence of cytoplasmic granule components associated with translation repression, including the fragile X mental retardation protein FMRP. The results thus reveal an autism-disrupted mechanism by which alternative splicing specializes neuronal translation to control higher order cognitive functioning.

A novel protein domain in an ancestral splicing factor drove the evolution of neural microexons.

The mechanisms by which entire programmes of gene regulation emerged during evolution are poorly understood. Neuronal microexons represent the most conserved class of alternative splicing in vertebrates, and are critical for proper brain development and function. Here, we discover neural microexon programmes in non-vertebrate species and trace their origin to bilaterian ancestors through the emergence of a previously uncharacterized 'enhancer of microexons' (eMIC) protein domain. The eMIC domain originated as an alternative, neural-enriched splice isoform of the pan-eukaryotic Srrm2/SRm300 splicing factor gene, and subsequently became fixed in the vertebrate and neuronal-specific splicing regulator Srrm4/nSR100 and its paralogue Srrm3. Remarkably, the eMIC domain is necessary and sufficient for microexon splicing, and functions by interacting with the earliest components required for exon recognition. The emergence of a novel domain with restricted expression in the nervous system thus resulted in the evolution of splicing programmes that qualitatively expanded the neuronal molecular complexity in bilaterians.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

Evolutionary recruitment of flexible Esrp-dependent splicing programs into diverse embryonic morphogenetic processes.

Epithelial-mesenchymal interactions are crucial for the development of numerous animal structures. Thus, unraveling how molecular tools are recruited in different lineages to control interplays between these tissues is key to understanding morphogenetic evolution. Here, we study Esrp genes, which regulate extensive splicing programs and are essential for mammalian organogenesis. We find that Esrp homologs have been independently recruited for the development of multiple structures across deuterostomes. Although Esrp is involved in a wide variety of ontogenetic processes, our results suggest ancient roles in non-neural ectoderm and regulating specific mesenchymal-to-epithelial transitions in deuterostome ancestors. However, consistent with the extensive rewiring of Esrp-dependent splicing programs between phyla, most developmental defects observed in vertebrate mutants are related to other types of morphogenetic processes. This is likely connected to the origin of an event in Fgfr, which was recruited as an Esrp target in stem chordates and subsequently co-opted into the development of many novel traits in vertebrates.

An atlas of alternative splicing profiles and functional associations reveals new regulatory programs and genes that simultaneously express multiple major isoforms.

Alternative splicing (AS) generates remarkable regulatory and proteomic complexity in metazoans. However, the functions of most AS events are not known, and programs of regulated splicing remain to be identified. To address these challenges, we describe the Vertebrate Alternative Splicing and Transcription Database (VastDB), the largest resource of genome-wide, quantitative profiles of AS events assembled to date. VastDB provides readily accessible quantitative information on the inclusion levels and functional associations of AS events detected in RNA-seq data from diverse vertebrate cell and tissue types, as well as developmental stages. The VastDB profiles reveal extensive new intergenic and intragenic regulatory relationships among different classes of AS and previously unknown and conserved landscapes of tissue-regulated exons. Contrary to recent reports concluding that nearly all human genes express a single major isoform, VastDB provides evidence that at least 48% of multiexonic protein-coding genes express multiple splice variants that are highly regulated in a cell/tissue-specific manner, and that >18% of genes simultaneously express multiple major isoforms across diverse cell and tissue types. Isoforms encoded by the latter set of genes are generally coexpressed in the same cells and are often engaged by translating ribosomes. Moreover, they are encoded by genes that are significantly enriched in functions associated with transcriptional control, implying they may have an important and wide-ranging role in controlling cellular activities. VastDB thus provides an unprecedented resource for investigations of AS function and regulation.

Complex transcriptional regulation and independent evolution of fungal-like traits in a relative of animals.

Cell-type specification through differential genome regulation is a hallmark of complex multicellularity. However, it remains unclear how this process evolved during the transition from unicellular to multicellular organisms. To address this question, we investigated transcriptional dynamics in the ichthyosporean Creolimax fragrantissima, a relative of animals that undergoes coenocytic development. We find that Creolimax utilizes dynamic regulation of alternative splicing, long inter-genic non-coding RNAs and co-regulated gene modules associated with animal multicellularity in a cell-type specific manner. Moreover, our study suggests that the different cell types of the three closest animal relatives (ichthyosporeans, filastereans and choanoflagellates) are the product of lineage-specific innovations. Additionally, a proteomic survey of the secretome reveals adaptations to a fungal-like lifestyle. In summary, the diversity of cell types among protistan relatives of animals and their complex genome regulation demonstrates that the last unicellular ancestor of animals was already capable of elaborate specification of cell types.

The origin of patterning systems in bilateria-insights from the Hox and ParaHox genes in Acoelomorpha.

Hox and ParaHox genes constitute two families of developmental regulators that pattern the Anterior-Posterior body axis in all bilaterians. The members of these two groups of genes are usually arranged in genomic clusters and work in a coordinated fashion, both in space and in time. While the mechanistic aspects of their action are relatively well known, it is still unclear how these systems evolved. For instance, we still need a proper model of how the Hox and ParaHox clusters were assembled over time. This problem is due to the shortage of information on gene complements for many taxa (mainly basal metazoans) and the lack of a consensus phylogenetic model of animal relationships to which we can relate our new findings. Recently, several studies have shown that the Acoelomorpha most probably represent the first offshoot of the Bilateria. This finding has prompted us, and others, to study the Hox and ParaHox complements in these animals, as well as their activity during development. In this review, we analyze how the current knowledge of Hox and ParaHox genes in the Acoelomorpha is shaping our view of bilaterian evolution.

High transcriptional complexity of the retinitis pigmentosa CERKL gene in human and mouse.

PURPOSE: To shed light on the pathogenicity of the mutations in the retinitis pigmentosa gene CERKL, the authors aimed to characterize its transcriptional repertoire and focused on the use of distinct promoters and alternative splicing in human and mouse tissues. METHODS: In silico genomic and transcriptomic computational customized analysis, combined with experimental RT-PCRs on different human and murine tissues and cell lines and immunohistochemistry, have been used to characterize the transcriptional spectrum of CERKL. In the mouse retina, Cerkl is detected primarily in ganglion cells and cones but can also be observed in rods. Cerkl is mainly cytosolic. It localizes in the outer segments of photoreceptors and in the perinuclear regions of some cells. RESULTS: An unexpected multiplicity of CERKL transcriptional start sites (four in each species) plus a high variety of alternative splicing events primarily affecting the 5' half of the gene generate >20 fully validated mRNA isoforms in human and 23 in mouse. Moreover, several translational start sites, compatible with a wide display of functional domains, contribute to the final protein complexity. CONCLUSIONS: This combined approach of in silico and experimental characterization of the CERKL gene provides a comprehensive picture of the species-specific transcriptional products in the retina, underscores highly tuned gene regulation in different tissues, and establishes a framework for the study of CERKL genotype-phenotype correlations.

Autosomal recessive retinitis pigmentosa with early macular affectation caused by premature truncation in PROM1.

PURPOSE: To identify the genetic basis of a large consanguineous Spanish pedigree affected with autosomal recessive retinitis pigmentosa (arRP) with premature macular atrophy and myopia. METHODS: After a high-throughput cosegregation gene chip was used to exclude all known RP and Leber congenital amaurosis (LCA) candidates, genome-wide screening and linkage analysis were performed. Direct mutational screening identified the pathogenic mutation, and primers were designed to obtain the RT-PCR products for isoform characterization. RESULTS: Mutational analysis detected a novel homozygous PROM1 mutation, c.869delG in exon 8 cosegregating with the disease. This variant causes a frameshift that introduces a premature stop codon, producing truncation of approximately two-thirds of the protein. Analysis of PROM1 expression in the lymphocytes of patients, carriers, and control subjects revealed an aberrant transcript that is degraded by the nonsense-mediated decay pathway, suggesting that the disease is caused by the absence of the PROM1 protein. Three (s2, s11 and s12) of the seven alternatively spliced isoforms reported in humans, accounted for 98% of the transcripts in the retina. Given that these three contained exon 8, no PROM1 isoform is expected in the affected retinas. CONCLUSIONS: A remarkable clinical finding in the affected family is early macular atrophy with concentric spared areas. The authors propose that the hallmark of PROM1 truncating mutations is early and severe progressive degeneration of both rods and cones and highlight this gene as a candidate of choice to prioritize in the molecular genetic study of patients with noncanonical clinical peripheral and macular affectation.

Comprehensive SNP-chip for retinitis pigmentosa-Leber congenital amaurosis diagnosis: new mutations and detection of mutational founder effects.

Fast and efficient high-throughput techniques are essential for the molecular diagnosis of highly heterogeneous hereditary diseases, such as retinitis pigmentosa (RP). We had previously approached RP genetic testing by devising a chip based on co-segregation analysis for the autosomal recessive forms. In this study, we aimed to design a diagnostic tool for all the known genes (40 up to now) responsible for the autosomal dominant and recessive RP and Leber congenital amaurosis (LCA). This new chip analyzes 240 single nucleotide polymorphisms (SNPs) (6 per gene) on a high-throughput genotyping platform (SNPlex, Applied Biosystems), and genetic diagnosis is based on the co-segregation analysis of SNP haplotypes in independent families. In a single genotyping step, the number of RP candidates to be screened for mutations is considerably reduced, and in the most informative families, all the candidates are ruled out at once. In a panel of RP Spanish pedigrees, the disease chip became a crucial tool for selecting those suitable for genome-wide RP gene search, and saved the burdensome direct mutational screening of every known RP gene. In a large adRP family, the chip allowed ruling out of all but the causative gene, and identification of an unreported null mutation (E181X) in PRPF31. Finally, on the basis of the conservation of the SNP haplotype linked to this pathogenic variant, we propose that the E181X mutation spread through a cohort of geographically isolated families by a founder effect.

Getting closer to a pre-vertebrate genome: the non-LTR retrotransposons of Branchiostoma floridae.

Non-LTR retrotransposons are common in vertebrate genomes and although present in invertebrates they appear at a much lower frequency. The cephalochordate amphioxus is the closest living relative to vertebrates and has been considered a good model for comparative analyses of genome expansions during vertebrate evolution. With the aim to assess the involvement of transposable elements in these events, we have analysed the non-LTR retrotransposons of Branchiostoma floridae. In silico searches have allowed to reconstruct non-LTR elements of six different clades (CR1, I, L1, L2, NeSL and RTE) and assess their structural features. According to the estimated copy number of these elements they account for less than 1% of the haploid genome, which reminds of the low abundance also encountered in the urochordate Ciona intestinalis. Amphioxus (B. floridae) and Ciona share a pre-vertebrate-like organization for the non-LTR retrotransposons (<150 copies, < 1% of the genome) versus the complexity associated to higher vertebrates (Homo sapiens >1.3.10(6) copies, > 20% of the genome).

Evolutionary genomics of the recently duplicated amphioxus Hairy genes.

Amphioxus Hairy genes have gone through a number of lineage-specific duplications, resulting in eight members, some of which are differentially expressed in the embryo. In order to gain insights into the evolution and function of this gene family we have compared their genomic structure and searched for conserved non-coding sequence elements. We have found that introns have been lost independently from these genes at least twice and after the duplication events. By carrying out phylogenetic footprinting between paralogues expressed in the embryo, we have found a differential distribution of conserved elements that could explain the limited overlap in expression patterns of Hairy genes in the amphioxus embryo. Furthermore, clustering of RBP-Jk binding sites in these conserved elements suggests that amphioxus Hairy genes are downstream targets of the Notch signaling pathway, as occurs in vertebrates. All of this evidence suggests that amphioxus Hairy genes have gone through a process of subfunctionalization shortly after their duplication, representing an extreme and rapid case of the duplication-degeneration-complementation model.

Preliminary observations on the spawning conditions of the European amphioxus (Branchiostoma lanceolatum) in captivity.

Members of the subphylum Cephalochordata, which include the genus Branchiostoma (i.e. amphioxus), represent the closest living invertebrate relatives of the vertebrates. To date, developmental studies have been carried out on three amphioxus species (the European Branchiostoma lanceolatum, the East Asian B. belcheri, and Floridian-Caribbean B. floridae). In most instances, adult animals have been collected from the field during their ripe season and allowed (or stimulated) to spawn in the laboratory. In any given year, dates of laboratory pawning have been limited by two factors. First, natural populations of these three most studied species of amphioxus are ripe, at most, for only a couple of months each year and, second, even when apparently ripe, animals spawn only at unpredictable intervals of every several days. This limited supply of living material hinders the development of amphioxus as a model system because this limitation makes it more difficult to work out protocols for new laboratory techniques. Therefore we are developing laboratory methods for increasing the number of amphioxus spawning dates per year. The present study found that a Mediterranean population of B. lanceolatum living near the Franco-Spanish border spawned naturally at the end of May and again at the end of June in 2003. Re-feeding experiments in the laboratory demonstrated that the gonads emptied at the end of May refilled with gametes by the end of June. We also found that animals with large gonads (both, obtained from the field and kept and fed at the laboratory during several weeks) could be induced to spawn in the laboratory out of phase with the field population if they were temperature shocked (spawning occurred 36 hours after a sustained increase in water temperature from 19 degrees C to 25 degrees C).

The non-LTR retrotransposons in Ciona intestinalis: new insights into the evolution of chordate genomes.

BACKGROUND: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. RESULTS: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. CONCLUSIONS: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.