Comparative analysis of embryonic cell lineage between Caenorhabditis briggsae and Caenorhabditis elegans

Divergence in neuronal signaling pathways despite conserved neuronal identity among Caenorhabditis species

One avenue to better understand brain evolution is to map molecular patterns of evolutionary changes in neuronal cell types across entire nervous systems of distantly related species. Generating whole-animal single-cell transcriptomes of three nematode species from the Caenorhabditis genus, we observed a remarkable stability of neuronal-cell-type identities over more than 45 million years of evolution. Conserved patterns of combinatorial expression of homeodomain transcription factors are among the best classifiers of homologous neuron classes. Unexpectedly, we discover an extensive divergence in neuronal signaling pathways. Although identities of neurotransmitter-producing neurons (glutamate, acetylcholine, γ-aminobutyric acid [GABA], and several monoamines) remain stable, expression of ionotropic and metabotropic receptors for all these neurotransmitter systems shows substantial divergence, resulting in more than half of all neuron classes changing their capacity to be receptive to specific neurotransmitters. Neuropeptidergic signaling is also remarkably divergent, both at the level of neuropeptide expression and receptor expression, yet the overall dense network topology of the wireless neuropeptidergic connectome remains stable. Novel neuronal signaling pathways are suggested by our discovery of small secreted proteins that show no obvious hallmarks of conventional neuropeptides but show similar patterns of highly neuron-type-specific and highly evolvable expression profiles. In conclusion, by investigating the evolution of entire nervous systems at the resolution of single-neuron classes, we uncover patterns that may reflect basic principles governing evolutionary novelty in neuronal circuits.

Variation in inbreeding depression within and among Caenorhabditis species

Outbreeding populations harbor large numbers of recessive deleterious alleles that reduce the fitness of inbred individuals, and this inbreeding depression potentially shapes the evolution of mating systems, acting as a counterweight to the inherent selective advantage of self-fertilization. The population biological factors that influence inbreeding depression are numerous and often difficult to disentangle. We investigated the utility of obligately-outcrossing (gonochoristic) Caenorhabditis nematodes as models for inbreeding depression. By systematically inbreeding lines from ten populations and tracking line extinction, we found that inbreeding depression is universal but highly variable among species and populations. Inbreeding depression was detected across the life cycle, from mating to embryo production to embryonic viability and larval growth, and reciprocal crosses implicated female-biased effects. In most cases, the surviving inbred lines have dramatically reduced fitness, but the variance among inbred lines is substantial and compatible with the idea that inbreeding depression need not be an obstacle to the evolution of selfing in these worms. Populations of some species, including Caenorhabditis becei, exhibited modest inbreeding depression and could be tractable laboratory models for gonochoristic Caenorhabditis.

Molecular patterns of evolutionary changes throughout the whole nervous system of multiple nematode species

One avenue to better understand brain evolution is to map molecular patterns of evolutionary changes in neuronal cell types across entire nervous systems of distantly related species. Generating whole-animal single-cell transcriptomes of three nematode species from the Caenorhabditis genus, we observed a remarkable stability of neuronal cell type identities over more than 45 million years of evolution. Conserved patterns of combinatorial expression of homeodomain transcription factors are among the best classifiers of homologous neuron classes. Unexpectedly, we discover an extensive divergence in neuronal signaling pathways. While identities of neurotransmitter-producing neurons (glutamate, acetylcholine, GABA and several monoamines) remain stable, ionotropic and metabotropic receptors for all these neurotransmitter systems show substantial divergence, resulting in more than half of all neuron classes changing their capacity to be receptive to specific neurotransmitters. Neuropeptidergic signaling is also remarkably divergent, both at the level of neuropeptide expression and receptor expression, yet the overall dense network topology of the wireless neuropeptidergic connectome remains stable. Novel neuronal signaling pathways are suggested by our discovery of small secreted proteins that show no obvious hallmarks of conventional neuropeptides, but show similar patterns of highly neuron-type-specific and highly evolvable expression profiles. In conclusion, by investigating the evolution of entire nervous systems at the resolution of single neuron classes, we uncover patterns that may reflect basic principles governing evolutionary novelty in neuronal circuits.

Neurogenesis in Caenorhabditis elegans

Animals rely on their nervous systems to process sensory inputs, integrate these with internal signals, and produce behavioral outputs. This is enabled by the highly specialized morphologies and functions of neurons. Neuronal cells share multiple structural and physiological features, but they also come in a large diversity of types or classes that give the nervous system its broad range of functions and plasticity. This diversity, first recognized over a century ago, spurred classification efforts based on morphology, function, and molecular criteria. Caenorhabditis elegans, with its precisely mapped nervous system at the anatomical level, an extensive molecular description of most of its neurons, and its genetic amenability, has been a prime model for understanding how neurons develop and diversify at a mechanistic level. Here, we review the gene regulatory mechanisms driving neurogenesis and the diversification of neuron classes and subclasses in C. elegans. We discuss our current understanding of the specification of neuronal progenitors and their differentiation in terms of the transcription factors involved and ensuing changes in gene expression and chromatin landscape. The central theme that has emerged is that the identity of a neuron is defined by modules of gene batteries that are under control of parallel yet interconnected regulatory mechanisms. We focus on how, to achieve these terminal identities, cells integrate information along their developmental lineages. Moreover, we discuss how neurons are diversified postembryonically in a time-, genetic sex-, and activity-dependent manner. Finally, we discuss how the understanding of neuronal development can provide insights into the evolution of neuronal diversity.

Gene duplication and evolutionary plasticity of lin-12/Notch gene function in Caenorhabditis

Gene duplication is an important substrate for the evolution of new gene functions, but the impacts of gene duplicates on their own activities and on the developmental networks in which they act are poorly understood. Here, we use a natural experiment of lin-12/Notch gene duplication within the nematode genus Caenorhabditis, combined with characterization of loss- and gain-of-function mutations, to uncover functional distinctions between the duplicate genes in 1 species (Caenorhabditis briggsae) and their single-copy ortholog in Caenorhabditis elegans. First, using improved genomic sequence and gene model characterization, we confirm that the C. briggsae genome includes 2 complete lin-12 genes, whereas most other genes encoding proteins that participate in the LIN-12 signaling pathway retain a one-to-one orthology with C. elegans. We use CRISPR-mediated genome editing to introduce alleles predicted to cause gain-of-function (gf) or loss-of-function (lf) into each C. briggsae gene and find that the gf mutations uncover functional distinctions not apparent from the lf alleles. Specifically, Cbr-lin-12.1(gf), but not Cbr-lin-12.2(gf), causes developmental defects similar to those observed in Cel-lin-12(gf). In contrast to Cel-lin-12(gf), however, the Cbr-lin-12.1(gf) alleles do not cause dominant phenotypes as compared to the wild type, and the mutant phenotype is observed only when 2 gf alleles are present. Our results demonstrate that gene duplicates can exhibit differential capacities to compensate for each other and to interfere with normal development, and uncover coincident gene duplication and evolution of developmental sensitivity to LIN-12/Notch activity.

