Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome

Molecular identification of a peroxidase gene controlling body size in the entomopathogenic nematode Steinernema hermaphroditum

The entomopathogenic nematode Steinernema hermaphroditum was recently rediscovered and is being developed as a genetically tractable experimental system for the study of previously unexplored biology, including parasitism of its insect hosts and mutualism with its bacterial endosymbiont Xenorhabdus griffiniae. Through whole-genome re-sequencing and genetic mapping we have for the first time molecularly identified the gene responsible for a mutationally defined phenotypic locus in an entomopathogenic nematode. In the process we observed an unexpected mutational spectrum following ethyl methansulfonate mutagenesis in this species. We find that the ortholog of the essential Caenorhabditis elegans peroxidase gene skpo-2 controls body size and shape in S. hermaphroditum. We confirmed this identification by generating additional loss-of-function mutations in the gene using CRISPR-Cas9. We propose that the identification of skpo-2 will accelerate gene targeting in other Steinernema entomopathogenic nematodes used commercially in pest control, as skpo-2 is X-linked and males hemizygous for loss of its function can mate, making skpo-2 an easily recognized and maintained marker for use in co-CRISPR.

CRISPR-Cas9 genome editing in Steinernema entomopathogenic nematodes

Molecular tool development in traditionally non-tractable animals opens new avenues to study gene functions in the relevant ecological context. Entomopathogenic nematodes (EPN) Steinernema and their symbiotic bacteria of Xenorhabdus spp are a valuable experimental system in the laboratory and are applicable in the field to promote agricultural productivity. The infective juvenile (IJ) stage of the nematode packages mutualistic symbiotic bacteria in the intestinal pocket and invades insects that are agricultural pests. The lack of consistent and heritable genetics tools in EPN targeted mutagenesis severely restricted the study of molecular mechanisms underlying both parasitic and mutualistic interactions. Here, I report a protocol for CRISPR-Cas9 based genome-editing that is successful in two EPN species, S. carpocapsae and S. hermaphroditum . I adapted a gonadal microinjection technique in S. carpocapsae , which created on-target modifications of a homologue Sc-dpy-10 (cuticular collagen) by homology-directed repair. A similar delivery approach was used to introduce various alleles in S. hermaphroditum including Sh-dpy-10 and Sh-unc-22 (a muscle gene), resulting in visible and heritable phenotypes of dumpy and twitching, respectively. Using conditionally dominant alleles of Sh-unc-22 as a co-CRISPR marker, I successfully modified a second locus encoding Sh-Daf-22 (a homologue of human sterol carrier protein SCPx), predicted to function as a core enzyme in the biosynthesis of nematode pheromone that is required for IJ development. As a proof of concept, Sh-daf-22 null mutant showed IJ developmental defects in vivo ( in insecta) . This research demonstrates that Steinernema spp are highly tractable for targeted mutagenesis and has great potential in the study of gene functions under controlled laboratory conditions within the relevant context of its ecological niche.

Low-temperature exposure has immediate and lasting effects on the stress tolerance, chemotaxis and proteome of entomopathogenic nematodes

Temperature is one of the most important factors affecting soil organisms, including the infective stages of parasites and entomopathogenic nematodes, which are important biological control agents. We investigated the response of 2 species of entomopathogenic nematodes to different storage regimes: cold (9°C), culture temperature (20°C) and temperature swapped from 9 to 20°C. For Steinernema carpocapsae, cold storage had profound effects on chemotaxis, stress tolerance and protein expression that were retained in temperature-swapped individuals. These effects included reversal of chemotactic response for 3 (prenol, methyl salicylate and hexanol) of the 4 chemicals tested, and enhanced tolerance to freezing (−10°C) and desiccation (75% RH). Label-free quantitative proteomics showed that cold storage induced widespread changes in S. carpocapsae, including an increase in heat-shock proteins and late embryogenesis abundant proteins. For Heterorhabditis megidis, cold storage had a less dramatic effect on chemotaxis (as previously shown for proteomic expression) and changes were not maintained on return to 20°C. Thus, cold temperature exposure has significant effects on entomopathogenic nematodes, but the nature of the change depends on the species. Steinernema carpocapsae, in particular, displays significant plasticity, and its behaviour and stress tolerance may be manipulated by brief exposure to low temperatures, with implications for its use as a biological control agent.

The effect of temperature conditioning (9°C and 20°C) on the proteome of entomopathogenic nematode infective juveniles

