Developmental milestones punctuate gene expression in the Caenorhabditis embryo

Widespread changes in gene expression accompany body size evolution in nematodes

Body size is a fundamental trait that drives multiple evolutionary and ecological patterns. Caenorhabditis inopinata is a fig-associated nematode that is exceptionally large relative to other members of the genus, including Caenorhabditis elegans. We previously showed that C. inopinata is large primarily due to postembryonic cell size expansion that occurs during the larval-to-adult transition. Here, we describe gene expression patterns in C. elegans and C. inopinata throughout this developmental period to understand the transcriptional basis of body size change. We performed RNA-seq in both species across the L3, L4, and adult stages. Most genes are differentially expressed across all developmental stages, consistent with C. inopinata's divergent ecology and morphology. We also used a model comparison approach to identify orthologues with divergent dynamics across this developmental period between the 2 species. This included genes connected to neurons, behavior, stress response, developmental timing, and small RNA/chromatin regulation. Multiple hypodermal collagens were also observed to harbor divergent developmental dynamics across this period, and genes important for molting and body morphology were also detected. Genes associated with transforming growth factor β signaling revealed idiosyncratic and unexpected transcriptional patterns given their role in body size regulation in C. elegans. This widespread transcriptional divergence between these species is unexpected and maybe a signature of the ecological and morphological divergence of C. inopinata. Alternatively, transcriptional turnover may be the rule in the Caenorhabditis genus, indicative of widespread developmental system drift among species. This work lays the foundation for future functional genetic studies interrogating the bases of body size evolution in this group.

High-resolution transcriptome datasets during embryogenesis of plant-parasitic nematodes

Understanding the transcriptional regulatory characteristics throughout the embryogenesis of plant-parasitic nematodes is crucial for elucidating their developmental processes' uniqueness. However, a challenge arises due to the lack of suitable technical methods for synchronizing the age of plant-parasitic nematodes embryo, it is difficult to collect detailed transcriptome data at each stage of embryonic development. Here, we recorded the 11 embryonic developmental time-points of endophytic nematode Meloidogyne incognita (isolated from Wuhan, China), Heterodera glycines (isolated from Wuhan, China), and Ditylenchus destructor (isolated from Jinan, China) species, and constructed transcriptome datasets of single embryos of these three species utilizing low-input smart-seq2 technology. The datasets encompassed 11 complete embryonic development stages, including Zygote, 2-cell, 4-cell, 8-cell, 24-44 cell, 64-78 cell, Comma, 1.5-fold, 2-fold, Moving, and L1, each stage generated four to five replicates, resulting in a total of 162 high-resolution transcriptome libraries. This high-resolution cross-species dataset serves as a crucial resource for comprehending the embryonic developmental properties of plant-parasitic nematodes and for identifying functional regulatory genes during embryogenesis.

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.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

An adapted MS2-MCP system to visualize endogenous cytoplasmic mRNA with live imaging in Caenorhabditis elegans

Live imaging of RNA molecules constitutes an invaluable means to track the dynamics of mRNAs, but live imaging in Caenorhabditis elegans has been difficult to achieve. Endogenous transcripts have been observed in nuclei, but endogenous mRNAs have not been detected in the cytoplasm, and functional mRNAs have not been generated. Here, we have adapted live imaging methods to visualize mRNA in embryonic cells. We have tagged endogenous transcripts with MS2 hairpins in the 3' untranslated region (UTR) and visualized them after adjusting MS2 Coat Protein (MCP) expression. A reduced number of these transcripts accumulates in the cytoplasm, leading to loss-of-function phenotypes. In addition, during epithelial morphogenesis, MS2-tagged mRNAs for dlg-1 fail to associate with the adherens junction, as observed for untagged, endogenous mRNAs. These defects are reversed by inactivating the nonsense-mediated decay pathway. RNA accumulates in the cytoplasm, mutant phenotypes are rescued, and dlg-1 RNA associates with the adherens junction. These data suggest that MS2 repeats can induce the degradation of endogenous RNAs and alter their cytoplasmic distribution. Although our focus is RNAs expressed in epithelial cells during morphogenesis, we find that this method can be applied to other cell types and stages.

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.

Integrating Complex Life Cycles in Comparative Developmental Biology

The goal of comparative developmental biology is identifying mechanistic differences in embryonic development between different taxa and how these evolutionary changes have led to morphological and organizational differences in adult body plans. Much of this work has focused on direct-developing species in which the adult forms straight from the embryo and embryonic modifications have direct effects on the adult. However, most animal lineages are defined by indirect development, in which the embryo gives rise to a larval body plan and the adult forms by transformation of the larva. Historically, much of our understanding of complex life cycles is viewed through the lenses of ecology and zoology. In this review, we discuss the importance of establishing developmental rather than morphological or ecological criteria for defining developmental mode and explicitly considering the evolutionary implications of incorporating complex life cycles into broad developmental comparisons of embryos across metazoans.

