Biology and genome of a newly discovered sibling species of Caenorhabditis elegans

Bioinformatic Annotation of Transposon DNA Processing Genes on the Long-Read Genome Assembly of Caenorhabditis elegans

Transposable elements (TEs) or transposons are thought to play roles in animal physiological processes, such as germline, early embryonic, and brain development, as well as aging. However, their roles have not been systematically investigated through experimental studies. In this study, we created a catalog of genes directly involved in replication, excision, or integration of transposon-coding DNA, which we refer to as transposon DNA processing genes (TDPGs). Specifically, to bridge the gap to experimental studies, we sought potentially functional TDPGs which maintain intact open reading frames and the amino acids at their catalytic cores on the latest long-read genome assembly of Caenorhabditis elegans, VC2010. Among 52 519 TE loci, we identified 145 potentially functional TDPGs encoded in long terminal repeat elements, long interspersed nuclear elements, terminal inverted repeat elements, Helitrons, and Mavericks/Polintons. Our TDPG catalog, which contains a feasible number of genes, allows for the experimental manipulation of TE mobility in vivo, regardless of whether the TEs are autonomous or non-autonomous, thereby potentially promoting the study of the physiological functions of TE mobility.

A ubiquinone precursor analogue does not clearly increase the growth rate of Caenorhabditis inopinata

The evolution of developmental rates may drive morphological change. Caenorhabditis inopinata develops nearly twice as slowly as Caenorhabditis elegans . clk-1 encodes a hydroxylase required for synthesizing ubiquinone, and mutant clk-1 slow growth phenotypes can be rescued by supplying animals with a ubiquinone precursor analogue, 2,4-dihydroxybenzoate. RNA-seq data showing low clk-1 expression raised the possibility that C. inopinata grows slowly because of reduced ubiquinone biosynthesis. C. inopinata did not reveal a clear reduction in the age of maturation when reared on 2,4-dihydroxybenzoate. Further scrutiny of RNA-seq results revealed multiple ubiquinone metabolism genes have low expression in C. inopinata . Divergent clk-1 expression alone may not be a major driver of the evolution of slow development in this species.

The Bacteria of a Fig Microcommunity

Understanding the biotic drivers of diversity is a major goal of microbial ecology. One approach towards tackling this issue is to interrogate relatively simple communities that are easy to observe and perturb. Figs (syconia) of the genus Ficus represent such a system. Here, we describe the microbial communities of Ficus septica figs, which are associated with the nematode Caenorhabditis inopinata (the sister species of the C. elegans genetic model system). In 2019, 38 Ficus septica figs (across 12 plants in Taiwan) were dissected, and metadata such as foundress wasp number and nematode occupancy were collected for each fig. Suspensions derived from interior fig material and fig surface washes were prepared for 16S microbial metabarcoding. Over 3,000 OTUs were detected, and microbial communities were dominated by members of Proteobacteria , Bacteroidota , and Actinobacteriota . Although microbial communities of fig exteriors and interiors can be distinguished, levels of microbial alpha diversity were comparable across these areas of the fig. Nematodes likewise had no detectable impact on microbial alpha diversity, although nematodes were associated with a modest change in microbial community composition. A handful of OTUs (associated with the genera Kosokonia , Ochobactrum , and Stenotrophomonas ) revealed potential differential abundance among figs varying in nematode occupancy. Additionally, foundress wasp number was negatively correlated with microbial alpha diversity. These findings set the stage for future studies that directly test the role of nematode and wasp occupancy on microbial communities, as well as investigations that probe nematode-microbe interactions through laboratory experiments. Taken together, these results constitute a fundamental step in characterizing the natural microbial communities of figs and Caenorhabditis nematodes.

Importance

Unraveling why different species live in different places is a longstanding open question in ecology, and it is clear that interspecific interactions among species are a major contributor to species distributions. Ficus figs are a useful system for ecological studies because they are relatively simple microcosms where characterizing animal community composition of multiple samples is straightforward. Additionally, Caenorhabditis inopinata , a close relative of the C. elegans genetic model system, thrives in Ficus septica figs. Here, we tie 16S microbial metabarcoding to nematode and wasp occupancy data to understand the causes of bacterial community composition in F. septica figs. We found that microbial composition, but not total diversity, varies among fig surface and interiors. Likewise, we found that nematode occupancy impacts microbial composition but not alpha diversity. Moreover, we show that as the number of foundress wasps increases, the microbial alpha diversity decreases. Finally, we identified OTUs that are potentially associated with nematode occupancy. Taken together, these results represent a key step in describing a microcommunity wherein ecological genetic hypotheses can be tested, as well as one that can potentially reveal the roles of uncharacterized genes in established model systems.

Life history in Caenorhabditis elegans: from molecular genetics to evolutionary ecology

