Chromosomal scale assembly reveals localized structural variants in avian caecal coccidian parasite Eimeria tenella.

Eimeria tenella is a major cause of caecal coccidiosis in commercial poultry chickens worldwide. Here, we report chromosomal scale assembly of Eimeria tenella strain APU2, a strain isolated from commercial broiler chickens in the U.S. We obtained 100× sequencing Oxford Nanopore Technology (ONT) and more than 800× Coverage of Illumina Next-Seq. We created the assembly using the hybrid approach implemented in MaSuRCA, achieving a contiguous 51.34 Mb chromosomal-scale scaffolding enabling identification of structural variations. The AUGUSTUS pipeline predicted 8060 genes, and BUSCO deemed the genomes 99% complete; 6278 (78%) genes were annotated with Pfam domains, and 1395 genes were assigned GO-terms. Comparing E. tenella strains (APU2, US isolate and Houghton, UK isolate) derived Houghton strain of E. tenella revealed 62,905 high stringency differences, of which 45,322 are single nucleotide polymorphisms (SNPs) (0.088%). The rate of transitions/transversions among the SNPs are 1.63 ts/tv. The strains possess conserved gene order but have profound sequence heterogeneity in a several chromosomal segments (chr 2, 11 and 15). Genic and intergenic variation in defined gene families was evaluated between the two strains to possibly identify sequences under selection. The average genic nucleotide diversity of 2.8 with average 2 kb gene length (0.145%) at genic level. We examined population structure using available E. tenella sequences in NCBI, revealing that the two E. tenella isolates from the U.S. (E. tenella APU2 and Wisconsin, "ERR296879") share a common maternal inheritance with the E. tenella Houghton. Our chromosomal level assembly promotes insight into Eimeria biology and evolution, hastening drug discovery and vaccine development.

Genetic evidence substantiates transmission of Trichinella spiralis from one swine farm to another.

BACKGROUND: Trichinella spiralis ranks seventh in the risk posed by foodborne parasites. It causes most human cases of trichinellosis and is the most frequent cause of Trichinella outbreaks on pig farms and in wild boar, worldwide. Veterinary inspectors seek the source of outbreaks in hopes of limiting the spread. Established molecular tools are inadequate for distinguishing among potential T. spiralis infection sources because genetic variability in these zoonotic pathogens is limited in Europe. Microsatellite markers proved successful in tracing an outbreak of T. britovi, a related parasite harboring much more genetic variation. Here, we successfully employed microsatellite markers to determine the genetic structure of T. spiralis isolates from two pig outbreaks, discovering notable uniformity among parasites within each farm and discovering an epidemiological link between these two outbreaks. METHODS: The individual larvae from five isolates of T. spiralis from two pig farms and from ten wild boars were genotyped using nine microsatellite markers to examine their genetic structure. RESULTS: Notably uniform parasite populations constituted each farm outbreak, and the parasites from the first and second outbreaks resembled each other to a notable degree, indicating an epidemiological link between them. Wild boar harbored more genetically variable larval cohorts, distinguishing them from parasites isolated from domestic pigs. CONCLUSIONS: Microsatellite markers succeeded in distinguishing isolates of the highly homogeneous T. spiralis, aiding efforts to track transmission. Each outbreak was composed of a homogenous group of parasites, suggesting a point source of contamination.

Divergence at mitochondrial and ribosomal loci indicates the split between Asian and European populations of Trichinella spiralis occurred prior to swine domestication.

Available evidence suggests that Trichinella spiralis first originated in Asia and subsequently spread to the rest of the world. Notably limited genetic diversity in European T. spiralis isolates indicates that the parasite went through a dramatic genetic bottleneck at some point in its history. Did this genetic bottleneck result from the transport of a limited number of T. spiralis infected pigs from Asian centers of domestication, or was the parasite resident in Europe far earlier than the domestication of pigs there? In order to explore this hypothesis, we generated complete mitochondrial genomes and ribosomal DNAs from seventeen European T. spiralis isolates, six North American isolates and seven Asian isolates using next generation sequencing. A total of 13,858 base pairs of mitochondrial DNA and 7431 nucleotides of the nuclear ribosomal DNA sequence from each isolate were aligned and subjected to phylogenetic analysis using T. nelsoni as an outgroup. We confirmed that North American and European isolates were tightly clustered within a single "western clade" and all Chinese T. spiralis isolates were placed within a well-supported sister clade. These results indicate that European T. spiralis did not directly descend from extant Chinese parasite populations. Furthermore, the amount of nucleotide divergence between the two clades suggests that they diverged before pigs were domesticated. Over evolutionary time periods, Chinese and European T. spiralis were likely maintained as separate populations. The data presented here indicates the genetic bottleneck observed in European T. spiralis did not result from a small number of founders introduced with Chinese pigs in the recent past, but derives from an earlier bottleneck in host populations associated with the end of the last glacial maximum.