Modulating cell proliferation by asymmetric division: A conserved pattern in the early embryogenesis of nematode species

In the early stage of the nematode Caenorhabditis elegans embryogenesis, the zygote divides asymmetrically into a symmetric fast lineage and an asymmetric slow lineage, producing 16 and 8 cells respectively almost at the same time, followed by the onset of gastrulation. It was recently reported that this cell division pattern is optimal for rapid cell proliferation. In this work, we compare the cell lineages of 9 nematode species, revealing that this pattern is conserved for >60 million years. It further suggests that such lineage design has an important functional role and it might speed up embryonic development in the nematode kingdom, not limited to C. elegans , and independent of the maternal-zygotic transition dynamics.

Control of successive unequal cell divisions by neural cell fate regulators determines embryonic neuroblast cell size

Asymmetric cell divisions often generate daughter cells of unequal size in addition to different fates. In some contexts, daughter cell size asymmetry is thought to be a key input to specific binary cell fate decisions. An alternative possibility is that unequal division is a mechanism by which a variety of cells of different sizes are generated during embryonic development. We show here that two unequal cell divisions precede neuroblast formation in the C lineage of Caenorhabditis elegans. The equalisation of these divisions in a pig-1/MELK mutant background has little effect on neuroblast specification. Instead, we demonstrate that let-19/MDT13 is a regulator of the proneural basic helix-loop-helix transcription factor hlh-14/ASCL1 and find that both are required to concomitantly regulate the acquisition of neuroblast identity and neuroblast cell size. Thus, embryonic neuroblast cell size in this lineage is progressively regulated in parallel with identity by key neural cell fate regulators. We propose that key cell fate determinants have a previously unappreciated function in regulating unequal cleavage, and therefore cell size, of the progenitor cells whose daughter cell fates they then go on to specify.

Orthologs of the Caenorhabditis elegans heterochronic genes have divergent functions in Caenorhabditis briggsae

The heterochronic genes of Caenorhabditis elegans comprise the best-studied pathway controlling the timing of tissue and organ formation in an animal. To begin to understand the evolution of this pathway and the significance of the relationships among its components, we characterized 11 Caenorhabditis briggsae orthologs of C. elegans heterochronic genes. Using CRISPR/Cas9, we made a variety of alleles and found that several mutant phenotypes differ in significant ways from those of C. elegans. Although most mutant orthologs displayed defects in developmental timing, their phenotypes could differ in which stages were affected, the penetrance and expressivity of the phenotypes, or by having additional pleiotropies that were not obviously connected to developmental timing. However, when examining pairwise epistasis and synergistic relationships, we found those paralleled the known relationships between their C. elegans orthologs, suggesting that the arrangements of these genes in functional modules are conserved, but the modules' relationships to each other and/or to their targets has drifted since the time of the species' last common ancestor. Furthermore, our investigation has revealed a relationship between this pathway to other aspects of the animal's growth and development, including gonad development, which is relevant to both species.

Mechanisms of lineage specification in Caenorhabditis elegans

The studies of cell fate and lineage specification are fundamental to our understanding of the development of multicellular organisms. Caenorhabditis elegans has been one of the premiere systems for studying cell fate specification mechanisms at single cell resolution, due to its transparent nature, the invariant cell lineage, and fixed number of somatic cells. We discuss the general themes and regulatory mechanisms that have emerged from these studies, with a focus on somatic lineages and cell fates. We next review the key factors and pathways that regulate the specification of discrete cells and lineages during embryogenesis and postembryonic development; we focus on transcription factors and include numerous lineage diagrams that depict the expression of key factors that specify embryonic founder cells and postembryonic blast cells, and the diverse somatic cell fates they generate. We end by discussing some future perspectives in cell and lineage specification.

Evolutionary Change in Gut Specification in Caenorhabditis Centers on the GATA Factor ELT-3 in an Example of Developmental System Drift

Cells in a developing animal embryo become specified by the activation of cell-type-specific gene regulatory networks. The network that specifies the gut in the nematode Caenorhabditis elegans has been the subject of study for more than two decades. In this network, the maternal factors SKN-1/Nrf and POP-1/TCF activate a zygotic GATA factor cascade consisting of the regulators MED-1,2 → END-1,3 → ELT-2,7, leading to the specification of the gut in early embryos. Paradoxically, the MED, END, and ELT-7 regulators are present only in species closely related to C. elegans, raising the question of how the gut can be specified without them. Recent work found that ELT-3, a GATA factor without an endodermal role in C. elegans, acts in a simpler ELT-3 → ELT-2 network to specify gut in more distant species. The simpler ELT-3 → ELT-2 network may thus represent an ancestral pathway. In this review, we describe the elucidation of the gut specification network in C. elegans and related species and propose a model by which the more complex network might have formed. Because the evolution of this network occurred without a change in phenotype, it is an example of the phenomenon of Developmental System Drift.

Highly efficient transgenesis with miniMos in Caenorhabditis briggsae

Caenorhabditis briggsae as a companion species for Caenorhabditis elegans has played an increasingly important role in study of evolution of development and genome and gene regulation. Aided by the isolation of its sister spices, it has recently been established as a model for speciation study. To take full advantage of the species for comparative study, an effective transgenesis method especially those with single-copy insertion is important for functional comparison. Here, we improved a transposon-based transgenesis methodology that had been originally developed in C. elegans but worked marginally in C. briggsae. By incorporation of a heat shock step, the transgenesis efficiency in C. briggsae with a single-copy insertion is comparable to that in C. elegans. We used the method to generate 54 independent insertions mostly consisting of a mCherry tag over the C. briggsae genome. We demonstrated the use of the tags in identifying interacting loci responsible for hybrid male sterility between C. briggsae and Caenorhabditis nigoni when combined with the GFP tags we generated previously. Finally, we demonstrated that C. briggsae tolerates the C. elegans toxin, PEEL-1, but not SUP-35, making the latter a potential negative selection marker against extrachromosomal array.

Delineating the mechanisms and design principles of Caenorhabditis elegans embryogenesis using in toto high-resolution imaging data and computational modeling

