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.

Pestiferous slugs and their associated nematodes in agricultural fields, greenhouses, and nurseries in Alberta, Canada

Some slug species are considered a nuisance in agriculture and horticulture worldwide, causing economic losses to growers. Phasmarhabditis is a genus of bacteria-feeding nematodes that can parasitize slugs and snails and thus potentially serve as a biological control agent. Canada had no record of Phasmarhabditis until a survey conducted in 2019 reported a Canadian strain of Phasmarhabditis californica from a single Arion rufus slug. To build on this discovery, we surveyed three major agricultural sites, ten greenhouses, and nurseries in Alberta from June to September 2021 to collect pest slug species and investigate their associated nematodes, specifically P. californica. Slugs were collected from the field and returned to the laboratory to check for emerging nematodes on White traps. We collected 1331 slugs belonging to nine species, with Deroceras reticulatum being the most common. Only 45 (3.38%) slug samples were positive for nematodes, and the majority were identified to species level: Alloionema appendiculatum, Caenorhabditis briggsae, Caenorhabditis elegans, Panagrolaimus subelongatus, and Mesorhabditis spiculigera. We did not isolate P. californica from any of the slugs collected from these survey sites, which included the original site where P. californica was discovered. However, four D. reticulatum slugs retrieved from a residential garden sample were infected with P. californica. These findings suggest the possibility of a fragmented distribution of P. californica across Alberta. Future research should focus on extensively surveying agriculture and horticulture sites and residential gardens in different provinces across Canada.

First Report of the Gastropod-Killing Nematode, Phasmarhabditis californica, in Washington State, U.S.A

Phasmarhabditis californica, a commercially available biological control agent in England, Scotland, and Wales (Nemaslug 2.0 ®), was discovered for the first time in Washington State during 2022. Nematodes were recovered from the invasive gastropods Arion hortensis, Deroceras reticulatum, and Oxychilus sp. in garden centers in both Vancouver and Spokane. The 18S ribosomal RNA gene was used to identify the nematodes. This discovery builds upon previous reports of P. californica in California and Oregon and suggests that the species is widespread throughout the west coast of the U.S. Future research directions with P. californica are proposed.

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.

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.

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.

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.

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.

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.

Analysis of nematode-endosymbiont coevolution in the Xiphinema americanum species complex using molecular markers of variable evolutionary rates

Bacterial symbioses play important roles in shaping diverse biological processes in nematodes, and serve as targets in nematode biocontrol strategies. Focusing on the Xiphinema americanum species complex, we expanded upon recent research investigating patterns of coevolution between Xiphinema spp. and Xiphinematobacter spp., utilising two symbiont genetic markers of varying evolutionary rates. Phylogenetic analysis of nematode mitochondrial DNA (mtDNA) revealed five strongly supported major clades. Analysis of slow-evolving 16S rDNA in bacterial symbionts resulted in a phylogenetic topology composed of four major clades that grouped taxa highly congruent with the nematode mtDNA topology. A faster evolving protein-coding symbiont gene (nad) provided more phylogenetic resolution with seven well-supported clades, also congruent with the nematode mtDNA tree topology. Our results reinforce recent studies suggesting extensive coevolution between Xiphinema spp. and their vertically transmitted endosymbionts Xiphinematobacter spp. and illustrate the advantages of including genetic markers of varying evolutionary rates in coevolutionary and phylogenetic studies.

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.

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.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3047-x) contains supplementary material, which is available to authorized users.

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.

Genome-wide polymorphism and signatures of selection in the symbiotic sea anemone Aiptasia

Background

Coral reef ecosystems are declining in response to global climate change and anthropogenic impacts. Yet patterns of standing genetic variation within cnidarian species, a major determinant of adaptive potential, are virtually unknown at genome-scale resolution. We explore patterns of genome-wide polymorphism and identify candidate loci under selection in the sea anemone Aiptasia, an important laboratory model system for studying the symbiosis between corals and dinoflagellate algae of the genus Symbiodinium.

Results

Low coverage genome sequencing revealed large genetic distances among globally widespread lineages, novel candidate targets of selection, and considerably higher heterozygosity than previously reported for Aiptasia. More than 670,000 single nucleotide polymorphisms were identified among 10 Aiptasia individuals including two pairs of genetic clones. Evolutionary relationships based on genome-wide polymorphism supported the current paradigm of a genetically distinct population from the US South Atlantic that harbors diverse Symbiodinium clades. However, anemones from the US South Atlantic demonstrated a striking lack of shared derived polymorphism. Heterozygosity was an important feature shaping nucleotide diversity patterns: at any given SNP site, more than a third of individuals genotyped were heterozygotes, and heterozygosity within individual genomes ranged from 0.37–0.58 %. Analysis of nonsynonymous and synonymous sites suggested that highly heterozygous regions are evolving under relaxed purifying selection compared to the rest of the Aiptasia genome. Genes previously identified as having elevated evolutionary rates in Aiptasia compared to other cnidarians were found in our study to be under strong purifying selection within Aiptasia. Candidate targets of selection, including lectins and genes involved in Rho GTPase signalling, were identified based on unusual signatures of nucleotide diversity, Tajima’s D, and heterozygosity compared to genome-wide averages.

