Constrained vertebrate evolution by pleiotropic genes

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

Deciphering the origin of developmental stability: The role of intracellular expression variability in evolutionary conservation

Progress in evolutionary developmental biology (evo-devo) has deepened our understanding of how intrinsic properties of embryogenesis, along with natural selection and population genetics, shape phenotypic diversity. A focal point of recent empirical and theoretical research is the idea that highly developmentally stable phenotypes are more conserved in evolution. Previously, we demonstrated that in Japanese medaka (Oryzias latipes), embryonic stages and genes with high stability, estimated through whole-embryo RNA-seq, are highly conserved in subsequent generations. However, the precise origin of the stability of gene expression levels evaluated at the whole-embryo level remained unclear. Such stability could be attributed to two distinct sources: stable intracellular expression levels or spatially stable expression patterns. Here we demonstrate that stability observed in whole-embryo RNA-seq can be attributed to stability at the cellular level (low variability in gene expression at the cellular levels). We quantified the intercellular variations in expression levels and spatial gene expression patterns for seven key genes involved in patterning dorsoventral and rostrocaudal regions during early development in medaka. We evaluated intracellular variability by counting transcripts and found its significant correlation with variation observed in whole-embryo RNA-seq data. Conversely, variation in spatial gene expression patterns, assessed through intraindividual left-right asymmetry, showed no correlation. Given the previously reported correlation between stability and conservation of expression levels throughout embryogenesis, our findings suggest a potential general trend: the stability or instability of developmental systems-and the consequent evolutionary diversity-may be primarily anchored in intrinsic fundamental elements such as the variability of intracellular states.

The Phylotypic Brain of Vertebrates, from Neural Tube Closure to Brain Diversification

BACKGROUND

The phylotypic or intermediate stages are thought to be the most evolutionary conserved stages throughout embryonic development. The contrast with divergent early and later stages derived from the concept of the evo-devo hourglass model. Nonetheless, this developmental constraint has been studied as a whole embryo process, not at organ level. In this review, we explore brain development to assess the existence of an equivalent brain developmental hourglass. In the specific case of vertebrates, we propose to split the brain developmental stages into: (1) Early: Neurulation, when the neural tube arises after gastrulation. (2) Intermediate: Brain patterning and segmentation, when the neuromere identities are established. (3) Late: Neurogenesis and maturation, the stages when the neurons acquire their functionality. Moreover, we extend this analysis to other chordates brain development to unravel the evolutionary origin of this evo-devo constraint.

SUMMARY

Based on the existing literature, we hypothesise that a major conservation of the phylotypic brain might be due to the pleiotropy of the inductive regulatory networks, which are predominantly expressed at this stage. In turn, earlier stages such as neurulation are rather mechanical processes, whose regulatory networks seem to adapt to environment or maternal geometries. The later stages are also controlled by inductive regulatory networks, but their effector genes are mostly tissue-specific and functional, allowing diverse developmental programs to generate current brain diversity. Nonetheless, all stages of the hourglass are highly interconnected: divergent neurulation must have a vertebrate shared end product to reproduce the vertebrate phylotypic brain, and the boundaries and transcription factor code established during the highly conserved patterning will set the bauplan for the specialised and diversified adult brain.

KEY MESSAGES

The vertebrate brain is conserved at phylotypic stages, but the highly conserved mechanisms that occur during these brain mid-development stages (Inducing Regulatory Networks) are also present during other stages. Oppositely, other processes as cell interactions and functional neuronal genes are more diverse and majoritarian in early and late stages of development, respectively. These phenomena create an hourglass of transcriptomic diversity during embryonic development and evolution, with a really conserved bottleneck that set the bauplan for the adult brain around the phylotypic stage.

Developmental hourglass: Verification by numerical evolution and elucidation by dynamical-systems theory

Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.

Critical Evaluation of Embedding Media for Histological Studies of Early Stages of Chick Embryo Development

A histological examination is an important tool in embryology, developmental biology, and correlated areas. Despite the amount of information available about tissue embedding and different media, there is a lack of information regarding best practices for embryonic tissues. Embryonic tissues are considered fragile structures, usually small in size, and frequently challenging to position correctly in media for the subsequent histological steps. Here, we discuss the embedding media and procedures that provided us with appropriate preservation of tissue and easier orientation of embryos at early development. Fertilized Gallus gallus eggs were incubated for 72 h, collected, fixed, processed, and embedded with paraplast, polyethylene glycol (PEG), or historesin. These resins were compared by the precision of tissue orientation, the preview of the embryos in the blocks, microtomy, contrast in staining, preservation, average time, and cost. Paraplast and PEG did not allow correct embryo orientation, even with agar–gelatin pre-embedded samples. Additionally, structural maintenance was hindered and did not allow detailed morphological assessment, presenting tissue shrinkage and disruption. Historesin provided precise tissue orientation and excellent preservation of structures. Assessing the performance of the embedding media contributes significantly to future developmental research, optimizing the processing of embryo specimens and improving results.

Weak gene-gene interaction facilitates the evolution of gene expression plasticity

BACKGROUND

Individual organisms may exhibit phenotypic plasticity when they acclimate to different conditions. Such plastic responses may facilitate or constrain the adaptation of their descendant populations to new environments, complicating their evolutionary trajectories beyond the genetic blueprint. Intriguingly, phenotypic plasticity itself can evolve in terms of its direction and magnitude during adaptation. However, we know little about what determines the evolution of phenotypic plasticity, including gene expression plasticity. Recent laboratory-based studies suggest dominance of reversing gene expression plasticity-plastic responses that move the levels of gene expression away from the new optima. Nevertheless, evidence from natural populations is still limited.

RESULTS

Here, we studied gene expression plasticity and its evolution in the montane and lowland populations of an elevationally widespread songbird-the Rufous-capped Babbler (Cyanoderma ruficeps)-with reciprocal transplant experiments and transcriptomic analyses; we set common gardens at altitudes close to these populations' native ranges. We confirmed the prevalence of reversing plasticity in genes associated with altitudinal adaptation. Interestingly, we found a positive relationship between magnitude and degree of evolution in gene expression plasticity, which was pertinent to not only adaptation-associated genes but also the whole transcriptomes from multiple tissues. Furthermore, we revealed that genes with weaker expressional interactions with other genes tended to exhibit stronger plasticity and higher degree of plasticity evolution, which explains the positive magnitude-evolution relationship.

