The biocontrol nematode Phasmarhabditis hermaphrodita infects and increases mortality of Monadenia fidelis, a non-target terrestrial gastropod species endemic to the Pacific Northwest of North America, in laboratory conditions.

Inundative biological control (biocontrol) efforts in pest management lead to the mass distribution of commercialized biocontrol agents. Many 'biocontrol gone awry' incidents have resulted in disastrous biodiversity impacts, leading to increased scrutiny of biocontrol efforts. The nematode Phasmarhabditis hermaphrodita is sold as a biocontrol agent on three continents and targets pest gastropods such as Deroceras reticulatum, the Grey Field Slug; P. hermaphrodita is not presently approved for use in the United States. Investigations into the potential for P. hermaphrodita to infect non-target gastropod species of conservation relevance, however, are limited. We examined the effects of three strains of P. hermaphrodita on mortality in Monadenia fidelis, the Pacific Sideband, a snail species endemic to the Pacific Northwest of North America, in laboratory conditions. Across a 71-day laboratory infectivity assay, snails exposed to each of the three nematode strains, each analyzed at two doses, experienced a mean 50% mortality by days 20-42. All nematode-treated snails were dead by the end of the study. By contrast, 30/30 water-control snails experienced no mortality. Nematodes killed smaller, juvenile-stage snails significantly faster than those in larger and more developmentally advanced hosts. Our results provide direct evidence that the biocontrol nematode P. hermaphrodita infects and kills M. fidelis, a non-target gastropod species endemic to the Pacific Northwest, in laboratory conditions. This study suggests that introduction of P. hermaphrodita to new ecosystems might negatively impact endemic gastropod biodiversity and advocates for further investigation of non-target effects, including in conditions closer to the natural environments of non-target species.

Genomic Analyses of Globodera pallida, A Quarantine Agricultural Pathogen in Idaho.

Globodera pallida is among the most significant plant-parasitic nematodes worldwide, causing major damage to potato production. Since it was discovered in Idaho in 2006, eradication efforts have aimed to contain and eradicate G. pallida through phytosanitary action and soil fumigation. In this study, we investigated genome-wide patterns of G. pallida genetic variation across Idaho fields to evaluate whether the infestation resulted from a single or multiple introduction(s) and to investigate potential evolutionary responses since the time of infestation. A total of 53 G. pallida samples (~1,042,000 individuals) were collected and analyzed, representing five different fields in Idaho, a greenhouse population, and a field in Scotland that was used for external comparison. According to genome-wide allele frequency and fixation index (Fst) analyses, most of the genetic variation was shared among the G. pallida populations in Idaho fields pre-fumigation, indicating that the infestation likely resulted from a single introduction. Temporal patterns of genome-wide polymorphisms involving (1) pre-fumigation field samples collected in 2007 and 2014 and (2) pre- and post-fumigation samples revealed nucleotide variants (SNPs, single-nucleotide polymorphisms) with significantly differentiated allele frequencies indicating genetic differentiation. This study provides insights into the genetic origins and adaptive potential of G. pallida invading new environments.

Phylogenetic evidence for the invasion of a commercialized European Phasmarhabditis hermaphrodita lineage into North America and New Zealand.

Biological control (biocontrol) as a component of pest management strategies reduces reliance on synthetic chemicals, and seemingly offers a natural approach that minimizes environmental impact. However, introducing a new organism to new environments as a classical biocontrol agent can have broad and unanticipated biodiversity effects and conservation consequences. Nematodes are currently used in a variety of commercial biocontrol applications, including the use of Phasmarhabditis hermaphrodita as an agent targeting pest slug and snail species. This species was originally discovered in Germany, and is generally thought to have European origins. P. hermaphrodita is sold under the trade name Nemaslug®, and is available only in European markets. However, this nematode species was discovered in New Zealand and the western United States, though its specific origins remained unclear. In this study, we analyzed 45 nematode strains representing eight different Phasmarhabditis species, collected from nine countries around the world. A segment of nematode mitochondrial DNA (mtDNA) was sequenced and subjected to phylogenetic analyses. Our mtDNA phylogenies were overall consistent with previous analyses based on nuclear ribosomal RNA (rRNA) loci. The recently discovered P. hermaphrodita strains in New Zealand and the United States had mtDNA haplotypes nearly identical to that of Nemaslug®, and these were placed together in an intraspecific monophyletic clade with high support in maximum likelihood and Bayesian analyses. We also examined bacteria that co-cultured with the nematode strains isolated in Oregon, USA, by analyzing 16S rRNA sequences. Eight different bacterial genera were found to associate with these nematodes, though Moraxella osloensis, the bacteria species used in the Nemaslug® formulation, was not detected. This study provided evidence that nematodes deriving from the Nemaslug® biocontrol product have invaded countries where its use is prohibited by regulatory agencies and not commercially available.

