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.     

Codon usage in Caenorhabditis elegans: delineation of translational selection and mutational biases.

Synonymous codon usage varies considerably among Caenorhabditis elegans genes. Multivariate statistical analyses reveal a single major trend among genes. At one end of the trend lie genes with relatively unbiased codon usage. These genes appear to be lowly expressed, and their patterns of codon usage are consistent with mutational biases influenced by the neighbouring nucleotide. At the other extreme lie genes with extremely biased codon usage. These genes appear to be highly expressed, and their codon usage seems to have been shaped by selection favouring a limited number of translationally optimal codons. Thus, the frequency of these optimal codons in a gene appears to be correlated with the level of gene expression, and may be a useful indicator in the case of genes (or open reading frames) whose expression levels (or even function) are unknown. A second, relatively minor trend among genes is correlated with the frequency of G at synonymously variable sites. It is not yet clear whether this trend reflects variation in base composition (or mutational biases) among regions of the C.elegans genome, or some other factor. Sequence divergence between C.elegans and C.briggsae has also been studied.     

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.     

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.     

Comparative genetics of sex determination: masculinizing mutations in Caenorhabditis briggsae

The nematodes Caenorhabditis elegans and C. briggsae independently evolved self-fertile hermaphroditism from gonochoristic ancestors. C. briggsae has variably divergent orthologs of nearly all genes in the C. elegans sex determination pathway. Their functional characterization has generally relied on reverse genetic approaches, such as RNA interference and cross-species transgene rescue and more recently on deletion mutations. We have taken an unbiased forward mutagenesis approach to isolating zygotic mutations that masculinize all tissues of C. briggsae hermaphrodites. The screens identified loss-of-function mutations in the C. briggsae orthologs of tra-1, tra-2, and tra-3. The somatic and germline phenotypes of these mutations are largely identical to those of their C. elegans homologs, including the poorly understood germline feminization of tra-1(lf) males. This overall conservation of Cb-tra phenotypes is in contrast to the fem genes, with which they directly interact and which are significantly divergent in germline function. In addition, we show that in both C. briggsae and C. elegans large C-terminal truncations of TRA-1 that retain the DNA-binding domain affect sex determination more strongly than somatic gonad development. Beyond these immediate results, this collection of mutations provides an essential foundation for further comparative genetic analysis of the Caenorhabditis sex determination pathway.     

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.     

Pattern of selective constraint in C. elegans and C. briggsae genomes.

Similarity between related genomes may carry information on selective constraint in each of them. We analysed patterns of similarity between several homologous regions of Caenorhabditis elegans and C. briggsae genomes. All homologous exons are quite similar. Alignments of introns and of intergenic sequences contain long gaps, segments where similarity is low and close to that between random sequences aligned using the same parameters, and segments of high similarity. Conservative estimates of the fractions of selectively constrained nucleotides are 72%, 17% and 18% for exons, introns and intergenic sequences, respectively. This implies that the total number of constrained nucleotides within non-coding sequences is comparable to that within coding sequences, so that at least one-third of nucleotides in C. elegans and C. briggsae genomes are under strong stabilizing selection.     

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.     

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.     

Comparative peptidomics of Caenorhabditis elegans versus C. briggsae by LC-MALDI-TOF MS

Neuropeptides are important signaling molecules that function in cell-cell communication as neurotransmitters or hormones to orchestrate a wide variety of physiological conditions and behaviors. These endogenous peptides can be monitored by high throughput peptidomics technologies from virtually any tissue or organism. The neuropeptide complement of the soil nematode Caenorhabditis elegans has been characterized by on-line two-dimensional liquid chromatography and quadrupole time-of-flight tandem mass spectrometry (2D-nanoLC Q-TOF MS/MS). Here, we use an alternative peptidomics approach combining liquid chromatography (LC) with matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry to map the peptide content of C. elegans and another Caenorhabditis species, Caenorhabditis briggsae. This study allows a better annotation of neuropeptide-encoding genes from the C. briggsae genome and provides a promising basis for further evolutionary comparisons.     

Prediction of structured non-coding RNAs in the genomes of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae

We present a survey for non-coding RNAs and other structured RNA motifs in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae using the RNAz program. This approach explicitly evaluates comparative sequence information to detect stabilizing selection acting on RNA secondary structure. We detect 3,672 structured RNA motifs, of which only 678 are known non-translated RNAs (ncRNAs) or clear homologs of known C. elegans ncRNAs. Most of these signals are located in introns or at a distance from known protein-coding genes. With an estimated false positive rate of about 50% and a sensitivity on the order of 50%, we estimate that the nematode genomes contain between 3,000 and 4,000 RNAs with evolutionary conserved secondary structures. Only a small fraction of these belongs to the known RNA classes, including tRNAs, snoRNAs, snRNAs, or microRNAs. A relatively small class of ncRNA candidates is associated with previously observed RNA-specific upstream elements.     

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.     

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.     

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.     

Hill-Robertson interference is a minor determinant of variations in codon bias across Drosophila melanogaster and Caenorhabditis elegans genomes

According to population genetics models, genomic regions with lower crossing-over rates are expected to experience less effective selection because of Hill-Robertson interference (HRi). The effect of genetic linkage is thought to be particularly important for a selection of weak intensity such as selection affecting codon usage. Consistent with this model, codon bias correlates positively with recombination rate in Drosophila melanogaster and Caenorhabditis elegans. However, in these species, the G+C content of both noncoding DNA and synonymous sites correlates positively with recombination, which suggests that mutation patterns and recombination are associated. To remove this effect of mutation patterns on codon bias, we used the synonymous sites of lowly expressed genes that are expected to be effectively neutral sites. We measured the differences between codon biases of highly expressed genes and their lowly expressed neighbors. In D. melanogaster we find that HRi weakly reduces selection on codon usage of genes located in regions of very low recombination; but these genes only comprise 4% of the total. In C. elegans we do not find any evidence for the effect of recombination on selection for codon bias. Computer simulations indicate that HRi poorly enhances codon bias if the local recombination rate is greater than the mutation rate. This prediction of the model is consistent with our data and with the current estimate of the mutation rate in D. melanogaster. The case of C. elegans, which is highly self-fertilizing, is discussed. Our results suggest that HRi is a minor determinant of variations in codon bias across the genome.