Sarcocystis cruzi (Hasselmann, 1923) Wenyon, 1926: redescription, molecular characterization and deposition of life cycle stages specimens in the Smithsonian Museum.

Currently, 7 named Sarcocystis species infect cattle: Sarcocystis hirsuta, S. cruzi, S. hominis, S. bovifelis, S. heydorni, S. bovini and S. rommeli; other, unnamed species also infect cattle. Of these parasites of cattle, a complete life cycle description is known only for S. cruzi, the most pathogenic species in cattle. The life cycle of S. cruzi was completed experimentally in 1982, before related parasite species were structurally characterized, and before the advent of molecular diagnostics; to our knowledge, no archived frozen tissues from the cattle employed in the original descriptions remain for DNA characterization. Here, we isolated DNA from a paraffin-embedded kidney of a calf experimentally infected with S. cruzi in 1980; we then sequenced portions of 18S rRNA, 28S rRNA, COX1 and Acetyl CoA genes and verified that each shares 99­100% similarity to other available isolates attributed to S. cruzi from naturally infected cattle. We also reevaluated histological sections of tissues of calves experimentally infected with S. cruzi in the original description, exploiting improvements in photographic technology to render clearer morphological detail. Finally, we reviewed all available studies of the life cycle of S. cruzi, noting that S. cruzi was transmitted between bison (Bison bison) and cattle (Bos taurus) and that the strain of parasite derived from bison appeared more pathogenic than the cattle strain. Based on these newfound molecular, morphological and physiological data, we thereby redescribed S. cruzi and deposited reference material in the Smithsonian Museum for posterity.

Global selective sweep of a highly inbred genome of the cattle parasite Neospora caninum.

Neospora caninum, a cyst-forming apicomplexan parasite, is a leading cause of neuromuscular diseases in dogs as well as fetal abortion in cattle worldwide. The importance of the domestic and sylvatic life cycles of Neospora, and the role of vertical transmission in the expansion and transmission of infection in cattle, is not sufficiently understood. To elucidate the population genomics of Neospora, we genotyped 50 isolates collected worldwide from a wide range of hosts using 19 linked and unlinked genetic markers. Phylogenetic analysis and genetic distance indices resolved a single genotype of N. caninum Whole-genome sequencing of 7 isolates from 2 different continents identified high linkage disequilibrium, significant structural variation, but only limited polymorphism genome-wide, with only 5,766 biallelic single nucleotide polymorphisms (SNPs) total. Greater than half of these SNPs (∼3,000) clustered into 6 distinct haploblocks and each block possessed limited allelic diversity (with only 4 to 6 haplotypes resolved at each cluster). Importantly, the alleles at each haploblock had independently segregated across the strains sequenced, supporting a unisexual expansion model that is mosaic at 6 genomic blocks. Integrating seroprevalence data from African cattle, our data support a global selective sweep of a highly inbred livestock pathogen that originated within European dairy stock and expanded transcontinentally via unisexual mating and vertical transmission very recently, likely the result of human activities, including recurrent migration, domestication, and breed development of bovid and canid hosts within similar proximities.

CD8α(+) dendritic cells are the critical source of interleukin-12 that controls acute infection by Toxoplasma gondii tachyzoites.

CD8α(+) dendritic cells (DCs) are important in vivo for cross-presentation of antigens derived from intracellular pathogens and tumors. Additionally, secretion of interleukin-12 (IL-12) by CD8α(+) DCs suggests a role for these cells in response to Toxoplasma gondii antigens, although it remains unclear whether these cells are required for protection against T. gondii infection. Toward this goal, we examined T. gondii infection of Batf3(-/-) mice, which selectively lack only lymphoid-resident CD8α(+) DCs and related peripheral CD103(+) DCs. Batf3(-/-) mice were extremely susceptible to T. gondii infection, with decreased production of IL-12 and interferon-γ. IL-12 administration restored resistance in Batf3(-/-) mice, and mice in which IL-12 production was ablated only from CD8α(+) DCs failed to control infection. These results reveal that the function of CD8α(+) DCs extends beyond a role in cross-presentation and includes a critical role for activation of innate immunity through IL-12 production during T. gondii infection.

