Hakuna Nematoda: genetic and phenotypic diversity in African isolates of Caenorhabditis elegans and C. briggsae

Caenorhabditis elegans and C. briggsae have many parallels in terms of morphology, life history and breeding system. Both species also share similar low levels of molecular diversity, although the global sampling of natural populations has been limited and geographically biased. In this study, we describe the first cultured isolates of C. elegans and C. briggsae from sub-Saharan Africa. We characterize these samples for patterns of nucleotide polymorphism and vulva precursor cell lineage, and conduct a series of hybrid crosses in C. briggsae to test for genetic incompatibilities. The distribution of genetic diversity confirms a lack of geographic structure to C. elegans sequences but shows genetic differentiation of C. briggsae into three distinct clades that may correspond to three latitudinal ranges. Despite low levels of molecular diversity, we find considerable variation in cell division frequency in African C. elegans for the P3.p vulva precursor cell, and in African C. briggsae for P4.p, a variation that was not previously observed in this species. Hybrid crosses did not reveal major incompatibilities between C. briggsae strains from Africa and elsewhere, and there was some evidence of inbreeding depression. These new African isolates suggest that important ecological factors may be shaping the patterns of diversity in C. briggsae, and that despite many similarities between C. elegans and C. briggsae, there may be more subtle differences in their natural histories than previously appreciated.     

Caenorhabditis briggsae recombinant inbred line genotypes reveal inter-strain incompatibility and the evolution of recombination

The nematode Caenorhabditis briggsae is an emerging model organism that allows evolutionary comparisons with C. elegans and exploration of its own unique biological attributes. To produce a high-resolution C. briggsae recombination map, recombinant inbred lines were generated from reciprocal crosses between two strains and genotyped at over 1,000 loci. A second set of recombinant inbred lines involving a third strain was also genotyped at lower resolution. The resulting recombination maps exhibit discrete domains of high and low recombination, as in C. elegans, indicating these are a general feature of Caenorhabditis species. The proportion of a chromosome's physical size occupied by the central, low-recombination domain is highly correlated between species. However, the C. briggsae intra-species comparison reveals striking variation in the distribution of recombination between domains. Hybrid lines made with the more divergent pair of strains also exhibit pervasive marker transmission ratio distortion, evidence of selection acting on hybrid genotypes. The strongest effect, on chromosome III, is explained by a developmental delay phenotype exhibited by some hybrid F2 animals. In addition, on chromosomes IV and V, cross direction-specific biases towards one parental genotype suggest the existence of cytonuclear epistatic interactions. These interactions are discussed in relation to surprising mitochondrial genome polymorphism in C. briggsae, evidence that the two strains diverged in allopatry, the potential for local adaptation, and the evolution of Dobzhansky-Muller incompatibilities. The genetic and genomic resources resulting from this work will support future efforts to understand inter-strain divergence as well as facilitate studies of gene function, natural variation, and the evolution of recombination in Caenorhabditis nematodes.     

Population genetics of Caenorhabditis elegans: the paradox of low polymorphism in a widespread species

Caenorhabditis elegans has become one of the most widely used model research organisms, yet we have little information on evolutionary processes and recent evolutionary history of this widespread species. We examined patterns of variation at 20 microsatellite loci in a sample of 23 natural isolates of C. elegans from various parts of the world. One-half of the loci were monomorphic among all strains, and overall genetic variation at microsatellite loci was low, relative to most other species. Some population structure was detected, but there was no association between the genetic and geographic distances among different natural isolates. Thus, despite the nearly worldwide occurrence of C. elegans, little evidence was found for local adaptation in strains derived from different parts of the world. The low levels of genetic variation within and among populations suggest that recent colonization and population expansion might have occurred. However, the patterns of variation are not consistent with population expansion. A possible explanation for the observed patterns is the action of background selection to reduce polymorphism, coupled with ongoing gene flow among populations worldwide.     

