Selfish conflict underlies RNA-mediated parent-of-origin effects.

Genomic imprinting-the non-equivalence of maternal and paternal genomes-is a critical process that has evolved independently in many plant and mammalian species1,2. According to kinship theory, imprinting is the inevitable consequence of conflictive selective forces acting on differentially expressed parental alleles3,4. Yet, how these epigenetic differences evolve in the first place is poorly understood3,5,6. Here we report the identification and molecular dissection of a parent-of-origin effect on gene expression that might help to clarify this fundamental question. Toxin-antidote elements (TAs) are selfish elements that spread in populations by poisoning non-carrier individuals7-9. In reciprocal crosses between two Caenorhabditis tropicalis wild isolates, we found that the slow-1/grow-1 TA is specifically inactive when paternally inherited. This parent-of-origin effect stems from transcriptional repression of the slow-1 toxin by the PIWI-interacting RNA (piRNA) host defence pathway. The repression requires PIWI Argonaute and SET-32 histone methyltransferase activities and is transgenerationally inherited via small RNAs. Remarkably, when slow-1/grow-1 is maternally inherited, slow-1 repression is halted by a translation-independent role of its maternal mRNA. That is, slow-1 transcripts loaded into eggs-but not SLOW-1 protein-are necessary and sufficient to counteract piRNA-mediated repression. Our findings show that parent-of-origin effects can evolve by co-option of the piRNA pathway and hinder the spread of selfish genes that require sex for their propagation.

Natural variation in infection specificity of Caenorhabditis briggsae isolates by two RNA viruses.

Antagonistic relationships such as host-virus interactions potentially lead to rapid evolution and specificity in interactions. The Orsay virus is so far the only horizontal virus naturally infecting the nematode C. elegans. In contrast, several related RNA viruses infect its congener C. briggsae, including Santeuil (SANTV) and Le Blanc (LEBV) viruses. Here we focus on the host's intraspecific variation in sensitivity to these two intestinal viruses. Many temperate-origin C. briggsae strains, including JU1264 and JU1498, are sensitive to both, while many tropical strains, such as AF16, are resistant to both. Interestingly, some C. briggsae strains exhibit a specific resistance, such as the HK104 strain, specifically resistant to LEBV. The viral sensitivity pattern matches the strains' geographic and genomic relationships. The heavily infected strains mount a seemingly normal small RNA response that is insufficient to suppress viral infection, while the resistant strains show no small RNA response, suggesting an early block in viral entry or replication. We use a genetic approach from the host side to map genomic regions participating in viral resistance polymorphisms. Using Advanced Intercrossed Recombinant Inbred Lines (RILs) between virus-resistant AF16 and SANTV-sensitive HK104, we detect Quantitative Trait Loci (QTLs) on chromosomes IV and III. Building RILs between virus-sensitive JU1498 and LEBV-resistant HK104 followed by bulk segregant analysis, we identify a chromosome II QTL. In both cases, further introgressions of the regions confirmed the QTLs. This diversity provides an avenue for studying virus entry, replication, and exit mechanisms, as well as host-virus specificity and the host response to a specific virus infection.

Structural basis of the 9-fold symmetry of centrioles.

The centriole, and the related basal body, is an ancient organelle characterized by a universal 9-fold radial symmetry and is critical for generating cilia, flagella, and centrosomes. The mechanisms directing centriole formation are incompletely understood and represent a fundamental open question in biology. Here, we demonstrate that the centriolar protein SAS-6 forms rod-shaped homodimers that interact through their N-terminal domains to form oligomers. We establish that such oligomerization is essential for centriole formation in C. elegans and human cells. We further generate a structural model of the related protein Bld12p from C. reinhardtii, in which nine homodimers assemble into a ring from which nine coiled-coil rods radiate outward. Moreover, we demonstrate that recombinant Bld12p self-assembles into structures akin to the central hub of the cartwheel, which serves as a scaffold for centriole formation. Overall, our findings establish a structural basis for the universal 9-fold symmetry of centrioles.

CaeNDR, the Caenorhabditis Natural Diversity Resource.

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.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes.

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Rewiring the Sex-Determination Pathway During the Evolution of Self-Fertility.