PopNet: A Markov Clustering Approach to Study Population Genetic Structure.

With the advent of low cost, high-throughput genome sequencing technology, population genomic data sets are being generated for hundreds of species of pathogenic, industrial, and agricultural importance. The challenge is how best to analyze and visually display these complex data sets to yield intuitive representations capable of capturing complex evolutionary relationships. Here we present PopNet, a novel computational method that identifies regions of shared ancestry in the chromosomes of related strains through clustering patterns of genetic variation. These relationships are subsequently visualized within a network by a novel implementation of chromosome painting. We apply PopNet to three diverse populations that feature differential rates of recombination and demonstrate its ability to capture evolutionary relationships as well as associate traits to specific loci. Compared with existing tools, PopNet provides substantial advances by both removing the need to predefine a single reference genome that can bias interpretation of population structure, as well as its ability to visualize multiple evolutionary relationships, such as recombination events and shared ancestry, across hundreds of strains.

Rhoptry Proteins ROP5 and ROP18 Are Major Murine Virulence Factors in Genetically Divergent South American Strains of Toxoplasma gondii.

Toxoplasma gondii has evolved a number of strategies to evade immune responses in its many hosts. Previous genetic mapping of crosses between clonal type 1, 2, and 3 strains of T. gondii, which are prevalent in Europe and North America, identified two rhoptry proteins, ROP5 and ROP18, that function together to block innate immune mechanisms activated by interferon gamma (IFNg) in murine hosts. However, the contribution of these and other virulence factors in more genetically divergent South American strains is unknown. Here we utilized a cross between the intermediately virulent North American type 2 ME49 strain and the highly virulent South American type 10 VAND strain to map the genetic basis for differences in virulence in the mouse. Quantitative trait locus (QTL) analysis of this new cross identified one peak that spanned the ROP5 locus on chromosome XII. CRISPR-Cas9 mediated deletion of all copies of ROP5 in the VAND strain rendered it avirulent and complementation confirmed that ROP5 is the major virulence factor accounting for differences between type 2 and type 10 strains. To extend these observations to other virulent South American strains representing distinct genetic populations, we knocked out ROP5 in type 8 TgCtBr5 and type 4 TgCtBr18 strains, resulting in complete loss of virulence in both backgrounds. Consistent with this, polymorphisms that show strong signatures of positive selection in ROP5 were shown to correspond to regions known to interface with host immunity factors. Because ROP5 and ROP18 function together to resist innate immune mechanisms, and a significant interaction between them was identified in a two-locus scan, we also assessed the role of ROP18 in the virulence of South American strains. Deletion of ROP18 in South American type 4, 8, and 10 strains resulted in complete attenuation in contrast to a partial loss of virulence seen for ROP18 knockouts in previously described type 1 parasites. These data show that ROP5 and ROP18 are conserved virulence factors in genetically diverse strains from North and South America, suggesting they evolved to resist innate immune defenses in ancestral T. gondii strains, and they have subsequently diversified under positive selection.

Genetic mapping reveals that sinefungin resistance in Toxoplasma gondii is controlled by a putative amino acid transporter locus that can be used as a negative selectable marker.

Quantitative trait locus (QTL) mapping studies have been integral in identifying and understanding virulence mechanisms in the parasite Toxoplasma gondii. In this study, we interrogated a different phenotype by mapping sinefungin (SNF) resistance in the genetic cross between type 2 ME49-FUDR(r) and type 10 VAND-SNF(r). The genetic map of this cross was generated by whole-genome sequencing of the progeny and subsequent identification of single nucleotide polymorphisms (SNPs) inherited from the parents. Based on this high-density genetic map, we were able to pinpoint the sinefungin resistance phenotype to one significant locus on chromosome IX. Within this locus, a single nonsynonymous SNP (nsSNP) resulting in an early stop codon in the TGVAND_290860 gene was identified, occurring only in the sinefungin-resistant progeny. Using CRISPR/CAS9, we were able to confirm that targeted disruption of TGVAND_290860 renders parasites sinefungin resistant. Because disruption of the SNR1 gene confers resistance, we also show that it can be used as a negative selectable marker to insert either a positive drug selection cassette or a heterologous reporter. These data demonstrate the power of combining classical genetic mapping, whole-genome sequencing, and CRISPR-mediated gene disruption for combined forward and reverse genetic strategies in T. gondii.