Assessing the evolutionary persistence of ecological relationships: A review and preview.

Species interactions, such as pollination, parasitism and predation, form the basis of functioning ecosystems. The origins and resilience of such interactions therefore merit attention. However, fossils only occasionally document ancient interactions, and phylogenetic methods are blind to recent interactions. Is there some other way to track shared species experiences? "Comparative demography" examines when pairs of species jointly thrived or declined. By forging links between ecology, epidemiology, and evolutionary biology, this method sheds light on biological adaptation, species resilience, and ecosystem health. Here, we describe how this method works, discuss examples, and suggest future directions in hopes of inspiring interest, imitators, and critics.

Hiding in plain sight: discovery and phylogeography of a cryptic species of Trichinella (Nematoda: Trichinellidae) in wolverine (Gulo gulo).

Understanding parasite diversity and distribution is essential in managing the potential impact of parasitic diseases in animals and people. Imperfect diagnostic methods, however, may conceal cryptic species. Here, we report the discovery and phylogeography of a previously unrecognized species of Trichinella in wolverine (Gulo gulo) from northwestern Canada that was indistinguishable from T. nativa using the standard multiplex PCR assay based on the expansion segment 5 (ESV) of ribosomal DNA. The novel genotype, designated as T13, was discovered when sequencing the mitochondrial genome. Phylogenetic analyses of the mitochondrial genome and of 15 concatenated single copy orthologs of nuclear DNA indicated a common ancestor for the encapsulated clade is shared by a subclade containing Trichinella spiralis and Trichinella nelsoni, and a subclade containing T13 and remaining taxa: T12 + (T2 + T6) + [(T5 + T9) + (T3 + T8)]. Of 95 individual hosts from 12 species of mammalian carnivores from northwestern Canada from which larvae were identified as T. nativa on multiplex PCR, only wolverines were infected with T13 (14 of 42 individuals). These infections were single or mixed with T. nativa and/or T6. Visual examination and motility testing confirmed that T13 is encapsulated and likely freeze-tolerant. We developed a new Polymerase Chain Reaction-Restriction Fragment Length Polymorphism which unequivocally distinguishes between T13 and T. nativa. We propose Trichinella chanchalensis n. sp. for T13, based on significant genetic divergence from other species of Trichinella and broad-based sampling of the Trichinella genome. Exploration of Alaskan and Siberian isolates may contribute to further resolution of a phylogeographically complex history for species of Trichinella across Beringia, including Trichinella chanchalensis n. sp. (T13).

Comparative demography elucidates the longevity of parasitic and symbiotic relationships.

Parasitic and symbiotic relationships govern vast nutrient and energy flows, yet controversy surrounds their longevity. Enduring relationships may engender parallel phylogenies among hosts and parasites, but so may ephemeral relationships when parasites colonize related hosts. An understanding of whether symbiont and host populations have grown and contracted in concert would be useful when considering the temporal durability of these relationships. Here, we devised methods to compare demographic histories derived from genomic data. We compared the historical growth of the agent of severe human malaria, Plasmodium falciparum, and its mosquito vector, Anopheles gambiae, to human and primate histories, thereby discerning long-term parallels and anthropogenic population explosions. The growth history of Trichinella spiralis, a zoonotic parasite disseminated by swine, proved regionally specific, paralleling distinctive growth histories for wild boar in Asia and Europe. Parallel histories were inferred for an anemone and its algal symbiont (Exaiptasia pallida and Symbiodinium minutum). Concerted growth in potatoes and the agent of potato blight (Solanum tuberosum and Phytophthora infestans) did not commence until the age of potato domestication. Through these examples, we illustrate the utility of comparative historical demography as a new exploratory tool by which to interrogate the origins and durability of myriad ecological relationships. To facilitate future use of this approach, we introduce a tool called C-PSMC to align and evaluate the similarity of demographic history curves.

Histopathological, morphological, and molecular characterization of Sarcocystis species in elk (Cervus elaphus) from Pennsylvania, USA.

