The other side of phenotypic plasticity: a developmental system that generates an invariant phenotype despite environmental variation

Understanding how the environment impacts development is of central interest in developmental and evolutionary biology. On the one hand, we would like to understand how the environment induces phenotypic changes (the study of phenotypic plasticity). On the other hand, we may ask how a development system maintains a stable and precise phenotypic output despite the presence of environmental variation. We study such developmental robustness to environmental variation using vulval cell fate patterning in the nematode Caenorhabditis elegans as a study system. Here we review both mechanistic and evolutionary aspects of these studies, focusing on recently obtained experimental results. First, we present evidence indicating that vulval formation is under stabilizing selection. Second, we discusss quantitative data on the precision and variability in the output of the vulval developmental system in different environments and different genetic backgrounds. Third, we illustrate how environmental and genetic variation modulate the cellular and molecular processes underlying the formation of the vulva. Fourth, we discuss the evolutionary significance of environmental sensitivity of this developmental system.     

Evolution of a system sensitive to stochastic noise: P3.p cell fate in Caenorhabditis

The C. elegans cell lineage is overall invariant. One rare instance of variability concerns P3.p, the most anterior vulva precursor cell, which may either fuse with the epidermis without dividing, or remain competent to form vulval tissue and divide. Here we examine the evolutionary properties of this stochastic variation in P3.p fate. In the Caenorhabditis genus, high P3.p competence is ancestral and reduction in P3.p competence and division frequency occurred in C. sp. 14 and in a clade of nine species. Within this clade, the frequency of P3.p division further varies within and among species, being intermediate in C. elegans and low in C. briggsae. P3.p fate frequency is sensitive to random mutation accumulation, suggesting that this trait may evolve rapidly because of its sensitivity to mutational impact. P3.p fate depends on LIN-39/Hox5 expression and we find that the peak of LIN-39/Hox5 protein level is displaced posteriorly in C. briggsae compared to C. elegans. However, P3.p fate specification is most sensitive to the dose of EGL-20 and CWN-1, two Wnts that are secreted in a long-range gradient from the posterior end of C. elegans larvae (accompanying article). A half-dose of either of these Wnts is sufficient to affect division frequency in C. elegans N2 to levels similar to those in C. briggsae. Symmetrically, we show that an increase in Wnt dose rescues anterior competence in C. briggsae. We propose that evolutionary variation in the concentration or interpretation of the long-range Wnt gradient may be involved in the rapid evolution of P3.p fate in Caenorhabditis.     

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.     

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

Comparative genomic analysis of important signaling pathways in Caenorhabditis briggsae and Caenorhabditis elegans reveals both conserved features and also differences. To build a framework to address the significance of these features we determined the C. briggsae embryonic cell lineage, using the tools StarryNite and AceTree. We traced both cell divisions and cell positions for all cells through all but the last round of cell division and for selected cells through the final round. We found the lineage to be remarkably similar to that of C. elegans. Not only did the founder cells give rise to similar numbers of progeny, the relative cell division timing and positions were largely maintained. These lineage similarities appear to give rise to similar cell fates as judged both by the positions of lineally equivalent cells and by the patterns of cell deaths in both species. However, some reproducible differences were seen, e.g., the P4 cell cycle length is more than 40% longer in C. briggsae than that in C. elegans (p < 0.01). The extensive conservation of embryonic development between such divergent species suggests that substantial evolutionary distance between these two species has not altered these early developmental cellular events, although the developmental defects of transpecies hybrids suggest that the details of the underlying molecular pathways have diverged sufficiently so as to not be interchangeable.     

Mutational bias for body size in rhabditid nematodes

Mutational bias is a potentially important agent of evolution, but it is difficult to disentangle the effects of mutation from those of natural selection. Mutation-accumulation experiments, in which mutations are allowed to accumulate at very small population size, thus minimizing the efficiency of natural selection, are the best way to separate the effects of mutation from those of selection. Body size varies greatly among species of nematode in the family rhabditidae; mutational biases are both a potential cause and a consequence of that variation. We report data on the cumulative effects of mutations that affect body size in three species of rhabditid nematode that vary fivefold in adult size. Results are very consistent with previous studies of mutations underlying fitness in the same strains: two strains of Caenorhabditis briggsae decline in body size about twice as fast as two strains of C. elegans, with a concomitant higher point estimate of the genomic mutation rate; the confamilial Oscheius myriophila is intermediate. There is an overall mutational bias, such that mutations reduce size on average, but the bias appears consistent between species. The genetic correlation between mutations that affect size and those underlying fitness is large and positive, on average.     

Predicting Phenotypic Diversity and the Underlying Quantitative Molecular Transitions

