Sex-biased gene expression in mammals.

Sex differences in gene expression start at puberty and vary across species and organs.

Distinct gene expression dynamics in developing and regenerating crustacean limbs.

Regenerating animals have the ability to reproduce body parts that were originally made in the embryo and subsequently lost due to injury. Understanding whether regeneration mirrors development is an open question in most regenerative species. Here, we take a transcriptomics approach to examine whether leg regeneration shows similar temporal patterns of gene expression as leg development in the embryo, in the crustacean Parhyale hawaiensis. We find that leg development in the embryo shows stereotypic temporal patterns of gene expression. In contrast, the dynamics of gene expression during leg regeneration show a higher degree of variation related to the physiology of individual animals. A major driver of this variation is the molting cycle. We dissect the transcriptional signals of individual physiology and regeneration to obtain clearer temporal signals marking distinct phases of leg regeneration. Comparing the transcriptional dynamics of development and regeneration we find that, although the two processes use similar sets of genes, the temporal patterns in which these genes are deployed are different and cannot be systematically aligned.

Tissue-specific patterns of regulatory changes underlying gene expression differences among Ficedula flycatchers and their naturally occurring F1 hybrids.

Changes in interacting cis- and trans-regulatory elements are important candidates for Dobzhansky-Muller hybrid incompatibilities and may contribute to hybrid dysfunction by giving rise to misexpression in hybrids. To gain insight into the molecular mechanisms and determinants of gene expression evolution in natural populations, we analyzed the transcriptome from multiple tissues of two recently diverged Ficedula flycatcher species and their naturally occurring F1 hybrids. Differential gene expression analysis revealed that the extent of differentiation between species and the set of differentially expressed genes varied across tissues. Common to all tissues, a higher proportion of Z-linked genes than autosomal genes showed differential expression, providing evidence for a fast-Z effect. We further found clear signatures of hybrid misexpression in brain, heart, kidney, and liver. However, while testis showed the highest divergence of gene expression among tissues, it showed no clear signature of misexpression in F1 hybrids, even though these hybrids were found to be sterile. It is therefore unlikely that incompatibilities between cis-trans regulatory changes explain the observed sterility. Instead, we found evidence that cis-regulatory changes play a significant role in the evolution of gene expression in testis, which illustrates the tissue-specific nature of cis-regulatory evolution bypassing constraints associated with pleiotropic effects of genes.

Developmental variability channels mouse molar evolution.

Do developmental systems preferentially produce certain types of variation that orient phenotypic evolution along preferred directions? At different scales, from the intra-population to the interspecific, the murine first upper molar shows repeated anterior elongation. Using a novel quantitative approach to compare the development of two mouse strains with short or long molars, we identified temporal, spatial and functional differences in tooth signaling center activity, that arise from differential tuning of the activation-inhibition mechanisms underlying tooth patterning. By tracing their fate, we could explain why only the upper first molar reacts via elongation of its anterior part. Despite a lack of genetic variation, individuals of the elongated strain varied in tooth length and the temporal dynamics of their signaling centers, highlighting the intrinsic instability of the upper molar developmental system. Collectively, these results reveal the variational properties of murine molar development that drive morphological evolution along a line of least resistance.

Over time species develop random mutations in their genetic sequence that causes their form to change. If this new form increases the survival of a species it will become favored through natural selection and is more likely to get passed on to future generations. But, the evolution of these new traits also depends on what happens during development. Developmental mechanisms control how an embryo progresses from a single cell to an adult organism made of many cells. Mutations that alter these processes can influence the physical outcome of development, and cause a new trait to form. This means that if many different mutations alter development in a similar way, this can lead to the same physical change, making it 'easy' for a new trait to repeatedly occur. Most of the research has focused on finding the mutations that underlie repeated evolution, but rarely on identifying the role of the underlying developmental mechanisms. To bridge this gap, Hayden et al. investigated how changes during development influence the shape and size of molar teeth in mice. In some wild species of mice, the front part of the first upper molar is longer than in other species. This elongation, which is repeatedly found in mice from different islands, likely came from developmental mechanisms. Tooth development in mice has been well-studied in the laboratory, and Hayden et al. started by identifying two strains of laboratory mice that mimic the teeth seen in their wild cousins, one with elongated upper first molars and another with short ones. Comparing how these two strains of mice developed their elongated or short teeth revealed key differences in the embryonic structures that form the upper molar and cause it to elongate. Further work showed that variations in these embryonic structures can even cause mice that are genetically identical to have longer or shorter upper first molars. These findings show how early differences during development can lead to small variations in form between adult species of mice. This study highlights how studying developmental differences as well as genetic sequences can further our understanding of how different species evolved.