Although evolution is driven by changes in how regulatory pathways control development, we know little about the molecular details underlying these transitions. The TRA-2 domain that mediates contact with TRA-1 is conserved in Caenorhabditis. By comparing the interaction of these proteins in two species, we identified a striking change in how sexual development is controlled. Identical mutations in this domain promote oogenesis in Caenorhabditis elegans but promote spermatogenesis in Caenorhabditis briggsae. Furthermore, the effects of these mutations involve the male-promoting gene fem-3 in C. elegans but are independent of fem-3 in C. briggsae. Finally, reciprocal mutations in these genes show that C. briggsae TRA-2 binds TRA-1 to prevent expression of spermatogenesis regulators. By contrast, in C. elegans TRA-1 sequesters TRA-2 in the germ line, allowing FEM-3 to initiate spermatogenesis. Thus, we propose that the flow of information within the sex determination pathway has switched directions during evolution. This result has important implications for how evolutionary change can occur.

Genetic exchange with an outcrossing sister species causes severe genome-wide dysregulation in a selfing Caenorhabditis nematode.

Different modes of reproduction evolve rapidly, with important consequences for genome composition. Selfing species often occupy a similar niche as their outcrossing sister species with which they are able to mate and produce viable hybrid progeny, raising the question of how they maintain genomic identity. Here, we investigate this issue by using the nematode Caenorhabditis briggsae, which reproduces as a hermaphrodite, and its outcrossing sister species Caenorhabditis nigoni We hypothesize that selfing species might develop some barriers to prevent gene intrusions through gene regulation. We therefore examined gene regulation in the hybrid F2 embryos resulting from reciprocal backcrosses between F1 hybrid progeny and C. nigoni or C. briggsae F2 hybrid embryos with ∼75% of their genome derived from C. briggsae (termed as bB2) were inviable, whereas those with ∼75% of their genome derived from C. nigoni (termed as nB2) were viable. Misregulation of transposable elements, coding genes, and small regulatory RNAs was more widespread in the bB2 compared with the nB2 hybrids, which is a plausible explanation for the differential phenotypes between the two hybrids. Our results show that regulation of the C. briggsae genome is strongly affected by genetic exchanges with its outcrossing sister species, C. nigoni, whereas regulation of the C. nigoni genome is more robust on genetic exchange with C. briggsae The results provide new insights into how selfing species might maintain their identity despite genetic exchanges with closely related outcrossing species.

The GATA factor ELT-3 specifies endoderm in Caenorhabditis angaria in an ancestral gene network.

Endoderm specification in Caenorhabditis elegans occurs through a network in which maternally provided SKN-1/Nrf, with additional input from POP-1/TCF, activates the GATA factor cascade MED-1,2→END-1,3→ELT-2,7. Orthologues of the MED, END and ELT-7 factors are found only among nematodes closely related to C. elegans, raising the question of how gut is specified in their absence in more distant species in the genus. We find that the C. angaria, C. portoensis and C. monodelphis orthologues of the GATA factor gene elt-3 are expressed in the early E lineage, just before their elt-2 orthologues. In C. angaria, Can-pop-1(RNAi), Can-elt-3(RNAi) and a Can-elt-3 null mutation result in a penetrant 'gutless' phenotype. Can-pop-1 is necessary for Can-elt-3 activation, showing that it acts upstream. Forced early E lineage expression of Can-elt-3 in C. elegans can direct the expression of a Can-elt-2 transgene and rescue an elt-7 end-1 end-3; elt-2 quadruple mutant strain to viability. Our results demonstrate an ancestral mechanism for gut specification and differentiation in Caenorhabditis involving a simpler POP-1→ELT-3→ELT-2 gene network.

Increasing Complexity of the N-Glycome During Caenorhabditis Development.

Caenorhabditis elegans is a frequently employed genetic model organism and has been the object of a wide range of developmental, genetic, proteomic, and glycomic studies. Here, using an off-line MALDI-TOF-MS approach, we have analyzed the N-glycans of mixed embryos and liquid- or plate-grown L4 larvae. Of the over 200 different annotatable N-glycan structures, variations between the stages as well as the mode of cultivation were observed. While the embryonal N-glycome appears less complicated overall, the liquid- and plate-grown larvae differ especially in terms of methylation of bisecting fucose, α-galactosylation of mannose, and di-ß-galactosylation of core α1,6-fucose. Furthermore, we analyzed the O-glycans by LC-electrospray ionization-MS following ß-elimination; especially the embryonal O-glycomes included a set of phosphorylcholine-modified structures, previously not shown to exist in nematodes. However, the set of glycan structures cannot be clearly correlated with levels of glycosyltransferase transcripts in developmental RNA-Seq datasets, but there is an indication for coordinated expression of clusters of potential glycosylation-relevant genes. Thus, there are still questions to be answered in terms of how and why a simple nematode synthesizes such a diverse glycome.

