Sequencing of Enteric Bacteria: Library Preparation Procedure Matters for Accurate Identification and Characterization.

Enzymatic library preparation kits are increasingly used for bacterial whole genome sequencing. While they offer a rapid workflow, the transposases used in the kits are recognized to be somewhat biased. The aim of this study was to optimize and validate a protocol for the Illumina DNA Prep kit (formerly Nextera DNA Flex) for sequencing enteric pathogens and compare its performance against the Nextera XT kit. One hundred forty-three strains of Campylobacter, Escherichia, Listeria, Salmonella, Shigella, and Vibrio were prepared with both methods and sequenced on the Illumina MiSeq using 300 and/or 500 cycle chemistries. Sequences were compared using core genome multilocus sequence typing (cgMLST), 7-gene multilocus sequence typing (MLST), and detection of markers encoding serotype, virulence, and antimicrobial resistance. Sequences for one Escherichia strain were downsampled to determine the minimum coverage required for the analyses. While organism-specific differences were observed, the Prep libraries generated longer average read lengths and less fragmented assemblies compared to the XT libraries. In downstream analysis, the most notable difference between the kits was observed for Escherichia, particularly for the 300 cycle sequences. The O group was not predicted in 32% and 4% of XT sequences when using blast and kmer algorithms, respectively, while the O group was predicted from all Prep sequences regardless of the algorithm. In addition, the ehxA gene was not detected in 6% of XT sequences and 34% were missing one or more of the type III secretion systems and/or plasmid-associated genes, which were detected in the Prep sequences. The coverage downsampling revealed that acceptable assembly quality and allele detection was achieved at 30 × coverage with the Prep libraries, whereas 40-50 × coverage was required for the XT libraries. The better performance of the Prep libraries was attributed to more even coverage, particularly in genome regions low in GC content.

ROCker Models for Reliable Detection and Typing of Short-Read Sequences Carrying ß-Lactamase Genes.

Identification of genes encoding ß-lactamases (BLs) from short-read sequences remains challenging due to the high frequency of shared amino acid functional domains and motifs in proteins encoded by BL genes and related non-BL gene sequences. Divergent BL homologs can be frequently missed during similarity searches, which has important practical consequences for monitoring antibiotic resistance. To address this limitation, we built ROCker models that targeted broad classes (e.g., class A, B, C, and D) and individual families (e.g., TEM) of BLs and challenged them with mock 150-bp- and 250-bp-read data sets of known composition. ROCker identifies most-discriminant bit score thresholds in sliding windows along the sequence of the target protein sequence and hence can account for nondiscriminative domains shared by unrelated proteins. BL ROCker models showed a 0% false-positive rate (FPR), a 0% to 4% false-negative rate (FNR), and an up-to-50-fold-higher F1 score [2 × precision × recall/(precision + recall)] compared to alternative methods, such as similarity searches using BLASTx with various e-value thresholds and BL hidden Markov models, or tools like DeepARG, ShortBRED, and AMRFinder. The ROCker models and the underlying protein sequence reference data sets and phylogenetic trees for read placement are freely available through http://enve-omics.ce.gatech.edu/data/rocker-bla. Application of these BL ROCker models to metagenomics, metatranscriptomics, and high-throughput PCR gene amplicon data should facilitate the reliable detection and quantification of BL variants encoded by environmental or clinical isolates and microbiomes and more accurate assessment of the associated public health risk, compared to the current practice. IMPORTANCE Resistance genes encoding ß-lactamases (BLs) confer resistance to the widely prescribed antibiotic class ß-lactams. Therefore, it is important to assess the prevalence of BL genes in clinical or environmental samples for monitoring the spreading of these genes into pathogens and estimating public health risk. However, detecting BLs in short-read sequence data is technically challenging. Our ROCker model-based bioinformatics approach showcases the reliable detection and typing of BLs in complex data sets and thus contributes toward solving an important problem in antibiotic resistance surveillance. The ROCker models developed substantially expand the toolbox for monitoring antibiotic resistance in clinical or environmental settings.

Comparison of four enzymatic library preparation kits for sequencing Shiga toxin-producing Escherichia coli for surveillance and outbreak detection.

Four enzymatic DNA library preparation kits were compared for sequencing Shiga toxin-producing E. coli. All kits produced high quality sequence data which performed equally well in the downstream analyses for surveillance and outbreak detection. Important differences were noted in the workflow user-friendliness and per sample cost.

Detection of Rickettsia Species, and Coxiella-Like and Francisella-Like Endosymbionts in Amblyomma americanum and Amblyomma maculatum from a Shared Field Site in Georgia, United States of America.

