Patterns of nucleotide polymorphism distinguish temperate and tropical wild isolates of Caenorhabditis briggsae

Caenorhabditis briggsae provides a natural comparison species for the model nematode C. elegans, given their similar morphology, life history, and hermaphroditic mode of reproduction. Despite C. briggsae boasting a published genome sequence and establishing Caenorhabditis as a model genus for genetics and development, little is known about genetic variation across the geographic range of this species. In this study, we greatly expand the collection of natural isolates and characterize patterns of nucleotide variation for six loci in 63 strains from three continents. The pattern of polymorphisms reveals differentiation between C. briggsae strains found in temperate localities in the northern hemisphere from those sampled near the Tropic of Cancer, with diversity within the tropical region comparable to what is found for C. elegans in Europe. As in C. elegans, linkage disequilibrium is pervasive, although recombination is evident among some variant sites, indicating that outcrossing has occurred at a low rate in the history of the sample. In contrast to C. elegans, temperate regions harbor extremely little variation, perhaps reflecting colonization and recent expansion of C. briggsae into northern latitudes. We discuss these findings in relation to their implications for selection, demographic history, and the persistence of self-fertilization.     

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.     

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.     

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.     

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.     

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.     

Inbreeding and outbreeding depression in Caenorhabditis nematodes

The nematode Caenorhabditis elegans reproduces primarily by self-fertilization of hermaphrodites, yet males are present at low frequencies in natural populations (androdioecy). The ancestral state of C. elegans was probably gonochorism (separate males and females), as in its relative C. remanei. Males may be maintained in C. elegans because outcrossed individuals escape inbreeding depression. The level of inbreeding depression is, however, expected to be low in such a highly selfing species, compared with an outcrosser like C. remanei. To investigate these issues, we measured life-history traits in the progeny of inbred versus outcrossed C. elegans and C. remanei individuals derived from recently isolated natural populations. In addition, we maintained inbred lines of C. remanei through 13 generations of full-sibling mating. Highly inbred C. remanei showed dramatic reductions in brood size and relative fitness compared to outcrossed individuals, with evidence of both direct genetic and maternal-effect inbreeding depression. This decline in fitness accumulated over time, causing extinction of nearly 90% of inbred lines, with no evidence of purging of deleterious mutations from the remaining lines. In contrast, pure strains of C. elegans performed better than crosses between strains, indicating outbreeding depression. The results are discussed in relation to the evolution of androdioecy and the effect of mating system on the level of inbreeding depression.     

Conservation of the C.elegans tra-2 3'UTR translational control.

The Caenorhabditis elegans sex-determination gene, tra-2, is translationally regulated by two 28 nt elements (DREs) located in the 3'UTR that bind a factor called DRF. This regulation requires the laf-1 gene activity. We demonstrate that the nematode Caenorhabditis briggsae tra-2 gene and the human oncogene GLI are translationally regulated by elements that are functionally equivalent to DREs. Here, we rename the DREs to TGEs (tra-2 and GLI elements). Similarly to the C.elegans tra-2 TGEs, the C.briggsae tra-2 and GLI TGEs repress translation of a reporter transgene in a laf-1 dependent manner. Furthermore, they regulate poly(A) tail length and bind DRF. We also find that the C.elegans TGEs control translation and poly(A) tail length in C.briggsae and rodent cells. Moreover, these same organisms contain a factor that specifically associates with the C.elegans TGEs. These findings are consistent with the TGE control being present in C.briggsae and rodent cells. Three lines of evidence indicate that C.briggsae tra-2 and GLI are translationally controlled in vivo by TGEs. First, like C.elegans tra-2 TGEs, the C.briggsae tra-2 and GLI TGEs control translation and poly(A) tail lengths in C.briggsae and rodent cells, respectively. Second, the same factor in C.briggsae and mammalian cells that binds to the C.elegans tra-2 TGEs binds the C.briggsae tra-2 and GLI TGEs. Third, deletion of the GLI TGE increases GLI's ability to transform cells. These findings suggest that TGE control is conserved and regulates the expression of other mRNAs.     

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.     

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.     

