Search

Studies of model organisms have provided important insights into how natural genetic differences shape trait variation. These discoveries are driven by the growing availability of genomes and the expansive experimental toolkits afforded to researchers using these species. For example, Caenorhabditis elegans is increasingly being used to identify and measure the effects of natural genetic variants on traits using quantitative genetics. Since 2016, the C. elegans Natural Diversity Resource (CeNDR) has facilitated many of these studies by providing an archive of wild strains, genome-wide sequence and variant data for each strain, and a genome-wide association (GWA) mapping portal for the C. elegans community. Here, we present an updated platform, the Caenorhabditis Natural Diversity Resource (CaeNDR), that enables quantitative genetics and genomics studies across the three Caenorhabditis species: C. elegans, C. briggsae and C. tropicalis. The CaeNDR platform hosts several databases that are continually updated by the addition of new strains, whole-genome sequence data and annotated variants. Additionally, CaeNDR provides new interactive tools to explore natural variation and enable GWA mappings. All CaeNDR data and tools are accessible through a freely available web portal located at caendr.org., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared to insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. F-box genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Moreover, using both neutrality tests and Extended Haplotype Homozygosity analysis, we identified signatures of strong positive selection in the F-box and csGPCR genes among the wild isolates, especially in the non-Hawaiian population. Accumulation of high-frequency derived alleles in these genes was found in non-Hawaiian population, leading to divergence from the ancestral genotype. In summary, we found that F-box and csGPCR genes harbour a large pool of natural variants, which may be subjected to positive selection. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular and intracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Caenorhabditis auriculariae, which was morphologically described in 1999, was re-isolated from a Platydema mushroom-associated beetle. Based on the re-isolated materials, some morphological characteristics were re-examined and ascribed to the species. In addition, to clarify phylogenetic relationships with other Caenorhabditis species and biological features of the nematode, the whole genome was sequenced and assembled into 109.5 Mb with 16,279 predicted protein-coding genes. Molecular phylogenetic analyses based on ribosomal RNA and 269 single-copy genes revealed the species is closely related to C. sonorae and C. monodelphis placing them at the most basal clade of the genus. C. auriculariae has morphological characteristics clearly differed from those two species and harbours a number of species-specific gene families, indicating its usefulness as a new outgroup species for Caenorhabditis evolutionary studies. A comparison of carbohydrate-active enzyme (CAZy) repertoires in genomes, which we found useful to speculate about the lifestyle of Caenorhabditis nematodes, suggested that C. auriculariae likely has a life-cycle with tight-association with insects.

Toxin-antidote elements (TAs) are selfish genetic dyads that spread in populations by selectively killing non-carriers. TAs are common in prokaryotes, but very few examples are known in animals. Here, we report the discovery of maternal-effect TAs in both C. tropicalis and C. briggsae, two distant relatives of C. elegans. In C. tropicalis, multiple TAs combine to cause a striking degree of intraspecific incompatibility: five elements reduce the fitness of >70% of the F 2 hybrid progeny of two Caribbean isolates. We identified the genes underlying one of the novel TAs, slow-1/grow-1, and found that its toxin, slow-1, is homologous to nuclear hormone receptors. Remarkably, although previously known TAs act during embryonic development, maternal loading of slow-1 in oocytes specifically slows down larval development, delaying the onset of reproduction by several days. Finally, we found that balancing selection acting on linked, conflicting TAs hampers their ability to spread in populations, leading to more stable genetic incompatibilities. Our findings indicate that TAs are widespread in Caenorhabditis species and target a wide range of developmental processes and that antagonism between them may cause lasting incompatibilities in natural populations. We expect that similar phenomena exist in other animal species., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

The nematode Caenorhabditis elegans is one of the key model systems in biology, including possessing the first fully assembled animal genome. Whereas C. elegans is a self-reproducing hermaphrodite with fairly limited within-population variation, its relative C. remanei is an outcrossing species with much more extensive genetic variation, making it an ideal parallel model system for evolutionary genetic investigations. Here, we greatly improve on previous assemblies by generating a chromosome-level assembly of the entire C. remanei genome (124.8 Mb of total size) using long-read sequencing and chromatin conformation capture data. Like other fully assembled genomes in the genus, we find that the C. remanei genome displays a high degree of synteny with C. elegans despite multiple within-chromosome rearrangements. Both genomes have high gene density in central regions of chromosomes relative to chromosome ends and the opposite pattern for the accumulation of repetitive elements. C. elegans and C. remanei also show similar patterns of interchromosome interactions, with the central regions of chromosomes appearing to interact with one another more than the distal ends. The new C. remanei genome presented here greatly augments the use of the Caenorhabditis as a platform for comparative genomics and serves as a basis for molecular population genetics within this highly diverse species., (Copyright © 2020 by the Genetics Society of America.)

The free-living nematode Caenorhabditis elegans is a key laboratory model for metazoan biology. C. elegans has also become a model for parasitic nematodes despite being only distantly related to most parasitic species. All of the ∼65 Caenorhabditis species currently in culture are free-living, with most having been isolated from decaying plant or fungal matter. Caenorhabditis bovis is a particularly unusual species that has been isolated several times from the inflamed ears of Zebu cattle in Eastern Africa, where it is associated with the disease bovine parasitic otitis. C. bovis is therefore of particular interest to researchers interested in the evolution of nematode parasitism. However, as C. bovis is not in laboratory culture, it remains little studied. Here, by sampling livestock markets and slaughterhouses in Western Kenya, we successfully reisolated C. bovis from the ear of adult female Zebu. We sequenced the genome of C. bovis using the Oxford Nanopore MinION platform in a nearby field laboratory and used the data to generate a chromosome-scale draft genome sequence. We exploited this draft genome sequence to reconstruct the phylogenetic relationships of C. bovis to other Caenorhabditis species and reveal the changes in genome size and content that have occurred during its evolution. We also identified expansions in several gene families that have been implicated in parasitism in other nematode species. The high-quality draft genome and our analyses thereof represent a significant advancement in our understanding of this unusual Caenorhabditis species., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

UDP- glucose: glycoprotein glucosyltransferase (UGGT) is a protein that operates as the gatekeeper for the endoplasmic reticulum (ER) quality control mechanism of glycoprotein folding. It is known that vertebrates and Caenorhabditis genomes harbor two uggt gene copies that exhibit differences in their properties.Bayesian phylogenetic inference based on 195 UGGT and UGGT-like protein sequences of an ample spectrum of eukaryotic species showed that uggt genes went through independent duplications in Caenorhabditis and vertebrates. In both lineages, the catalytic domain of the duplicated genes was subjected to a strong purifying selective pressure, while the recognition domain was subjected to episodic positive diversifying selection. Selective relaxation in the recognition domain was more pronounced in Caenorhabditis uggt-b than in vertebrates uggt-2 Structural bioinformatics analysis revealed that Caenorhabditis UGGT-b protein lacks essential sequences proposed to be involved in the recognition of unfolded proteins. When we assayed glucosyltrasferase activity of a chimeric protein composed by Caenorhabditis uggt-b recognition domain fused to S. pombe catalytic domain expressed in yeast, no activity was detected.The present results support the conservation of the UGGT activity in the catalytic domain and a putative divergent function of the recognition domain for the UGGT2 protein in vertebrates, which would have gone through a specialization process. In Caenorhabditis , uggt-b evolved under different constraints compared to uggt-a which, by means of a putative neofunctionalization process, resulted in a non-redundant paralog. The non-canonical function of uggt-b in the worm lineage highlights the need to take precautions before generalizing gene functions in model organisms., (Copyright © 2020 Caraballo et al.)