Two abundant species of aggressive ticks commonly feed on humans in Georgia: the Gulf Coast tick (Amblyomma maculatum) and the Lone Star tick (A. americanum). A. maculatum is the primary host of Rickettsia parkeri, "Candidatus Rickettsia andeanae," and a Francisella-like endosymbiont (AmacFLE), whereas A. americanum is the primary host for R. amblyommatis, Ehrlichia chaffeensis, E. ewingii, and a Coxiella-like endosymbiont (AamCLE). Horizontal transmission of R. parkeri from A. maculatum to A. americanum by co-feeding has been described, and R. amblyommatis has been found infrequently in A. maculatum ticks. We assessed the prevalence of these agents and whether exchange of tick-associated bacteria is common between A. maculatum and A. americanum collected from the same field site. Unengorged ticks were collected May-August 2014 in west-central Georgia from a 4.14 acre site by flagging and from humans and canines traversing that site. All DNA samples were screened with quantitative PCR assays for the bacteria found in both ticks, and the species of any Rickettsia detected was identified by species-specific TaqMan assays or sequencing of the rickettsial ompA gene. Only R. amblyommatis (15) and AamCLE (39) were detected in 40 A. americanum, while the 74 A. maculatum only contained R. parkeri (30), "Candidatus Rickettsia andeanae" (3), and AmacFLE (74). Neither tick species had either Ehrlichia species. Consequently, we obtained no evidence for the frequent exchange of these tick-borne agents in a natural setting despite high levels of carriage of each agent and the common observance of infestation of both ticks on both dogs and humans at this site. Based on these data, exchange of these Rickettsia, Coxiella, and Francisella agents between A. maculatum and A. americanum appears to be an infrequent event.

Detection of Rickettsia asembonensis in Fleas (Siphonaptera: Pulicidae, Ceratophyllidae) Collected in Five Counties in Georgia, United States.

We conducted a molecular survey of Rickettsia in fleas collected from opossums, road-killed and live-trapped in peridomestic and rural settings, state parks, and from pet cats and dogs in Georgia, United States during 1992-2014. The cat flea, Ctenocephalides felis (Bouché) was the predominant species collected from cats and among the archival specimens from opossums found in peridomestic settings. Polygenis gwyni (Fox) was more prevalent on opossums and a single cotton rat trapped in sylvatic settings. Trapped animals were infested infrequently with the squirrel flea, Orchopeas howardi (Baker) and C. felis. TaqMan assays targeting the BioB gene of Rickettsia felis and the OmpB gene of Rickettsia typhi were used to test 291 flea DNAs for Rickettsia. A subset of 53 C. felis collected from a cat in 2011 was tested in 18 pools which were all bioB TaqMan positive (34% minimum infection prevalence). Of 238 fleas tested individually, 140 (58.8%, 95% confidence interval [CI]: 52.5-64.9%) DNAs were bioB positive. Detection of bioB was more prevalent in individual C. felis (91%) compared to P. gwyni (13.4%). Twenty-one (7.2%) were ompB TaqMan positive, including 18 C. felis (9.5%) and 3 P. gwyni (3.2%). Most of these fleas were also positive with bioB TaqMan; however, sequencing of gltA amplicons detected only DNA of Rickettsia asembonensis. Furthermore, only the R. asembonensis genotype was identified based on NlaIV restriction analysis of a larger ompB fragment. These findings contribute to understanding the diversity of Rickettsia associated with fleas in Georgia and emphasize the need for development of more specific molecular tools for detection and field research on rickettsial pathogens.

Metagenomic Approaches for Public Health Surveillance of Foodborne Infections: Opportunities and Challenges.

Foodborne disease surveillance in the United States is at a critical point. Clinical and diagnostic laboratories are using culture-independent diagnostic tests (CIDTs) to identify the pathogen causing foodborne illness from patient specimens. CIDTs are molecular tests that allow doctors to rapidly identify the bacteria causing illness within hours. CIDTs, unlike previous gold standard methods such as bacterial culture, do not produce an isolate that can be subtyped as part of the national molecular subtyping network for foodborne disease surveillance, PulseNet. Without subtype information, cases can no longer be linked using molecular data to identify potentially related cases that are part of an outbreak. In this review, we discuss the public health needs for a molecular subtyping approach directly from patient specimen and highlight different approaches, including amplicon and shotgun metagenomic sequencing.

Seroprevalence and risk factors for infection with Bartonella bacilliformis in Loja province, Ecuador.

