Character trees from transcriptome data: Origin and individuation of morphological characters and the so-called "species signal"

Comparative transcriptomics in serial organs uncovers early and pan-organ developmental changes associated with organ-specific morphological adaptation

Mice have evolved a new dental plan with two additional cusps on the upper molar, while hamsters were retaining the ancestral plan. By comparing the dynamics of molar development with transcriptome time series, we found at least three early changes in mouse upper molar development. Together, they redirect spatio-temporal dynamics to ultimately form two additional cusps. The mouse lower molar has undergone much more limited phenotypic evolution. Nevertheless, its developmental trajectory evolved as much as that of the upper molar and co-evolved with it. Among the coevolving changes, some are clearly involved in the new upper molar phenotype. We found a similar level of coevolution in bat limbs. In conclusion, our study reveals how serial organ morphology has adapted through organ-specific developmental changes, as expected, but also through shared changes that have organ-specific effects on the final phenotype. This highlights the important role of developmental system drift in one organ to accommodate adaptation in another.

Cell type and cell signaling innovations underlying mammalian pregnancy

How fetal and maternal cell types have co-evolved to enable mammalian placentation poses a unique evolutionary puzzle. Here, we present a multi-species atlas integrating single-cell transcriptomes from six species bracketing therian mammal diversity. We find that invasive trophoblasts share a gene-expression signature across eutherians, and evidence that endocrine decidual cells evolved stepwise from an immunomodulatory cell type retained in Tenrec with affinity to human decidua of menstruation. We recover evolutionary patterns in ligand-receptor signaling: fetal and maternal cells show a pronounced tendency towards disambiguation, but a predicted arms race dynamic between them is limited. We reconstruct cell communication networks of extinct mammalian ancestors, finding strong integration of fetal trophoblast into maternal networks. Together, our results reveal a dynamic history of cell type and signaling evolution.

Synopsis

The fetal-maternal interface is one of the most intense loci of cell-cell signaling in the human body. Invasion of cells from the fetal placenta into the uterus, and the corresponding transformation of maternal tissues called decidualization, first evolved in the stem lineage of eutherian mammals( 1 , 2 ). Single-cell studies of the human fetal-maternal interface have provided new insight into the cell type diversity and cell-cell interactions governing this chimeric organ( 3-5 ). However, the fetal-maternal interface is also one of the most rapidly evolving, and hence most diverse, characters among mammals( 6 ), and an evolutionary analysis is missing. Here, we present and compare single-cell data from the fetal-maternal interface of species bracketing key events in mammal phylogeny: a marsupial (opossum, Monodelphis domestica ), the afrotherian Tenrec ecaudatus, and four Euarchontoglires - guinea pig and mouse (Rodentia) together with recent macaque and human data (primates) ( 4 , 5 , 7 ). We infer cell type homologies, identify a gene-expression signature of eutherian invasive trophoblast conserved over 99 million years, and discover a predecidual cell in the tenrec which suggests stepwise evolution of the decidual stromal cell. We reconstruct ancestral cell signaling networks, revealing the integration of fetal cell types into the interface. Finally, we test two long-standing theoretical predictions, the disambiguation hypothesis( 8 ) and escalation hypothesis( 9 ), at transcriptome-wide scale, finding divergence between fetal and maternal signaling repertoires but arms race dynamics restricted to a small subset of ligand-receptor pairs. In so doing, we trace the co-evolutionary history of cell types and their signaling across mammalian viviparity.

Comparative chloroplast genomes study of five officinal Ardisia Species: Unraveling interspecific diversity and evolutionary insights in Ardisia

Ardisia S.W. (Primulaceae), naturally distributed in tropical and subtropical regions, has edible and medicinal values and is prevalent in clinical and daily use in China. More genetic information for distinct species delineation is needed to support the development and utilization of the genus Ardisia. We sequenced, annotated, and compared the chloroplast genomes of five Ardisia species: A. brunnescens, A. pusilla, A. squamulosa, A. crenata, and A. brevicaulis in this study. We found a typical quadripartite structure in all five chloroplast genomes, with lengths ranging from 155,045 to 156,943 bp. Except for A. pusilla, which lacked the ycf15 gene, the other four Ardisia species contained 114 unique genes, including 79 protein-coding genes, 30 tRNAs, and four rRNAs. In addition, the rps19 pseudogene gene was present only in A. brunnescens. Five highly variable DNA barcodes were identified for five Ardisia species, including trnT-GGU-psbD, trnT-UGU-trnL-UAA, rps4-trnT-UGU, rpl32-trnL-UAG, and rpoB-trnC-GAA. The RNA editiing sites of protein-coding genes in the five Ardisia plastome were characterized and compared, and 274 (A. crenata)-288 (A. brevicaulis) were found. The results of the phylogenetic analysis were consistent with the morphological classification. Sequence alignment and phylogenetic analysis showed that ycf15 genes were highly divergent in Primulaceae. Reconstructions of ancestral character states indicated that leaf margin morphology is critical for classifying the genus Ardisia, with a rodent-like character being the most primitive. These results provide valuable information on the taxonomy and evolution of Ardisia plants.