REDHORSE-REcombination and Double crossover detection in Haploid Organisms using next-geneRation SEquencing data.

BACKGROUND: Next-generation sequencing technology provides a means to study genetic exchange at a higher resolution than was possible using earlier technologies. However, this improvement presents challenges as the alignments of next generation sequence data to a reference genome cannot be directly used as input to existing detection algorithms, which instead typically use multiple sequence alignments as input. We therefore designed a software suite called REDHORSE that uses genomic alignments, extracts genetic markers, and generates multiple sequence alignments that can be used as input to existing recombination detection algorithms. In addition, REDHORSE implements a custom recombination detection algorithm that makes use of sequence information and genomic positions to accurately detect crossovers. REDHORSE is a portable and platform independent suite that provides efficient analysis of genetic crosses based on Next-generation sequencing data. RESULTS: We demonstrated the utility of REDHORSE using simulated data and real Next-generation sequencing data. The simulated dataset mimicked recombination between two known haploid parental strains and allowed comparison of detected break points against known true break points to assess performance of recombination detection algorithms. A newly generated NGS dataset from a genetic cross of Toxoplasma gondii allowed us to demonstrate our pipeline. REDHORSE successfully extracted the relevant genetic markers and was able to transform the read alignments from NGS to the genome to generate multiple sequence alignments. Recombination detection algorithm in REDHORSE was able to detect conventional crossovers and double crossovers typically associated with gene conversions whilst filtering out artifacts that might have been introduced during sequencing or alignment. REDHORSE outperformed other commonly used recombination detection algorithms in finding conventional crossovers. In addition, REDHORSE was the only algorithm that was able to detect double crossovers. CONCLUSION: REDHORSE is an efficient analytical pipeline that serves as a bridge between genomic alignments and existing recombination detection algorithms. Moreover, REDHORSE is equipped with a recombination detection algorithm specifically designed for Next-generation sequencing data. REDHORSE is portable, platform independent Java based utility that provides efficient analysis of genetic crosses based on Next-generation sequencing data. REDHORSE is available at http://redhorse.sourceforge.net/ .

Sarcocystis caninum and Sarcocystis svanai n. spp. (Apicomplexa: Sarcocystidae) Associated with Severe Myositis and Hepatitis in the Domestic Dog (Canis familiaris).

There are several reports of Sarcocystis sarcocysts in muscles of dogs, but these species have not been named. Additionally, there are two reports of Sarcocystis neurona in dogs. Here, we propose two new names, Sarcocystis caninum, and Sarcocystis svanai for sarcocysts associated with clinical muscular sarcocystosis in four domestic dogs (Canis familiaris), one each from Montana and Colorado in the USA, and two from British Columbia, Canada. Only the sarcocyst stage was identified. Most of the sarcocysts identified were S. caninum. Sarcocysts were studied using light microscopy, transmission electron microscopy (TEM), and polymerase chain reaction. Based on collective results two new species, S. caninum and S. svanai were designated. Sarcocystis caninum and S. svanai were structurally distinct. Sarcocystis caninum sarcocysts were up to 1.2 mm long and up to 75 µm wide. By light microscopy, the sarcocyst wall was relatively thin and smooth. By TEM, the sarcocyst wall was "type 9", 1-2 µm thick, and contained villar protrusions that lacked microtubules. Bradyzoites in sections were 7-9 µm long. Sarcocysts of S. svanai were few and were identified by TEM. Sarcocystis svanai sarcocysts were "type 1", thin walled (< 0.5 µm), and the wall lacked villar protrusions but had tiny blebs that did not invaginate. DNA was extracted either from infected frozen muscle biopsies or formalin-fixed paraffin-embedded sections. Dogs were either singly infected with S. caninum or multiply co-infected with S. caninum and S. svanai (the result of a mixed infection) based on multilocus DNA sequencing and morphology. BLASTn analysis established that the sarcocysts identified in these dogs were similar to, but not identical to Sarcocystis canis or Sarcocystis arctosi, parasites found to infect polar bears (Ursus maritimus) or brown bears (Ursus arctosi), respectively. However, the S. caninum sequence showed 100% identify over the 18S rRNA region sequenced to that of S. arctica, a parasite known to infect Arctic foxes (Vulpes lagopus).

