Widespread changes in gene expression accompany body size evolution in nematodes

Body size is a fundamental trait that drives multiple evolutionary and ecological patterns. Caenorhabditis inopinata is a fig-associated nematode that is exceptionally large relative to other members of the genus, including C. elegans. We previously showed that C. inopinata is large primarily due to postembryonic cell size expansion that occurs during the larval-to-adult transition. Here, we describe gene expression patterns in C. elegans and C. inopinata throughout this developmental period to understand the transcriptional basis of body size change. We performed RNA-seq in both species across the L3, L4, and adult stages. Most genes are differentially expressed across all developmental stages, consistent with C. inopinata's divergent ecology and morphology. We also used a model comparison approach to identify orthologs with divergent dynamics across this developmental period between the two species. This included genes connected to neurons, behavior, stress response, developmental timing, and small RNA/chromatin regulation. Multiple hypodermal collagens were also observed to harbor divergent developmental dynamics across this period, and genes important for molting and body morphology were also detected. Genes associated with TGF-β signaling revealed idiosyncratic and unexpected transcriptional patterns given their role in body size regulation in C. elegans. Widespread transcriptional divergence between these species is unexpected and may be a signature of the ecological and morphological divergence of C. inopinata. Alternatively, transcriptional turnover may be the rule in the Caenorhabditis genus, indicative of widespread developmental system drift among species. This work lays the foundation for future functional genetic studies interrogating the bases of body size evolution in this group.

An Emerging Animal Model for Querying the Role of Whole Genome Duplication in Development, Evolution, and Disease

Whole genome duplication (WGD) or polyploidization can occur at the cellular, tissue, and organismal levels. At the cellular level, tetraploidization has been proposed as a driver of aneuploidy and genome instability and correlates strongly with cancer progression, metastasis, and the development of drug resistance. WGD is also a key developmental strategy for regulating cell size, metabolism, and cellular function. In specific tissues, WGD is involved in normal development (e.g., organogenesis), tissue homeostasis, wound healing, and regeneration. At the organismal level, WGD propels evolutionary processes such as adaptation, speciation, and crop domestication. An essential strategy to further our understanding of the mechanisms promoting WGD and its effects is to compare isogenic strains that differ only in their ploidy. Caenorhabditis elegans (C. elegans) is emerging as an animal model for these comparisons, in part because relatively stable and fertile tetraploid strains can be produced rapidly from nearly any diploid strain. Here, we review the use of Caenorhabditis polyploids as tools to understand important developmental processes (e.g., sex determination, dosage compensation, and allometric relationships) and cellular processes (e.g., cell cycle regulation and chromosome dynamics during meiosis). We also discuss how the unique characteristics of the C. elegans WGD model will enable significant advances in our understanding of the mechanisms of polyploidization and its role in development and disease.

Widespread sex ratio polymorphism in Caenorhabditis nematodes

Although equal sex ratio is ubiquitous and represents an equilibrium in evolutionary theory, biased sex ratios are predicted for certain local conditions. Cases of sex ratio bias have been mostly reported for single species, but little is known about its evolution above the species level. Here, we surveyed progeny sex ratios in 23 species of the nematode genus Caenorhabditis, including 19 for which we tested multiple strains. For the species with multiple strains, five species had female-biased and two had non-biased sex ratios in all strains, respectively. The other 12 species showed polymorphic sex ratios across strains. Female-biased sex ratios could be due to sperm competition whereby X-bearing sperm outcompete nullo-X sperm during fertilization. In this model, when sperm are limited allowing all sperm to be used, sex ratios are expected to be equal. However, in assays limiting mating to a few hours, most strains showed similarly biased sex ratios compared with unlimited mating experiments, except that one C. becei strain showed significantly reduced female bias compared with unlimited mating. Our study shows frequent polymorphism in sex ratios within Caenorhabditis species and that sperm competition alone cannot explain the sex ratio bias.

Diet can alter the cost of resistance to a natural parasite in Caenorhabditis elegans

Resistance to parasites confers a fitness advantage, yet hosts show substantial variation in resistance in natural populations. Evolutionary theory indicates that resistant and susceptible genotypes can coexist if resistance is costly, but there is mixed evidence that resistant individuals have lower fitness in the absence of parasites. One explanation for this discrepancy is that the cost of resistance varies with environmental context. We tested this hypothesis using Caenorhabditis elegans and its natural microsporidian parasite, Nematocida ironsii. We used multiple metrics to compare the fitness of two near-isogenic host genotypes differing at regions associated with resistance to N. ironsii. To quantify the effect of the environment on the cost associated with these known resistance regions, we measured fitness on three microbial diets. We found that the cost of resistance varied with both diet and the measure of fitness. We detected no cost to resistance, irrespective of diet, when fitness was measured as fecundity. However, we detected a cost when fitness was measured in terms of population growth, and the magnitude of this cost varied with diet. These results provide a proof of concept that, by mediating the cost of resistance, environmental context may govern the rate and nature of resistance evolution in heterogeneous environments.

Developing an empirical model for spillover and emergence: Orsay virus host range in Caenorhabditis

A lack of tractable experimental systems in which to test hypotheses about the ecological and evolutionary drivers of disease spillover and emergence has limited our understanding of these processes. Here we introduce a promising system: Caenorhabditis hosts and Orsay virus, a positive-sense single-stranded RNA virus that naturally infects C. elegans. We assayed species across the Caenorhabditis tree and found Orsay virus susceptibility in 21 of 84 wild strains belonging to 14 of 44 species. Confirming patterns documented in other systems, we detected effects of host phylogeny on susceptibility. We then tested whether susceptible strains were capable of transmitting Orsay virus by transplanting exposed hosts and determining whether they transmitted infection to conspecifics during serial passage. We found no evidence of transmission in 10 strains (virus undetectable after passaging in all replicates), evidence of low-level transmission in 5 strains (virus lost between passage 1 and 5 in at least one replicate) and evidence of sustained transmission in 6 strains (including all three experimental C. elegans strains) in at least one replicate. Transmission was strongly associated with viral amplification in exposed populations. Variation in Orsay virus susceptibility and transmission among Caenorhabditis strains suggests that the system could be powerful for studying spillover and emergence.

First report of Caenorhabditis brenneri (Nematoda: Rhabditida) isolated from the cadaver of Philippinella moellendorffi (Stylommatophora: Ariophantidae), a terrestrial slug in the Philippines

Gastropod-associated nematodes have been previously studied and documented worldwide, with some species forming host-specific association as obligate parasites of molluscs while others form intermediate and temporary association. Philippinella moellendorffi from Imelda, Zamboanga Sibugay, Philippines, were collected, washed and maintained in the laboratory until death. Cadavers were placed on nutrient agar to allow nematode proliferation. Nematode pure culture was obtained using one gravid female for propagation. Morphology and molecular analyses (18S ribosomal DNA (rDNA) and D2-D3 expansion segments of 28S rDNA) were employed as diagnostic tools in identifying the nematode species isolated from P. moellendorffi. The newly isolated nematode was identified as Caenorhabditis brenneri, thus designated as C. brenneri strain IZSP from the Philippines. This is the first record of C. brenneri isolated from the terrestrial slug P. moellendorffi.

Caenorhabditis nematodes colonize ephemeral resource patches in neotropical forests

