Targeted next generation sequencing of Cyclospora cayetanensis mitochondrial genomes from seeded fresh produce and other seeded food samples.

Cyclospora cayetanensis, a coccidian apicomplexan parasite, causes large outbreaks of foodborne diarrheal disease globally. Tracking the source of C. cayetanensis oocyst contamination in food items is essential to reduce, even prevent outbreaks. We previously showed that a genotyping method based on mitochondrial single nucleotide polymorphism (SNP) profiles had discriminatory power in classifying C. cayetanensis clinical isolates. In food specimens, low level contamination by oocysts and difficulties in DNA extraction present significant challenges in genotyping method development. Here, we report the development of a highly sensitive, custom-designed, targeted sequencing method based on the Illumina AmpliSeq platform; our method was capable of consistently generating near-complete mitochondrial genome sequences of C. cayetanensis from foods with low levels of contamination. To simulate environmentally observed contamination levels in foods, we seeded various food matrices, such as fresh produce and prepared dishes, with known quantities of oocysts, and isolated genomic DNA from washed food samples. Using the Ampliseq Targeted Sequencing method, we obtained near-complete mitochondrial genome sequences of C. cayetanensis from food samples seeded with as low as five to ten oocysts and used the data in downstream analysis. The flexibility of the AmpliSeq platform could potentially allow for more genomic targets to be added to achieve higher discriminatory power. This level of sensitivity in capturing high resolution genome data from contaminated food samples is a critical milestone towards the potential development of a comprehensive genotyping method for C. cayetanensis.

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.

Assessment of Commercial DNA Cleanup Kits for Elimination of Real-Time PCR Inhibitors in the Detection of Cyclospora cayetanensis in Cilantro.

ABSTRACT: Inhibited reactions have occasionally been observed when cilantro samples were processed for the detection of Cyclospora cayetanensis using quantitative real-time PCR (qPCR). Partial or total inhibition of PCR reactions, including qPCR, can occur, leading to decreased sensitivity or false-negative results. If inhibition occurs, this implies the need for additional purification or cleanup treatments of the extracted DNA to remove inhibitors prior to molecular detection. Our objective was to evaluate the performance of five commercial DNA cleanup kits (QIAquick purification kit from Qiagen [kit 1], OneStep PCR inhibitor removal by Zymo Research [kit 2], NucleoSpin genomic DNA cleanup XS from Macherey-Nagel [kit 3], DNA IQ system by Promega [kit 4], and DNeasy PowerPlant pro kit from Qiagen [5]) to minimize qPCR inhibition using the U.S. Food and Drug Administration-validated Bacteriological Analytical Manual (BAM) Chapter 19b method for detection of C. cayetanensis in cilantro samples containing soil. Each of the five commercial DNA cleanup kits evaluated was able to reduce the qPCR internal amplification control cycle threshold values to those considered to be normal for noninhibited samples, allowing unambiguous interpretation of results in cilantro samples seeded at both a high oocyst level (200 oocysts) and a low oocyst level (10 oocysts). Of the five kits compared, kits 1, 2, and 3 did not show significant differences in the detection of C. cayetanensis, while significantly higher cycle threshold values, indicating lower recovery of the target DNA, were observed from kits 4 and/or 5 in samples seeded with 200 and 10 oocysts (P < 0.05). This comparative study provides recommendations on the use of commercial cleanup kits which could be implemented when inhibition is observed in the detection of C. cayetanensis in cilantro samples using the BAM Chapter 19b method.

Molecular typing of Cyclospora cayetanensis in produce and clinical samples using targeted enrichment of complete mitochondrial genomes and next-generation sequencing.

