Epidemiologic Utility of a Framework for Partition Number Selection When Dissecting Hierarchically Clustered Genetic Data Evaluated on the Intestinal Parasite Cyclospora cayetanensis.

Comparing parasite genotypes to inform parasitic disease outbreak investigations involves computation of genetic distances that are typically analyzed by hierarchical clustering to identify related isolates, indicating a common source. A limitation of hierarchical clustering is that hierarchical clusters are not discrete; they are nested. Consequently, small groups of similar isolates exist within larger groups that get progressively larger as relationships become increasingly distant. Investigators must dissect hierarchical trees at a partition number ensuring grouped isolates belong to the same strain; a process typically performed subjectively, introducing bias into resultant groupings. We describe an unbiased, probabilistic framework for partition number selection that ensures partitions comprise isolates that are statistically likely to belong to the same strain. We computed distances and established a normalized distribution of background distances that we used to demarcate a threshold below which the closeness of relationships is unlikely to be random. Distances are hierarchically clustered and the dendrogram dissected at a partition number where most within-partition distances fall below the threshold. We evaluated this framework by partitioning 1,137 clustered Cyclospora cayetanensis genotypes, including 552 isolates epidemiologically linked to various outbreaks. The framework was 91% sensitive and 100% specific in assigning epidemiologically linked isolates to the same partition.

Retrospective evaluation of an integrated molecular-epidemiological approach to cyclosporiasis outbreak investigations - United States, 2021.

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.

Cyclospora cayetanensis comprises at least 3 species that cause human cyclosporiasis.

The apicomplexan parasite Cyclospora cayetanensis causes seasonal foodborne outbreaks of the gastrointestinal illness cyclosporiasis. Prior to the coronavirus disease-2019 pandemic, annually reported cases were increasing in the USA, leading the US Centers for Disease Control and Prevention to develop a genotyping tool to complement cyclosporiasis outbreak investigations. Thousands of US isolates and 1 from China (strain CHN_HEN01) were genotyped by Illumina amplicon sequencing, revealing 2 lineages (A and B). The allelic composition of isolates was examined at each locus. Two nuclear loci (CDS3 and 360i2) distinguished lineages A and B. CDS3 had 2 major alleles: 1 almost exclusive to lineage A and the other to lineage B. Six 360i2 alleles were observed ­ 2 exclusive to lineage A (alleles A1 and A2), 2 to lineage B (B1 and B2) and 1 (B4) was exclusive to CHN_HEN01 which shared allele B3 with lineage B. Examination of heterozygous genotypes revealed that mixtures of A- and B-type 360i2 alleles occurred rarely, suggesting a lack of gene flow between lineages. Phylogenetic analysis of loci from whole-genome shotgun sequences, mitochondrial and apicoplast genomes, revealed that CHN_HEN01 represents a distinct lineage (C). Retrospective examination of epidemiologic data revealed associations between lineage and the geographical distribution of US infections plus strong temporal associations. Given the multiple lines of evidence for speciation within human-infecting Cyclospora, we provide an updated taxonomic description of C. cayetanensis, and describe 2 novel species as aetiological agents of human cyclosporiasis: Cyclospora ashfordi sp. nov. and Cyclospora henanensis sp. nov. (Apicomplexa: Eimeriidae).

Simultaneous targeted amplicon deep sequencing and library preparation for a time and cost-effective universal parasite diagnostic sequencing approach.