Widespread changes in gene expression accompany body size evolution in nematodes

Body size is a fundamental trait that drives multiple evolutionary and ecological patterns. Caenorhabditis inopinata is a fig-associated nematode that is exceptionally large relative to other members of the genus, including C. elegans. We previously showed that C. inopinata is large primarily due to postembryonic cell size expansion that occurs during the larval-to-adult transition. Here, we describe gene expression patterns in C. elegans and C. inopinata throughout this developmental period to understand the transcriptional basis of body size change. We performed RNA-seq in both species across the L3, L4, and adult stages. Most genes are differentially expressed across all developmental stages, consistent with C. inopinata's divergent ecology and morphology. We also used a model comparison approach to identify orthologs with divergent dynamics across this developmental period between the two species. This included genes connected to neurons, behavior, stress response, developmental timing, and small RNA/chromatin regulation. Multiple hypodermal collagens were also observed to harbor divergent developmental dynamics across this period, and genes important for molting and body morphology were also detected. Genes associated with TGF-β signaling revealed idiosyncratic and unexpected transcriptional patterns given their role in body size regulation in C. elegans. Widespread transcriptional divergence between these species is unexpected and may be a signature of the ecological and morphological divergence of C. inopinata. Alternatively, transcriptional turnover may be the rule in the Caenorhabditis genus, indicative of widespread developmental system drift among species. This work lays the foundation for future functional genetic studies interrogating the bases of body size evolution in this group.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

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.

A novel time-lapse imaging method for studying developing bacterial biofilms

In nature, bacteria prevailingly reside in the form of biofilms. These elaborately organized surface-bound assemblages of bacterial cells show numerous features of multicellular organization. We recently showed that biofilm growth is a true developmental process, which resembles developmental processes in multicellular eukaryotes. To study the biofilm growth in a fashion of eukaryotic ontogeny, it is essential to define dynamics and critical transitional phases of this process. The first step in this endeavor is to record the gross morphological changes of biofilm ontogeny under standardized conditions. This visual information is instrumental in guiding the sampling strategy for the later omics analyses of biofilm ontogeny. However, none of the currently available visualizations methods is specifically tailored for recording gross morphology across the whole biofilm development. To address this void, here we present an affordable Arduino-based approach for time-lapse visualization of complete biofilm ontogeny using bright field stereomicroscopy with episcopic illumination. The major challenge in recording biofilm development on the air-solid interphase is water condensation, which compromises filming directly through the lid of a Petri dish. To overcome these trade-offs, we developed an Arduino microcontroller setup which synchronizes a robotic arm, responsible for opening and closing the Petri dish lid, with the activity of a stereomicroscope-mounted camera and lighting conditions. We placed this setup into a microbiological incubator that maintains temperature and humidity during the biofilm growth. As a proof-of-principle, we recorded biofilm development of five Bacillus subtilis strains that show different morphological and developmental dynamics.

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.

Phylotranscriptomic patterns of network stochasticity and pathway dynamics during embryogenesis

MOTIVATION

The hourglass model is a popular evo-devo model depicting that the developmental constraints in the middle of a developmental process are higher, and hence the phenotypes are evolutionarily more conserved, than those that occur in early and late ontogeny stages. Although this model has been supported by studies analyzing developmental gene expression data, the evolutionary explanation and molecular mechanism behind this phenomenon are not fully understood yet. To approach this problem, Raff proposed a hypothesis and claimed that higher interconnectivity among elements in an organism during organogenesis resulted in the larger constraints at the mid-developmental stage. By employing stochastic network analysis and gene-set pathway analysis, we aim to demonstrate such changes of interconnectivity claimed in Raff's hypothesis.

RESULTS

We first compared the changes of network randomness among developmental processes in different species by measuring the stochasticity within the biological network in each developmental stage. By tracking the network entropy along each developmental process, we found that the network stochasticity follows an anti-hourglass trajectory, and such a pattern supports Raff's hypothesis in dynamic changes of interconnections among biological modules during development. To understand which biological functions change during the transition of network stochasticity, we sketched out the pathway dynamics along the developmental stages and found that species may activate similar groups of biological processes across different stages. Moreover, higher interspecies correlations are found at the mid-developmental stages.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Neuronal control of maternal provisioning in response to social cues