First report of molecular characterization and phylogeny of Trichuris fossor Hall, 1916 (Nematoda: Trichuridae).

Because species of Trichuris are morphologically similar and ranges of host preference are variable, using molecular data to evaluate species delineations is essential for properly quantifying biodiversity of and relationships within Trichuridae. Trichuris fossor has been reported from Thomomys spp. (Rodentia: Geomyidae, 'pocket gophers') hosts based on morphological features alone. Partial 18S rRNA sequences for specimens identified as T. fossor based on morphology, along with sequences from 26 additional taxa, were used for a phylogenetic analysis. Evolutionary histories were constructed using maximum likelihood and Bayesian inference. In both analyses, the specimens fell within the Trichuris clade with 100% support and formed a distinct subclade with 100% support. These results confirm that T. fossor is a distinct species and represent the first molecular report for it. Relatedness among species within the family were well resolved in the BI tree. This study represents an initial effort to obtain a more comprehensive view of Trichuridae by including a new clade member, T. fossor. A better understanding of Trichuridae phylogeny could contribute to further characterization of host-associations, including species that infect livestock and humans.Because species of Trichuris are morphologically similar and ranges of host preference are variable, using molecular data to evaluate species delineations is essential for properly quantifying biodiversity of and relationships within Trichuridae. Trichuris fossor has been reported from Thomomys spp. (Rodentia: Geomyidae, 'pocket gophers') hosts based on morphological features alone. Partial 18S rRNA sequences for specimens identified as T. fossor based on morphology, along with sequences from 26 additional taxa, were used for a phylogenetic analysis. Evolutionary histories were constructed using maximum likelihood and Bayesian inference. In both analyses, the specimens fell within the Trichuris clade with 100% support and formed a distinct subclade with 100% support. These results confirm that T. fossor is a distinct species and represent the first molecular report for it. Relatedness among species within the family were well resolved in the BI tree. This study represents an initial effort to obtain a more comprehensive view of Trichuridae by including a new clade member, T. fossor. A better understanding of Trichuridae phylogeny could contribute to further characterization of host-associations, including species that infect livestock and humans.

Mitochondrial DNA Variation and Selfish Propagation Following Experimental Bottlenecking in Two Distantly Related Caenorhabditis briggsae Isolates.

Understanding mitochondrial DNA (mtDNA) evolution and inheritance has broad implications for animal speciation and human disease models. However, few natural models exist that can simultaneously represent mtDNA transmission bias, mutation, and copy number variation. Certain isolates of the nematode Caenorhabditis briggsae harbor large, naturally-occurring mtDNA deletions of several hundred basepairs affecting the NADH dehydrogenase subunit 5 (nduo-5) gene that can be functionally detrimental. These deletion variants can behave as selfish DNA elements under genetic drift conditions, but whether all of these large deletion variants are transmitted in the same preferential manner remains unclear. In addition, the degree to which transgenerational mtDNA evolution profiles are shared between isolates that differ in their propensity to accumulate the nduo-5 deletion is also unclear. We address these knowledge gaps by experimentally bottlenecking two isolates of C. briggsae with different nduo-5 deletion frequencies for up to 50 generations and performing total DNA sequencing to identify mtDNA variation. We observed multiple mutation profile differences and similarities between C. briggsae isolates, a potentially species-specific pattern of copy number dysregulation, and some evidence for genetic hitchhiking in the deletion-bearing isolate. Our results further support C. briggsae as a practical model for characterizing naturally-occurring mtgenome variation and contribute to the understanding of how mtgenome variation persists in animal populations and how it presents in mitochondrial disease states.