Factors shaping the distribution and abundance of species include life-history traits, population structure, and stochastic colonization-extinction dynamics. Field studies of model species groups help reveal the roles of these factors. Species of Caenorhabditis nematodes are highly divergent at the sequence level but exhibit highly conserved morphology, and many of these species live in sympatry on microbe-rich patches of rotten material. Here, we use field experiments and large-scale opportunistic collections to investigate species composition, abundance, and colonization efficiency of Caenorhabditis species in two of the world's best-studied lowland tropical field sites: Barro Colorado Island in Panamá and La Selva in Sarapiquí, Costa Rica. We observed seven species of Caenorhabditis, four of them known only from these collections. We formally describe two species and place them within the Caenorhabditis phylogeny. While these localities contain species from many parts of the phylogeny, both localities were dominated by globally distributed androdiecious species. We found that Caenorhabditis individuals were able to colonize baits accessible only through phoresy and preferentially colonized baits that were in direct contact with the ground. We estimate the number of colonization events per patch to be low.

Mutation, selection, and the prevalence of the Caenorhabditis elegans heat-sensitive mortal germline phenotype

Caenorhabditis elegans strains with the heat-sensitive mortal germline phenotype become progressively sterile over the course of a few tens of generations when maintained at temperatures near the upper range of C. elegans' tolerance. Mortal germline is transgenerationally heritable, and proximately under epigenetic control. Previous studies have suggested that mortal germline presents a relatively large mutational target and that mortal germline is not uncommon in natural populations of C. elegans. The mortal germline phenotype is not monolithic. Some strains exhibit a strong mortal germline phenotype, in which individuals invariably become sterile over a few generations, whereas other strains show a weaker (less penetrant) phenotype in which the onset of sterility is slower and more stochastic. We present results in which we (1) quantify the rate of mutation to the mortal germline phenotype and (2) quantify the frequency of mortal germline in a collection of 95 wild isolates. Over the course of ∼16,000 meioses, we detected one mutation to a strong mortal germline phenotype, resulting in a point estimate of the mutation rate UMrt≈ 6×10-5/genome/generation. We detected no mutations to a weak mortal germline phenotype. Six out of 95 wild isolates have a strong mortal germline phenotype, and although quantification of the weak mortal germline phenotype is inexact, the weak mortal germline phenotype is not rare in nature. We estimate a strength of selection against mutations conferring the strong mortal germline phenotype s¯≈0.1%, similar to selection against mutations affecting competitive fitness. The appreciable frequency of weak mortal germline variants in nature combined with the low mutation rate suggests that mortal germline may be maintained by balancing selection.

Small-Molecule Mn Antioxidants in Caenorhabditis elegans and Deinococcus radiodurans Supplant MnSOD Enzymes during Aging and Irradiation

Denham Harman's oxidative damage theory identifies superoxide (O₂^•-) radicals as central agents of aging and radiation injury, with Mn²⁺-dependent superoxide dismutase (MnSOD) as the principal O₂^•--scavenger. However, in the radiation-resistant nematode Caenorhabditis elegans, the mitochondrial antioxidant enzyme MnSOD is dispensable for longevity, and in the model bacterium Deinococcus radiodurans, it is dispensable for radiation resistance. Many radiation-resistant organisms accumulate small-molecule Mn²⁺-antioxidant complexes well-known for their catalytic ability to scavenge O₂^•-, along with MnSOD, as exemplified by D. radiodurans. Here, we report experiments that relate the MnSOD and Mn-antioxidant content to aging and oxidative stress resistances and which indicate that C. elegans, like D. radiodurans, may rely on Mn-antioxidant complexes as the primary defense against reactive oxygen species (ROS). Wild-type and ΔMnSOD D. radiodurans and C. elegans were monitored for gamma radiation sensitivities over their life spans while gauging Mn²⁺-antioxidant content by electron paramagnetic resonance (EPR) spectroscopy, a powerful new approach to determining the in vivo Mn-antioxidant content of cells as they age. As with D. radiodurans, MnSOD is dispensable for radiation survivability in C. elegans, which hyperaccumulates Mn-antioxidants exceptionally protective of proteins. Unexpectedly, ΔMnSOD mutants of both the nematodes and bacteria exhibited increased gamma radiation survival compared to the wild-type. In contrast, the loss of MnSOD renders radiation-resistant bacteria sensitive to atmospheric oxygen during desiccation. Our results support the concept that the disparate responses to oxidative stress are explained by the accumulation of Mn-antioxidant complexes which protect, complement, and can even supplant MnSOD. IMPORTANCE The current theory of cellular defense against oxidative damage identifies antioxidant enzymes as primary defenders against ROS, with MnSOD being the preeminent superoxide (O₂^•-) scavenger. However, MnSOD is shown to be dispensable both for radiation resistance and longevity in model organisms, the bacterium Deinococcus radiodurans and the nematode Caenorhabditis elegans. Measured by electron paramagnetic resonance (EPR) spectroscopy, small-molecule Mn-antioxidant content was shown to decline in unison with age-related decreases in cell proliferation and radioresistance, which again are independent of MnSOD presence. Most notably, the Mn-antioxidant content of C. elegans drops precipitously in the last third of its life span, which links with reports that the steady-state level of oxidized proteins increases exponentially during the last third of the life span in animals. This leads us to propose that global responses to oxidative stress must be understood through an extended theory that includes small-molecule Mn-antioxidants as potent O₂^•--scavengers that complement, and can even supplant, MnSOD.

Imaging Glycosaminoglycan Modification Patterns In Vivo

Glycosaminoglycans (GAGs) such as heparan sulfates (HS) or chondroitin sulfates (CS) are long unbranched polymers of a disaccharide comprised of hexuronic acid and hexosamine. Attached to a protein backbone via a characteristic tetrasaccharide, the GAG chains are non-uniformly modified by sulfations, epimerizations, and deacetylations. The resultant glycan chains contain highly modified domains, separated by sections of sparse or no modifications. These GAG domains are central to the role of glycans in binding to proteins and mediating protein-protein interactions. Since HS and CS domains are not genetically encoded, they cannot be visualized and studied with conventional methods in vivo. We describe a transgenic approach using single chain variable fragment (scFv) antibodies that bind HS or CS. By transgenically expressing fluorescently tagged scFv antibodies, we can directly visualize both HS and CS domains in live Caenorhabditis elegans revealing unprecedented cellular specificity and evolutionary conservation (Attreed et al., Nat Methods 9(5): 477-479, 2012; Attreed et al., Glycobiology 26(8): 862-870, 2016) (unpublished). The approach allows concomitant co-labeling of multiple GAG domains, the study of GAG dynamics, and could lend itself to a genetic analysis of GAG domain biosynthesis or function.

Genome-wide characterization, evolution, structure, and expression analysis of the F-box genes in Caenorhabditis

BACKGROUND

F-box proteins represent a diverse class of adaptor proteins of the ubiquitin-proteasome system (UPS) that play critical roles in the cell cycle, signal transduction, and immune response by removing or modifying cellular regulators. Among closely related organisms of the Caenorhabditis genus, remarkable divergence in F-box gene copy numbers was caused by sizeable species-specific expansion and contraction. Although F-box gene number expansion plays a vital role in shaping genomic diversity, little is known about molecular evolutionary mechanisms responsible for substantial differences in gene number of F-box genes and their functional diversification in Caenorhabditis. Here, we performed a comprehensive evolution and underlying mechanism analysis of F-box genes in five species of Caenorhabditis genus, including C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei.

RESULTS