Conclusions

This study represents the first genome-wide analysis of Tajima’s D in a cnidarian. Our results shed light on patterns of intraspecific genome-wide polymorphism in a model for studies of coral-algae symbiosis and present genetic targets for future research on evolutionary and cellular processes in early-diverging metazoans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2488-6) contains supplementary material, which is available to authorized users.

Equine Protozoal Myeloencephalitis: An Updated Consensus Statement with a Focus on Parasite Biology, Diagnosis, Treatment, and Prevention

Equine protozoal myeloencephalitis (EPM) remains an important neurologic disease of horses. There are no pathognomonic clinical signs for the disease. Affected horses can have focal or multifocal central nervous system (CNS) disease. EPM can be difficult to diagnose antemortem. It is caused by either of 2 parasites, Sarcocystis neurona and Neospora hughesi, with much less known about N. hughesi. Although risk factors such as transport stress and breed and age correlations have been identified, biologic factors such as genetic predispositions of individual animals, and parasite‐specific factors such as strain differences in virulence, remain largely undetermined. This consensus statement update presents current published knowledge of the parasite biology, host immune response, disease pathogenesis, epidemiology, and risk factors. Importantly, the statement provides recommendations for EPM diagnosis, treatment, and prevention.

Comparative Genomics of a Plant-Parasitic Nematode Endosymbiont Suggest a Role in Nutritional Symbiosis

Bacterial mutualists can modulate the biochemical capacity of animals. Highly coevolved nutritional mutualists do this by synthesizing nutrients missing from the host’s diet. Genomics tools have advanced the study of these partnerships. Here we examined the endosymbiont Xiphinematobacter (phylum Verrucomicrobia) from the dagger nematode Xiphinema americanum, a migratory ectoparasite of numerous crops that also vectors nepovirus. Previously, this endosymbiont was identified in the gut, ovaries, and eggs, but its role was unknown. We explored the potential role of this symbiont using fluorescence in situ hybridization, genome sequencing, and comparative functional genomics. We report the first genome of an intracellular Verrucomicrobium and the first exclusively intracellular non-Wolbachia nematode symbiont. Results revealed that Xiphinematobacter had a small 0.916-Mb genome with only 817 predicted proteins, resembling genomes of other mutualist endosymbionts. Compared with free-living relatives, conserved proteins were shorter on average, and there was large-scale loss of regulatory pathways. Despite massive gene loss, more genes were retained for biosynthesis of amino acids predicted to be essential to the host. Gene ontology enrichment tests showed enrichment for biosynthesis of arginine, histidine, and aromatic amino acids, as well as thiamine and coenzyme A, diverging from the profiles of relatives Akkermansia muciniphilia (in the human colon), Methylacidiphilum infernorum, and the mutualist Wolbachia from filarial nematodes. Together, these features and the location in the gut suggest that Xiphinematobacter functions as a nutritional mutualist, supplementing essential nutrients that are depleted in the nematode diet. This pattern points to evolutionary convergence with endosymbionts found in sap-feeding insects.

Selfish Mitochondrial DNA Proliferates and Diversifies in Small, but not Large, Experimental Populations of Caenorhabditis briggsae

Evolutionary interactions across levels of biological organization contribute to a variety of fundamental processes including genome evolution, reproductive mode transitions, species diversification, and extinction. Evolutionary theory predicts that so-called “selfish” genetic elements will proliferate when the host effective population size (N_e) is small, but direct tests of this prediction remain few. We analyzed the evolutionary dynamics of deletion-containing mitochondrial DNA (ΔmtDNA) molecules, previously characterized as selfish elements, in six different natural strains of the nematode Caenorhabditis briggsae allowed to undergo experimental evolution in a range of population sizes (N = 1, 10, 100, and 1,000) for a maximum of 50 generations. Mitochondrial DNA (mtDNA) was analyzed for replicate lineages at each five-generation time point. Ten different ΔmtDNA molecule types were observed and characterized across generations in the experimental populations. Consistent with predictions from evolutionary theory, lab lines evolved in small-population sizes (e.g., nematode N = 1) were more susceptible to accumulation of high levels of preexisting ΔmtDNA compared with those evolved in larger populations. New ΔmtDNA elements were observed to increase in frequency and persist across time points, but almost exclusively at small population sizes. In some cases, ΔmtDNA levels decreased across generations when population size was large (nematode N = 1,000). Different natural strains of C. briggsae varied in their susceptibilities to ΔmtDNA accumulation, owing in part to preexisting compensatory mtDNA alleles in some strains that prevent deletion formation. This analysis directly demonstrates that the evolutionary trajectories of ΔmtDNA elements depend upon the population-genetic environments and molecular-genetic features of their hosts.