Detecting adaptive convergent amino acid evolution.

In evolutionary genomics, researchers have taken an interest in identifying substitutions that subtend convergent phenotypic adaptations. This is a difficult question that requires distinguishing foreground convergent substitutions that are involved in the convergent phenotype from background convergent substitutions. Those may be linked to other adaptations, may be neutral or may be the consequence of mutational biases. Furthermore, there is no generally accepted definition of convergent substitutions. Various methods that use different definitions have been proposed in the literature, resulting in different sets of candidate foreground convergent substitutions. In this article, we first describe the processes that can generate foreground convergent substitutions in coding sequences, separating adaptive from non-adaptive processes. Second, we review methods that have been proposed to detect foreground convergent substitutions in coding sequences and expose the assumptions that underlie them. Finally, we examine their power on simulations of convergent changes-including in the presence of a change in the efficacy of selection-and on empirical alignments. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Developmental and comparative transcriptomic identification of iridophore contribution to white barring in clownfish.

Actinopterygian fishes harbor at least eight distinct pigment cell types, leading to a fascinating diversity of colors. Among this diversity, the cellular origin of the white color appears to be linked to several pigment cell types such as iridophores or leucophores. We used the clownfish Amphiprion ocellaris, which has a color pattern consisting of white bars over a darker body, to characterize the pigment cells that underlie the white hue. We observe by electron microscopy that cells in white bars are similar to iridophores. In addition, the transcriptomic signature of clownfish white bars exhibits similarities with that of zebrafish iridophores. We further show by pharmacological treatments that these cells are necessary for the white color. Among the top differentially expressed genes in white skin, we identified several genes (fhl2a, fhl2b, saiyan, gpnmb, and apoD1a) and show that three of them are expressed in iridophores. Finally, we show by CRISPR/Cas9 mutagenesis that these genes are critical for iridophore development in zebrafish. Our analyses provide clues to the genomic underpinning of color diversity and allow identification of new iridophore genes in fish.

CAARS: comparative assembly and annotation of RNA-Seq data.

MOTIVATION: RNA sequencing (RNA-Seq) is a widely used approach to obtain transcript sequences in non-model organisms, notably for performing comparative analyses. However, current bioinformatic pipelines do not take full advantage of pre-existing reference data in related species for improving RNA-Seq assembly, annotation and gene family reconstruction. RESULTS: We built an automated pipeline named CAARS to combine novel data from RNA-Seq experiments with existing multi-species gene family alignments. RNA-Seq reads are assembled into transcripts by both de novo and assisted assemblies. Then, CAARS incorporates transcripts into gene families, builds gene alignments and trees and uses phylogenetic information to classify the genes as orthologs and paralogs of existing genes. We used CAARS to assemble and annotate RNA-Seq data in rodents and fishes using distantly related genomes as reference, a difficult case for this kind of analysis. We showed CAARS assemblies are more complete and accurate than those assembled by a standard pipeline consisting of de novo assembly coupled with annotation by sequence similarity on a guide species. In addition to annotated transcripts, CAARS provides gene family alignments and trees, annotated with orthology relationships, directly usable for downstream comparative analyses. AVAILABILITY AND IMPLEMENTATION: CAARS is implemented in Python and Ocaml and is freely available at https://github.com/carinerey/caars. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Accurate Detection of Convergent Amino-Acid Evolution with PCOC.

