Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes.

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Post-insemination selection dominates pre-insemination selection in driving rapid evolution of male competitive ability.

Sexual reproduction is a complex process that contributes to differences between the sexes and divergence between species. From a male's perspective, sexual selection can optimize reproductive success by acting on the variance in mating success (pre-insemination selection) as well as the variance in fertilization success (post-insemination selection). The balance between pre- and post-insemination selection has not yet been investigated using a strong hypothesis-testing framework that directly quantifies the effects of post-insemination selection on the evolution of reproductive success. Here we use experimental evolution of a uniquely engineered genetic system that allows sperm production to be turned off and on in obligate male-female populations of Caenorhabditis elegans. We show that enhanced post-insemination competition increases the efficacy of selection and surpasses pre-insemination sexual selection in driving a polygenic response in male reproductive success. We find that after 10 selective events occurring over 30 generations post-insemination selection increased male reproductive success by an average of 5- to 7-fold. Contrary to expectation, enhanced pre-insemination competition hindered selection and slowed the rate of evolution. Furthermore, we found that post-insemination selection resulted in a strong polygenic response at the whole-genome level. Our results demonstrate that post-insemination sexual selection plays a critical role in the rapid optimization of male reproductive fitness. Therefore, explicit consideration should be given to post-insemination dynamics when considering the population effects of sexual selection.

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations.

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.

A large close relative of C. elegans is slow-developing but not long-lived.

BACKGROUND: Variation in body size is thought to be a major driver of a wide variety of ecological and evolutionary patterns, including changes in development, reproduction, and longevity. Additionally, drastic changes in natural context often have profound effects on multiple fitness-related traits. Caenorhabditis inopinata is a recently-discovered fig-associated nematode that is unusually large relative to other members of the genus, including the closely related model system C. elegans. Here we test whether the dramatic increase in body size and shift in ecological context has led to correlated changes in key life history and developmental parameters within this species. RESULTS: Using four developmental milestones, C. inopinata was found to have a slower rate of development than C. elegans across a range of temperatures. Despite this, C. inopinata did not reveal any differences in adult lifespan from C. elegans after accounting for differences in developmental timing and reproductive mode. C. inopinata fecundity was generally lower than that of C. elegans, but fitness improved under continuous-mating, consistent with sperm-limitation under gonochoristic (male/female) reproduction. C. inopinata also revealed greater fecundity and viability at higher temperatures. CONCLUSION: Consistent with observations in other ectotherms, slower growth in C. inopinata indicates a potential trade-off between body size and developmental timing, whereas its unchanged lifespan suggests that longevity is largely uncoupled from its increase in body size. Additionally, temperature-dependent patterns of fitness in C. inopinata are consistent with its geographic origins in subtropical Okinawa. Overall, these results underscore the extent to which changes in ecological context and body size can shape life history traits.

Proteomic and evolutionary analyses of sperm activation identify uncharacterized genes in Caenorhabditis nematodes.

BACKGROUND: Nematode sperm have unique and highly diverged morphology and molecular biology. In particular, nematode sperm contain subcellular vesicles known as membranous organelles that are necessary for male fertility, yet play a still unknown role in overall sperm function. Here we take a novel proteomic approach to characterize the functional protein complement of membranous organelles in two Caenorhabditis species: C. elegans and C. remanei. RESULTS: We identify distinct protein compositions between membranous organelles and the activated sperm body. Two particularly interesting and undescribed gene families-the Nematode-Specific Peptide family, group D and the here designated Nematode-Specific Peptide family, group F-localize to the membranous organelle. Both multigene families are nematode-specific and exhibit patterns of conserved evolution specific to the Caenorhabditis clade. These data suggest gene family dynamics may be a more prevalent mode of evolution than sequence divergence within sperm. Using a CRISPR-based knock-out of the NSPF gene family, we find no evidence of a male fertility effect of these genes, despite their high protein abundance within the membranous organelles. CONCLUSIONS: Our study identifies key components of this unique subcellular sperm component and establishes a path toward revealing their underlying role in reproduction.

Comparative genomic analysis of upstream miRNA regulatory motifs in Caenorhabditis.

