Evolution of discrete Notch-like receptors from a distant gene duplication in Caenorhabditis

SUMMARY Caenorhabditis elegans possesses two Notch-like receptors, LIN-12 and GLP-1, which have both overlapping and individual biological functions. We examined the lin-12 and glp-1 genes in closely related nematodes to learn about their evolution. Here we report molecular and functional analyses of lin-12 orthologs from two related nematodes, C. briggsae (Cb) and C. remanei (Cr). In addition, we compare these lin-12 findings with similar studies of Cb-glp-1 and Cr-glp-1 orthologs. Cb-LIN-12 and Cr-LIN-12 retain the same number and order of motifs as Ce-LIN-12. Intriguingly, we find that LIN-12 conservation differs from that of GLP-1 in two respects. First, individual motifs are conserved to a different degree for the two receptors. For example, the transmembrane domain is 16-32% identical among LIN-12 orthologs but 65-70% identical among GLP-1 orthologs. Second, certain amino acids are conserved in a receptor-specific manner, a phenomenon most prevalent in the CC-linker. We suggest that LIN-12 and GLP-1 have been molded by selective constraints that are receptor specific and that the two proteins may not be entirely interchangeable. To analyze the functions of the lin-12 orthologs, we used RNA-mediated interference (RNAi). Cb-lin-12(RNAi) or Cr-lin-12(RNAi) progeny are nearly 100% Lag, a larval lethality typical of C. elegans lin-12 glp-1 double mutants, but not the primary defect observed in Ce-lin-12 null mutants or Ce-lin-12(RNAi). Therefore, LIN-12 functions are similar, but not identical, among the Caenorhabditis species. We suggest that ancestral functions may have been divided between LIN-12 and GLP-1 receptors in a process contributing to the retention of both genes after gene duplication (i.e., subfunctionalization).     

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.     

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.     

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.     

Evolution of regulatory elements producing a conserved gene expression pattern in Caenorhabditis

Natural selection acts at the level of function, not at the logistical level of how organisms achieve a particular function. Consequently, significant DNA sequence and regulatory differences can achieve the same function, such as a particular gene expression pattern. To investigate how regulatory features underlying a conserved function can evolve, we compared the regulation of a conserved gene expression pattern in the related species Caenorhabditis elegans and C. briggsae. We find that both C. elegans and C. briggsae express the ovo-related zinc finger gene lin-48 in the same pattern in hindgut cells. However, the regulation of this gene by the Pax-2/5/8 protein EGL-38 differs in two important ways. First, specific differences in the regulatory sequences of lin-48 result in the presence of two redundant EGL-38 response elements in C. elegans, whereas the redundancy is absent in C. briggsae. Second, there is a single egl-38 gene in C. briggsae. In contrast, the gene is duplicated in C. elegans, with only one copy retaining the ability to regulate lin-48 in vivo. These results illustrate molecular changes that can occur despite maintenance of conserved gene function in different 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.     

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.     

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.     

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.     

Evolutionarily conserved regions in Caenorhabditis transposable elements deduced by sequence comparison.

In this paper we present the sequence of an intact Caenorhabditis briggsae transposable element, Tcb2. Tcb2 is 1606 base pairs in length and contains 80 base pair imperfect terminal repeats and a single open reading frame. We have identified blocks of T-rich repeats in the regions 150-200 and 1421-1476 of this element which are conserved in the Caenorhabditis elegans element Tc1. The sequence conservation of these regions in elements from different Caenorhabditis species suggests that they are of functional importance. A single open reading frame corresponding to the major open reading frame of Tc1 is conserved among Tc1, Tcb1, and Tcb2. Comparison of the first 550 nucleotides of the sequence among the three elements has allowed the evaluation of a model proposing an extension of the major open reading frame. Our data support the suggestion that Tc1 is capable of producing a 335 amino acid protein. A comparison of the sequence coding for the amino and carboxy termini of the 273 amino acid transposase from Caenorhabditis Tc1-like elements and Drosophila HB1 showed different amounts of divergence for each of these regions, indicating that the two functional domains have undergone different amounts of selection. Our data are not compatible with the proposal that Tc1-related sequences have been acquired via horizontal transmission. The divergence of Tc1 from the two C. briggsae elements, Tcb1 and Tcb2, indicated that all three elements have been diverging from each other for approximately the same amount of time as the genomes of the two species.     

The genome of Oscheius tipulae: determination of size, complexity, and structure by DNA reassociation using fluorescent dye.

This work describes the physicochemical characterization of the genome and telomere structure from the nematode Oscheius tipulae CEW1. Oscheius tipulae is a free-living nematode belonging to the family Rhabditidae and has been used as a model system for comparative genetic studies. A new protocol that combines fluorescent detection of double-stranded DNA and S1 nuclease was used to determine the genome size of O. tipulae as 100.8 Mb (approximately 0.1 pg DNA/haploid nucleus). The genome of this nematode is made up of 83.4% unique copy sequences, 9.4% intermediate repetitive sequences, and 7.2% highly repetitive sequences, suggesting that its structure is similar to those of other nematodes of the genus Caenorhabditis. We also showed that O. tipulae has the same telomere repeats already found in Caenorhabditis elegans at the ends and in internal regions of the chromosomes. Using a cassette-ligation-mediated PCR protocol we were able to obtain 5 different putative subtelomeric sequences of O. tipulae, which show no similarity to C. elegans or C. briggsae subtelomeric regions. DAPI staining of hermaphrodite gonad cells show that, as detected in C. elegans and other rhabditids, O. tipulae have a haploid complement of 6 chromosomes.     

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.     

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.