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

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.

Sapp SG, Peterson A, Qvarnstrom Y. Comparison of two novel one-tube nested real-time qPCR assays to detect human-infecting Cyclospora spp

IMPORTANCE

Human-infecting Cyclospora spp. cause gastrointestinal distress among healthy individuals contributing to morbidity and putting stress on the economics of countries and companies in the form of produce recalls. Accessible and easy-to-use diagnostic tools available to a wide variety of laboratories would aid in the early detection of possible outbreaks of cyclosporiasis. This, in turn, will assist in the timely traceback investigation to the suspected source of an outbreak by informing the smallest possible recall and protecting consumers from contaminated produce. This manuscript describes two novel detection methods with improved performance for the causative agents of cyclosporiasis when compared to the currently used 18S assay.

Evaluation of nuclear and mitochondrial phylogenetics for the subtyping of Cyclospora cayetanensis

Cyclospora cayetanensis is an enteric coccidian parasite responsible for gastrointestinal disease transmitted through contaminated food and water. It has been documented in several countries, mostly with low-socioeconomic levels, although major outbreaks have hit developed countries. Detection methods based on oocyst morphology, staining, and molecular testing have been developed. However, the current MLST panel offers an opportunity for enhancement, as amplification of all molecular markers remains unfeasible in the majority of samples. This study aims to address this challenge by evaluating two approaches for analyzing the genetic diversity of C. cayetanensis and identifying reliable markers for subtyping: core homologous genes and mitochondrial genome analysis. A pangenome was constructed using 36 complete genomes of C. cayetanensis, and a haplotype network and phylogenetic analysis were conducted using 33 mitochondrial genomes. Through the analysis of the pangenome, 47 potential markers were identified, emphasizing the need for more sequence data to achieve comprehensive characterization. Additionally, the analysis of mitochondrial genomes revealed 19 single-nucleotide variations that can serve as characteristic markers for subtyping this parasite. These findings not only contribute to the selection of molecular markers for C. cayetanensis subtyping, but they also drive the knowledge toward the potential development of a comprehensive genotyping method for this parasite.

The limit of detection of the BioFire® FilmArray® gastrointestinal panel for the foodborne parasite Cyclospora cayetanensis

Cyclosporiasis is a foodborne diarrheal illness caused by the parasite Cyclospora cayetanensis. The BioFire® FilmArray® gastrointestinal (FilmArray GI) panel is a common method for diagnosing cyclosporiasis from clinical stool samples. The currently published limit of detection (LOD) of this panel is in genome equivalents; however, it is unclear how this relates to the number of C. cayetanensis oocysts in a clinical sample. In this study, we developed a technique to determine the LOD in terms of oocysts, using a cell sorter to sort 1 to 50 C. cayetanensis oocyst(s) previously purified from three human stool sources. We found the FilmArray GI panel detected samples with ≥20 C. cayetanensis oocysts in 100% of replicates, with varying detection among samples with 1, 5, or 10 C. cayetanensis oocysts. This method provides a parasitologically relevant LOD that should enable comparison among C. cayetanensis detection techniques, including the FilmArray GI panel.

Novel insights on the genetic population structure of human-infecting Cyclospora spp. and evidence for rapid subtype selection among isolates from the USA

Human-infecting Cyclospora was recently characterized as three species, two of which (C. cayetanensis and C. ashfordi) are currently responsible for all known human infections in the USA, yet much remains unknown about the genetic structure within these two species. Here, we investigate Cyclospora genotyping data from 2018 through 2022 to ascertain if there are temporal patterns in the genetic structure of Cyclospora parasites that cause infections in US residents from year to year. First, we investigate three levels of genetic characterization: species, subpopulation, and strain, to elucidate annual trends in Cyclospora infections. Next, we determine if shifts in genetic diversity can be linked to any of the eight loci used in our Cyclospora genotyping approach. We observed fluctuations in the abundance of Cyclospora types at the species and subpopulation levels, but no significant temporal trends were identified; however, we found recurrent and sporadic strains within both C. ashfordi and C. cayetanensis. We also uncovered major shifts in the mitochondrial genotypes in both species, where there was a universal increase in abundance of a specific mitochondrial genotype that was relatively abundant in 2018 but reached near fixation (was observed in over 96% of isolates) in C. ashfordi by 2022. Similarly, this allele jumped from 29% to 82% relative abundance of isolates belonging to C. cayetanensis. Overall, our analysis uncovers previously unknown temporal-genetic patterns in US Cyclospora types from 2018 through 2022 and is an important step to presenting a clearer picture of the factors influencing cyclosporiasis outbreaks in the USA.