Sarcocyst Development in Raccoons (Procyon lotor) Inoculated with Different Strains of Sarcocystis neurona Culture-Derived Merozoites

Sarcocystis neurona is considered the major etiologic agent of equine protozoal myeloencephalitis (EPM), a neurological disease in horses. Raccoon ( Procyon lotor ) is considered the most important intermediate host in the life cycle of S. neurona in the United States; S. neurona sarcocysts do mature in raccoon muscles, and raccoons also develop clinical signs simulating EPM. The focus of this study was to determine if sarcocysts would develop in raccoons experimentally inoculated with different host-derived strains of in vitro-cultivated S. neurona merozoites. Four raccoons were inoculated with strains derived from a raccoon, a sea otter, a cat, and a horse. Raccoon tissues were fed to laboratory-raised opossums ( Didelphis virginiana ), the definitive host of S. neurona . Intestinal scraping revealed sporocysts in opossums who received muscle tissue from raccoons inoculated with the raccoon-derived or the sea otter-derived isolates. These results demonstrate that sarcocysts can mature in raccoons inoculated with in vitro-derived S. neurona merozoites. In contrast, the horse and cat-derived isolates did not produce microscopically or biologically detected sarcocysts. Immunoblot analysis revealed both antigenic and antibody differences when testing the inoculated raccoons. Immunohistochemical staining indicated differences in staining between the merozoite and sarcocyst stages. The successful infections achieved in this study indicates that the life cycle can be manipulated in the laboratory without affecting subsequent stage development, thereby allowing further purification of strains and artificial maintenance of the life cycle.

An update on Sarcocystis neurona infections in animals and equine protozoal myeloencephalitis (EPM)

Equine protozoal myeloencephalitis (EPM) is a serious disease of horses, and its management continues to be a challenge for veterinarians. The protozoan Sarcocystis neurona is most commonly associated with EPM. S. neurona has emerged as a common cause of mortality in marine mammals, especially sea otters (Enhydra lutris). EPM-like illness has also been recorded in several other mammals, including domestic dogs and cats. This paper updates S. neurona and EPM information from the last 15 years on the advances regarding life cycle, molecular biology, epidemiology, clinical signs, diagnosis, treatment and control.

In vivo quantification reveals extensive natural variation in mitochondrial form and function in Caenorhabditis briggsae

We have analyzed natural variation in mitochondrial form and function among a set of Caenorhabditis briggsae isolates known to harbor mitochondrial DNA structural variation in the form of a heteroplasmic nad5 gene deletion (nad5Δ) that correlates negatively with organismal fitness. We performed in vivo quantification of 24 mitochondrial phenotypes including reactive oxygen species level, membrane potential, and aspects of organelle morphology, and observed significant among-isolate variation in 18 traits. Although several mitochondrial phenotypes were non-linearly associated with nad5Δ levels, most of the among-isolate phenotypic variation could be accounted for by phylogeographic clade membership. In particular, isolate-specific mitochondrial membrane potential was an excellent predictor of clade membership. We interpret this result in light of recent evidence for local adaptation to temperature in C. briggsae. Analysis of mitochondrial-nuclear hybrid strains provided support for both mtDNA and nuclear genetic variation as drivers of natural mitochondrial phenotype variation. This study demonstrates that multicellular eukaryotic species are capable of extensive natural variation in organellar phenotypes and highlights the potential of integrating evolutionary and cell biology perspectives.

Variation in base-substitution mutation in experimental and natural lineages of Caenorhabditis nematodes

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species.

Exposure to Sarcocystis spp. in horses from Spain determined by Western blot analysis using Sarcocystis neurona merozoites as heterologous antigen

Horses serve as an intermediate host for several species of Sarcocystis, all of which utilize canids as the definitive host. Sarcocystis spp. infection and formation of latent sarcocysts in horses often appears to be subclinical, but morbidity can occur, especially when the parasite burden is large. A serological survey was conducted to determine the presence of antibodies against Sarcocystis spp. in seemingly healthy horses from the Galicia region of Spain. Western blot analyses using Sarcocystis neurona merozoites as heterologous antigen suggested greater than 80% seroprevalance of Sarcocystis spp. in a sample set of 138 horses. The serum samples were further tested with enzyme-linked immunosorbent assays (ELISAs) based on recombinant S. neurona-specific surface antigens (rSnSAGs). As expected for horses from the Eastern Hemisphere, less than 4% of the serum samples were positive when analyzed with either the rSnSAG2 or the rSnSAG4/3 ELISAs. An additional 246 horses were tested using the rSnSAG2 ELISA, which revealed that less than 3% of the 384 samples were seropositive. Collectively, the results of this serologic study suggested that a large proportion of horses from this region of Spain are exposed to Sarcocystis spp. Furthermore, the anti-Sarcocystis seroreactivity in these European horses could be clearly distinguished from anti-S. neurona antibodies using the rSnSAG2 and rSnSAG4/3 ELISAs.