Molecular evolution and quantitative variation for chemosensory behaviour in the nematode genus Caenorhabditis

Caenorhabditis elegans is a model organism in biology, yet despite the tremendous information generated from genetic, genomic and functional analyses, C. elegans has rarely been used to address questions in ecological genetics. Here, we analyse genetic variation for chemosensory behaviour, an ecologically important trait that is also genetically well characterized, at both the phenotypic and molecular levels within three species of the genus Caenorhabditis. We show that the G-protein ODR-3 plays an important role in chemosensory avoidance behaviour and identify orthologues of odr-3 in C. briggsae and C. remanei. Both quantitative genetic analysis of chemosensory behaviour and molecular population genetic analysis of odr-3 show that there is little genetic variation among a worldwide collection of isolates of the primarily selfing C. elegans, whereas there is substantially more variation within a single population of the outcrossing C. remanei. Although there are a large number of substitutions at silent sites within odr-3 among the three species, molecular evolution at the protein level is extremely conserved, suggesting that odr-3 plays an important role in cell signalling during chemosensation and/or neuronal cilia development in C. remanei and in C. briggsae as it does in C. elegans. Our results suggest that C. remanei may be a more suitable subject for ecological and evolutionary genetic studies than C. elegans.     

High nucleotide polymorphism and rapid decay of linkage disequilibrium in wild populations of Caenorhabditis remanei

The common ancestor of the self-fertilizing nematodes Caenorhabditis elegans and C. briggsae must have reproduced by obligate outcrossing, like most species in this genus. However, we have only a limited understanding about how genetic variation is patterned in such male-female (gonochoristic) Caenorhabditis species. Here, we report results from surveying nucleotide variation of six nuclear loci in a broad geographic sample of wild isolates of the gonochoristic C. remanei. We find high levels of diversity in this species, with silent-site diversity averaging 4.7%, implying an effective population size close to 1 million. Additionally, the pattern of polymorphisms reveals little evidence for population structure or deviation from neutral expectations, suggesting that the sampled C. remanei populations approximate panmixis and demographic equilibrium. Combined with the observation that linkage disequilibrium between pairs of polymorphic sites decays rapidly with distance, this suggests that C. remanei will provide an excellent system for identifying the genetic targets of natural selection from deviant patterns of polymorphism and linkage disequilibrium. The patterns revealed in this obligately outcrossing species may provide a useful model of the evolutionary circumstances in C. elegans' gonochoristic progenitor. This will be especially important if self-fertilization evolved recently in C. elegans history, because most of the evolutionary time separating C. elegans from its known relatives would have occurred in a state of obligate outcrossing.     

A comparative genomic analysis of the small heat shock proteins in Caenorhabditis elegans and briggsae

The small heat shock proteins (sHSPs) are a ubiquitous family of molecular chaperones. We have identified 18 sHSPs in the Caenorhabditis elegans genome and 20 sHSPs in the Caenorhabditis briggsae genome. Analysis of phylogenetic relationships and evolutionary dynamics of the sHSPs in these two genomes reveals a very complex pattern of evolution. The sHSPs in C. elegans and C. briggsae do not display clear orthologous relationships with other invertebrate sHSPs. But many sHSPs in C. elegans have orthologs in C. briggsae. One group of sHSPs, the HSP16s, has a very unusual evolutionary history. Although there are a number of HSP16s in both the C. elegans and C. briggsae genomes, none of the HSP16s display orthologous relationships across these two species. The HSP16s have an unusual gene pair structure and a complex evolutionary history shaped by gene duplication, gene conversion, and purifying selection. We found no evidence of recent positive selection acting on any of the sHSPs in C. elegans or in C. briggsae. There is also no evidence of functional divergence within the pairs of orthologous C. elegans and C. briggsae sHSPs. However, the evolutionary patterns do suggest that functional divergence has occurred between the sHSPs in C. elegans and C. briggsae and the sHSPs in more distantly related invertebrates.     