Mothers contribute cytoplasmic components to their progeny in a process called maternal provisioning. Provisioning is influenced by the parental environment, but the molecular pathways that transmit environmental cues between generations are not well understood. Here, we show that, in Caenorhabditis elegans, social cues modulate maternal provisioning to regulate gene silencing in offspring. Intergenerational signal transmission depends on a pheromone-sensing neuron and neuronal FMRFamide (Phe-Met-Arg-Phe)-like peptides. Parental FMRFamide-like peptide signaling dampens oxidative stress resistance and promotes the deposition of mRNAs for translational components in progeny, which, in turn, reduces gene silencing. This study identifies a previously unknown pathway for intergenerational communication that links neuronal responses to maternal provisioning. We suggest that loss of social cues in the parental environment represents an adverse environment that stimulates stress responses across generations.

Single worm transcriptomics identifies a developmental core network of oscillating genes with deep conservation across nematodes

High-resolution spatial and temporal maps of gene expression have facilitated a comprehensive understanding of animal development and evolution. In nematodes, the small body size represented a major challenge for such studies, but recent advancements have helped overcome this limitation. Here, we have implemented single worm transcriptomics (SWT) in the nematode model organism Pristionchus pacificus to provide a high-resolution map of the developmental transcriptome. We selected 38 time points from hatching of the J2 larvae to young adults to perform transcriptome analysis over 60 h of postembryonic development. A mean sequencing depth of 4.5 million read pairs allowed the detection of more than 23,135 (80%) of all genes. Nearly 3000 (10%) genes showed oscillatory expression with discrete expression levels, phases, and amplitudes. Gene age analysis revealed an overrepresentation of ancient gene classes among oscillating genes, and around one-third of them have 1:1 orthologs in C. elegans. One important gene family overrepresented among oscillating genes is collagens. Several of these collagen genes are regulated by the developmental switch gene eud-1, indicating a potential function in the regulation of mouth-form plasticity, a key developmental process in this facultative predatory nematode. Together, our analysis provides (1) an updated protocol for SWT in nematodes that is applicable to many microscopic species, (2) a 1- to 2-h high-resolution catalog of P. pacificus gene expression throughout postembryonic development, and (3) a comparative analysis of oscillatory gene expression between the two model organisms P. pacificus and C. elegans and associated evolutionary dynamics.

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.

Comparative Transcriptomics across Nematode Life Cycles Reveal Gene Expression Conservation and Correlated Evolution in Adjacent Developmental Stages

Nematodes are highly abundant animals with diverse habitats and lifestyles. Some are free living whereas others parasitize animals or plants, and among the latter, infection abilities change across developmental stages to infect hosts and complete life cycles. To determine the relationship between transcriptome evolution and morphological divergences among nematodes, we compared 48 transcriptomes of different developmental stages across eight nematode species. The transcriptomes were clustered broadly into embryo, larva, and adult stages, with the developmental plastic stages were separated from common larval stages within the larval branch. This suggests that development was the major determining factor after lifestyle changes, such as parasitism, during transcriptome evolution. Such patterns were partly accounted for by tissue-specific genes—such as those in oocytes and the hypodermis—being expressed at different proportions. Although nematodes typically have 3–5 larval stages, the transcriptomes for these stages were found to be highly correlated within each species, suggesting high similarity among larval stages across species. For the Caenorhabditis elegans–Caenorhabditis briggsae and Strongyloides stercoralis–Strongyloides venezuelensis comparisons, we found that ∼50% of genes were expressed at multiple stages, whereas half of their orthologs were also expressed in multiple but different stages. Such frequent changes in expression have resulted in concerted transcriptome evolution across adjacent stages, thus generating species-specific transcriptomes over the course of nematode evolution. Our study provides a first insight into the evolution of nematode transcriptomes beyond embryonic development.

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.

The tweety Gene Family: From Embryo to Disease

The tweety genes encode gated chloride channels that are found in animals, plants, and even simple eukaryotes, signifying their deep evolutionary origin. In vertebrates, the tweety gene family is highly conserved and consists of three members—ttyh1, ttyh2, and ttyh3—that are important for the regulation of cell volume. While research has elucidated potential physiological functions of ttyh1 in neural stem cell maintenance, proliferation, and filopodia formation during neural development, the roles of ttyh2 and ttyh3 are less characterized, though their expression patterns during embryonic and fetal development suggest potential roles in the development of a wide range of tissues including a role in the immune system in response to pathogen-associated molecules. Additionally, members of the tweety gene family have been implicated in various pathologies including cancers, particularly pediatric brain tumors, and neurodegenerative diseases such as Alzheimer’s and Parkinson’s disease. Here, we review the current state of research using information from published articles and open-source databases on the tweety gene family with regard to its structure, evolution, expression during development and adulthood, biochemical and cellular functions, and role in human disease. We also identify promising areas for further research to advance our understanding of this important, yet still understudied, family of genes.

Systematic screening of long intergenic noncoding RNAs expressed during chicken embryogenesis