Globally diverse Toxoplasma gondii isolates comprise six major clades originating from a small number of distinct ancestral lineages.

Marked phenotypic variation characterizes isolates of Toxoplasma gondii, a ubiquitous zoonotic parasite that serves as an important experimental model for studying apicomplexan parasites. Progress in identifying the heritable basis for clinically and epidemiologically significant differences requires a robust system for describing and interpreting evolutionary subdivisions in this prevalent pathogen. To develop such a system, we have examined more than 950 isolates collected from around the world and genotyped them using three independent sets of polymorphic DNA markers, sampling 30 loci distributed across all nuclear chromosomes as well as the plastid genome. Our studies reveal a biphasic pattern consisting of regions in the Northern Hemisphere where a few, highly clonal and abundant lineages predominate; elsewhere, and especially in portions of South America are characterized by a diverse assemblage of less common genotypes that show greater evidence of recombination. Clustering methods were used to organize the marked genetic diversity of 138 unique genotypes into 15 haplogroups that collectively define six major clades. Analysis of gene flow indicates that a small number of ancestral lineages gave rise to the existing diversity through a process of limited admixture. Identification of reference strains for these major groups should facilitate future studies on comparative genomics and identification of genes that control important biological phenotypes including pathogenesis and transmission.

NextGen sequencing reveals short double crossovers contribute disproportionately to genetic diversity in Toxoplasma gondii.

BACKGROUND: Toxoplasma gondii is a widespread protozoan parasite of animals that causes zoonotic disease in humans. Three clonal variants predominate in North America and Europe, while South American strains are genetically diverse, and undergo more frequent recombination. All three northern clonal variants share a monomorphic version of chromosome Ia (ChrIa), which is also found in unrelated, but successful southern lineages. Although this pattern could reflect a selective advantage, it might also arise from non-Mendelian segregation during meiosis. To understand the inheritance of ChrIa, we performed a genetic cross between the northern clonal type 2 ME49 strain and a divergent southern type 10 strain called VAND, which harbors a divergent ChrIa. RESULTS: NextGen sequencing of haploid F1 progeny was used to generate a genetic map revealing a low level of conventional recombination, with an unexpectedly high frequency of short, double crossovers. Notably, both the monomorphic and divergent versions of ChrIa were isolated with equal frequency. As well, ChrIa showed no evidence of being a sex chromosome, of harboring an inversion, or distorting patterns of segregation. Although VAND was unable to self fertilize in the cat, it underwent successful out-crossing with ME49 and hybrid survival was strongly associated with inheritance of ChrIII from ME49 and ChrIb from VAND. CONCLUSIONS: Our findings suggest that the successful spread of the monomorphic ChrIa in the wild has not been driven by meiotic drive or related processes, but rather is due to a fitness advantage. As well, the high frequency of short double crossovers is expected to greatly increase genetic diversity among progeny from genetic crosses, thereby providing an unexpected and likely important source of diversity.

Virulence differences in Toxoplasma mediated by amplification of a family of polymorphic pseudokinases.

The population structure of Toxoplasma gondii includes three highly prevalent clonal lineages referred to as types I, II, and III, which differ greatly in virulence in the mouse model. Previous studies have implicated a family of serine/threonine protein kinases found in rhoptries (ROPs) as important in mediating virulence differences between strain types. Here, we explored the genetic basis of differences in virulence between the highly virulent type I lineage and moderately virulent type II based on successful genetic cross between these lineages. Genome-wide association revealed that a single quantitative trait locus controls the dramatic difference in lethality between these strain types. Neither ROP16 nor ROP18, previously implicated in virulence of T. gondii, was found to contribute to differences between types I and II. Instead, the major virulence locus contained a tandem cluster of polymorphic alleles of ROP5, which showed similar protein expression between strains. ROP5 contains a conserved serine/threonine protein kinase domain that includes only part of the catalytic triad, and hence, all members are considered to be pseudokinases. Genetic disruption of the entire ROP5 locus in the type I lineage led to complete attenuation of acute virulence, and complementation with ROP5 restored lethality to WT levels. These findings reveal that a locus of polymorphic pseudokinases plays an important role in pathogenesis of toxoplasmosis in the mouse model.