The nematode (roundworm) Caenorhabditis elegans is one of the most popular animal models for the study of developmental biology, as its invariant development and transparent body enable in toto cellular-resolution fluorescence microscopy imaging of developmental processes at 1-min intervals. This has led to the development of various computational tools for the systematic and automated analysis of imaging data to delineate the molecular and cellular processes throughout the embryogenesis of C. elegans, such as those associated with cell lineage, cell migration, cell morphology, and gene activity. In this review, we first introduce C. elegans embryogenesis and the development of techniques for tracking cell lineage and reconstructing cell morphology during this process. We then contrast the developmental modes of C. elegans and the customized technologies used for studying them with the ones of other animal models, highlighting its advantage for studying embryogenesis with exceptional spatial and temporal resolution. This is followed by an examination of the physical models that have been devised-based on accurate determinations of developmental processes afforded by analyses of imaging data-to interpret the early embryonic development of C. elegans from subcellular to intercellular levels of multiple cells, which focus on two key processes: cell polarization and morphogenesis. We subsequently discuss how quantitative data-based theoretical modeling has improved our understanding of the mechanisms of C. elegans embryogenesis. We conclude by summarizing the challenges associated with the acquisition of C. elegans embryogenesis data, the construction of algorithms to analyze them, and the theoretical interpretation.

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.

Neuronal specification in C. elegans: combining lineage inheritance with intercellular signaling

The nervous system is composed of a high diversity of neuronal types. How this diversity is generated during development is a key question in neurobiology. Addressing this question is one of the reasons that led Sydney Brenner to develop the nematode C. elegans as a model organism. While there was initially a debate on whether the neuronal specification follows a 'European' model (determined by ancestry) or an 'American' model (determined by intercellular communication), several decades of research have established that the truth lies somewhere in between. Neurons are specified by the combination of transcription factors inherited from the ancestor cells and signaling between neighboring cells (especially Wnt and Notch signaling). This converges to the activation in newly generated postmitotic neurons of a specific set of terminal selector transcription factors that initiate and maintain the differentiation of the neuron. In this review, we also discuss the evolution of these specification mechanisms in other nematodes and beyond.

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.

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.

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.

Specific Interactions Between Autosome and X Chromosomes Cause Hybrid Male Sterility in Caenorhabditis Species

Hybrid male progeny from interspecies crosses are more prone to sterility or inviability than hybrid female progeny, and the male sterility and inviability often demonstrate parent-of-origin asymmetry. However, the underlying genetic mechanism of asymmetric sterility or inviability remains elusive. We previously established a genome-wide hybrid incompatibility (HI) landscape between Caenorhabditis briggsae and C. nigoni by phenotyping a large collection of C. nigoni strains each carrying a C. briggsae introgression. In this study, we systematically dissect the genetic mechanism of asymmetric sterility and inviability in both hybrid male and female progeny between the two species. Specifically, we performed reciprocal crosses between C . briggsae and different C. nigoni strains that each carry a GFP-labeled C. briggsae genomic fragment referred to as introgression, and scored the HI phenotypes in the F1 progeny. The aggregated introgressions cover 94.6% of the C. briggsae genome, including 100% of the X chromosome. Surprisingly, we observed that two C. briggsae X fragments that produce C. nigoni male sterility as an introgression rescued hybrid F1 sterility in males fathered by C. briggsae Subsequent backcrossing analyses indicated that a specific interaction between the X-linked interaction and one autosome introgression is required to rescue the hybrid male sterility. In addition, we identified another two C. briggsae genomic intervals on chromosomes II and IV that can rescue the inviability, but not the sterility, of hybrid F1 males fathered by C. nigoni, suggesting the involvement of differential epistatic interactions in the asymmetric hybrid male fertility and inviability. Importantly, backcrossing of the rescued sterile males with C. nigoni led to the isolation of a 1.1-Mb genomic interval that specifically interacts with an X-linked introgression, which is essential for hybrid male fertility. We further identified three C. briggsae genomic intervals on chromosome I, II, and III that produced inviability in all F1 progeny, dependent on or independent of the parent-of-origin. Taken together, we identified multiple independent interacting loci that are responsible for asymmetric hybrid male and female sterility, and inviability, which lays a foundation for their molecular characterization.

From "the Worm" to "the Worms" and Back Again: The Evolutionary Developmental Biology of Nematodes

Since the earliest days of research on nematodes, scientists have noted the developmental and morphological variation that exists within and between species. As various cellular and developmental processes were revealed through intense focus on Caenorhabditis elegans, these comparative studies have expanded. Within the genus Caenorhabditis, they include characterization of intraspecific polymorphisms and comparisons of distinct species, all generally amenable to the same laboratory culture methods and supported by robust genomic and experimental tools. The C. elegans paradigm has also motivated studies with more distantly related nematodes and animals. Combined with improved phylogenies, this work has led to important insights about the evolution of nematode development. First, while many aspects of C. elegans development are representative of Caenorhabditis, and of terrestrial nematodes more generally, others vary in ways both obvious and cryptic. Second, the system has revealed several clear examples of developmental flexibility in achieving a particular trait. This includes developmental system drift, in which the developmental control of homologous traits has diverged in different lineages, and cases of convergent evolution. Overall, the wealth of information and experimental techniques developed in C. elegans is being leveraged to make nematodes a powerful system for evolutionary cellular and developmental biology.

Dramatic evolution of body length due to postembryonic changes in cell size in a newly discovered close relative of Caenorhabditis elegans

Understanding morphological diversity-and morphological constraint-has been a central question in evolutionary biology since its inception. Nematodes of the genus Caenorhabditis, which contains the well-studied model organism C. elegans, display remarkable morphological consistency in the face of extensive genetic divergence. Here, we provide a description of the broad developmental patterns of a newly discovered species, C. sp. 34, which was isolated from fresh figs in Okinawa and which is among the closest known relatives of C. elegans. C. sp. 34 displays an extremely large body size; it can grow to be nearly twice as long as C. elegans and all other known members of the genus. Observations of the timing of developmental milestones reveal that C. sp. 34 develops about twice as slowly as C. elegans. Measurements of embryonic and larval size show that the size difference between C. sp. 34 and C. elegans is largely due to postembryonic events, particularly during the transition from larval to adult stages. This difference in size is not attributable to differences in germ line chromosome number or the number of somatic cells. The overall difference in body size is therefore largely attributable to changes in cell size via increased cytoplasmic volume. Because of its close relationship to C. elegans, the distinctness of C. sp. 34 provides an ideal system for the detailed analysis of evolutionary diversification. The context of over 40 years of C. elegans developmental genetics also reveals clues into how natural selection and developmental constraint act jointly to promote patterns of morphological stasis and divergence in this group.

Establishment of Signaling Interactions with Cellular Resolution for Every Cell Cycle of Embryogenesis