The SnSAG merozoite surface antigens of Sarcocystis neurona are expressed differentially during the bradyzoite and sporozoite life cycle stages

Sarcocystis neurona is a two-host coccidian parasite whose complex life cycle progresses through multiple developmental stages differing at morphological and molecular levels. The S. neurona merozoite surface is covered by multiple, related glycosylphosphatidylinositol-linked proteins, which are orthologous to the surface antigen (SAG)/SAG1-related sequence (SRS) gene family of Toxoplasma gondii. Expression of the SAG/SRS proteins in T. gondii and another related parasite Neospora caninum is life-cycle stage specific and seems necessary for parasite transmission and persistence of infection. In the present study, the expression of S. neurona merozoite surface antigens (SnSAGs) was evaluated in the sporozoite and bradyzoite stages. Western blot analysis was used to compare SnSAG expression in merozoites versus sporozoites, while immunocytochemistry was performed to examine expression of the SnSAGs in merozoites versus bradyzoites. These analyses revealed that SnSAG2, SnSAG3 and SnSAG4 are expressed in sporozoites, while SnSAG5 was appeared to be downregulated in this life cycle stage. In S. neurona bradyzoites, it was found that SnSAG2, SnSAG3, SnSAG4 and SnSAG5 were either absent or expression was greatly reduced. As shown for T. gondii, stage-specific expression of the SnSAGs may be important for the parasite to progress through its developmental stages and complete its life cycle successfully. Thus, it is possible that the SAG switching mechanism by these parasites could be exploited as a point of intervention. As well, the alterations in surface antigen expression during different life cycle stages may need to be considered when designing prospective approaches for protective vaccination.

Detection of antibodies against Sarcocystis neurona, Neospora spp., and Toxoplasma gondii in horses from Costa Rica

Serum samples from 315 horses from Costa Rica, Central America, were examined for the presence of antibodies against Sarcocystis neurona, Neospora spp., and Toxoplasma gondii by using the surface antigen (SAG) SnSAG2 enzyme-linked immunosorbent assay (ELISA), the NhSAG1 ELISA, and the modified agglutination test, respectively. Anti- S. neurona antibodies were found in 42.2% of the horses by using the SnSAG2 ELISA. Anti- Neospora spp. antibodies were found in only 3.5% of the horses by using the NhSAG1 ELISA, and only 1 of these horses was confirmed seropositive by Western blot. Antibodies to T. gondii were found in 34.0% of the horses tested, which is higher than in previous reports from North and South America. The finding of anti- S. neurona antibodies in horses from geographical areas where Didelphis marsupialis has wide distribution suggests that D. marsupialis is a potential definitive host for this parasite and a source of infection for these horses.

Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans

Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ∼ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.

High rate of large deletions in Caenorhabditis briggsae mitochondrial genome mutation processes

Mitochondrial DNA (mtDNA) mutations underlie a variety of human genetic disorders and are associated with the aging process. mtDNA polymorphisms are widely used in a variety of evolutionary applications. Although mtDNA mutation spectra are known to differ between distantly related model organisms, the extent to which mtDNA mutation processes vary between more closely related species and within species remains enigmatic. We analyzed mtDNA divergence in two sets of 250-generation Caenorhabditis briggsae mutation-accumulation (MA) lines, each derived from a different natural isolate progenitor: strain HK104 from Okayama, Japan, and strain PB800 from Ohio, United States. Both sets of C. briggsae MA lines accumulated numerous large heteroplasmic mtDNA deletions, whereas only one similar event was observed in a previous analysis of Caenorhabditis elegans MA line mtDNA. Homopolymer length change mutations were frequent in both sets of C. briggsae MA lines and occurred in both intergenic and protein-coding gene regions. The spectrum of C. briggsae mtDNA base substitution mutations differed from the spectrum previously observed in C. elegans. In C. briggsae, the HK104 MA lines experienced many different base substitution types, whereas the PB800 lines displayed only C:G --> T:A transitions, although the difference was not significant. Over half of the mtDNA base substitutions detected in the C. briggsae MA lines were in a heteroplasmic state, whereas all those previously characterized in C. elegans MA line mtDNA were fixed changes, indicating a narrower mtDNA bottleneck in C. elegans as compared with C. briggsae. Our results show that C. briggsae mtDNA is highly susceptible to large deletions and that the mitochondrial mutation process varies between Caenorhabditis nematode species.

Evolution of Caenorhabditis mitochondrial genome pseudogenes and Caenorhabditis briggsae natural isolates