The Caenorhabditis elegans ARIP-4 DNA helicase couples mitochondrial surveillance to immune, detoxification, and antiviral pathways.

Surveillance of Caenorhabditis elegans mitochondrial status is coupled to defense responses such as drug detoxification, immunity, antiviral RNA interference (RNAi), and regulation of life span. A cytochrome p540 detoxification gene, cyp-14A4, is specifically activated by mitochondrial dysfunction. The nuclear hormone receptor NHR-45 and the transcriptional Mediator component MDT-15/MED15 are required for the transcriptional activation of cyp-14A4 by mitochondrial mutations, gene inactivations, or toxins. A genetic screen for mutations that fail to activate this cytochrome p450 gene upon drug or mutation-induced mitochondrial dysfunction identified a DNA helicase ARIP-4 that functions in concert with the NHR-45 transcriptional regulatory cascade. In response to mitochondrial dysfunction, ARIP-4 and NHR-45 protein interaction is enhanced, and they relocalize from the nuclear periphery to the interior of intestinal nuclei. NHR-45/ARIP-4 also regulates the transcriptional activation of the eol-1 gene that encodes a decapping enzyme required for enhanced RNAi and transgene silencing of mitochondrial mutants. In the absence of arip-4, animals were more susceptible to the mitochondrial inhibitor antimycin. Thus, ARIP-4 serves as a transcriptional coactivator of NHR-45 to promote this defense response. A null mutation in arip-4 extends the life span and health span of both wild type and a mitochondrial mutant, suggesting that the activation of detoxification pathways is deleterious to health when the mitochondrial dysfunction is caused by mutation that cannot be cytochrome p450-detoxified. Thus, arip-4 acts in a pathway that couples mitochondrial surveillance to the activation of downstream immunity, detoxification, and RNAi responses.

Genome Size Changes by Duplication, Divergence, and Insertion in Caenorhabditis Worms.

Genome size has been measurable since the 1940s but we still do not understand genome size variation. Caenorhabditis nematodes show strong conservation of chromosome number but vary in genome size between closely related species. Androdioecy, where populations are composed of males and self-fertile hermaphrodites, evolved from outcrossing, female-male dioecy, three times in this group. In Caenorhabditis, androdioecious genomes are 10-30% smaller than dioecious species, but in the nematode Pristionchus, androdioecy evolved six times and does not correlate with genome size. Previous hypotheses include genome size evolution through: 1) Deletions and "genome shrinkage" in androdioecious species; 2) Transposable element (TE) expansion and DNA loss through large deletions (the "accordion model"); and 3) Differing TE dynamics in androdioecious and dioecious species. We analyzed nematode genomes and found no evidence for these hypotheses. Instead, nematode genome sizes had strong phylogenetic inertia with increases in a few dioecious species, contradicting the "genome shrinkage" hypothesis. TEs did not explain genome size variation with the exception of the DNA transposon Mutator which was twice as abundant in dioecious genomes. Across short and long evolutionary distances Caenorhabditis genomes evolved through small structural mutations including gene-associated duplications and insertions. Seventy-one protein families had significant, parallel decreases across androdioecious Caenorhabditis including genes involved in the sensory system, regulatory proteins and membrane-associated immune responses. Our results suggest that within a dynamic landscape of frequent small rearrangements in Caenorhabditis, reproductive mode mediates genome evolution by altering the precise fates of individual genes, proteins, and the phenotypes they underlie.

Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution.

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a "faster male" model for Haldane's rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a "large-X effect" for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process.

Efficient targeted recombination with CRISPR/Cas9 in hybrids of Caenorhabditis nematodes with suppressed recombination.