Long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes, including embryogenesis and development. To provide a systematic analysis of lncRNAs expressed during chicken embryogenesis, we used Iso-Seq and RNA-Seq to identify potential lncRNAs at embryonic stages from d 1 to d 8 of incubation: sequential stages covering gastrulation, somitogenesis, and organogenesis. The data characterized an expanded landscape of lncRNAs, yielding 45,410 distinct lncRNAs (31,282 genes). Amongst these, a set of 13,141 filtered intergenic lncRNAs (lincRNAs) transcribed from 9803 lincRNA gene loci, of which, 66.5% were novel, were further analyzed. These lincRNAs were found to share many characteristics with mammalian lincRNAs, including relatively short lengths, fewer exons, lower expression levels, and stage-specific expression patterns. Functional studies motivated by "guilt-by-association" associated individual lincRNAs with specific GO functions, providing an important resource for future studies of lincRNA function. Most importantly, a weighted gene co-expression network analysis suggested that genes of the brown module were specifically associated with the day 2 stage. LincRNAs within this module were co-expressed with proteins involved in hematopoiesis and lipid metabolism. This study presents the systematic identification of lincRNAs in developing chicken embryos and will serve as a powerful resource for the study of lincRNA functions.

The biogenetic law and the Gastraea theory: From Ernst Haeckel's discoveries to contemporary views

More than 150 years ago, in 1866, Ernst Haeckel published a book in two volumes called Generelle Morphologie der Organismen (General Morphology of Organisms) in the first volume of which he formulated his biogenetic law, famously stating that ontogeny recapitulates phylogeny. Here, we describe Haeckel's original idea as first formulated in the Generelle Morphologie der Organismen and later further developed in other publications until the present situation in which molecular data are used to test the "hourglass model," which can be seen as a modern version of the biogenetic law. We also tell the story about his discovery, while traveling in Norway, of an unknown organism, Magosphaera planula, that was important in that it helped to precipitate his ideas into what was to become the Gastraea theory. We also follow further development and reformulations of the Gastraea theory by other scientists, notably the Russian school. Elias Metchnikoff developed the Phagocytella hypothesis for the origin of metazoans based on studies of a colonial flagellate. Alexey Zakhvatin focused on deducing the ancestral life cycle and the cell types of the last common ancestor of all metazoans, and Kirill V. Mikhailov recently pursued this line of research further.

Dynamical systems approach to evolution-development congruence: Revisiting Haeckel's recapitulation theory

It is acknowledged that embryonic development has a tendency to proceed from common toward specific. Ernst Haeckel raised the question of why that tendency prevailed through evolution, and the question remains unsolved. Here, we revisit Haeckel's recapitulation theory, that is, the parallelism between evolution and development through numerical evolution and dynamical systems theory. By using intracellular gene expression dynamics with cell‐to‐cell interaction over spatially aligned cells to represent the developmental process, gene regulation networks (GRN) that govern these dynamics evolve under the selection pressure to achieve a prescribed spatial gene expression pattern. For most numerical evolutionary experiments, the evolutionary pattern changes over generations, as well as the developmental pattern changes governed by the evolved GRN exhibit remarkable similarity. Changes in both patterns consisted of several epochs where stripes are formed in a short time, whereas for other temporal regimes, the pattern hardly changes. In evolution, these quasi‐stationary generations are needed to achieve relevant mutations, whereas, in development, they are due to some gene expressions that vary slowly and control the pattern change. These successive epochal changes in development and evolution are represented as common bifurcations in dynamical systems theory, regulating working network structure from feedforward subnetwork to those containing feedback loops. The congruence is the correspondence between successive acquisitions of subnetworks through evolution and changes in working subnetworks in development. Consistency of the theory with the segmentation gene‐expression dynamics is discussed. Novel outlook on recapitulation and heterochrony are provided, testable experimentally by the transcriptome and network analysis.

We confirmed recapitulation, parallelism between evolution and development, by evolutionary simulation of pattern formation. Its origin is analyzed in terms of bifurcation in dynamical systems theory, slow‐mode control, and network structures.

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.

Syndecan receptors: pericellular regulators in development and inflammatory disease

The syndecans are the major family of transmembrane proteoglycans, usually bearing multiple heparan sulfate chains. They are present on virtually all nucleated cells of vertebrates and are also present in invertebrates, indicative of a long evolutionary history. Genetic models in both vertebrates and invertebrates have shown that syndecans link to the actin cytoskeleton and can fine-tune cell adhesion, migration, junction formation, polarity and differentiation. Although often associated as co-receptors with other classes of receptors (e.g. integrins, growth factor and morphogen receptors), syndecans can nonetheless signal to the cytoplasm in discrete ways. Syndecan expression levels are upregulated in development, tissue repair and an array of human diseases, which has led to the increased appreciation that they may be important in pathogenesis not only as diagnostic or prognostic agents, but also as potential targets. Here, their functions in development and inflammatory diseases are summarized, including their potential roles as conduits for viral pathogen entry into cells.

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.