Although most metazoan mitochondrial genomes are highly streamlined and encode little noncoding DNA outside of the "AT" region, the accumulation of mitochondrial pseudogenes and other types of noncoding DNA has been observed in a growing number of animal groups. The nematode species Caenorhabditis briggsae harbors two mitochondrial DNA (mtDNA) pseudogenes, named Psinad5-1 and Psinad5-2, presumably derived from the nad5 protein-coding gene. Here, we provide an in-depth analysis of mtDNA pseudogene evolution in C. briggsae natural isolates and related Caenorhabditis species. Mapping the observed presence and absence of the pseudogenes onto phylogenies suggests that Psinad5-1 originated in the ancestor to C. briggsae and its recently discovered outcrossing relative species Caenorhabditis sp. 5 and Caenorhabditis sp. 9. However, Psinad5-1 was not detected in Caenorhabditis sp. 9 natural isolates, suggesting a lineage-specific loss of this pseudogene in this species. Our results corroborated the previous finding that Psinad5-2 originated within C. briggsae. The observed pattern of mitochondrial pseudogene gain and loss in Caenorhabditis was inconsistent with predictions of the tandem duplication-random loss model of mitochondrial genome evolution and suggests that intralineage recombination-like mechanisms might play a major role in Caenorhabditis mtDNA evolution. Natural variation was analyzed at the pseudogenes and flanking mtDNA sequences in 141 geographically diverse C. briggsae natural isolates. Although phylogenetic analysis placed the majority of isolates into the three previously established major intraspecific clades of C. briggsae, two new and unexpected haplotypes fell outside of these conventional groupings. Psinad5-2 copy number variation was observed among C. briggsae isolates collected from the same geographic site. Patterns of nucleotide diversity were analyzed in Psinad5-1 and Psinad5-2, and confidence intervals were found to overlap values from synonymous sites in protein-coding genes, consistent with neutral expectations. Our findings provide new insights into the mode and tempo of mitochondrial genome and pseudogene evolution both within and between Caenorhabditis nematode species.

Muller's Ratchet and compensatory mutation in Caenorhabditis briggsae mitochondrial genome evolution

Background

The theory of Muller' Ratchet predicts that small asexual populations are doomed to accumulate ever-increasing deleterious mutation loads as a consequence of the magnified power of genetic drift and mutation that accompanies small population size. Evidence for Muller's Ratchet and knowledge on its underlying molecular mechanisms, however, are lacking for natural populations.

Results

We characterized mitochondrial genome evolutionary processes in Caenorhabditis briggsae natural isolates to show that numerous lineages experience a high incidence of nonsynonymous substitutions in protein-coding genes and accumulate unusual deleterious noncoding DNA stretches with associated heteroplasmic function-disrupting genome deletions. Isolate-specific deletion proportions correlated negatively with nematode fecundity, suggesting that these deletions might negatively affect C. briggsae fitness. However, putative compensatory mutations were also observed that are predicted to reduce heteroplasmy levels of deleterious deletions. Paradoxically, compensatory mutations were observed in one major intraspecific C. briggsae clade where population sizes are estimated to be very small (and selection is predicted to be relatively weak), but not in a second major clade where population size estimates are much larger and selection is expected to be more efficient.

Conclusion

This study provides evidence that the mitochondrial genomes of animals evolving in nature are susceptible to Muller's Ratchet, suggests that context-dependent compensatory mutations can accumulate in small populations, and predicts that Muller's Ratchet can affect fundamental evolutionary forces such as the rate of mutation.

Soybean bacterial artificial chromosome contigs anchored with RFLPs: insights into genome duplication and gene clustering

Surveying the soybean genome with 683 bacterial artificial chromosome (BAC) contiguous groups (contigs) anchored by restriction fragment length polymorphisms (RFLPs) enabled us to explore microsyntenic relationships among duplicated regions and also to examine the physical organization of hypomethylated (and presumably gene-rich) genomic regions. Numerous cases where nonhomologous RFLPs hybridized to common BAC clones indicated that RFLPs were physically clustered in soybean, apparently in less than 25% of the genome. By extension, we speculate that most of the genes are clustered in less than 275 M of the soybean genome. Approximately 40%-45% of this gene-rich portion is associated with the RFLP-anchored contigs described in this study. Similarities in genome organization among BAC contigs from duplicate genomic regions were also examined. Homoeologous BAC contigs often exhibited extensive microsynteny. Furthermore, paralogs recovered from duplicate contigs shared 86%-100% sequence identity.

Redescription of Neospora caninum and its differentiation from related coccidia

Neospora caninum is a protozoan parasite of animals, which before 1984 was misidentified as Toxoplasma gondii. Infection by this parasite is a major cause of abortion in cattle and causes paralysis in dogs. Since the original description of N. caninum in 1988, considerable progress has been made in the understanding of its life cycle, biology, genetics and diagnosis. In this article, the authors redescribe the parasite, distinguish it from related coccidia, and provide accession numbers to its type specimens deposited in museums.

Characterization of the Oregon isolate of Neospora hughesi from a horse