During development, signaling networks control the formation of multicellular patterns. To what extent quantitative fluctuations in these complex networks may affect multicellular phenotype remains unclear. Here, we describe a computational approach to predict and analyze the phenotypic diversity that is accessible to a developmental signaling network. Applying this framework to vulval development in C. elegans, we demonstrate that quantitative changes in the regulatory network can render ~500 multicellular phenotypes. This phenotypic capacity is an order-of-magnitude below the theoretical upper limit for this system but yet is large enough to demonstrate that the system is not restricted to a select few outcomes. Using metrics to gauge the robustness of these phenotypes to parameter perturbations, we identify a select subset of novel phenotypes that are the most promising for experimental validation. In addition, our model calculations provide a layout of these phenotypes in network parameter space. Analyzing this landscape of multicellular phenotypes yielded two significant insights. First, we show that experimentally well-established mutant phenotypes may be rendered using non-canonical network perturbations. Second, we show that the predicted multicellular patterns include not only those observed in C. elegans, but also those occurring exclusively in other species of the Caenorhabditis genus. This result demonstrates that quantitative diversification of a common regulatory network is indeed demonstrably sufficient to generate the phenotypic differences observed across three major species within the Caenorhabditis genus. Using our computational framework, we systematically identify the quantitative changes that may have occurred in the regulatory network during the evolution of these species. Our model predictions show that significant phenotypic diversity may be sampled through quantitative variations in the regulatory network without overhauling the core network architecture. Furthermore, by comparing the predicted landscape of phenotypes to multicellular patterns that have been experimentally observed across multiple species, we systematically trace the quantitative regulatory changes that may have occurred during the evolution of the Caenorhabditis genus.     

Why the short face? Developmental disintegration of the neurocranium drives convergent evolution in neotropical electric fishes

Convergent evolution is widely viewed as strong evidence for the influence of natural selection on the origin of phenotypic design. However, the emerging evo‐devo synthesis has highlighted other processes that may bias and direct phenotypic evolution in the presence of environmental and genetic variation. Developmental biases on the production of phenotypic variation may channel the evolution of convergent forms by limiting the range of phenotypes produced during ontogeny. Here, we study the evolution and convergence of brachycephalic and dolichocephalic skull shapes among 133 species of Neotropical electric fishes (Gymnotiformes: Teleostei) and identify potential developmental biases on phenotypic evolution. We plot the ontogenetic trajectories of neurocranial phenotypes in 17 species and document developmental modularity between the face and braincase regions of the skull. We recover a significant relationship between developmental covariation and relative skull length and a significant relationship between developmental covariation and ontogenetic disparity. We demonstrate that modularity and integration bias the production of phenotypes along the brachycephalic and dolichocephalic skull axis and contribute to multiple, independent evolutionary transformations to highly brachycephalic and dolichocephalic skull morphologies.     

Molecular population genetics and phenotypic sensitivity to ethanol for a globally diverse sample of the nematode Caenorhabditis briggsae

New genomic resources and genetic tools of the past few years have advanced the nematode genus Caenorhabditis as a model for comparative biology. However, understanding of natural genetic variation at molecular and phenotypic levels remains rudimentary for most species in this genus, and for C. briggsae in particular. Here we characterize phenotypic variation in C. briggsae’s sensitivity to the potentially important and variable environmental toxin, ethanol, for globally diverse strains. We also quantify nucleotide variation in a new sample of 32 strains from four continents, including small islands, and for the closest‐known relative of this species (C. sp. 9). We demonstrate that C. briggsae exhibits little heritable variation for the effects of ethanol on the norm of reaction for survival and reproduction. Moreover, C. briggsae does not differ significantly from C. elegans in our assays of its response to this substance that both species likely encounter regularly in habitats of rotting fruit and vegetation. However, we uncover drastically more molecular genetic variation than was known previously for this species, despite most strains, including all island strains, conforming to the broad biogeographic patterns described previously. Using patterns of sequence divergence between populations and between species, we estimate that the self‐fertilizing mode of reproduction by hermaphrodites in C. briggsae likely evolved sometime between 0.9 and 10 million generations ago. These insights into C. briggsae’s natural history and natural genetic variation greatly expand the potential of this organism as an emerging model for studies in molecular and quantitative genetics, the evolution of development, and ecological genetics.     

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.     

The Genome Sequence of Caenorhabditis briggsae: A Platform for Comparative Genomics

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes.     

Autosomal genes of autosomal/X-linked duplicated gene pairs and germ-line proliferation in Caenorhabditis elegans