CONCLUSIONS

Our experimental evidence demonstrates that species may initiate their adaptation to new habitats with genes exhibiting strong expression plasticity. We also highlight the role of expression interdependence among genes in regulating the magnitude and evolution of expression plasticity. This study illuminates how the evolution of phenotypic plasticity in gene expression facilitates the adaptation of species to challenging environments in nature.

Stability in gene expression and body-plan development leads to evolutionary conservation

BACKGROUND

Phenotypic evolution is mainly explained by selection for phenotypic variation arising from factors including mutation and environmental noise. Recent theoretical and experimental studies have suggested that phenotypes with greater developmental stability tend to have a constant phenotype and gene expression level within a particular genetic and environmental condition, and this positively correlates with stronger evolutionary conservation, even after the accumulation of genetic changes. This could reflect a novel mechanism that contributes to evolutionary conservation; however, it remains unclear whether developmental stability is the cause, or whether at least it contributes to their evolutionary conservation. Here, using Japanese medaka lines, we tested experimentally whether developmental stages and gene expression levels with greater stability led to their evolutionary conservation.

RESULTS

We first measured the stability of each gene expression level and developmental stage (defined here as the whole embryonic transcriptome) in the inbred F0 medaka population. We then measured their evolutionary conservation in the F3 generation by crossing the F0 line with the distantly related Japanese medaka line (Teradomori), followed by two rounds of intra-generational crossings. The results indicated that the genes and developmental stages that had smaller variations in the F0 generation showed lower diversity in the hybrid F3 generation, which implies a causal relationship between stability and evolutionary conservation.

CONCLUSIONS

These findings suggest that the stability in phenotypes, including the developmental stages and gene expression levels, leads to their evolutionary conservation; this most likely occurs due to their low potential to generate phenotypic variation. In addition, since the highly stable developmental stages match with the body-plan-establishment stage, it also implies that the developmental stability potentially contributed to the strict conservation of animal body plan.

Unraveling signatures of chicken genetic diversity and divergent selection in breed-specific patterns of early myogenesis, nitric oxide metabolism and post-hatch growth

Introduction: Due to long-term domestication, breeding and divergent selection, a vast genetic diversity in poultry currently exists, with various breeds being characterized by unique phenotypic and genetic features. Assuming that differences between chicken breeds divergently selected for economically and culturally important traits manifest as early as possible in development and growth stages, we aimed to explore breed-specific patterns and interrelations of embryo myogenesis, nitric oxide (NO) metabolism and post-hatch growth rate (GR). Methods: These characteristics were explored in eight breeds of different utility types (meat-type, dual purpose, egg-type, game, and fancy) by incubating 70 fertile eggs per breed. To screen the differential expression of seven key myogenesis associated genes (MSTN, GHR, MEF2C, MYOD1, MYOG, MYH1, and MYF5), quantitative real-time PCR was used. Results: We found that myogenesis associated genes expressed in the breast and thigh muscles in a coordinated manner showing breed specificity as a genetic diversity signature among the breeds studied. Notably, coordinated ("accord") expression patterns of MSTN, GHR, and MEFC2 were observed both in the breast and thigh muscles. Also, associated expression vectors were identified for MYOG and MYOD1 in the breast muscles and for MYOG and MYF5 genes in the thigh muscles. Indices of NO oxidation and post-hatch growth were generally concordant with utility types of breeds, with meat-types breeds demonstrating higher NO oxidation levels and greater GR values as compared to egg-type, dual purpose, game and fancy breeds. Discussion: The results of this study suggest that differences in early myogenesis, NO metabolism and post-hatch growth are breed-specific; they appropriately reflect genetic diversity and accurately capture the evolutionary history of divergently selected chicken breeds.

Potential contribution of intrinsic developmental stability toward body plan conservation

Background

Despite the morphological diversity of animals, their basic anatomical patterns—the body plans in each animal phylum—have remained highly conserved over hundreds of millions of evolutionary years. This is attributed to conservation of the body plan-establishing developmental period (the phylotypic period) in each lineage. However, the evolutionary mechanism behind this phylotypic period conservation remains under debate. A variety of hypotheses based on the concept of modern synthesis have been proposed, such as negative selection in the phylotypic period through its vulnerability to embryonic lethality. Here we tested a new hypothesis that the phylotypic period is developmentally stable; it has less potential to produce phenotypic variations than the other stages, and this has most likely led to the evolutionary conservation of body plans.

Results

By analyzing the embryos of inbred Japanese medaka embryos raised under the same laboratory conditions and measuring the whole embryonic transcriptome as a phenotype, we found that the phylotypic period has greater developmental stability than other stages. Comparison of phenotypic differences between two wild medaka populations indicated that the phylotypic period and its genes in this period remained less variational, even after environmental and mutational modifications accumulated during intraspecies evolution. Genes with stable expression levels were enriched with those involved in cell-cell signalling and morphological specification such as Wnt and Hox, implying possible involvement in body plan development of these genes.

Conclusions

This study demonstrated the correspondence between the developmental stage with low potential to produce phenotypic variations and that with low diversity in micro- and macroevolution, namely the phylotypic period. Whereas modern synthesis explains evolution as a process of shaping of phenotypic variations caused by mutations, our results highlight the possibility that phenotypic variations are readily limited by the intrinsic nature of organisms, namely developmental stability, thus biasing evolutionary outcomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01276-5.

Distinguishing Evolutionary Conservation from Derivedness

While the concept of “evolutionary conservation” has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or “derivedness”. We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, “conservation” mainly considers genes or traits which species have in common, while “derivedness” additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, “transcriptomic derivedness index”, for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits.

The toggle switch model for gene expression change during the prenatal-to-postnatal transition in mammals

The prenatal-to-postnatal transition is a pivotal process in the life cycle whereby an organism shifts from responding to intrauterine cues to undergoing extrauterine stresses with many physiological adaptations. However, the molecular basis underlying the evolutionarily conserved physiological adaptations remains elusive. Here, we analyze the transcriptomes of seven organs across developmental time points from five mammalian species by constructing computational coexpression networks and report a developmental shift of gene expression at the perinatal stage. The low-to-high and high-to-low expressed genes tightly coalesce in the functional categories and gene regulatory pathways that implicate the physiological adaptions during the prenatal-to-postnatal transition, including lipid metabolism, circadian rhythm, immune response, cell cycle, and cell division. The low-to-high and high-to-low expressed genes around the perinatal stage tend to form the mutually inhibitory toggle switch gene pairs linking the gene regulatory networks in response to the environmental changes. We thus propose the toggle switch model for the developmental shift of gene expression as a mechanic framework to investigate how the physiological adaptations occur during the prenatal-to-postnatal transition.