Intercellular signaling interactions play a key role in breaking fate symmetry during animal development. Identification of signaling interactions at cellular resolution is technically challenging, especially in a developing embryo. Here, we develop a platform that allows automated inference and validation of signaling interactions for every cell cycle of Caenorhabditis elegans embryogenesis. This is achieved by the generation of a systems-level cell contact map, which consists of 1114 highly confident intercellular contacts, by modeling analysis and is validated through cell membrane labeling coupled with cell lineage analysis. We apply the map to identify cell pairs between which a Notch signaling interaction takes place. By generating expression patterns for two ligands and two receptors of the Notch signaling pathway with cellular resolution using the automated expression profiling technique, we are able to refine existing and identify novel Notch interactions during C. elegans embryogenesis. Targeted cell ablation followed by cell lineage analysis demonstrates the roles of signaling interactions during cell division in breaking fate symmetry. Finally, we describe the development of a website that allows online access to the cell-cell contact map for mapping of other signaling interactions by the community. The platform can be adapted to establish cellular interactions from any other signaling pathway.

Comparative proteome analysis between C . briggsae embryos and larvae reveals a role of chromatin modification proteins in embryonic cell division

Caenorhabditis briggsae has emerged as a model for comparative biology against model organism C. elegans. Most of its cell fate specifications are completed during embryogenesis whereas its cell growth is achieved mainly in larval stages. The molecular mechanism underlying the drastic developmental changes is poorly understood. To gain insights into the molecular changes between the two stages, we compared the proteomes between the two stages using iTRAQ. We identified a total of 2,791 proteins in the C. briggsae embryos and larvae, 247 of which undergo up- or down-regulation between the two stages. The proteins that are upregulated in the larval stages are enriched in the Gene Ontology categories of energy production, protein translation, and cytoskeleton; whereas those upregulated in the embryonic stage are enriched in the categories of chromatin dynamics and posttranslational modification, suggesting a more active chromatin modification in the embryos than in the larva. Perturbation of a subset of chromatin modifiers followed by cell lineage analysis suggests their roles in controlling cell division pace. Taken together, we demonstrate a general molecular switch from chromatin modification to metabolism during the transition from C. briggsae embryonic to its larval stages using iTRAQ approach. The switch might be conserved across metazoans.

Differing roles for sur-2/MED23 in C. elegans and C. briggsae vulval development

Normal vulval development in the nematode Caenorhabditis briggsae is identical to that in the related Caenorhabditis elegans. However, several experiments suggest that there are differences between the two species with respect to the contribution of EGF/Ras signaling. To investigate these differences genetically, we have characterized a C. briggsae mutant strain that phenocopies the effect observed when C. briggsae animals are treated with U0126, an inhibitor of the EGF pathway component MEK. We identify that the gene affected in the mutant strain is Cbr-sur-2, which encodes a MED23 mediator complex protein that acts downstream of EGF signaling in C. elegans and other organisms, such as mammals. When Cbr-sur-2 and Cel-sur-2 mutants are compared, we find that the production of additional vulval cells from P5.p and P7.p in C. elegans is dependent on proper development of P6.p, while C. briggsae does not have a similar requirement. Combined chemical and genetic interference with the EGF pathway completely eliminates vulval development in C. elegans but not in C. briggsae. Our results provide genetic evidence for the differing requirements for EGF signaling in the two species.

Conservation of anatomically restricted glycosaminoglycan structures in divergent nematode species

Heparan sulfates (HS) are glycosaminoglycans of the extracellular matrices and characterized by complex modification patterns owing to sulfations, epimerization, and acetylation. Distinct HS modification patterns have been shown to modulate protein-protein interactions during development in general and of the nervous system in particular. This has led to the heparan sulfate code hypothesis, which posits that specifically modified HS epitopes are distributed in a tissue and cell-specific fashion to orchestrate neural circuit formation. Whether an HS code exists in vivo, how specific or how evolutionarily conserved the anatomical distribution of an HS code may be has remained unknown. Here we conduct a systematic comparison of HS modification patterns in the nematode Caenorhabditis elegans using transgenic expression of 33 different HS-specific single chain variable fragment antibodies. We find that some HS modification patterns are widely distributed in the nervous system. In contrast, other HS modification patterns appear highly cell-specific in both non-neuronal and neuronal cells. Some patterns can be as restricted in their localization as to single neurites or synaptic connections between two neurons. This restricted anatomical localization of specific HS patterns can be evolutionarily conserved over a span of 80-100 million years in the divergent nematode species Caenorhabditis briggsae suggesting structural and, possibly functional conservation of glycosaminoglycan structures similar to proteins. These findings suggest a HS code with subcellularly localized, unique glycan identities in the nervous system.

SSBD: a database of quantitative data of spatiotemporal dynamics of biological phenomena

Motivation: Rapid advances in live-cell imaging analysis and mathematical modeling have produced a large amount of quantitative data on spatiotemporal dynamics of biological objects ranging from molecules to organisms. There is now a crucial need to bring these large amounts of quantitative biological dynamics data together centrally in a coherent and systematic manner. This will facilitate the reuse of this data for further analysis.

Results: We have developed the Systems Science of Biological Dynamics database (SSBD) to store and share quantitative biological dynamics data. SSBD currently provides 311 sets of quantitative data for single molecules, nuclei and whole organisms in a wide variety of model organisms from Escherichia coli to Mus musculus. The data are provided in Biological Dynamics Markup Language format and also through a REST API. In addition, SSBD provides 188 sets of time-lapse microscopy images from which the quantitative data were obtained and software tools for data visualization and analysis.

Availability and Implementation: SSBD is accessible at http://ssbd.qbic.riken.jp.

Contact:sonami@riken.jp

Developmental robustness in the Caenorhabditis elegans embryo

Developmental robustness is the ability of an embryo to develop normally despite many sources of variation, from differences in the environment to stochastic cell-to-cell differences in gene expression. The nematode Caenorhabditis elegans exhibits an additional level of robustness: Unlike most other animals, the embryonic pattern of cell divisions is nearly identical from animal to animal. The endoderm (gut) lineage is an ideal model for studying such robustness as the juvenile gut has a simple anatomy, consisting of 20 cells that are derived from a single cell, E, and the gene regulatory network that controls E specification shares features with developmental regulatory networks in many other systems, including genetic redundancy, parallel pathways, and feed-forward loops. Early studies were initially concerned with identifying the genes in the network, whereas recent work has focused on understanding how the endoderm produces a robust developmental output in the face of many sources of variation. Genetic control exists at three levels of endoderm development: Progenitor specification, cell divisions within the developing gut, and maintenance of gut differentiation. Recent findings show that specification genes regulate all three of these aspects of gut development, and that mutant embryos can experience a "partial" specification state in which some, but not all, E descendants adopt a gut fate. Ongoing studies using newer quantitative and genome-wide methods promise further insights into how developmental gene-regulatory networks buffer variation.