A tradeoff evolution between acoustic fat bodies and skull muscles in toothed whales

Toothed whales have developed specialized echolocation abilities that are crucial for underwater activities. Acoustic fat bodies, including the melon, extramandibular fat body, and intramandibular fat body, are vital for echolocation. This study explores the transcriptome of acoustic fat bodies in toothed whales, revealing some insight into their evolutionary origins and ecological significance. Comparative transcriptome analysis of acoustic fat bodies and related tissues in a harbor porpoise and a Pacific white-sided dolphin reveals that acoustic fat bodies possess characteristics of both muscle and adipose tissue, occupying an intermediate position. The melon and extramandibular fat body exhibit specific muscle-related functions, implying an evolutionary connection between acoustic fat bodies and muscle tissue. Furthermore, we suggested that the melon and extramandibular fat body originate from intramuscular adipose tissue, a component of white adipose tissue. The extramandibular fat body has been identified as an evolutionary homolog of the masseter muscle, supported by the specific expression of MYH16, a pivotal protein in masticatory muscles. The intramandibular fat body, located within the mandibular foramen, shows possibilities of the presence of several immune-related functions, likely due to its proximity to bone marrow. Furthermore, this study sheds light on leucine modification in the catabolic pathway, which leads to the accumulation of isovaleric acid in acoustic fat bodies. Swallowing without chewing, a major toothed whale feeding ecology adaptation, makes the masticatory muscle redundant and leads to the formation of the extramandibular fat body. We propose that the intramuscular fat enlargement in facial muscles, which influences acoustic fat body development, is potentially related to the substantial reorganization of head morphology in toothed whales during aquatic adaptation.

Experimental and phylogenetic evidence for correlated gene expression evolution in endometrial and skin fibroblasts

Gene expression change is a dominant mode of evolution. Mutations, however, can affect gene expression in multiple cell types. Therefore, gene expression evolution in one cell type can lead to similar gene expression changes in another cell type. Here, we test this hypothesis by investigating dermal skin fibroblasts (SFs) and uterine endometrial stromal fibroblasts (ESFs). The comparative dataset consists of transcriptomes from cultured SF and ESF of nine mammalian species. We find that evolutionary changes in gene expression in SF and ESF are highly correlated. The experimental dataset derives from a SCID mouse strain selected for slow cancer growth leading to substantial gene expression changes in SFs. We compared the gene expression profiles of SF with that of ESF and found a significant correlation between them. We discuss the implications of these findings for the evolutionary correlation between placental invasiveness and vulnerability to metastatic cancer.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Rebuilding ships while at sea-Character individuality, homology, and evolutionary innovation

How novel traits originate in evolution is still one of the most perplexing questions in Evolutionary Biology. Building on a previous account of evolutionary innovation, I here propose that evolutionary novelties are those individualized characters that are not homologous to any characters in the ancestor. To clarify this definition, I here provide a detailed analysis of the concepts of "character individuality" and "homology" first, before addressing their role for our understanding of evolutionary innovation. I will argue (1) that functional as well as structural considerations are important for character individualization; and (2) that compositional (structural) and positional homology need to be clearly distinguished to properly describe the evolutionary transformations of hierarchically structured characters. My account will therefore integrate functional and structural perspectives and put forward a new multi-level view of character identity and transformation.

Cell-type profiling in salamanders identifies innovations in vertebrate forebrain evolution

The evolution of advanced cognition in vertebrates is associated with two independent innovations in the forebrain: the six-layered neocortex in mammals and the dorsal ventricular ridge (DVR) in sauropsids (reptiles and birds). How these innovations arose in vertebrate ancestors remains unclear. To reconstruct forebrain evolution in tetrapods, we built a cell-type atlas of the telencephalon of the salamander Pleurodeles waltl. Our molecular, developmental, and connectivity data indicate that parts of the sauropsid DVR trace back to tetrapod ancestors. By contrast, the salamander dorsal pallium is devoid of cellular and molecular characteristics of the mammalian neocortex yet shares similarities with the entorhinal cortex and subiculum. Our findings chart the series of innovations that resulted in the emergence of the mammalian six-layered neocortex and the sauropsid DVR.

Cell types as species: Exploring a metaphor

The concept of "cell type," though fundamental to cell biology, is controversial. Cells have historically been classified into types based on morphology, physiology, or location. More recently, single cell transcriptomic studies have revealed fine-scale differences among cells with similar gross phenotypes. Transcriptomic snapshots of cells at various stages of differentiation, and of cells under different physiological conditions, have shown that in many cases variation is more continuous than discrete, raising questions about the relationship between cell type and cell state. Some researchers have rejected the notion of fixed types altogether. Throughout the history of discussions on cell type, cell biologists have compared the problem of defining cell type with the interminable and often contentious debate over the definition of arguably the most important concept in systematics and evolutionary biology, "species." In the last decades, systematics, like cell biology, has been transformed by the increasing availability of molecular data, and the fine-grained resolution of genetic relationships have generated new ideas about how that variation should be classified. There are numerous parallels between the two fields that make exploration of the "cell types as species" metaphor timely. These parallels begin with philosophy, with discussion of both cell types and species as being either individuals, groups, or something in between (e.g., homeostatic property clusters). In each field there are various different types of lineages that form trees or networks that can (and in some cases do) provide criteria for grouping. Developing and refining models for evolutionary divergence of species and for cell type differentiation are parallel goals of the two fields. The goal of this essay is to highlight such parallels with the hope of inspiring biologists in both fields to look for new solutions to similar problems outside of their own field.