Derivedness Index for Estimating Degree of Phenotypic Evolution of Embryos: A Study of Comparative Transcriptomic Analyses of Chordates and Echinoderms

Species retaining ancestral features, such as species called living fossils, are often regarded as less derived than their sister groups, but such discussions are usually based on qualitative enumeration of conserved traits. This approach creates a major barrier, especially when quantifying the degree of phenotypic evolution or degree of derivedness, since it focuses only on commonly shared traits, and newly acquired or lost traits are often overlooked. To provide a potential solution to this problem, especially for inter-species comparison of gene expression profiles, we propose a new method named "derivedness index" to quantify the degree of derivedness. In contrast to the conservation-based approach, which deals with expressions of commonly shared genes among species being compared, the derivedness index also considers those that were potentially lost or duplicated during evolution. By applying our method, we found that the gene expression profiles of penta-radial phases in echinoderm tended to be more highly derived than those of the bilateral phase. However, our results suggest that echinoderms may not have experienced much larger modifications to their developmental systems than chordates, at least at the transcriptomic level. In vertebrates, we found that the mid-embryonic and organogenesis stages were generally less derived than the earlier or later stages, indicating that the conserved phylotypic period is also less derived. We also found genes that potentially explain less derivedness, such as Hox genes. Finally, we highlight technical concerns that may influence the measured transcriptomic derivedness, such as read depth and library preparation protocols, for further improvement of our method through future studies. We anticipate that this index will serve as a quantitative guide in the search for constrained developmental phases or processes.

From Cell Types to an Integrated Understanding of Brain Evolution: The Case of the Cerebral Cortex

With the discovery of the incredible diversity of neurons, Cajal and coworkers laid the foundation of modern neuroscience. Neuron types are not only structural units of nervous systems but also evolutionary units, because their identities are encoded in the genome. With the advent of high-throughput cellular transcriptomics, neuronal identities can be characterized and compared systematically across species. The comparison of neurons in mammals, reptiles, and birds indicates that the mammalian cerebral cortex is a mosaic of deeply conserved and recently evolved neuron types. Using the cerebral cortex as a case study, this review illustrates how comparing neuron types across species is key to reconciling observations on neural development, neuroanatomy, circuit wiring, and physiology for an integrated understanding of brain evolution.

New Genes Interacted with Recent Whole Genome Duplicates in the Fast Stem Growth of Bamboos

As drivers of evolutionary innovations, new genes allow organisms to explore new niches. However, clear examples of this process remain scarce. Bamboos, the unique grass lineage diversifying into the forest, have evolved with a key innovation of fast growth of woody stem, reaching up to 1 m/day. Here, we identify 1,622 bamboo-specific orphan genes that appeared in recent 46 million years, and 19 of them evolved from noncoding ancestral sequences with entire de novo origination process reconstructed. The new genes evolved gradually in exon−intron structure, protein length, expression specificity, and evolutionary constraint. These new genes, whether or not from de novo origination, are dominantly expressed in the rapidly developing shoots, and make transcriptomes of shoots the youngest among various bamboo tissues, rather than reproductive tissue in other plants. Additionally, the particularity of bamboo shoots has also been shaped by recent whole-genome duplicates (WGDs), which evolved divergent expression patterns from ancestral states. New genes and WGDs have been evolutionarily recruited into coexpression networks to underline fast-growing trait of bamboo shoot. Our study highlights the importance of interactions between new genes and genome duplicates in generating morphological innovation.

The developmental hourglass model is applicable to the spinal cord based on single-cell transcriptomes and non-conserved cis-regulatory elements

The developmental hourglass model predicts that embryonic morphology is most conserved at the mid‐embryonic stage and diverges at the early and late stages. To date, this model has been verified by examining the anatomical features or gene expression profiles at the whole embryonic level. Here, by data mining approach utilizing multiple genomic and transcriptomic datasets from different species in combination, and by experimental validation, we demonstrate that the hourglass model is also applicable to a reduced element, the spinal cord. In the middle of spinal cord development, dorsoventrally arrayed neuronal progenitor domains are established, which are conserved among vertebrates. By comparing the publicly available single‐cell transcriptome datasets of mice and zebrafish, we found that ventral subpopulations of post‐mitotic spinal neurons display divergent molecular profiles. We also detected the non‐conservation of cis‐regulatory elements located around the progenitor fate determinants, indicating that the cis‐regulatory elements contributing to the progenitor specification are evolvable. These results demonstrate that, despite the conservation of the progenitor domains, the processes before and after the progenitor domain specification diverged. This study will be helpful to understand the molecular basis of the developmental hourglass model.

The hourglass model of evolutionary conservation during embryogenesis extends to developmental enhancers with signatures of positive selection

Inter-species comparisons of both morphology and gene expression within a phylum have revealed a period in the middle of embryogenesis with more similarity between species compared to earlier and later time-points. This "developmental hourglass" pattern has been observed in many phyla, yet the evolutionary constraints on gene expression, and underlying mechanisms of how this is regulated, remains elusive. Moreover, the role of positive selection on gene regulation in the more diverged earlier and later stages of embryogenesis remains unknown. Here, using DNase-seq to identify regulatory regions in two distant Drosophila species (D. melanogaster and D. virilis), we assessed the evolutionary conservation and adaptive evolution of enhancers throughout multiple stages of embryogenesis. This revealed a higher proportion of conserved enhancers at the phylotypic period, providing a regulatory basis for the hourglass expression pattern. Using an in silico mutagenesis approach, we detect signatures of positive selection on developmental enhancers at early and late stages of embryogenesis, with a depletion at the phylotypic period, suggesting positive selection as one evolutionary mechanism underlying the hourglass pattern of animal evolution.

Osteological characters of birds and reptiles are more congruent with molecular phylogenies than soft characters are