BACKGROUND: Outbreaks of cyclosporiasis, a diarrheal illness caused by Cyclospora cayetanensis, have been a public health issue in the USA since the mid 1990's. In 2018, 2299 domestically acquired cases of cyclosporiasis were reported in the USA as a result of multiple large outbreaks linked to different fresh produce commodities. Outbreak investigations are hindered by the absence of standardized molecular epidemiological tools for C. cayetanensis. For other apicomplexan coccidian parasites, multicopy organellar DNA such as mitochondrial genomes have been used for detection and molecular typing. METHODS: We developed a workflow to obtain complete mitochondrial genome sequences from cilantro samples and clinical samples for typing of C. cayetanensis isolates. The 6.3 kb long C. cayetanensis mitochondrial genome was amplified by PCR in four overlapping amplicons from genomic DNA extracted from cilantro, seeded with oocysts, and from stool samples positive for C. cayetanensis by diagnostic methods. DNA sequence libraries of pooled amplicons were prepared and sequenced via next-generation sequencing (NGS). Sequence reads were assembled using a custom bioinformatics pipeline. RESULTS: This approach allowed us to sequence complete mitochondrial genomes from the samples studied. Sequence alterations, such as single nucleotide polymorphism (SNP) profiles and insertion and deletions (InDels), in mitochondrial genomes of 24 stool samples from patients with cyclosporiasis diagnosed in 2014, exhibited discriminatory power. The cluster dendrogram that was created based on distance matrices of the complete mitochondrial genome sequences, indicated distinct strain-level diversity among the 2014 C. cayetanensis outbreak isolates analyzed in this study. CONCLUSIONS: Our results suggest that genomic analyses of mitochondrial genome sequences may help to link outbreak cases to the source.

Evaluation of the U.S. Food and Drug Administration validated molecular method for detection of Cyclospora cayetanensis oocysts on fresh and frozen berries.

Outbreaks and sporadic cases of Cyclospora cayetanensis have been linked to consumption of berries. The efficacy of the U.S. Food and Drug Administration (FDA) method for detection of C. cayetanensis was evaluated in fresh berries (blackberries, strawberries, blueberries and mixed berries) and in frozen mixed berries. The protocol included seeding with C. cayetanensis oocysts, produce washing, DNA extraction and a dual TaqMan assay. As few as five oocysts were detected in every type of fresh berry analyzed. All berry samples seeded with 200 oocysts were positive and all unseeded berry samples were negative. No significant differences were observed among any of the berry types analyzed in detection rates, CT values and estimated oocyst recovery percentages. Mixed berries were seeded and frozen for up to seven weeks. As few as five oocysts were also detected. No significant differences were observed in C. cayetanensis CT values between fresh and frozen mixed berries at any seeding level. In conclusion, the FDA BAM Chapter 19B method for the detection of Cyclospora was robust, consistent, and showed high sensitivity in all types of berries analyzed. Evaluation of the FDA detection method in berries will provide reliable laboratory support for surveillance programs and for outbreak investigations.

Evaluation of the U.S. Food and Drug Administration validated method for detection of Cyclospora cayetanensis in high-risk fresh produce matrices and a method modification for a prepared dish.

The performance of the U.S. Food and Drug Administration (FDA) validated method for regulatory detection of Cyclospora cayetanensis in leafy greens and berries was evaluated in additional high-risk fresh produce items and in a dish prepared with these produce commodities. The method was robust and reproducible in basil, parsley, shredded carrots, shredded cabbage and carrot mix, and could detect as few as 5 oocysts in 25 g samples. Some differences in C. cayetanensis detection were found among the fresh produce analyzed. Significantly lower target gene copy numbers per reaction were obtained with shredded carrots, and shredded cabbage and carrot mix compared to leafy greens, which highlights the importance of evaluating the performance characteristics of validated methods in different food matrices. In the prepared dish, coleslaw with dressing, the method was optimized to detect 5 oocysts in a 25 g sample by using 1.0% Alconox® in the washing solution instead of 0.1% as originally described. These data are important to assess the prevalence of C. cayetanensis in different produce items and to support outbreak investigations.

Interlaboratory validation of an improved method for detection of Cyclospora cayetanensis in produce using a real-time PCR assay.

A collaborative validation study was performed to evaluate the performance of a new U.S. Food and Drug Administration method developed for detection of the protozoan parasite, Cyclospora cayetanensis, on cilantro and raspberries. The method includes a sample preparation step in which oocysts are recovered from produce using an enhanced produce washing solution containing 0.1% Alconox and a commercially available method to disrupt the C. cayetanensis oocysts and extract DNA. A real-time PCR assay targeting the C. cayetanensis 18S rDNA gene with an internal amplification control to monitor PCR inhibition provides species-specific identification. Five laboratories blindly analyzed a total of 319 samples consisting of 25 g of cilantro or 50 g of raspberries which were either uninoculated or artificially contaminated with C. cayetanensis oocysts. Detection rates for cilantro inoculated with 200, 10, and 5 oocysts, were 100%, 80%, and 31%, respectively. For raspberries, the detection rates for samples inoculated with 200, 10, and 5 oocysts were 100%, 90% and 50%, respectively. All uninoculated samples, DNA blank extracts, and no-template PCR controls were negative. Reproducibility between laboratories and analysts was high and the method was shown to be an effective analytical tool for detection of C. cayetanensis in produce.