"Evaluating the Benefits and Limits of Multiple Displacement Amplification with Whole-Genome Oxford Nanopore Sequencing".

Multiple Displacement Amplification (MDA) outperforms conventional PCR in long fragment and whole genome amplification which makes it attractive to couple with long-read sequencing of samples with limited quantities of DNA to obtain improved genome assemblies. Here, we explore the efficacy and limits of MDA for genome sequence assembly using Oxford Nanopore Technologies (ONT) rapid library preparations and minION sequencing. We successfully generated almost complete genome sequences for all organisms examined, including Cryptosporidium meleagridis, Staphylococcus aureus, Enterococcus faecium, and Escherichia coli, with the ability to generate high-quality data from samples starting with only 0.025 ng of total DNA. Controlled sheared DNA samples exhibited a distinct pattern of size-increase after MDA, which may be associated with the amplification of long, low-abundance fragments present in the assay, as well as generating concatemeric sequences during amplification. To address concatemers, we developed a computational pipeline (CADECT: Concatemer Detection Tool) to identify and remove putative concatemeric sequences. This study highlights the efficacy of MDA in generating high-quality genome assemblies from limited amounts of input DNA. Also, the CADECT pipeline effectively mitigated the impact of concatemeric sequences, enabling the assembly of contiguous sequences even in cases where the input genomic DNA was degraded. These results have significant implications for the study of organisms that are challenging to culture in vitro, such as Cryptosporidium, and for expediting critical results in clinical settings with limited quantities of available genomic DNA.

Comparative genomics of Giardia duodenalis sub-assemblage AI beaver (Be-2) and human (WB-C6) strains show remarkable homozygosity, sequence similarity, and conservation of VSP genes.

Giardia duodenalis, a major cause of waterborne infection, infects a wide range of mammalian hosts and is subdivided into eight genetically well-defined assemblages named A through H. However, fragmented genomes and a lack of comparative analysis within and between the assemblages render unclear the molecular mechanisms controlling host specificity and differential disease outcomes. To address this, we generated a near-complete de novo genome of AI assemblage using the Oxford Nanopore platform by sequencing the Be-2 genome. We generated 148,144 long-reads with quality scores of > 7. The final genome assembly consists of only nine contigs with an N50 of 3,045,186 bp. This assembly agrees closely with the assembly of another strain in the AI assemblage (WB-C6). However, a critical difference is that a region previously placed in the five-prime region of Chr5 belongs to Chr4 of Be-2. We find a high degree of conservation in the ploidy, homozygosity, and the presence of cysteine-rich variant-specific surface proteins (VSPs) within the AI assemblage. Our assembly provides a nearly complete genome of a member of the AI assemblage of G. duodenalis, aiding population genomic studies capable of elucidating Giardia transmission, host range, and pathogenicity.

Evolutionarily divergent, unstable filamentous actin is essential for gliding motility in apicomplexan parasites.

Apicomplexan parasites rely on a novel form of actin-based motility called gliding, which depends on parasite actin polymerization, to migrate through their hosts and invade cells. However, parasite actins are divergent both in sequence and function and only form short, unstable filaments in contrast to the stability of conventional actin filaments. The molecular basis for parasite actin filament instability and its relationship to gliding motility remain unresolved. We demonstrate that recombinant Toxoplasma (TgACTI) and Plasmodium (PfACTI and PfACTII) actins polymerized into very short filaments in vitro but were induced to form long, stable filaments by addition of equimolar levels of phalloidin. Parasite actins contain a conserved phalloidin-binding site as determined by molecular modeling and computational docking, yet vary in several residues that are predicted to impact filament stability. In particular, two residues were identified that form intermolecular contacts between different protomers in conventional actin filaments and these residues showed non-conservative differences in apicomplexan parasites. Substitution of divergent residues found in TgACTI with those from mammalian actin resulted in formation of longer, more stable filaments in vitro. Expression of these stabilized actins in T. gondii increased sensitivity to the actin-stabilizing compound jasplakinolide and disrupted normal gliding motility in the absence of treatment. These results identify the molecular basis for short, dynamic filaments in apicomplexan parasites and demonstrate that inherent instability of parasite actin filaments is a critical adaptation for gliding motility.