Despite increased use of genomic data in phylogenetics, morphological information remains vital for resolving evolutionary relationships, particularly for fossil taxa. The properties and models of evolution of molecular sequence data are well characterized and mature, relative to those of morphological data. Furthermore, heterogeneity, integration and relative homoplasy of empirical morphological data could prove problematic for phylogenetic reconstruction. Here we compare osteological and non-osteological characters of 28 morphological datasets of extant saurians in terms of their homoplasy relative to molecular trees. Analysis of individual avian datasets finds osteological characters to be significantly more consistent with molecular data than soft characters are. Significant differences between morphological partitions were also observed in the age at which characters resolved on molecular trees. Osteological character changes occur relatively earlier in deep branches, whilst soft-tissue character transitions are more recent in shallow branches. The combined results demonstrate differences in evolutionary dynamics between morphological partitions. This may reflect evolutionary constraints acting on osteological characters, compared with the relative lability of soft characters. Furthermore, it provides some support to phylogenetic interpretations of fossil data, including dinosaurs, which are predominately osteological. Recent advances in amphibian and mammal phylogenetics may make these patterns possible to test for all tetrapods.

Comparative Transcriptomics Reveals Distinct Patterns of Gene Expression Conservation through Vertebrate Embryogenesis

Despite life's diversity, studies of variation often remind us of our shared evolutionary past. Abundant genome sequencing and analyses of gene regulatory networks illustrate that genes and entire pathways are conserved, reused, and elaborated in the evolution of diversity. Predating these discoveries, 19th-century embryologists observed that though morphology at birth varies tremendously, certain stages of vertebrate embryogenesis appear remarkably similar across vertebrates. In the mid to late 20th century, anatomical variability of early and late-stage embryos and conservation of mid-stages embryos (the "phylotypic" stage) was named the hourglass model of diversification. This model has found mixed support in recent analyses comparing gene expression across species possibly owing to differences in species, embryonic stages, and gene sets compared. We compare 186 microarray and RNA-seq data sets covering embryogenesis in six vertebrate species. We use an unbiased clustering approach to group stages of embryogenesis by transcriptomic similarity and ask whether gene expression similarity of clustered embryonic stages deviates from a null expectation. We characterize expression conservation patterns of each gene at each evolutionary node after correcting for phylogenetic nonindependence. We find significant enrichment of genes exhibiting early conservation, hourglass, late conservation patterns in both microarray and RNA-seq data sets. Enrichment of genes showing patterned conservation through embryogenesis indicates diversification of embryogenesis may be temporally constrained. However, the circumstances under which each pattern emerges remain unknown and require both broad evolutionary sampling and systematic examination of embryogenesis across species.

Cross-species comparisons and in vitro models to study tempo in development and homeostasis

Time is inherent to biological processes. It determines the order of events and the speed at which they take place. However, we still need to refine approaches to measure the course of time in biological systems and understand what controls the pace of development. Here, we argue that the comparison of biological processes across species provides molecular insight into the timekeeping mechanisms in biology. We discuss recent findings and the open questions in the field and highlight the use of in vitro systems as tools to investigate cell-autonomous control as well as the coordination of temporal mechanisms within tissues. Further, we discuss the relevance of studying tempo for tissue transplantation, homeostasis and lifespan.

Polygenic Adaptation: Integrating Population Genetics and Gene Regulatory Networks

The adaptation of populations to local environments often relies on the selection of optimal values for polygenic traits. Here, we first summarize the results obtained from different quantitative genetics and population genetics models, about the genetic architecture of polygenic traits and their response to directional selection. We then highlight the contribution of systems biology to the understanding of the molecular bases of polygenic traits and the evolution of gene regulatory networks involved in these traits. Finally, we discuss the need for a unifying framework merging the fields of population genetics, quantitative genetics and systems biology to better understand the molecular bases of polygenic traits adaptation.

On the sequence heterochrony of cranial ossification of bats in light of Haeckel's recapitulation theory

Haeckel's recapitulation theory has been a controversial topic in evolutionary biology. However, we have seen some recent cases applying Haeckel's view to interpret the interspecific variation of prenatal ontogeny. To revisit the validity of Haeckel's recapitulation theory, we take bats that have undergone drastic morphological changes and possess a characteristic ecology as a case study. All members of Rhinolophoidea and Yangochiroptera can generate an ultrasonic pulse from the larynx to interpret surrounding objects (laryngeal echolocation) whereas Pteropodidae lacks such ability. It is known that the petrosal bone is particularly derived in shape and expanded in laryngeal echolocators. If Haeckel's recapitulation theory holds, the formation of this derived trait should occur later than those of other bones. Therefore, we compared the prenatal ossification timing of the petrosal in 15 bat species and five outgroup species. We found that the ossification of the petrosal is accelerated in laryngeal echolocators while it is the last bone to ossify in non-laryngeal echolocating bats and non-volant mammals, which runs counter to the prediction generated by Haeckel's recapitulation theory. We point out the evolutionarily labile nature of trait developmental timing and emphasize that Haeckel's recapitulation theory does not hold in many cases. We caution that generating predictions on ancestral conditions and evolutionary history leading from Haeckel's recapitulation theory is not well supported.

Measuring potential effects of the developmental burden associated with the vertebrate notochord

The notochord functions primarily as a supporting tissue to maintain the anteroposterior axis of primitive chordates, a function that is replaced entirely by the vertebral column in many vertebrates. The notochord still appears during vertebrate embryogenesis and plays a crucial role in the developmental pattern formation of surrounding structures, such as the somites and neural tube, providing the basis for the vertebrate body plan. The indispensable role of the notochord has often been referred to as the developmental burden and used to explain the evolutionary conservation of notochord; however, the existence of this burden has not been successfully exemplified so far. Since the adaptive value of target tissues appears to result in the evolutionary conservation of upstream structures through the developmental burden, we performed comparative gene expression profiling of the notochord, somites, and neural tube during the mid‐embryonic stages in turtles and chicken to measure their evolutionary conservation. When compared with the somites and neural tube, overall gene expression profiles in the notochord showed significantly lower or merely comparable levels of conservation. However, genes involved in inductive signalings, such as the sonic hedgehog (Shh) cascade and the formation of functional primary cilia, showed relatively higher levels of conservation in all the three structures analyzed. Collectively, these results suggest that shh signals are critical as the inductive source and receiving structures, possibly constituting the inter‐dependencies of developmental burden.