Herein, we identified and characterized 594, 192, 377, 39, 1426 F-box homologs encoding putative F-box proteins in the genome of C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. Our work suggested that extensive species-specific tandem duplication followed by a small amount of gene loss was the primary mechanism responsible for F-box gene number divergence in Caenorhabditis genus. After F-box gene duplication events occurred, multiple mechanisms have contributed to gene structure divergence, including exon/intron gain/loss, exonization/pseudoexonization, exon/intron boundaries alteration, exon splits, and intron elongation by tandem repeats. Based on high-throughput RNA sequencing data analysis, we proposed that F-box gene functions have diversified by sub-functionalization through highly divergent stage-specific expression patterns in Caenorhabditis species.

CONCLUSIONS

Massive species-specific tandem duplications and occasional gene loss drove the rapid evolution of the F-box gene family in Caenorhabditis, leading to complex gene structural variation and diversified functions affecting growth and development within and among Caenorhabditis species. In summary, our findings outline the evolution of F-box genes in the Caenorhabditis genome and lay the foundation for future functional studies.

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

Cross-Predicting Essential Genes between Two Model Eukaryotic Species Using Machine Learning

Experimental studies of Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular and cellular processes in metazoans at large. Since the publication of their genomes, functional genomic investigations have identified genes that are essential or non-essential for survival in each species. Recently, a range of features linked to gene essentiality have been inferred using a machine learning (ML)-based approach, allowing essentiality predictions within a species. Nevertheless, predictions between species are still elusive. Here, we undertake a comprehensive study using ML to discover and validate features of essential genes common to both C. elegans and D. melanogaster. We demonstrate that the cross-species prediction of gene essentiality is possible using a subset of features linked to nucleotide/protein sequences, protein orthology and subcellular localisation, single-cell RNA-seq, and histone methylation markers. Complementary analyses showed that essential genes are enriched for transcription and translation functions and are preferentially located away from heterochromatin regions of C. elegans and D. melanogaster chromosomes. The present work should enable the cross-prediction of essential genes between model and non-model metazoans.

Competitive fitness analysis using Convolutional Neural Network

We developed a procedure for estimating competitive fitness by using Caenorhabditis elegans as a model organism and a Convolutional Neural Network (CNN) as a tool. Competitive fitness is usually the most informative fitness measure, and competitive fitness assays often rely on green fluorescent protein (GFP) marker strains. CNNs are a class of deep learning neural networks, which are well suited for image analysis and object classification. Our model analyses involved image classification of nematodes as wild-type vs. GFP-expressing, and counted both categories. The performance was analyzed with (i) precision and recall parameters, and (ii) comparison of the wild-type frequency calculated from the model against that obtained by visual scoring of the same images. The average precision and recall varied from 0.79 to 0.87 and from 0.84 to 0.92, respectively, depending on worm density in the images. Compared with manual counting, the model decreased counting time at least 20-fold while preventing human errors. Given the rapid development in the field of CNN, the model, which is fully available on GitHub, can be further optimized and adapted for other image-based uses.

Cultivation of Caenorhabditis elegans on new cheap monoxenic media without peptone

The study of species biodiversity within the Caenorhabditis genus of nematodes would be facilitated by the isolation of as many species as possible. So far, over 50 species have been found, usually associated with decaying vegetation or soil samples, with many from Africa, South America and Southeast Asia. Scientists based in these regions can contribute to Caenorhabditis sampling and their proximity would allow intensive sampling, which would be useful for understanding the natural history of these species. However, severely limited research budgets are often a constraint for these local scientists. In this study, we aimed to find a more economical, alternative growth media to rear Caenorhabditis and related species. We tested 25 media permutations using cheaper substitutes for the reagents found in the standard nematode growth media (NGM) and found three media combinations that performed comparably to NGM with respect to the reproduction and longevity of C. elegans. These new media should facilitate the isolation and characterization of Caenorhabditis and other free-living nematodes for the researchers in the poorer regions such as Africa, South America, and Southeast Asia where nematode diversity appears high.

Ubiquitous Selfish Toxin-Antidote Elements in Caenorhabditis Species

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₂ 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.

A larger target leads to faster evolution

The speed at which a cell fate decision in nematode worms evolves is due to the number of genes that control the decision, rather than to a high mutation rate.

Environmental variation mediates the evolution of anticipatory parental effects

Theory maintains that when future environment is predictable, parents should adjust the phenotype of their offspring to match the anticipated environment. The plausibility of positive anticipatory parental effects is hotly debated and the experimental evidence for the evolution of such effects is currently lacking. We experimentally investigated the evolution of anticipatory maternal effects in a range of environments that differ drastically in how predictable they are. Populations of the nematode Caenorhabditis remanei, adapted to 20°C, were exposed to a novel temperature (25°C) for 30 generations with either positive or zero correlation between parent and offspring environment. We found that populations evolving in novel environments that were predictable across generations evolved a positive anticipatory maternal effect, because they required maternal exposure to 25°C to achieve maximum reproduction in that temperature. In contrast, populations evolving under zero environmental correlation had lost this anticipatory maternal effect. Similar but weaker patterns were found if instead rate‐sensitive population growth was used as a fitness measure. These findings demonstrate that anticipatory parental effects evolve in response to environmental change so that ill‐fitting parental effects can be rapidly lost. Evolution of positive anticipatory parental effects can aid population viability in rapidly changing but predictable environments.

Evolution and Developmental System Drift in the Endoderm Gene Regulatory Network of Caenorhabditis and Other Nematodes

Developmental gene regulatory networks (GRNs) underpin metazoan embryogenesis and have undergone substantial modification to generate the tremendous variety of animal forms present on Earth today. The nematode Caenorhabditis elegans has been a central model for advancing many important discoveries in fundamental mechanistic biology and, more recently, has provided a strong base from which to explore the evolutionary diversification of GRN architecture and developmental processes in other species. In this short review, we will focus on evolutionary diversification of the GRN for the most ancient of the embryonic germ layers, the endoderm. Early embryogenesis diverges considerably across the phylum Nematoda. Notably, while some species deploy regulative development, more derived species, such as C. elegans, exhibit largely mosaic modes of embryogenesis. Despite the relatively similar morphology of the nematode gut across species, widespread variation has been observed in the signaling inputs that initiate the endoderm GRN, an exemplar of developmental system drift (DSD). We will explore how genetic variation in the endoderm GRN helps to drive DSD at both inter- and intraspecies levels, thereby resulting in a robust developmental system. Comparative studies using divergent nematodes promise to unveil the genetic mechanisms controlling developmental plasticity and provide a paradigm for the principles governing evolutionary modification of an embryonic GRN.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

Diversification of the Caenorhabditis heat shock response by Helitron transposable elements

Heat Shock Factor 1 (HSF-1) is a key regulator of the heat shock response (HSR). Upon heat shock, HSF-1 binds well-conserved motifs, called Heat Shock Elements (HSEs), and drives expression of genes important for cellular protection during this stress. Remarkably, we found that substantial numbers of HSEs in multiple Caenorhabditis species reside within Helitrons, a type of DNA transposon. Consistent with Helitron-embedded HSEs being functional, upon heat shock they display increased HSF-1 and RNA polymerase II occupancy and up-regulation of nearby genes in C. elegans. Interestingly, we found that different genes appear to be incorporated into the HSR by species-specific Helitron insertions in C. elegans and C. briggsae and by strain-specific insertions among different wild isolates of C. elegans. Our studies uncover previously unidentified targets of HSF-1 and show that Helitron insertions are responsible for rewiring and diversifying the Caenorhabditis HSR.

Artificial selection for increased dispersal results in lower fitness

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.

Vertical transmission in Caenorhabditis nematodes of RNA molecules encoding a viral RNA-dependent RNA polymerase