Neospora hughesi was isolated in cell cultures inoculated with homogenate of spinal cord from a horse in Oregon. Tachyzoites of this Oregon isolate of N. hughesi were maintained continuously by cell culture passage and tachyzoites were infective to immunosuppressed mice. Gamma interferon gene knockout (KO) mice injected with tachyzoites developed fatal myocarditis and numerous tachyzoites were seen in lesions. Gerbils (Meriones unguiculatus) inoculated with tachyzoites developed antibodies (> or = 1:500) as indicated by the Neospora caninum agglutination test but did not develop clinical signs, and Neospora organisms were not demonstrable in their tissues. Tissue cysts were not found in gerbils, nude mice, KO mice, immunosuppressed outbred Swiss Webster mice, or BALB/c mice injected with the Oregon isolate of N. hughesi. Ultrastructurally, tachyzoites of the Oregon isolate from the myocardium of infected KO mice and from cell culture were similar to N. caninum tachyzoites. Western blot analysis using NcSAG1 and NcSRS2 polyclonal and monoclonal antibodies and characterization of the internal transcribed spacer 1 sequences from the equine isolates and different isolates of N. caninum from dogs and cattle indicated that the Oregon isolate of N. hughesi is distinct from N. caninum isolates from cattle and dogs.

Initiation of a Sarcocystis neurona expressed sequence tag (EST) sequencing project: a preliminary report

To accelerate genetic and molecular characterization of Sarcocystis neurona, the primary causative agent of equine protozoal myeloencephalitis (EPM), a sequencing project has been initiated that will generate approximately 7000-8000 expressed sequence tags (ESTs) from this apicomplexan parasite. Poly(A)(+) RNA was isolated from culture-derived S. neurona merozoites, and a cDNA library was constructed in a unidirectional lambda phage cloning vector. Sixty phage clones were randomly picked from the library, and the cDNA inserts were amplified from these clones using the T3 and T7 primers that flank the multi-cloning site of the lambda vector. This analysis demonstrated that 100% (60/60) of the clones selected from this library contained recombinant cDNA inserts ranging in size from 0.4 to 4.0 kilobases (kb) with an average size of 1.23kb. Single-pass sequencing from the 5' end of the 60 amplified cDNAs produced high-quality nucleotide sequence from 53 of the clones. Comparison of these ESTs to the current gene databases revealed significant matches for 10 of the ESTs, six of which are similar to sequences from other Apicomplexa (i.e., Toxoplasma gondii). Importantly, none of the ESTs were of obvious mammalian origin, thus indicating that the cDNAs in this library were derived primarily from parasite mRNA and not from mRNA of the bovine turbinate host cells. Collectively, these data indicate that the described cDNA library will provide an excellent substrate for generating a portion of the ESTs that are planned from S. neurona. This sequencing project will greatly hasten gene discovery for this protozoan pathogen thereby enhancing efforts towards the development of improved diagnostics, treatments, and preventatives for EPM. In addition, the S. neurona ESTs will represent a significant contribution to the extensive database of sequences from the Apicomplexa. Comparative analyses of these apicomplexan sequences will likely offer a multitude of important information about the biology and evolutionary history of this phylogenetic grouping of parasites.

Biological and molecular characterizations of Toxoplasma gondii strains obtained from southern sea otters (Enhydra lutris nereis)

Toxoplasma gondii was isolated from brain or heart tissue from 15 southern sea otters (Enhydra lutris nereis) in cell cultures. These strains were used to infect mice that developed antibodies to T. gondii as detected in the modified direct agglutination test and had T. gondii tissue cysts in their brains at necropsy. Mouse brains containing tissue cysts from 4 of the strains were fed to 4 cats. Two of the cats excreted T. gondii oocysts in their feces that were infectious for mice. Molecular analyses of 13 strains indicated that they were all type II strains, but that they were genetically distinct from one another.

Molecular characterization of a thrombospondin-related anonymous protein homologue in Neospora caninum

Thrombospondin-related anonymous protein (TRAP) family members participate in attachment and invasion of host cells by apicomplexan parasites. A TRAP homologue in Neospora caninum strain Nc-1 (NcMIC2) was cloned, sequenced and found to be 61% identical (75% similar) at the amino acid level to Toxoplasma gondii MIC2 (TgMIC2). Similar to TgMIC2, the predicted amino acid sequence of NcMIC2 contains one integrin-like domain (I or A domain), five thrombospondin (TSP) repeats, a putative transmembrane spanning region and intracellular C-terminus, and was localized to micronemes by cryo-immunoelectron microscopy. The secretion of NcMIC2 was temperature dependent and was induced at or above 25 degrees C. The secreted form of NcMIC2 released into the medium was found to be proteolytically processed such that it lacked the C-terminal domain. Secretion of NcMIC2 was regulated by calcium, since several agents which raise intracellular calcium levels were shown to promote NcMIC2 secretion and chelation of [Ca(2+)](i) abrogated release. As a member of the growing family of apicomplexan TRAP proteins, NcMIC2 may play an important role in attachment and invasion by N. caninum into host cells.