Entomopathogenic nematodes (EPN) of the genera Steinernema and Heterorhabditis are parasites which kill and reproduce within insects. While both have life cycles centred around their developmentally arrested, nonfeeding and stress tolerant infective juvenile (IJ) stage, they are relatively distantly related. These IJs are promising biocontrol agents, and their shelf life and stress tolerance may be enhanced by storage at low temperatures. The purpose of this study was to investigate how the proteome of the IJs of two distantly related EPN species is affected by storage at 9°C (for up to 9 weeks) and 20°C (for up to 6 weeks), using label-free quantitative proteomics. Overall, more proteins were detected in S. carpocapsae (2422) than in H. megidis (1582). The S. carpocapsae proteome was strongly affected by temperature, while the H. megidis proteome was affected by both time and temperature. The proteins which increased in abundance to the greatest extent in S. carpocapsae IJs after conditioning at 9°C were chaperone proteins, and proteins related to stress. The proteins which increased in abundance the most after storage at 20°C were proteins related to the cytoskeleton, cell signalling, proteases and their inhibitors, which may have roles in infection. The proteins which decreased in abundance to the greatest extent in S. carpocapsae after both 9°C and 20°C storage were those associated with metabolism, stress and the cytoskeleton. After storage at both temperatures, the proteins increased to the greatest extent in H. megidis IJs were those associated with the cytoskeleton, cell signalling and carbon metabolism, and the proteins decreased in abundance to the greatest extent were heat shock and ribosomal proteins, and those associated with metabolism. As the longest-lived stage of the EPN life cycle, IJs may be affected by proteostatic stress, caused by the accumulation of misfolded proteins and toxic aggregates. The substantial increase of chaperone proteins in S. carpocapsae, and to a greater extent at 9°C, and the general decrease in ribosomal and chaperone proteins in H. megidis may represent species-specific proteostasis mechanisms. Similarly, organisms accumulate reactive oxygen species (ROS) over time and both species exhibited a gradual increase in proteins which enhance ROS tolerance, such as catalase. The species-specific responses of the proteome in response to storage temperature, and over time, may reflect the phylogenetic distance and/or different ecological strategies.

The entomopathogenic nematode Steinernema hermaphroditum is a self-fertilizing hermaphrodite and a genetically tractable system for the study of parasitic and mutualistic symbiosis

Entomopathogenic nematodes (EPNs), including Heterorhabditis and Steinernema, are parasitic to insects and contain mutualistically symbiotic bacteria in their intestines (Photorhabdus and Xenorhabdus, respectively) and therefore offer opportunities to study both mutualistic and parasitic symbiosis. The establishment of genetic tools in EPNs has been impeded by limited genetic tractability, inconsistent growth in vitro, variable cryopreservation, and low mating efficiency. We obtained the recently described Steinernema hermaphroditum strain CS34 and optimized its in vitro growth, with a rapid generation time on a lawn of its native symbiotic bacteria Xenorhabdus griffiniae. We developed a simple and efficient cryopreservation method. Previously, S. hermaphroditum isolated from insect hosts was described as producing hermaphrodites in the first generation. We discovered that CS34, when grown in vitro, produced consecutive generations of autonomously reproducing hermaphrodites accompanied by rare males. We performed mutagenesis screens in S. hermaphroditum that produced mutant lines with visible and heritable phenotypes. Genetic analysis of the mutants demonstrated that this species reproduces by self-fertilization rather than parthenogenesis and that its sex is determined chromosomally. Genetic mapping has thus far identified markers on the X chromosome and three of four autosomes. We report that S. hermaphroditum CS34 is the first consistently hermaphroditic EPN and is suitable for genetic model development to study naturally occurring mutualistic symbiosis and insect parasitism.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Genome announcement of Steinernema khuongi and its associated symbiont from Florida

Citrus root weevil (Diaprepes abbreviates) causes significant yield loss in citrus, especially in Florida. A promising source of control for this pest is biological control agents, namely, native entomopathogenic nematodes (EPNs) within the genus Steinernema. Two species of endemic EPN in Florida are S. diaparepesi, abundant within the central ridge, and S. khuongi, dominating the flatwood regions of the state. These citrus-growing regions differ significantly in their soil habitats, which impacts the potential success of biological control measures. Although the genome sequence of S. diaprepesi is currently available, the genome sequence of S. khuongi and identity of the symbiotic bacteria is still unknown. Understanding the genomic differences between these two nematodes and their favored habitats can inform successful biological control practices. Here, MiSeq libraries were used to simultaneously sequence and assemble the draft genome of S. khuongi and its associated symbionts. The final draft genome for S. khuongi has 8,794 contigs with a total length of ∼82 Mb, a largest contig of 428,226 bp, and N50 of 46 kb; its BUSCO scores indicate that it is > 86% complete. An associated bacterial genome was assembled with a total length of ∼3.5 Mb, a largest contig at 116,532 bp, and N50 of 17,487 bp. The bacterial genome encoded 3,721 genes, similar to other Xenorhabdus genomes. Comparative genomics identified the symbiotic bacteria of S. khuongi as Xenorhabdus poinarii. These new draft genomes of a host and symbiont can be used as a valuable tool for comparative genomics with other EPNs and its symbionts to understand host range and habitat suitability.

Transcriptome Profile Changes Associated With Heat Shock Reaction in the Entomopathogenic Nematode, Steinernema carpocapsae

The entomopathogenic nematode Steinernema carpocapsae has been used for control of soil insects. However, S. carpocapse is sensitive to environmental factors, particularly temperature. We studied an S. carpocapse group that was shocked with high temperature. We also studied the transcriptome-level responses associated with temperature stress using a BGIseq sequencing platform. We de novo assembled the reads from the treatment and control groups into one transcriptome consisting of 43.9 and 42.9 million clean reads, respectively. Based on the genome database, we aligned the clean reads to the Nr, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases and analyzed the differentially expressed genes (DEGs). Compared with the control, the heat-shocked group had significant differential expression of the heat shock protein (HSP) family, antioxidase [glutathione S-transferases (GSTs) and superoxide dismutase (SOD)], monooxygenase (P450), and transcription factor genes (DAF-16 and DAF-2). These DEGs were demonstrated to be part of the Longevity pathway and insulin/insulin-like signaling pathway. The results revealed the potential mechanisms, at the transcriptional level, of S. carpocapsae under thermal stress.