A novel regulatory gene promotes novel cell fate by suppressing ancestral fate in the sea anemone Nematostella vectensis

In this study, we demonstrate how a new cell type can arise through duplication of an ancestral cell type followed by functional divergence of the new daughter cell. Specifically, we show that stinging cells in a cnidarian (namely, a sea anemone) emerged by duplication of an ancestral neuron followed by inhibition of the RFamide neuropeptide it once secreted. This finding is evidence that stinging cells evolved from a specific subtype of neurons and suggests other neuronal subtypes may have been coopted for other novel secretory functions.

Cnidocytes (i.e., stinging cells) are an unequivocally novel cell type used by cnidarians (i.e., corals, jellyfish, and their kin) to immobilize prey. Although they are known to share a common evolutionary origin with neurons, the developmental program that promoted the emergence of cnidocyte fate is not known. Using functional genomics in the sea anemone, Nematostella vectensis, we show that cnidocytes develop by suppression of neural fate in a subset of neurons expressing RFamide. We further show that a single regulatory gene, a C₂H₂-type zinc finger transcription factor (ZNF845), coordinates both the gain of novel (cnidocyte-specific) traits and the inhibition of ancestral (neural) traits during cnidocyte development and that this gene arose by domain shuffling in the stem cnidarian. Thus, we report a mechanism by which a truly novel regulatory gene (ZNF845) promotes the development of a truly novel cell type (cnidocyte) through duplication of an ancestral cell lineage (neuron) and inhibition of its ancestral identity (RFamide).

Butterfly eyespots evolved via cooption of an ancestral gene-regulatory network that also patterns antennae, legs, and wings

Where do butterfly eyespots come from? One of the long-standing questions in the field of evolution concerns addressing where novel complex traits come from. Here we show that butterfly eyespots, a novel complex trait, likely originated from the redeployment of a preexisting gene-regulatory network regulating antennae, legs, and wings, to novel locations on the wing.

Butterfly eyespots are beautiful novel traits with an unknown developmental origin. Here we show that eyespots likely originated via cooption of parts of an ancestral appendage gene-regulatory network (GRN) to novel locations on the wing. Using comparative transcriptome analysis, we show that eyespots cluster most closely with antennae, relative to multiple other tissues. Furthermore, three genes essential for eyespot development, Distal-less (Dll), spalt (sal), and Antennapedia (Antp), share similar regulatory connections as those observed in the antennal GRN. CRISPR knockout of cis-regulatory elements (CREs) for Dll and sal led to the loss of eyespots, antennae, legs, and also wings, demonstrating that these CREs are highly pleiotropic. We conclude that eyespots likely reused an ancient GRN for their development, a network also previously implicated in the development of antennae, legs, and wings.

Single-nucleus transcriptomes reveal evolutionary and functional properties of cell types in the Drosophila accessory gland

Many traits responsible for male reproduction evolve quickly, including gene expression phenotypes in germline and somatic male reproductive tissues. Rapid male evolution in polyandrous species is thought to be driven by competition among males for fertilizations and conflicts between male and female fitness interests that manifest in postcopulatory phenotypes. In Drosophila, seminal fluid proteins secreted by three major cell types of the male accessory gland and ejaculatory duct are required for female sperm storage and use, and influence female postcopulatory traits. Recent work has shown that these cell types have overlapping but distinct effects on female postcopulatory biology, yet relatively little is known about their evolutionary properties. Here, we use single-nucleus RNA-Seq of the accessory gland and ejaculatory duct from Drosophila melanogaster and two closely related species to comprehensively describe the cell diversity of these tissues and their transcriptome evolution for the first time. We find that seminal fluid transcripts are strongly partitioned across the major cell types, and expression of many other genes additionally defines each cell type. We also report previously undocumented diversity in main cells. Transcriptome divergence was found to be heterogeneous across cell types and lineages, revealing a complex evolutionary process. Furthermore, protein adaptation varied across cell types, with potential consequences for our understanding of selection on male postcopulatory traits.

Evolution of higher mesenchymal CD44 expression in the human lineage