Life history is defined by traits that reflect key components of fitness, especially those relating to reproduction and survival. Research in life history seeks to unravel the relationships among these traits and understand how life history strategies evolve to maximize fitness. As such, life history research integrates the study of the genetic and developmental mechanisms underlying trait determination with the evolutionary and ecological context of Darwinian fitness. As a leading model organism for molecular and developmental genetics, Caenorhabditis elegans is unmatched in the characterization of life history-related processes, including developmental timing and plasticity, reproductive behaviors, sex determination, stress tolerance, and aging. Building on recent studies of natural populations and ecology, the combination of C. elegans' historical research strengths with new insights into trait variation now positions it as a uniquely valuable model for life history research. In this review, we summarize the contributions of C. elegans and related species to life history and its evolution. We begin by reviewing the key characteristics of C. elegans life history, with an emphasis on its distinctive reproductive strategies and notable life cycle plasticity. Next, we explore intraspecific variation in life history traits and its underlying genetic architecture. Finally, we provide an overview of how C. elegans has guided research on major life history transitions both within the genus Caenorhabditis and across the broader phylum Nematoda. While C. elegans is relatively new to life history research, significant progress has been made by leveraging its distinctive biological traits, establishing it as a highly cross-disciplinary system for life history studies.

Visualizing genomic evolution in Caenorhabditis through WormSynteny

Understanding the syntenic relationships among genomes is crucial to elucidate the genomic mechanisms that drive the evolution of species. The nematode Caenorhabditis is a good model for studying genomic evolution due to the well-established biology of Caenorhabditis elegans and the availability of > 50 genomes in the genus. However, effective alignment of more than ten species in Caenorhabditis has not been conducted before, and there is currently no tool to visualize the synteny of more than two species. In this study, we used Progressive Cactus, a recently developed multigenome aligner, to align the genomes of eleven Caenorhabditis species. Through the progressive alignment, we reconstructed nine ancestral genomes, analyzed the mutational types that cause genomic rearrangement during speciation, and found that insertion and duplication are the major driving forces for genome expansion. Dioecious species appear to expand their genomes more than androdioecious species. We then built an online interactive app called WormSynteny to visualize the syntenic relationship among the eleven species. Users can search the alignment dataset using C. elegans query sequences, construct synteny plots at different genomic scales, and use a set of options to control alignment output and plot presentation. We showcased the use of WormSynteny to visualize the syntenic conservation of one-to-one orthologues among species, tandem and dispersed gene duplication in C. elegans, and the evolution of exon and intron structures. Importantly, the integration of orthogroup information with synteny linkage in WormSynteny allows the easy visualization of conserved genomic blocks and disruptive rearrangement. In conclusion, WormSynteny provides immediate access to the syntenic relationships among the most widely used Caenorhabditis species and can facilitate numerous comparative genomics studies. This pilot study with eleven species also serves as a proof-of-concept to a more comprehensive larger-scale analysis using hundreds of nematode genomes, which is expected to reveal mechanisms that drive genomic evolution in the Nematoda phylum. Finally, the WormSynteny software provides a generalizable solution for visualizing the output of Progressive Cactus with interactive graphics, which would be useful for a broad community of genome researchers.

Global analysis of neuropeptide receptor conservation across phylum Nematoda

BACKGROUND

The phylum Nematoda is incredibly diverse and includes many parasites of humans, livestock, and plants. Peptide-activated G protein-coupled receptors (GPCRs) are central to the regulation of physiology and numerous behaviors, and they represent appealing pharmacological targets for parasite control. Efforts are ongoing to characterize the functions and define the ligands of nematode GPCRs, with already most peptide GPCRs known or predicted in Caenorhabditis elegans. However, comparative analyses of peptide GPCR conservation between C. elegans and other nematode species are limited, and many nematode GPCRs remain orphan. A phylum-wide perspective on peptide GPCR profiles will benefit functional and applied studies of nematode peptide GPCRs.

RESULTS

We constructed a pan-phylum resource of C. elegans peptide GPCR orthologs in 125 nematode species using a semi-automated pipeline for analysis of predicted proteome datasets. The peptide GPCR profile varies between nematode species of different phylogenetic clades and multiple C. elegans peptide GPCRs have orthologs across the phylum Nematoda. We identified peptide ligands for two highly conserved orphan receptors, NPR-9 and NPR-16, that belong to the bilaterian galanin/allatostatin A (Gal/AstA) and somatostatin/allatostatin C (SST/AstC) receptor families. The AstA-like NLP-1 peptides activate NPR-9 in cultured cells and are cognate ligands of this receptor in vivo. In addition, we discovered an AstC-type peptide, NLP-99, that activates the AstC-type receptor NPR-16. In our pan-phylum resource, the phylum-wide representation of NPR-9 and NPR-16 resembles that of their cognate ligands more than those of allatostatin-like peptides that do not activate these receptors.

CONCLUSIONS

The repertoire of C. elegans peptide GPCR orthologs varies across phylogenetic clades and several peptide GPCRs show broad conservation in the phylum Nematoda. Our work functionally characterizes the conserved receptors NPR-9 and NPR-16 as the respective GPCRs for the AstA-like NLP-1 peptides and the AstC-related peptide NLP-99. NLP-1 and NLP-99 are widely conserved in nematodes and their representation matches that of their receptor in most species. These findings demonstrate the conservation of a functional Gal/AstA and SST/AstC signaling system in nematodes. Our dataset of C. elegans peptide GPCR orthologs also lays a foundation for further functional studies of peptide GPCRs in the widely diverse nematode phylum.

Transposon-mediated genic rearrangements underlie variation in small RNA pathways