In organisms composed of a single cell, RNAs of viral origin may be transmitted to daughter cells at cell division without passing through an extracellular virion stage. These RNAs usually encode an RNA-dependent RNA polymerase that enables their replication. For some of these agents, such as Narnaviruses, no capsid protein is expressed, and thus, they are called capsidless viruses. Here, we identify putative capsidless viral RNAs in animals, in nematodes closely related to the model organism Caenorhabditis elegans. We show that these RNAs are transmitted vertically through the host germline. Our work provides evidence that animal cells harbor capsidless viruses.

Here, we report on the discovery in Caenorhabditis nematodes of multiple vertically transmitted RNAs coding for putative RNA-dependent RNA polymerases. Their sequences share similarity to distinct RNA viruses, including bunyaviruses, narnaviruses, and sobemoviruses. The sequences are present exclusively as RNA and are not found in DNA form. The RNAs persist in progeny after bleach treatment of adult animals, indicating vertical transmission of the RNAs. We tested one of the infected strains for transmission to an uninfected strain and found that mating of infected animals with uninfected animals resulted in infected progeny. By in situ hybridization, we detected several of these RNAs in the cytoplasm of the male and female germline of the nematode host. The Caenorhabditis hosts were found defective in degrading exogenous double-stranded RNAs, which may explain retention of viral-like RNAs. Strikingly, one strain, QG551, harbored three distinct virus-like RNA elements. Specific patterns of small RNAs complementary to the different viral-like RNAs were observed, suggesting that the different RNAs are differentially recognized by the RNA interference (RNAi) machinery. While vertical transmission of viruses in the family Narnaviridae, which are known as capsidless viruses, has been described in fungi, these observations provide evidence that multicellular animal cells harbor similar viruses.

Natural Variation and Genetic Determinants of Caenorhabditis elegans Sperm Size

The diversity in sperm shape and size represents a powerful paradigm to understand how selection drives the evolutionary diversification of cell morphology. Experimental work on the sperm biology of the male-hermaphrodite nematode Caenorhabditis elegans has elucidated diverse factors important for sperm fertilization success, including the competitive superiority of larger sperm. Yet despite extensive research, the molecular mechanisms regulating C. elegans sperm size and the genetic basis underlying natural variation in sperm size remain unknown. To address these questions, we quantified male sperm size variation of a worldwide panel of 97 genetically distinct C. elegans strains, allowing us to uncover significant genetic variation in male sperm size. Aiming to characterize the molecular genetic basis of C. elegans male sperm size variation using a genome-wide association study, we did not detect any significant quantitative trait loci. We therefore focused on the genetic analysis of pronounced sperm size differences observed between recently diverged laboratory strains (N2 vs. LSJ1/2). Using mutants and quantitative complementation tests, we demonstrate that variation in the gene nurf-1 underlies the evolution of small sperm in the LSJ lineage. Given the previous discovery that this same nurf-1 variation was central for hermaphrodite laboratory adaptation, the evolution of reduced male sperm size in LSJ strains likely reflects a pleiotropic consequence. Together, our results provide a comprehensive quantification of natural variation in C. elegans sperm size and first insights into the genetic determinants of Caenorhabditis sperm size, pointing at an involvement of the NURF chromatin remodeling complex.

Genome structure predicts modular transcriptome responses to genetic and environmental conditions

Understanding the plasticity, robustness and modularity of transcriptome expression to genetic and environmental conditions is crucial to deciphering how organisms adapt in nature. To test how genome architecture influences transcriptome profiles, we quantified expression responses for distinct temperature-adapted genotypes of the nematode Caenorhabditis briggsae when exposed to chronic temperature stresses throughout development. We found that 56% of the 8,795 differentially expressed genes show genotype-specific changes in expression in response to temperature (genotype-by-environment interactions, GxE). Most genotype-specific responses occur under heat stress, indicating that cold vs. heat stress responses involve distinct genomic architectures. The 22 co-expression modules that we identified differ in their enrichment of genes with genetic vs. environmental vs. interaction effects, as well as their genomic spatial distributions, functional attributes and rates of molecular evolution at the sequence level. Genes in modules enriched for simple effects of either genotype or temperature alone tend to evolve especially rapidly, consistent with disproportionate influence of adaptation or weaker constraint on these subsets of loci. Chromosome-scale heterogeneity in nucleotide polymorphism, however, rather than the scale of individual genes predominates as the source of genetic differences among expression profiles, and natural selection regimes are largely decoupled between coding sequences and noncoding flanking sequences that contain cis-regulatory elements. These results illustrate how the form of transcriptome modularity and genome structure contribute to predictable profiles of evolutionary change.

Genetically Distinct Behavioral Modules Underlie Natural Variation in Thermal Performance Curves

Thermal reaction norms pervade organismal traits as stereotyped responses to temperature, a fundamental environmental input into sensory and physiological systems. Locomotory behavior represents an especially plastic read-out of animal response, with its dynamic dependence on environmental stimuli presenting a challenge for analysis and for understanding the genomic architecture of heritable variation. Here we characterize behavioral reaction norms as thermal performance curves for the nematode Caenorhabditis briggsae, using a collection of 23 wild isolate genotypes and 153 recombinant inbred lines to quantify the extent of genetic and plastic variation in locomotory behavior to temperature changes. By reducing the dimensionality of the multivariate phenotypic response with a function-valued trait framework, we identified genetically distinct behavioral modules that contribute to the heritable variation in the emergent overall behavioral thermal performance curve. Quantitative trait locus mapping isolated regions on Chromosome II associated with locomotory activity at benign temperatures and Chromosome V loci related to distinct aspects of sensitivity to high temperatures, with each quantitative trait locus explaining up to 28% of trait variation. These findings highlight how behavioral responses to environmental inputs as thermal reaction norms can evolve through independent changes to genetically distinct modular components of such complex phenotypes.

Rapid population-wide declines in stem cell number and activity during reproductive aging in C. elegans

C. elegans hermaphrodites display dramatic age-related decline of reproduction early in life, while somatic functions are still robust. To understand reproductive aging, we analyzed the assembly line of oocyte production that generates fertilized eggs. Aging germlines displayed both sporadic and population-wide changes. A small fraction of aging animals displayed endomitotic oocytes in the germline and other defects. By contrast, all animals displayed age-related decreases in germline size and function. As early as day 3 of adulthood, animals displayed fewer stem cells and a slower cell cycle, which combine to substantially decrease progenitor zone output. The C. elegans germline is the only adult tissue that contains stem cells, allowing the analysis of stem cells in aging. To investigate the mechanism of the decrease in stem cell number, we analyzed the Notch signaling pathway. The Notch effectors LST-1 and SYGL-1 displayed age-related decreases in expression domains, suggesting a role for Notch signaling in germline aging. The results indicate that although sporadic defects account for the sterility of some animals, population-wide changes account for the overall pattern of reproductive aging.

N-glycomic Complexity in Anatomical Simplicity: Caenorhabditis elegans as a Non-model Nematode?

Caenorhabditis elegans is a genetically well-studied model nematode or "worm"; however, its N-glycomic complexity is actually baffling and still not completely unraveled. Some features of its N-glycans are, to date, unique and include bisecting galactose and up to five fucose residues associated with the asparagine-linked Man2-3GlcNAc2 core; the substitutions include galactosylation of fucose, fucosylation of galactose and methylation of mannose or fucose residues as well as phosphorylcholine on antennal (non-reducing) N-acetylglucosamine. Only some of these modifications are shared with various other nematodes, while others have yet to be detected in any other species. Thus, C. elegans can be used as a model for some aspects of N-glycan function, but its glycome is far from identical to those of other organisms and is actually far from simple. Possibly the challenges of its native environment, which differ from those of parasitic or necromenic species, led to an anatomically simple worm possessing a complex glycome.