The seroprevalence and epidemiology of Bartonella bacilliformis infection in the Andean highlands of Ecuador is largely unknown. We conducted a sero-epidemiologic survey of 319 healthy children aged 1-15 years living in six rural, mountain communities in Loja Province, Ecuador. Blood was collected by finger stick onto filter paper and dried, and the eluted sera analyzed for antibodies to B. bacilliformis by rPap31 ELISA. Demographic, entomologic, and household variables were assessed to investigate associated risk factors for antibody seropositivity to B. bacilliformis. Seroprevalence of 28% was found among children in the study communities. Increased risk of seropositivity was associated with the presence of lumber piles near houses. Decreased risk of seropositivity was observed with the presence of animal waste and incremental 100 meter increases in elevation. Although investigation of clinical cases of Carrion's disease was not within the scope of this study, our serology data suggest that infection of children with B. bacilliformis is prevalent in this region of Ecuador and is largely unrecognized and undiagnosed. This study highlights the need to further investigate the prevalence, pathogenesis, epidemiology, and disease impact of this pathogen in Ecuador.

A Comparative Analysis of the Lyve-SET Phylogenomics Pipeline for Genomic Epidemiology of Foodborne Pathogens.

Modern epidemiology of foodborne bacterial pathogens in industrialized countries relies increasingly on whole genome sequencing (WGS) techniques. As opposed to profiling techniques such as pulsed-field gel electrophoresis, WGS requires a variety of computational methods. Since 2013, United States agencies responsible for food safety including the CDC, FDA, and USDA, have been performing whole-genome sequencing (WGS) on all Listeria monocytogenes found in clinical, food, and environmental samples. Each year, more genomes of other foodborne pathogens such as Escherichia coli, Campylobacter jejuni, and Salmonella enterica are being sequenced. Comparing thousands of genomes across an entire species requires a fast method with coarse resolution; however, capturing the fine details of highly related isolates requires a computationally heavy and sophisticated algorithm. Most L. monocytogenes investigations employing WGS depend on being able to identify an outbreak clade whose inter-genomic distances are less than an empirically determined threshold. When the difference between a few single nucleotide polymorphisms (SNPs) can help distinguish between genomes that are likely outbreak-associated and those that are less likely to be associated, we require a fine-resolution method. To achieve this level of resolution, we have developed Lyve-SET, a high-quality SNP pipeline. We evaluated Lyve-SET by retrospectively investigating 12 outbreak data sets along with four other SNP pipelines that have been used in outbreak investigation or similar scenarios. To compare these pipelines, several distance and phylogeny-based comparison methods were applied, which collectively showed that multiple pipelines were able to identify most outbreak clusters and strains. Currently in the US PulseNet system, whole genome multi-locus sequence typing (wgMLST) is the preferred primary method for foodborne WGS cluster detection and outbreak investigation due to its ability to name standardized genomic profiles, its central database, and its ability to be run in a graphical user interface. However, creating a functional wgMLST scheme requires extended up-front development and subject-matter expertise. When a scheme does not exist or when the highest resolution is needed, SNP analysis is used. Using three Listeria outbreak data sets, we demonstrated the concordance between Lyve-SET SNP typing and wgMLST. Availability: Lyve-SET can be found at https://github.com/lskatz/Lyve-SET.

The mitochondrial genome of the lone star tick (Amblyomma americanum).

Amblyomma americanum is an abundant tick in the southeastern, midwestern, and northeastern United States. It is a vector of multiple diseases, but limited genomic resources are available for it. We sequenced the complete mitochondrial genome of a single female A. americanum collected in Georgia using the Illumina platform. The consensus sequence was 14,709 bp long, and the mean coverage across the assembly was >12,000×. All expected tick genomic features were present, including two "Tick-Box" motifs, and in the expected order for the Metastriata. Heteroplasmy rates were low compared to the most closely related tick for which data are available, Amblyomma cajennense. The phylogeny derived from the concatenated protein coding and rRNA genes from the 33 available tick mitochondrial genomes was consistent with those previously proposed for the Acari. This is the first complete mitochondrial sequence for A. americanum, which provides a useful reference for future studies of A. americanum population genetics and tick phylogeny.

Seasonal weather, nutrients, and conspecific presence impacts on the southern house mosquito oviposition dynamics in combined sewage overflows.

Combined sewage overflows have created favorable conditions for the establishment of the southern house mosquito, Culex quinquefasciatus Say (Diptera: Culicidae), larvae in natural creeks that would otherwise be unsuitable for the development of this mosquito species. Here, we show the results from a seminatural experiment carried over the three seasons of mosquito activity (spring, summer, and fall) in Tanyard Creek, Atlanta, GA. In this study we manipulated the amount of nutrients by further enriching combined sewage overflow water, and tracked weather variables, organic nutrient concentration, exposure time to conspecifics, and the number of egg rafts collected in experimental containers. We found season and nutrient enrichment to be the most important variables explaining the differences in egg rafts counts. Further analyses suggest that temperature may also play a role in seasonal oviposition patterns. The results from this study suggest that nutrient enrichment and adequate temperatures are important factors shaping Cx. quinquefasciatus oviposition seasonality in combined sewage overflows.