The Homeobox Genes of Caenorhabditis elegans and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis

Homeobox genes play crucial roles for the development of multicellular eukaryotes. We have generated a revised list of all homeobox genes for Caenorhabditis elegans and provide a nomenclature for the previously unnamed ones. We show that, out of 103 homeobox genes, 70 are co-orthologous to human homeobox genes. 14 are highly divergent, lacking an obvious ortholog even in other Caenorhabditis species. One of these homeobox genes encodes 12 homeodomains, while three other highly divergent homeobox genes encode a novel type of double homeodomain, termed HOCHOB. To understand how transcription factors regulate cell fate during development, precise spatio-temporal expression data need to be obtained. Using a new imaging framework that we developed, Endrov, we have generated spatio-temporal expression profiles during embryogenesis of over 60 homeobox genes, as well as a number of other developmental control genes using GFP reporters. We used dynamic feedback during recording to automatically adjust the camera exposure time in order to increase the dynamic range beyond the limitations of the camera. We have applied the new framework to examine homeobox gene expression patterns and provide an analysis of these patterns. The methods we developed to analyze and quantify expression data are not only suitable for C. elegans, but can be applied to other model systems or even to tissue culture systems.

Full-genome evolutionary histories of selfing, splitting, and selection in Caenorhabditis

The nematode Caenorhabditis briggsae is a model for comparative developmental evolution with C. elegans. Worldwide collections of C. briggsae have implicated an intriguing history of divergence among genetic groups separated by latitude, or by restricted geography, that is being exploited to dissect the genetic basis to adaptive evolution and reproductive incompatibility; yet, the genomic scope and timing of population divergence is unclear. We performed high-coverage whole-genome sequencing of 37 wild isolates of the nematode C. briggsae and applied a pairwise sequentially Markovian coalescent (PSMC) model to 703 combinations of genomic haplotypes to draw inferences about population history, the genomic scope of natural selection, and to compare with 40 wild isolates of C. elegans. We estimate that a diaspora of at least six distinct C. briggsae lineages separated from one another approximately 200,000 generations ago, including the “Temperate” and “Tropical” phylogeographic groups that dominate most samples worldwide. Moreover, an ancient population split in its history approximately 2 million generations ago, coupled with only rare gene flow among lineage groups, validates this system as a model for incipient speciation. Low versus high recombination regions of the genome give distinct signatures of population size change through time, indicative of widespread effects of selection on highly linked portions of the genome owing to extreme inbreeding by self-fertilization. Analysis of functional mutations indicates that genomic context, owing to selection that acts on long linkage blocks, is a more important driver of population variation than are the functional attributes of the individually encoded genes.

Comprehensive single cell-resolution analysis of the role of chromatin regulators in early C. elegans embryogenesis

Chromatin regulators are widely expressed proteins with diverse roles in gene expression, nuclear organization, cell cycle regulation, pluripotency, physiology and development, and are frequently mutated in human diseases such as cancer. Their inhibition often results in pleiotropic effects that are difficult to study using conventional approaches. We have developed a semi-automated nuclear tracking algorithm to quantify the divisions, movements and positions of all nuclei during the early development of Caenorhabditis elegans and have used it to systematically study the effects of inhibiting chromatin regulators. The resulting high dimensional datasets revealed that inhibition of multiple regulators, including F55A3.3 (encoding FACT subunit SUPT16H), lin-53 (RBBP4/7), rba-1 (RBBP4/7), set-16 (MLL2/3), hda-1 (HDAC1/2), swsn-7 (ARID2), and let-526 (ARID1A/1B) affected cell cycle progression and caused chromosome segregation defects. In contrast, inhibition of cir-1 (CIR1) accelerated cell division timing in specific cells of the AB lineage. The inhibition of RNA polymerase II also accelerated these division timings, suggesting that normal gene expression is required to delay cell cycle progression in multiple lineages in the early embryo. Quantitative analyses of the dataset suggested the existence of at least two functionally distinct SWI/SNF chromatin remodeling complex activities in the early embryo, and identified a redundant requirement for the egl-27 and lin-40 MTA orthologs in the development of endoderm and mesoderm lineages. Moreover, our dataset also revealed a characteristic rearrangement of chromatin to the nuclear periphery upon the inhibition of multiple general regulators of gene expression. Our systematic, comprehensive and quantitative datasets illustrate the power of single cell-resolution quantitative tracking and high dimensional phenotyping to investigate gene function. Furthermore, the results provide an overview of the functions of essential chromatin regulators during the early development of an animal.

Gonad morphogenesis defects drive hybrid male sterility in asymmetric hybrid breakdown of Caenorhabditis nematodes

Determining the causes and evolution of reproductive barriers to gene flow between populations, speciation, is the key to understanding the origin of diversity in nature. Many species manifest hybrid breakdown when they intercross, characterized by increasingly exacerbated problems in later generations of hybrids. Recently, Caenorhabditis nematodes have emerged as a genetic model for studying speciation, and here we investigate the nature and causes of hybrid breakdown between Caenorhabditis remanei and C. latens. We quantify partial F1 hybrid inviability and extensive F2 hybrid inviability; the ~75% F2 embryonic arrest occurs primarily during gastrulation or embryonic elongation. Moreover, F1 hybrid males exhibit Haldane's rule asymmetrically for both sterility and inviability, being strongest when C. remanei serves as maternal parent. We show that the mechanism by which sterile hybrid males are incapable of transferring sperm or a copulatory plug involves defective gonad morphogenesis, which we hypothesize results from linker cell defects in migration and/or cell death during development. This first documented case of partial hybrid male sterility in Caenorhabditis follows expectations of Darwin's corollary to Haldane's rule for asymmetric male fitness, providing a powerful foundation for molecular dissection of intrinsic reproductive barriers and divergence of genetic pathways controlling organ morphogenesis.

Pervasive divergence of transcriptional gene regulation in Caenorhabditis nematodes

Because there is considerable variation in gene expression even between closely related species, it is clear that gene regulatory mechanisms evolve relatively rapidly. Because primary sequence conservation is an unreliable proxy for functional conservation of cis-regulatory elements, their assessment must be carried out in vivo. We conducted a survey of cis-regulatory conservation between C. elegans and closely related species C. briggsae, C. remanei, C. brenneri, and C. japonica. We tested enhancers of eight genes from these species by introducing them into C. elegans and analyzing the expression patterns they drove. Our results support several notable conclusions. Most exogenous cis elements direct expression in the same cells as their C. elegans orthologs, confirming gross conservation of regulatory mechanisms. However, the majority of exogenous elements, when placed in C. elegans, also directed expression in cells outside endogenous patterns, suggesting functional divergence. Recurrent ectopic expression of different promoters in the same C. elegans cells may reflect biases in the directions in which expression patterns can evolve due to shared regulatory logic of coexpressed genes. The fact that, despite differences between individual genes, several patterns repeatedly emerged from our survey, encourages us to think that general rules governing regulatory evolution may exist and be discoverable.