Contrasting patterns of molecular evolution in metazoan germ line genes

BACKGROUND

Germ lines are the cell lineages that give rise to the sperm and eggs in animals. The germ lines first arise from primordial germ cells (PGCs) during embryogenesis: these form from either a presumed derived mode of preformed germ plasm (inheritance) or from an ancestral mechanism of inductive cell-cell signalling (induction). Numerous genes involved in germ line specification and development have been identified and functionally studied. However, little is known about the molecular evolutionary dynamics of germ line genes in metazoan model systems.

RESULTS

Here, we studied the molecular evolution of germ line genes within three metazoan model systems. These include the genus Drosophila (N=34 genes, inheritance), the fellow insect Apis (N=30, induction), and their more distant relative Caenorhabditis (N=23, inheritance). Using multiple species and established phylogenies in each genus, we report that germ line genes exhibited marked variation in the constraint on protein sequence divergence (dN/dS) and codon usage bias (CUB) within each genus. Importantly, we found that de novo lineage-specific inheritance (LSI) genes in Drosophila (osk, pgc) and in Caenorhabditis (pie-1, pgl-1), which are essential to germ plasm functions under the derived inheritance mode, displayed rapid protein sequence divergence relative to the other germ line genes within each respective genus. We show this may reflect the evolution of specialized germ plasm functions and/or low pleiotropy of LSI genes, features not shared with other germ line genes. In addition, we observed that the relative ranking of dN/dS and of CUB between genera were each more strongly correlated between Drosophila and Caenorhabditis, from different phyla, than between Drosophila and its insect relative Apis, suggesting taxonomic differences in how germ line genes have evolved.

CONCLUSIONS

Taken together, the present results advance our understanding of the evolution of animal germ line genes within three well-known metazoan models. Further, the findings provide insights to the molecular evolution of germ line genes with respect to LSI status, pleiotropy, adaptive evolution as well as PGC-specification mode.

Demographic consequences of reproductive interference in multi-species communities

BACKGROUND

Reproductive interference can mediate interference competition between species through sexual interactions that reduce the fitness of one species by another. Theory shows that the positive frequency-dependent effects of such costly errors in mate recognition can dictate species coexistence or exclusion even with countervailing resource competition differences between species. While usually framed in terms of pre-mating or post-zygotic costs, reproductive interference manifests between individual Caenorhabditis nematodes from negative interspecies gametic interactions: sperm cells from interspecies matings can migrate ectopically to induce female sterility and premature death. The potential for reproductive interference to exert population level effects on Caenorhabditis trait evolution and community structure, however, remains unknown.

RESULTS

Here we test whether a species that is superior in individual-level reproductive interference (C. nigoni) can exact negative demographic effects on competitor species that are superior in resource competition (C. briggsae and C. elegans). We observe coexistence over six generations and find evidence of demographic reproductive interference even under conditions unfavorable to its influence. C. briggsae and C. elegans show distinct patterns of reproductive interference in competitive interactions with C. nigoni.

CONCLUSIONS

These results affirm that individual level negative effects of reproductive interference mediated by gamete interactions can ramify to population demography, with the potential to influence patterns of species coexistence separately from the effects of direct resource competition.

Physiological Starvation Promotes Caenorhabditis elegans Vulval Induction

Studying how molecular pathways respond to ecologically relevant environmental variation is fundamental to understand organismal development and its evolution. Here we characterize how starvation modulates Caenorhabditis elegans vulval cell fate patterning – an environmentally sensitive process, with a nevertheless robust output. Past research has shown many vulval mutants affecting EGF-Ras-MAPK, Delta-Notch and Wnt pathways to be suppressed by environmental factors, such as starvation. Here we aimed to resolve previous, seemingly contradictory, observations on how starvation modulates levels of vulval induction. Using the strong starvation suppression of the Vulvaless phenotype of lin-3/egf reduction-of-function mutations as an experimental paradigm, we first tested for a possible involvement of the sensory system in relaying starvation signals to affect vulval induction: mutation of various sensory inputs, DAF-2/Insulin or DAF-7/TGF-β signaling did not abolish lin-3(rf) starvation suppression. In contrast, nutrient deprivation induced by mutation of the intestinal peptide transporter gene pept-1 or the TOR pathway component rsks-1 (the ortholog of mammalian P70S6K) very strongly suppressed lin-3(rf) mutant phenotypes. Therefore, physiologically starved animals induced by these mutations tightly recapitulated the effects of external starvation on vulval induction. While both starvation and pept-1 RNAi were sufficient to increase Ras and Notch pathway activities in vulval cells, the highly penetrant Vulvaless phenotype of a tissue-specific null allele of lin-3 was not suppressed by either condition. This and additional results indicate that partial lin-3 expression is required for starvation to affect vulval induction. These results suggest a cross-talk between nutrient deprivation, TOR-S6K and EGF-Ras-MAPK signaling during C. elegans vulval induction.

Field studies reveal a close relative of C. elegans thrives in the fresh figs of Ficus septica and disperses on its Ceratosolen pollinating wasps

BACKGROUND

Biotic interactions are ubiquitous and require information from ecology, evolutionary biology, and functional genetics in order to be understood. However, study systems that are amenable to investigations across such disparate fields are rare. Figs and fig wasps are a classic system for ecology and evolutionary biology with poor functional genetics; Caenorhabditis elegans is a classic system for functional genetics with poor ecology. In order to help bridge these disciplines, here we describe the natural history of a close relative of C. elegans, Caenorhabditis inopinata, that is associated with the fig Ficus septica and its pollinating Ceratosolen wasps.

RESULTS

To understand the natural context of fig-associated Caenorhabditis, fresh F. septica figs from four Okinawan islands were sampled, dissected, and observed under microscopy. C. inopinata was found in all islands where F. septica figs were found. C.i nopinata was routinely found in the fig interior and almost never observed on the outside surface. C. inopinata was only found in pollinated figs, and C. inopinata was more likely to be observed in figs with more foundress pollinating wasps. Actively reproducing C. inopinata dominated early phase figs, whereas late phase figs with emerging wasp progeny harbored C. inopinata dauer larvae. Additionally, C. inopinata was observed dismounting from Ceratosolen pollinating wasps that were placed on agar plates. C. inopinata was not found on non-pollinating, parasitic Philotrypesis wasps. Finally, C. inopinata was only observed in F. septica figs among five Okinawan Ficus species sampled.

CONCLUSION

These are the first detailed field observations of C. inopinata, and they suggest a natural history where this species proliferates in early phase F. septica figs and disperses from late phase figs on Ceratosolen pollinating fig wasps. While consistent with other examples of nematode diversification in the fig microcosm, the fig and wasp host specificity of C. inopinata is highly divergent from the life histories of its close relatives and frames hypotheses for future investigations. This natural co-occurrence of the fig/fig wasp and C. inopinata study systems sets the stage for an integrated research program that can help to explain the evolution of interspecific interactions.

X exceptionalism in Caenorhabditis speciation