High-Temporal-Resolution smFISH Method for Gene Expression Studies in Caenorhabditis elegans Embryos

Recent development in fluorescence-based molecular tools has contributed significantly to developmental studies, including embryogenesis. Many of these tools rely on multiple steps of sample manipulation, so obtaining large sample sizes presents a major challenge as it can be labor-intensive and time-consuming. However, large sample sizes are required to uncover critical aspects of embryogenesis, for example, subtle phenotypic differences or gene expression dynamics. This problem is particularly relevant for single-molecule fluorescence in situ hybridization (smFISH) studies in Caenorhabditis elegans embryogenesis. Microfluidics can help address this issue by allowing a large number of samples and parallelization of experiments. However, performing efficient reagent exchange on chip for large numbers of embryos remains a bottleneck. Here, we present a microfluidic pipeline for large-scale smFISH imaging of C. elegans embryos with minimized labor. We designed embryo traps and engineered a protocol allowing for efficient chemical exchange for hundreds of C. elegans embryos simultaneously. Furthermore, the device design and small footprint optimize imaging throughput by facilitating spatial registration and enabling minimal user input. We conducted the smFISH protocol on chip and demonstrated that image quality is preserved. With one device replacing the equivalent of 10 glass slides of embryos mounted manually, our microfluidic approach greatly increases throughput. Finally, to highlight the capability of our platform to perform longitudinal studies with high temporal resolution, we conducted a temporal analysis of par-1 gene expression in early C. elegans embryos. The method demonstrated here paves the way for systematic high-temporal-resolution studies that will benefit large-scale RNAi and drug screens and in systems beyond C. elegans embryos.

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.

In Vivo Simultaneous Analysis of Gene Expression by Dual-Color Luciferases in Caenorhabditis elegans

Both fluorescent and luminescent observation are widely used to examine real-time gene expression patterns in living organisms. Several fluuorescent and luminescent proteins with specific optical properties have been developed and applied for simultaneous, multi-color observation of more than two gene expression profiles. Compared to fluorescent proteins, however, the application of multi-color luminescent imaging in living organisms is still limited. In this study, we introduced two-color luciferases into the soil nematode C. elegans and performed simultaneous analysis of two gene expression profiles. Using a green-emitting luciferase Eluc (emerald luciferase) and red-emitting luciferase SLR (stable luciferase red), the expression patterns of two genes were simultaneously observed in single animals from embryonic to adult stages over its whole life span. In addition, dual gene activities were observed at the single embryo level, with the simultaneous observation of morphological changes. These are the first application of a two-color luciferase system into a whole animal and suggest that precise relationship of expression patterns of multiple genes of interest can be analyzed over the whole life of the animal, dependent on the changes in genetic and/or environmental conditions.

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.

Speciation and the developmental alarm clock

New species arise as the genomes of populations diverge. The developmental 'alarm clock' of speciation sounds off when sufficient divergence in genetic control of development leads hybrid individuals to infertility or inviability, the world awoken to the dawn of new species with intrinsic post-zygotic reproductive isolation. Some developmental stages will be more prone to hybrid dysfunction due to how molecular evolution interacts with the ontogenetic timing of gene expression. Considering the ontogeny of hybrid incompatibilities provides a profitable connection between 'evo-devo' and speciation genetics to better link macroevolutionary pattern, microevolutionary process, and molecular mechanisms. Here, we explore speciation alongside development, emphasizing their mutual dependence on genetic network features, fitness landscapes, and developmental system drift. We assess models for how ontogenetic timing of reproductive isolation can be predictable. Experiments and theory within this synthetic perspective can help identify new rules of speciation as well as rules in the molecular evolution of development.

Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms

Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.

Conserved Patterns in Developmental Processes and Phases, Rather than Genes, Unite the Highly Divergent Bilateria