Genetic analysis of dauer formation in Caenorhabditis briggsae

Molecular changes that underlie evolutionary changes in behavior and physiology are not well understood. Dauer formation in Caenorhabditis elegans is a temperature-sensitive process controlled through a network of signaling pathways associated with sensory neurons and is potentially an excellent system in which to investigate molecular changes in neuronal function during evolution. To begin to investigate the evolution of dauer formation in the genus Caenorhabditis at the molecular level, we isolated dauer-formation mutations in C. briggsae, a species closely related to the model organism C. elegans. We identified mutations in orthologs of C. elegans genes daf-2 (insulin receptor), daf-3 (Smad), and daf-4 (TGF-beta type 2 receptor), as well as genes required for formation of sensory cilia. Phenotypic analyses revealed that functions of these genes are conserved between C. elegans and C. briggsae. Analysis of C. briggsae mutations also revealed a significant difference between the two species in their responses to high temperatures (>26 degrees). C. elegans is strongly induced to form dauers at temperatures above 26 degrees, near the upper limit for growth of C. elegans. In contrast, C. briggsae, which is capable of growth at higher temperatures than C. elegans, lacks this response.     

Molecular relationships between closely related strains and species of nematodes

Summary Electrophoretic comparisons have been made for 24 enzymes in theBergerac andBristol strains ofCaenorhabditis elegans and the related species,Caenorhabditis briggsae. No variation was detected between the two strains ofC. elegans. In contrast, the two species,C. elegans andC. briggsae exhibited electrophoretic differences in 22 of 24 enzymes. A consensus 5S rRNA sequence was determined forC. elegans and found to be identical to that fromC. briggsae. By analogy with other species with relatively well established fossil records it can be inferred that the time of divergence between the two nematode species is probably in the tens of millions of years.The limited anatomical evolution during a time period in which proteins undergo extensive changes supports the hypothesis that anatomical evolution is not dependent on overall protein changes.     

Diversity in mating behavior of hermaphroditic and male-female Caenorhabditis nematodes

In this study, we addressed why Caenorhabditis elegans males are inefficient at fertilizing their hermaphrodites. During copulation, hermaphrodites generally move away from males before they become impregnated. C. elegans hermaphrodites reproduce by internal self-fertilization, so that copulation with males is not required for species propagation. The hermaphroditic mode of reproduction could potentially relax selection for genes that optimize male mating behavior. We examined males from hermaphroditic and gonochoristic (male-female copulation) Caenorhabditis species to determine if they use different sensory and motor mechanisms to control their mating behavior. Instead, we found through laser ablation analysis and behavioral observations that hermaphroditic C. briggsae and gonochoristic C. remanei and Caenorhabditis species 4, PB2801 males produce a factor that immobilizes females during copulation. This factor also stimulates the vulval slit to widen, so that the male copulatory spicules can easily insert. C. elegans and C. briggsae hermaphrodites are not affected by this factor. We suggest that sensory and motor execution of mating behavior have not significantly changed among males of different Caenorhabditis species; however, during the evolution of internal self-fertilization, hermaphrodites have lost the ability to respond to the male soporific-inducing factor.     

Smoke without fire: most reported cases of intron gain in nematodes instead reflect intron losses

Identification of recently gained spliceosomal introns would provide crucial evidence in the continuing debate concerning the age and evolutionary significance of introns. A previously published genomic analysis reported to have identified 122 introns that had been gained since the divergence of the nematodes Caenorhabidits elegans and Caenorhabditis briggsae approximately 100 MYA. However, using newly available genomic sequence from additional Caenorhabditis species, we show that 74% (60/81) of the reported gains in C. elegans are present in a C. briggsae relative. This pattern indicates that these introns represent losses in C. briggsae, not gains in C. elegans. In addition, 61% (25/41) of the reported gains in C. briggsae are present in the more distant C. briggsae relative, in a pattern suggesting that additional reported gains in C. elegans and/or C. briggsae may in fact represent unrecognized losses. These results underscore the dominance of intron loss over intron gain in recent eukaryotic evolution, the pitfalls associated with parsimony in inferring intron gains, and the importance of genomic sequencing of clusters of closely related species for drawing accurate inferences about genome evolution.     