Experimental approaches to understanding virulence in toxoplasmosis

Toxoplasma gondii is a widespread protozoan parasite that causes severe disease only in immunocompromised individuals. Equipped with excellent animal models and relatively advanced systems for genetics, T. gondii provides an excellent system for understanding pathogenesis. Resistance to toxoplasmosis is governed by rapid innate and adaptive immunity that is characterized by a Th1 type profile of cytokines. Despite this effective response, acute infections can cause considerable damage and the parasite effectively establishes a long-term chronic infection that predisposes the host to reactivation and provides a means of eventual transmission. This complex interaction is brought about by the differentiation of the parasite from a rapidly replicating, lytic form (known as the tachyzoite) to a slow-growing form (known as the bradyzoite) that gives rise to chronic infection. The population structure of T. gondii is remarkably clonal, consisting of just three predominant lineages that are geographically widespread and found in a variety of hosts including humans. Acute virulence is strongly associated with the type I genotype which exhibits an enhanced replication rate in vitro and higher tissue burdens in vivo relative to non-virulent lineages. The pathology associated with acute infection appears to be due to excessive production of acute inflammatory mediators, suggesting that disease is partly due to over-response of the host immune system. A combination of refined animal models and newly developed genetic tools for establishing the relative contribution of genes to pathogenesis will enable a comprehensive analysis of the molecular basis of virulence in toxoplasmosis.

Differentiation of Neospora hughesi from Neospora caninum based on their immunodominant surface antigen, SAG1 and SRS2

Neospora hughesi is a newly recognised parasite that is closely related to Neospora caninum, and is a cause of equine protozoal myeloencephalitis. We have characterised two N. hughesi immunodominant tachyzoite antigens which exhibit antigenic and molecular differences from the homologous tachyzoite antigens on N. caninum. These antigens on N. hughesi are referred to as NhSAG1 and NhSRS2, using the same mnemonics as used for the N. caninum antigens (NcSAG1 and NcSRS2), and are homologous to Toxoplasma gondii surface antigen 1 (SAG1) and SAG1-related sequence 2 (SRS2). The NcSAG1 and NcSRS2 were antigenically conserved in six different N. caninum isolates from cattle and dogs. The two equine-derived Neospora isolates, one designated as N. hughesi, were similar to each other but different from N. caninum. There was 6% difference in amino acid identity between NcSAG1 and NhSAG1, whereas there was a 9% difference when NcSRS2 and NhSRS2 were compared. The polymorphism of these genes and their corresponding proteins provide additional markers which can be used to distinguish N. caninum from N. hughesi.

Comparison of the major antigens of Neospora caninum and Toxoplasma gondii

The Apicomplexa are a diverse group of parasitic protozoa with very ancient phylogenetic roots. Consistent with their phylogeny, the extant species share conserved proteins and traits that were found in their apicomplexan progenitor, but at the same time they have diverged to occupy different biological niches (e.g. host-range and cell type). Characterisation of gene and protein diversity is important for distinguishing between related parasites, for determining their phylogeny, and for providing insight into factors that determine host restriction, cell preference, and virulence. The value of molecular characterisations and comparisons between species is well illustrated by the close phylogenetic relationship between Neospora caninum and Toxoplasma gondii. These two organisms have nearly identical morphology and can cause similar pathology and disease. Consequently, N. caninum has often been incorrectly identified as T. gondii, thus demonstrating the need for studies addressing the molecular and antigenic composition of Neospora. In this review, we describe the major antigenic proteins that have been characterised in N. caninum. These show homology to T. gondii proteins, yet possess unique antigenic characteristics that distinguish them from their homologues and enable their use for specific serological diagnoses and parasite identification.

Neospora caninum: tachyzoites express a potent type-I nucleoside triphosphate hydrolase

We have identified type I nucleoside triphosphate hydrolase (NTPase; EC 3.6. 1.3) activity, previously thought to be restricted to the virulent strains of Toxoplasma gondii, in the cell extracts of Neospora caninum tachyzoites. Sequence analysis of a complete cDNA from Nc-1 strain indicated that N. caninum NTPases shared approximately 69% identity to the NTPases of T. gondii and are most similar to the NTPase-I isozyme. Southern blot analysis of genomic DNA and sequence analysis of two independent NTP clones from the Nc-1 strain revealed the presence of multiple genes, at least two of which are transcribed. Substrate specificity and Km values for MgATP2- and MgADP- hydrolysis for recombinant or partially purified native NcNTPase were the same as those for the type I isozyme (NTPase-I). Significantly, no type II enzyme (NTPase-II) activity for NDP hydrolysis was detected in cell extracts of N. caninum, although it is universally present in all T. gondii strains that have been tested. This intriguing difference between these two closely related apicomplexan parasites may provide insight into the function of the NTPases during intracellular parasitism.

The p29 and p35 immunodominant antigens of Neospora caninum tachyzoites are homologous to the family of surface antigens of Toxoplasma gondii