BACKGROUND: Homology-based recombination (HR) is the cornerstone of genetic mapping. However, a lack of sufficient sequence homology or the presence of a genomic rearrangement prevents HR through crossing, which inhibits genetic mapping in relevant genomic regions. This is particularly true in species hybrids whose genomic sequences are highly divergent along with various genome arrangements, making the mapping of genetic loci, such as hybrid incompatibility (HI) loci, through crossing impractical. We previously mapped tens of HI loci between two nematodes, Caenorhabditis briggsae and C. nigoni, through the repeated backcrossing of GFP-linked C. briggsae fragments into C. nigoni. However, the median introgression size was over 7 Mb, indicating apparent HR suppression and preventing the subsequent cloning of the causative gene underlying a given HI phenotype. Therefore, a robust method that permits recombination independent of sequence homology is desperately desired. RESULTS: Here, we report a method of highly efficient targeted recombination (TR) induced by CRISPR/Cas9 with dual guide RNAs (gRNAs), which circumvents the HR suppression in hybrids between the two species. We demonstrated that a single gRNA was able to induce efficient TR between highly homologous sequences only in the F1 hybrids but not in the hybrids that carry a GFP-linked C. briggsae fragment in an otherwise C. nigoni background. We achieved highly efficient TR, regardless of sequence homology or genetic background, when dual gRNAs were used that each specifically targeted one parental chromosome. We further showed that dual gRNAs were able to induce efficient TR within genomic regions that had undergone inversion, in which HR-based recombination was expected to be suppressed, supporting the idea that dual-gRNA-induced TR can be achieved through nonhomology-based end joining between two parental chromosomes. CONCLUSIONS: Recombination suppression can be circumvented through CRISPR/Cas9 with dual gRNAs, regardless of sequence homology or the genetic background of the species hybrid. This method is expected to be applicable to other situations in which recombination is suppressed in interspecies or intrapopulation hybrids.

Novel and improved Caenorhabditis briggsae gene models generated by community curation.

BACKGROUND: The nematode Caenorhabditis briggsae has been used as a model in comparative genomics studies with Caenorhabditis elegans because of their striking morphological and behavioral similarities. However, the potential of C. briggsae for comparative studies is limited by the quality of its genome resources. The genome resources for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans. The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of software-derived gene predictions that contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 gene models and underlying transcriptomic data to repair software-derived errors. RESULTS: We designed a detailed workflow to train a team of nine students to manually curate gene models using RNA read alignments. We manually inspected the gene models, proposed corrections to the coding sequences of over 8,000 genes, and modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality and showed that manual curation led to substantial improvements in the protein sequence length accuracy of QX1410 genes. Additionally, collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome. CONCLUSIONS: Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. Our manual curation efforts have brought the QX1410 gene models to a comparable level of quality as the extensively curated AF16 gene models. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

Temperature-dependent Small RNA Expression Depends on Wild Genetic Backgrounds of Caenorhabditis briggsae.

Geographically distinct populations can adapt to the temperature conditions of their local environment, leading to temperature-dependent fitness differences between populations. Consistent with local adaptation, phylogeographically distinct Caenorhabditis briggsae nematodes show distinct fitness responses to temperature. The genetic mechanisms underlying local adaptation, however, remain unresolved. To investigate the potential role of small noncoding RNAs in genotype-specific responses to temperature, we quantified small RNA expression using high-throughput sequencing of C. briggsae nematodes from tropical and temperate strain genotypes reared under three temperature conditions (14 °C, 20 °C, and 30 C). Strains representing both tropical and temperate regions showed significantly lower expression of PIWI-interacting RNAs (piRNAs) at high temperatures, primarily mapping to a large ∼7 Mb long piRNA cluster on chromosome IV. We also documented decreased expression of 22G-RNAs antisense to protein-coding genes and other genomic features at high rearing temperatures for the thermally-intolerant temperate strain genotype, but not for the tropical strain genotype. Reduced 22G-RNA expression was widespread along chromosomes and among feature types, indicative of a genome-wide response. Targets of the EGO-1/CSR-1 22G-RNA pathway were most strongly impacted compared with other 22G-RNA pathways, implicating the CSR-1 Argonaute and its RNA-dependent RNA polymerase EGO-1 in the genotype-dependent modulation of C. briggsae 22G-RNAs under chronic thermal stress. Our work suggests that gene regulation via small RNAs may be an important contributor to the evolution of local adaptations.