CD44 is an extracellular matrix receptor implicated in cancer progression. CD44 increases the invasibility of skin (SF) and endometrial stromal fibroblasts (ESF) by cancer and trophoblast cells. We reasoned that the evolution of CD44 expression can affect both, the fetal–maternal interaction through CD44 in ESF as well as vulnerability to malignant cancer through expression in SF. We studied the evolution of CD44 expression in mammalian SF and ESF and demonstrate that in the human lineage evolved higher CD44 expression. Isoform expression in cattle and human is very similar suggesting that differences in invasibility are not due to the nature of expressed isoforms. We then asked whether the concerted gene expression increase in both cell types is due to shared regulatory mechanisms or due to cell type-specific factors. Reporter gene experiments with cells and cis-regulatory elements from human and cattle show that the difference of CD44 expression is due to cis effects as well as cell type-specific trans effects. These results suggest that the concerted expression increase is likely due to selection acting on both cell types because the evolutionary change in cell type-specific factors requires selection on cell type-specific functions. This scenario implies that the malignancy enhancing effects of elevated CD44 expression in humans likely evolved as a side-effect of positive selection on a yet unidentified other function of CD44. A possible candidate is the anti-fibrotic effect of CD44 but there are no reliable data showing that humans and primates are less fibrotic than other mammals.

Assessing evolutionary and developmental transcriptome dynamics in homologous cell types

Background

During development, complex organ patterns emerge through the precise temporal and spatial specification of different cell types. On an evolutionary timescale, these patterns can change, resulting in morphological diversification. It is generally believed that homologous anatomical structures are built—largely—by homologous cell types. However, whether a common evolutionary origin of such cell types is always reflected in the conservation of their intrinsic transcriptional specification programs is less clear.

Results

Here, we developed a user‐friendly bioinformatics workflow to detect gene co‐expression modules and test for their conservation across developmental stages and species boundaries. Using a paradigm of morphological diversification, the tetrapod limb, and single‐cell RNA‐sequencing data from two distantly related species, chicken and mouse, we assessed the transcriptional dynamics of homologous cell types during embryonic patterning. With mouse limb data as reference, we identified 19 gene co‐expression modules with varying tissue or cell type‐restricted activities. Testing for co‐expression conservation revealed modules with high evolutionary turnover, while others seemed maintained—to different degrees, in module make‐up, density or connectivity—over developmental and evolutionary timescales.

Conclusions

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

Extracellular matrix gene expression signatures as cell type and cell state identifiers

Transcriptomic signatures based on cellular mRNA expression profiles can be used to categorize cell types and states. Yet whether different functional groups of genes perform better or worse in this process remains largely unexplored. Here we test the core matrisome - that is, all genes coding for structural proteins of the extracellular matrix - for its ability to delineate distinct cell types in embryonic single-cell RNA-sequencing (scRNA-seq) data. We show that even though expressed core matrisome genes correspond to less than 2% of an entire cellular transcriptome, their RNA expression levels suffice to recapitulate essential aspects of cell type-specific clustering. Notably, using scRNA-seq data from the embryonic limb, we demonstrate that core matrisome gene expression outperforms random gene subsets of similar sizes and can match and exceed the predictive power of transcription factors. While transcription factor signatures generally perform better in predicting cell types at early stages of chicken and mouse limb development, i.e., when cells are less differentiated, the information content of the core matrisome signature increases in more differentiated cells. Moreover, using cross-species analyses, we show that these cell type-specific signatures are evolutionarily conserved. Our findings suggest that each cell type produces its own unique extracellular matrix, or matreotype, which becomes progressively more refined and cell type-specific as embryonic tissues mature.

Evolutionary Cell Type Mapping with Single-Cell Genomics

A fundamental characteristic of animal multicellularity is the spatial coexistence of functionally specialized cell types that are all encoded by a single genome sequence. Cell type transcriptional programs are deployed and maintained by regulatory mechanisms that control the asymmetric, differential access to genomic information in each cell. This genome regulation ultimately results in specific cellular phenotypes. However, the emergence, diversity, and evolutionary dynamics of animal cell types remain almost completely unexplored beyond a few species. Single-cell genomics is emerging as a powerful tool to build comprehensive catalogs of cell types and their associated gene regulatory programs in non-traditional model species. We review the current state of sampling efforts across the animal tree of life and challenges ahead for the comparative study of cell type programs. We also discuss how the phylogenetic integration of cell atlases can lead to the development of models of cell type evolution and a phylogenetic taxonomy of cells.

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Evolution of Embryo Implantation Was Enabled by the Origin of Decidual Stromal Cells in Eutherian Mammals

Mammalian pregnancy evolved in the therian stem lineage, that is, before the common ancestor of marsupials and eutherian (placental) mammals. Ancestral therian pregnancy likely involved a brief phase of attachment between the fetal and maternal tissues followed by parturition-similar to the situation in most marsupials including the opossum. In all eutherians, however, embryo attachment is followed by implantation, allowing for a stable fetal-maternal interface and an extended gestation. Embryo attachment induces an attachment reaction in the uterus that is homologous to an inflammatory response. Here, we elucidate the evolutionary mechanism by which the ancestral inflammatory response was transformed into embryo implantation in the eutherian lineage. We performed a comparative uterine transcriptomic and immunohistochemical study of three eutherians, armadillo (Dasypus novemcinctus), hyrax (Procavia capensis), and rabbit (Oryctolagus cuniculus); and one marsupial, opossum (Monodelphis domestica). Our results suggest that in the eutherian lineage, the ancestral inflammatory response was domesticated by suppressing one of its modules detrimental to pregnancy, namely, neutrophil recruitment by cytokine IL17A. Further, we propose that this suppression was mediated by decidual stromal cells, a novel cell type in eutherian mammals. We tested a prediction of this model in vitro and showed that decidual stromal cells can suppress the production of IL17A from helper T cells. Together, these results provide a mechanistic understanding of early stages in the evolution of eutherian pregnancy.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