Sarcocystis sarcocysts are common in many species of domestic and wild animals. Here, we report sarcocystosis in muscles from 91 free range elk (Cervus elaphus) from Pennsylvania, USA, tested by histopathology, transmission electron microscopy (TEM), and DNA sequencing. Sarcocysts were detected in hematoxylin and eosin (HE)-stained sections from 83 of 91 (91.2%) elk, including 83/91 (91.2%) tongues and 15/17 (88.2%) hearts. With respect to age, sarcocysts were found in 0/5 calves, 8/9 (88.8%) yearlings, and 75/77 (97.4%) adults. Sarcocysts were identified in 62/69 (89.4%) females and 21/22 (91.2%) males. Associated lesions were mild and consisted of inflammatory foci around degenerate sarcocysts. There were two morphologically distinct sarcocysts based on wall thickness, thin (< 0.5 µm) and thick-walled (> 4.0 µm). Thin-walled sarcocysts had a TEM "type 2" and villar protrusions (vps), identical to Sarcocystis wapiti previously described from elk in western USA. This species was present both in tongue and heart samples and was detected in all infected elk. Thick-walled sarcocysts consisted of three morphologic variants, referred to herein as subkinds A, B, C. Subkind A sarcocysts were rare; only four sarcocysts were found in three elk. Histologically, they had a 5-8-µm thick wall with tufted vp. By TEM, the sarcocyst wall was "type 12" and appeared similar to Sarcocystis sybillensis, previously described from elk in USA. Subkind B, Sarcocystis sp.1 sarcocysts were also rare, found in only 1 elk. These sarcocysts had 6.7-7.3-µm-thick wall with TEM "type 15b" vp. Subkind C Sarcocystis sp.2 sarcocysts were more common (22/91). Morphologically, the sarcocyst wall was 6.1-6.8 µm thick and contained "type 10b" vp. Comparisons of ribosomal DNA loci with published sequences indicated all sarcocysts were similar to what has previously been isolated from cervid hosts across the northern hemisphere. Phylogenetic analysis placed the thin-walled S. wapiti within a strongly supported clade with S. linearis and S. taeniata, while the thick-walled cysts were very closely related to S. truncata, S. elongata, S. silva, and S. tarandi. Further sequencing is needed to produce molecular diagnostics to distinguish among these species. North American elk are hosts to multiple Sarcocystis species with diverse morphology, deriving from two separate evolutionary lineages.

Sarcocystis cymruensis: discovery in Western Hemisphere in the Brown rat (Rattus norvegicus) from Grenada, West Indies: redescription, molecular characterization, and transmission to IFN-γ gene knockout mice via sporocysts from experimentally infected domestic cat (Felis catus).

Rodents are intermediate hosts for many species of Sarcocystis. Little is known of Sarcocystis cymruensis that uses the Brown rat (Rattus norvegicus) as intermediate hosts and the domestic cat (Felis catus) as experimental definitive host. Here, we identified and described Sarcocystis cymruensis in naturally infected R. norvegicus from Grenada, West Indies. Rats (n = 167) were trapped in various locations in two parishes (St. George and St. David). Microscopic, thin (< 1 µm) walled, slender sarcocysts were found in 11 of 156 (7.0%) rats skeletal muscles by squash examination. A laboratory-raised cat fed naturally infected rat tissues excreted sporocysts that were infectious for interferon gamma gene knockout (KO) mice, but not to Swiss Webster outbred albino mice. All inoculated mice remained asymptomatic, and microscopic S. cymruensis-like sarcocysts were found in the muscles of KO mice euthanized on day 70, 116, and 189 post inoculation (p.i.). Sarcocysts from infected KO mice were infective for cats at day 116 but not at 70 days p.i. By transmission electron microscopy, the sarcocyst wall was "type 1a." Detailed morphological description of the cyst wall, metrocytes, and bradyzoites is given for the first time. Additionally, molecular data on S. cymruensis are presented also for the first time. Molecular characterization of sarcocysts 18S rDNA and 28S rDNA, ITS-1, and cox1 loci showed the highest similarity with S. rodentifelis and S. muris. In conclusion, the present study described the natural infection of S. cymruensis in Brown rat for the first time in a Caribbean country and provided its molecular characteristics.

Long-read sequencing improves assembly of Trichinella genomes 10-fold, revealing substantial synteny between lineages diverged over 7 million years.

