Catherine D. Carrillo

Large-Scale Comparative Genomics Meta-Analysis of Campylobacter jejuni Isolates Reveals Low Level of Genome Plasticity

ABSTRACT We have used comparative genomic hybridization (CGH) on a full-genome Campylobacter jejuni microarray to examine genome-wide gene conservation patterns among 51 strains isolated from food and clinical sources. These data have been integrated with data from three previous C. jejuni CGH studies to perform a meta-analysis that included 97 strains from the four separate data sets. Although many genes were found to be divergent across multiple strains (n = 350), many genes (n = 249) were uniquely variable in single strains. Thus, the strains in each data set comprise strains with a unique genetic diversity not found in the strains in the other data sets. Despite the large increase in the collective number of variable C. jejuni genes (n = 599) found in the meta-analysis data set, nearly half of these (n = 276) mapped to previously defined variable loci, and it therefore appears that large regions of the C. jejuni genome are genetically stable. A detailed analysis of the microarray data revealed that divergent genes could be differentiated on the basis of the amplitudes of their differential microarray signals. Of 599 variable genes, 122 could be classified as highly divergent on the basis of CGH data. Nearly all highly divergent genes (117 of 122) had divergent neighbors and showed high levels of intraspecies variability. The approach outlined here has enabled us to distinguish global trends of gene conservation in C. jejuni and has enabled us to define this group of genes as a robust set of variable markers that can become the cornerstone of a new generation of genotyping methods that use genome-wide C. jejuni gene variability data.

Modulation of Toxin Production by the Flagellar Regulon in Clostridium difficile

ABSTRACT We show in this study that toxin production in Clostridium difficile is altered in cells which can no longer form flagellar filaments. The impact of inactivation of fliC, CD0240, fliF, fliG, fliM, and flhB-fliR flagellar genes upon toxin levels in culture supernatants was assessed using cell-based cytotoxicity assay, proteomics, immunoassay, and immunoblotting approaches. Each of these showed that toxin levels in supernatants were significantly increased in a fliC mutant compared to that in the C. difficile 630 parent strain. In contrast, the toxin levels in supernatants secreted from other flagellar mutants were significantly reduced compared with that in the parental C. difficile 630 strain. Transcriptional analysis of the pathogenicity locus genes (tcdR, tcdB, tcdE, and tcdA) revealed a significant increase of all four genes in the fliC mutant strain, while transcription of all four genes was significantly reduced in fliM, fliF, fliG, and flhB-fliR mutants. These results demonstrate that toxin transcription in C. difficile is modulated by the flagellar regulon. More significantly, mutant strains showed a corresponding change in virulence compared to the 630 parent strain when tested in a hamster model of C. difficile infection. This is the first demonstration of differential flagellum-related transcriptional regulation of toxin production in C. difficile and provides evidence for elaborate regulatory networks for virulence genes in C. difficile.

Current methods for molecular typing of Campylobacter species

Campylobacter remains one of the most common bacterial causes of gastroenteritis worldwide. Tracking sources of this organism is challenging due to the large numbers of human cases, and the prevalence of this organism throughout the environment due to growth in a wide range of animal species. Many molecular subtyping methods have been developed to characterize Campylobacter species, but only a few are commonly used in molecular epidemiology studies. This review examines the applicability of these methods, as well as the role that emerging whole genome sequencing technologies will play in tracking sources of Campylobacter spp. infection.

Comparative genomics profiling of clinical isolates of Aeromonas salmonicida using DNA microarrays

BackgroundAeromonas salmonicida has been isolated from numerous fish species and shows wide variation in virulence and pathogenicity. As part of a larger research program to identify virulence genes and candidates for vaccine development, a DNA microarray was constructed using a subset of 2024 genes from the draft genome sequence of A. salmonicida subsp. salmonicida strain A449. The microarray included genes encoding known virulence-associated factors in A. salmonicida and homologs of virulence genes of other pathogens. We used microarray-based comparative genomic hybridizations (M-CGH) to compare selected A. salmonicida sub-species and other Aeromonas species from different hosts and geographic locations.ResultsResults showed variable carriage of virulence-associated genes and generally increased variation in gene content across sub-species and species boundaries. The greatest variation was observed among genes associated with plasmids and transposons. There was little correlation between geographic region and degree of variation for all isolates tested.ConclusionWe have used the M-CGH technique to identify subsets of conserved genes from amongst this set of A. salmonicida virulence genes for further investigation as potential vaccine candidates. Unlike other bacterial characterization methods that use a small number of gene or DNA-based functions, M-CGH examines thousands of genes and/or whole genomes and thus is a more comprehensive analytical tool for veterinary or even human health research.

Occurrence of Campylobacter spp. in Raw and Ready-to-Eat Foods and in a Canadian Food Service Operation