Comparative RNAi screens in C. elegans and C. briggsae reveal the impact of developmental system drift on gene function

Although two related species may have extremely similar phenotypes, the genetic networks underpinning this conserved biology may have diverged substantially since they last shared a common ancestor. This is termed Developmental System Drift (DSD) and reflects the plasticity of genetic networks. One consequence of DSD is that some orthologous genes will have evolved different in vivo functions in two such phenotypically similar, related species and will therefore have different loss of function phenotypes. Here we report an RNAi screen in C. elegans and C. briggsae to identify such cases. We screened 1333 genes in both species and identified 91 orthologues that have different RNAi phenotypes. Intriguingly, we find that recently evolved genes of unknown function have the fastest evolving in vivo functions and, in several cases, we identify the molecular events driving these changes. We thus find that DSD has a major impact on the evolution of gene function and we anticipate that the C. briggsae RNAi library reported here will drive future studies on comparative functional genomics screens in these nematodes.

Automated cellular annotation for high-resolution images of adult Caenorhabditis elegans

Motivation: Advances in high-resolution microscopy have recently made possible the analysis of gene expression at the level of individual cells. The fixed lineage of cells in the adult worm Caenorhabditis elegans makes this organism an ideal model for studying complex biological processes like development and aging. However, annotating individual cells in images of adult C.elegans typically requires expertise and significant manual effort. Automation of this task is therefore critical to enabling high-resolution studies of a large number of genes.

Results: In this article, we describe an automated method for annotating a subset of 154 cells (including various muscle, intestinal and hypodermal cells) in high-resolution images of adult C.elegans. We formulate the task of labeling cells within an image as a combinatorial optimization problem, where the goal is to minimize a scoring function that compares cells in a test input image with cells from a training atlas of manually annotated worms according to various spatial and morphological characteristics. We propose an approach for solving this problem based on reduction to minimum-cost maximum-flow and apply a cross-entropy–based learning algorithm to tune the weights of our scoring function. We achieve 84% median accuracy across a set of 154 cell labels in this highly variable system. These results demonstrate the feasibility of the automatic annotation of microscopy-based images in adult C.elegans.

Contact:saerni@cs.stanford.edu

Boolean genetic network model for the control of C. elegans early embryonic cell cycles

Background

In Caenorhabditis elegans early embryo, cell cycles only have two phases: DNA synthesis and mitosis, which are different from the typical 4-phase cell cycle. Modeling this cell-cycle process into network can fill up the gap in C. elegans cell-cycle study and provide a thorough understanding on the cell-cycle regulations and progressions at the network level.

Methods

In this paper, C. elegans early embryonic cell-cycle network has been constructed based on the knowledge of key regulators and their interactions from literature studies. A discrete dynamical Boolean model has been applied in computer simulations to study dynamical properties of this network. The cell-cycle network is compared with random networks and tested under several perturbations to analyze its robustness. To investigate whether our proposed network could explain biological experiment results, we have also compared the network simulation results with gene knock down experiment data.

Results

With the Boolean model, this study showed that the cell-cycle network was stable with a set of attractors (fixed points). A biological pathway was observed in the simulation, which corresponded to a whole cell-cycle progression. The C. elegans network was significantly robust when compared with random networks of the same size because there were less attractors and larger basins than random networks. Moreover, the network was also robust under perturbations with no significant change of the basin size. In addition, the smaller number of attractors and the shorter biological pathway from gene knock down network simulation interpreted the shorter cell-cycle lengths in mutant from the RNAi gene knock down experiment data. Hence, we demonstrated that the results in network simulation could be verified by the RNAi gene knock down experiment data.

Conclusions

A C. elegans early embryonic cell cycles network was constructed and its properties were analyzed and compared with those of random networks. Computer simulation results provided biologically meaningful interpretations of RNAi gene knock down experiment data.

Evolutionarily dynamic roles of a PUF RNA-binding protein in the somatic development of Caenorhabditis briggsae

Gene duplication and divergence has emerged as an important aspect of developmental evolution. The genomes of Caenorhabditis nematodes encode an ancient family of PUF RNA-binding proteins. Most have been implicated in germline development, and are often redundant with paralogs of the same sub-family. An exception is Cbr-puf-2 (one of three Caenorhabditis briggsae PUF-2 sub-family paralogs), which is required for development past the second larval stage. Here, we provide a detailed functional characterization of Cbr-puf-2. The larval arrest of Cbr-puf-2 mutant animals is caused by inefficient breakdown of bacterial food, which leads to starvation. Cbr-puf-2 is required for the normal grinding cycle of the muscular terminal bulb during early larval stages, and is transiently expressed in this tissue. In addition, rescue of larval arrest reveals that Cbr-puf-2 also promotes normal vulval development. It is expressed in the anchor cell (which induces vulval fate) and vulval muscles, but not in the vulva precursor cells (VPCs) themselves. This contrasts with the VPC-autonomous repression of vulval development described for the Caenorhabditis elegans homologs fbf-1/2. These different roles for PUF proteins occur even as the vulva and pharynx maintain highly conserved anatomies across Caenorhabditis, indicating pervasive developmental system drift (DSD). Because Cbr-PUF-2 shares RNA-binding specificity with its paralogs and with C. elegans FBF, we suggest that functional novelty of RNA-binding proteins evolves through changes in the site of their expression, perhaps in concert with cis-regulatory evolution in target mRNAs.

Genetic control of vulval development in Caenorhabditis briggsae

The nematode Caenorhabditis briggsae is an excellent model organism for the comparative analysis of gene function and developmental mechanisms. To study the evolutionary conservation and divergence of genetic pathways mediating vulva formation, we screened for mutations in C. briggsae that cause the egg-laying defective (Egl) phenotype. Here, we report the characterization of 13 genes, including three that are orthologs of Caenorhabditis elegans unc-84 (SUN domain), lin-39 (Dfd/Scr-related homeobox), and lin-11 (LIM homeobox). Based on the morphology and cell fate changes, the mutants were placed into four different categories. Class 1 animals have normal-looking vulva and vulva-uterine connections, indicating defects in other components of the egg-laying system. Class 2 animals frequently lack some or all of the vulval precursor cells (VPCs) due to defects in the migration of P-cell nuclei into the ventral hypodermal region. Class 3 animals show inappropriate fusion of VPCs to the hypodermal syncytium, leading to a reduced number of vulval progeny. Finally, class 4 animals exhibit abnormal vulval invagination and morphology. Interestingly, we did not find mutations that affect VPC induction and fates. Our work is the first study involving the characterization of genes in C. briggsae vulva formation, and it offers a basis for future investigations of these genes in C. elegans.