In the history of life, some phenotypes have been acquired several times independently, through convergent evolution. Recently, lots of genome-scale studies have been devoted to identify nucleotides or amino acids that changed in a convergent manner when the convergent phenotypes evolved. These efforts have had mixed results, probably because of differences in the detection methods, and because of conceptual differences about the definition of a convergent substitution. Some methods contend that substitutions are convergent only if they occur on all branches where the phenotype changed toward the exact same state at a given nucleotide or amino acid position. Others are much looser in their requirements and define a convergent substitution as one that leads the site at which they occur to prefer a phylogeny in which species with the convergent phenotype group together. Here, we suggest to look for convergent shifts in amino acid preferences instead of convergent substitutions to the exact same amino acid. We define as convergent shifts substitutions that occur on all branches where the phenotype changed and such that they correspond to a change in the type of amino acid preferred at this position. We implement the corresponding model into a method named PCOC. We show on simulations that PCOC better recovers convergent shifts than existing methods in terms of sensitivity and specificity. We test it on a plant protein alignment where convergent evolution has been studied in detail and find that our method recovers several previously identified convergent substitutions and proposes credible new candidates.

Evolution of mitotic spindle behavior during the first asymmetric embryonic division of nematodes.

Asymmetric cell division is essential to generate cellular diversity. In many animal cells, the cleavage plane lies perpendicular to the mitotic spindle, and it is the spindle positioning that dictates the size of the daughter cells. Although some properties of spindle positioning are conserved between distantly related model species and different cell types, little is known of the evolutionary robustness of the mechanisms underlying this event. We recorded the first embryonic division of 42 species of nematodes closely related to Caenorhabditis elegans, which is an excellent model system to study the biophysical properties of asymmetric spindle positioning. Our recordings, corresponding to 128 strains from 27 Caenorhabditis and 15 non-Caenorhabditis species (accessible at http://www.ens-lyon.fr/LBMC/NematodeCell/videos/), constitute a powerful collection of subcellular phenotypes to study the evolution of various cellular processes across species. In the present work, we analyzed our collection to the study of asymmetric spindle positioning. Although all the strains underwent an asymmetric first cell division, they exhibited large intra- and inter-species variations in the degree of cell asymmetry and in several parameters controlling spindle movement, including spindle oscillation, elongation, and displacement. Notably, these parameters changed frequently during evolution with no apparent directionality in the species phylogeny, with the exception of spindle transverse oscillations, which were an evolutionary innovation at the base of the Caenorhabditis genus. These changes were also unrelated to evolutionary variations in embryo size. Importantly, spindle elongation, displacement, and oscillation each evolved independently. This finding contrasts starkly with expectations based on C. elegans studies and reveals previously unrecognized evolutionary changes in spindle mechanics. Collectively, these data demonstrate that, while the essential process of asymmetric cell division has been conserved over the course of nematode evolution, the underlying spindle movement parameters can combine in various ways. Like other developmental processes, asymmetric cell division is subject to system drift.

Recombination, meiotic expression and human codon usage.

Synonymous codon usage (SCU) varies widely among human genes. In particular, genes involved in different functional categories display a distinct codon usage, which was interpreted as evidence that SCU is adaptively constrained to optimize translation efficiency in distinct cellular states. We demonstrate here that SCU is not driven by constraints on tRNA abundance, but by large-scale variation in GC-content, caused by meiotic recombination, via the non-adaptive process of GC-biased gene conversion (gBGC). Expression in meiotic cells is associated with a strong decrease in recombination within genes. Differences in SCU among functional categories reflect differences in levels of meiotic transcription, which is linked to variation in recombination and therefore in gBGC. Overall, the gBGC model explains 70% of the variance in SCU among genes. We argue that the strong heterogeneity of SCU induced by gBGC in mammalian genomes precludes any optimization of the tRNA pool to the demand in codon usage.

Transcriptomic signatures shaped by cell proportions shed light on comparative developmental biology.

BACKGROUND: Comparative transcriptomics can answer many questions in developmental and evolutionary developmental biology. Most transcriptomic studies start by showing global patterns of variation in transcriptomes that differ between species or organs through developmental time. However, little is known about the kinds of expression differences that shape these patterns. RESULTS: We compared transcriptomes during the development of two morphologically distinct serial organs, the upper and lower first molars of the mouse. We found that these two types of teeth largely share the same gene expression dynamics but that three major transcriptomic signatures distinguish them, all of which are shaped by differences in the relative abundance of different cell types. First, lower/upper molar differences are maintained throughout morphogenesis and stem from differences in the relative abundance of mesenchyme and from constant differences in gene expression within tissues. Second, there are clear time-shift differences in the transcriptomes of the two molars related to cusp tissue abundance. Third, the transcriptomes differ most during early-mid crown morphogenesis, corresponding to exaggerated morphogenetic processes in the upper molar involving fewer mitotic cells but more migrating cells. From these findings, we formulate hypotheses about the mechanisms enabling the two molars to reach different phenotypes. We also successfully applied our approach to forelimb and hindlimb development. CONCLUSIONS: Gene expression in a complex tissue reflects not only transcriptional regulation but also abundance of different cell types. This knowledge provides valuable insights into the cellular processes underpinning differences in organ development. Our approach should be applicable to most comparative developmental contexts.

Diversity in Morphology and Locomotory Behavior Is Associated with Niche Expansion in the Semi-aquatic Bugs.

Acquisition of new ecological opportunities is a major driver of adaptation and species diversification [1-4]. However, how groups of organisms expand their habitat range is often unclear [3]. We study the Gerromorpha, a monophyletic group of heteropteran insects that occupy a large variety of water surface-associated niches, from small puddles to open oceans [5, 6]. Due to constraints related to fluid dynamics [7-9] and exposure to predation [5, 10], we hypothesize that selection will favor high speed of locomotion in the Gerromorpha that occupy water-air interface niches relative to the ancestral terrestrial life style. Through biomechanical assays and phylogenetic reconstruction, we show that only species that occupy water surface niches can generate high maximum speeds. Basally branching lineages with ancestral mode of locomotion, consisting of tripod gait, achieved increased speed on the water through increasing midleg length, stroke amplitude, and stroke frequency. Derived lineages evolved rowing as a novel mode of locomotion through simultaneous sculling motion almost exclusively of the midlegs. We demonstrate that this change in locomotory behavior significantly reduced the requirement for high stroke frequency and energy expenditure. Furthermore, we show how the evolution of rowing, by reducing stroke frequency, may have eliminated the constraint on body size, which may explain the evolution of larger Gerromorpha. This correlation between the diversity in locomotion behaviors and niche specialization suggests that changes in morphology and behavior may facilitate the invasion and diversification in novel environments.

The genome of the crustacean Parhyale hawaiensis, a model for animal development, regeneration, immunity and lignocellulose digestion.

The amphipod crustacean Parhyale hawaiensis is a blossoming model system for studies of developmental mechanisms and more recently regeneration. We have sequenced the genome allowing annotation of all key signaling pathways, transcription factors, and non-coding RNAs that will enhance ongoing functional studies. Parhyale is a member of the Malacostraca clade, which includes crustacean food crop species. We analysed the immunity related genes of Parhyale as an important comparative system for these species, where immunity related aquaculture problems have increased as farming has intensified. We also find that Parhyale and other species within Multicrustacea contain the enzyme sets necessary to perform lignocellulose digestion ('wood eating'), suggesting this ability may predate the diversification of this lineage. Our data provide an essential resource for further development of Parhyale as an experimental model. The first malacostracan genome will underpin ongoing comparative work in food crop species and research investigating lignocellulose as an energy source.

Tinkering signaling pathways by gain and loss of protein isoforms: the case of the EDA pathway regulator EDARADD.

BACKGROUND: Only a handful of signaling pathways are major actors of development and responsible for both the conservation and the diversification of animal morphologies. To explain this twofold nature, gene duplication and enhancer evolution were predominantly put forth as tinkering mechanisms whereas the evolution of alternative isoforms has been, so far, overlooked. We investigate here the role of gain and loss of isoforms using Edaradd, a gene of the Ecodysplasin pathway, implicated in morphological evolution. A previous study had suggested a scenario of isoform gain and loss with an alternative isoform (A) newly gained in mammals but secondarily lost in mouse lineage. RESULTS: For a comprehensive view of A and B Edaradd isoforms history during mammal evolution, we obtained sequences for both isoforms in representative mammals and performed in vitro translations to support functional predictions. We showed that the ancestral B isoform is well conserved, whereas the mammal-specific A isoform was lost at least 7 times independently in terminal lineages throughout mammal phylogeny. Then, to gain insights into the functional relevance of this evolutionary pattern, we compared the biological function of these isoforms: i) In cellulo promoter assays showed that they are transcribed from two alternative promoters, only B exhibiting feedback regulation. ii) RT-PCR in various tissues and ENCODE data suggested that B isoform is systematically expressed whereas A isoform showed a more tissue-specific expression. iii) Both isoforms activated the NF-κB pathway in an in cellulo reporter assay, albeit at different levels and with different dynamics since A isoform exhibited feedback regulation at the protein level. Finally, only B isoform could rescue a zebrafish edaradd knockdown. CONCLUSIONS: These results suggest that the newly evolved A isoform enables modulating EDA signaling in specific conditions and with different dynamics. We speculate that during mammal diversification, A isoform regulation may have evolved rapidly, accompanying and possibly supporting the diversity of ectodermal appendages, while B isoform may have ensured essential roles. This study makes the case to pay greater attention to mosaic loss of evolutionarily speaking "young" isoforms as an important mechanism underlying phenotypic diversity and not simply as a manifestation of neutral evolution.