Bilateria are the predominant clade of animals on Earth. Despite having evolved a wide variety of body plans and developmental modes, they are characterized by common morphological traits. By default, researchers have tried to link clade-specific genes to these traits, thus distinguishing bilaterians from non-bilaterians, by their gene content. Here we argue that it is rather biological processes that unite Bilateria and set them apart from their non-bilaterian sisters, with a less complex body morphology. To test this hypothesis, we compared proteomes of bilaterian and non-bilaterian species in an elaborate computational pipeline, aiming to search for a set of bilaterian-specific genes. Despite the limited confidence in their bilaterian specificity, we nevertheless detected Bilateria-specific functional and developmental patterns in the sub-set of genes conserved in distantly related Bilateria. Using a novel multi-species GO-enrichment method, we determined the functional repertoire of genes that are widely conserved among Bilateria. Analyzing expression profiles in three very distantly related model species-D. melanogaster, D. rerio and C. elegans-we find characteristic peaks at comparable stages of development and a delayed onset of expression in embryos. In particular, the expression of the conserved genes appears to peak at the phylotypic stage of different bilaterian phyla. In summary, our study illustrate how development connects distantly related Bilateria after millions of years of divergence, pointing to processes potentially separating them from non-bilaterians. We argue that evolutionary biologists should return from a purely gene-centric view of evolution and place more focus on analyzing and defining conserved developmental processes and periods.

Evolution and Developmental System Drift in the Endoderm Gene Regulatory Network of Caenorhabditis and Other Nematodes

Developmental gene regulatory networks (GRNs) underpin metazoan embryogenesis and have undergone substantial modification to generate the tremendous variety of animal forms present on Earth today. The nematode Caenorhabditis elegans has been a central model for advancing many important discoveries in fundamental mechanistic biology and, more recently, has provided a strong base from which to explore the evolutionary diversification of GRN architecture and developmental processes in other species. In this short review, we will focus on evolutionary diversification of the GRN for the most ancient of the embryonic germ layers, the endoderm. Early embryogenesis diverges considerably across the phylum Nematoda. Notably, while some species deploy regulative development, more derived species, such as C. elegans, exhibit largely mosaic modes of embryogenesis. Despite the relatively similar morphology of the nematode gut across species, widespread variation has been observed in the signaling inputs that initiate the endoderm GRN, an exemplar of developmental system drift (DSD). We will explore how genetic variation in the endoderm GRN helps to drive DSD at both inter- and intraspecies levels, thereby resulting in a robust developmental system. Comparative studies using divergent nematodes promise to unveil the genetic mechanisms controlling developmental plasticity and provide a paradigm for the principles governing evolutionary modification of an embryonic GRN.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

Distinct functions and temporal regulation of methylated histone H3 during early embryogenesis

During the first hours of embryogenesis, formation of higher-order heterochromatin coincides with the loss of developmental potential. Here, we examine the relationship between these two events, and we probe the processes that contribute to the timing of their onset. Mutations that disrupt histone H3 lysine 9 (H3K9) methyltransferases reveal that the methyltransferase MET-2 helps terminate developmental plasticity, through mono- and di-methylation of H3K9 (me1/me2), and promotes heterochromatin formation, through H3K9me3. Although loss of H3K9me3 perturbs formation of higher-order heterochromatin, embryos are still able to terminate plasticity, indicating that the two processes can be uncoupled. Methylated H3K9 appears gradually in developing C. elegans embryos and depends on nuclear localization of MET-2. We find that the timing of H3K9me2 and nuclear MET-2 is sensitive to rapid cell cycles, but not to zygotic genome activation or cell counting. These data reveal distinct roles for different H3K9 methylation states in the generation of heterochromatin and loss of developmental plasticity by MET-2, and identify the cell cycle as a crucial parameter of MET-2 regulation.

C. elegans synMuv B proteins regulate spatial and temporal chromatin compaction during development

Tissue-specific establishment of repressive chromatin through creation of compact chromatin domains during development is necessary to ensure proper gene expression and cell fate. Caenorhabditis elegans synMuv B proteins are important for the soma/germline fate decision and mutants demonstrate ectopic germline gene expression in somatic tissue, especially at high temperature. We show that C. elegans synMuv B proteins regulate developmental chromatin compaction and that the timing of chromatin compaction is temperature sensitive in both wild type and synMuv B mutants. Chromatin compaction in mutants is delayed into developmental time periods when zygotic gene expression is upregulated and demonstrates an anterior-to-posterior pattern. Loss of this patterned compaction coincides with the developmental time period of ectopic germline gene expression, which leads to a developmental arrest in synMuv B mutants. Finally, accelerated cell division rates at elevated temperature may contribute to a lack of coordination between expression of tissue specific transcription programs and chromatin compaction at high temperature. Thus, chromatin organization during development is regulated both spatially and temporally by synMuv B proteins to establish repressive chromatin in a tissue-specific manner to ensure proper gene expression.

Reconstructing the Transcriptional Ontogeny of Maize and Sorghum Supports an Inverse Hourglass Model of Inflorescence Development