MicroRNAs (miRNAs) comprise a class of short noncoding RNA molecules that play diverse developmental and physiological roles by controlling mRNA abundance and protein output of the vast majority of transcripts. Despite the importance of miRNAs in regulating gene function, we still lack a complete understanding of how miRNAs themselves are transcriptionally regulated. To fill this gap, we predicted regulatory sequences by searching for abundant short motifs located upstream of miRNAs in eight species of Caenorhabditis nematodes. We identified three conserved motifs across the Caenorhabditis phylogeny that show clear signatures of purifying selection from comparative genomics, patterns of nucleotide changes in motifs of orthologous miRNAs, and correlation between motif incidence and miRNA expression. We then validated our predictions with transgenic green fluorescent protein reporters and site-directed mutagenesis for a subset of motifs located in an enhancer region upstream of let-7 We demonstrate that a CT-dinucleotide motif is sufficient for proper expression of GFP in the seam cells of adult C. elegans, and that two other motifs play incremental roles in combination with the CT-rich motif. Thus, functional tests of sequence motifs identified through analysis of molecular evolutionary signatures provide a powerful path for efficiently characterizing the transcriptional regulation of miRNA genes.

Field studies reveal a close relative of C. elegans thrives in the fresh figs of Ficus septica and disperses on its Ceratosolen pollinating wasps.

BACKGROUND: Biotic interactions are ubiquitous and require information from ecology, evolutionary biology, and functional genetics in order to be understood. However, study systems that are amenable to investigations across such disparate fields are rare. Figs and fig wasps are a classic system for ecology and evolutionary biology with poor functional genetics; Caenorhabditis elegans is a classic system for functional genetics with poor ecology. In order to help bridge these disciplines, here we describe the natural history of a close relative of C. elegans, Caenorhabditis inopinata, that is associated with the fig Ficus septica and its pollinating Ceratosolen wasps. RESULTS: To understand the natural context of fig-associated Caenorhabditis, fresh F. septica figs from four Okinawan islands were sampled, dissected, and observed under microscopy. C. inopinata was found in all islands where F. septica figs were found. C.i nopinata was routinely found in the fig interior and almost never observed on the outside surface. C. inopinata was only found in pollinated figs, and C. inopinata was more likely to be observed in figs with more foundress pollinating wasps. Actively reproducing C. inopinata dominated early phase figs, whereas late phase figs with emerging wasp progeny harbored C. inopinata dauer larvae. Additionally, C. inopinata was observed dismounting from Ceratosolen pollinating wasps that were placed on agar plates. C. inopinata was not found on non-pollinating, parasitic Philotrypesis wasps. Finally, C. inopinata was only observed in F. septica figs among five Okinawan Ficus species sampled. CONCLUSION: These are the first detailed field observations of C. inopinata, and they suggest a natural history where this species proliferates in early phase F. septica figs and disperses from late phase figs on Ceratosolen pollinating fig wasps. While consistent with other examples of nematode diversification in the fig microcosm, the fig and wasp host specificity of C. inopinata is highly divergent from the life histories of its close relatives and frames hypotheses for future investigations. This natural co-occurrence of the fig/fig wasp and C. inopinata study systems sets the stage for an integrated research program that can help to explain the evolution of interspecific interactions.

Reproductive Mode and the Evolution of Genome Size and Structure in Caenorhabditis Nematodes.

The self-fertile nematode worms Caenorhabditis elegans, C. briggsae, and C. tropicalis evolved independently from outcrossing male-female ancestors and have genomes 20-40% smaller than closely related outcrossing relatives. This pattern of smaller genomes for selfing species and larger genomes for closely related outcrossing species is also seen in plants. We use comparative genomics, including the first high quality genome assembly for an outcrossing member of the genus (C. remanei) to test several hypotheses for the evolution of genome reduction under a change in mating system. Unlike plants, it does not appear that reductions in the number of repetitive elements, such as transposable elements, are an important contributor to the change in genome size. Instead, all functional genomic categories are lost in approximately equal proportions. Theory predicts that self-fertilization should equalize the effective population size, as well as the resulting effects of genetic drift, between the X chromosome and autosomes. Contrary to this, we find that the self-fertile C. briggsae and C. elegans have larger intergenic spaces and larger protein-coding genes on the X chromosome when compared to autosomes, while C. remanei actually has smaller introns on the X chromosome than either self-reproducing species. Rather than being driven by mutational biases and/or genetic drift caused by a reduction in effective population size under self reproduction, changes in genome size in this group of nematodes appear to be caused by genome-wide patterns of gene loss, most likely generated by genomic adaptation to self reproduction per se.