Comparative transcriptomics of primary cells in vertebrates

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Gene expression of functionally-related genes coevolves across fungal species: detecting coevolution of gene expression using phylogenetic comparative methods

BACKGROUND

Researchers often measure changes in gene expression across conditions to better understand the shared functional roles and regulatory mechanisms of different genes. Analogous to this is comparing gene expression across species, which can improve our understanding of the evolutionary processes shaping the evolution of both individual genes and functional pathways. One area of interest is determining genes showing signals of coevolution, which can also indicate potential functional similarity, analogous to co-expression analysis often performed across conditions for a single species. However, as with any trait, comparing gene expression across species can be confounded by the non-independence of species due to shared ancestry, making standard hypothesis testing inappropriate.

RESULTS

We compared RNA-Seq data across 18 fungal species using a multivariate Brownian Motion phylogenetic comparative method (PCM), which allowed us to quantify coevolution between protein pairs while directly accounting for the shared ancestry of the species. Our work indicates proteins which physically-interact show stronger signals of coevolution than randomly-generated pairs. Interactions with stronger empirical and computational evidence also showing stronger signals of coevolution. We examined the effects of number of protein interactions and gene expression levels on coevolution, finding both factors are overall poor predictors of the strength of coevolution between a protein pair. Simulations further demonstrate the potential issues of analyzing gene expression coevolution without accounting for shared ancestry in a standard hypothesis testing framework. Furthermore, our simulations indicate the use of a randomly-generated null distribution as a means of determining statistical significance for detecting coevolving genes with phylogenetically-uncorrected correlations, as has previously been done, is less accurate than PCMs, although is a significant improvement over standard hypothesis testing. These methods are further improved by using a phylogenetically-corrected correlation metric.

CONCLUSIONS

Our work highlights potential benefits of using PCMs to detect gene expression coevolution from high-throughput omics scale data. This framework can be built upon to investigate other evolutionary hypotheses, such as changes in transcription regulatory mechanisms across species.

Evolution of placental invasion and cancer metastasis are causally linked

Among mammals, placental invasion is correlated with vulnerability to malignancy. Animals with more invasive placentation (e.g. humans) are more vulnerable to malignancy. To explain this correlation, we propose the hypothesis of Evolved Levels of Invasibility: the evolution of invasibility of stromal tissue affects both, placental and cancer invasion. We provide evidence for this hypothesis using an in vitro model. We find that bovine endometrial and skin fibroblasts are more resistant to invasion than their human counterparts. Gene expression profiling identified genes with high expression in human but not in bovine fibroblasts. Knocking down a subset of them in human fibroblasts leads to stronger resistance to cancer cell invasion. Identifying the evolutionary determinants of stromal invasibility can provide significant insights to develop rational anti-metastatic therapeutics.

The Evolution of Gene Expression Underlying Vision Loss in Cave Animals

Dissecting the evolutionary genetic processes underlying eye reduction and vision loss in obligate cave-dwelling organisms has been a long-standing challenge in evolutionary biology. Independent vision loss events in related subterranean organisms can provide critical insight into these processes as well as into the nature of convergent loss of complex traits. Advances in evolutionary developmental biology have illuminated the significant role of heritable gene expression variation in the evolution of new forms. Here, we analyze gene expression variation in adult eye tissue across the freshwater crayfish, representing four independent vision-loss events in caves. Species and individual expression patterns cluster by eye function rather than phylogeny, suggesting convergence in transcriptome evolution in independently blind animals. However, this clustering is not greater than what is observed in surface species with conserved eye function after accounting for phylogenetic expectations. Modeling expression evolution suggests that there is a common increase in evolutionary rates in the blind lineages, consistent with a relaxation of selective constraint maintaining optimal expression levels. This is evidence for a repeated loss of expression constraint in the transcriptomes of blind animals and that convergence occurs via a similar trajectory through genetic drift.

Phototransduction Gene Expression and Evolution in Cave and Surface Crayfishes