Coevolution within and between regulatory loci can preserve promoter function despite evolutionary rate acceleration

Phenotypes that appear to be conserved could be maintained not only by strong purifying selection on the underlying genetic systems, but also by stabilizing selection acting via compensatory mutations with balanced effects. Such coevolution has been invoked to explain experimental results, but has rarely been the focus of study. Conserved expression driven by the unc-47 promoters of Caenorhabditis elegans and C. briggsae persists despite divergence within a cis-regulatory element and between this element and the trans-regulatory environment. Compensatory changes in cis and trans are revealed when these promoters are used to drive expression in the other species. Functional changes in the C. briggsae promoter, which has experienced accelerated sequence evolution, did not lead to alteration of gene expression in its endogenous environment. Coevolution among promoter elements suggests that complex epistatic interactions within cis-regulatory elements may facilitate their divergence. Our results offer a detailed picture of regulatory evolution in which subtle, lineage-specific, and compensatory modifications of interacting cis and trans regulators together maintain conserved gene expression patterns.

Some phenotypes, including gene expression patterns, are conserved between distantly related species. However, the molecular bases of those phenotypes are not necessarily conserved. Instead, regulatory DNA sequences and the proteins with which they interact can change over time with balanced effects, preserving expression patterns and concealing regulatory divergence. Coevolution between interacting molecules makes gene regulation highly species-specific, and it can be detected when the cis-regulatory DNA of one species is used to drive expression in another species. In this way, we identified regions of the C. elegans and C. briggsae unc-47 promoters that have coevolved with the lineage-specific trans-regulatory environments of these organisms. The C. briggsae promoter experienced accelerated sequence change relative to related species. All of this evolution occurred without changing the expression pattern driven by the promoter in its endogenous environment.

A method for rapid and simultaneous mapping of genetic loci and introgression sizes in nematode species

Caenorhabditis briggsae is emerging as an attractive model organism not only in studying comparative biology against C. elegans, but also in developing novel experimentation avenues. In particular, recent identification of a new Caenorhabditis species, C. sp.9 with which it can mate and produce viable progeny provides an opportunity for studying the genetics of hybrid incompatibilities (HI) between the two. Mapping of a specific HI locus demands repeated backcrossing to get hold of the specific genomic region underlying an observed phenotype. To facilitate mapping of HI loci between C. briggsae and C. sp.9, an efficient mapping method and a genetic map ideally consisting of dominant markers are required for systematic introgression of genomic fragments between the two species. We developed a fast and cost-effective method for high throughput mapping of dominant loci with resolution up to 1 million bps in C. briggsae. The method takes advantage of the introgression between C. briggsae and C. sp.9 followed by PCR genotyping using C. briggsae specific primers. Importantly, the mapping results can not only serve as an effective way for estimating the chromosomal position of a genetic locus in C. briggsae, but also provides size information for the introgression fragment in an otherwise C. sp.9 background. In addition, it also helps generate introgression line as a side-product that is invaluable for the subsequent mapping of HI loci. The method will greatly facilitate the construction of a genetic map consisting of dominant markers and pave the way for systematic isolation of HI loci between C. briggsae and C. sp.9 which has so far not been attempted between nematode species. The method is designed for mapping of a dominant allele, but can be easily adapted for mapping of any other type of alleles in any other species if introgression between a sister species pair is feasible.

RNA-seq analysis of the C. briggsae transcriptome

Curation of a high-quality gene set is the critical first step in genome research, enabling subsequent analyses such as ortholog assignment, cis-regulatory element finding, and synteny detection. In this project, we have reannotated the genome of Caenorhabditis briggsae, the best studied sister species of the model organism Caenorhabditis elegans. First, we applied a homology-based gene predictor genBlastG to annotate the C. briggsae genome. We then validated and further improved the C. briggsae gene annotation through RNA-seq analysis of the C. briggsae transcriptome, which resulted in the first validated C. briggsae gene set (23,159 genes), among which 7347 genes (33.9% of all genes with introns) have all of their introns confirmed. Most genes (14,812, or 68.3%) have at least one intron validated, compared with only 3.9% in the most recent WormBase release (WS228). Of all introns in the revised gene set (103,083), 61,503 (60.1%) have been confirmed. Additionally, we have identified numerous trans-splicing leaders (SL1 and SL2 variants) in C. briggsae, leading to the first genome-wide annotation of operons in C. briggsae (1105 operons). The majority of the annotated operons (564, or 51.0%) are perfectly conserved in C. elegans, with an additional 345 operons (or 31.2%) somewhat divergent. Additionally, RNA-seq analysis revealed over 10 thousand small-size assembly errors in the current C. briggsae reference genome that can be readily corrected. The revised C. briggsae genome annotation represents a solid platform for comparative genomics analysis and evolutionary studies of Caenorhabditis species.

Plectus - a stepping stone in embryonic cell lineage evolution of nematodes

BACKGROUND

Recent studies have challenged the widespread view that the pattern of embryogenesis found in Caenorhabditis elegans (clade 9) is characteristic of nematodes in general. To understand this still largely unexplored landscape of developmental events, we set out to examine more distantly related nematodes in detail for temporospatial differences in pattern formation and cell specification. Members of the genus Plectus (clade 6) seem to be suitable candidates to show variety, with certain idiosyncratic features during early development and the convenient availability of cultivatable species.

METHODS

The study was conducted using 4-D lineage analysis, 3-D modeling of developing embryos and laser-induced ablation of individual blastomeres.

RESULTS

Detailed cell lineage studies of several Plectus species reveal that pattern formation and cell fate assignment differ markedly from C. elegans. Descendants of the first somatic founder cell S1 (AB) - but not the progeny of other founder cells - demonstrate extremely variable spatial arrangements illustrating that here distinct early cell-cell interactions between invariant partners, as found in C. elegans, cannot take place. Different from C. elegans, in Plectus alternative positional variations among early S1 blastomeres resulting in a 'situs inversus' pattern, nevertheless give rise to adults with normal left-right asymmetries. In addition, laser ablations of early blastomeres uncover inductions between variable cell partners.

CONCLUSIONS

Our results suggest that embryonic cell specification in Plectus is not correlated with cell lineage but with position. With this peculiarity, Plectus appears to occupy an intermediate position between basal nematodes displaying a variable early development and the C. elegans-like invariant pattern. We suggest that indeterminate pattern formation associated with late, position-dependent fate assignment represents a plesiomorphic character among nematodes predominant in certain basal clades but lost in derived clades. Thus, the behavior of S1 cells in Plectus can be considered an evolutionary relict in a transition phase between two different developmental strategies.