Oscheius tipulae, a widespread hermaphroditic soil nematode, displays a higher genetic diversity and geographical structure than Caenorhabditis elegans

The nematode Oscheius tipulae belongs to the same family (Rhabditidae) as the model species Caenorhabditis elegans . Both species reproduce through self-fertilizing hermaphrodites and facultative males. Recent studies have shown that the self-fertile C. elegans and C. briggsae displayed a 20-fold lower genetic diversity than the male–female species C. remanei . Several explanations have been put forward to account for this difference, including their mode of reproduction and dynamic population structure. Here, we present the results of extensive worldwide sampling of O. tipulae , which we previously used as a laboratory organism for developmental genetics. We found that O. tipulae is much more widespread and common in soil throughout the world than Caenorhabditis species. We analysed 63 O. tipulae isolates from several continents using amplified fragment length polymorphism (AFLP). We found that O. tipulae harbours a 5-fold higher genetic diversity than C. elegans and C. briggsae . As in C. elegans , a high proportion of this diversity was found locally. Yet, we detected significant geographical differentiation, both at the worldwide scale with a latitudinal structure and between three localities in France. In summary, O. tipulae exhibited significantly higher levels of genetic diversity and large-scale geographical structure than C. elegans , despite their shared mode of reproduction. This species difference in genetic diversity may be explained by a number of other differences, such as population size, distribution, migration and dynamics. Due to its widespread occurrence and relatively high genetic diversity, O. tipulae may be a promising study species for evolutionary studies.     

Conservation of gene organization and trans-splicing in the glyceraldehyde-3-phosphate dehydrogenase-encoding genes of Caenorhabditis briggsae.

The genes encoding body-wall-specific glyceraldehyde-3-phosphate dehydrogenase from Caenorhabditis briggsae were sequenced and compared to the homologous genes from Caenorhabditis elegans. The direct tandem organization of these genes, gpd-2 and gpd-3, and the size and location of the two introns in each gene are the same in C. elegans and C. briggsae. Primer-extension studies demonstrated that the two genes in C. briggsae are trans-splice differentially with the same splice leader (SL) RNAs as are observed in C. elegans. The gdp-2 gene is trans-spliced with SL1 while gdp-3 is trans-spliced with SL2. Significant sequence conservation was observed within the promoter regions of each species and may indicate those regions responsible for body-wall-muscle-specific gene expression and/or differential trans-splicing. Comparisons of the sequences suggest that the tandem repeat of the genes has been subjected to concerted evolution and that C. briggsae and C. elegans diverged much earlier than would be anticipated based on morphological similarities alone. Finally, an open reading frame found several hundred nucleotides upstream from gpd-2, in both species, appears to be homologous to the ATP synthase subunit, ATPase inhibitor protein, from bovine mitochondria.     

Genetic variation for postzygotic reproductive isolation between Caenorhabditis briggsae and Caenorhabditis sp. 9