In the absence of light in caves, animals have repeatedly evolved reduced eyes and visual systems. Whether the underlying genetic components remain intact in blind species remains unanswered across taxa. The freshwater crayfish have evolved to live in caves multiple times throughout their history; therefore, this system provides an opportunity to probe the genetic patterns and processes underlying repeated vision loss. Using transcriptomic data from the eyes of 14 species of cave and surface crayfishes, we identify the expression of 17 genes putatively related to visual phototransduction. We find a similarly complete repertoire of phototransduction gene families expressed in cave and surface species, but that the expression levels of those transcripts are consistently lower in cave species. We find statistical support for episodic positive selection, increased and decreased selection strength in caves, depending on the gene family. Analyses of gene expression evolution suggest convergent and possibly adaptive downregulation of these genes across eye-reduction events. Our results reveal a combination of evolutionary processes acting on the sequences and gene expression levels of vision-related genes underlying the loss of vision in caves.

Pervasive Correlated Evolution in Gene Expression Shapes Cell and Tissue Type Transcriptomes

The evolution and diversification of cell types is a key means by which animal complexity evolves. Recently, hierarchical clustering and phylogenetic methods have been applied to RNA-seq data to infer cell type evolutionary history and homology. A major challenge for interpreting this data is that cell type transcriptomes may not evolve independently due to correlated changes in gene expression. This nonindependence can arise for several reasons, such as common regulatory sequences for genes expressed in multiple tissues, that is, pleiotropic effects of mutations. We develop a model to estimate the level of correlated transcriptome evolution (LCE) and apply it to different data sets. The results reveal pervasive correlated transcriptome evolution among different cell and tissue types. In general, tissues related by morphology or developmental lineage exhibit higher LCE than more distantly related tissues. Analyzing new data collected from bird skin appendages suggests that LCE decreases with the phylogenetic age of tissues compared, with recently evolved tissues exhibiting the highest LCE. Furthermore, we show correlated evolution can alter patterns of hierarchical clustering, causing different tissue types from the same species to cluster together. To identify genes that most strongly contribute to the correlated evolution signal, we performed a gene-wise estimation of LCE on a data set with ten species. Removing genes with high LCE allows for accurate reconstruction of evolutionary relationships among tissue types. Our study provides a statistical method to measure and account for correlated gene expression evolution when interpreting comparative transcriptome data.

NeuroSystematics and Periodic System of Neurons: Model vs Reference Species at Single-Cell Resolution

There is more than one way to develop neuronal complexity, and animals frequently use different molecular toolkits to achieve similar functional outcomes (=convergent evolution). Neurons are different not only because they have different functions, but also because neurons and circuits have different genealogies, and perhaps independent origins at the broadest scale from ctenophores and cnidarians to cephalopods and primates. By combining modern phylogenomics, single-neuron sequencing (scRNA-seq), machine learning, single-cell proteomics, and metabolomic across Metazoa, it is possible to reconstruct the evolutionary histories of neurons tracing them to ancestral secretory cells. Comparative data suggest that neurons, and perhaps synapses, evolved at least 2-3 times (in ctenophore, cnidarian and bilateral lineages) during ∼600 million years of animal evolution. There were also several independent events of the nervous system centralization either from a common bilateral/cnidarian ancestor without the bona fide neurons or from the urbilaterian with diffuse, nerve-net type nervous system. From the evolutionary standpoint, (i) a neuron should be viewed as a functional rather than a genetic character, and (ii) any given neural system might be chimeric and composed of different cell lineages with distinct origins and evolutionary histories. The identification of distant neural homologies or examples of convergent evolution among 34 phyla will not only allow the reconstruction of neural systems' evolution but together with single-cell "omic" approaches the proposed synthesis would lead to the "Periodic System of Neurons" with predictive power for neuronal phenotypes and plasticity. Such a phylogenetic classification framework of Neuronal Systematics (NeuroSystematics) might be a conceptual analog of the Periodic System of Chemical Elements. scRNA-seq profiling of all neurons in an entire brain or Brain-seq is now fully achievable in many nontraditional reference species across the entire animal kingdom. Arguably, marine animals are the most suitable for the proposed tasks because the world oceans represent the greatest taxonomic and body-plan diversity.

The Origin of Animal Multicellularity and Cell Differentiation

Over 600 million years ago, animals evolved from a unicellular or colonial organism whose cell(s) captured bacteria with a collar complex, a flagellum surrounded by a microvillar collar. Using principles from evolutionary cell biology, we reason that the transition to multicellularity required modification of pre-existing mechanisms for extracellular matrix synthesis and cytokinesis. We discuss two hypotheses for the origin of animal cell types: division of labor from ancient plurifunctional cells and conversion of temporally alternating phenotypes into spatially juxtaposed cell types. Mechanistic studies in diverse animals and their relatives promise to deepen our understanding of animal origins and cell biology.

Brownian model of transcriptome evolution and phylogenetic network visualization between tissues

While phylogenetic analysis of transcriptomes of the same tissue is usually congruent with the species tree, the controversy emerges when multiple tissues are included, that is, whether species from the same tissue are clustered together, or different tissues from the same species are clustered together. Recent studies have suggested that phylogenetic network approach may shed some lights on our understanding of multi-tissue transcriptome evolution; yet the underlying evolutionary mechanism remains unclear. In this paper we develop a Brownian-based model of transcriptome evolution under the phylogenetic network that can statistically distinguish between the patterns of species-clustering and tissue-clustering. Our model can be used as a null hypothesis (neutral transcriptome evolution) for testing any correlation in tissue evolution, can be applied to cancer transcriptome evolution to study whether two tumors of an individual appeared independently or via metastasis, and can be useful to detect convergent evolution at the transcriptional level.