Genomic Signatures of Sexual Conflict.

Sexual conflict is a specific class of intergenomic conflict that describes the reciprocal sex-specific fitness costs generated by antagonistic reproductive interactions. The potential for sexual conflict is an inherent property of having a shared genome between the sexes and, therefore, is an extreme form of an environment-dependent fitness effect. In this way, many of the predictions from environment-dependent selection can be used to formulate expected patterns of genome evolution under sexual conflict. However, the pleiotropic and transmission constraints inherent to having alleles move across sex-specific backgrounds from generation to generation further modulate the anticipated signatures of selection. We outline methods for detecting candidate sexual conflict loci both across and within populations. Additionally, we consider the ability of genome scans to identify sexually antagonistic loci by modeling allele frequency changes within males and females due to a single generation of selection. In particular, we highlight the need to integrate genotype, phenotype, and functional information to truly distinguish sexual conflict from other forms of sexual differentiation.

Intrinsic differences between males and females determine sex-specific consequences of inbreeding.

BACKGROUND: Inbreeding increases homozygosity and exposes deleterious recessive alleles, generally decreasing the fitness of inbred individuals. Interestingly, males and females are usually affected differently by inbreeding, though the more vulnerable sex depends on the species and trait measured. RESULTS: We used the soil-dwelling nematode Caenorhabditis remanei to examine sex-specific inbreeding depression across nine lineages, five levels of inbreeding, and hundreds of thousands of progeny. Female nematodes consistently suffered greater fitness losses than their male counterparts, especially at high levels of inbreeding. CONCLUSIONS: These results suggest that females experience stronger selection on genes contributing to reproductive traits. Inbreeding depression in males may be further reduced by sex chromosome hemizygosity, which affects the dominance of some mutations, as well as by the absence of sexual selection. Determining the relative contributions of sex-specific expression, genes on the sex chromosomes, and the environment they are filtered through-including opportunities for sexual selection-may explain the frequent though inconsistent records of sex differences in inbreeding depression, along with their implications for conservation and the evolution of mating systems.

High-specificity detection of rare alleles with Paired-End Low Error Sequencing (PELE-Seq).

BACKGROUND: Polymorphic loci exist throughout the genomes of a population and provide the raw genetic material needed for a species to adapt to changes in the environment. The minor allele frequencies of rare Single Nucleotide Polymorphisms (SNPs) within a population have been difficult to track with Next-Generation Sequencing (NGS), due to the high error rate of standard methods such as Illumina sequencing. RESULTS: We have developed a wet-lab protocol and variant-calling method that identifies both sequencing and PCR errors, called Paired-End Low Error Sequencing (PELE-Seq). To test the specificity and sensitivity of the PELE-Seq method, we sequenced control E. coli DNA libraries containing known rare alleles present at frequencies ranging from 0.2-0.4 % of the total reads. PELE-Seq had higher specificity and sensitivity than standard libraries. We then used PELE-Seq to characterize rare alleles in a Caenorhabditis remanei nematode worm population before and after laboratory adaptation, and found that minor and rare alleles can undergo large changes in frequency during lab-adaptation. CONCLUSION: We have developed a method of rare allele detection that mitigates both sequencing and PCR errors, called PELE-Seq. PELE-Seq was evaluated using control E. coli populations and was then used to compare a wild C. remanei population to a lab-adapted population. The PELE-Seq method is ideal for investigating the dynamics of rare alleles in a broad range of reduced-representation sequencing methods, including targeted amplicon sequencing, RAD-Seq, ddRAD, and GBS. PELE-Seq is also well-suited for whole genome sequencing of mitochondria and viruses, and for high-throughput rare mutation screens.

Natural and experimental evolution of sexual conflict within Caenorhabditis nematodes.