The process of speciation is key to the origins of biodiversity, and yet the Caenorhabditis nematode model system has contributed little to this topic. Genetic studies of speciation in the genus are now feasible, owing to crosses between the recently discovered Caenorhabditis sp. 9 and the well-known C. briggsae producing fertile F(1) hybrid females. We dissected patterns of postzygotic reproductive isolation between these species by crossing eight isogenic strains of C. briggsae reciprocally with six strains of C. sp. 9. We determined that overall patterns of reproductive isolation are robust across these genetic backgrounds. However, we also quantified significant heritable variation within each species for interspecific hybrid incompatibilities for total adult progeny, egg-to-adult viability, and the percentage of male progeny. This demonstrates that intraspecific variation for interspecific hybrid incompatibility occurs despite extensive, albeit incomplete, reproductive isolation. Therefore, this emerging general phenomenon of variable reproductive isolation is not restricted to highly interfertile, early-stage incipient species, but also applies to species in the latest stages of the speciation process. Furthermore, we confirm Haldane's rule and demonstrate strongly asymmetric parent-of-origin effects (Darwin's corollary) that consistently manifest more extremely when hermaphroditic C. briggsae serves as maternal parent. These findings highlight Caenorhabditis as an emerging system for understanding the genetics of general patterns of reproductive isolation.     

Clustered organization of reproductive genes in the C. elegans genome

Defining the forces that sculpt genome organization is fundamental for understanding the origin, persistence, and diversification of species. The genomic sequences of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae provide an excellent opportunity to explore the dynamics of chromosome evolution. Extensive chromosomal rearrangement has accompanied divergence from their common ancestor, an event occurring roughly 100 million years ago (Mya); yet, morphologically, these species are nearly indistinguishable and both reproduce primarily by self-fertilization. Here, we show that genes expressed during spermatogenesis (sperm genes) are nonrandomly distributed across the C. elegans genome into three large clusters located on two autosomes. In addition to sperm genes, these chromosomal regions are enriched for genes involved in the hermaphrodite sperm/oocyte switch and in the reception of sperm signals that control fertilization. Most loci are present in single copy, suggesting that cluster formation is largely due to gene aggregation and not to tandem duplication. Comparative mapping indicates that the C. briggsae genome differs dramatically from the C. elegans genome in clustering. Because clustered genes have a direct role in reproduction and thus fitness, their aggregated pattern might have been shaped by natural selection, perhaps as hermaphroditism evolved.     

Patterns of molecular evolution in Caenorhabditis preclude ancient origins of selfing

The evolution of self-fertilization can mediate pronounced changes in genomes as a by-product of a drastic reduction in effective population size and the concomitant accumulation of slightly deleterious mutations by genetic drift. In the nematode genus Caenorhabditis, a highly selfing lifestyle has evolved twice independently, thus permitting an opportunity to test for the effects of mode of reproduction on patterns of molecular evolution on a genomic scale. Here we contrast rates of nucleotide substitution and codon usage bias among thousands of orthologous groups of genes in six species of Caenorhabditis, including the classic model organism Caenorhabditis elegans. Despite evidence that weak selection on synonymous codon usage is pervasive in the history of all species in this genus, we find little difference among species in the patterns of codon usage bias and in replacement-site substitution. Applying a model of relaxed selection on codon usage to the C. elegans and C. briggsae lineages suggests that self-fertilization is unlikely to have evolved more than approximately 4 million years ago, which is less than a quarter of the time since they shared a common ancestor with outcrossing species. We conclude that the profound changes in mating behavior, physiology, and developmental mechanisms that accompanied the transition from an obligately outcrossing to a primarily selfing mode of reproduction evolved in the not-too-distant past.     

Comparative genomics and adaptive selection of the ATP-binding-cassette gene family in caenorhabditis species

ABC transporters constitute one of the largest gene families in all species. They are mostly involved in transport of substrates across membranes. We have previously demonstrated that the Caenorhabditis elegans ABC family shows poor one-to-one gene orthology with other distant model organisms. To address the evolution dynamics of this gene family among closely related species, we carried out a comparative analysis of the ABC family among the three nematode species C. elegans, C. briggsae, and C. remanei. In contrast to the previous observations, the majority of ABC genes in the three species were found in orthologous trios, including many tandemly duplicated ABC genes, indicating that the gene duplication took place before speciation. Species-specific expansions of ABC members are rare and mostly observed in subfamilies A and B. C. briggsae and C. remanei orthologous ABC genes tend to cluster on trees, with those of C. elegans as an outgroup, consistent with their proposed species phylogeny. Comparison of intron/exon structures of the highly conserved ABCE subfamily members also indicates a closer relationship between C. briggsae and C. remanei than between either of these species and C. elegans. A comparison between insect and mammalian species indicates lineage-specific duplications or deletions of ABC genes, while the family size remains relatively constant. Sites undergoing positive selection within subfamily D, which are implicated in very-long-chain fatty acid transport, were identified. The evolution of these sites might be driven by the changes in food source with time.     