We recently described a targeted amplicon deep sequencing (TADS) strategy that utilizes a nested PCR targeting the 18S rDNA gene of blood-borne parasites. The assay facilitates selective digestion of host DNA by targeting enzyme restriction sites present in vertebrates but absent in parasites. This enriching of parasite-derived amplicon drastically reduces the proportion of host-derived reads during sequencing and results in the sensitive detection of several clinically important blood parasites including Plasmodium spp., Babesia spp., kinetoplastids, and filarial nematodes. Despite these promising results, high costs and the laborious nature of metagenomics sequencing are prohibitive to the routine use of this assay in most laboratories. We describe and evaluate a new metagenomic approach that utilizes a set of primers modified from our original assay that incorporates Illumina barcodes and adapters during the PCR steps. This modification makes amplicons immediately compatible with sequencing on the Illumina MiSeq platform, removing the need for a separate library preparation, which is expensive and time-consuming. We compared this modified assay to our previous nested TADS assay in terms of preparation speed, limit of detection (LOD), and cost. Our modifications reduced assay turnaround times from 7 to 5 days. The cost decreased from approximately $40 per sample to $11 per sample. The modified assay displayed comparable performance in the detection and differentiation of human-infecting Plasmodium spp., Babesia spp., kinetoplastids, and filarial nematodes in clinical samples. The LOD of this modified approach was determined for malaria parasites and remained similar to that previously reported for our earlier assay (0.58 Plasmodium falciparum parasites/µL of blood). These modifications markedly reduced costs and turnaround times, making the assay more amenable to routine diagnostic applications.

Genetic characterization and description of Leishmania (Leishmania) ellisi sp. nov.: a new human-infecting species from the USA.

In a 2018 report, an unusual case of cutaneous leishmaniasis was described in a 72-year-old female patient residing in Arizona, United States of America (USA). Preliminary analysis of the 18S rDNA and glyceraldehyde-3-phosphate dehydrogenase genes supported the conclusion that the Leishmania strain (strain 218-L139) isolated from this case was a novel species, though a complete taxonomic description was not provided. Identification of Leishmania at the species level is critical for clinical management and epidemiologic investigations so it is important that novel human-infecting species are characterized taxonomically and assigned a unique scientific name compliant with the ICZN code. Therefore, we sought to provide a complete taxonomic description of Leishmania strain 218-L139. Phylogenetic analysis of several nuclear loci and partial maxicircle genome sequences supported its position within the subgenus Leishmania and further clarified the distinctness of this new species. Morphological characterization of cultured promastigotes and amastigotes from the original case material is also provided. Thus, we conclude that Leishmania (Leishmania) ellisi is a new cause of autochthonous cutaneous leishmaniasis in the USA.

Genotyping Cyclospora cayetanensis From Multiple Outbreak Clusters With An Emphasis on a Cluster Linked to Bagged Salad Mix-United States, 2020.

Cyclosporiasis is a diarrheal illness caused by the foodborne parasite Cyclospora cayetanensis. Annually reported cases have been increasing in the United States prompting development of genotyping tools to aid cluster detection. A recently developed Cyclospora genotyping system based on 8 genetic markers was applied to clinical samples collected during the cyclosporiasis peak period of 2020, facilitating assessment of its epidemiologic utility. While the system performed well and helped inform epidemiologic investigations, inclusion of additional markers to improve cluster detection was supported. Consequently, investigations have commenced to identify additional markers to enhance performance.

Evaluation of various distance computation methods for construction of haplotype-based phylogenies from large MLST datasets.

Multi-locus sequence typing (MLST) is widely used to investigate genetic relationships among eukaryotic taxa, including parasitic pathogens. MLST analysis workflows typically involve construction of alignment-based phylogenetic trees - i.e., where tree structures are computed from nucleotide differences observed in a multiple sequence alignment (MSA). Notably, alignment-based phylogenetic methods require that all isolates/taxa are represented by a single sequence. When multiple loci are sequenced these sequences may be concatenated to produce one tree that includes information from all loci. Alignment-based phylogenetic techniques are robust and widely used yet possess some shortcomings, including how heterozygous sites are handled, intolerance for missing data (i.e., partial genotypes), and differences in the way insertions-deletions (indels) are scored/treated during tree construction. In certain contexts, 'haplotype-based' methods may represent a viable alternative to alignment-based techniques, as they do not possess the aforementioned limitations. This is namely because haplotype-based methods assess genetic similarity based on numbers of shared (i.e., intersecting) haplotypes as opposed to similarities in nucleotide composition observed in an MSA. For haplotype-based comparisons, choosing an appropriate distance statistic is fundamental, and several statistics are available to choose from. However, a comprehensive assessment of various available statistics for their ability to produce a robust haplotype-based phylogenetic reconstruction has not yet been performed. We evaluated seven distance statistics by applying them to extant MLST datasets from the gastrointestinal parasite Cyclospora cayetanensis and two species of pathogenic nematode of the genus Strongyloides. We compare the genetic relationships identified using each statistic to epidemiologic, geographic, and host metadata. We show that Barratt's heuristic definition of genetic distance was the most robust among the statistics evaluated. Consequently, it is proposed that Barratt's heuristic represents a useful approach for use in the context of challenging MLST datasets possessing features (i.e., high heterozygosity, partial genotypes, and indel or repeat-based polymorphisms) that confound or preclude the use of alignment-based methods.

