The genus Leishmania includes approximately 53 species, 20 of which cause human leishmaniais; a significant albeit neglected tropical disease. Leishmaniasis has afflicted humans for millennia, but how ancient is Leishmania and where did it arise? These questions have been hotly debated for decades and several theories have been proposed. One theory suggests Leishmania originated in the Palearctic, and dispersed to the New World via the Bering land bridge. Others propose that Leishmania evolved in the Neotropics. The Multiple Origins theory suggests that separation of certain Old World and New World species occurred due to the opening of the Atlantic Ocean. Some suggest that the ancestor of the dixenous genera Leishmania, Endotrypanum and Porcisia evolved on Gondwana between 90 and 140 million years ago. In the present study a detailed molecular and morphological characterisation was performed on a novel Australian trypanosomatid following its isolation in Australia's tropics from the native black fly, Simulium (Morops) dycei Colbo, 1976. Phylogenetic analyses were conducted and confirmed this parasite as a sibling to Zelonia costaricensis, a close relative of Leishmania previously isolated from a reduviid bug in Costa Rica. Consequently, this parasite was assigned the name Zelonia australiensis sp. nov. Assuming Z. costaricensis and Z. australiensis diverged when Australia and South America became completely separated, their divergence occurred between 36 and 41 million years ago at least. Using this vicariance event as a calibration point for a phylogenetic time tree, the common ancestor of the dixenous genera Leishmania, Endotrypanum and Porcisia appeared in Gondwana approximately 91 million years ago. Ultimately, this study contributes to our understanding of trypanosomatid diversity, and of Leishmania origins by providing support for a Gondwanan origin of dixenous parasitism in the Leishmaniinae.

Hierarchical clustering of pathogen genotypes is widely used to complement epidemiologic investigations of outbreaks. Investigators must dissect trees to obtain genetic partitions that provide epidemiologists with meaningful information. Statistical approaches to tree dissection often require a user-defined parameter to predict the optimal partition number and augmenting this parameter can drastically impact resultant partition memberships. Here, we demonstrate how to optimize a given tree dissection parameter to maximize accuracy irrespective of the tree dissection method used. We hierarchically clustered 1,873 genotypes of the foodborne pathogen Cyclospora spp., including 587 possessing links to historic outbreaks. We dissected the resulting tree using a statistical method requiring users to select the value of a 'stringency parameter' (s), with a recommended value of 95% to 99.5%. We dissected this hierarchical tree across s-values from 94% to 99.5% (at increments of 0.25%), to identify a value that maximized partitioning accuracy, defined as the degree to which genetic partitions conform to known epidemiologic groupings. We show that s-values of 96.5% and 96.75% yield the highest accuracy (> 99.9%) when clustering Cyclospora sp. isolates with known epidemiologic linkages. In practice, the optimized s-value will generate robust genetic partitions comprising isolates likely derived from a common food source, even when the epidemiologic grouping is not known prior to genetic clustering. While the s-value is specific to the tree dissection method used here, the optimization approach described could be applied to any parameter/method used to dissect hierarchical trees.

Cyclospora cayetanensis is an intestinal parasite responsible for the diarrheal illness, cyclosporiasis. Molecular genotyping, using targeted amplicon sequencing, provides a complementary tool for outbreak investigations, especially when epidemiological data are insufficient for linking cases and identifying clusters. The goal of this study was to identify candidate genotyping markers using a novel workflow for detection of segregating single nucleotide polymorphisms (SNPs) in C. cayetanensis genomes. Four whole C. cayetanensis genomes were compared using this workflow and four candidate markers were selected for evaluation of their genotyping utility by PCR and Sanger sequencing. These four markers covered 13 SNPs and resolved parasites from 57 stool specimens, differentiating C. cayetanensis into 19 new unique genotypes.

Previously, it was suggested that haemadipsid leeches represent an important vector of trypanosomes amongst native animals in Australia. Consequently, Chtonobdella bilineata leeches were investigated for the presence of trypanosome species by polymerase chain reaction (PCR), DNA sequencing and in vitro isolation. Phylogenetic analysis ensued to further define the populations present. PCR targeting the 28S rDNA demonstrated that over 95% of C. bilineata contained trypanosomes; diversity profiling by deep amplicon sequencing of 18S rDNA indicated the presence of four different clusters related to the Trypanosoma (Megatrypanum) theileri. Novy–MacNeal–Nicolle slopes with liquid overlay were used to isolate trypanosomes into culture that proved similar in morphology to Trypanosoma cyclops in that they contained a large numbers of acidocalcisomes. Phylogeny of 18S rDNA/GAPDH/ND5 DNA sequences from primary cultures and subclones showed the trypanosomes were monophyletic, with T. cyclops as a sister group. Blood-meal analysis of leeches showed that leeches primarily contained blood from swamp wallaby (Wallabia bicolour), human (Homo sapiens) or horse (Equus sp.). The leech C. bilineata is a host for at least five lineages of Trypanosoma sp. and these are monophyletic with T. cyclops; we propose Trypanosoma cyclops australiensis as a subspecies of T. cyclops based on genetic similarity and biogeography considerations.

Targeted amplicon deep sequencing (TADS) of the 16S rRNA gene is commonly used to explore and characterize bacterial microbiomes. Meanwhile, attempts to apply TADS to the detection and characterization of entire parasitic communities have been hampered since conserved regions of many conserved parasite genes, such as the 18S rRNA gene, are also conserved in their eukaryotic hosts. As a result, targeted amplification of 18S rRNA from clinical samples using universal primers frequently results in competitive priming and preferential amplification of host DNA. Here, we describe a novel method that employs a single pair of universal primers to capture all blood-borne parasites while reducing host 18S rRNA template and enhancing the amplification of parasite 18S rRNA for TADS. This was achieved using restriction enzymes to digest the 18S rRNA gene at cut sites present only in the host sequence prior to PCR amplification.