Comparative sequence analysis of Cyclospora cayetanensis apicoplast genomes originating from diverse geographical regions.

BACKGROUND: Cyclospora cayetanensis is an emerging coccidian parasite that causes endemic and epidemic diarrheal disease called cyclosporiasis, and this infection is associated with consumption of contaminated produce or water in developed and developing regions. Food-borne outbreaks of cyclosporiasis have occurred almost every year in the USA since the 1990s. Investigations of these outbreaks are currently hampered due to lack of molecular epidemiological tools for trace back analysis. The apicoplast of C. cayetanensis, a relict non-photosynthetic plastid with an independent genome, provides an attractive target to discover sequence polymorphisms useful as genetic markers for detection and trace back analysis of the parasite. Distinct differences in the apicoplast genomes of C. cayetanensis could be useful in designing advanced molecular methods for rapid detection and, subtyping and geographical source attribution, which would aid outbreak investigations and surveillance studies. METHODS: To obtain the genome sequence of the C. cayetanensis apicoplast, we sequenced the C. cayetanensis genomic DNA extracted from clinical stool samples, assembled and annotated a 34,146 bp-long circular sequence, and used this sequence as a reference genome in this study. We compared the genome and the predicted proteome to the data available from other apicomplexan parasites. To initialize the search for genetic markers, we mapped the raw sequence reads from an additional 11 distinct clinical stool samples originating from Nepal, New York, Texas, and Indonesia to the apicoplast reference genome. RESULTS: We identified several high quality single nucleotide polymorphisms (SNPs) and small insertion/deletions spanning the apicoplast genome supported by extensive sequencing reads data, and a 30 bp sequence repeat at the terminal spacer region in a Nepalese sample. The predicted proteome consists of 29 core apicomplexan peptides found in most of the apicomplexans. Cluster analysis of these C. cayetanensis apicoplast genomes revealed a familiar pattern of tight grouping with Eimeria and Toxoplasma, separated from distant species such as Plasmodium and Babesia. CONCLUSIONS: SNPs and sequence repeats identified in this study may be useful as genetic markers for identification and differentiation of C. cayetanensis isolates found and could facilitate outbreak investigations.

The Complete Mitochondrial Genome of the Foodborne Parasitic Pathogen Cyclospora cayetanensis.

Cyclospora cayetanensis is a human-specific coccidian parasite responsible for several food and water-related outbreaks around the world, including the most recent ones involving over 900 persons in 2013 and 2014 outbreaks in the USA. Multicopy organellar DNA such as mitochondrion genomes have been particularly informative for detection and genetic traceback analysis in other parasites. We sequenced the C. cayetanensis genomic DNA obtained from stool samples from patients infected with Cyclospora in Nepal using the Illumina MiSeq platform. By bioinformatically filtering out the metagenomic reads of non-coccidian origin sequences and concentrating the reads by targeted alignment, we were able to obtain contigs containing Eimeria-like mitochondrial, apicoplastic and some chromosomal genomic fragments. A mitochondrial genomic sequence was assembled and confirmed by cloning and sequencing targeted PCR products amplified from Cyclospora DNA using primers based on our draft assembly sequence. The results show that the C. cayetanensis mitochondrion genome is 6274 bp in length, with 33% GC content, and likely exists in concatemeric arrays as in Eimeria mitochondrial genomes. Phylogenetic analysis of the C. cayetanensis mitochondrial genome places this organism in a tight cluster with Eimeria species. The mitochondrial genome of C. cayetanensis contains three protein coding genes, cytochrome (cytb), cytochrome C oxidase subunit 1 (cox1), and cytochrome C oxidase subunit 3 (cox3), in addition to 14 large subunit (LSU) and nine small subunit (SSU) fragmented rRNA genes.

Genomic analysis of stress response against arsenic in Caenorhabditis elegans.