Deep homology in the age of next-generation sequencing

The principle of homology is central to conceptualizing the comparative aspects of morphological evolution. The distinctions between homologous or non-homologous structures have become blurred, however, as modern evolutionary developmental biology (evo-devo) has shown that novel features often result from modification of pre-existing developmental modules, rather than arising completely de novo. With this realization in mind, the term 'deep homology' was coined, in recognition of the remarkably conserved gene expression during the development of certain animal structures that would not be considered homologous by previous strict definitions. At its core, it can help to formulate an understanding of deeper layers of ontogenetic conservation for anatomical features that lack any clear phylogenetic continuity. Here, we review deep homology and related concepts in the context of a gene expression-based homology discussion. We then focus on how these conceptual frameworks have profited from the recent rise of high-throughput next-generation sequencing. These techniques have greatly expanded the range of organisms amenable to such studies. Moreover, they helped to elevate the traditional gene-by-gene comparison to a transcriptome-wide level. We will end with an outlook on the next challenges in the field and how technological advances might provide exciting new strategies to tackle these questions.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

The evolutionary origin of bilaterian smooth and striated myocytes

The dichotomy between smooth and striated myocytes is fundamental for bilaterian musculature, but its evolutionary origin is unsolved. In particular, interrelationships of visceral smooth muscles remain unclear. Absent in fly and nematode, they have not yet been characterized molecularly outside vertebrates. Here, we characterize expression profile, ultrastructure, contractility and innervation of the musculature in the marine annelid Platynereis dumerilii and identify smooth muscles around the midgut, hindgut and heart that resemble their vertebrate counterparts in molecular fingerprint, contraction speed and nervous control. Our data suggest that both visceral smooth and somatic striated myocytes were present in the protostome-deuterostome ancestor and that smooth myocytes later co-opted the striated contractile module repeatedly – for example, in vertebrate heart evolution. During these smooth-to-striated myocyte conversions, the core regulatory complex of transcription factors conveying myocyte identity remained unchanged, reflecting a general principle in cell type evolution.

DOI:http://dx.doi.org/10.7554/eLife.19607.001

The origin and evolution of cell types

Cell types are the basic building blocks of multicellular organisms and are extensively diversified in animals. Despite recent advances in characterizing cell types, classification schemes remain ambiguous. We propose an evolutionary definition of a cell type that allows cell types to be delineated and compared within and between species. Key to cell type identity are evolutionary changes in the 'core regulatory complex' (CoRC) of transcription factors, that make emergent sister cell types distinct, enable their independent evolution and regulate cell type-specific traits termed apomeres. We discuss the distinction between developmental and evolutionary lineages, and present a roadmap for future research.

The Transcriptomic Evolution of Mammalian Pregnancy: Gene Expression Innovations in Endometrial Stromal Fibroblasts

The endometrial stromal fibroblast (ESF) is a cell type present in the uterine lining of therian mammals. In the stem lineage of eutherian mammals, ESF acquired the ability to differentiate into decidual cells in order to allow embryo implantation. We call the latter cell type "neo-ESF" in contrast to "paleo-ESF" which is homologous to eutherian ESF but is not able to decidualize. In this study, we compare the transcriptomes of ESF from six therian species: Opossum (Monodelphis domestica; paleo-ESF), mink, rat, rabbit, human (all neo-ESF), and cow (secondarily nondecidualizing neo-ESF). We find evidence for strong stabilizing selection on transcriptome composition suggesting that the expression of approximately 5,600 genes is maintained by natural selection. The evolution of neo-ESF from paleo-ESF involved the following gene expression changes: Loss of expression of genes related to inflammation and immune response, lower expression of genes opposing tissue invasion, increased markers for proliferation as well as the recruitment of FOXM1, a key gene transiently expressed during decidualization. Signaling pathways also evolve rapidly and continue to evolve within eutherian lineages. In the bovine lineage, where invasiveness and decidualization were secondarily lost, we see a re-expression of genes found in opossum, most prominently WISP2, and a loss of gene expression related to angiogenesis. The data from this and previous studies support a scenario, where the proinflammatory paleo-ESF was reprogrammed to express anti-inflammatory genes in response to the inflammatory stimulus coming from the implanting conceptus and thus paving the way for extended, trans-cyclic gestation.

New features of desiccation tolerance in the lichen photobiont Trebouxia gelatinosa are revealed by a transcriptomic approach