The occurrence of Campylobacter spp. in a variety of foods from Ottawa, Ontario, Canada, and raw milk samples from across Canada was determined over a 2-year period. The samples consisted of 55 raw foods (chicken, pork, and beef), 126 raw milk samples from raw milk cheese manufacturers, and 135 ready-to-eat foods (meat products, salads, and raw milk cheeses). Campylobacter jejuni was detected in 4 of the 316 samples analyzed: 1 raw beef liver sample and 3 raw chicken samples. An isolation rate of 9.7% was observed among the raw chicken samples tested. This study also investigated the role of cross-contamination in disseminating Campylobacter from raw poultry within a food service operation specializing in poultry dishes. Accordingly, kitchen surfaces within a restaurant in Ottawa, Ontario, were sampled between March and August 2001. Tests of the sampling method indicated that as few as 100 Campylobacter cells could be detected if sampling was done within 45 min of inoculation; however, Campylobacter spp. were not detected in 125 swabs of surfaces within the kitchens of this food service operation. Despite the reported high prevalence of Campylobacter spp. in raw poultry, this organism was not detected on surfaces within a kitchen of a restaurant specializing in poultry dishes.

A framework for assessing the concordance of molecular typing methods and the true strain phylogeny of Campylobacter jejuni and C. coli using draft genome sequence data

Tracking of sources of sporadic cases of campylobacteriosis remains challenging, as commonly used molecular typing methods have limited ability to unambiguously link genetically related strains. Genomics has become increasingly prominent in the public health response to enteric pathogens as methods enable characterization of pathogens at an unprecedented level of resolution. However, the cost of sequencing and expertise required for bioinformatic analyses remains prohibitive, and these comprehensive analyses are limited to a few priority strains. Although several molecular typing methods are currently widely used for epidemiological analysis of campylobacters, it is not clear how accurately these methods reflect true strain relationships. To address this, we have developed a framework and associated computational tools to rapidly analyze draft genome sequence data for the assessment of molecular typing methods against a “gold standard” based on the phylogenetic analysis of highly conserved core (HCC) genes with high sequence quality. We analyzed 104 publicly available whole genome sequences (WGS) of C. jejuni and C. coli. In addition to in silico determination of multi-locus sequence typing (MLST), flaA, and porA type, as well as comparative genomic fingerprinting (CGF) type, we inferred a “reference” phylogeny based on 389 HCC genes. Molecular typing data were compared to the reference phylogeny for concordance using the adjusted Wallace coefficient (AWC) with confidence intervals. Although MLST targets the sequence variability in core genes and CGF targets insertions/deletions of accessory genes, both methods are based on multi-locus analysis and provided better estimates of true phylogeny than methods based on single loci (porA, flaA). A more comprehensive WGS dataset including additional genetically related strains, both epidemiologically linked and unlinked, will be necessary to more comprehensively assess the performance of subtyping methods for outbreak investigations and surveillance activities. Analyses of the strengths and weaknesses of widely used typing methodologies in inferring true strain relationships will provide guidance in the interpretation of this data for epidemiological purposes.

An enhanced technique combining pre-enrichment and passive filtration increases the isolation efficiency of Campylobacter jejuni and Campylobacter coli from water and animal fecal samples.

Improved isolation techniques from environmental water and animal samples are vital to understanding Campylobacter epidemiology. In this study, the efficiency of selective enrichment in Bolton Broth (BB) followed by plating on charcoal cefoperazone deoxycholate agar (CCDA) (conventional method) was compared with an approach combining BB enrichment and passive filtration (membrane method) adapted from a method previously developed for testing of broiler meat, in the isolation of thermophilic campylobacters from surface water and animal fecal samples. The conventional method led to recoveries of Campylobacter from 36.7% of the water samples and 78.0% of the fecal samples and similar numbers, 38.3% and 76.0%, respectively, were obtained with the membrane method. To investigate the genetic diversity of Campylobacter jejuni and Campylobacter coli obtained by these two methods, isolates were analyzed using Comparative Genomic Fingerprinting, a high-resolution subtyping technique. The conventional and membrane methods yielded similar numbers of Campylobacter subtypes from water (25 and 28, respectively) and fecal (15 and 17, respectively) samples. Although there was no significant difference in recovery rates between the conventional and membrane methods, a significant improvement in isolation efficiency was obtained by using the membrane method, with a false-positive rate of 1.6% compared with 30.7% obtained using the conventional method. In conclusion, although the two methods are comparable in sensitivity, the membrane method had higher specificity, making it a cost-effective procedure for the enhanced isolation of C. jejuni and C. coli from water and animal fecal samples.

GeneSippr: a rapid whole-genome approach for the identification and characterization of foodborne pathogens such as priority Shiga toxigenic Escherichia coli.