Speciation genetics research in diverse organisms shows the X-chromosome to be exceptional in how it contributes to "rules" of speciation. Until recently, however, the nematode phylum has been nearly silent on this issue, despite the model organism Caenorhabditis elegans having touched most other topics in biology. Studies of speciation with Caenorhabditis accelerated with the recent discovery of species pairs showing partial interfertility. The resulting genetic analyses of reproductive isolation in nematodes demonstrate key roles for the X-chromosome in hybrid male sterility and inviability, opening up new understanding of the genetic causes of Haldane's rule, Darwin's corollary to Haldane's rule, and enabling tests of the large-X effect hypothesis. Studies to date implicate improper chromatin regulation of the X-chromosome by small RNA pathways as integral to hybrid male dysfunction. Sexual transitions in reproductive mode to self-fertilizing hermaphroditism inject distinctive molecular evolutionary features into the speciation process for some species. Caenorhabditis also provides unique opportunities for analysis in a system with XO sex determination that lacks a Y-chromosome, sex chromosome-dependent sperm competition differences and mechanisms of gametic isolation, exceptional accessibility to the development process and rapid experimental evolution. As genetic analysis of reproductive isolation matures with investigation of multiple pairs of Caenorhabditis species and new species discovery, nematodes will provide a powerful complement to more established study organisms for deciphering the genetic basis of and rules to speciation.

Studying placebo effects in model organisms will help us understand them in humans

The placebo effect is widely recognized but important questions remain, for example whether the capacity to respond to a placebo is an evolved, and potentially ubiquitous trait, or an unpredictable side effect of another evolved process. Understanding this will determine the degree to which the physiology underlying placebo effects might be manipulated or harnessed to optimize medical treatments. We argue that placebo effects are cases of phenotypic plasticity where once predictable cues are now unpredictable. Importantly, this explains why placebo-like effects are observed in less complex organisms such as worms and flies. Further, this indicates that such species present significant opportunities to test hypotheses that would be ethically or pragmatically impossible in humans. This paradigm also suggests that data informative of human placebo effects pre-exist in studies of model organisms.

A Model for Evolutionary Ecology of Disease: The Case for Caenorhabditis Nematodes and Their Natural Parasites

Many of the outstanding questions in disease ecology and evolution call for combining observation of natural host-parasite populations with experimental dissection of interactions in the field and the laboratory. The "rewilding" of model systems holds great promise for this endeavor. Here, we highlight the potential for development of the nematode Caenorhabditis elegans and its close relatives as a model for the study of disease ecology and evolution. This powerful laboratory model was disassociated from its natural habitat in the 1960s. Today, studies are uncovering that lost natural history, with several natural parasites described since 2008. Studies of these natural Caenorhabditis-parasite interactions can reap the benefits of the vast array of experimental and genetic tools developed for this laboratory model. In this review, we introduce the natural parasites of C. elegans characterized thus far and discuss resources available to study them, including experimental (co)evolution, cryopreservation, behavioral assays, and genomic tools. Throughout, we present avenues of research that are interesting and feasible to address with caenorhabditid nematodes and their natural parasites, ranging from the maintenance of outcrossing to the community dynamics of host-associated microbes. In combining natural relevance with the experimental power of a laboratory supermodel, these fledgling host-parasite systems can take on fundamental questions in evolutionary ecology of disease.

Differences in the genetic control of early egg development and reproduction between C. elegans and its parthenogenetic relative D. coronatus

BACKGROUND

The free-living nematode Diploscapter coronatus is the closest known relative of Caenorhabditis elegans with parthenogenetic reproduction. It shows several developmental idiosyncracies, for example concerning the mode of reproduction, embryonic axis formation and early cleavage pattern (Lahl et al. in Int J Dev Biol 50:393-397, 2006). Our recent genome analysis (Hiraki et al. in BMC Genomics 18:478, 2017) provides a solid foundation to better understand the molecular basis of developmental idiosyncrasies in this species in an evolutionary context by comparison with selected other nematodes. Our genomic data also yielded indications for the view that D. coronatus is a product of interspecies hybridization.

RESULTS

In a genomic comparison between D. coronatus, C. elegans, other representatives of the genus Caenorhabditis and the more distantly related Pristionchus pacificus and Panagrellus redivivus, certain genes required for central developmental processes in C. elegans like control of meiosis and establishment of embryonic polarity were found to be restricted to the genus Caenorhabditis. The mRNA content of early D. coronatus embryos was sequenced and compared with similar stages in C. elegans and Ascaris suum. We identified 350 gene families transcribed in the early embryo of D. coronatus but not in the other two nematodes. Looking at individual genes transcribed early in D. coronatus but not in C. elegans and A. suum, we found that orthologs of most of these are present in the genomes of the latter species as well, suggesting heterochronic shifts with respect to expression behavior. Considerable genomic heterozygosity and allelic divergence lend further support to the view that D. coronatus may be the result of an interspecies hybridization. Expression analysis of early acting single-copy genes yields no indication for silencing of one parental genome.

CONCLUSIONS

Our comparative cellular and molecular studies support the view that the genus Caenorhabditis differs considerably from the other studied nematodes in its control of development and reproduction. The easy-to-culture parthenogenetic D. coronatus, with its high-quality draft genome and only a single chromosome when haploid, offers many new starting points on the cellular, molecular and genomic level to explore alternative routes of nematode development and reproduction.

A Behavioral Survey of the Effects of Kavalactones on Caenorhabditis elegans Neuromuscular Transmission

Kava is a plant root extract that is widely consumed by Pacific Islanders. Kava contains a class of lactone compounds called kavalactones. The sedative and anxiolytic effects of kava are likely attributed to the efficacies of kavalactones on the nervous system. Although some studies have implicated the potencies of certain kavalactone species on γ-aminobutyric acid transmission, evidence supporting the action of kavalactones on the eukaryotic neuromuscular junction (NMJ) and acetylcholine (ACh) transmission is scant. Here, we used behavioral assays to demonstrate the effects of kavalactones at the Caenorhabditis elegans NMJ. Our results suggest that kavalactones disrupt the inhibitory-excitatory balance at the NMJ. Such perturbation of NMJ activity is likely due to excess or prolonged ACh transmission. In addition, we found that kavain, a major constituent of kava, induced worm paralysis but not convulsions. Hence, the modulatory action of kavain could be distinct from the other kavalactone species.

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

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

Mating and male pheromone kill Caenorhabditis males through distinct mechanisms

Differences in longevity between sexes is a mysterious yet general phenomenon across great evolutionary distances. To test the roles of responses to environmental cues and sexual behaviors in longevity regulation, we examined Caenorhabditis male lifespan under solitary, grouped, and mated conditions. We find that neurons and the germline are required for male pheromone-dependent male death. Hermaphrodites with a masculinized nervous system secrete male pheromone and are susceptible to male pheromone killing. Male pheromone-mediated killing is unique to androdioecious Caenorhabditis, and may reduce the number of males in hermaphroditic populations; neither males nor females of gonochoristic species are susceptible to male pheromone killing. By contrast, mating-induced death, which is characterized by germline-dependent shrinking, glycogen loss, and ectopic vitellogenin expression, utilizes distinct molecular pathways and is shared between the sexes and across species. The study of sex- and species-specific regulation of aging reveals deeply conserved mechanisms of longevity and population structure regulation.

DOI:http://dx.doi.org/10.7554/eLife.23493.001

In many animals, different sexes have different life expectancies. This holds true for a roundworm species called Caenorhabditis elegans that has commonly been used to study aging and lifespan. Unlike some related Caenorhabditis roundworm species (which consist of male and female worms), C. elegans worms are predominantly hermaphrodites and reproduce by self-fertilization. C. elegans males are normally rare. However, under stressful conditions the number of males increases to reduce inbreeding and so help the worm population to adapt to the environment.