Complex Transmission Patterns and Age-Related Dynamics of a Selfish mtDNA Deletion.

Despite wide-ranging implications of selfish mitochondrial DNA (mtDNA) elements for human disease and topics in evolutionary biology (e.g., speciation), the forces controlling their formation, age-related accumulation, and offspring transmission remain largely unknown. Selfish mtDNA poses a significant challenge to genome integrity, mitochondrial function, and organismal fitness. For instance, numerous human diseases are associated with mtDNA mutations; however, few genetic systems can simultaneously represent pathogenic mitochondrial genome evolution and inheritance. The nematode Caenorhabditis briggsae is one such system. Natural C. briggsae isolates harbor varying levels of a large-scale deletion affecting the mitochondrial nduo-5 gene, termed nad5Δ. A subset of these isolates contains putative compensatory mutations that may reduce the risk of deletion formation. We studied the dynamics of nad5Δ heteroplasmy levels during animal development and transmission from mothers to offspring in genetically diverse C. briggsae natural isolates. Results support previous work demonstrating that nad5Δ is a selfish element and that heteroplasmy levels of this deletion can be quite plastic, exhibiting high degrees of inter-family variability and divergence between generations. The latter is consistent with a mitochondrial bottleneck effect, and contrasts with previous findings from a laboratory-derived model uaDf5 mtDNA deletion in C. elegans. However, we also found evidence for among-isolate differences in the ability to limit nad5Δ accumulation, the pattern of which suggested that forces other than the compensatory mutations are important in protecting individuals and populations from rampant mtDNA deletion expansion over short time scales.

Variable Abundance and Distribution of Wolbachia and Cardinium Endosymbionts in Plant-Parasitic Nematode Field Populations.

The bacterial endosymbiont Wolbachia interacts with different invertebrate hosts, engaging in diverse symbiotic relationships. Wolbachia is often a reproductive parasite in arthropods, but an obligate mutualist in filarial nematodes. Wolbachia was recently discovered in plant-parasitic nematodes, and, is thus far known in just two genera Pratylenchus and Radopholus, yet the symbiont's function remains unknown. The occurrence of Wolbachia in these economically important plant pests offers an unexplored biocontrol strategy. However, development of Wolbachia-based biocontrol requires an improved understanding of symbiont-host functional interactions and the symbiont's prevalence among nematode field populations. This study used a molecular-genetic approach to assess the prevalence of a Wolbachia lineage (wPpe) in 32 field populations of Pratylenchus penetrans. Populations were examined from eight different plant species in Washington, Oregon, and California. Nematodes were also screened for the endosymbiotic bacterium Cardinium (cPpe) that was recently shown to co-infect P. penetrans. Results identified wPpe in 9/32 and cPpe in 1/32 of P. penetrans field populations analyzed. No co-infection was observed in field populations. Wolbachia was detected in nematodes from 4/8 plant-hosts examined (raspberry, strawberry, clover, and lily), and in all three states surveyed. Cardinium was detected in nematodes from mint in Washington. In the wPpe-infected P. penetrans populations collected from raspberry, the prevalence of wPpe infection ranged from 11 to 58%. This pattern is unlike that in filarial nematodes where Wolbachia is an obligate mutualist and occurs in 100% of the host. Further analysis of wPpe-infected populations revealed female-skewed sex ratios (up to 96%), with the degree of skew positively correlating with wPpe prevalence. Uninfected nematode populations had approximately equal numbers of males and females. Comparisons of 54 wPpe 16S ribosomal RNA sequences revealed high similarity across the geographic isolates, with 45 of 54 isolates being identical at this locus. The complete absence of wPpe among some populations and low prevalence in others suggest that this endosymbiont is not an obligate mutualist of P. penetrans. The observed sex ratio bias in wPpe-infected nematode populations is similar to that observed in arthropods where Wolbachia acts as a reproductive manipulator, raising the question of a similar role in plant-parasitic nematodes.