Dispersal often covaries with other traits, and this covariation was shown to have a genetic basis. Here, we wanted to explore to what extent genetic constraints and correlational selection can explain patterns of covariation between dispersal and key life-history traits-lifespan and reproduction. A prediction from the fitness-associated dispersal hypothesis was that lower genetic quality is associated with higher dispersal propensity as driven by the benefits of genetic mixing. We wanted to contrast it with a prediction from a different model that individuals putting more emphasis on current rather than future reproduction disperse more, as they are expected to be more risk-prone and exploratory. However, if dispersal has inherent costs, this will also result in a negative genetic correlation between higher rates of dispersal and some aspects of performance. To explore this issue, we used the dioecious nematode Caenorhabditis remanei and selected for increased and decreased dispersal propensity for 10 generations, followed by five generations of relaxed selection. Dispersal propensity responded to selection, and females from high-dispersal lines dispersed more than females from low-dispersal lines. Females selected for increased dispersal propensity produced fewer offspring and were more likely to die from matricide, which is associated with a low physiological condition in Caenorhabditis nematodes. There was no evidence for differences in age-specific reproductive effort between high- and low-dispersal females. Rather, reproductive output of high-dispersal females was consistently reduced. We argue that our data provide support for the fitness-associated dispersal hypothesis., (© 2019 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2019 European Society For Evolutionary Biology.)

Three RNA viruses related to nodaviruses were previously described to naturally infect the nematode Caenorhabditis elegans and its relative, Caenorhabditis briggsae Here, we report on a collection of more than 50 viral variants from wild-caught Caenorhabditis. We describe the discovery of a new related virus, the Mělník virus, infecting C. briggsae , which similarly infects intestinal cells. In France, a frequent pattern of coinfection of C. briggsae by the Santeuil virus and Le Blanc virus was observed at the level of an individual nematode and even a single cell. We do not find evidence of reassortment between the RNA1 and RNA2 molecules of Santeuil and Le Blanc viruses. However, by studying patterns of evolution of each virus, reassortments of RNA1 and RNA2 among variants of each virus were identified. We develop assays to test the relative infectivity and competitive ability of the viral variants and detect an interaction between host genotype and Santeuil virus genotype, such that the result depends on the host strain. IMPORTANCE The roundworm Caenorhabditis elegans is a laboratory model organism in biology. We study natural populations of this small animal and its relative, C. briggsae , and the viruses that infect them. We previously discovered three RNA viruses related to nodaviruses and here describe a fourth one, called the Mělník virus. These viruses have a genome composed of two RNA molecules. We find that two viruses may infect the same animal and the same cell. The two RNA molecules may be exchanged between variants of a given viral species. We study the diversity of each viral species and devise an assay of their infectivity and competitive ability. Using this assay, we show that the outcome of the competition also depends on the host., (Copyright © 2019 Frézal et al.)

Telomeric DNA is composed of simple tandem repeat sequences and has a G-rich strand that runs 5' to 3' toward the chromosome terminus. Small RNAs with homology to telomeres have been observed in several organisms and could originate from telomeres or from interstitial telomere sequences (ITSs), which are composites of degenerate and perfect telomere repeat sequences found on chromosome arms. We identified Caenorhabditis elegans small RNAs composed of the Caenorhabditis telomere sequence (TTAGGC) n with up to three mismatches, which might interact with telomeres. We rigorously defined ITSs for genomes of C. elegans and for two closely related nematodes, Caenorhabditis briggsae and Caenorhabditis remanei Most telomeric small RNAs with mismatches originated from ITSs, which were depleted from mRNAs but were enriched in introns whose genes often displayed hallmarks of genomic silencing. C. elegans small RNAs composed of perfect telomere repeats were very rare but their levels increased by several orders of magnitude in C. briggsae and C. remanei Major small RNA species in C. elegans begin with a 5' guanine nucleotide, which was strongly depleted from perfect telomeric small RNAs of all three Caenorhabditis species. Perfect G-rich or C-rich telomeric small RNAs commonly began with 5' UAGGCU and 5' UUAGGC or 5' CUAAGC, respectively. In contrast, telomeric small RNAs with mismatches had a mixture of all four 5' nucleotides. We suggest that perfect telomeric small RNAs have a mechanism of biogenesis that is distinct from known classes of small RNAs and that a dramatic change in their regulation occurred during recent Caenorhabditis evolution., (© 2019 Frenk et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Long non-coding RNAs (lncRNAs) are a heterogeneous class of genes that do not code for proteins. Since lncRNAs (or a fraction thereof) are expected to be functional, many efforts have been dedicated to catalog lncRNAs in numerous organisms, but our knowledge of lncRNAs in non vertebrate species remains very limited. Here, we annotated lncRNAs using transcriptomic data from the same larval stage of four Caenorhabditis species. The number of annotated lncRNAs in self-fertile nematodes was lower than in out-crossing species. We used a combination of approaches to identify putatively homologous lncRNAs: synteny, sequence conservation, and structural conservation. We classified a total of 1,532 out of 7,635 genes from the four species into families of lncRNAs with conserved synteny and expression at the larval stage, suggesting that a large fraction of the predicted lncRNAs may be species specific. Despite both sequence and local secondary structure seem to be poorly conserved, sequences within families frequently shared BLASTn hits and short sequence motifs, which were more likely to be unpaired in the predicted structures. We provide the first multi-species catalog of lncRNAs in nematodes and identify groups of lncRNAs with conserved synteny and expression, that share exposed motifs.

Hybrid male sterility often evolves before female sterility or inviability of hybrids, implying that the accumulation of divergence between separated lineages should lead hybrid male sterility to have a more polygenic basis. However, experimental evidence is mixed. Here, we use the nematodes Caenorhabditis remanei and C. latens to characterize the underlying genetic basis of asymmetric hybrid male sterility and hybrid inviability. We demonstrate that hybrid male sterility is consistent with a simple genetic basis, involving a single X-autosome incompatibility. We also show that hybrid inviability involves more genomic compartments, involving diverse nuclear-nuclear incompatibilities, a mito-nuclear incompatibility, and maternal effects. These findings demonstrate that male sensitivity to genetic perturbation may be genetically simple compared to hybrid inviability in Caenorhabditis and motivates tests of generality for the genetic architecture of hybrid incompatibility across the breadth of phylogeny.

DNA recombination is required for effective segregation and diversification of genomes and for the successful completion of meiosis. Recent studies in various species hybrids have demonstrated a genetic link between DNA recombination and speciation. Consistent with this, we observed a striking suppression of recombination in the hybrids between two nematodes, the hermaphroditic Caenorhabditis briggsae and the gonochoristic C. nigoni. To unravel the molecular basis underlying the recombination suppression in their hybrids, we generated a C. nigoni genome with chromosome-level contiguity and produced an improved C. briggsae genome with resolved gaps up to 2.8 Mb. The genome alignment reveals not only high sequence divergences but also pervasive intra- and inter-chromosomal sequence re-arrangements between the two species, which are plausible culprits for the observed suppression. Comparison of recombination boundary sequences suggests that recombination in the hybrid requires extensive sequence homology, which is rarely seen between the two genomes. The new genomes and genomic libraries form invaluable resources for studying genome evolution, hybrid incompatibilities and sex evolution for this pair of model species.