Investigations into the factors that affect the lifespan of C. elegans have mostly studied hermaphrodites. For example, one recent study showed that mating shortens the lifespan of hermaphrodites. Another study showed that pheromones – hormones that change the behavior of other worms – also shorten hermaphrodite lifespan. The male pheromone is produced by males and sensed by both males and hermaphrodites. But does mating and male pheromone affect the lifespan of male roundworms?

Shi et al. have now studied Caenorhabditis worms of different species and sexes to investigate how sexual behaviors and male pheromone regulate the lifespan of male roundworms. The results of the experiments revealed two distinct mechanisms of male death. Firstly, mating caused the males of many different Caenorhabditis species to shrink and die, and also killed females and hermaphrodites. Secondly, the males of hermaphroditic species – and only these males – could also be killed by male pheromone.

The results suggest that death from mating may be an unavoidable cost of reproducing that is seen across all sexes and species of roundworm. In contrast, death by male pheromone may be a way of culling the male population in hermaphroditic species, for example, after stressful conditions have caused a sudden increase in the number of male worms.

Further work is now needed to investigate the finer details of the mechanisms by which mating and male pheromone cause death. Ultimately, this work in Caenorhabditis could be extended to help us to understand how other animals regulate their lifespan and maintain an optimum ratio of the sexes.

DOI:http://dx.doi.org/10.7554/eLife.23493.002

Cytoplasmic-Nuclear Incompatibility Between Wild Isolates of Caenorhabditis nouraguensis

How species arise is a fundamental question in biology. Species can be defined as populations of interbreeding individuals that are reproductively isolated from other such populations. Therefore, understanding how reproductive barriers evolve between populations is essential for understanding the process of speciation. Hybrid incompatibility (for example, hybrid sterility or lethality) is a common and strong reproductive barrier in nature. Here we report a lethal incompatibility between two wild isolates of the nematode Caenorhabditis nouraguensis Hybrid inviability results from the incompatibility between a maternally inherited cytoplasmic factor from each strain and a recessive nuclear locus from the other. We have excluded the possibility that maternally inherited endosymbiotic bacteria cause the incompatibility by treating both strains with tetracycline and show that hybrid death is unaffected. Furthermore, cytoplasmic-nuclear incompatibility commonly occurs between other wild isolates, indicating that this is a significant reproductive barrier within C. nouraguensis We hypothesize that the maternally inherited cytoplasmic factor is the mitochondrial genome and that mitochondrial dysfunction underlies hybrid death. This system has the potential to shed light on the dynamics of divergent mitochondrial-nuclear coevolution and its role in promoting speciation.

Host-Specific Functional Significance of Caenorhabditis Gut Commensals

The gut microbiota is an important contributor to host health and fitness. Given its importance, microbiota composition should not be left to chance. However, what determines this composition is far from clear, with results supporting contributions of both environmental factors and host genetics. To gauge the relative contributions of host genetics and environment, specifically the microbial diversity, we characterized the gut microbiotas of Caenorhabditis species spanning 200-300 million years of evolution, and raised on different composted soil environments. Comparisons were based on 16S rDNA deep sequencing data, as well as on functional evaluation of gut isolates. Worm microbiotas were distinct from those in their respective soil environment, and included bacteria previously identified as part of the C. elegans core microbiota. Microbiotas differed between experiments initiated with different soil communities, but within each experiment, worm microbiotas clustered according to host identity, demonstrating a dominant contribution of environmental diversity, but also a significant contribution of host genetics. The dominance of environmental contributions hindered identification of host-associated microbial taxa from 16S data. Characterization of gut isolates from C. elegans and C. briggsae, focusing on the core family Enterobacteriaceae, were also unable to expose phylogenetic distinctions between microbiotas of the two species. However, functional evaluation of the isolates revealed host-specific contributions, wherein gut commensals protected their own host from infection, but not a non-host. Identification of commensal host-specificity at the functional level, otherwise overlooked in standard sequence-based analyses, suggests that the contribution of host genetics to shaping of gut microbiotas may be greater than previously realized.

Chemical activation of a food deprivation signal extends lifespan

Model organisms subject to dietary restriction (DR) generally live longer. Accompanying this lifespan extension are improvements in overall health, based on multiple metrics. This indicates that pharmacological treatments that mimic the effects of DR could improve health in humans. To find new chemical structures that extend lifespan, we screened 30 000 synthetic, diverse drug‐like chemicals in Caenorhabditis elegans and identified several structurally related compounds that acted through DR mechanisms. The most potent of these NP1 impinges upon a food perception pathway by promoting glutamate signaling in the pharynx. This results in the overriding of a GPCR pathway involved in the perception of food and which normally acts to decrease glutamate signals. Our results describe the activation of a dietary restriction response through the pharmacological masking of a novel sensory pathway that signals the presence of food. This suggests that primary sensory pathways may represent novel targets for human pharmacology.

Trichinella spiralis: Adaptation and parasitism

Publication of the genome from the clade I organism, Trichinella spiralis, has provided us an avenue to address more holistic problems in parasitology; namely the processes of adaptation and the evolution of parasitism. Parasitism among nematodes has evolved in multiple, independent events. Deciphering processes that drive species diversity and adaptation are keys to understanding parasitism and advancing control strategies. Studies have been put forth on morphological and physiological aspects of parasitism and adaptation in nematodes; however, data is now coming available to investigate adaptation, host switching and parasitism at the genomic level. Herein we compare proteomic data from the clade I parasite, Trichinella spiralis with data from Brugia malayi (clade III), Meloidogyne hapla and Meloidogyne incognita (clade IV), and free-living nematodes belonging to the genera Caenorhabditis and Pristionchus (clade V). We explore changes in protein family birth/death and expansion/reduction over the course of metazoan evolution using Homo sapiens, Drosophila melanogaster and Saccharomyces cerevisiae as outgroups for the phylum Nematoda. We further examine relationships between these changes and the ability and/or result of nematodes adapting to their environments. Data are consistent with gene loss occurring in conjunction with nematode specialization resulting from parasitic worms acclimating to well-defined, environmental niches. We observed evidence for independent, lateral gene transfer events involving conserved genes that may have played a role in the evolution of nematode parasitism. In general, parasitic nematodes gained proteins through duplication and lateral gene transfer, and lost proteins through random mutation and deletions. Data suggest independent acquisition rather than ancestral inheritance among the Nematoda followed by selective gene loss over evolutionary time. Data also show that parasitism and adaptation affected a broad range of proteins, especially those involved in sensory perception, metabolism, and transcription/translation. New protein gains with functions related to regulating transcription and translation, and protein family expansions with functions related to morphology and body development have occurred in association with parasitism. Further gains occurred as a result of lateral gene transfer and in particular, with the cyanase protein family In contrast, reductions and/or losses have occurred in protein families with functions related to metabolic process and signal transduction. Taking advantage of the independent occurrences of parasitism in nematodes, which enabled us to distinguish changes associated with parasitism from species specific niche adaptation, our study provides valuable insights into nematode parasitism at a proteome level using T. spiralis as a benchmark for early adaptation to or acquisition of parasitism.

Comparative genomic analysis of upstream miRNA regulatory motifs in Caenorhabditis