Sex and Mitonuclear Adaptation in Experimental Caenorhabditis elegans Populations.

To reveal phenotypic and functional genomic patterns of mitonuclear adaptation, a laboratory adaptation study with Caenorhabditis elegans nematodes was conducted in which independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain (ETC) mutant, gas-1 Following 60 generations of evolution in large population sizes with competition for food resources, two distinct classes of lines representing different degrees of adaptive response emerged: a low-fitness class that exhibited minimal or no improvement compared to the gas-1 mutant ancestor, and a high-fitness class containing lines that exhibited partial recovery of wild-type fitness. Many lines that achieved higher reproductive and competitive fitness levels were also noted to evolve high frequencies of males during the experiment, consistent with adaptation in these lines having been facilitated by outcrossing. Whole-genome sequencing and analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome and could be classified into a small number of functional genomic categories. A highly nonrandom pattern of mitochondrial DNA mutation was observed within high-fitness gas-1 lines, with parallel fixations of nonsynonymous base substitutions within genes encoding NADH dehydrogenase subunits I and VI. These mitochondrial gene products reside within ETC complex I alongside the nuclear-encoded GAS-1 protein, suggesting that rapid adaptation of select gas-1 recovery lines was driven by fixation of compensatory mitochondrial mutations.

Genome Announcement: The Draft Genomes of Two Radopholus similis populations from Costa Rica.

Radopholus similis is an economically important pest of both banana and citrus in tropical regions. Here we present draft genomes from two populations of R. similis from Costa Rica that were created and assembled using short read libraries from Illumina HiSeq technology.

Comparative Genomics of Wolbachia-Cardinium Dual Endosymbiosis in a Plant-Parasitic Nematode.

Wolbachia and Cardinium are among the most important and widespread of all endosymbionts, occurring in nematodes and more than half of insect and arachnid species, sometimes as coinfections. These symbionts are of significant interest as potential biocontrol agents due to their abilities to cause major effects on host biology and reproduction through cytoplasmic incompatibility, sex ratio distortion, or obligate mutualism. The ecological and metabolic effects of coinfections are not well understood. This study examined a Wolbachia-Cardinium coinfection in the plant-parasitic nematode (PPN), Pratylenchus penetrans, producing the first detailed study of such a coinfection using fluorescence in situ hybridization (FISH), polymerase chain reaction (PCR), and comparative genomic analysis. Results from FISH and single-nematode PCR showed 123/127 individuals in a focal population carried Cardinium (denoted strain cPpe), and 48% were coinfected with Wolbachia strain wPpe. Both endosymbionts showed dispersed tissue distribution with highest densities in the anterior intestinal walls and gonads. Phylogenomic analyses confirmed an early place of cPpe and long distance from a sister strain in another PPN, Heterodera glycines, supporting a long history of both Cardinium and Wolbachia in PPNs. The genome of cPpe was 1.36 Mbp with 35.8% GC content, 1,131 predicted genes, 41% having no known function, and missing biotin and lipoate synthetic capacity and a plasmid present in other strains, despite having a slightly larger genome compared to other sequenced Cardinium. The larger genome revealed expansions of gene families likely involved in host-cellular interactions. More than 2% of the genes of cPpe and wPpe were identified as candidate horizontally transferred genes, with some of these from eukaryotes, including nematodes. A model of the possible Wolbachia-Cardinium interaction is proposed with possible complementation in function for pathways such as methionine and fatty acid biosynthesis and biotin transport.

Identification of candidate effector genes of Pratylenchus penetrans.