Trebouxia is the most common lichen-forming genus of aero-terrestrial green algae and all its species are desiccation tolerant (DT). The molecular bases of this remarkable adaptation are, however, still largely unknown. We applied a transcriptomic approach to a common member of the genus, T. gelatinosa, to investigate the alteration of gene expression occurring after dehydration and subsequent rehydration in comparison to cells kept constantly hydrated. We sequenced, de novo assembled and annotated the transcriptome of axenically cultured T. gelatinosa by using Illumina sequencing technology. We tracked the expression profiles of over 13,000 protein-coding transcripts. During the dehydration/rehydration cycle c. 92 % of the total protein-coding transcripts displayed a stable expression, suggesting that the desiccation tolerance of T. gelatinosa mostly relies on constitutive mechanisms. Dehydration and rehydration affected mainly the gene expression for components of the photosynthetic apparatus, the ROS-scavenging system, Heat Shock Proteins, aquaporins, expansins, and desiccation related proteins (DRPs), which are highly diversified in T. gelatinosa, whereas Late Embryogenesis Abundant Proteins were not affected. Only some of these phenomena were previously observed in other DT green algae, bryophytes and resurrection plants, other traits being distinctive of T. gelatinosa, and perhaps related to its symbiotic lifestyle. Finally, the phylogenetic inference extended to DRPs of other chlorophytes, embryophytes and bacteria clearly pointed out that DRPs of chlorophytes are not orthologous to those of embryophytes: some of them were likely acquired through horizontal gene transfer from extremophile bacteria which live in symbiosis within the lichen thallus.

A robust (re-)annotation approach to generate unbiased mapping references for RNA-seq-based analyses of differential expression across closely related species

BACKGROUND

RNA-seq based on short reads generated by next generation sequencing technologies has become the main approach to study differential gene expression. Until now, the main applications of this technique have been to study the variation of gene expression in a whole organism, tissue or cell type under different conditions or at different developmental stages. However, RNA-seq also has a great potential to be used in evolutionary studies to investigate gene expression divergence in closely related species.

RESULTS

We show that the published genomes and annotations of the three closely related Drosophila species D. melanogaster, D. simulans and D. mauritiana have limitations for inter-specific gene expression studies. This is due to missing gene models in at least one of the genome annotations, unclear orthology assignments and significant gene length differences in the different species. A comprehensive evaluation of four statistical frameworks (DESeq2, DESeq2 with length correction, RPKM-limma and RPKM-voom-limma) shows that none of these methods sufficiently accounts for inter-specific gene length differences, which inevitably results in false positive candidate genes. We propose that published reference genomes should be re-annotated before using them as references for RNA-seq experiments to include as many genes as possible and to account for a potential length bias. We present a straight-forward reciprocal re-annotation pipeline that allows to reliably compare the expression for nearly all genes annotated in D. melanogaster.

CONCLUSIONS

We conclude that our reciprocal re-annotation of previously published genomes facilitates the analysis of significantly more genes in an inter-specific differential gene expression study. We propose that the established pipeline can easily be applied to re-annotate other genomes of closely related animals and plants to improve comparative expression analyses.

The Lungfish Transcriptome: A Glimpse into Molecular Evolution Events at the Transition from Water to Land

Lungfish and coelacanths are the only living sarcopterygian fish. The phylogenetic relationship of lungfish to the last common ancestor of tetrapods and their close morphological similarity to their fossil ancestors make this species uniquely interesting. However their genome size, the largest among vertebrates, is hampering the generation of a whole genome sequence. To provide a partial solution to the problem, a high-coverage lungfish reference transcriptome was generated and assembled. The present findings indicate that lungfish, not coelacanths, are the closest relatives to land-adapted vertebrates. Whereas protein-coding genes evolve at a very slow rate, possibly reflecting a “living fossil” status, transposable elements appear to be active and show high diversity, suggesting a role for them in the remarkable expansion of the lungfish genome. Analyses of single genes and gene families documented changes connected to the water to land transition and demonstrated the value of the lungfish reference transcriptome for comparative studies of vertebrate evolution.

Inferring cell type innovations by phylogenetic methods-concepts, methods, and limitations

Multicellular organisms are composed of distinct cell types that have specific roles in the body. Each cell type is a product of two kinds of historical processes-development and evolution. Although the concept of a cell type is difficult to define, the cell type concept based on the idea of the core regulatory network (CRN), a gene regulatory network that determines the identity of a cell type, illustrates the essential aspects of the cell type concept. The first step toward elucidating cell type evolution is to reconstruct the evolutionary relationships of cell types, or the cell type tree. The sister cell type model assumes that a new cell type evolves through divergence from a multifunctional ancestral cell type, creating tree-like evolutionary relationships between cell types. The process of generating a cell type tree can also be understood as the sequential addition of a new branching point on an ancestral cell differentiation hierarchy in evolution. A cell type tree thus represents an intertwined history of cell type evolution and development. Cell type trees can be reconstructed from high-throughput sequencing data, and the reconstruction of a cell type tree leads to the discovery of genes that are functionally important for a cell type. Although many issues including the lack of cross-species comparisons and the lack of a proper model for cell type evolution remain, the study of the origin of a new cell type using phylogenetic methods offers a promising new research avenue in developmental evolution. J. Exp. Zool. (Mol. Dev. Evol.) 324B: 653-661, 2015. © 2015 Wiley Periodicals, Inc.