Potential evolutionary effects toward notochord by developmental burden was evaluated by Laser Micro Dissection RNAseq (LMDseq). Notochord was less conserved than neural tube and somites; however, genes in sonic hedgehog (shh) signaling cascade was found to be evolutionarily conserved (not only in notochord but also in somites and neural tube). These results suggest that Shh signals are critical as the inductive source and receiving structures, possibly constituting the inter‐dependencies of developmental burden. Further studies that directly measure the burden required to verify the hypothesis are awaited.

Embryonic evidence uncovers convergent origins of laryngeal echolocation in bats

Bats are the second-most speciose group of mammals, comprising 20% of species diversity today. Their global explosion, representing one of the greatest adaptive radiations in mammalian history, is largely attributed to their ability of laryngeal echolocation and powered flight, which enabled them to conquer the night sky, a vast and hitherto unoccupied ecological niche. While there is consensus that powered flight evolved only once in the lineage, whether laryngeal echolocation has a single origin in bats or evolved multiple times independently remains disputed. Here, we present developmental evidence in support of laryngeal echolocation having multiple origins in bats. This is consistent with a non-echolocating bat ancestor and independent gain of echolocation in Yinpterochiroptera and Yangochiroptera, as well as the gain of primitive echolocation in the bat ancestor, followed by convergent evolution of laryngeal echolocation in Yinpterochiroptera and Yangochiroptera, with loss of primitive echolocation in pteropodids. Our comparative embryological investigations found that there is no developmental difference in the hearing apparatus between non-laryngeal echolocating bats (pteropodids) and terrestrial non-bat mammals. In contrast, the echolocation system is developed heterotopically and heterochronically in the two phylogenetically distant laryngeal echolocating bats (rhinolophoids and yangochiropterans), providing the first embryological evidence that the echolocation system evolved independently in these bats.

Developmental hourglass and heterochronic shifts in fin and limb development

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as “developmental constraint”. These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved.

The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb.

To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint.

Darwin’s theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

The developmental hourglass model and recapitulation: An attempt to integrate the two models

Recapitulation is a hypothetical concept that assumes embryogenesis of an animal parallels its own phylogenetic history, sequentially developing from more ancestral features to more derived ones. This concept predicts that the earliest developmental stage of various animals should represent the most evolutionarily conserved patterns. Recent transcriptome‐based studies, on the other hand, have reported that mid‐embryonic, organogenetic periods show the highest level of conservation (the developmental hourglass model). This, however, does not rule out the possibility that recapitulation would still be detected after the mid‐embryonic period. In accordance with this, recapitulation‐like morphological features are enriched in late developmental stages. Moreover, our recent chromatin accessibility‐based study provided molecular evidence for recapitulation in the mid‐to‐late embryogenesis of vertebrates, as newly evolved gene regulatory elements tended to be activated at late embryonic stages. In this review, we revisit the recapitulation hypothesis, together with recent molecular‐based studies that support the developmental hourglass model. We contend that the recapitulation hypothesis does not entirely contradict the developmental hourglass model and that these two may even coexist in later embryonic stages of vertebrates. Finally, we review possible mechanisms underlying the recapitulation pattern of chromatin accessibility together with the hourglass‐like evolutionary conservation in vertebrate embryogenesis.

Recapitulation pattern has been reported for chromatin accessibility during the mid‐to‐late embryogenesis.

The observed recapitulation pattern and the developmental hourglass model may coexist.

The possible evolutionary mechanisms underlying tendencies of embryonic evolution were discussed.

Sexual and asexual development: two distinct programs producing the same tunicate

Colonial tunicates are the only chordate that possess two distinct developmental pathways to produce an adult body: either sexually through embryogenesis or asexually through a stem cell-mediated renewal termed blastogenesis. Using the colonial tunicate Botryllus schlosseri, we combine transcriptomics and microscopy to build an atlas of the molecular and morphological signatures at each developmental stage for both pathways. The general molecular profiles of these processes are largely distinct. However, the relative timing of organogenesis and ordering of tissue-specific gene expression are conserved. By comparing the developmental pathways of B. schlosseri with other chordates, we identify hundreds of putative transcription factors with conserved temporal expression. Our findings demonstrate that convergent morphology need not imply convergent molecular mechanisms but that it showcases the importance that tissue-specific stem cells and transcription factors play in producing the same mature body through different pathways.

MeDAS: a Metazoan Developmental Alternative Splicing database

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Embryo-Like Features in Developing Bacillus subtilis Biofilms

Correspondence between evolution and development has been discussed for more than two centuries. Recent work reveals that phylogeny-ontogeny correlations are indeed present in developmental transcriptomes of eukaryotic clades with complex multicellularity. Nevertheless, it has been largely ignored that the pervasive presence of phylogeny-ontogeny correlations is a hallmark of development in eukaryotes. This perspective opens a possibility to look for similar parallelisms in biological settings where developmental logic and multicellular complexity are more obscure. For instance, it has been increasingly recognized that multicellular behavior underlies biofilm formation in bacteria. However, it remains unclear whether bacterial biofilm growth shares some basic principles with development in complex eukaryotes. Here we show that the ontogeny of growing Bacillus subtilis biofilms recapitulates phylogeny at the expression level. Using time-resolved transcriptome and proteome profiles, we found that biofilm ontogeny correlates with the evolutionary measures, in a way that evolutionary younger and more diverged genes were increasingly expressed toward later timepoints of biofilm growth. Molecular and morphological signatures also revealed that biofilm growth is highly regulated and organized into discrete ontogenetic stages, analogous to those of eukaryotic embryos. Together, this suggests that biofilm formation in Bacillus is a bona fide developmental process comparable to organismal development in animals, plants, and fungi. Given that most cells on Earth reside in the form of biofilms and that biofilms represent the oldest known fossils, we anticipate that the widely adopted vision of the first life as a single-cell and free-living organism needs rethinking.

How can recapitulation be reconciled with modern concepts of evolution?