Transposable elements (TEs) can alter host gene structure and expression, whereas host organisms develop mechanisms to repress TE activities. In the nematode Caenorhabditis elegans, a small interfering RNA pathway dependent on the helicase ERI-6/7 primarily silences retrotransposons and recent genes of likely viral origin. By studying gene expression variation among wild C. elegans strains, we found that structural variants and transposon remnants likely underlie expression variation in eri-6/7 and the pathway targets. We further found that multiple insertions of the DNA transposons, Polintons, reshuffled the eri-6/7 locus and induced inversion of eri-6 in some wild strains. In the inverted configuration, gene function was previously shown to be repaired by unusual trans-splicing mediated by direct repeats. We identified that these direct repeats originated from terminal inverted repeats of Polintons. Our findings highlight the role of host-transposon interactions in driving rapid host genome diversification among natural populations and shed light on evolutionary novelty in genes and splicing mechanisms.

Regulatory logic and transposable element dynamics in nematode worm genomes

Genome sequencing has revealed a tremendous diversity of transposable elements (TEs) in eukaryotes but there is little understanding of the evolutionary processes responsible for TE diversity. Non-autonomous TEs have lost the machinery necessary for transposition and rely on closely related autonomous TEs for critical proteins. We studied two mathematical models of TE regulation, one assuming that both autonomous tranposons and their non-autonomous relatives operate under the same regulatory logic, competing for transposition resources, and one assuming that autonomous TEs self-attenuate transposition while non-autonomous transposons continually increase, parasitizing their autonomous relatives. We implemented these models in stochastic simulations and studied how TE regulatory relationships influence transposons and populations. We found that only outcrossing populations evolving with Parasitic TE regulation resulted in stable maintenance of TEs. We tested our model predictions in Caenorhabditis genomes by annotating TEs in two focal families, autonomous LINEs and their non-autonomous SINE relatives and the DNA transposon Mutator. We found broad variation in autonomous - non-autonomous relationships and rapid mutational decay in the sequences that allow non-autonomous TEs to transpose. Together, our results suggest that individual TE families evolve according to disparate regulatory rules that are relevant in the early, acute stages of TE invasion.

Evolutionary changes of noncoding elements associated with transition of sexual mode in Caenorhabditis nematodes

The transition of the sexual mode occurs widely in animal evolution. In Caenorhabditis nematodes, androdioecy, a sexual polymorphism composed of males and hermaphrodites having the ability to self-fertilize, has evolved independently multiple times. While the modification of noncoding regulatory elements likely contributed to the evolution of hermaphroditism, little is known about these changes. Here, we conducted a genome-wide analysis of conserved noncoding elements (CNEs) focusing on the evolution of hermaphroditism in Caenorhabditis nematodes. We found that, in androdioecious nematodes, mutations rapidly accumulated in CNEs' neighboring genes associated with sexual traits. Expression analysis indicate that the identified CNEs are involved in spermatogenesis in hermaphrodites and associated with the transition of gene expression from dioecious to androdioecious nematodes. Last, genome editing of a CNE neighboring laf-1 resulted in a change in its expression in the gonadal region undergoing spermatogenesis. Our bioinformatic and experimental analyses highlight the importance of CNEs in gene regulation associated with the development of hermaphrodites.

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 Caenorhabditis 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 orthologues with divergent dynamics across this developmental period between the 2 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 transforming growth factor β signaling revealed idiosyncratic and unexpected transcriptional patterns given their role in body size regulation in C. elegans. This widespread transcriptional divergence between these species is unexpected and maybe 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.

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

The impact of differential transposition activities of autonomous and nonautonomous hAT transposable elements on genome architecture and gene expression in Caenorhabditis inopinata

Transposable elements are DNA sequences capable of moving within genomes and significantly influence genomic evolution. The nematode Caenorhabditis inopinata exhibits a much higher transposable element copy number than its sister species, Caenorhabditis elegans. In this study, we identified a novel autonomous transposable element belonging to the hAT superfamily from a spontaneous transposable element-insertion mutant in C. inopinata and named this transposon Ci-hAT1. Further bioinformatic analyses uncovered 3 additional autonomous hAT elements-Ci-hAT2, Ci-hAT3, and Ci-hAT4-along with over 1,000 copies of 2 nonautonomous miniature inverted-repeat transposable elements, mCi-hAT1 and mCi-hAT4, likely derived from Ci-hAT1 and Ci-hAT4 through internal deletion. We tracked at least 3 sequential transpositions of Ci-hAT1 over several years. However, the transposition rates of the other 3 autonomous hAT elements were lower, suggesting varying activity levels. Notably, the distribution patterns of the 2 miniature inverted-repeat transposable element families differed significantly: mCi-hAT1 was primarily located in the chromosome arms, a pattern observed in the transposable elements of other Caenorhabditis species, whereas mCi-hAT4 was more evenly distributed across chromosomes. Additionally, interspecific transcriptome analysis indicated that C. inopinata genes with upstream or intronic these miniature inverted-repeat transposable element insertions tend to be more highly expressed than their orthologous genes in C. elegans. These findings highlight the significant role of de-silenced transposable elements in driving the evolution of genomes and transcriptomes, leading to species-specific genetic diversity.

Behavioral evolution: No sex please, we're hermaphrodites