Relationship between Eimeria oocyst infectivity for chickens and in vitro excystation of E. acervulina, E. maxima, and E. tenella oocyst during long-term storage.

Vaccination of chickens against avian coccidiosis in chickens often involves storing Eimeria oocysts for months after oocyst propagation and sporulation. The purpose of this study was to determine how long E. acervulina, E. maxima, and E. tenella oocysts remained viable when stored at refrigeration (4°C) or egg room (20°C) temperatures. Separate tubes containing E. acervulina, E. maxima, or E. tenella oocysts were stored at these temperatures and a sample removed every 3 mo for inoculating chickens for evidence of a patent infection. Also, an aliquot of each Eimeria species at each time-temperature combination was subjected to in vitro excystation to quantify the relative number of released sporozoites to intact (nonexcysted) sporocysts. Eimeria tenella appeared to be most susceptible to storage in that no oocyst production was observed at 9 mo at either temperature. Although E. maxima oocysts were viable at 9 mo, no oocyst production was observed at 12 mo storage at these 2 temperatures. Quite unexpected was that E. acervulina was much more stable than E. tenella and E. maxima remaining viable up to and including 27 mo at 4°C and up to and including 12 mo at 20°C. No consistent correlation was observed between in vivo oocyst production and in vitro excystation arising from these 2 respective temperatures (E. acervulina r = 0.58, r = 0.54; E. maxima r = 0.90, r = 0.54; E. tenella r = 0.38, r = 0.90). These data indicate that attention must be paid to time and temperature of Eimeria oocyst storage, and that sporozoite excystation may not be a good indicator of oocyst viability, particularly at later timepoints in incubation.

The same genotype of Sarcocystis neurona responsible for mass mortality in marine mammals induced a clinical outbreak in raccoons (Procyon lotor) 10 years later.

Here, we report the first known outbreak of clinical protozoal myeloencephalitis in naturally infected raccoons by the parasite Sarcocystis neurona. The North American opossum (Didelphis virginiana) and the South American opossum (Didelphis albiventris) are its known definitive hosts. Several other animal species are its intermediate or aberrant hosts. The raccoon (Procyon lotor) is considered the most important intermediate host for S. neurona in the USA. More than 50% of raccoons in the USA have sarcocysts in their muscles, however clinical sarcocystosis in raccoons is rare. In 2014, 38 free-living raccoons were found dead or moribund on the grounds of the Saint Louis Zoo, Missouri, USA. Moribund individuals were weak, lethargic, and mildly ataxic; several with oculo-nasal discharge. Seven raccoons were found dead and 31 were humanely euthanized. Postmortem examinations were conducted on nine raccoons. Neural lesions compatible with acute sarcocystosis were detected in eight raccoons. The predominant lesions were meningoencephalitis and perivascular mononuclear cells. Histologic evidence for the Canine Distemper Virus was found in one raccoon. Schizonts and merozoites were present in the encephalitic lesions of four raccoons. Mature sarcocysts were present within myocytes of five raccoons. In six raccoons, S. neurona schizonts and merozoites were confirmed by immunohistochemical staining with S. neurona-specific polyclonal antibodies. Viable S. neurona was isolated from the brains of two raccoons by bioassay in interferon gamma gene knockout mice and in cell cultures seeded directly with raccoon brain homogenate. Molecular characterization was based on raccoon no. 68. Molecular characterization based on multi-locus typing at five surface antigens (SnSAG1-5-6, SnSAG3 and SnSAG4) and the ITS-1 marker within the ssrRNA locus, using DNA isolated from bradyzoites released from sarcocysts in a naturally infected raccoon (no. 68), confirmed the presence of S. neurona antigen type I, the same genotype that caused a mass mortality event in which 40 southern sea otters stranded dead or dying within a 3 week period in April 2004 with S. neurona-associated disease. An expanded set of genotyping markers was next applied. This study reports the following new genotyping markers at 18S rRNA, 28S rRNA, COX1, ITS-1, RON1, RON2, GAPDH1, ROP20, SAG2, SnSRS21 and TUBA1 markers. The identity of Sarcocystis spp. infecting raccoons is discussed.

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.