To understand Haeckel's idea of recapitulation with modern evolutionary biology, one has to realize how evolutionarily conserved embryonic stages appear sequentially in developmental processes as chains of causality. Whether the idea of evolution was accepted or not, Haeckel and von Baer commonly saw an importance of a particularly conserved mid-embryonic stage in biphasic development of metazoans, the phylotype, that defines an animal phylum as the developmental source of a basic body plan. In an evolutionary context, the phylotypic stage was once understood by Haeckel to reflect the common ancestor of animal phyla, which went through hypermorphosis independently into various phyla. Recent comprehensive molecular studies, however, accumulated data to refute this idea. The conserved embryonic pattern does not reflect an ancestral adult morphology but appears to have arisen primarily as an embodiment of developmental constraints established through evolutionary processes. How the developmental burden results in a nested series of constraints will solve the recapitulative tendency of developmental programs.

Network and Evolutionary Analysis of Human Epigenetic Regulators to Unravel Disease Associations

We carried out a system-level analysis of epigenetic regulators (ERs) and detailed the protein–protein interaction (PPI) network characteristics of disease-associated ERs. We found that most diseases associated with ERs can be clustered into two large groups, cancer diseases and developmental diseases. ER genes formed a highly interconnected PPI subnetwork, indicating a high tendency to interact and agglomerate with one another. We used the disease module detection (DIAMOnD) algorithm to expand the PPI subnetworks into a comprehensive cancer disease ER network (CDEN) and developmental disease ER network (DDEN). Using the transcriptome from early mouse developmental stages, we identified the gene co-expression modules significantly enriched for the CDEN and DDEN gene sets, which indicated the stage-dependent roles of ER-related disease genes during early embryonic development. The evolutionary rate and phylogenetic age distribution analysis indicated that the evolution of CDEN and DDEN genes was mostly constrained, and these genes exhibited older evolutionary age. Our analysis of human polymorphism data revealed that genes belonging to DDEN and Seed-DDEN were more likely to show signs of recent positive selection in human history. This finding suggests a potential association between positive selection of ERs and risk of developmental diseases through the mechanism of antagonistic pleiotropy.

Variations of telencephalic development that paved the way for neocortical evolution

Charles Darwin stated, "community in embryonic structure reveals community of descent". Thus, to understand how the neocortex emerged during mammalian evolution we need to understand the evolution of the development of the pallium, the source of the neocortex. In this article, we review the variations in the development of the pallium that enabled the production of the six-layered neocortex. We propose that an accumulation of subtle modifications from very early brain development accounted for the diversification of vertebrate pallia and the origin of the neocortex. Initially, faint differences of expression of secretable morphogens promote a wide variety in the proportions and organization of sectors of the early pallium in different vertebrates. It prompted different sectors to host varied progenitors and distinct germinative zones. These cells and germinative compartments generate diverse neuronal populations that migrate and mix with each other through radial and tangential migrations in a taxon-specific fashion. Together, these early variations had a profound influence on neurogenetic gradients, lamination, positioning, and connectivity. Gene expression, hodology, and physiological properties of pallial neurons are important features to suggest homologies, but the origin of cells and their developmental trajectory are fundamental to understand evolutionary changes. Our review compares the development of the homologous pallial sectors in sauropsids and mammals, with a particular focus on cell lineage, in search of the key changes that led to the appearance of the mammalian neocortex.

Toward understanding of evolutionary constraints: experimental and theoretical approaches

Although organisms have diversified remarkably through evolution, they do not exhibit unlimited variability. During evolution, the phenotypic changes do not occur at random; instead, they are directional and restricted by the constraints imposed on them. Despite the perceived importance of characterizing the unevenness of these changes, studies on evolutionary constraints have been primarily qualitative in nature. In this review, we focus on the recent studies of evolutionary constraints, which are based on the quantification of high-dimensional phenotypic and genotypic data. Furthermore, we present a theoretical analysis that enables us to predict evolutionary constraints on the basis of phenotypic fluctuation, modeled on the fluctuation-response relationship in statistical physics. The review lays emphasis on the tight interactions between experimental and theoretical analyses in evolutionary biology that will contribute to a better understanding of evolutionary constraints.

Patterns of Aging Biomarkers, Mortality, and Damaging Mutations Illuminate the Beginning of Aging and Causes of Early-Life Mortality

An increase in the probability of death has been a defining feature of aging, yet human perinatal mortality starts high and decreases with age. Previous evolutionary models suggested that organismal aging begins after the onset of reproduction. However, we find that mortality and incidence of diseases associated with aging follow a U-shaped curve with the minimum before puberty, whereas quantitative biomarkers of aging, including somatic mutations and DNA methylation, do not, revealing that aging starts early but is masked by early-life mortality. Moreover, our genetic analyses point to the contribution of damaging mutations to early mortality. We propose that mortality patterns are governed, in part, by negative selection against damaging mutations in early life, manifesting after the corresponding genes are first expressed. Deconvolution of mortality patterns suggests that deleterious changes rather than mortality are the defining characteristic of aging and that aging begins in very early life.

Inter-embryo gene expression variability recapitulates the hourglass pattern of evo-devo

BACKGROUND

The evolution of embryological development has long been characterized by deep conservation. In animal development, the phylotypic stage in mid-embryogenesis is more conserved than either early or late stages among species within the same phylum. Hypotheses to explain this hourglass pattern have focused on purifying the selection of gene regulation. Here, we propose an alternative-genes are regulated in different ways at different stages and have different intrinsic capacities to respond to perturbations on gene expression.

RESULTS

To eliminate the influence of natural selection, we quantified the expression variability of isogenetic single embryo transcriptomes throughout fly Drosophila melanogaster embryogenesis. We found that the expression variability is lower at the phylotypic stage, supporting that the underlying regulatory architecture in this stage is more robust to stochastic variation on gene expression. We present evidence that the phylotypic stage is also robust to genetic variations on gene expression. Moreover, chromatin regulation appears to play a key role in the variation and evolution of gene expression.

CONCLUSIONS

We suggest that a phylum-level pattern of embryonic conservation can be explained by the intrinsic difference of gene regulatory mechanisms in different stages.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Li et al. investigate the evolution and genetic basis of the adult pentameral body plan in echinoderms using genomic, transcriptomic, and proteomic data. They determine that the last common ancestor of echinoderms contained an organized Hox cluster and endoskeleton genes, and suggest that cooption of bilateral development genes was involved in evolution of the pentameric body plan.

Deeply conserved synteny resolves early events in vertebrate evolution