Many 'hard-wired', innate animal behaviors are related to reproduction. So what happens when reproductive systems evolve? New research in nematodes has identified principles underlying the co-evolution of reproductive strategy and sexual behavior, revealing some surprises and raising intriguing new questions.

Evolution remodels olfactory and mating-receptive behaviors in the transition from female to hermaphrodite reproduction

Male/hermaphrodite species have arisen multiple times from a male/female ancestral state in nematodes, providing a model to study behavioral adaptations to different reproductive strategies. Here, we examined the mating behaviors of male/female (gonochoristic) Caenorhabditis species in comparison with male/hermaphrodite (androdiecious) close relatives. We find that females from two species in the Elegans group chemotax to volatile odor from males, but hermaphrodites do not. Females, but not hermaphrodites, also display known mating-receptive behaviors such as sedation when male reproductive structures contact the vulva. Focusing on the male/female species C. nigoni, we show that female chemotaxis to males is limited to adult females approaching adult or near-adult males and relies upon the AWA neuron-specific transcription factor ODR-7, as does male chemotaxis to female odor as previously shown in C. elegans. However, female receptivity during mating contact is odr-7 independent. All C. nigoni female behaviors are suppressed by mating and all are absent in young hermaphrodites from the sister species C. briggsae. However, latent receptivity during mating contact can be uncovered in mutant or aged C. briggsae hermaphrodites that lack self-sperm. These results reveal two mechanistically distinct components of the shift from female to hermaphrodite behavior: the loss of female-specific odr-7-dependent chemotaxis and a sperm-dependent state of reduced receptivity to mating contact. Hermaphrodites from a second androdioecious species, C. tropicalis, recover all female behaviors upon aging, including chemotaxis to males. Regaining mating receptivity after sperm depletion could maximize hermaphrodite fitness across their lifespan.

Fantastic beasts and how to study them: rethinking experimental animal behavior

Humans have been trying to understand animal behavior at least since recorded history. Recent rapid development of new technologies has allowed us to make significant progress in understanding the physiological and molecular mechanisms underlying behavior, a key goal of neuroethology. However, there is a tradeoff when studying animal behavior and its underlying biological mechanisms: common behavior protocols in the laboratory are designed to be replicable and controlled, but they often fail to encompass the variability and breadth of natural behavior. This Commentary proposes a framework of 10 key questions that aim to guide researchers in incorporating a rich natural context into their experimental design or in choosing a new animal study system. The 10 questions cover overarching experimental considerations that can provide a template for interspecies comparisons, enable us to develop studies in new model organisms and unlock new experiments in our quest to understand behavior.

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

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.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

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.

Evolution of fertilization ability in obligatorily outcrossing populations of Caenorhabditis elegans

In species reproducing by selfing, the traits connected with outcrossing typically undergo degeneration, a phenomenon called selfing syndrome. In Caenorhabditis elegans nematodes, selfing syndrome affects many traits involved in mating, rendering cross-fertilization highly inefficient. In this study, we investigated the evolution of cross-fertilization efficiency in populations genetically modified to reproduce by obligatory outcrossing. Following the genetic modification, replicate obligatorily outcrossing were maintained for over 100 generations, at either optimal (20 °C) or elevated (24 °C) temperatures, as a part of a broader experimental evolution program. Subsequently, fertilization rates were assayed in the evolving populations, as well as their ancestors who had the obligatory outcrossing introduced but did not go through experimental evolution. Fertilization effectivity was measured by tracking the fractions of fertilized females in age-synchronized populations, through 8 h since reaching adulthood. In order to check the robustness of our measurements, each evolving population was assayed in two or three independent replicate blocks. Indeed, we found high levels of among-block variability in the fertilization trajectories, and in the estimates of divergence between evolving populations and their ancestors. We also identified five populations which appear to have evolved increased fertilization efficiency, relative to their ancestors. However, due to the abovementioned high variability, this set of populations should be treated as candidate, with further replications needed to either confirm or disprove their divergence from ancestors. Furthermore, we also discuss additional observations we have made concerning fertilization trajectories.

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.

daf-42 is an evolutionarily young gene essential for dauer development in Caenorhabditis elegans

Under adverse environmental conditions, nematodes arrest into dauer, an alternative developmental stage for diapause. Dauer endures unfavorable environments and interacts with host animals to access favorable environments, thus playing a critical role in survival. Here, we report that in Caenorhabditis elegans, daf-42 is essential for development into the dauer stage, as the null mutant of daf-42 exhibited a "no viable dauer" phenotype in which no viable dauers were obtained in any dauer-inducing conditions. Long-term time lapse microscopy of synchronized larvae revealed that daf-42 is involved in developmental changes from the pre-dauer L2d stage to the dauer stage. daf-42 encodes large, disordered proteins of various sizes that are expressed in and secreted from the seam cells within a narrow time window shortly before the molt into dauer stage. Transcriptome analysis showed that the transcription of genes involved in larval physiology and dauer metabolism is highly affected by the daf-42 mutation. Contrary to the notion that essential genes that control the life and death of an organism may be well conserved across diverse species, daf-42 is an evolutionarily young gene conserved only in the Caenorhabditis genus. Our study shows that dauer formation is a vital process that is controlled not only by conserved genes but also by newly emerged genes, providing important insights into evolutionary mechanisms.