Genome assemblies can form the basis of comparative analyses fostering insight into the evolutionary genetics of a parasite's pathogenicity, host-pathogen interactions, environmental constraints and invasion biology; however, the length and complexity of many parasite genomes has hampered the development of well-resolved assemblies. In order to improve Trichinella genome assemblies, the genome of the sylvatic encapsulated species Trichinella murrelli was sequenced using third-generation, long-read technology and, using syntenic comparisons, scaffolded to a reference genome assembly of Trichinella spiralis, markedly improving both. A high-quality draft assembly for T. murrelli was achieved that totalled 63·2 Mbp, half of which was condensed into 26 contigs each longer than 571 000 bp. When compared with previous assemblies for parasites in the genus, ours required 10-fold fewer contigs, which were five times longer, on average. Better assembly across repetitive regions also enabled resolution of 8 Mbp of previously indeterminate sequence. Furthermore, syntenic comparisons identified widespread scaffold misassemblies in the T. spiralis reference genome. The two new assemblies, organized for the first time into three chromosomal scaffolds, will be valuable resources for future studies linking phenotypic traits within each species to their underlying genetic bases.

Morphological and molecular characterization of Sarcocystis arctica-like sarcocysts from the Arctic fox (Vulpes lagopus) from Alaska, USA.

The muscles of herbivores commonly harbor sarcocysts of parasites belonging to species in the genus Sarcocystis, but such muscle parasites are rare in carnivores. Here, we report Sarcocystis arctica-like sarcocysts in muscles of Arctic foxes (Vulpes lagopus) from Alaska, USA, for the first time. The tongues of 56 foxes were examined for Sarcocystis infection using several methods. Sarcocystis bradyzoites were detected in pepsin digests of 13 (23.2%), and sarcocysts were found in histological sections stained with hematoxylin and eosin (HE) of 9 (16.0%). By light microscopy, sarcocysts were up to 4 mm long and up to 245 µm wide. In HE-stained sections, the sarcocyst wall appeared smooth and up to 1.5 µm thick without visible protrusions. By transmission electron microscopy, the sarcocyst wall had a wavy parasitophorous vacuolar membrane (pvm) folded as pleomorphic villar protrusions (vp), sometimes with anastomoses of villar tips. The vp and the ground substance (gs) layer were smooth and without microtubules. The gs was up to 2.0 µm thick. The total width of the wall including vp and the gs was up to 4.0 µm. The vp were up to 3.0 µm long and most closely resembled "type 9c." All sarcocysts were mature and contained numerous 8.1 × 2.1 µm sized bradyzoites. Molecular characterization (at 18S rDNA, 28S rDNA, ITS-1, and cox1) showed the highest affinity for S. arctica of the Arctic fox (V. lagopus) from Norway. In the present investigation, we provide evidence that sarcocysts are common in tongues of Alaskan Arctic foxes suggesting that these carnivores are serving as intermediate hosts, and we also provide ultrastructure of S. arctica from the Arctic fox for the first time.

Hybridization is limited between two lineages of freeze-resistant Trichinella during coinfection in a mouse model.

Hybridization between two closely related but distinct genetic lineages may lead to homogenization of the two lineages with potentially novel phenotypes, or selective pressure to avoid hybridization if the two lineages are truly distinct. Trichinella nativa and Trichinella T6 are zoonotic nematode parasites which can be distinguished genetically despite occasional hybridization. Here, using an experimental murine model, we attempt to determine whether there are barriers to hybridization when sizeable numbers of each lineage are allowed to coinfect a host. Two mice were independently infected with equal numbers of T. nativa and T6. The offspring of these coinfections were genotyped at two microsatellite loci and one mitochondrial locus capable of distinguishing T. nativa from T6 genotypes. Among larvae in the F1 generation, offspring of every possible mating were encountered. Most larvae (63.6%) derived from T. nativa×T. nativa matings, while 21.1% of offspring were the product of T6×T6 matings, and only 15.3% were hybrid offspring of T. nativa×T6 crosses, differing markedly from null expectations. In this experimental model, T. nativa and Trichinella T6 were able to mate, but ratios of offspring indicated pre- or post-zygotic barriers to hybridization that may include assortative mating, genetic incompatibilities, and/or differences in the fitness of offspring. These barriers would limit gene flow between these two lineages in a natural setting, serving as a barrier to their homogenization and promoting their persistence as distinct and separate entities.

Hybridization between previously isolated ancestors may explain the persistence of exactly two ancient lineages in the genome of the oyster parasite Perkinsus marinus.