Arsenic, a known human carcinogen, is widely distributed around the world and found in particularly high concentrations in certain regions including Southwestern US, Eastern Europe, India, China, Taiwan and Mexico. Chronic arsenic poisoning affects millions of people worldwide and is associated with increased risk of many diseases including arthrosclerosis, diabetes and cancer. In this study, we explored genome level global responses to high and low levels of arsenic exposure in Caenorhabditis elegans using Affymetrix expression microarrays. This experimental design allows us to do microarray analysis of dose-response relationships of global gene expression patterns. High dose (0.03%) exposure caused stronger global gene expression changes in comparison with low dose (0.003%) exposure, suggesting a positive dose-response correlation. Biological processes such as oxidative stress, and iron metabolism, which were previously reported to be involved in arsenic toxicity studies using cultured cells, experimental animals, and humans, were found to be affected in C. elegans. We performed genome-wide gene expression comparisons between our microarray data and publicly available C. elegans microarray datasets of cadmium, and sediment exposure samples of German rivers Rhine and Elbe. Bioinformatics analysis of arsenic-responsive regulatory networks were done using FastMEDUSA program. FastMEDUSA analysis identified cancer-related genes, particularly genes associated with leukemia, such as dnj-11, which encodes a protein orthologous to the mammalian ZRF1/MIDA1/MPP11/DNAJC2 family of ribosome-associated molecular chaperones. We analyzed the protective functions of several of the identified genes using RNAi. Our study indicates that C. elegans could be a substitute model to study the mechanism of metal toxicity using high-throughput expression data and bioinformatics tools such as FastMEDUSA.

Genomic analysis of immune response against Vibrio cholerae hemolysin in Caenorhabditis elegans.

Vibrio cholerae cytolysin (VCC) is among the accessory V. cholerae virulence factors that may contribute to disease pathogenesis in humans. VCC, encoded by hlyA gene, belongs to the most common class of bacterial toxins, known as pore-forming toxins (PFTs). V. cholerae infects and kills Caenorhabditis elegans via cholerae toxin independent manner. VCC is required for the lethality, growth retardation and intestinal cell vacuolation during the infection. However, little is known about the host gene expression responses against VCC. To address this question we performed a microarray study in C. elegans exposed to V. cholerae strains with intact and deleted hlyA genes.Many of the VCC regulated genes identified, including C-type lectins, Prion-like (glutamine [Q]/asparagine [N]-rich)-domain containing genes, genes regulated by insulin/IGF-1-mediated signaling (IIS) pathway, were previously reported as mediators of innate immune response against other bacteria in C. elegans. Protective function of the subset of the genes up-regulated by VCC was confirmed using RNAi. By means of a machine learning algorithm called FastMEDUSA, we identified several putative VCC induced immune regulatory transcriptional factors and transcription factor binding motifs. Our results suggest that VCC is a major virulence factor, which induces a wide variety of immune response- related genes during V. cholerae infection in C. elegans.

Vibrio cholerae hemolysin is required for lethality, developmental delay, and intestinal vacuolation in Caenorhabditis elegans.

BACKGROUND: Cholera toxin (CT) and toxin-co-regulated pili (TCP) are the major virulence factors of Vibrio cholerae O1 and O139 strains that contribute to the pathogenesis of disease during devastating cholera pandemics. However, CT and TCP negative V. cholerae strains are still able to cause severe diarrheal disease in humans through mechanisms that are not well understood. METHODOLOGY/PRINCIPAL FINDINGS: To determine the role of other virulence factors in V. cholerae pathogenesis, we used a CT and TCP independent infection model in the nematode Caenorhabditis elegans and identified the hemolysin A (hlyA) gene as a factor responsible for animal death and developmental delay. We demonstrated a correlation between the severity of infection in the nematode and the level of hemolytic activity in the V. cholerae biotypes. At the cellular level, V. cholerae infection induces formation of vacuoles in the intestinal cells in a hlyA dependent manner, consistent with the previous in vitro observations. CONCLUSIONS/SIGNIFICANCE: Our data strongly suggest that HlyA is a virulence factor in C. elegans infection leading to lethality and developmental delay presumably through intestinal cytopathic changes.

A method to rank order water soluble compounds according to their toxicity using Caenorhabditis elegans, a Complex Object Parametric Analyzer and Sorter, and axenic liquid media.