The timely identification and characterization of foodborne bacteria for risk assessment purposes is a key operation in outbreak investigations. Current methods require several days and/or provide low-resolution characterization. Here we describe a whole-genome-sequencing (WGS) approach (GeneSippr) enabling same-day identification of colony isolates recovered from investigative food samples. The identification of colonies of priority Shiga-toxigenic Escherichia coli (STEC) (i.e., serogroups O26, O45, O103, O111, O121, O145 and O157) served as a proof of concept. Genomic DNA was isolated from single colonies and sequencing was conducted on the Illumina MiSeq instrument with raw data sampling from the instrument following 4.5 hrs of sequencing. Modeling experiments indicated that datasets comprised of 21-nt reads representing approximately 4-fold coverage of the genome were sufficient to avoid significant gaps in sequence data. A novel bioinformatic pipeline was used to identify the presence of specific marker genes based on mapping of the short reads to reference sequence libraries, along with the detection of dispersed conserved genomic markers as a quality control metric to assure the validity of the analysis. STEC virulence markers were correctly identified in all isolates tested, and single colonies were identified within 9 hrs. This method has the potential to produce high-resolution characterization of STEC isolates, and whole-genome sequence data generated following the GeneSippr analysis could be used for isolate identification in place of lengthy biochemical characterization and typing methodologies. Significant advantages of this procedure include ease of adaptation to the detection of any gene marker of interest, as well as to the identification of other foodborne pathogens for which genomic markers have been defined.

Local genes for local bacteria: Evidence of allopatry in the genomes of transatlantic Campylobacter populations

The genetic structure of bacterial populations can be related to geographical locations of isolation. In some species, there is a strong correlation between geographical distance and genetic distance, which can be caused by different evolutionary mechanisms. Patterns of ancient admixture in Helicobacter pylori can be reconstructed in concordance with past human migration, whereas in Mycobacterium tuberculosis it is the lack of recombination that causes allopatric clusters. In Campylobacter, analyses of genomic data and molecular typing have been successful in determining the reservoir host species, but not geographical origin. We investigated biogeographical variation in highly recombining genes to determine the extent of clustering between genomes from geographically distinct Campylobacter populations. Whole‐genome sequences from 294 Campylobacter isolates from North America and the UK were analysed. Isolates from within the same country shared more recently recombined DNA than isolates from different countries. Using 15 UK/American closely matched pairs of isolates that shared ancestors, we identify regions that have frequently and recently recombined to test their correlation with geographical origin. The seven genes that demonstrated the greatest clustering by geography were used in an attribution model to infer geographical origin which was tested using a further 383 UK clinical isolates to detect signatures of recent foreign travel. Patient records indicated that in 46 cases, travel abroad had occurred <2 weeks prior to sampling, and genomic analysis identified that 34 (74%) of these isolates were of a non‐UK origin. Identification of biogeographical markers in Campylobacter genomes will contribute to improved source attribution of clinical Campylobacter infection and inform intervention strategies to reduce campylobacteriosis.

Genomic Tools for Customized Recovery and Detection of Foodborne Shiga Toxigenic Escherichia coli

Genomic antimicrobial resistance (AMR) prediction tools have the potential to support foodborne illness outbreak investigations through their application in the analysis of bacterial genomes from causative strains. The AMR marker profile of a strain of interest, initially identified in outbreak-associated clinical samples, may serve as the basis for customization of selective enrichment media, facilitating its recovery from samples in a food safety investigation. Different possibilities for AMR analyses include the use of comprehensive AMR gene databases such as the Comprehensive Antibiotic Resistance Database, which can be mined with in-house bioinformatics alignment tools (e.g., Antimicrobial Resistance Marker Identifier), or publicly available tools based on clinically relevant acquired AMR gene databases (e.g., ResFinder). In combination with a previously reported pipeline (SigSeekr) designed to identify specific DNA sequences associated with a particular strain for its rapid identification by PCR, it should be possible to deploy custom recovery and identification tools for the efficient detection of priority pathogens such as Shiga toxigenic Escherichia coli (STEC) outbreak strains within the time frame of an active investigation. Using a laboratory STEC strain as a model, trimethoprim resistance identified by both Antimicrobial Resistance Marker Identifier and ResFinder was used as the basis for its selective recovery against a background of commensal E. coli bacteria in ground beef samples. Enrichment in modified tryptic soy broth containing trimethoprim greatly enhanced the recovery of low numbers of model strain cells inoculated in ground beef samples, as verified by the enumeration of colonies on plating media using a strain-specific PCR method to determine the recovery efficiency for the target strain. We discuss the relative merits of different AMR marker prediction tools for this purpose and describe how such tools can be utilized to good effect in a typical outbreak investigation scenario.