Although it is widely believed that early vertebrate evolution was shaped by ancient whole-genome duplications, the number, timing and mechanism of these events remain elusive. Here, we infer the history of vertebrates through genomic comparisons with a new chromosome-scale sequence of the invertebrate chordate amphioxus. We show how the karyotypes of amphioxus and diverse vertebrates are derived from 17 ancestral chordate linkage groups (and 19 ancestral bilaterian groups) by fusion, rearrangement and duplication. We resolve two distinct ancient duplications based on patterns of chromosomal conserved synteny. All extant vertebrates share the first duplication, which occurred in the mid/late Cambrian by autotetraploidization (that is, direct genome doubling). In contrast, the second duplication is found only in jawed vertebrates and occurred in the mid-late Ordovician by allotetraploidization (that is, genome duplication following interspecific hybridization) from two now-extinct progenitors. This complex genomic history parallels the diversification of vertebrate lineages in the fossil record.

Identification and characterization of male reproduction-related genes in pig (Sus scrofa) using transcriptome analysis

BACKGROUND

The systematic interrogation of reproduction-related genes was key to gain a comprehensive understanding of the molecular mechanisms underlying male reproductive traits in mammals. Here, based on the data collected from the NCBI SRA database, this study first revealed the genes involved in porcine male reproduction as well their uncharacterized transcriptional characteristics.

RESULTS

Results showed that the transcription of porcine genome was more widespread in testis than in other organs (the same for other mammals) and that testis had more tissue-specific genes (1210) than other organs. GO and GSEA analyses suggested that the identified test is-specific genes (TSGs) were associated with male reproduction. Subsequently, the transcriptional characteristics of porcine TSGs, which were conserved across different mammals, were uncovered. Data showed that 195 porcine TSGs shared similar expression patterns with other mammals (cattle, sheep, human and mouse), and had relatively higher transcription abundances and tissue specificity than low-conserved TSGs. Additionally, further analysis of the results suggested that alternative splicing, transcription factors binding, and the presence of other functionally similar genes were all involved in the regulation of porcine TSGs transcription.

CONCLUSIONS

Overall, this analysis revealed an extensive gene set involved in the regulation of porcine male reproduction and their dynamic transcription patterns. Data reported here provide valuable insights for a further improvement of the economic benefits of pigs as well as future treatments for male infertility.

pmar1/phb homeobox genes and the evolution of the double-negative gate for endomesoderm specification in echinoderms

In several model animals, the earliest phases of embryogenesis are regulated by lineage-specific genes, such as Drosophila bicoid Sea urchin (echinoid) embryogenesis is initiated by zygotic expression of pmar1, a paired-class homeobox gene that has been considered to be present only in the lineage of modern urchins (euechinoids). In euechinoids, Pmar1 promotes endomesoderm specification by repressing the hairy and enhancer of split C (hesC) gene. Here, we have identified the basal echinoid (cidaroid) pmar1 gene, which also promotes endomesoderm specification but not by repressing hesC A further search for related genes demonstrated that other echinoderms have pmar1-related genes named phb Functional analyses of starfish Phb proteins indicated that, similar to cidaroid Pmar1, they promote activation of endomesoderm regulatory gene orthologs via an unknown repressor that is not HesC. Based on these results, we propose that Pmar1 may have recapitulated the regulatory function of Phb during the early diversification of echinoids and that the additional repressor HesC was placed under the control of Pmar1 in the euechinoid lineage. This case provides an exceptional model for understanding how early developmental processes diverge.

Adaptation and Conservation throughout the Drosophila melanogaster Life-Cycle

Previous studies of the evolution of genes expressed at different life-cycle stages of Drosophila melanogaster have not been able to disentangle adaptive from nonadaptive substitutions when using nonsynonymous sites. Here, we overcome this limitation by combining whole-genome polymorphism data from D. melanogaster and divergence data between D. melanogaster and Drosophila yakuba. For the set of genes expressed at different life-cycle stages of D. melanogaster, as reported in modENCODE, we estimate the ratio of substitutions relative to polymorphism between nonsynonymous and synonymous sites (α) and then α is discomposed into the ratio of adaptive (ω_a) and nonadaptive (ω_na) substitutions to synonymous substitutions.

We find that the genes expressed in mid- and late-embryonic development are the most conserved, whereas those expressed in early development and postembryonic stages are the least conserved. Importantly, we found that low conservation in early development is due to high rates of nonadaptive substitutions (high ω_na), whereas in postembryonic stages it is due, instead, to high rates of adaptive substitutions (high ω_a).

By using estimates of different genomic features (codon bias, average intron length, exon number, recombination rate, among others), we also find that genes expressed in mid- and late-embryonic development show the most complex architecture: they are larger, have more exons, more transcripts, and longer introns. In addition, these genes are broadly expressed among all stages. We suggest that all these genomic features are related to the conservation of mid- and late-embryonic development. Globally, our study supports the hourglass pattern of conservation and adaptation over the life-cycle.

Adaptive Evolution of Animal Proteins over Development: Support for the Darwin Selection Opportunity Hypothesis of Evo-Devo

A driving hypothesis of evolutionary developmental biology is that animal morphological diversity is shaped both by adaptation and by developmental constraints. Here, we have tested Darwin’s “selection opportunity” hypothesis, according to which high evolutionary divergence in late development is due to strong positive selection. We contrasted it to a “developmental constraint” hypothesis, according to which late development is under relaxed negative selection. Indeed, the highest divergence between species, both at the morphological and molecular levels, is observed late in embryogenesis and postembryonically. To distinguish between adaptation and relaxation hypotheses, we investigated the evidence of positive selection on protein-coding genes in relation to their expression over development, in fly Drosophila melanogaster, zebrafish Danio rerio, and mouse Mus musculus. First, we found that genes specifically expressed in late development have stronger signals of positive selection. Second, over the full transcriptome, genes with evidence for positive selection trend to be expressed in late development. Finally, genes involved in pathways with cumulative evidence of positive selection have higher expression in late development. Overall, there is a consistent signal that positive selection mainly affects genes and pathways expressed in late embryonic development and in adult. Our results imply that the evolution of embryogenesis is mostly conservative, with most adaptive evolution affecting some stages of postembryonic gene expression, and thus postembryonic phenotypes. This is consistent with the diversity of environmental challenges to which juveniles and adults are exposed.