Developmental milestones punctuate gene expression in the Caenorhabditis embryo

A fundamental question in developmental biology relates to the connection between morphological stages and their underlying molecular activity. Here we demonstrate that, at the molecular level, embryonic development in five Caenorhabditis species proceeds through two distinct milestones in which the transcriptome is resistant to differences in species-specific developmental timings. By comparing the complete protein-coding transcriptomes of individually timed embryos across ten morphological markers, we found that developmental milestones can be characterized by their expression dynamics and activation of key developmental regulators. This approach led us to discover the nematode phylotypic stage and to show that in chordates and arthropods it is represented as two distinct stages, suggesting that animal body plans might evolve by uncoupling and elaboration on formerly synchronous processes.

Evolution of embryonic development in nematodes

BACKGROUND

Nematodes can be subdivided into basal Enoplea (clades 1 and 2) and more derived Chromadorea (clades 3 to 12). Embryogenesis of Caenorhabditis elegans (clade 9) has been analyzed in most detail. Their establishment of polarity and asymmetric cleavage requires the differential localization of PAR proteins. Earlier studies on selected other nematodes revealed that embryonic development of nematodes is more diverse than the essentially invariant development of C. elegans and the classic study object Ascaris had suggested. To obtain a more detailed picture of variations and evolutionary trends we compared embryonic cell lineages and pattern formation in embryos of all 12 nematode clades.

METHODS

The study was conducted using 4-D microscopy and 3-D modeling of developing embryos.

RESULTS

We found dramatic differences compared to C. elegans in Enoplea but also considerable variations among Chromadorea. We discovered 'Polarity Organizing Centers' (POCs) that orient cleavage spindles along the anterior-posterior axis in distinct cells over consecutive cell generations. The resulting lineally arranged blastomeres represent a starting point for the establishment of bilateral symmetry within individual lineages. We can discern six different early cleavage types and suggest that these variations are due to modifications in the activity of the POCs in conjunction with changes in the distribution of PAR proteins. In addition, our studies indicate that lineage complexity advanced considerably during evolution, that is we observe trends towards an increase of somatic founder cells, from monoclonal to polyclonal lineages and from a variable (position-dependent) to an invariable (lineage-dependent) way of cell fate specification. In contrast to the early phase of embryogenesis, the second half ('morphogenesis') appears similar in all studied nematodes. Comparison of early cleavage between the basal nematode Tobrilus stefanskii and the tardigrade Hypsibius dujardini revealed surprising similarities indicating that the presence of POCs is not restricted to nematode embryos.

CONCLUSIONS

The pattern of cleavage, spatial arrangement and differentiation of cells diverged dramatically during the history of the phylum Nematoda without corresponding changes in the phenotype. While in all studied representatives the same distinctive developmental steps need to be taken, cell behavior leading to these is not conserved.

Internalization of multiple cells during C. elegans gastrulation depends on common cytoskeletal mechanisms but different cell polarity and cell fate regulators

Understanding the links between developmental patterning mechanisms and force-producing cytoskeletal mechanisms is a central goal in studies of morphogenesis. Gastrulation is the first morphogenetic event in the development of many organisms. Gastrulation involves the internalization of surface cells, often driven by the contraction of actomyosin networks that are deployed with spatial precision-both in specific cells and in a polarized manner within each cell. These cytoskeletal mechanisms rely on different cell fate and cell polarity regulators in different organisms. Caenorhabditis elegans gastrulation presents an opportunity to examine the extent to which diverse mechanisms may be used by dozens of cells that are internalized at distinct times within a single organism. We identified 66 cells that are internalized in C. elegans gastrulation, many of which were not known previously to gastrulate. To gain mechanistic insights into how these cells internalize, we genetically manipulated cell fate, cell polarity and cytoskeletal regulators and determined the effects on cell internalization. We found that cells of distinct lineages depend on common actomyosin-based mechanisms to gastrulate, but different cell fate regulators, and, surprisingly, different cell polarity regulators. We conclude that diverse cell fate and cell polarity regulators control common mechanisms of morphogenesis in C. elegans. The results highlight the variety of developmental patterning mechanisms that can be associated with common cytoskeletal mechanisms in the morphogenesis of an animal embryo.

Natural variation in life history and aging phenotypes is associated with mitochondrial DNA deletion frequency in Caenorhabditis briggsae

BACKGROUND

Mutations that impair mitochondrial functioning are associated with a variety of metabolic and age-related disorders. A barrier to rigorous tests of the role of mitochondrial dysfunction in aging processes has been the lack of model systems with relevant, naturally occurring mitochondrial genetic variation. Toward the goal of developing such a model system, we studied natural variation in life history, metabolic, and aging phenotypes as it relates to levels of a naturally-occurring heteroplasmic mitochondrial ND5 deletion recently discovered to segregate among wild populations of the soil nematode, Caenorhabditis briggsae. The normal product of ND5 is a central component of the mitochondrial electron transport chain and integral to cellular energy metabolism.

RESULTS

We quantified significant variation among C. briggsae isolates for all phenotypes measured, only some of which was statistically associated with isolate-specific ND5 deletion frequency. We found that fecundity-related traits and pharyngeal pumping rate were strongly inversely related to ND5 deletion level and that C. briggsae isolates with high ND5 deletion levels experienced a tradeoff between early fecundity and lifespan. Conversely, oxidative stress resistance was only weakly associated with ND5 deletion level while ATP content was unrelated to deletion level. Finally, mean levels of reactive oxygen species measured in vivo showed a significant non-linear relationship with ND5 deletion level, a pattern that may be driven by among-isolate variation in antioxidant or other compensatory mechanisms.

CONCLUSIONS

Our findings suggest that the ND5 deletion may adversely affect fitness and mitochondrial functioning while promoting aging in natural populations, and help to further establish this species as a useful model for explicit tests of hypotheses in aging biology and mitochondrial genetics.

Cell identification and cell lineage analysis

Caenorhabditis elegans is uniquely suited to the analysis of cell lineage patterns. C. elegans has a small number of somatic cells whose position and morphology are almost invariant from animal to animal. Because C. elegans is virtually transparent, cells can be identified in live animals using a simple bright-field microscopy technique, Nomarski differential interference contrast (DIC), or by expression of transgenic fluorescent reporter genes. The small size and rapid development of C. elegans mean that animals can develop while under continuous observation, allowing cell lineages to be analyzed throughout embryonic and postembryonic development. Embryonic cell lineages can also be traced semiautomatically using timelapse imaging of GFP-labeled nuclei. Analysis of mutant cell lineages remains important for defining the roles of developmental control genes.