Assembling meaningful comparisons between species is a major limitation in studying the evolution of organismal form. To understand development in maize and sorghum, closely related species with architecturally distinct inflorescences, we collected RNA-seq profiles encompassing inflorescence body-plan specification in both species. We reconstructed molecular ontogenies from 40 B73 maize tassels and 47 BTx623 sorghum panicles and separated them into transcriptional stages. To discover new markers of inflorescence development, we used random forest machine learning to determine stage by RNA-seq. We used two descriptions of transcriptional conservation to identify hourglass-like stages during inflorescence development. Despite a relatively short 12 million years since their last common ancestor, we found maize and sorghum inflorescences are most different during their hourglass-like stages of development, following an inverse-hourglass model of development. We discuss whether agricultural selection may account for the rapid divergence signatures in these species and the observed separation of evolutionary pressure and developmental reprogramming.

How Weird is The Worm? Evolution of the Developmental Gene Toolkit in Caenorhabditis elegans

Comparative developmental biology and comparative genomics are the cornerstones of evolutionary developmental biology. Decades of fruitful research using nematodes have produced detailed accounts of the developmental and genomic variation in the nematode phylum. Evolutionary developmental biologists are now utilising these data as a tool with which to interrogate the evolutionary basis for the similarities and differences observed in Nematoda. Nematodes have often seemed atypical compared to the rest of the animal kingdom-from their totally lineage-dependent mode of embryogenesis to their abandonment of key toolkit genes usually deployed by bilaterians for proper development-worms are notorious rule breakers of the bilaterian handbook. However, exploring the nature of these deviations is providing answers to some of the biggest questions about the evolution of animal development. For example, why is the evolvability of each embryonic stage not the same? Why can evolution sometimes tolerate the loss of genes involved in key developmental events? Lastly, why does natural selection act to radically diverge toolkit genes in number and sequence in certain taxa? In answering these questions, insight is not only being provided about the evolution of nematodes, but of all metazoans.

Ancestral roles of the Fam20C family of secreted protein kinases revealed in C. elegans

Gerson-Gurwitz et al. show that FAMK-1, the sole Fam20-like secreted protein kinase in C. elegans, contributes to embryonic development and fertility, acts in the late secretory pathway, phosphorylates lectins, and is important in embryonic and tissue contexts where cells are subjected to mechanical strain.

Fam20C is a secreted protein kinase mutated in Raine syndrome, a human skeletal disorder. In vertebrates, bone and enamel proteins are major Fam20C substrates. However, Fam20 kinases are conserved in invertebrates lacking bone and enamel, suggesting other ancestral functions. We show that FAMK-1, the Caenorhabditis elegans Fam20C orthologue, contributes to fertility, embryogenesis, and development. These functions are not fulfilled when FAMK-1 is retained in the early secretory pathway. During embryogenesis, FAMK-1 maintains intercellular partitions and prevents multinucleation; notably, temperature elevation or lowering cortical stiffness reduces requirement for FAMK-1 in these contexts. FAMK-1 is expressed in multiple adult tissues that undergo repeated mechanical strain, and selective expression in the spermatheca restores fertility. Informatic, biochemical, and functional analysis implicate lectins as FAMK-1 substrates. These findings suggest that FAMK-1 phosphorylation of substrates, including lectins, in the late secretory pathway is important in embryonic and tissue contexts where cells are subjected to mechanical strain.

Dynamic Transcriptional Landscape of the Early Chick Embryo

Defining the dynamic transcriptome of the early embryo at high resolution would assist greatly in understanding vertebrate development. Here, we describe the dynamic transcription landscape of early chick embryo development using advanced single-molecule long-read isoform sequencing (Iso-Seq) and RNA-Seq technology. Our transcriptomic profiling reflected the time course of chicken embryonic development from day 1 to day 8 of incubation, a period encompassing gastrulation, somitogenesis, and organogenesis. This analysis identified transcriptional isoforms, alternative splicing (AS) events, fusion transcripts, alternative polyadenylation (APA) sites, and novel genes. Our results showed that intron retention (IR) represented the most abundant AS type and displayed distinct features and dynamic modulation during development. Moreover, we constructed a high-resolution expression profile across embryonic development. Our combined expression dataset correlates distinct gene clusters with specific morphological changes, and provides the first framework for the molecular basis of early chicken embryogenesis. Analysis of gene expression in the developing chicken embryo highlighted the dynamic nature and complexity of the chicken transcriptome and demonstrated that dramatically increased IR events are associated with distinct gene sets.

Extensive intraspecies cryptic variation in an ancient embryonic gene regulatory network