MicroRNAs (miRNAs) comprise a class of short noncoding RNA molecules that play diverse developmental and physiological roles by controlling mRNA abundance and protein output of the vast majority of transcripts. Despite the importance of miRNAs in regulating gene function, we still lack a complete understanding of how miRNAs themselves are transcriptionally regulated. To fill this gap, we predicted regulatory sequences by searching for abundant short motifs located upstream of miRNAs in eight species of Caenorhabditis nematodes. We identified three conserved motifs across the Caenorhabditis phylogeny that show clear signatures of purifying selection from comparative genomics, patterns of nucleotide changes in motifs of orthologous miRNAs, and correlation between motif incidence and miRNA expression. We then validated our predictions with transgenic green fluorescent protein reporters and site-directed mutagenesis for a subset of motifs located in an enhancer region upstream of let-7 We demonstrate that a CT-dinucleotide motif is sufficient for proper expression of GFP in the seam cells of adult C. elegans, and that two other motifs play incremental roles in combination with the CT-rich motif. Thus, functional tests of sequence motifs identified through analysis of molecular evolutionary signatures provide a powerful path for efficiently characterizing the transcriptional regulation of miRNA genes.

Increased Protein Stability and Decreased Protein Turnover in the Caenorhabditis elegans Ins/IGF-1 daf-2 Mutant

In Caenorhabditis elegans, cellular proteostasis is likely essential for longevity. Autophagy has been shown to be essential for lifespan extension of daf-2 insulin/IGF mutants. Therefore, it can be hypothesized that daf-2 mutants achieve this phenotype by increasing protein turnover. However, such a mechanism would exert a substantial energy cost. By using classical ³⁵S pulse-chase labeling, we observed that protein synthesis and degradation rates are decreased in young adults of the daf-2 insulin/IGF mutants. Although reduction of protein turnover may be energetically favorable, it may lead to accumulation and aggregation of damaged proteins. As this has been shown not to be the case in daf-2 mutants, another mechanism must exist to maintain proteostasis in this strain. We observed that proteins isolated from daf-2 mutants are more soluble in acidic conditions due to increased levels of trehalose. This suggests that trehalose may decrease the potential for protein aggregation and increases proteostasis in the daf-2 mutants. We postulate that daf-2 mutants save energy by decreasing protein turnover rates and instead stabilize their proteome by trehalose.

Metabolome and proteome changes with aging in Caenorhabditis elegans

To expand the understanding of aging in the model organism Caenorhabditis elegans, global quantification of metabolite and protein levels in young and aged nematodes was performed using mass spectrometry. With age, there was a decreased abundance of proteins functioning in transcription termination, mRNA degradation, mRNA stability, protein synthesis, and proteasomal function. Furthermore, there was altered S-adenosyl methionine metabolism as well as a decreased abundance of the S-adenosyl methionine synthetase (SAMS-1) protein. Other aging-related changes included alterations in free fatty acid levels and composition, decreased levels of ribosomal proteins, decreased levels of NADP-dependent isocitrate dehydrogenase (IDH1), a shift in the cellular redox state, an increase in sorbitol content, alterations in free amino acid levels, and indications of altered muscle function and sarcoplasmic reticulum Ca(2+) homeostasis. There were also decreases in pyrimidine and purine metabolite levels, most markedly nitrogenous bases. Supplementing the culture medium with cytidine (a pyrimidine nucleoside) or hypoxanthine (a purine base) increased lifespan slightly, suggesting that aging-induced alterations in ribonucleotide metabolism affect lifespan. An age-related increase in body size, lipotoxicity from ectopic yolk lipoprotein accumulation, a decline in NAD(+) levels, and mitochondrial electron transport chain dysfunction may explain many of these changes. In addition, dietary restriction in aged worms resulting from sarcopenia of the pharyngeal pump likely decreases the abundance of SAMS-1, possibly leading to decreased phosphatidylcholine levels, larger lipid droplets, and ER and mitochondrial stress. The complementary use of proteomics and metabolomics yielded unique insights into the molecular processes altered with age in C. elegans.

Gametic selection, developmental trajectories, and extrinsic heterogeneity in Haldane's rule

Deciphering the genetic and developmental causes of the disproportionate rarity, inviability, and sterility of hybrid males, Haldane's rule, is important for understanding the evolution of reproductive isolation between species. Moreover, extrinsic and prezygotic factors can contribute to the magnitude of intrinsic isolation experienced between species with partial reproductive compatibility. Here, we use the nematodes Caenorhabditis briggsae and C. nigoni to quantify the sensitivity of hybrid male viability to extrinsic temperature and developmental timing, and test for a role of mito-nuclear incompatibility as a genetic cause. We demonstrate that hybrid male inviability manifests almost entirely as embryonic, not larval, arrest and is maximal at the lowest rearing temperatures, indicating an intrinsic-by-extrinsic interaction to hybrid inviability. Crosses using mitochondrial substitution strains that have reciprocally introgressed mitochondrial and nuclear genomes show that mito-nuclear incompatibility is not a dominant contributor to postzygotic isolation and does not drive Haldane's rule in this system. Crosses also reveal that competitive superiority of X-bearing sperm provides a novel means by which postmating prezygotic factors exacerbate the rarity of hybrid males. These findings highlight the important roles of gametic, developmental, and extrinsic factors in modulating the manifestation of Haldane's rule.

The cell survival pathways of the primordial RNA-DNA complex remain conserved in the extant genomes and may function as proto-oncogenes

Malignantly transformed (cancer) cells of multicellular hosts, including human cells, operate activated biochemical pathways that recognizably derived from unicellular ancestors. The descendant heat shock proteins of thermophile archaea now chaperon oncoproteins. The ABC cassettes of toxin-producer zooxantella Symbiodinia algae pump out the cytoplasmic toxin molecules; malignantly transformed cells utilize the derivatives of these cassettes to get rid of chemotherapeuticals. High mobility group helix-loop-helix proteins, protein arginine methyltransferases, proliferating cell nuclear antigens, and Ki-67 nuclear proteins, that protect and repair DNA in unicellular life forms, support oncogenes in transformed cells. The cell survival pathways of Wnt-β-catenin, Hedgehog, PI3K, MAPK-ERK, STAT, Ets, JAK, Pak, Myb, achaete scute, circadian rhythms, Bruton kinase and others, which are physiological in uni- and early multicellular eukaryotic life forms, are constitutively encoded in complex oncogenic pathways in selected single cells of advanced multicellular eukaryotic hosts. Oncogenes and oncoproteins in advanced multicellular hosts recreate selected independently living and immortalized unicellular life forms, which are similar to extinct and extant protists. These unicellular life forms are recognized at the clinics as autologous "cancer cells".

Mitochondrial genome of Caenorhabditis nigoni (Rhabditida: Rhabditidae)

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.

Microevolution of nematode miRNAs reveals diverse modes of selection

Micro-RNA (miRNA) genes encode abundant small regulatory RNAs that play key roles during development and in homeostasis by fine tuning and buffering gene expression. This layer of regulatory control over transcriptional networks is preserved by selection across deep evolutionary time, yet selection pressures on individual miRNA genes in contemporary populations remain poorly characterized in any organism. Here, we quantify nucleotide variability for 129 miRNAs in the genome of the nematode Caenorhabditis remanei to understand the microevolution of this important class of regulatory genes. Our analysis of three population samples and C. remanei's sister species revealed ongoing natural selection that constrains evolution of all sequence domains within miRNA hairpins. We also show that new miRNAs evolve faster than older miRNAs but that selection nevertheless favors their persistence. Despite the ongoing importance of purging of new mutations, we discover a trove of >400 natural miRNA sequence variants that include single nucleotide polymorphisms in seed motifs, indels that ablate miRNA functional domains, and origination of new miRNAs by duplication. Moreover, we demonstrate substantial nucleotide divergence of pre-miRNA hairpin alleles between populations and sister species. These findings from the first global survey of miRNA microevolution in Caenorhabditis support the idea that changes in gene expression, mediated through divergence in miRNA regulation, can contribute to phenotypic novelty and adaptation to specific environments in the present day as well as the distant past.