Pratylenchus penetrans is one of the most important species of root lesion nematodes (RLNs) because of its detrimental and economic impact in a wide range of crops. Similar to other plant-parasitic nematodes (PPNs), P. penetrans harbours a significant number of secreted proteins that play key roles during parasitism. Here, we combined spatially and temporally resolved next-generation sequencing datasets of P. penetrans to select a list of candidate genes aimed at the identification of a panel of effector genes for this species. We determined the spatial expression of transcripts of 22 candidate effectors within the oesophageal glands of P. penetrans by in situ hybridization. These comprised homologues of known effectors of other PPNs with diverse putative functions, as well as novel pioneer effectors specific to RLNs. It is noteworthy that five of the pioneer effectors encode extremely proline-rich proteins. We then combined in situ localization of effectors with available genomic data to identify a non-coding motif enriched in promoter regions of a subset of P. penetrans effectors, and thus a putative hallmark of spatial expression. Expression profiling analyses of a subset of candidate effectors confirmed their expression during plant infection. Our current results provide the most comprehensive panel of effectors found for RLNs. Considering the damage caused by P. penetrans, this information provides valuable data to elucidate the mode of parasitism of this nematode and offers useful suggestions regarding the potential use of P. penetrans-specific target effector genes to control this important pathogen.

Adaptive Evolution under Extreme Genetic Drift in Oxidatively Stressed Caenorhabditis elegans.

A mutation-accumulation (MA) experiment with Caenorhabditis elegans nematodes was conducted in which replicate, independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain mutant, gas-1. The original intent of the study was to assess the effect of electron transport chain dysfunction involving elevated reactive oxygen species production on patterns of spontaneous germline mutation. In contrast to results of standard MA experiments, gas-1 MA lines evolved slightly higher mean fitness alongside reduced among-line genetic variance compared with their ancestor. Likewise, the gas-1 MA lines experienced partial recovery to wildtype reactive oxygen species levels. Whole-genome sequencing and analysis revealed that the molecular spectrum but not the overall rate of nuclear DNA mutation differed from wildtype patterns. Further analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome, and could be classified into a small number of functional-genomic categories. Characterization of a backcrossed four-mutation set isolated from one gas-1 MA line revealed this combination to be beneficial on both gas-1 mutant and wildtype genetic backgrounds. Our combined results suggest that selection favoring beneficial mutations can be powerful even under unfavorable population genetic conditions, and agree with fitness landscape theory predicting an inverse relationship between population fitness and the likelihood of adaptation.

The Draft Genome of Globodera ellingtonae.

Globodera ellingtonae is a newly described potato cyst nematode (PCN) found in Idaho, Oregon, and Argentina. Here, we present a genome assembly for G. ellingtonae, a relative of the quarantine nematodes G. pallida and G. rostochiensis, produced using data from Illumina and Pacific Biosciences DNA sequencing technologies.

Genomic evidence for plant-parasitic nematodes as the earliest Wolbachia hosts.

Wolbachia, one of the most widespread endosymbionts, is a target for biological control of mosquito-borne diseases (malaria and dengue virus), and antibiotic elimination of infectious filarial nematodes. We sequenced and analyzed the genome of a new Wolbachia strain (wPpe) in the plant-parasitic nematode Pratylenchus penetrans. Phylogenomic analyses placed wPpe as the earliest diverging Wolbachia, suggesting two evolutionary invasions into nematodes. The next branches comprised strains in sap-feeding insects, suggesting Wolbachia may have first evolved as a nutritional mutualist. Genome size, protein content, %GC, and repetitive DNA allied wPpe with mutualistic Wolbachia, whereas gene repertoire analyses placed it between parasite (A, B) and mutualist (C, D, F) groups. Conservation of iron metabolism genes across Wolbachia suggests iron homeostasis as a potential factor in its success. This study enhances our understanding of this globally pandemic endosymbiont, highlighting genetic patterns associated with host changes. Combined with future work on this strain, these genomic data could help provide potential new targets for plant-parasitic nematode control.

Paternal Mitochondrial Transmission in Intra-Species Caenorhabditis briggsae Hybrids.

To study mitochondrial-nuclear genetic interactions in the nematode Caenorhabditis briggsae, our three laboratories independently created 38 intra-species cytoplasmic-nuclear hybrid (cybrid) lines. Although the cross design combines maternal mitotypes with paternal nuclear genotypes, eight lines (21%) unexpectedly contained paternal mitotypes. All eight share in common ancestry of one of two genetically related strains. This unexpected parallel observation of paternal mitochondrial transmission, undesirable given our intent of creating cybrids, provides a serendipitous experimental model and framework to study the molecular and evolutionary basis of uniparental mitochondrial inheritance.