Neospora caninum is an apicomplexan parasite that is closely related to Toxoplasma gondii and has been found to be associated with neurological disorders in dogs and congenital infections and abortions in cattle. We have identified two surface proteins of 29 and 35 kDa (designated Ncp29 and Ncp35, respectively) from N. caninum tachyzoites that are the predominant antigens recognized by antisera from Neospora-infected animals. Monoclonal antibodies against Ncp29 and Ncp35 were used to analyze several independent and diverse N. caninum isolates; both antigens were recognized in all isolates, suggesting that they are well conserved. Localization studies and surface labeling with biotin demonstrated that Ncp29 and Ncp35 are membrane associated and displayed on the surface of the parasite. After treatment with phosphatidylinositol-specific phospholipase C, parasite lysates were analyzed with antibodies against the cross-reacting determinant of glycosylphosphatidylinositol anchors. Approximately six glycolipid-anchored surface proteins were identified, with the two most prominent corresponding to Ncp29 and Ncp35. Sequence comparisons of Ncp29 and Ncp35 with GenBank indicated that they are most similar to the T. gondii surface antigen 1 (SAG1) and surface antigen 1-related sequence 2 (SRS2), respectively. Consequently, Ncp29 has been designated NcSAG1 and Ncp35 has been designated NcSRS2. Both NcSAG1 and NcSRS2 contain a tandemly duplicated motif and 12 absolutely conserved cysteines which are also found in all of the SAG and SRS proteins of T. gondii. Maintenance of these motifs and the 12 cysteine residues suggests that these surface antigens share a similar secondary and tertiary structure that is presumably important for a conserved function that these antigens serve during infection.

Targeted disruption of the GRA2 locus in Toxoplasma gondii decreases acute virulence in mice

Following invasion into the host cell, the protozoan Toxoplasma gondii secretes a variety of proteins that modify the parasitophorous vacuole. Within the vacuole, the 28-kDa dense granule protein known as GRA2 is specifically targeted to the tubulovesicular network which forms connections with the vacuolar membrane. To investigate the importance of GRA2, we derived from strain RH a mutant T. gondii line in which GRA2 was disrupted by replacement with the marker Ble (selecting for phleomycin resistance). The Deltagra2 mutant invaded and grew normally in both fibroblasts and macrophages in vitro; however, it was less virulent during acute infection in mice. The survival rate of mice inoculated with Deltagra2 was significantly higher; some infected mice survived the acute infection, whereas all mice infected with the wild-type strain RH succumbed to early death. Chronic infection by Deltagra2 was detected by positive serology, immunohistochemical detection of parasites and cysts in the brain, and reisolation of parasites by bioassay at 6 weeks postinfection. Thus, absence of GRA2 partially attenuates the virulence of T. gondii during the acute phase of infection and allows for establishment of chronic infection by the otherwise highly virulent RH strain. These results establish that GRA2 plays an important role during in vivo infection and provide a potential model for examining acute pathogenesis by T. gondii.

Genotypic analysis of Toxoplasma gondii isolates from pigs

To determine the prevalence of the 3 primary clonal lineages of Toxoplasma gondii (strain types I, II, and III) in a potential food source of infection for humans, we analyzed 43 isolates of T. gondii that had been collected from pigs at an abattoir in Iowa. Parasites were harvested as in vitro-grown tachyzoites, and their genotypes were determined at the SAG1 and SAG2 loci. On the basis of the allele identified at the SAG2 locus, isolates were grouped into 1 of the 3 primary lineages. Type II strains were by far the most prevalent, accounting for 83.7% of the isolates. The type III genotype was identified in only 16.3% of the isolates. These prevalences differ significantly from a previous sampling of isolates from animals but are similar to the frequencies with which they occur in human disease cases. Similar to the previously characterized strain P89, strains P62 and P105 appeared to have recombinant genotypes. The type I genotype was not identified in the isolates from pigs although these strains have previously been shown to account for approximately 10-25% of toxoplasmosis cases in humans.

Development of molecular genetics for Neospora caninum: A complementary system to Toxoplasma gondii

The development of molecular genetics has greatly enhanced the study of Toxoplasma gondii, and investigations into the biology and pathology associated with neosporosis will be similarly benefited by the development of molecular tools for Neospora caninum. We have demonstrated the feasibility of using the existing DNA vectors developed for T. gondii to transfect and transform the Nc-1 strain of Neospora. We have also shown that T. gondii proteins are faithfully expressed and targeted in N. caninum, indicating the suitability of using Neospora as a heterologous expression system for studying T. gondii. These studies provide the basis for initiating molecular genetic studies on N. caninum and will allow for a number of molecular comparisons of these two closely related, though phenotypically distinct, parasites. Here we describe the methods and reagents used to perform genetic manipulations of N. caninum, and we present some of the principles and potential utilities of these molecular studies, including the use of N. caninum as a heterologous system for the study of T. gondii proteins.