To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsae C. nigoni 's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsae C. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short ( mss ) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss -null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors., (Copyright © 2018, American Association for the Advancement of Science.)

Normal vulval development in the nematode Caenorhabditis briggsae is identical to that in the related Caenorhabditis elegans. However, several experiments suggest that there are differences between the two species with respect to the contribution of EGF/Ras signaling. To investigate these differences genetically, we have characterized a C. briggsae mutant strain that phenocopies the effect observed when C. briggsae animals are treated with U0126, an inhibitor of the EGF pathway component MEK. We identify that the gene affected in the mutant strain is Cbr-sur-2, which encodes a MED23 mediator complex protein that acts downstream of EGF signaling in C. elegans and other organisms, such as mammals. When Cbr-sur-2 and Cel-sur-2 mutants are compared, we find that the production of additional vulval cells from P5.p and P7.p in C. elegans is dependent on proper development of P6.p, while C. briggsae does not have a similar requirement. Combined chemical and genetic interference with the EGF pathway completely eliminates vulval development in C. elegans but not in C. briggsae. Our results provide genetic evidence for the differing requirements for EGF signaling in the two species.

Limiting the debilitating consequences of ageing is a major medical challenge of our time. Robust pharmacological interventions that promote healthy ageing across diverse genetic backgrounds may engage conserved longevity pathways. Here we report results from the Caenorhabditis Intervention Testing Program in assessing longevity variation across 22 Caenorhabditis strains spanning 3 species, using multiple replicates collected across three independent laboratories. Reproducibility between test sites is high, whereas individual trial reproducibility is relatively low. Of ten pro-longevity chemicals tested, six significantly extend lifespan in at least one strain. Three reported dietary restriction mimetics are mainly effective across C. elegans strains, indicating species and strain-specific responses. In contrast, the amyloid dye ThioflavinT is both potent and robust across the strains. Our results highlight promising pharmacological leads and demonstrate the importance of assessing lifespans of discrete cohorts across repeat studies to capture biological variation in the search for reproducible ageing interventions.

To facilitate comparative genomic study in the Caenorhabditis species, the mitochondrial genome (mitogenome) of a nematode species Caenorhabditis nigoni (previous name: Caenorhabditis sp. 9) was generated using next-generation sequencing. The mitogenome length is 13,413 bp, containing 12 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs) and 2 non-coding regions (NCR). The genome organization and nucleotide composition is very similar to that of the mitogenome of C. elegans and C. briggsae. Mitogenome of C. nigoni shows higher sequence similarity to C. briggsae than to C. elegans, which is consistent with the fact that C. nigoni is a sister species of C. briggsae. However, as in C. elegans, two NCRs present in the mitogenome of C. briggsae are missing in C. nigoni. The mitogenome sequence of C. nigoni plays an important role in further studies of phylogenetics, population genetics and evolutionary genetics in nematode species.