Functional elements and domains inferred from sequence comparisons of a heat shock gene in two nematodes

Caenorhabditis elegans and Caenorhabditis briggsae are two closely related nematode species that are nearly identical morphologically. Interspecific cross-hybridizing DNA appears to be restricted primarily to coding regions. We compared portions of the hsp-3 homologs, two grp 78-like genes, from C. elegans and C. briggsae and detected regions of DNA identity in the coding region, the 5' flanking DNAs, and the introns. The hsp-3 homologs share approximately 98% and 93% identity at the amino acid and nucleotide levels, respectively. Using the nucleotide substitution rate at the silent third position of the codons, we have estimated a lower limit for the date of divergence between C. elegans and C. briggsae to be approximately 23-32 million years ago. The 5' flanking DNAs and one of the introns contain elements that are highly conserved between C. elegans and C. briggsae. Some of the regions of nucleotide identity in the 5' flanking DNAs correspond to previously detected identities including viral enhancer sequences, a heat shock element, and an element present in the regulatory regions of mammalian grp78 and grp94 genes. We propose that a comparison of C. elegans and C. briggsae sequences will be useful in the detection of potential regulatory and structural elements.     

Conservation of sequence and function of the pag-3 genes from C. elegans and C. briggsae

The Caenorhabditis briggsae homologue of the Caenorhabditis elegans pag-3 gene was cloned and sequenced. When transformed into a C. elegans pag-3 mutant, the C. briggsae pag-3 gene rescued the pag-3 reverse kinker and lethargic phenotypes. The C. elegans pag-3 gene fused to lacZ was expressed in the same pattern in C. elegans and C. briggsae. Unlike many gene homologues compared between C. elegans and C. briggsae, extensive sequence conservation was found in the non-coding regions upstream of the pag-3 exons, in several of the introns and in the downstream non-coding region. Furthermore, the splice acceptor and splice donor sites were conserved, and the size of the introns and exons was surprisingly similar. The predicted protein sequence of C. briggsae PAG-3 was 85% identical to the protein sequence of C. elegans PAG-3. Because so much of the non-coding region of pag-3 was conserved, the control of pag-3 may be quite complex, involving the binding of many trans-acting factors. These results suggest the evolutionary conservation of the pag-3 gene sequence, its expression and function.     

Levels of DNA polymorphism vary with mating system in the nematode genus caenorhabditis

Self-fertilizing species often harbor less genetic variation than cross-fertilizing species, and at least four different models have been proposed to explain this trend. To investigate further the relationship between mating system and genetic variation, levels of DNA sequence polymorphism were compared among three closely related species in the genus Caenorhabditis: two self-fertilizing species, Caenorhabditis elegans and C. briggsae, and one cross-fertilizing species, C. remanei. As expected, estimates of silent site nucleotide diversity were lower in the two self-fertilizing species. For the mitochondrial genome, diversity in the selfing species averaged 42% of diversity in C. remanei. Interestingly, the reduction in genetic variation was much greater for the nuclear than for the mitochondrial genome. For two nuclear genes, diversity in the selfing species averaged 6 and 13% of diversity in C. remanei. We argue that either population bottlenecks or the repeated action of natural selection, coupled with high levels of selfing, are likely to explain the observed reductions in species-wide genetic diversity.