Investigation of US Cyclospora cayetanensis outbreaks in 2019 and evaluation of an improved Cyclospora genotyping system against 2019 cyclosporiasis outbreak clusters.

Cyclosporiasis is an illness characterised by watery diarrhoea caused by the food-borne parasite Cyclospora cayetanensis. The increase in annual US cyclosporiasis cases led public health agencies to develop genotyping tools that aid outbreak investigations. A team at the Centers for Disease Control and Prevention (CDC) developed a system based on deep amplicon sequencing and machine learning, for detecting genetically-related clusters of cyclosporiasis to aid epidemiologic investigations. An evaluation of this system during 2018 supported its robustness, indicating that it possessed sufficient utility to warrant further evaluation. However, the earliest version of CDC's system had some limitations from a bioinformatics standpoint. Namely, reliance on proprietary software, the inability to detect novel haplotypes and absence of a strategy to select an appropriate number of discrete genetic clusters would limit the system's future deployment potential. We recently introduced several improvements that address these limitations and the aim of this study was to reassess the system's performance to ensure that the changes introduced had no observable negative impacts. Comparison of epidemiologically-defined cyclosporiasis clusters from 2019 to analogous genetic clusters detected using CDC's improved system reaffirmed its excellent sensitivity (90%) and specificity (99%), and confirmed its high discriminatory power. This C. cayetanensis genotyping system is robust and with ongoing improvement will form the basis of a US-wide C. cayetanensis genotyping network for clinical specimens.

A new subspecies of Trypanosoma cyclops found in the Australian terrestrial leech Chtonobdella bilineata.

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.

Evaluation of an ensemble-based distance statistic for clustering MLST datasets using epidemiologically defined clusters of cyclosporiasis.

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.

Machine learning-based analyses support the existence of species complexes for Strongyloides fuelleborni and Strongyloides stercoralis.

Human strongyloidiasis is a serious disease mostly attributable to Strongyloides stercoralis and to a lesser extent Strongyloides fuelleborni, a parasite mainly of non-human primates. The role of animals as reservoirs of human-infecting Strongyloides is ill-defined, and whether dogs are a source of human infection is debated. Published multi-locus sequence typing (MLST) studies attempt to elucidate relationships between Strongyloides genotypes, hosts, and distributions, but typically examine relatively few worms, making it difficult to identify population-level trends. Combining MLST data from multiple studies is often impractical because they examine different combinations of loci, eliminating phylogeny as a means of examining these data collectively unless hundreds of specimens are excluded. A recently-described machine learning approach that facilitates clustering of MLST data may offer a solution, even for datasets that include specimens sequenced at different combinations of loci. By clustering various MLST datasets as one using this procedure, we sought to uncover associations among genotype, geography, and hosts that remained elusive when examining datasets individually. Multiple datasets comprising hundreds of S. stercoralis and S. fuelleborni individuals were combined and clustered. Our results suggest that the commonly proposed 'two lineage' population structure of S. stercoralis (where lineage A infects humans and dogs, lineage B only dogs) is an over-simplification. Instead, S. stercoralis seemingly represents a species complex, including two distinct populations over-represented in dogs, and other populations vastly more common in humans. A distinction between African and Asian S. fuelleborni is also supported here, emphasizing the need for further resolving these taxonomic relationships through modern investigations.

Diversity profiling of xenic cultures of Dientamoeba fragilis following systematic antibiotic treatment and prospects for genome sequencing.