Heparan sulfates (HS) are glycosaminoglycans of the extracellular matrices and characterized by complex modification patterns owing to sulfations, epimerization, and acetylation. Distinct HS modification patterns have been shown to modulate protein-protein interactions during development in general and of the nervous system in particular. This has led to the heparan sulfate code hypothesis, which posits that specifically modified HS epitopes are distributed in a tissue and cell-specific fashion to orchestrate neural circuit formation. Whether an HS code exists in vivo, how specific or how evolutionarily conserved the anatomical distribution of an HS code may be has remained unknown. Here we conduct a systematic comparison of HS modification patterns in the nematode Caenorhabditis elegans using transgenic expression of 33 different HS-specific single chain variable fragment antibodies. We find that some HS modification patterns are widely distributed in the nervous system. In contrast, other HS modification patterns appear highly cell-specific in both non-neuronal and neuronal cells. Some patterns can be as restricted in their localization as to single neurites or synaptic connections between two neurons. This restricted anatomical localization of specific HS patterns can be evolutionarily conserved over a span of 80-100 million years in the divergent nematode species Caenorhabditis briggsae suggesting structural and, possibly functional conservation of glycosaminoglycan structures similar to proteins. These findings suggest a HS code with subcellularly localized, unique glycan identities in the nervous system., (© The Author 2016. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: The high density of tandem repeat sequences (satellites) in nematode genomes and the availability of genome sequences from several species in the group offer a unique opportunity to better understand the evolutionary dynamics and the functional role of these sequences. We take advantage of the previously developed SATFIND program to study the satellites in four Caenorhabditis species and investigate these questions., Methods: The identification and comparison of satellites is carried out in three steps. First we find all the satellites present in each species with the SATFIND program. Each satellite is defined by its length, number of repeats, and repeat sequence. Only satellites with at least ten repeats are considered. In the second step we build satellite families with a newly developed alignment program. Satellite families are defined by a consensus sequence and the number of satellites in the family. Finally we compare the consensus sequence of satellite families in different species., Results: We give a catalog of individual satellites in each species. We have also identified satellite families with a related sequence and compare them in different species. We analyze the turnover of satellites: they increased in size through duplications of fragments of 100-300 bases. It appears that in many cases they have undergone an explosive expansion. In C. elegans we have identified a subset of large satellites that have strong affinity for the centromere protein CENP-A. We have also compared our results with those obtained from other species, including one nematode and three mammals., Conclusions: Most satellite families found in Caenorhabditis are species-specific; in particular those with long repeats. A subset of these satellites may facilitate the formation of kinetochores in mitosis. Other satellite families in C. elegans are either related to Helitron transposons or to meiotic pairing centers.

Background: How do very small animals with limited long-distance dispersal abilities move between locations, especially if they prefer ephemeral micro-habitats that are only available for short periods of time? The free-living model nematode Caenorhabditis elegans and several congeneric taxa appear to be common in such short-lived environments, for example decomposing fruits or other rotting plant material. Dispersal is usually assumed to depend on animal vectors, yet all current data is based on only a limited number of studies. In our project we performed three comprehensive field surveys on possible invertebrate vectors in North German locations containing populations of C. elegans and two related species, especially C. remanei, and combined these screens with an experimental analysis of persistence in one of the vector taxa., Results: Our field survey revealed that Caenorhabditis nematodes are commonly found in slugs, isopods, and chilopods, but are not present in the remaining taxonomic groups examined. Surprisingly, the nematodes were frequently isolated from the intestines of slugs, even if slugs were not collected in close association with suitable substrates for Caenorhabditis proliferation. This suggests that the nematodes are able to enter the slug intestines and persist for certain periods of time. Our experimental analysis confirmed the ability of C. elegans to invade slug intestines and subsequently be excreted alive with the slug feces, although only for short time periods under laboratory conditions., Conclusions: We conclude that three invertebrate taxonomic groups represent potential vectors of Caenorhabditis nematodes. The nematodes appear to have evolved specific adaptations to enter and persist in the harsh environment of slug intestines, possibly indicating first steps towards a parasitic life-style.

Evolutionary interactions across levels of biological organization contribute to a variety of fundamental processes including genome evolution, reproductive mode transitions, species diversification, and extinction. Evolutionary theory predicts that so-called "selfish" genetic elements will proliferate when the host effective population size (Ne) is small, but direct tests of this prediction remain few. We analyzed the evolutionary dynamics of deletion-containing mitochondrial DNA (ΔmtDNA) molecules, previously characterized as selfish elements, in six different natural strains of the nematode Caenorhabditis briggsae allowed to undergo experimental evolution in a range of population sizes (N = 1, 10, 100, and 1,000) for a maximum of 50 generations. Mitochondrial DNA (mtDNA) was analyzed for replicate lineages at each five-generation time point. Ten different ΔmtDNA molecule types were observed and characterized across generations in the experimental populations. Consistent with predictions from evolutionary theory, lab lines evolved in small-population sizes (e.g., nematode N = 1) were more susceptible to accumulation of high levels of preexisting ΔmtDNA compared with those evolved in larger populations. New ΔmtDNA elements were observed to increase in frequency and persist across time points, but almost exclusively at small population sizes. In some cases, ΔmtDNA levels decreased across generations when population size was large (nematode N = 1,000). Different natural strains of C. briggsae varied in their susceptibilities to ΔmtDNA accumulation, owing in part to preexisting compensatory mtDNA alleles in some strains that prevent deletion formation. This analysis directly demonstrates that the evolutionary trajectories of ΔmtDNA elements depend upon the population-genetic environments and molecular-genetic features of their hosts., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The roundworm Caenorhabditis elegans has risen to the status of a top model organism for biological research in the last fifty years. Among laboratory animals, this tiny nematode is one of the simplest and easiest organisms to handle. And its life outside the laboratory is beginning to be unveiled. Like other model organisms, C. elegans has a boom-and-bust lifestyle. It feasts on ephemeral bacterial blooms in decomposing fruits and stems. After resource depletion, its young larvae enter a migratory diapause stage, called the dauer. Organisms known to be associated with C. elegans include migration vectors (such as snails, slugs and isopods) and pathogens (such as microsporidia, fungi, bacteria and viruses). By deepening our understanding of the natural history of C. elegans, we establish a broader context and improved tools for studying its biology.

Studies of vulval development in the nematode C. elegans have identified many genes that are involved in cell division and differentiation processes. Some of these encode components of conserved signal transduction pathways mediated by EGF, Notch, and Wnt. To understand how developmental mechanisms change during evolution, we are doing a comparative analysis of vulva formation in C. briggsae, a species that is closely related to C. elegans. Here, we report 14 mutations in 7 Multivulva (Muv) genes in C. briggsae that inhibit inappropriate division of vulval precursors. We have developed a new efficient and cost-effective gene mapping method to localize Muv mutations to small genetic intervals on chromosomes, thus facilitating cloning and functional studies. We demonstrate the utility of our method by determining molecular identities of three of the Muv genes that include orthologs of Cel-lin-1 (ETS) and Cel-lin-31 (Winged-Helix) of the EGF-Ras pathway and Cel-pry-1 (Axin), of the Wnt pathway. The remaining four genes reside in regions that lack orthologs of known C. elegans Muv genes. Inhibitor studies demonstrate that the Muv phenotype of all four new genes is dependent on the activity of the EGF pathway kinase, MEK. One of these, Cbr-lin(gu167), shows modest increase in the expression of Cbr-lin-3/EGF compared to wild type. These results argue that while Cbr-lin(gu167) may act upstream of Cbr-lin-3/EGF, the other three genes influence the EGF pathway downstream or in parallel to Cbr-lin-3. Overall, our findings demonstrate that the genetic program underlying a conserved developmental process includes both conserved and divergent functional contributions., (© 2014 Wiley Periodicals, Inc.)

Determining the causes and evolution of reproductive barriers to gene flow between populations, speciation, is the key to understanding the origin of diversity in nature. Many species manifest hybrid breakdown when they intercross, characterized by increasingly exacerbated problems in later generations of hybrids. Recently, Caenorhabditis nematodes have emerged as a genetic model for studying speciation, and here we investigate the nature and causes of hybrid breakdown between Caenorhabditis remanei and C. latens. We quantify partial F1 hybrid inviability and extensive F2 hybrid inviability; the ~75% F2 embryonic arrest occurs primarily during gastrulation or embryonic elongation. Moreover, F1 hybrid males exhibit Haldane's rule asymmetrically for both sterility and inviability, being strongest when C. remanei serves as maternal parent. We show that the mechanism by which sterile hybrid males are incapable of transferring sperm or a copulatory plug involves defective gonad morphogenesis, which we hypothesize results from linker cell defects in migration and/or cell death during development. This first documented case of partial hybrid male sterility in Caenorhabditis follows expectations of Darwin's corollary to Haldane's rule for asymmetric male fitness, providing a powerful foundation for molecular dissection of intrinsic reproductive barriers and divergence of genetic pathways controlling organ morphogenesis., (© 2014 Wiley Periodicals, Inc.)

The major families of chromatin remodelers have been conserved throughout eukaryotic evolution. Because they play broad, pleiotropic roles in gene regulation, it was not known if their functions could change rapidly. Here, we show that major alterations in the use of chromatin remodelers are possible, because the nucleosome remodeling factor (NURF) complex has acquired a unique role in the sperm/oocyte decision of the nematode Caenorhabditis briggsae. First, lowering the activity of C. briggsae NURF-1 or ISW-1, the core components of the NURF complex, causes germ cells to become oocytes rather than sperm. This observation is based on the analysis of weak alleles and null mutations that were induced with TALENs and on RNA interference. Second, qRT-polymerase chain reaction data show that the C. briggsae NURF complex promotes the expression of Cbr-fog-1 and Cbr-fog-3, two genes that control the sperm/oocyte decision. This regulation occurs in the third larval stage and affects the expression of later spermatogenesis genes. Third, double mutants reveal that the NURF complex and the transcription factor TRA-1 act independently on Cbr-fog-1 and Cbr-fog-3. TRA-1 binds both promoters, and computer analyses predict that these binding sites are buried in nucleosomes, so we suggest that the NURF complex alters chromatin structure to allow TRA-1 access to Cbr-fog-1 and Cbr-fog-3. Finally, lowering NURF activity by mutation or RNA interference does not affect this trait in other nematodes, including the sister species C. nigoni, so it must have evolved recently. We conclude that altered chromatin remodeling could play an important role in evolutionary change., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Understanding the robustness of developmental systems requires insights into the sensitivity of underlying molecular and cellular parameters to perturbations, and how such sensitivity evolves. We address these issues using vulval cell fate determination--a reproducible and robust patterning process regulated by a cross-talk of EGF-Ras-MAPK and Delta-Notch pathways. Although the final vulval cell fate pattern is identical in all Caenorhabditis species, the patterning process underlies extensive cryptic genetic variation between and within species. Here, we tested whether this cryptic genetic variation translates into variation in developmental sensitivity to environmental perturbations. We disrupted vulval patterning using thermal perturbations to quantify and compare environmental sensitivity of different system parameters between distinct genotypes of C. elegans and C. briggsae. Thermal perturbations globally debuffered vulval development, triggering diverse pattering variants, whose frequency and spectra were strongly species- and genotype-dependent. This condition-dependent variation indicates that environmental sensitivity of different system properties, such as vulval competence or vulval induction, is subject to evolutionary change. High temperature induced a genotype-specific decrease of secondary fate induction and corresponding Notch pathway activity in the C. elegans N2 strain; in contrast, hypoinduction of the primary cell fate was never observed. Vulval precursor cells therefore differ in temperature sensitivity and such cell-specific sensitivity shows evolutionary variation. We further compared spectra of temperature-induced vulval variants to the ones induced by mutation accumulation in the same genotypes. In response to either perturbation, we observed similar genotype-dependence of variant production, allowing identification of distinct system features most sensitive to both mutation and environment. Taken together, we show how sensitivity of system parameters regulating Caenorhabditis vulval development depends on subtle interactions between perturbations and genetic background. Our results imply that cryptic genetic variation may reflect evolutionary variation in developmental robustness, therefore potentially contributing to the maintenance of phenotypic precision when facing perturbations., (© 2014 Wiley Periodicals, Inc.)

The rapid pace of species discovery outstrips the rate of species description in many taxa. This problem is especially acute for Caenorhabditis nematodes, where the naming of distinct species would greatly improve their visibility and usage for biological research, given the thousands of scientists studying Caenorhabditis. Species description and naming has been hampered in Caenorhabditis, in part due to the presence of morphologically cryptic species despite complete biological reproductive isolation and often enormous molecular divergence. With the aim of expediting species designations, here we propose and apply a revised framework for species diagnosis and description in this group. Our solution prioritizes reproductive isolation over traditional morphological characters as the key feature in delineating and diagnosing new species, reflecting both practical considerations and conceptual justifications. DNA sequence divergence criteria help prioritize crosses for establishing patterns of reproductive isolation among the many species of Caenorhabditis known to science, such as with the ribosomal internal transcribed spacer-2 (ITS2) DNA barcode. By adopting this approach, we provide new species name designations for 15 distinct biological species, thus increasing the number of named Caenorhabditis species in laboratory culture by nearly 3-fold. We anticipate that the improved accessibility of these species to the research community will expand the opportunities for study and accelerate our understanding of diverse biological phenomena.

In this study, a modified miRsearch program was developed in C++ for the detection of miRNAs. All the mature miRNA sequences of Caenorhabditis elegans, Caenorhabditis briggsae, Caenorhabditis remanei, Drosophila melanogaster, Homo sapiens and Rattus norvegicus available in miRbase was searched by this program for homologous sequences with a maximum of 3-mismatches in the chromosomes of C. elegans excluding the miRNAs of C. elegans. The same strategy was repeated for C. briggsae excluding C. briggsae miRNAs. The probable pre-miRNA sequences with stem loop secondary structures were assessed by implementation of longest-common-subsequence (LCS) algorithm with appropriate scoring system. As miRNA genes could be on either strand, each sequence was searched in both forward and reverse strands. The putative miRNAs were viewed through Mapviewer to identify their intronic or intergenic location in C. elegans and C. briggsae genomes. Further, the quality of stem-loop formation of the remaining pre-miRNA sequences was assessed through RNAFOLD. This algorithm will be helpful in detection of potential miRNAs in future sequencing data, making this an invaluable tool for miRNA prediction.

Background: Gene duplicates often exhibit asymmetric rates of molecular evolution in their early evolutionary existence. This asymmetry in rates is thought to signify the maintenance of the ancestral function by one copy and the removal of functional constraint on the other copy, enabling it to embark on a novel evolutionary trajectory. Here I focused on a large population of evolutionarily young gene duplicates (KS ≤ 0.14) in the Caenorhabditis elegans genome in order to conduct the first combined analysis of four predictors (evolutionary age, chromosomal location, structural resemblance between duplicates, and duplication span) which may be implicated in the asymmetric sequence divergence of paralogs at the nucleotide and amino acid level. In addition, I investigate if either paralog is equally likely to embark on a trajectory of accelerated sequence evolution or whether the derived paralog is more likely to exhibit faster sequence evolution., Results: Three predictors (evolutionary age of duplicates, chromosomal location and duplication span) serve as major determinants of sequence asymmetry between C. elegans paralogs. Paralogs diverge asymmetrically in sequence with increasing evolutionary age, the relocation of one copy to a different chromosome and attenuated duplication spans that likely fail to capture the entire ancestral repertoire of coding sequence and regulatory elements. Furthermore, for paralogs residing on the same chromosome, opposite transcriptional orientation and increased genomic distance do not increase sequence asymmetry between paralogs. For a subset of duplicate pairs wherein the ancestral versus derived paralog could be distinguished, the derived paralogs are more likely to evolve at accelerated rates., Conclusions: This genome-wide study of evolutionarily young duplicates stemming primarily from DNA-mediated small-scale duplication events demonstrates that genomic relocation to a new chromosome has important consequences for asymmetric divergence of paralogs, akin to paralogs arising from RNA-mediated duplication events. Additionally, the duplication span is negatively correlated with sequence rate asymmetry among paralogs, suggesting that attenuated duplication spans stemming from incomplete duplication of the ORF and/or ancestral regulatory elements further accelerate sequence divergence between paralogs. Cumulatively, derived copies exhibit accelerated rates of sequence evolution suggesting that they are primed for a divergent evolutionary trajectory by changes in structure and genomic context at inception.

The nematode Caenorhabditis elegans contains each of the broad classes of eukaryotic small RNAs, including microRNAs (miRNAs), endogenous small-interfering RNAs (endo-siRNAs), and piwi-interacting RNAs (piRNAs). To better understand the evolution of these regulatory RNAs, we deep-sequenced small RNAs from C. elegans and three closely related nematodes: C. briggsae, C. remanei, and C. brenneri. The results reveal a fluid landscape of small RNA pathways with essentially no conservation of individual sequences aside from a subset of miRNAs. We identified 54 miRNA families that are conserved in each of the four species, as well as numerous miRNAs that are species-specific or shared between only two or three species. Despite a lack of conservation of individual piRNAs and siRNAs, many of the features of each pathway are conserved between the different species. We show that the genomic distribution of 26G siRNAs and the tendency for piRNAs to cluster is conserved between C. briggsae and C. elegans. We also show that, in each species, 26G siRNAs trigger stage-specific secondary siRNA formation. piRNAs in each species also trigger secondary siRNA formation from targets containing up to three mismatches. Finally, we show that the production of male- and female-specific piRNAs is conserved in all four species, suggesting distinct roles for piRNAs in male and female germlines.

Background: In stark contrast to the wealth of detail about C. elegans developmental biology and molecular genetics, biologists lack basic data for understanding the abundance and distribution of Caenorhabditis species in natural areas that are unperturbed by human influence., Methods: Here we report the analysis of dense sampling from a small, remote site in the Amazonian rain forest of the Nouragues Natural Reserve in French Guiana., Results: Sampling of rotting fruits and flowers revealed proliferating populations of Caenorhabditis, with up to three different species co-occurring within a single substrate sample, indicating remarkable overlap of local microhabitats. We isolated six species, representing the highest local species richness for Caenorhabditis encountered to date, including both tropically cosmopolitan and geographically restricted species not previously isolated elsewhere. We also documented the structure of within-species molecular diversity at multiple spatial scales, focusing on 57 C. briggsae isolates from French Guiana. Two distinct genetic subgroups co-occur even within a single fruit. However, the structure of C. briggsae population genetic diversity in French Guiana does not result from strong local patterning but instead presents a microcosm of global patterns of differentiation. We further integrate our observations with new data from nearly 50 additional recently collected C. briggsae isolates from both tropical and temperate regions of the world to re-evaluate local and global patterns of intraspecific diversity, providing the most comprehensive analysis to date for C. briggsae population structure across multiple spatial scales., Conclusions: The abundance and species richness of Caenorhabditis nematodes is high in a Neotropical rainforest habitat that is subject to minimal human interference. Microhabitat preferences overlap for different local species, although global distributions include both cosmopolitan and geographically restricted groups. Local samples for the cosmopolitan C. briggsae mirror its pan-tropical patterns of intraspecific polymorphism. It remains an important challenge to decipher what drives Caenorhabditis distributions and diversity within and between species.

Mating system transitions dramatically alter the evolutionary trajectories of genomes that can be revealed by contrasts of species with disparate modes of reproduction. For such transitions in Caenorhabditis nematodes, some major causes of genome variation in selfing species have been discerned. And yet, we have only limited understanding of species-wide population genetic processes for their outcrossing relatives, which represent the reproductive state of the progenitors of selfing species. Multilocus-multipopulation sequence polymorphism data provide a powerful means to uncover the historical demography and evolutionary processes that shape genomes. Here we survey nucleotide polymorphism across the X chromosome for three populations of the outcrossing nematode Caenorhabditis remanei and demonstrate its divergence from a fourth population describing a closely related new species from China, C. sp. 23. We find high genetic variation globally and within each local population sample. Despite geographic barriers and moderate genetic differentiation between Europe and North America, considerable gene flow connects C. remanei populations. We discovered C. sp. 23 while investigating C. remanei, observing strong genetic differentiation characteristic of reproductive isolation that was confirmed by substantial F2 hybrid breakdown in interspecific crosses. That C. sp. 23 represents a distinct biological species provides a cautionary example of how standard practice can fail for mating tests of species identity in this group. This species pair permits full application of divergence population genetic methods to obligately outcrossing species of Caenorhabditis and also presents a new focus for interrogation of the genetics and evolution of speciation with the Caenorhabditis model system.

Molecular hyperdiversity has been documented in viruses, prokaryotes and eukaryotes. Such organisms undermine the assumptions of the infinite-sites mutational model, because multiple mutational events at a site comprise a non-negligible portion of polymorphisms. Moreover, different sampling schemes of individuals from species with subdivided populations can profoundly influence resulting patterns and interpretations of molecular variation. Inspired by molecular hyperdiversity in the nematode Caenorhabditis sp. 5, which exhibits average pairwise differences among synonymous sites of >5% as well as modest population structure, we investigated via coalescent simulation the joint effects of a finite-sites mutation (FSM) process and population subdivision on the variant frequency spectrum. From many demes interconnected through a stepping-stone migration model, we constructed local samples from a single deme, pooled samples from several demes and scattered samples of a single individual from numerous demes. Compared with a single panmictic population at equilibrium, we find that high population mutation rates induce a deficit of rare variants (positive Tajima's D) under a FSM model. Population structure also induces such a skew for local samples when migration is high and for pooled samples when migration is low. Contrasts of sampling schemes for C. sp. 5 imply high mutational input coupled with high migration. We propose that joint analysis of local, pooled and scattered samples for species with subdivided populations provides a means of improving inference of demographic history, by virtue of the partially distinct patterns of polymorphism that manifest when sequences are analyzed according to differing sampling schemes., (© 2012 Blackwell Publishing Ltd.)

It is important to understand the forces that shape the size and evolutionary histories of gene families. Here, we investigated the evolution of non-protein-coding RNA genes in the genomes of Caenorhabditis nematodes. We specifically focused on nested arrangements, that is, cases in which an RNA gene is entirely contained in an intron of another gene. Comparing these arrangements between species simplifies the inference of orthology and, therefore, of evolutionary fates of nested genes. Two distinct patterns are evident in the data. Genes encoding small nuclear RNAs (snRNAs) and transfer RNAs form large families, which have persisted since before the common ancestor of Metazoa. Yet, individual genes die relatively rapidly, with few orthologs having survived since the divergence of Caenorhabditis elegans and Caenorhabditis briggsae. In contrast, genes encoding small nucleolar RNAs (snoRNAs) are either single-copy or form small families. Individual snoRNAs turn over at a relatively slow rate-most C. elegans genes have clearly identifiable orthologs in C. briggsae. We also found that in Drosophila, genes from larger snRNA families die at a faster rate than their counterparts from single-gene families. These results suggest that a relationship between family size and the rate of gene turnover may be a general feature of genome evolution.

Pleiotropic developmental regulators have been repeatedly linked to the evolution of anatomical novelties. Known mechanisms include cis-regulatory DNA changes that alter regulator transcription patterns or modify target-gene linkages. Here, we examine the role of another form of regulation, translational control, in the repeated evolution of self-fertile hermaphroditism in Caenorhabditis nematodes. Caenorhabditis elegans hermaphrodites initiate spermatogenesis in an otherwise female body through translational repression of the gene tra-2. This repression is mediated by GLD-1, an RNA-binding protein also required for oocyte meiosis and differentiation. By contrast, we show that in the convergently hermaphroditic Caenorhabditis briggsae, GLD-1 acts to promote oogenesis. The opposite functions of gld-1 in these species are not gene-intrinsic, but instead result from the unique contexts for its action that evolved in each. In C. elegans, GLD-1 became essential for promoting XX spermatogenesis via changes in the tra-2 mRNA and evolution of the species-specific protein FOG-2. C. briggsae GLD-1 became an essential repressor of sperm-promoting genes, including Cbr-puf-8, and did not evolve a strong association with tra-2. Despite its variable roles in sex determination, the function of gld-1 in female meiotic progression is ancient and conserved. This conserved role may explain why gld-1 is repeatedly recruited to regulate hermaphroditism. We conclude that, as with transcription factors, spatially localized translational regulators play important roles in the evolution of anatomical novelties.

Background: The nematode Caenorhabditis elegans is a major laboratory model in biology. Only ten Caenorhabditis species were available in culture at the onset of this study. Many of them, like C. elegans, were mostly isolated from artificial compost heaps, and their more natural habitat was unknown., Results: Caenorhabditis nematodes were found to be proliferating in rotten fruits, flowers and stems. By collecting a large worldwide set of such samples, 16 new Caenorhabditis species were discovered. We performed mating tests to establish biological species status and found some instances of semi-fertile or sterile hybrid progeny. We established barcodes for all species using ITS2 rDNA sequences. By obtaining sequence data for two rRNA and nine protein-coding genes, we determined the likely phylogenetic relationships among the 26 species in culture. The new species are part of two well-resolved sister clades that we call the Elegans super-group and the Drosophilae super-group. We further scored phenotypic characters such as reproductive mode, mating behavior and male tail morphology, and discuss their congruence with the phylogeny. A small space between rays 2 and 3 evolved once in the stem species of the Elegans super-group; a narrow fan and spiral copulation evolved once in the stem species of C. angaria, C. sp. 8 and C. sp. 12. Several other character changes occurred convergently. For example, hermaphroditism evolved three times independently in C. elegans, C. briggsae and C. sp. 11. Several species can co-occur in the same location or even the same fruit. At the global level, some species have a cosmopolitan distribution: C. briggsae is particularly widespread, while C. elegans and C. remanei are found mostly or exclusively in temperate regions, and C. brenneri and C. sp. 11 exclusively in tropical zones. Other species have limited distributions, for example C. sp. 5 appears to be restricted to China, C. sp. 7 to West Africa and C. sp. 8 to the Eastern United States., Conclusions: Caenorhabditis are "fruit worms", not soil nematodes. The 16 new species provide a resource and their phylogeny offers a framework for further studies into the evolution of genomic and phenotypic characters.

Self-fertile hermaphrodites have evolved independently several times in the genus Caenorhabditis [1, 2]. These XX hermaphrodites make smaller sperm than males [3, 4], which they use to fertilize their own oocytes. Because larger sperm outcompete smaller sperm in nematodes [3-5], it had been assumed that this dimorphism evolved in response to sperm competition. However, we show that it was instead caused by a developmental bias. When we transformed females of the species Caenorhabditis remanei into hermaphrodites [6], their sperm were significantly smaller than those of males. Because this species never makes hermaphrodites in the wild, this dimorphism cannot be due to selection. Instead, analyses of the related nematode Caenorhabditis elegans suggest that this dimorphism might reflect the development of sperm within the distinct physiological environment of the hermaphrodite gonad. These results reveal a new mechanism for some types of developmental bias-the effects of a novel physical location alter the development of ectopic cells in predictable ways., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Background: Hybridization of heterologous (non-specific) nucleic acids onto arrays designed for model-organisms has been proposed as a viable genomic resource for estimating sequence variation and gene expression in non-model organisms. However, conventional methods of normalization that assume equivalent distributions (such as quantile normalization) are inappropriate when applied to non-specific (heterologous) hybridization. We propose an algorithm for normalizing and centering intensity data from heterologous hybridization that makes no prior assumptions of distribution, reduces the false appearance of homology, and provides a way for researchers to confirm whether heterologous hybridization is suitable., Results: Data are normalized by adjusting for Gibbs free energy binding, and centered by adjusting for the median of a common set of control probes assumed to be equivalently dissimilar for all species. This procedure was compared to existing approaches and found to be as successful as Loess normalization at detecting sequence variations (deletions) and even more successful than quantile normalization at reducing the accumulation of false positive probe matches between two related nematode species, Caenorhabditis elegans and C. briggsae. Despite the improvements, we still found that probe fluorescence intensity was too poorly correlated with sequence similarity to result in reliable detection of matching probe sequence., Conclusions: Cross-species hybridizations can be a way to adapt genome-enabled tools for closely related non-model organisms, but data must be appropriately normalized and centered in a way that accommodates hybridization of nucleic acids with diverged sequence. For short, 25-mer probes, hybridization intensity alone may be insufficiently correlated with sequence similarity to allow reliable inference of homology at the probe level.

Background: Mutations that impair mitochondrial functioning are associated with a variety of metabolic and age-related disorders. A barrier to rigorous tests of the role of mitochondrial dysfunction in aging processes has been the lack of model systems with relevant, naturally occurring mitochondrial genetic variation. Toward the goal of developing such a model system, we studied natural variation in life history, metabolic, and aging phenotypes as it relates to levels of a naturally-occurring heteroplasmic mitochondrial ND5 deletion recently discovered to segregate among wild populations of the soil nematode, Caenorhabditis briggsae. The normal product of ND5 is a central component of the mitochondrial electron transport chain and integral to cellular energy metabolism., Results: We quantified significant variation among C. briggsae isolates for all phenotypes measured, only some of which was statistically associated with isolate-specific ND5 deletion frequency. We found that fecundity-related traits and pharyngeal pumping rate were strongly inversely related to ND5 deletion level and that C. briggsae isolates with high ND5 deletion levels experienced a tradeoff between early fecundity and lifespan. Conversely, oxidative stress resistance was only weakly associated with ND5 deletion level while ATP content was unrelated to deletion level. Finally, mean levels of reactive oxygen species measured in vivo showed a significant non-linear relationship with ND5 deletion level, a pattern that may be driven by among-isolate variation in antioxidant or other compensatory mechanisms., Conclusions: Our findings suggest that the ND5 deletion may adversely affect fitness and mitochondrial functioning while promoting aging in natural populations, and help to further establish this species as a useful model for explicit tests of hypotheses in aging biology and mitochondrial genetics.

Caenorhabditis elegans is a premier model organism upon which considerable knowledge of basic cell and developmental biology has been built. Yet, as is true for many traditional model systems, we have limited knowledge of the ecological context in which these systems evolved, severely limiting our understanding of gene function. A better grasp of the ecology of model systems would help us immensely in understanding the functionality of genes and evolution of genomes in an environmental context. Consequently, there are ongoing efforts to uncover natural populations of this model system globally. Here, we describe the discovery of a Caenorhabditis briggsae strain and its bacterial associate (Serratia sp.) that form an entomopathogenic complex in the wild. Laboratory experiments confirm that this nematode and its natural bacterial associate can penetrate, kill and reproduce in an insect host and that the bacterial associate can induce this insect pathogenic life cycle in other Caenorhabditis species, including C. elegans. Our findings suggest that this life history may be widespread in nature and critical to the understanding of the biology of this important model organism. Caenorhabditis-insect interaction could be a key factor in our quest for a better grasp of gene functionality in this important model species. The discovered association, consequently, would provide an ecological framework for functional genomics of Caenorhabditis.

We have developed a nematode transformation vector carrying the bacterial neomycin resistance gene (NeoR) and shown that it could confer resistance to G-418 on both wild-type Caenorhabditis elegans and C. briggsae. This selection system allows hands-off maintenance and enrichment of transgenic worms carrying non-integrated transgenes on selective plates. We also show that this marker can be used for Mos1-mediated single-copy insertion in wild-type genetic backgrounds (MosSCI-biotic).

Although the nematode Caenorhabditis elegans produces self-fertile hermaphrodites, it descended from a male/female species, so hermaphroditism provides a model for the origin of novel traits. In the related species C. remanei, which has only male and female sexes, lowering the activity of tra-2 by RNA interference created XX animals that made spermatids as well as oocytes, but their spermatids could not activate without the addition of male seminal fluid. However, by lowering the expression of both tra-2 and swm-1, a gene that regulates sperm activation in C. elegans, we produced XX animals with active sperm that were self-fertile. Thus, the evolution of hermaphroditism in Caenorhabditis probably required two steps: a mutation in the sex-determination pathway that caused XX spermatogenesis and a mutation that allowed these spermatids to self-activate.

We present a highly accurate gene-prediction system for eukaryotic genomes, called mGene. It combines in an unprecedented manner the flexibility of generalized hidden Markov models (gHMMs) with the predictive power of modern machine learning methods, such as Support Vector Machines (SVMs). Its excellent performance was proved in an objective competition based on the genome of the nematode Caenorhabditis elegans. Considering the average of sensitivity and specificity, the developmental version of mGene exhibited the best prediction performance on nucleotide, exon, and transcript level for ab initio and multiple-genome gene-prediction tasks. The fully developed version shows superior performance in 10 out of 12 evaluation criteria compared with the other participating gene finders, including Fgenesh++ and Augustus. An in-depth analysis of mGene's genome-wide predictions revealed that approximately 2200 predicted genes were not contained in the current genome annotation. Testing a subset of 57 of these genes by RT-PCR and sequencing, we confirmed expression for 24 (42%) of them. mGene missed 300 annotated genes, out of which 205 were unconfirmed. RT-PCR testing of 24 of these genes resulted in a success rate of merely 8%. These findings suggest that even the gene catalog of a well-studied organism such as C. elegans can be substantially improved by mGene's predictions. We also provide gene predictions for the four nematodes C. briggsae, C. brenneri, C. japonica, and C. remanei. Comparing the resulting proteomes among these organisms and to the known protein universe, we identified many species-specific gene inventions. In a quality assessment of several available annotations for these genomes, we find that mGene's predictions are most accurate.

Background: MicroRNAs (miRNAs) are small, noncoding RNA molecules that act as post-transcriptional regulators of gene expression. Studies concerning transcriptional regulation of miRNAs have so far concentrated on those located within the intergenic region of the genome and the search for putative promoters, thus leaving open the question of the existence of possible regulatory elements common to all miRNAs including those located in introns of protein coding genes., Results: In this study, we initially searched for motifs occurring in the area 1000 bp upstream from all miRNAs independent of their genomic location. We discovered a previously unknown sequence motif GANNNNGA that displayed a conserved distribution in the nematode worms Caenorhabditis elegans and Caenorhabditis briggsae. This motif had a peak occurrence at 500 bp upstream, with a sharp drop-off toward the miRNA start site. Further analysis indicated that this motif was locally restricted and not enriched 1000-5000 bp upstream or 0-2000 bp downstream of the miRNA start site. In addition, this motif was observed to be most abundant in the upstream sequences of two important miRNAs, mir-1 and mir-124. This abundance was also conserved in phylogenetically distant species including human and mouse., Conclusion: The results show that the motif GANNNNGA is conserved close to miRNA precursor start sites, suggesting that it may be involved in miRNA sequence recognition or regulation. This data provides important knowledge for the identification and computational prediction of miRNA sequences.

Background: The phylum Nematoda is biologically diverse, including parasites of plants and animals as well as free-living taxa. Underpinning this diversity will be commensurate diversity in expressed genes, including gene sets associated specifically with evolution of parasitism., Methods and Findings: Here we have analyzed the extensive expressed sequence tag data (available for 37 nematode species, most of which are parasites) and define over 120,000 distinct putative genes from which we have derived robust protein translations. Combined with the complete proteomes of Caenorhabditis elegans and Caenorhabditis briggsae, these proteins have been grouped into 65,000 protein families that in turn contain 40,000 distinct protein domains. We have mapped the occurrence of domains and families across the Nematoda and compared the nematode data to that available for other phyla. Gene loss is common, and in particular we identify nearly 5,000 genes that may have been lost from the lineage leading to the model nematode C. elegans. We find a preponderance of novelty, including 56,000 nematode-restricted protein families and 26,000 nematode-restricted domains. Mapping of the latest time-of-origin of these new families and domains across the nematode phylogeny revealed ongoing evolution of novelty. A number of genes from parasitic species had signatures of horizontal transfer from their host organisms, and parasitic species had a greater proportion of novel, secreted proteins than did free-living ones., Conclusions: These classes of genes may underpin parasitic phenotypes, and thus may be targets for development of effective control measures.

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.

Dietary restriction is known to promote longevity in a variety of eukaryotic organisms. Most studies of dietary restriction have been performed on animals bred for many generations under conditions that differ substantially from their natural environment, raising the possibility that some apparent beneficial effects of dietary restriction are due to adaptation to laboratory conditions. To address this question in an invertebrate model, we determined the effect of dietary restriction by bacterial deprivation on life span in five different wild-derived Caenorhabditis elegans strains and two strains of the related species Caenorhabditis remanei. Longevity was enhanced in each of the wild-derived C. elegans strains, in most cases to a degree similar to that observed in N2, the standard laboratory strain. Both strains of C. remanei were substantially longer lived any of the C. elegans isolates, produced larger brood sizes, and retained the ability to produce offspring for a longer period of time. Dietary restriction failed to increase mean life span in one C. remanei isolate, but significantly increased the maximum life span of both C. remanei strains. Thus, we find no evidence that adaptation to laboratory conditions has significantly altered the aging process in C. elegans under either standard or food-restricted conditions.

We used network-level conservation between pairs of fly (Drosophila melanogaster/D. pseudoobscura) and worm (Caenorhabditis elegans/C. briggsae) genomes to detect highly conserved mRNA motifs in 3' untranslated regions. Many of these elements are complementary to the 5' extremity of known microRNAs (miRNAs), and likely correspond to their target sites. We also identify known targets of RNA-binding proteins, and many novel sites not yet known to be functional. Coherent sets of genes with similar function often bear the same conserved elements, providing new insights into their cellular functions. We also show that target sites for distinct miRNAs are often simultaneously conserved, suggesting combinatorial regulation by multiple miRNAs. A genome-wide search for conserved stem-loops, containing complementary sequences to the novel sites, revealed many new candidate miRNAs that likely target them. We also provide evidence that posttranscriptional networks have undergone extensive rewiring across distant phyla, despite strong conservation of regulatory elements themselves.

C. elegans is a member of a group of nematodes called rhabditids, which encompasses a large number of ecologically and genetically diverse species. A new, preliminary phylogenetic analysis is presented for concatenated sequences of three nuclear genes for 48 rhabditid and diplogastrid species (including 10 Caenorhabditis species), as well as four species representing the outgroup. Although many relationships are well-resolved, more data are still needed to resolve some key relationships, particularly near the base of the rhabditid tree. There is high confidence for two major clades: (1) a clade comprising Mesorhabditis Parasitorhabditis, Pelodera, Teratorhabditis plus a few other species; (2) a large clade (Eurhabditis) comprising most of the remaining rhabditid genera, including Caenorhabditis and its sistergroup Protorhabditis-Prodontorhabditis-Diploscapter. Eurhabditis also contains the parasitic strongylids, the entomopathogenic Heterorhabditis, and the monophyletic group Oscheius which includes the satellite model organism O. tipulae. The relationships within Caenorhabditis are well resolved. The analysis also suggests that rhabditids include diplogastrids, to which the second satellite model organism Pristionchus pacificus belongs. Genetic disparity within Caenorhabditis is as great as that across vertebrates, suggesting Caenorhabditis lineages are quickly evolving, ancient, or both. The phylogenetic tree can be used to reconstruct evolutionary events within rhabditids. For instance, the reproductive mode changed multiple times from gonochorism to hermaphroditism, but only once from hermaphroditism to gonochorism. Complete retraction of the male tail tip, leading to a blunt, peloderan tail, evolved at least once. Reversions to unretracted tail tips occurred within both major rhabditid groups. The phylogeny also provides a guide to species which would be good candidates for future genome projects and comparative studies.

The Caenorhabditis elegans pos-1 gene encodes a zinc-finger protein that is required for germline specification during embryogenesis. The maternally provided mRNA is translationally regulated both spatially and temporally during early development. We have cloned orthologs of pos-1 from C. briggsae and C. remanei, two Caenorhabditis species that have diverged from C. elegans by approximately 20-40 million years. Two regions in the 3' untranslated region are highly conserved among all three species. We find that the pos-1 RNA is expressed in the hermaphrodite and female gonads of C. briggsae and C. remanei but POS-1 protein is not detected at high levels in C. briggsae until the 2-cell stage of embryogenesis. The protein expression is restricted to the germline precursors of the embryo. We conclude that pos-1 appears to be translationally regulated in C. briggsae as it is in C. elegans and speculate the conserved 3' UTR sequences may be involved.