Theory predicts that neutral genetic variation accumulates within populations to a level determined by gains through mutation and losses by genetic drift. This balance results in a characteristic distribution of allelic variation with the maximum allelic difference determined by effective population size. Here, we report a striking departure from these expectations in the form of allelic dimorphism, observed at the majority of seven loci examined in Perkinsus marinus, an important oyster parasite that causes Dermo disease. DNA sequences were collected from five loci flanking microsatellite repeats and two loci coding for superoxide dismutase enzymes that may mediate the parasite's interaction with its host. Based on 474 sequences, sampled across 5000 km of the eastern United States coastline, no more than two alleles were observed at each locus (discounting singletons). Depending on the locus, the common allele ranged in overall frequency from 72% to 92%. At each locus the two alleles differed substantially (3.8% sequence difference, on average), and the among-locus variance in divergences was not sufficient to reject a simultaneous origin for all dimorphisms using approximate Bayesian methods. Dimorphic alleles were estimated to have diverged from a common ancestral allele at least 0.9 million years ago. Across these seven loci, only five other alleles were ever observed, always as singletons and differing from the dimorphic alleles by no more than two nucleotides. Free recombination could potentially have shuffled these dimorphisms into as many as 243 multilocus combinations, but the existence of only ten combinations among all samples strongly supports low recombination frequencies and is consistent with the observed absence of intragenic recombination. We consider several demographic and evolutionary hypotheses to explain these patterns. Few can be conclusively rejected with the present data, but we advance a recent hybridization of ancient divergent lineages scenario as the most parsimonious.

Use of leaf litter breakdown and macroinvertebrates to evaluate gradient of recovery in an acid mine impacted stream remediated with an active alkaline doser.

The spatial congruence of chemical and biological recovery along an 18-km acid mine impaired stream was examined to evaluate the efficacy of treatment with an alkaline doser. Two methods were used to evaluate biological recovery: the biological structure of the benthic macroinvertebrate community and several ecosystem processing measures (leaf litter breakdown, microbial respiration rates) along the gradient of improved water chemistry. We found that the doser successfully reduced the acidity and lowered dissolved metals (Al, Fe, and Mn), but downstream improvements were not linear. Water chemistry was more variable, and precipitated metals were elevated in a 3-5-km "mixing zone" immediately downstream of the doser, then stabilized into a "recovery zone" 10-18 km below the doser. Macroinvertebrate communities exhibited a longitudinal pattern of recovery, but it did not exactly match the water chemistry gradient Taxonomic richness (number of families) recovered about 6.5 km downstream of the doser, while total abundance and % EPT taxa recovery were incomplete except at the most downstream site, 18 km away. The functional measures of ecosystem processes (leaf litter breakdown, microbial respiration of conditioned leaves, and shredder biomass) closely matched the measures of community structure and also showed a more modest longitudinal trend of biological recovery than expected based on pH and alkalinity. The measures of microbial respiration had added diagnostic value and indicated that biological recovery downstream of the doser is limited by factors other than habitat and acidity/alkalinity, perhaps episodes of AMD and/or impaired energy/nutrient inputs. A better understanding of the factors that govern spatial and temporal variations in acid mine contaminants, especially episodic events, will improve our ability to predict biological recovery after remediation.

An evolutionary legacy of sex and clonal reproduction in the protistan oyster parasite Perkinsus marinus.

Perkinsus marinus, a protozoan parasite of the eastern oyster Crassostrea virginica, causes Dermo disease which limits fecundity and causes high mortality in host populations. The long-term efficacy of management strategies for suppressing this disease in both aquaculture and restoration settings depends on the potential rate of evolutionary response by P. marinus. Sexual reproduction has never been demonstrated in vitro or in previous population genetic studies. We developed high resolution microsatellite markers and amplified alleles directly from infected oyster genomic DNA. Of 336 infected oysters from four populations between Massachusetts and Florida, 129 (48%) appeared to be infected with a single parasite genotype and were subjected to population genetic analyses assuming diploidy. The great diversity of multilocus genotypes observed is incompatible with strictly clonal reproduction. Substantial heterozygote deficits in three populations suggest that sexual reproduction often involves inbreeding. At the same time, significant multilocus linkage disequilibrium occurred in most sampled populations, and several genotypes were sampled repeatedly in each of two populations, indicating that asexual reproduction also occurs in P. marinus populations. Interestingly, where this parasite has recently expanded its range, lower strain diversity, significant heterozygote excess, and highly heterozygous multilocus genotypes suggests clonal propagation of recent recombinants. Taken together, these data suggest that P. marinus employs multiple reproductive modes, and that over the short term, selection acts upon independent parasite lineages rather than upon individual loci in a cohesive, interbreeding population. Nevertheless, high genotypic diversity is the evolutionary legacy of sex in P. marinus. Anthropogenic movement of infected oysters may increase outcrossing opportunities, potentially facilitating rapid evolution of this parasite.