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.     

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.     

Regulatory elements required for development of caenorhabditis elegans hermaphrodites are conserved in the tra-2 homologue of C. remanei, a male/female sister species

The Caenorhabditis elegans hermaphrodite is essentially a female that produces sperm. In C. elegans, tra-2 promotes female fates and must be repressed to achieve hermaphrodite spermatogenesis. In an effort to learn how mating systems evolve, we have cloned tra-2 from C. remanei, the closest gonochoristic relative of C. elegans. We found its structure to be similar to that of Ce-tra-2 but its sequence to be divergent. RNA interference demonstrates that Cr-tra-2 promotes female fates. Two sites of tra-2 regulation are required for the onset of hermaphrodite spermatogenesis in C. elegans. One, the MX region of TRA-2, is as well conserved in C. remanei as it is in C. briggsae (another male/hermaphrodite species), suggesting that this control is not unique to hermaphrodites. Another, the DRE/TGE element of the tra-2 3' UTR, was not detected by sequence analysis. However, gel-shift assays demonstrate that a factor in C. remanei can bind specifically to the Cr-tra-2 3' UTR, suggesting that this translational control is also conserved. We propose that both controls are general and do not constitute a novel "switch" that enables sexual mosaicism in hermaphrodites. However, subtle quantitative or qualitative differences in their employment may underlie differences in mating system seen in Caenorhabditis.     

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.     

Conservation of glp-1 regulation and function in nematodes

The Caenorhabditis elegans (Ce) glp-1 gene encodes a Notch-like receptor. We have cloned glp-1 from C. briggsae (Cb) and C. remanei (Cr), two Caenorhabditis species that have diverged from C. elegans by roughly 20-40 million years. By sequence analysis, we find that the Cb-GLP-1 and Cr-GLP-1 proteins have retained the same motif architecture as Ce-GLP-1, including number of domains. In addition, two regions (CC-linker and regions flanking the ANK repeats) are as highly conserved as regions previously recognized as essential for signaling (e.g., ANK repeats). Phylogenetic analysis of glp-1 sequences suggests a C. briggsae/C. remanei clade with C. elegans as a sister taxon. Using RNAi to test biological functions, we find that Ce-glp-1, Cb-glp-1, and Cr-glp-1 are all required for proliferation of germline stem cells and for specifying blastomere fates in the embryo. In addition, certain biological roles of Cb-glp-1, e.g., in the vulva, have diverged from those of Ce-glp-1 and Cr-glp-1, suggesting a change in either regulation or function of the Cb-glp-1 gene during evolution. Finally, the regulation of glp-1 mRNA, previously analyzed for Ce-glp-1, is conserved in Cb-glp-1, and we identify conserved 3' UTR sequences that may serve as regulatory elements.     

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.     

Evidence for active maintenance of inverted repeat structures identified by a comparative genomic approach

Inverted repeats have been found to occur in both prokaryotic and eukaryotic genomes. Usually they are short and some have important functions in various biological processes. However, long inverted repeats are rare and can cause genome instability. Analyses of C. elegans genome identified long, nearly-perfect inverted repeat sequences involving both divergently and convergently oriented homologous gene pairs and complete intergenic sequences. Comparisons with the orthologous regions from the genomes of C. briggsae and C. remanei show that the inverted repeat structures are often far more conserved than the sequences. This observation implies that there is an active mechanism for maintaining the inverted repeat nature of the sequences.     

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.     

Rapid coevolution of the nematode sex-determining genes fem-3 and tra-2

Unlike many features of metazoan development, sex determination is not widely conserved among phyla. However, the recent demonstration that one gene family controls sexual development in Drosophila, C. elegans, and vertebrates suggests that sex determination mechanisms may have evolved from a common pathway that has diverged radically since the Cambrian. Sex determination gene sequences often evolve quickly, but it is not known how this relates to higher-order pathways or what selective or neutral forces are driving it. In such a rapidly evolving developmental pathway, the fate of functionally linked genes is of particular interest. To investigate a pair of such genes, we cloned orthologs of the key C. elegans male-promoting gene fem-3 from two sister species, C. briggsae and C. remanei. We employed RNA interference to show that in all three species, the male-promoting function of fem-3 and its epistatic relationship with its female-promoting upstream repressor, tra-2, are conserved. Consistent with this, the FEM-3 protein interacts with TRA-2 in each species, but in a strictly species-specific manner. Because FEM-3 is the most divergent protein yet described in Caenorhabditis and the FEM-3 binding domain of TRA-2 is itself hypervariable, a key protein-protein interaction is rapidly evolving in concert. Extrapolation of this result to larger phylogenetic scales helps explain the dissimilarity of the sex determination systems across phyla.     

A family of zinc-finger proteins is required for chromosome-specific pairing and synapsis during meiosis in C. elegans