RAR/RXR binding dynamics distinguish pluripotency from differentiation associated cis-regulatory elements.

In mouse embryonic cells, ligand-activated retinoic acid receptors (RARs) play a key role in inhibiting pluripotency-maintaining genes and activating some major actors of cell differentiation. To investigate the mechanism underlying this dual regulation, we performed joint RAR/RXR ChIP-seq and mRNA-seq time series during the first 48 h of the RA-induced Primitive Endoderm (PrE) differentiation process in F9 embryonal carcinoma (EC) cells. We show here that this dual regulation is associated with RAR/RXR genomic redistribution during the differentiation process. In-depth analysis of RAR/RXR binding sites occupancy dynamics and composition show that in undifferentiated cells, RAR/RXR interact with genomic regions characterized by binding of pluripotency-associated factors and high prevalence of the non-canonical DR0-containing RA response element. By contrast, in differentiated cells, RAR/RXR bound regions are enriched in functional Sox17 binding sites and are characterized with a higher frequency of the canonical DR5 motif. Our data offer an unprecedentedly detailed view on the action of RA in triggering pluripotent cell differentiation and demonstrate that RAR/RXR action is mediated via two different sets of regulatory regions tightly associated with cell differentiation status.

Transcriptomics of developing embryos and organs: A raising tool for evo-devo.