C. elegans ageing is accelerated by a self-destructive reproductive programme

In post-reproductive C. elegans, destructive somatic biomass repurposing supports production of yolk which, it was recently shown, is vented and can serve as a foodstuff for larval progeny. This is reminiscent of the suicidal reproductive effort (reproductive death) typical of semelparous organisms such as Pacific salmon. To explore the possibility that C. elegans exhibits reproductive death, we have compared sibling species pairs of the genera Caenorhabditis and Pristionchus with hermaphrodites and females. We report that yolk venting and constitutive, early pathology involving major anatomical changes occur only in hermaphrodites, which are also shorter lived. Moreover, only in hermaphrodites does germline removal suppress senescent pathology and markedly increase lifespan. This is consistent with the hypothesis that C. elegans exhibit reproductive death that is suppressed by germline ablation. If correct, this would imply a major difference in the ageing process between C. elegans and most higher organisms, and potentially explain the exceptional plasticity in C. elegans ageing.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

Repurposing degradation pathways for modular metabolite biosynthesis in nematodes

Recent studies have revealed that Caenorhabditis elegans and other nematodes repurpose products from biochemical degradation pathways for the combinatorial assembly of complex modular structures that serve diverse signaling functions. Building blocks from neurotransmitter, amino acid, nucleoside and fatty acid metabolism are attached to scaffolds based on the dideoxyhexose ascarylose or glucose, resulting in hundreds of modular ascarosides and glucosides. Genome-wide association studies have identified carboxylesterases as the key enzymes mediating modular assembly, enabling rapid compound discovery via untargeted metabolomics and suggesting that modular metabolite biosynthesis originates from the 'hijacking' of conserved detoxification mechanisms. Modular metabolites thus represent a distinct biosynthetic strategy for generating structural and functional diversity in nematodes, complementing the primarily polyketide synthase- and nonribosomal peptide synthetase-derived universe of microbial natural products. Although many aspects of modular metabolite biosynthesis and function remain to be elucidated, their identification demonstrates how phenotype-driven compound discovery, untargeted metabolomics and genomic approaches can synergize to facilitate the annotation of metabolic dark matter.

The CERV protein of Cer1, a C. elegans LTR retrotransposon, is required for nuclear export of viral genomic RNA and can form giant nuclear rods

Retroviruses and closely related LTR retrotransposons export full-length, unspliced genomic RNA (gRNA) for packaging into virions and to serve as the mRNA encoding GAG and POL polyproteins. Because gRNA often includes splice acceptor and donor sequences used to splice viral mRNAs, retroelements must overcome host mechanisms that retain intron-containing RNAs in the nucleus. Here we examine gRNA expression in Cer1, an LTR retrotransposon in C. elegans which somehow avoids silencing and is highly expressed in germ cells. Newly exported Cer1 gRNA associates rapidly with the Cer1 GAG protein, which has structural similarity with retroviral GAG proteins. gRNA export requires CERV (C. elegans regulator of viral expression), a novel protein encoded by a spliced Cer1 mRNA. CERV phosphorylation at S214 is essential for gRNA export, and phosphorylated CERV colocalizes with nuclear gRNA at presumptive sites of transcription. By electron microscopy, tagged CERV proteins surround clusters of distinct, linear fibrils that likely represent gRNA molecules. Single fibrils, or groups of aligned fibrils, also localize near nuclear pores. During the C. elegans self-fertile period, when hermaphrodites fertilize oocytes with their own sperm, CERV concentrates in two nuclear foci that are coincident with gRNA. However, as hermaphrodites cease self-fertilization, and can only produce cross-progeny, CERV undergoes a remarkable transition to form giant nuclear rods or cylinders that can be up to 5 microns in length. We propose a novel mechanism of rod formation, in which stage-specific changes in the nucleolus induce CERV to localize to the nucleolar periphery in flattened streaks of protein and gRNA; these streaks then roll up into cylinders. The rods are a widespread feature of Cer1 in wild strains of C. elegans, but their function is not known and might be limited to cross-progeny. We speculate that the adaptive strategy Cer1 uses for the identical self-progeny of a host hermaphrodite might differ for heterozygous cross-progeny sired by males. For example, mating introduces male chromosomes which can have different, or no, Cer1 elements.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

Background

The nematode Caenorhabditis briggsae has been used as a model for genomics studies compared to Caenorhabditis elegans because of its striking morphological and behavioral similarities. These studies yielded numerous findings that have expanded our understanding of nematode development and evolution. However, the potential of C. briggsae to study nematode biology is limited by the quality of its genome resources. The reference genome and gene models 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 protein-coding gene predictions generated from short- and long-read transcriptomic data. Because of the limitations of gene prediction software, the existing gene models for QX1410 contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 software-derived gene models and underlying transcriptomic data to improve the protein-coding gene models of the C. briggsae QX1410 genome.

Results