Neospora caninum is an apicomplexan parasite responsible for neosporosis, a disease causing hind limb paralysis in dogs and abortion in cattle, resulting in substantial economic losses to beef and dairy industries. Marked differences in pathogenicity exist between N. caninum strains suggesting that intrinsic genetic differences exist between them. These differences likely exist in genes expressed during the tachyzoite lifecycle stage which is responsible for the pathogenesis of neosporosis. An improved understanding of these genetic differences is essential to understanding N. caninum virulence, though such knowledge is scarce. Using a variant detection workflow we compared the tachyzoite transcriptomes of two N. caninum strains with different virulence properties: NC-Liverpool (virulent) and NC-Nowra (avirulent). This workflow identified 3130 SNPs and 6123 indels between the strains, and nine markers capturing 30 variants were Sanger sequenced for both strains. Sequencing of these loci was extended to an additional eight strains and subsequent phylogenetic analysis supported a genetic population structure comprised of two major clades with no geographical segregation. Sequence polymorphisms within coding regions of tachyzoite-associated genes were concentrated on chromosomes XI and XII, with 19 distinct tachyzoite-associated SNP hotspot regions identified within coding regions of the N. caninum nuclear genome. The variants were predominantly located in loci associated with protein binding, protein-protein interactions, transcription, and translation. Furthermore, 468 nonsynonymous SNPs identified within protein-coding genes were associated with protein kinase activity, protein binding, protein phosphorylation, and proteolysis. This work may implicate these processes and the specific proteins involved as novel effectors of N. caninum tachyzoite virulence.

Cyclosporiasis is an infection caused by Cyclospora cayetanensis, which is acquired by consumption of contaminated fresh food or water. In the United States, cases of cyclosporiasis are often associated with foodborne outbreaks linked to imported fresh produce or travel to disease-endemic countries. Epidemiologic investigation has been the primary method for linking outbreak cases. A molecular typing marker that can identify genetically related samples would be helpful in tracking outbreaks. We evaluated the mitochondrial junction region as a potential genotyping marker. We tested stool samples from 134 laboratory-confirmed cases in the United States by using PCR and Sanger sequencing. All but 2 samples were successfully typed and divided into 14 sequence types. Typing results were identical among samples within each epidemiologically defined case cluster for 7 of 10 clusters. These findings suggest that this marker can distinguish between distinct case clusters and might be helpful during cyclosporiasis outbreak investigations.

Targeted amplicon deep sequencing (TADS) has enabled characterization of diverse bacterial communities, yet the application of TADS to communities of parasites has been relatively slow to advance. The greatest obstacle to this has been the genetic diversity of parasitic agents, which include helminths, protozoa, arthropods, and some acanthocephalans. Meanwhile, universal amplification of conserved loci from all parasites without amplifying host DNA has proven challenging. Pan-eukaryotic PCRs preferentially amplify the more abundant host DNA, obscuring parasite-derived reads following TADS. Flaherty et al. (2018) described a pan-parasitic TADS method involving amplification of eukaryotic 18S rDNA regions possessing restriction sites only in vertebrates. Using this method, host DNA in total DNA extracts could be selectively digested prior to PCR using restriction enzymes, thereby increasing the number of parasite-derived reads obtained following NGS. This approach showed promise though was only as sensitive as conventional PCR.

Outbreaks of cyclosporiasis, a food-borne illness caused by the coccidian parasite Cyclospora cayetanensis have increased in the USA in recent years, with approximately 2300 laboratory-confirmed cases reported in 2018. Genotyping tools are needed to inform epidemiological investigations, yet genotyping Cyclospora has proven challenging due to its sexual reproductive cycle which produces complex infections characterized by high genetic heterogeneity. We used targeted amplicon deep sequencing and a recently described ensemble-based distance statistic that accommodates heterogeneous (mixed) genotypes and specimens with partial genotyping data, to genotype and cluster 648 C. cayetanensis samples submitted to CDC in 2018. The performance of the ensemble was assessed by comparing ensemble-identified genetic clusters to analogous clusters identified independently based on common food exposures. Using these epidemiologic clusters as a gold standard, the ensemble facilitated genetic clustering with 93.8% sensitivity and 99.7% specificity. Hence, we anticipate that this procedure will greatly complement epidemiologic investigations of cyclosporiasis.

Cyclosporiasis results from an infection of the small intestine by Cyclospora parasites after ingestion of contaminated food or water, often leading to gastrointestinal distress. Recent developments in temporally linking genetically related Cyclospora isolates demonstrated effectiveness in supporting epidemiological investigations. We used 'temporal-genetic clusters' (TGCs) to investigate reported cyclosporiasis cases in the United States during the 2021 peak-period (1 May - 31 August 2021). Our approach split 655 genotyped isolates into 55 genetic clusters and 31 TGCs. We linked two large multi-state epidemiological clusters (Epidemiologic Cluster 1 [n = 136 cases, 54 genotyped] and Epidemiologic Cluster 2 [n = 42 cases, 15 genotyped]) to consumption of lettuce varieties; however, product traceback did not identify a specific product for either cluster due to the lack of detailed product information. To evaluate the utility of TGCs, we performed a retrospective case study comparing investigation outcomes of outbreaks first detected using epidemiological methods with those of the same outbreaks had TGCs been used to first detect them. Our study results indicate that adjustments to routine epidemiological approaches could link additional cases to epidemiological clusters of cyclosporiasis. Overall, we show that CDC's integrated genotyping and epidemiological investigations provide valuable insights into cyclosporiasis outbreaks in the United States.