The presence of bacterial DNA in Dientamoeba fragilis DNA extracts from culture poses a substantial challenge to sequencing the D. fragilis genome. However, elimination of bacteria from D. fragilis cultures has proven difficult in the past, presumably due to its dependence on some unknown prokaryote/s. This study explored options for removal of bacteria from D. fragilis cultures and for the generation of genome sequence data from D. fragilis. DNA was extracted from human faecal samples and xenic D. fragilis cultures. Extracts were subjected to 16S ribosomal DNA bacterial diversity profiling. Xenic D. fragilis cultures were then subject to antibiotic treatment regimens that systematically removed bacterial species depending on their membrane structure (Gram-positive or Gram-negative) and aerobic requirements. The impact of these treatments on cultures was assessed by 16S amplicon sequencing. Prior to antibiotic treatment, the cultures were dominated by Gram-negative bacteria. Addition of meropenem to cultures eliminated anaerobic Gram-negative bacteria, but it also led to protozoan death after 5 days incubation. The seeding of meropenem resistant Klebsiella pneumoniae strain KPC-2 into cultures before treatment by meropenem prevented death of D. fragilis cells beyond this 5 day period, suggesting that one or more species of Gram-negative bacteria may be an essential nutritional requirement for D. fragilis. Gram-positive cells were completely eliminated using vancomycin without affecting trophozoite growth. Finally, this study shows that genome sequencing of D. fragilis is feasible following bacterial elimination from cultures as the result of the major advances occurring in bioinformatics. We provide evidence on this fact by successfully sequencing the D. fragilis 28S large ribosomal DNA subunit gene using culture-derived DNA.

Mitochondrial Junction Region as Genotyping Marker for Cyclospora cayetanensis.

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.

Genotyping genetically heterogeneous Cyclospora cayetanensis infections to complement epidemiological case linkage.

Sexually reproducing pathogens such as Cyclospora cayetanensis often produce genetically heterogeneous infections where the number of unique sequence types detected at any given locus varies depending on which locus is sequenced. The genotypes assigned to these infections quickly become complex when additional loci are analysed. This genetic heterogeneity confounds the utility of traditional sequence-typing and phylogenetic approaches for aiding epidemiological trace-back, and requires new methods to address this complexity. Here, we describe an ensemble of two similarity-based classification algorithms, including a Bayesian and heuristic component that infer the relatedness of C. cayetanensis infections. The ensemble requires a set of haplotypes as input and assigns arbitrary distances to specimen pairs reflecting their most likely relationships. The approach was applied to data generated from a test cohort of 88 human fecal specimens containing C. cayetanensis, including 30 from patients whose infections were associated with epidemiologically defined outbreak clusters of cyclosporiasis. The ensemble assigned specimens to plausible clusters of genetically related infections despite their complex haplotype composition. These relationships were corroborated by a significant number of epidemiological linkages (P < 0.0001) suggesting the ensemble's utility for aiding epidemiological trace-back investigations of cyclosporiasis.

Dientamoeba fragilis, the Neglected Trichomonad of the Human Bowel.

Dientamoeba fragilis is a protozoan parasite of the human bowel, commonly reported throughout the world in association with gastrointestinal symptoms. Despite its initial discovery over 100 years ago, arguably, we know less about this peculiar organism than any other pathogenic or potentially pathogenic protozoan that infects humans. The details of its life cycle and mode of transmission are not completely known, and its potential as a human pathogen is debated within the scientific community. Recently, several major advances have been made with respect to this organism's life cycle and molecular biology. While many questions remain unanswered, these and other recent advances have given rise to some intriguing new leads, which will pave the way for future research. This review encompasses a large body of knowledge generated on various aspects of D. fragilis over the last century, together with an update on the most recent developments. This includes an update on the latest diagnostic techniques and treatments, the clinical aspects of dientamoebiasis, the development of an animal model, the description of a D. fragilis cyst stage, and the sequencing of the first D. fragilis transcriptome.

Angiostrongylus cantonensis: a review of its distribution, molecular biology and clinical significance as a human pathogen.