Gene expression across mammalian organ development

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, rhesus macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified correspondences of developmental stages across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

A hypothesis on paradoxical privileged portal vein metastasis of hepatocellular carcinoma. Can organ evolution shed light on patterns of human pathology, and vice versa?

Unlike other carcinomas, hepatocellular carcinoma (HCC) metastasizes to distant organs relatively rarely. In contrast, it routinely metastasizes to liver vasculature/liver, affecting portal veins 3-10 times more often than hepatic veins. This portal metastatic predominance is traditionally rationalized within the model of a reverse portal flow, due to accompanying liver cirrhosis. However, this intuitive model is not coherent with facts: 1) reverse portal flow occurs in fewer than 10% of cirrhotic patients, while portal metastasis occurs in 30-100% of HCC cases, and 2) portal vein prevalence of HCC metastasis is also characteristic of HCC in non-cirrhotic livers. Therefore, we must assume that the route for HCC metastatic dissemination is the same as for other carcinomas: systemic dissemination via the draining vessel, i.e., via the hepatic vein. In this light, portal prevalence versus hepatic vein of HCC metastasis appears as a puzzling pattern, particularly in cases when portal HCC metastases have appeared as the sole manifestation of HCC. Considering that other GI carcinomas (colorectal, pancreatic, gastric and small bowel) invariably disseminate via portal vein, but very rarely form portal metastasis, portal prevalence of HCC metastasis appears as a paradox. However, nature does not contradict itself; it is rather our wrong assumptions that create paradoxes. The 'portal paradox' becomes a logical event within the hypothesis that the formation of the unique portal venous system preceded the appearance of liver in evolution of chordates. The analysis suggests that the appearance of the portal venous system, supplying hormones and growth factors of pancreatic family, which includes insulin, glucagon, somatostatin, and pancreatic polypeptide (HGFPF) to midgut diverticulum in the early evolution of chordates (in an Amphioxus-like ancestral animal), promoted differentiation of enterocytes into hepatocytes and their further evolution to the liver of vertebrates. These promotional-dependent interactions are conserved in the vertebrate lineage. I hypothesize that selective homing and proliferation of malignant hepatocytes (i.e., HCC cells) in the portal vein environment are due to a uniquely high concentration of HGFPF in portal blood. HGFPF are also necessary for liver function and renewal and are significantly extracted by hepatocytes from passing blood, creating a concentration gradient of HGFPF between the portal blood and hepatic vein outflow, making post-liver vasculature and remote organs less favorable spaces for HCC growth. It also suggested that the portal vein environment (i.e., HGFPF) promotes the differentiation of more aggressive HCC clones from already-seeded portal metastases, explaining the worse outcome of HCC with the portal metastatic pattern. The analysis also offers new hypothesis on the phylogenetic origin of the hepatic diverticulum of cephalochordates, with certain implications for the modeling of the chordate phylogeny.

Gene Expression Does Not Support the Developmental Hourglass Model in Three Animals with Spiralian Development

It has been proposed that animals have a pattern of developmental evolution resembling an hourglass because the most conserved development stage—often called the phylotypic stage—is always in midembryonic development. Although the topic has been debated for decades, recent studies using molecular data such as RNA-seq gene expression data sets have largely supported the existence of periods of relative evolutionary conservation in middevelopment, consistent with the phylotypic stage and the hourglass concepts. However, so far this approach has only been applied to a limited number of taxa across the tree of life. Here, using established phylotranscriptomic approaches, we found a surprising reverse hourglass pattern in two molluscs and a polychaete annelid, representatives of the Spiralia, an understudied group that contains a large fraction of metazoan body plan diversity. These results suggest that spiralians have a divergent midembryonic stage, with more conserved early and late development, which is the inverse of the pattern seen in almost all other organisms where these phylotranscriptomic approaches have been reported. We discuss our findings in light of proposed reasons for the phylotypic stage and hourglass model in other systems.

The phylum Vertebrata: a case for zoological recognition

The group Vertebrata is currently placed as a subphylum in the phylum Chordata, together with two other subphyla, Cephalochordata (lancelets) and Urochordata (ascidians). The past three decades, have seen extraordinary advances in zoological taxonomy and the time is now ripe for reassessing whether the subphylum position is truly appropriate for vertebrates, particularly in light of recent advances in molecular phylogeny, comparative genomics, and evolutionary developmental biology. Four lines of current research are discussed here. First, molecular phylogeny has demonstrated that Deuterostomia comprises Ambulacraria (Echinodermata and Hemichordata) and Chordata (Cephalochordata, Urochordata, and Vertebrata), each clade being recognized as a mutually comparable phylum. Second, comparative genomic studies show that vertebrates alone have experienced two rounds of whole-genome duplication, which makes the composition of their gene family unique. Third, comparative gene-expression profiling of vertebrate embryos favors an hourglass pattern of development, the most conserved stage of which is recognized as a phylotypic period characterized by the establishment of a body plan definitively associated with a phylum. This mid-embryonic conservation is supported robustly in vertebrates, but only weakly in chordates. Fourth, certain complex patterns of body plan formation (especially of the head, pharynx, and somites) are recognized throughout the vertebrates, but not in any other animal groups. For these reasons, we suggest that it is more appropriate to recognize vertebrates as an independent phylum, not as a subphylum of the phylum Chordata.

Amphioxus functional genomics and the origins of vertebrate gene regulation

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

Genomic, epigenomic and transcriptomic data derived from the Mediterranean amphioxus (Branchiostoma lanceolatum) provide insights into the evolution of the genomic regulatory landscape of chordates.

Development and Evolutionary Constraints in Animals

We review the evolutionary importance of developmental mechanisms in constraining evolutionary changes in animals—in other words, developmental constraints. We focus on hard constraints that can act on macroevolutionary timescales. In particular, we discuss the causes and evolutionary consequences of the ancient metazoan constraint that differentiated cells cannot divide and constraints against changes of phylotypic stages in vertebrates and other higher taxa. We conclude that in all cases these constraints are caused by complex and highly controlled global interactivity of development, the disturbance of which has grave consequences. Mutations that affect such global interactivity almost unavoidably have many deleterious pleiotropic effects, which will be strongly selected against and will lead to long-term evolutionary stasis. The discussed developmental constraints have pervasive consequences for evolution and critically restrict regeneration capacity and body plan evolution.