Recurrent but Short-Lived Duplications of Centromeric Proteins in Holocentric Caenorhabditis Species.

Centromeric histones (CenH3s) are essential for chromosome inheritance during cell division in most eukaryotes. CenH3 genes have rapidly evolved and undergone repeated gene duplications and diversification in many plant and animal species. In Caenorhabditis species, two independent duplications of CenH3 (named hcp-3 for HoloCentric chromosome-binding Protein 3) were previously identified in C. elegans and C. remanei. Using phylogenomic analyses in 32 Caenorhabditis species, we find strict retention of the ancestral hcp-3 gene and 10 independent duplications. Most hcp-3L (hcp-3-like) paralogs are only found in 1-2 species, are expressed in both males and females/hermaphrodites, and encode histone fold domains with 69-100% identity to ancestral hcp-3. We identified novel N-terminal protein motifs, including putative kinetochore protein-interacting motifs and a potential separase cleavage site, which are well conserved across Caenorhabditis HCP-3 proteins. Other N-terminal motifs vary in their retention across paralogs or species, revealing potential subfunctionalization or functional loss following duplication. An N-terminal extension in the hcp-3L gene of C. afra revealed an unprecedented protein fusion, where hcp-3L fused to duplicated segments from hcp-4 (nematode CENP-C). By extending our analyses beyond CenH3, we found gene duplications of six inner and outer kinetochore genes in Caenorhabditis, which appear to have been retained independent of hcp-3 duplications. Our findings suggest that centromeric protein duplications occur frequently in Caenorhabditis nematodes, are selectively retained for short evolutionary periods, then degenerate or are lost entirely. We hypothesize that unique challenges associated with holocentricity in Caenorhabditis may lead to this rapid "revolving door" of kinetochore protein paralogs.

The compact genome of Caenorhabditis niphades n. sp., isolated from a wood-boring weevil, Niphades variegatus.

BACKGROUND: The first metazoan genome sequenced, that of Caenorhabditis elegans, has motivated animal genome evolution studies. To date > 50 species from the genus Caenorhabditis have been sequenced, allowing research on genome variation. RESULTS: In the present study, we describe a new gonochoristic species, Caenorhabditis niphades n. sp., previously referred as C. sp. 36, isolated from adult weevils (Niphades variegatus), with whom they appear to be tightly associated during its life cycle. Along with a species description, we sequenced the genome of C. niphades n. sp. and produced a chromosome-level assembly. A genome comparison highlighted that C. niphades n. sp. has the smallest genome (59 Mbp) so far sequenced in the Elegans supergroup, despite being closely related to a species with an exceptionally large genome, C. japonica. CONCLUSIONS: The compact genome of C. niphades n. sp. can serve as a key resource for comparative evolutionary studies of genome and gene number expansions in Caenorhabditis species.

A role for worm cutl-24 in background- and parent-of-origin-dependent ER stress resistance.

BACKGROUND: Organisms in the wild can acquire disease- and stress-resistance traits that outstrip the programs endogenous to humans. Finding the molecular basis of such natural resistance characters is a key goal of evolutionary genetics. Standard statistical-genetic methods toward this end can perform poorly in organismal systems that lack high rates of meiotic recombination, like Caenorhabditis worms. RESULTS: Here we discovered unique ER stress resistance in a wild Kenyan C. elegans isolate, which in inter-strain crosses was passed by hermaphrodite mothers to hybrid offspring. We developed an unbiased version of the reciprocal hemizygosity test, RH-seq, to explore the genetics of this parent-of-origin-dependent phenotype. Among top-scoring gene candidates from a partial-coverage RH-seq screen, we focused on the neuronally-expressed, cuticlin-like gene cutl-24 for validation. In gene-disruption and controlled crossing experiments, we found that cutl-24 was required in Kenyan hermaphrodite mothers for ER stress tolerance in their inter-strain hybrid offspring; cutl-24 was also a contributor to the trait in purebred backgrounds. CONCLUSIONS: These data establish the Kenyan strain allele of cutl-24 as a determinant of a natural stress-resistant state, and they set a precedent for the dissection of natural trait diversity in invertebrate animals without the need for a panel of meiotic recombinants.

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.

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands.

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three-year period, measured various environmental parameters at each sampling site, and whole-genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500-1,500 m) where ambient air temperatures are cool (15-20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper-divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.