We designed a detailed workflow to train a team of nine students to manually curate genes using RNA read alignments and predicted gene models. We manually inspected the gene models using the genome annotation editor, Apollo, and proposed corrections to the coding sequences of over 8,000 genes. Additionally, we 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 before and after curation. Manual curation led to a substantial improvement in the protein sequence length accuracy of QX1410 genes. We also compared the curated QX1410 gene models against the existing AF16 gene models. The manual curation efforts yielded QX1410 gene models that are similar in quality to the extensively curated AF16 gene models in terms of protein-length accuracy and biological completeness scores. 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. Comparative genomic analysis using a related species with high-quality reference genome(s) and gene models can be used to quantify improvements in gene model quality in a newly sequenced genome. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. The chromosome-level reference genome for the C. briggsae strain QX1410 far surpasses the quality of the genome of the laboratory strain AF16, and our manual curation efforts have brought the QX1410 gene models to a comparable level of quality to the previous reference, AF16. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

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.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

Ionizing radiation affects the demography and the evolution of Caenorhabditis elegans populations

Ionizing radiation can reduce survival, reproduction and affect development, and lead to the extinction of populations if their evolutionary response is insufficient. However, demographic and evolutionary studies on the effects of ionizing radiation are still scarce. Using an experimental evolution approach, we analyzed population growth rate and associated change in life history traits across generations in Caenorhabditis elegans populations exposed to 0, 1.4, and 50.0 mGy.h^-1 of ionizing radiation (gamma external irradiation). We found a higher population growth rate in the 1.4 mGy.h^-1 treatment and a lower in the 50.0 mGy.h^-1 treatment compared to the control. Realized fecundity was lower in both 1.4 and 50.0 mGy.h^-1 than control treatment. High irradiation levels decreased brood size from self-fertilized hermaphrodites, specifically early brood size. Finally, high irradiation levels decreased hatching success compared to the control condition. In reciprocal-transplant experiments, we found that life in low irradiation conditions led to the evolution of higher hatching success and late brood size. These changes could provide better tolerance against ionizing radiation, investing more in self-maintenance than in reproduction. These evolutionary changes were with some costs of adaptation. This study shows that ionizing radiation has both demographic and evolutionary consequences on populations.

Holotype sequencing of Silvataresholzenthali Rázuri-Gonzales, Ngera & Pauls, 2022 (Trichoptera, Pisuliidae)

While DNA barcodes are increasingly provided in descriptions of new species, the whole mitochondrial and nuclear genomes are still rarely included. This is unfortunate because whole genome sequencing of holotypes allows perpetual genetic characterization of the most representative specimen for a given species. Thus, de novo genomes are invaluable additional diagnostic characters in species descriptions, provided the structural integrity of the holotype specimens remains intact. Here, we used a minimally invasive method to extract DNA of the type specimen of the recently described caddisfly species Silvataresholzenthali Rázuri-Gonzales, Ngera & Pauls, 2022 (Trichoptera: Pisuliidae) from the Democratic Republic of the Congo. A low-cost next generation sequencing strategy was used to generate the complete mitochondrial and draft nuclear genome of the holotype. The data in its current form is an important extension to the morphological species description and valuable for phylogenomic studies.

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.

Making "Sense" of Ecology from a Genetic Perspective: Caenorhabditis elegans, Microbes and Behavior

Our knowledge of animal and behavior in the natural ecology is based on over a century's worth of valuable field studies. In this post-genome era, however, we recognize that genes are the underpinning of ecological interactions between two organisms. Understanding how genes contribute to animal ecology, which is essentially the intersection of two genomes, is a tremendous challenge. The bacterivorous nematode Caenorhabditis elegans, one of the most well-known genetic animal model experimental systems, experiences a complex microbial world in its natural habitat, providing us with a window into the interplay of genes and molecules that result in an animal-microbial ecology. In this review, we will discuss C. elegans natural ecology, how the worm uses its sensory system to detect the microbes and metabolites that it encounters, and then discuss some of the fascinating ecological dances, including behaviors, that have evolved between the nematode and the microbes in its environment.

Opposing directions of stage-specific body shape change in a close relative of C. elegans

Background

Body size is a fundamental organismal trait. However, as body size and ecological contexts change across developmental time, evolutionary divergence may cause unexpected patterns of body size diversity among developmental stages. This may be particularly evident in polyphenic developmental stages specialized for dispersal. The dauer larva is such a stage in nematodes, and Caenorhabditis species disperse by traveling on invertebrate carriers. Here, we describe the morphology of a stress-resistant, dauer-like larval stage of the nematode Caenorhabditis inopinata, whose adults can grow to be nearly twice as long as its close relative, the model organism C. elegans.

Results

We find that a dauer-like, stress-resistant larval stage in two isolates of C. inopinata is on average 13% shorter and 30% wider than the dauer larvae of C. elegans, despite its much longer adult stage. Additionally, many C. inopinata dauer-like larvae were ensheathed, a possible novelty in this lineage reminiscent of the infective juveniles of parasitic nematodes. Variation in dauer-like larva formation frequency among twenty-four wild isolates of C. inopinata was also observed, although frequencies were low across all isolates (< 2%), with many isolates unable to produce dauer-like larvae under conventional laboratory conditions.

Conclusion