The mitochondrial genome of Globodera ellingtonae is composed of two circles with segregated gene content and differential copy numbers.

BACKGROUND: The evolution of animal mitochondrial (mt) genomes has resulted in a highly conserved structure: a single compact circular chromosome approximately 14 to 20 kb long. Within the last two decades exceptions to this conserved structure, such as the division of the genome into multiple chromosomes, have been reported in a diverse set of metazoans. We report on the two circle multipartite mt genome of a newly described cyst nematode, Globodera ellingtonae. RESULTS: The G. ellingtonae mt genome was found to be comprised of two circles, each larger than any other multipartite circular mt chromosome yet reported, and both were larger than the single mt circle of the model nematode Caenorhabditis elegans. The genetic content of the genome was disproportionately divided between the two circles, although they shared a ~6.5 kb non-coding region. The 17.8 kb circle (mtDNA-I) contained ten protein-coding genes and two tRNA genes, whereas the 14.4 kb circle (mtDNA-II) contained two protein-coding genes, 20 tRNA genes and both rRNA genes. Perhaps correlated with this division of genetic content, the copy number of mtDNA-II was more than four-fold that of mtDNA-I in individual nematodes. The difference in copy number increased between second-stage and fourth-stage juveniles. CONCLUSIONS: The segregation of gene types to different mt circles in G. ellingtonae could provide benefit by localizing gene functional types to independent transcriptional units. This is the first report of both two-circle and several-circle mt genomes within a single genus. The differential copy number associated with this multipartite mt organization could provide a model system for deconstructing mechanisms regulating mtDNA copy number both in somatic cells and during germline development.

Paths of Heritable Mitochondrial DNA Mutation and Heteroplasmy in Reference and gas-1 Strains of Caenorhabditis elegans.

Heteroplasmy-the presence of more than one mitochondrial DNA (mtDNA) sequence type in a cell, tissue, or individual-impacts human mitochondrial disease and numerous aging-related syndromes. Understanding the trans-generational dynamics of mtDNA is critical to understanding the underlying mechanisms of mitochondrial disease and evolution. We investigated mtDNA mutation and heteroplasmy using a set of wild-type (N2 strain) and mitochondrial electron transport chain (ETC) mutant (gas-1) mutant Caenorhabditis elegans mutation-accumulation (MA) lines. The N2 MA lines, derived from a previous experiment, were bottlenecked for 250 generations. The gas-1 MA lines were created for this study, and bottlenecked in the laboratory for up to 50 generations. We applied Illumina-MiSeq DNA sequencing to L1 larvae from five gas-1 MA lines and five N2 MA lines to detect and characterize mtDNA mutation and heteroplasmic inheritance patterns evolving under extreme drift. mtDNA copy number increased in both sets of MA lines: three-fold on average among the gas-1 MA lines and five-fold on average among N2 MA lines. Eight heteroplasmic single base substitution polymorphisms were detected in the gas-1 MA lines; only one was observed in the N2 MA lines. Heteroplasmy frequencies ranged broadly in the gas-1 MA lines, from as low as 2.3% to complete fixation (homoplasmy). An initially low-frequency (<5%) heteroplasmy discovered in the gas-1 progenitor was observed to fix in one gas-1 MA line, achieve higher frequency (37.4%) in another, and be lost in the other three lines. A similar low-frequency heteroplasmy was detected in the N2 progenitor, but was lost in all five N2 MA lines. We identified three insertion-deletion (indel) heteroplasmies in gas-1 MA lines and six indel variants in the N2 MA lines, most occurring at homopolymeric nucleotide runs. The observed bias toward accumulation of single nucleotide polymorphisms in gas-1 MA lines is consistent with the idea that impaired mitochondrial activity renders mtDNA more vulnerable to this type of mutation. The consistent increases in mtDNA copy number implies that extreme genetic drift provides a permissive environment for elevated organelle genome copy number in C. elegans reference and gas-1 strains. This study broadens our understanding of the heteroplasmic mitochondrial mutation process in a multicellular model organism.