Selection at a single locus leads to widespread expansion of Toxoplasma gondii lineages that are virulent in mice.

Pathogenicity differences among laboratory isolates of the dominant clonal North American and European lineages of Toxoplasma gondii are largely controlled by polymorphisms and expression differences in rhoptry secretory proteins (ROPs). However, the extent to which such differences control virulence in natural isolates of T. gondii, including those from more diverse genetic backgrounds, is uncertain. We elucidated the evolutionary history and functional consequences of diversification in the serine/threonine kinase ROP18, a major virulence determinant in the mouse model. We characterized the extent of sequence polymorphism and the evolutionary forces acting on ROP18 and several antigen-encoding genes within a large collection of natural isolates, comparing them to housekeeping genes and introns. Surprisingly, despite substantial genetic diversity between lineages, we identified just three principal alleles of ROP18, which had very ancient ancestry compared to other sampled loci. Expression and allelic differences between these three alleles of ROP18 accounted for much of the variation in acute mouse virulence among natural isolates. While the avirulent type III allele was the most ancient, intermediate virulent (type II) and highly virulent (type I) lineages predominated and showed evidence of strong selective pressure. Out-group comparison indicated that historical loss of an upstream regulatory element increased ROP18 expression, exposing it to newfound diversifying selection, resulting in greatly enhanced virulence in the mouse model and expansion of new lineages. Population sweeps are evident in many genomes, yet their causes and evolutionary histories are rarely known. Our results establish that up-regulation of expression and selection at ROP18 in T. gondii has resulted in three distinct alleles with widely different levels of acute virulence in the mouse model. Preservation of all three alleles in the wild indicates they are likely adaptations for different niches. Our findings demonstrate that sweeping changes in population structure can result from alterations in a single gene.

Life Cycle and Transmission of Cyclospora cayetanensis: Knowns and Unknowns.

Although infections with Cyclospora cayetanensis are prevalent worldwide, many aspects of this parasite's life cycle and transmission remain unknown. Humans are the only known hosts of this parasite. Existing information on its endogenous development has been derived from histological examination of only a few biopsy specimens. Its asexual and sexual stages occur in biliary-intestinal epithelium. In histological sections, its stages are less than 10 µm, making definitive identification difficult. Asexual (schizonts) and sexual (gamonts) are located in epithelial cells. Male microgamonts have two flagella; female macrogametes contain wall-forming bodies. Oocysts are excreted in feces unsporulated. Sporulation occurs in the environment, but there are many unanswered questions concerning dissemination and survival of C. cayetanensis oocysts. Biologically and phylogenetically, C. cayetanensis closely resembles Eimeria spp. that parastize chickens; among them, E. acervulina most closely resembles C. cayetanensis in size. Here, we review known and unknown aspects of its life cycle and transmission and discuss the appropriateness of surrogates best capable of hastening progress in understanding its biology and developing mitigating strategies.

Hastening Progress in Cyclospora Requires Studying Eimeria Surrogates.

Cyclospora cayetanensis is an enigmatic human parasite that sickens thousands of people worldwide. The scarcity of research material and lack of any animal model or cell culture system slows research, denying the produce industry, epidemiologists, and regulatory agencies of tools that might aid diagnosis, risk assessment, and risk abatement. Fortunately, related species offer a strong foundation when used as surrogates to study parasites of this type. Species of Eimeria lend themselves especially well as surrogates for C. cayetanensis. Those Eimeria that infect poultry can be produced in abundance, share many biological features with Cyclospora, pose no risk to the health of researchers, and can be studied in their natural hosts. Here, we overview the actual and potential uses of such surrogates to advance understanding of C. cayetanensis biology, diagnostics, control, and genomics, focusing on opportunities to improve prevention, surveillance, risk assessment, and risk reduction. Studying Eimeria surrogates accelerates progress, closing important research gaps and refining promising tools for producers and food safety regulators to monitor and ameliorate the food safety risks imposed by this emerging, enigmatic parasite.