We report molecular genetic studies of three genes involved in early germ-line proliferation in Caenorhabditis elegans that lend unexpected insight into a germ-line/soma functional separation of autosomal/X-linked duplicated gene pairs. In a genetic screen for germ-line proliferation-defective mutants, we identified mutations in rpl-11.1 (L11 protein of the large ribosomal subunit), pab-1 [a poly(A)-binding protein], and glp-3/eft-3 (an elongation factor 1-α homolog). All three are members of autosome/X gene pairs. Consistent with a germ-line-restricted function of rpl-11.1 and pab-1, mutations in these genes extend life span and cause gigantism. We further examined the RNAi phenotypes of the three sets of rpl genes (rpl-11, rpl-24, and rpl-25) and found that for the two rpl genes with autosomal/X-linked pairs (rpl-11 and rpl-25), zygotic germ-line function is carried by the autosomal copy. Available RNAi results for highly conserved autosomal/X-linked gene pairs suggest that other duplicated genes may follow a similar trend. The three rpl and the pab-1/2 duplications predate the divergence between C. elegans and C. briggsae, while the eft-3/4 duplication appears to have occurred in the lineage to C. elegans after it diverged from C. briggsae. The duplicated C. briggsae orthologs of the three C. elegans autosomal/X-linked gene pairs also display functional differences between paralogs. We present hypotheses for evolutionary mechanisms that may underlie germ-line/soma subfunctionalization of duplicated genes, taking into account the role of X chromosome silencing in the germ line and analogous mammalian phenomena.     

Species differentiation in Caenorhabditis briggsae and Caenorhabditis elegans

Identification of five laboratory strains (1-5) of putative Caenorhabditis briggsae was undertaken. Examination of the male bursal ray arrangement, mating tests with males of Caenorhabditis elegans, malate dehydrogenase zymograms, and SDS polyacrylamide electrophoresis demonstrated that strain 4 was C. briggsae and the others were C. elegans.     

Microevolutionary studies in nematodes: a beginning

Comparisons between related species often allow the detailed genetic analysis of evolutionary processes. Here we advocate the use of the nematode Caenorhabditis elegans (and several other rhabditid species) as model systems for microevolutionary studies. Compared to Drosophila species, which have been a mainstay of such studies, C. elegans has a self-fertilizing mode of reproduction, a shorter life cycle and a convenient cell-level analysis of phenotypic variation. Data concerning its population genetics and ecology are still scarce, however. We review molecular, behavioral and developmental intraspecific polymorphisms for populations of C. elegans, Oscheius sp. 1 and Pristionchus pacificus. Focusing on vulval development, which has been well characterized in several species, we discuss relationships between patterns of variations: (1) for a given genotype (developmental variants), (2) after mutagenesis (mutability), (3) in different populations of the same species (polymorphisms) and (4) between closely related species. These studies have revealed that evolutionary variations between sister species affect those characters that show phenotypic developmental variants, that are mutable and that are polymorphic within species.     

Nematode selenoproteome: the use of the selenocysteine insertion system to decode one codon in an animal genome?

Selenocysteine (Sec) is co-translationally inserted into selenoproteins in response to codon UGA with the help of the selenocysteine insertion sequence (SECIS) element. The number of selenoproteins in animals varies, with humans having 25 and mice having 24 selenoproteins. To date, however, only one selenoprotein, thioredoxin reductase, has been detected in Caenorhabditis elegans, and this enzyme contains only one Sec. Here, we characterize the selenoproteomes of C.elegans and Caenorhabditis briggsae with three independent algorithms, one searching for pairs of homologous nematode SECIS elements, another searching for Cys- or Sec-containing homologs of potential nematode selenoprotein genes and the third identifying Sec-containing homologs of annotated nematode proteins. These methods suggest that thioredoxin reductase is the only Sec-containing protein in the C.elegans and C.briggsae genomes. In contrast, we identified additional selenoproteins in other nematodes. Assuming that Sec insertion mechanisms are conserved between nematodes and other eukaryotes, the data suggest that nematode selenoproteomes were reduced during evolution, and that in an extreme reduction case Sec insertion systems probably decode only a single UGA codon in C.elegans and C.briggsae genomes. In addition, all detected genes had a rare form of SECIS element containing a guanosine in place of a conserved adenosine present in most other SECIS structures, suggesting that in organisms with small selenoproteomes SECIS elements may change rapidly.     

The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes.     

Differences in the selection response of serially repeated color pattern characters: Standing variation,development, and evolution

Background  

There is spectacular morphological diversity in nature but lineages typically display a limited range of phenotypes. Because developmental processes generate the phenotypic variation that fuels natural selection, they are a likely source of evolutionary biases, facilitating some changes and limiting others. Although shifts in developmental regulation are associated with morphological differences between taxa, it is unclear how underlying mechanisms affect the rate and direction of evolutionary change within populations under selection.     

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.