Innovations in metazoan development arise from evolutionary modification of gene regulatory networks (GRNs). We report widespread cryptic variation in the requirement for two key regulatory inputs, SKN-1/Nrf2 and MOM-2/Wnt, into the C. elegans endoderm GRN. While some natural isolates show a nearly absolute requirement for these two regulators, in others, most embryos differentiate endoderm in their absence. GWAS and analysis of recombinant inbred lines reveal multiple genetic regions underlying this broad phenotypic variation. We observe a reciprocal trend, in which genomic variants, or knockdown of endoderm regulatory genes, that result in a high SKN-1 requirement often show low MOM-2/Wnt requirement and vice-versa, suggesting that cryptic variation in the endoderm GRN may be tuned by opposing requirements for these two key regulatory inputs. These findings reveal that while the downstream components in the endoderm GRN are common across metazoan phylogeny, initiating regulatory inputs are remarkably plastic even within a single species.

Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant

Gene regulatory networks (GRNs) consist of different molecular interactions that closely work together to establish proper gene expression in time and space. Especially in higher eukaryotes, many questions remain on how these interactions collectively coordinate gene regulation. We study high quality GRNs consisting of undirected protein–protein, genetic and homologous interactions, and directed protein–DNA, regulatory and miRNA–mRNA interactions in the worm Caenorhabditis elegans and the plant Arabidopsis thaliana. Our data-integration framework integrates interactions in composite network motifs, clusters these in biologically relevant, higher-order topological network motif modules, overlays these with gene expression profiles and discovers novel connections between modules and regulators. Similar modules exist in the integrated GRNs of worm and plant. We show how experimental or computational methodologies underlying a certain data type impact network topology. Through phylogenetic decomposition, we found that proteins of worm and plant tend to functionally interact with proteins of a similar age, while at the regulatory level TFs favor same age, but also older target genes. Despite some influence of the duplication mode difference, we also observe at the motif and module level for both species a preference for age homogeneity for undirected and age heterogeneity for directed interactions. This leads to a model where novel genes are added together to the GRNs in a specific biological functional context, regulated by one or more TFs that also target older genes in the GRNs. Overall, we detected topological, functional and evolutionary properties of GRNs that are potentially universal in all species.

Molecular evolution across developmental time reveals rapid divergence in early embryogenesis

Ontogenetic development hinges on the changes in gene expression in time and space within an organism, suggesting that the demands of ontogenetic growth can impose or reveal predictable pattern in the molecular evolution of genes expressed dynamically across development. Here, we characterize coexpression modules of the Caenorhabditis elegans transcriptome, using a time series of 30 points from early embryo to adult. By capturing the functional form of expression profiles with quantitative metrics, we find fastest evolution in the distinctive set of genes with transcript abundance that declines through development from a peak in young embryos. These genes are highly enriched for oogenic function and transient early zygotic expression, are nonrandomly distributed in the genome, and correspond to a life stage especially prone to inviability in interspecies hybrids. These observations conflict with the "early conservation model" for the evolution of development, although expression-weighted sequence divergence analysis provides some support for the "hourglass model." Genes in coexpression modules that peak toward adulthood also evolve fast, being hyper-enriched for roles in spermatogenesis, implicating a history of sexual selection and relaxation of selection on sperm as key factors driving rapid change to ontogenetically distinguishable coexpression modules of genes. We propose that these predictable trends of molecular evolution for dynamically expressed genes across ontogeny predispose particular life stages, early embryogenesis in particular, to hybrid dysfunction in the speciation process.

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.

Repression of Germline Genes in Caenorhabditis elegans Somatic Tissues by H3K9 Dimethylation of Their Promoters

Repression of germline-promoting genes in somatic cells is critical for somatic development and function. To study how germline genes are repressed in somatic tissues, we analyzed key histone modifications in three Caenorhabditis elegans synMuv B mutants, lin-15B, lin-35, and lin-37-all of which display ectopic expression of germline genes in the soma. LIN-35 and LIN-37 are members of the conserved DREAM complex. LIN-15B has been proposed to work with the DREAM complex but has not been shown biochemically to be a member of the complex. We found that, in wild-type worms, synMuv B target genes and germline genes are enriched for the repressive histone modification dimethylation of histone H3 on lysine 9 (H3K9me2) at their promoters. Genes with H3K9me2 promoter localization are evenly distributed across the autosomes, not biased toward autosomal arms, as are the broad H3K9me2 domains. Both synMuv B targets and germline genes display a dramatic reduction of H3K9me2 promoter localization in lin-15B mutants, but much weaker reduction in lin-35 and lin-37 mutants. This difference between lin-15B and DREAM complex mutants likely represents a difference in molecular function for these synMuv B proteins. In support of the pivotal role of H3K9me2 in regulation of germline genes by LIN-15B, global loss of H3K9me2 but not H3K9me3 results in phenotypes similar to synMuv B mutants, high-temperature larval arrest, and ectopic expression of germline genes in the soma. We propose that LIN-15B-driven enrichment of H3K9me2 at promoters of germline genes contributes to repression of those genes in somatic tissues.