Homologous chromosome pairing and synapsis are prerequisite for accurate chromosome segregation during meiosis. Here, we show that a family of four related C2H2 zinc-finger proteins plays a central role in these events in C. elegans. These proteins are encoded within a tandem gene cluster. In addition to the X-specific HIM-8 protein, three additional paralogs collectively mediate the behavior of the five autosomes. Each chromosome relies on a specific member of the family to pair and synapse with its homolog. These "ZIM" proteins concentrate at special regions called meiotic pairing centers on the corresponding chromosomes. These sites are dispersed along the nuclear envelope during early meiotic prophase, suggesting a role analogous to the telomere-mediated meiotic bouquet in other organisms. To gain insight into the evolution of these components, we characterized homologs in C. briggsae and C. remanei, which revealed changes in copy number of this gene family within the nematode lineage.     

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.     

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.     

REPRODUCTIVE ISOLATION IN RHABDITIDAE (NEMATODA: SECERNENTEA); MECHANISMS THAT ISOLATE SIX SPECIES OF THREE GENERA

We have attempted interspecific hybridizations among six species of rhabditid nematodes: Caenorhabditis elegans, Caenorhabditis briggsae, Caenorhabditis remanei, Caenorhabditis sp. v, Rhabditis sp., and Pelodera teres. Copulation was observed in all crosses between Caenorhabditis species; however, none resulted in the generation of stable hybrid populations. No copulation was observed in crosses between Caenorhabditis males and Rhabditis or Pelodera females, even when congeneric females were present, suggesting that Caenorhabditis males are able to selectively recognize congeneric females by a short-range stimulus. All pairwise combinations of Caenorhabditis species were isolated to some degree by gametic mechanisms; 7 of 12 combinations were cross infertile and 5 of 12 were cross-fertile but had low brood sizes. In cross-fertile combinations, most hybrid embryos were inviable and arrested prior to gastrulation. Only in crosses of C. briggsae males to C. sp. v females did any hybrids survive embryogenesis. Most of these C. briggsae/C. sp. v hybrids arrested during larval development, and the few that reached adulthood invariably were female. These results are consistent with the presence of at least two lethal factors in the C. briggsae-C. sp. v combination: a maternal lethal factor in the cytoplasm of C. briggsae and a recessive lethal factor on the X chromosome of C. sp. v.     

The egg surface LDL receptor repeat-containing proteins EGG-1 and EGG-2 are required for fertilization in Caenorhabditis elegans

The molecular machinery that mediates sperm-egg interactions at fertilization is largely unknown. We identify two partially redundant egg surface LDL receptor repeat-containing proteins (EGG-1 and EGG-2) that are required for Caenorhabditis elegans fertility in hermaphrodites, but not males. Wild-type sperm cannot enter the morphologically normal oocytes produced by hermaphrodites that lack egg-1 and egg-2 function despite direct gamete contact. Furthermore, we find that levels of meiotic maturation/ovulation and sperm migratory behavior are altered in egg-1 mutants. These observations suggest an unexpected regulatory link between fertilization and other events necessary for reproductive success. egg-1 and egg-2 are the result of a gene duplication in the nematode lineage leading to C. elegans. The two closely related species C. briggsae and C. remanei encode only a single egg-1/egg-2 homolog that is required for hermaphrodite/female fertility. In addition to being the first identified egg components of the nematode fertilization machinery, the egg-1 and egg-2 gene duplication could be vital with regards to maximizing C. elegans fecundity and understanding the evolutionary differentiation of molecular function and speciation.     

Mode and Tempo of Molecular Evolution in the Nematode Caenorhabditis: Cytochrome Oxidase II and Calmodulin Sequences

Through direct sequencing methods, the mitochondrial gene for cytochrome oxidase subunit two (CO II) and the single-copy nuclear gene for calmodulin were compared among strains of Caenorhabidits elegans and two other Caenorhabditis species (C. remanei and C. briggsae). In addition the CO II sequence was determined from a distantly related nematode, Steinernema intermedii. Among the 11 strains of C. elegans tested, there are four types of CO II gene, arising from two major lineages. Levels of intraspecific difference in the CO II gene are low (less than 2.0%) compared to the extraordinary divergence between congeneric species, which is about 50% when corrected for multiple hits. Concordant with the increase in divergence between taxa is a change in the pattern of substitution from a strong transition bias (24 transitions compared to two transversions) within species to a substitution pattern that appears to reflect the base composition of the mitochondrial genome when more divergent nematodes are compared. The base composition of the Caenorhabditis CO II gene is strongly biased toward A + T at all three positions of codons and appears to constrain the amino acid composition of the protein. Both the CO II and calmodulin genes show extreme conservation of amino acid sequences. When the accumulation of changes at silent sites in the two genes is compared among strains, it becomes evident that the mitochondrial gene is changing faster than the nuclear gene.     

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.     

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.