An improved framework for detecting discrete epidemiologically meaningful partitions in hierarchically clustered genetic data

Motivation

Hierarchical clustering of microbial genotypes has the limitation that hierarchical clusters are nested, where smaller groups of related isolates exist within larger groups that get progressively larger as relationships become increasingly distant. In an epidemiologic context, investigators must dissect hierarchical trees into discrete groupings that are epidemiologically meaningful. We recently described a statistical framework (Method A) for dissecting hierarchical trees that attempts to minimize investigator bias. Here, we apply a modified version of that framework (Method B) to a hierarchical tree constructed from 2111 genotypes of the foodborne parasite Cyclospora, including 639 genotypes linked to epidemiologically defined outbreaks. To evaluate Method B's performance, we examined the concordance between these epidemiologically defined groupings and the genetic partitions identified. We also used the same epidemiologic clusters to evaluate the performance of Method A, plus two tree-dissection methods (cutreeHybrid and cutreeDynamic) available within the Dynamic Tree Cut R package, in addition to the TreeCluster method and PARNAS.

Results

Compared to the other methods, Method B, TreeCluster, and PARNAS were the most accurate (99.4%) in identifying genetic groups that reflected the epidemiologic groupings, noting that TreeCluster and PARNAS performed identically on our dataset. CutreeHybrid identified groups reflecting patterns in the wider Cyclospora population structure but lacked finer, strain-level discrimination (Simpson's D: cutreeHybrid=0.785). CutreeDynamic displayed good strain discrimination (Simpson's D = 0.933), though lacked sensitivity (77%). At two different threshold/radius settings TreeCluster/PARNAS displayed similar utility to Method B. However, Method B computes a tree-dissection threshold automatically, and the threshold/radius settings used when executing TreeCluster/PARNAS here were computed using Method B. Using a TreeCluster threshold of 0.045 as recommended in the TreeCluster documentation, epidemiologic utility dropped markedly below that of Method B.

Availability and implementation

Relevant code and data are publicly available. Source code (Method B) and instructions for its use are available here: https://github.com/Joel-Barratt/Hierarchical-tree-dissection-framework.

Cyclospora cayetanensis: A Perspective (2020-2023) with Emphasis on Epidemiology and Detection Methods

Cyclospora cayetanensis infections are prevalent worldwide, and the parasite has become a major public health and food safety concern. Although important efforts have been dedicated to advance toward preventing and reducing incidences of cyclosporiasis, there are still several knowledge gaps that hamper the implementation of effective measures to prevent the contamination of produce and water with Cyclospora oocysts. Some of these data gaps can be attributed to the fact that access to oocysts is a limiting factor in C. cayetanensis research. There are no animal models or in vivo or in vitro culture systems to propagate the oocysts needed to facilitate C. cayetanensis research. Thus, researchers must rely upon limited supplies of oocysts obtained from naturally infected human patients considerably restricting what can be learnt about this parasite. Despite the limited supply of C. cayetanensis oocysts, several important advances have happened in the past 3 years. Great progress has been made in the Cyclospora field in the areas of molecular characterization of strains and species, generation of genomes, and development of novel detection methods. This comprehensive perspective summarizes research published from 2020 to 2023 and evaluates what we have learnt and identifies those aspects in which further research is needed.

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.

Development of a targeted amplicon sequencing method for genotyping Cyclospora cayetanensis from fresh produce and clinical samples with enhanced genomic resolution and sensitivity

Outbreaks of cyclosporiasis, an enteric illness caused by the parasite Cyclospora cayetanensis, have been associated with consumption of various types of fresh produce. Although a method is in use for genotyping C. cayetanensis from clinical specimens, the very low abundance of C. cayetanensis in food and environmental samples presents a greater challenge. To complement epidemiological investigations, a molecular surveillance tool is needed for use in genetic linkage of food vehicles to cyclosporiasis illnesses, estimation of the scope of outbreaks or clusters of illness, and determination of geographical areas involved. We developed a targeted amplicon sequencing (TAS) assay that incorporates a further enrichment step to gain the requisite sensitivity for genotyping C. cayetanensis contaminating fresh produce samples. The TAS assay targets 52 loci, 49 of which are located in the nuclear genome, and encompasses 396 currently known SNP sites. The performance of the TAS assay was evaluated using lettuce, basil, cilantro, salad mix, and blackberries inoculated with C. cayetanensis oocysts. A minimum of 24 markers were haplotyped even at low contamination levels of 10 oocysts in 25 g leafy greens. The artificially contaminated fresh produce samples were included in a genetic distance analysis based on haplotype presence/absence with publicly available C. cayetanensis whole genome sequence assemblies. Oocysts from two different sources were used for inoculation, and samples receiving the same oocyst preparation clustered together, but separately from the other group, demonstrating the utility of the assay for genetically linking samples. Clinical fecal samples with low parasite loads were also successfully genotyped. This work represents a significant advance in the ability to genotype C. cayetanensis contaminating fresh produce along with greatly expanding the genomic diversity included for genetic clustering of clinical specimens.

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).

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.