BACKGROUND: Although males and females need one another in order to reproduce, they often have different reproductive interests, which can lead to conflict between the sexes. The intensity and frequency of male-male competition for fertilization opportunities is thought to be an important contributor to this conflict. The nematode genus Caenorhabditis provides an opportunity to test this hypothesis because the frequency of males varies widely among species with different mating systems. RESULTS: We find evidence that there is strong inter- and intra-sexual conflict within C. remanei, a dioecious species composed of equal frequencies of males and females. In particular, some C. remanei males greatly reduce female lifespan following mating, and their sperm have a strong competitive advantage over the sperm of other males. In contrast, our results suggest that both types of conflict have been greatly reduced within C. elegans, which is an androdioecious species that is composed of self-fertilizing hermaphrodites and rare males. Using experimental evolution in mutant C. elegans populations in which sperm production is blocked in hermaphrodites (effectively converting them to females), we find that the consequences of sexual conflict observed within C. remanei evolve rapidly within C. elegans populations experiencing high levels of male-male competition. CONCLUSIONS: Together, these complementary data sets support the hypothesis that the intensity of intersexual conflict varies with the intensity of competition among males, and that male-induced collateral damage to mates can evolve very rapidly within populations.

Modeling the evolution of complex genetic systems: the gene network family tree.

In 1994 and 1996, Andreas Wagner introduced a novel model in two papers addressing the evolution of genetic regulatory networks. This work, and a suite of papers that followed using similar models, helped integrate network thinking into biology and motivate research focused on the evolution of genetic networks. The Wagner network has its mathematical roots in the Ising model, a statistical physics model describing the activity of atoms on a lattice, and in neural networks. These models have given rise to two branches of applications, one in physics and biology and one in artificial intelligence and machine learning. Here, we review development along these branches, outline similarities and differences between biological models of genetic regulatory circuits and neural circuits models used in machine learning, and identify ways in which these models can provide novel insights into biological systems.

Mutation load and rapid adaptation favour outcrossing over self-fertilization.

The tendency of organisms to reproduce by cross-fertilization despite numerous disadvantages relative to self-fertilization is one of the oldest puzzles in evolutionary biology. For many species, the primary obstacle to the evolution of outcrossing is the cost of production of males, individuals that do not directly contribute offspring and thus diminish the long-term reproductive output of a lineage. Self-fertilizing ('selfing') organisms do not incur the cost of males and therefore should possess at least a twofold numerical advantage over most outcrossing organisms. Two competing explanations for the widespread prevalence of outcrossing in nature despite this inherent disadvantage are the avoidance of inbreeding depression generated by selfing and the ability of outcrossing populations to adapt more rapidly to environmental change. Here we show that outcrossing is favoured in populations of Caenorhabditis elegans subject to experimental evolution both under conditions of increased mutation rate and during adaptation to a novel environment. In general, fitness increased with increasing rates of outcrossing. Thus, each of the standard explanations for the maintenance of outcrossing are correct, and it is likely that outcrossing is the predominant mode of reproduction in most species because it is favoured under ecological conditions that are ubiquitous in natural environments.

The transgenerational effects of heat stress in the nematode Caenorhabditis remanei are negative and rapidly eliminated under direct selection for increased stress resistance in larvae.

Parents encountering stress environments can influence the phenotype of their offspring in a form of transgenerational phenotypic plasticity that has the potential to be adaptive if offspring are thereby better able to deal with future stressors. Here, we test for the existence of anticipatory parental effects in the heat stress response in the highly polymorphic nematode Caenorhabditis remanei. Rather providing an anticipatory response, parents subject to a prior heat stress actually produce offspring that are less able to survive a severe heat shock. Selection on heat shock resistance within the larvae via experimental evolution leads to a loss of sensitivity (robustness) to environmental variation during both the parental and larval periods. Whole genome transcriptional analysis of both ancestor and selected lines shows that there is weak correspondence between genetic pathways induced via temperature shifts during parental and larval periods. Parental effects can evolve very rapidly via selection acting directly on offspring.

Hermaphrodite life history and the maintenance of partial selfing in experimental populations of Caenorhabditis elegans.