Comparative transcriptomics has become an important tool for revisiting many evo-devo questions and exploring new ones, and its importance is likely to increase in the near future, partly because RNA-seq data open many new possibilities. The aim of this opinion piece is twofold. In the first section, we discuss the particularities of transcriptomic studies in evo-devo, focusing mainly on RNA-seq data. The preliminary processing steps (getting coding sequences as well as expression levels) are challenging, because many studied species do not have a sequenced genome. The next step (interpreting expression differences) is also challenging, due to several issues with interpreting expression levels in complex tissues, managing developmental stages and species heterochronies, and the problem of conceptualizing expression differences. In the second section, we discuss some past and possible future applications of transcriptomic approaches (using microarray or RNA-seq) to three major themes in evo-devo: the evolution of the developmental toolkit, the genetic and developmental basis for phenotypic changes, and the general rules of the evolution of development. We believe that conceptual and technical tools are necessary in order to fully exploit the richness of multispecies transcriptomic time-series data.

Programmed genome rearrangements: in lampreys, all cells are not equal.

How can organisms silence deleterious gene loci? A recent study has shed light on a very brute mechanism in a jawless vertebrate: the irreversible deletion of massive chunks of genomic DNA.

microRNA complements in deuterostomes: origin and evolution of microRNAs.

Although numerous studies have emphasized the role of microRNAs (miRNAs) in the control of many different cellular processes, they might also exert a profound effect on the macroevolution of animal body plans. It has been hypothesized that, because miRNAs increase genic precision and are continuously being added to metazoan genomes through geologic time, miRNAs might be instrumental for canalization of development and morphological evolution. Nonetheless, an outstanding question remains: how are new miRNAs constantly evolving? To address this question, we assessed the miRNA complements of four deuterostome species, chosen because of their sequenced genomes and well-resolved phylogeny. Our comparative analysis shows that each of these four species is characterized by a unique repertoire of miRNAs, with few instances of miRNA loss. Moreover, we find that almost half of the miRNAs identified in this study are located in intronic regions of protein coding genes, suggesting that new miRNAs might arise from intronic regions in a process we term intronic exaptation. We also show that miRNAs often occur within cotranscribed clusters, and describe the biological function of one of these conserved clusters, the miR-1/miR-133 cluster. Taken together, our work shows that miRNAs can easily emerge within already transcribed regions of DNA, whether it be introns or preexisting clusters of miRNAs and/or miRNAs and protein coding genes, and because of their regulatory roles, these novel players change the structure of gene regulatory networks, with potential macroevolutionary results.