Complex Object Parametric Analyzer and Sorter (COPAS) parameters Time of Flight (TOF) and Extinction (EXT) were utilized to assess growth and development in Caenorhabditis elegans exposed to (in order of decreasing toxicity) sodium arsenite, sodium fluoride, caffeine, valproic acid, sodium borate or DMSO in C. elegans Habitation Medium (CeHM) for 72h. Using multivariate statistical modeling and unique sub sampling procedures mean p-value ratios were calculated for each compound. Comparison of mean p-value ratios and/or the percent change in mean-p value ratios to controls were utilized to assess test compound toxicity. Using this assay 5 of the 6 compounds tested (83.3%) were correctly ranked according to their toxicity based on oral rat LD50 data. Test compounds were ranked from most toxic to least toxic as follows: sodium arsenite, sodium fluoride, sodium borate, valproic acid, caffeine and DMSO. Sodium borate was found to be more toxic than caffeine and valproic acid in this bioassay. This study suggests that axenic liquid culture may be used to expose large numbers of nematodes to water soluble toxicants and the COPAS parameters TOF and EXT may be used as functional biomarkers to assess a toxin's effect on growth and development in C. elegans.

N-ethylmaleimide sensitive factor is required for fusion of the C. elegans uterine anchor cell.

The fusion of the Caenorhabditis elegans uterine anchor cell (AC) with the uterine-seam cell (utse) is an excellent model system for studying cell-cell fusion, which is essential to animal development. We obtained an egg-laying defective (Egl) mutant in which the AC fails to fuse with the utse. This defect is highly specific: other aspects of utse development and other cell fusions appear to occur normally. We find that defect is due to a missense mutation in the nsf-1 gene, which encodes N-ethylmaleimide-sensitive factor (NSF), an intracellular membrane fusion factor. There are two NSF-1 isoforms, which are expressed in distinct tissues through two separate promoters. NSF-1L is expressed in the uterus, including the AC. We find that nsf-1 is required cell-autonomously in the AC for its fusion with the utse. Our results establish AC fusion as a paradigm for studying cell fusion at single cell resolution and demonstrate that the NSF ATPase is a key player in this process.

C. elegans Kallmann syndrome protein KAL-1 interacts with syndecan and glypican to regulate neuronal cell migrations.

The anosmin-1 protein family regulates cell migration, axon guidance, and branching, by mechanisms that are not well understood. We show that the C. elegans anosmin-1 ortholog KAL-1 promotes migrations of ventral neuroblasts prior to epidermal enclosure. KAL-1 does not modulate FGF signaling in neuroblast migration and acts in parallel to other neuroblast migration pathways. Defects in heparan sulfate (HS) synthesis or in specific HS modifications disrupt neuroblast migrations and affect the KAL-1 pathway. KAL-1 binds the cell surface HS proteoglycans syndecan/SDN-1 and glypican/GPN-1. This interaction is mediated via HS side chains and requires specific HS modifications. SDN-1 and GPN-1 are expressed in ventral neuroblasts and have redundant roles in KAL-1-dependent neuroblast migrations. Our findings suggest that KAL-1 interacts with multiple HSPGs to promote cell migration.

Neurosurgery: functional regeneration after laser axotomy.

Understanding how nerves regenerate is an important step towards developing treatments for human neurological disease, but investigation has so far been limited to complex organisms (mouse and zebrafish) in the absence of precision techniques for severing axons (axotomy). Here we use femtosecond laser surgery for axotomy in the roundworm Caenorhabditis elegans and show that these axons functionally regenerate after the operation. Application of this precise surgical technique should enable nerve regeneration to be studied in vivo in its most evolutionarily simple form.

Genetic analysis of the Caenorhabditis elegans pax-6 locus: roles of paired domain-containing and nonpaired domain-containing isoforms.

PAX-6 proteins are involved in eye and brain development in many animals. In the nematode Caenorhabditis elegans the pax-6 locus encodes multiple PAX-6 isoforms both with and without a paired domain. Mutations in the C. elegans pax-6 locus can be grouped into three classes. Mutations that affect paired domain-containing isoforms cause defects in epidermal morphogenesis, epidermal cell fates, and gonad cell migration and define the class I (vab-3) complementation group. The class II mutation mab-18(bx23) affects nonpaired domain-containing isoforms and transforms the fate of a sensory organ in the male tail. Class III mutations affect both paired domain and nonpaired domain isoforms; the most severe class III mutations are candidate null mutations in pax-6. Class III mutant phenotypes do not resemble a simple sum of class I and class II phenotypes. A comparison of class I and class III phenotypes indicates that PAX-6 isoforms can interact additively, synergistically, or antagonistically, depending on the cellular context.

The EGL-13 SOX domain transcription factor affects the uterine pi cell lineages in Caenorhabditis elegans.

We isolated egl-13 mutants in which the pi cells of the Caenorhabditis elegans uterus initially appeared to develop normally but then underwent an extra round of cell division. The data suggest that egl-13 is required for maintenance of the pi cell fate.