BACKGROUND: Classic population genetics theory predicts that mixed reproductive systems, where self reproduction (selfing) and outcrossing co-exist, should not be as common as they are in nature. One means of reconciling theory with observations is to recognize that sexual conflict between males and hermaphrodites and/or constraints in the allocation of resources towards sex functions in hermaphrodites can balance the fitness components of selfing and outcrossing. RESULTS: Using experimental evolution in Caenorhabditis elegans, we test whether the adaptive maintenance of partial selfing is due to sexual conflict and/or to the evolution of sex allocation towards male function in hermaphrodites. For this, we characterized the reproductive schedule and longevity patterns in hermaphrodites under selfing and under outcrossing with naïve males that did not have the opportunity to evolve with them. A shift in reproductive schedule towards earlier reproduction would be indicative of adaptation in our imposed life-cycle, while longevity is expected to evolve as a response to the harm that males impinge on hermaphrodites upon mating. To determine adaptation in the absence of constraints in sex allocation, we also characterized the life history of females that reproduced during experimental evolution through obligate mating with males. As expected with adaptation, we find that after 100 generations of experimental evolution, selfing hermaphrodites and females showed improved reproduction at earlier ages. We did not observe similar reproductive shifts in outcrossed hermaphrodites. We further find increased longevity in outcrossed females after evolution but not in outcrossed hermaphrodites, a result that indicates that sexual conflicts were likely more prevalent under male-female evolution than under male-hermaphrodite evolution. CONCLUSIONS: Taken together, our findings suggest that the adaptive maintenance of partial selfing during C. elegans experimental evolution resulted from the evolution of sex allocation towards male function in hermaphrodites.

Epistasis--the essential role of gene interactions in the structure and evolution of genetic systems.

Epistasis, or interactions between genes, has long been recognized as fundamentally important to understanding the structure and function of genetic pathways and the evolutionary dynamics of complex genetic systems. With the advent of high-throughput functional genomics and the emergence of systems approaches to biology, as well as a new-found ability to pursue the genetic basis of evolution down to specific molecular changes, there is a renewed appreciation both for the importance of studying gene interactions and for addressing these questions in a unified, quantitative manner.

Heritable epigenetic variation facilitates long-term maintenance of epigenetic and genetic variation.

How genetic and phenotypic variation are maintained has long been one of the fundamental questions in population and quantitative genetics. A variety of factors have been implicated to explain the maintenance of genetic variation in some contexts (e.g. balancing selection), but the potential role of epigenetic regulation to influence population dynamics has been understudied. It is well recognized that epigenetic regulation, including histone methylation, small RNA expression, and DNA methylation, helps to define differences between cell types and facilitate phenotypic plasticity. In recent years, empirical studies have shown the potential for epigenetic regulation to also be heritable for at least a few generations without selection, raising the possibility that differences in epigenetic regulation can act alongside genetic variation to shape evolutionary trajectories. Heritable differences in epigenetic regulation that arise spontaneously are termed "epimutations." Epimutations differ from genetic mutations in 2 key ways-they occur at a higher rate and the loci at which they occur often revert back to their original state within a few generations. Here, we present an extension of the standard population genetic model with selection to incorporate epigenetic variation arising via epimutation. Our model assumes a diploid, sexually reproducing population with random mating. In addition to spontaneous genetic mutation, we included parameters for spontaneous epimutation and back-epimutation, allowing for 4 potential epialleles at a single locus (2 genetic alleles, each with 2 epigenetic states), each of which affect fitness. We then analyzed the conditions under which stable epialleles were maintained. Our results show that highly reversible epialleles can be maintained in long-term equilibrium under neutral conditions in a manner that depends on the epimutation and back-epimutation rates, which we term epimutation-back-epimutation equilibrium. On the other hand, epialleles that compensate for deleterious mutations cause deviations from the expectations of mutation-selection balance by a simple factor that depends on the epimutation and back-epimutation rates. We also numerically analyze several sets of fitness parameters for which large deviations from mutation-selection balance occur. Together, these results demonstrate that transient epigenetic regulation may be an important factor in the maintenance of both epigenetic and genetic variation in populations.

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans.

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

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.