microRNAs reveal the interrelationships of hagfish, lampreys, and gnathostomes and the nature of the ancestral vertebrate.

Hagfish and lampreys are the only living representatives of the jawless vertebrates (agnathans), and compared with jawed vertebrates (gnathostomes), they provide insight into the embryology, genomics, and body plan of the ancestral vertebrate. However, this insight has been obscured by controversy over their interrelationships. Morphological cladistic analyses have identified lampreys and gnathostomes as closest relatives, whereas molecular phylogenetic studies recover a monophyletic Cyclostomata (hagfish and lampreys as closest relatives). Here, we show through deep sequencing of small RNA libraries, coupled with genomic surveys, that Cyclostomata is monophyletic: hagfish and lampreys share 4 unique microRNA families, 15 unique paralogues of more primitive microRNA families, and 22 unique substitutions to the mature gene products. Reanalysis of morphological data reveals that support for cyclostome paraphyly was based largely on incorrect character coding, and a revised dataset is not decisive on the mono- vs. paraphyly of cyclostomes. Furthermore, we show fundamental conservation of microRNA expression patterns among lamprey, hagfish, and gnathostome organs, implying that the role of microRNAs within specific organs is coincident with their appearance within the genome and is conserved through time. Together, these data support the monophyly of cyclostomes and suggest that the last common ancestor of all living vertebrates was a more complex organism than conventionally accepted by comparative morphologists and developmental biologists.

Localized hypermutation and associated gene losses in legume chloroplast genomes.

Point mutations result from errors made during DNA replication or repair, so they are usually expected to be homogeneous across all regions of a genome. However, we have found a region of chloroplast DNA in plants related to sweetpea (Lathyrus) whose local point mutation rate is at least 20 times higher than elsewhere in the same molecule. There are very few precedents for such heterogeneity in any genome, and we suspect that the hypermutable region may be subject to an unusual process such as repeated DNA breakage and repair. The region is 1.5 kb long and coincides with a gene, ycf4, whose rate of evolution has increased dramatically. The product of ycf4, a photosystem I assembly protein, is more divergent within the single genus Lathyrus than between cyanobacteria and other angiosperms. Moreover, ycf4 has been lost from the chloroplast genome in Lathyrus odoratus and separately in three other groups of legumes. Each of the four consecutive genes ycf4-psaI-accD-rps16 has been lost in at least one member of the legume "inverted repeat loss" clade, despite the rarity of chloroplast gene losses in angiosperms. We established that accD has relocated to the nucleus in Trifolium species, but were unable to find nuclear copies of ycf4 or psaI in Lathyrus. Our results suggest that, as well as accelerating sequence evolution, localized hypermutation has contributed to the phenomenon of gene loss or relocation to the nucleus.