Most Caenorhabditis species thrive on rotting plants and disperse on snails, slugs, or isopods (among others) whereas C. inopinata is ecologically divergent and thrives in fresh Ficus septica figs and disperses on their pollinating wasps. While there is some unknown factor of the fig environment that promotes elongated body size in C. inopinata adults, the small size or unique life history of its fig wasp carrier may be driving the divergent morphology of its stress-resistant larval stages. Further characterization of the behavior, development, and morphology of this stage will refine connections to homologous developmental stages in other species and determine whether ecological divergence across multiple developmental stages can promote unexpected and opposing changes in body size dimensions within a single species.

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.

Reference-Guided De Novo Genome Assembly of the Flour Beetle Tribolium freemani

The flour beetle Tribolium freemani is a sibling species of the model organism and important pest Tribolium castaneum. The two species are so closely related that they can produce hybrid progeny, but the genetic basis of their differences has not been revealed. In this work, we sequenced the T. freemani genome by applying PacBio HiFi technology. Using the well-assembled T. castaneum genome as a reference, we assembled 262 Mb of the T. freemani genomic sequence and anchored it in 10 linkage groups corresponding to nine autosomes and sex chromosome X. The assembly showed 99.8% completeness of conserved insect genes, indicating a high-quality reference genome. Comparison with the T. castaneum assembly revealed that the main differences in genomic sequence between the two sibling species come from repetitive DNA, including interspersed and tandem repeats. In this work, we also provided the complete assembled mitochondrial genome of T. freemani. Although the genome assembly needs to be ameliorated in tandemly repeated regions, the first version of the T. freemani reference genome and the complete mitogenome presented here represent useful resources for comparative evolutionary studies of related species and for further basic and applied research on different biological aspects of economically important pests.

Transgenesis of the gonochoristic nematode Caenorhabditis inopinata by microparticle bombardment with hygromycin B selection

The gonochoristic nematode Caenorhabditis inopinata is the phylogenetically closest species to the well-studied nematode Caenorhabditis elegans (Kanzaki et al. , 2018). While C. inopinata has been expected to be a useful comparative model for C. elegans , efficient transgenesis methods have not been available. Here, we established a method to integrate transgenes into the C. inopinata genome by microparticle bombardment with hygromycin B selection. C. elegans- derived genetic elements tested in this study, including universal and germline-specific promoters, ORFs, and 3’UTRs, were all functional in C. inopinata. Using this method, transgenic C. inopinata strains that express fluorescent subcellular markers were established.

Bursaphelenchus suri n. sp.: A second Bursaphelenchus syconial parasite of figs supports adaptive radiation among section Sycomorus figs

The nematode genus Bursaphelenchus is a highly divergent group. This genus mainly consists of mycophagous entomophilic species, but some species have specialized as obligate or facultative plant parasites, facultative insect parasites, or exhibit feeding dimorphism (phenotypic plasticity) leading to mycophagous and predatory forms. In the present study, a new Bursaphelenchus species, B. suri n. sp. was isolated from fresh syconia (figs) of Ficus sur and is described and illustrated based on its typological characters and molecular phylogenetic status. The new species is characterized by its highly derived feeding structures found in obligate plant parasites, lip possessing a labial disc and a long and thick stylet with a long conus and extremely well-developed basal swellings. In addition, slender body of both sexes is characteristic of the species. The new species is phylogenetically and typologically closely related to B. sycophilus, i.e., these two species share the characteristic feeding structures and form a well-supported clade within the B. fungivorus group in the genus. Biologically, these two species are both isolated from fresh figs of the section Sycomorus. However, the new species differs from B. sycophilus by the length of the female post-uterine sac and the shape of the male spicule, i.e., the new species has a long post-uterine sac and spicule condylus without dorsal recurvature. Thus, the new species is the second obligate fig parasite of the genus, and the evolutionary relationship between the B. suri n. sp. and B. sycophilus clade and section Sycomorus figs is hypothesized as an example of adaptive radiation with more species to be discovered.

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.

Chromosome-level reference genomes for two strains of Caenorhabditis briggsae: an improved platform for comparative genomics

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines (RILs) generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing (RNA-seq) data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Association of Fig Pollinating Wasps and Fig Nematodes inside Male and Female Figs of a Dioecious Fig Tree in Sumatra, Indonesia

Simple Summary

It has been known that fig pollinating wasps can transfer nematodes when they enter receptive figs to deposit their eggs (male figs in male trees) and pollinate the flowers (female figs in female trees) of a dioecious fig tree. However, the fate of nematodes transferred into female figs was unknown, since the pollinating wasps that enter female figs cannot reproduce. This study examined whether the nematodes transferred by pollinating wasps into female figs of Ficus hispida could develop and reproduce in the same way as in male figs. Three species of nematodes were found to develop within the male and female figs of F. hispida, with populations peaking at about the same time for both sexes of figs. Our findings showed that the female figs could support the growth and reproduction of the three nematodes; however, there was no pollinator offspring to transfer the nematodes out of the figs. The results provide a better understanding of the development of nematodes in male and female figs of a dioecious fig tree, which may also influence the biodiversity of the tropical ecosystem.

Abstract