Genetic characterization of Strongyloides fuelleborni infecting free-roaming African vervets (Chlorocebus aethiops sabaeus) on the Caribbean island of St. Kitts

Human strongyloidiasis is an important neglected tropical disease primarily caused by the nematode Strongyloides stercoralis, and to a lesser extent Strongyloides fuelleborni which mainly infects non-human primates. Zoonotic sources of infection have important implications for control and prevention of morbidity and mortality caused by strongyloidiasis. Recent molecular evidence suggests that for S. fuelleborni, primate host specificity is variable among genotypes across the Old World, and consequently that these types likely vary in their capacity for human spillover infections. Populations of free-roaming vervet monkeys (Chlorocebus aethiops sabaeus), introduced to the Caribbean Island of Staint Kitts from Africa, live in close contact with humans, and concern has arisen regarding their potential to serve as reservoirs of zoonotic infections. In this study, we sought to determine the genotypes of S. fuelleborni infecting St. Kitts vervets to explore whether they are potential reservoirs for human-infecting S. fuelleborni types. Fecal specimens were collected from St. Kitts vervets and S. fuelleborni infections were confirmed microscopically and by PCR. Strongyloides fuelleborni genotypes were determined from positive fecal specimens using an Illumina amplicon sequencing-based genotyping approach targeting the mitochondrial cox1 locus and 18S rDNA hypervariable regions I and IV of Strongyloides species. Phylogenetic analysis of resultant genotypes supported that S. fuelleborni from St. Kitts vervets is of an exclusively African variety, falling within the same monophyletic group as an isolate which has been detected previously in a naturally infected human from Guinea-Bissau. This observation highlights that St. Kitts vervets may serve as potential reservoirs for zoonotic S. fuelleborni infection, which warrants further exploration.

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.

Development of a Molecular Marker Based on the Mitochondrial Genome for Detection of Cyclospora cayetanensis in Food and Water Samples

Cyclospora cayetanensis is a coccidian parasite that causes diarrheal illness outbreaks worldwide. The development of new laboratory methods for detection of C. cayetanensis is of critical importance because of the high potential for environmental samples to be contaminated with a myriad of microorganisms, adversely impacting the specificity when testing samples from various sources using a single molecular assay. In this study, a new sequencing-based method was designed targeting a specific fragment of C. cayetanensis cytochrome oxidase gene and developed as a complementary method to the TaqMan qPCR present in the U.S. FDA BAM Chapter 19b and Chapter 19c. The comparative results between the new PCR protocol and the qPCR for detection of C. cayetanensis in food and water samples provided similar results in both matrices with the same seeding level. The target region and primers in the protocol discussed in this study contain sufficient Cyclospora-specific sequence fidelity as observed by sequence comparison with other Eimeriidae species. The sequence of the PCR product appears to represent a robust target for identifying C. cayetanensis on samples from different sources. Such a sensitive method for detection of C. cayetanensis would add to the target repertoire of qPCR-based screening strategies for food and water samples.

Genotyping Canadian Cyclospora cayetanensis Isolates to Supplement Cyclosporiasis Outbreak Investigations

Cyclospora cayetanensis is an emerging foodborne parasite that causes cyclosporiasis, an enteric disease of humans. Domestically acquired outbreaks have been reported in Canada every spring or summer since 2013. To date, investigations into the potential sources of infection have relied solely on epidemiological data. To supplement the epidemiological data with genetic information, we genotyped 169 Canadian cyclosporiasis cases from stool specimens collected from 2010 to 2021 using an existing eight-marker targeted amplicon deep (TADS) scheme specific to C. cayetanensis as previously described by the US Centers for Disease Control and Prevention (CDC). This is the first study to genotype Canadian Cyclospora cayetanensis isolates, and it focuses on evaluating the genotyping performance and genetic clustering. Genotyping information was successfully collected with at least part of one of the markers in the TADS assay for 97.9% of specimens, and 81.1% of cyclosporiasis cases met the minimum requirements to genetically cluster into 20 groups. The performance of the scheme suggests that examining cyclosporiasis cases genetically will be a valuable tool for supplementing epidemiological outbreak investigations and to minimize further infections. Further research is required to expand the number of discriminatory markers to improve genetic clustering.

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.