Angiostrongylus cantonensis is a metastrongyloid nematode found widely in the Asia-Pacific region, and the aetiological agent of angiostrongyliasis; a disease characterized by eosinophilic meningitis. Rattus rats are definitive hosts of A. cantonensis, while intermediate hosts include terrestrial and aquatic molluscs. Humans are dead-end hosts that usually become infected upon ingestion of infected molluscs. A presumptive diagnosis is often made based on clinical features, a history of mollusc consumption, eosinophilic pleocytosis in cerebral spinal fluid, and advanced imaging such as computed tomography. Serological tests are available for angiostrongyliasis, though many tests are still under development. While there is no treatment consensus, therapy often includes a combination of anthelmintics and corticosteroids. Angiostrongyliasis is relatively rare, but is often associated with morbidity and sometimes mortality. Recent reports suggest the parasites' range is increasing, leading to fatalities in regions previously considered Angiostrongylus-free, and sometimes, delayed diagnosis in newly invaded regions. Increased awareness of angiostrongyliasis would facilitate rapid diagnosis and improved clinical outcomes. This paper summarizes knowledge on the parasites' life cycle, clinical aspects and epidemiology. The molecular biology of Angiostrongylus spp. is also discussed. Attention is paid to the significance of angiostrongyliasis in Australia, given the recent severe cases reported from the Sydney region.

Evaluation of the Increased Genetic Resolution and Utility for Source Tracking of a Recently Developed Method for Genotyping Cyclospora cayetanensis.

Cyclospora cayetanensis is a foodborne parasite that causes cyclosporiasis, an enteric illness in humans. Genotyping methods are used to genetically discriminate between specimens from cyclosporiasis cases and can complement source attribution investigations if the method is sufficiently sensitive for application to food items. A very sensitive targeted amplicon sequencing (TAS) assay for genotyping C. cayetanensis encompassing 52 loci was recently designed. In this study, we analyzed 66 genetically diverse clinical specimens to assess the change in phylogenetic resolution between the TAS assay and a currently employed eight-marker scheme. Of the 52 markers, ≥50 were successfully haplotyped for all specimens, and these results were used to generate a hierarchical cluster dendrogram. Using a previously described statistical approach to dissect hierarchical trees, the 66 specimens resolved into 24 and 27 distinct genetic clusters for the TAS and an 8-loci scheme, respectively. Although the specimen composition of 15 clusters was identical, there were substantial differences between the two dendrograms, highlighting the importance of both inclusion of additional genome coverage and choice of loci to target for genotyping. To evaluate the ability to genetically link contaminated food samples with clinical specimens, C. cayetanensis was genotyped from DNA extracted from raspberries inoculated with fecal specimens. The contaminated raspberry samples were assigned to clusters with the corresponding clinical specimen, demonstrating the utility of the TAS assay for traceback efforts.

Sources and Prevalence of Cyclospora cayetanensis in Southeastern U.S. Growing Environments.

Recent cyclosporiasis outbreaks associated with fresh produce grown in the United States highlight the need to better understand Cyclospora cayetanensis prevalence in U.S. agricultural environments. In this study, C. cayetanensis occurrence was assessed in municipal wastewater sludge, on-farm portable toilets, irrigation pond water, and spent packing house dump tank water in a Southeastern Georgia growing region over two years. Detection of the C. cayetanensis 18S rRNA qPCR gene target in pond samples was 0%, 28%, and 42% (N = 217) depending on the detection definition used, and ≤1% in dump tank samples (N = 46). However, no qPCR detections were confirmed by sequencing, suggesting false detection occurred due to cross-reactions. C. cayetanensis qPCR detections were confirmed in 9% of wastewater sludge samples (N = 76). The human-specific fecal markers HF183 and crAssphage were detected in 33% and 6% of pond samples, respectively, and 4% and 0% of dump tank samples, respectively. Despite community Cyclospora shedding and evidence of human fecal contamination in irrigation water, there was no correlation between C. cayetanensis and HF183 qPCR detections, further supporting that 18S gene target qPCR amplifications were due to cross-reactions. When evaluating C. cayetanensis qPCR environmental detection data, the impact of assay specificity and detection criteria should be considered. Moreover, additional sequence-based testing may be needed to appropriately interpret Cyclospora qPCR environmental data.

Optimizing hierarchical tree dissection parameters using historic epidemiologic data as 'ground truth'.

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.