Differences in collagen prolyl 4-hydroxylase assembly between two Caenorhabditis nematode species despite high amino acid sequence identity of the enzyme subunits.

The collagen prolyl 4-hydroxylases (P4Hs) are essential for proper extracellular matrix formation in multicellular organisms. The vertebrate enzymes are alpha(2)beta(2) tetramers, in which the beta subunits are identical to protein disulfide isomerase (PDI). Unique P4H forms have been shown to assemble from the Caenorhabditis elegans catalytic alpha subunit isoforms PHY-1 and PHY-2 and the beta subunit PDI-2. A mixed PHY-1/PHY-2/(PDI-2)(2) tetramer is the major form, while PHY-1/PDI-2 and PHY-2/PDI-2 dimers are also assembled but less efficiently. Cloning and characterization of the orthologous subunits from the closely related nematode Caenorhabditis briggsae revealed distinct differences in the assembly of active P4H forms in spite of the extremely high amino acid sequence identity (92-97%) between the C. briggsae and C. elegans subunits. In addition to a PHY-1/PHY-2(PDI-2)(2) tetramer and a PHY-1/PDI-2 dimer, an active (PHY-2)(2)(PDI-2)(2) tetramer was formed in C. briggsae instead of a PHY-2/PDI-2 dimer. Site-directed mutagenesis studies and generation of inter-species hybrid polypeptides showed that the N-terminal halves of the Caenorhabditis PHY-2 polypeptides determine their assembly properties. Genetic disruption of C. briggsae phy-1 (Cb-dpy-18) via a Mos1 insertion resulted in a small (short) phenotype that is less severe than the dumpy (short and fat) phenotype of the corresponding C. elegans mutants (Ce-dpy-18). C. briggsae phy-2 RNA interference produced no visible phenotype in the wild type nematodes but produced a severe dumpy phenotype and larval arrest in phy-1 mutants. Genetic complementation of the C. briggsae and C. elegans phy-1 mutants was achieved by injection of a wild type phy-1 gene from either species.     

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.     

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.     

Cryptic quantitative evolution of the vulva intercellular signaling network in Caenorhabditis

BACKGROUND: The Caenorhabditis vulva is formed from a row of Pn.p precursor cells, which adopt a spatial cell-fate pattern-3 degrees 3 degrees 2 degrees 1 degrees 2 degrees 3 degrees -centered on the gonadal anchor cell. This pattern is robustly specified by an intercellular signaling network including EGF/Ras induction from the anchor cell and Delta/Notch signaling between the precursor cells. It is unknown how the roles and quantitative contributions of these signaling pathways have evolved in closely related Caenorhabditis species. RESULTS: Cryptic evolution in the network is uncovered by quantification of cell-fate-pattern frequencies obtained after displacement of the system out of its normal range, either by anchor-cell ablations or through LIN-3/EGF overexpression. Silent evolution in the Caenorhabditis genus covers a large neutral space of cell-fate patterns. Direct induction of the 1 degrees fate as in C. elegans appeared within the genus. C. briggsae displays a graded induction of 1 degrees and 2 degrees fates, with 1 degrees fate induction requiring a longer time than in C. elegans, and a reduced lateral inhibition of adjacent 1 degrees fates. C. remanei displays a strong lateral induction of 2 degrees fates relative to vulval-fate activation in the central cell. This evolution in cell-fate pattern space can be experimentally reconstituted by mild variations of Ras, Wnt, and Notch pathway activities in C. elegans and C. briggsae. CONCLUSIONS: Quantitative evolution in the roles of graded induction by LIN-3/EGF and Notch signaling is demonstrated for the Caenorhabditis vulva signaling network. This evolutionary system biology approach provides a quantitative view of the variational properties of this biological system.     

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.     

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.     

Haldane's rule by sexual transformation in Caenorhabditis

Haldane's rule in C. briggsae x C. remanei broods was caused by sexual transformation; XX and XO hybrids were female. C. briggsae and C. remanei variants that partially suppress hybrid sexual transformation were identified. Effects of variant strains were cumulative. Hence, aberrant sex determination is a reproductive isolation mechanism in Caenorhabditis.