Nematodes can grow within the inflorescences of many fig trees (Ficus spp., Moraceae); however, the feeding behaviour of most nematodes is not known. Fig pollinating wasps (Hymenoptera: Agaonidae) transfer nematodes into young figs upon the wasps’ entry into the figs to deposit their eggs. Most Asian fig trees, however, are functionally dioecious, and the pollinating wasps that enter female figs are unable to reproduce. They fail to produce the offspring required to carry the new generations of nematodes. We examined whether female figs of F. hispida can nonetheless support the development of phoretic nematode populations. Nematodes were extracted from male and female figs sampled in Sumatra, Indonesia, to compare the growth of their populations within the figs. We found three species of nematodes that grew within figs of male and female trees of F. hispida: Ficophagus cf. centerae (Aphelenchoididae), Martininema baculum (Aphelenchoididae) and Caenorhabditis sp (Rhabditidae). The latter species (Caenorhabditis sp.) has never been reported to be associated with F. hispida before. Nematode populations peaked at around 120–140 individuals in both sexes of figs, at the time when a succeeding generation of adult fig wasps appeared within male figs. The female figs could support the growth and reproduction of the three nematodes species; however, the absence of vectors meant that female figs remained as traps from which there could be no escape.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

The fifth Japanese meeting on biological function and evolution through interactions between hosts and transposable elements

The fifth Japanese meeting on host–transposon interactions, titled “Biological Function and Evolution through Interactions between Hosts and Transposable Elements (TEs),” was held online on August 26–27, 2021. The meeting was supported by National Institute of Genetics and aimed to bring together researchers studying the diverse roles of TEs in genome function and evolution, as well as host defense systems against TE mobility by chromatin and RNA modifications and protein-protein interactions. Here, we present the highlights of the talks.

The distribution of mutational effects on fitness in Caenorhabditis elegans inferred from standing genetic variation

The distribution of fitness effects (DFE) for new mutations is one of the most theoretically important but difficult to estimate properties in population genetics. A crucial challenge to inferring the DFE from natural genetic variation is the sensitivity of the site frequency spectrum to factors like population size change, population substructure, genome structure, and nonrandom mating. Although inference methods aim to control for population size changes, the influence of nonrandom mating remains incompletely understood, despite being a common feature of many species. We report the DFE estimated from 326 genomes of Caenorhabditis elegans, a nematode roundworm with a high rate of self-fertilization. We evaluate the robustness of DFE inferences using simulated data that mimics the genomic structure and reproductive life history of C. elegans. Our observations demonstrate how the combined influence of self-fertilization, genome structure, and natural selection on linked sites can conspire to compromise estimates of the DFE from extant polymorphisms with existing methods. These factors together tend to bias inferences toward weakly deleterious mutations, making it challenging to have full confidence in the inferred DFE of new mutations as deduced from standing genetic variation in species like C. elegans. Improved methods for inferring the DFE are needed to appropriately handle strong linked selection and selfing. These results highlight the importance of understanding the combined effects of processes that can bias our interpretations of evolution in natural populations.

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.

Combinatorial Assembly of Modular Glucosides via Carboxylesterases Regulates C. elegans Starvation Survival

The recently discovered modular glucosides (MOGLs) form a large metabolite library derived from combinatorial assembly of moieties from amino acid, neurotransmitter, and lipid metabolism in the model organism C. elegans. Combining CRISPR-Cas9 genome editing, comparative metabolomics, and synthesis, we show that the carboxylesterase homologue Cel-CEST-1.2 is responsible for specific 2-O-acylation of diverse glucose scaffolds with a wide variety of building blocks, resulting in more than 150 different MOGLs. We further show that this biosynthetic role is conserved for the closest homologue of Cel-CEST-1.2 in the related nematode species C. briggsae, Cbr-CEST-2. Expression of Cel-cest-1.2 and MOGL biosynthesis are strongly induced by starvation conditions in C. elegans, one of the premier model systems for mechanisms connecting nutrition and physiology. Cel-cest-1.2-deletion results in early death of adult animals under starvation conditions, providing first insights into the biological functions of MOGLs.

Balancing selection maintains hyper-divergent haplotypes in Caenorhabditis elegans

Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing decreases the effective population size, effective recombination rate and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis, selfing has evolved at least three times, and all three species, including the model organism Caenorhabditis elegans, show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long-read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example of how species can avoid the evolutionary dead end associated with selfing.

The double-stranded DNA-binding proteins TEBP-1 and TEBP-2 form a telomeric complex with POT-1

Telomeres are bound by dedicated proteins, which protect them from DNA damage and regulate telomere length homeostasis. In the nematode Caenorhabditis elegans, a comprehensive understanding of the proteins interacting with the telomere sequence is lacking. Here, we harnessed a quantitative proteomics approach to identify TEBP-1 and TEBP-2, two paralogs expressed in the germline and embryogenesis that associate to telomeres in vitro and in vivo. tebp-1 and tebp-2 mutants display strikingly distinct phenotypes: tebp-1 mutants have longer telomeres than wild-type animals, while tebp-2 mutants display shorter telomeres and a Mortal Germline. Notably, tebp-1;tebp-2 double mutant animals have synthetic sterility, with germlines showing signs of severe mitotic and meiotic arrest. Furthermore, we show that POT-1 forms a telomeric complex with TEBP-1 and TEBP-2, which bridges TEBP-1/-2 with POT-2/MRT-1. These results provide insights into the composition and organization of a telomeric protein complex in C. elegans.