SeekDeep: single-base resolution de novo clustering for amplicon deep sequencing

PCR amplicon deep sequencing continues to transform the investigation of genetic diversity in viral, bacterial, and eukaryotic populations. In eukaryotic populations such as Plasmodium falciparum infections, it is important to discriminate sequences differing by a single nucleotide polymorphism. In bacterial populations, single-base resolution can provide improved resolution towards species and strains. Here, we introduce the SeekDeep suite built around the qluster algorithm, which is capable of accurately building de novo clusters representing true, biological local haplotypes differing by just a single base. It outperforms current software, particularly at low frequencies and at low input read depths, whether resolving single-base differences or traditional OTUs. SeekDeep is open source and works with all major sequencing technologies, making it broadly useful in a wide variety of applications of amplicon deep sequencing to extract accurate and maximal biologic information

Next-Generation Sequencing and Comparative Analysis of Sequential Outbreaks Caused by Multidrug-Resistant Acinetobacter baumannii at a Large Academic Burn Center

Next-generation sequencing (NGS) analysis has emerged as a promising molecular epidemiological method for investigating health care-associated outbreaks. Here, we used NGS to investigate a 3-year outbreak of multidrug-resistant Acinetobacter baumannii (MDRAB) at a large academic burn center. A reference genome from the index case was generated using de novo assembly of PacBio reads. Forty-six MDRAB isolates were analyzed by pulsed-field gel electrophoresis (PFGE) and sequenced using an Illumina platform. After mapping to the index case reference genome, four samples were excluded due to low coverage, leaving 42 samples for further analysis. Multilocus sequence types (MLST) and the presence of acquired resistance genes were also determined from the sequencing data. A transmission network was inferred from genomic and epidemiological data using a Bayesian framework. Based on single-nucleotide variant (SNV) differences, this MDRAB outbreak represented three sequential outbreaks caused by distinct clones. The first and second outbreaks were caused by sequence type 2 (ST2), while the third outbreak was caused by ST79. For the second outbreak, the MLST and PFGE results were discordant. However, NGS-based SNV typing detected a recombination event and consequently enabled a more accurate phylogenetic analysis. The distribution of resistance genes varied among the three outbreaks. The first- and second-outbreak strains possessed a bla OXA-23-like group, while the third-outbreak strains harbored a bla OXA-40-like group. NGS-based analysis demonstrated the superior resolution of outbreak transmission networks for MDRAB and provided insight into the mechanisms of strain diversification between sequential outbreaks through recombination

A prolonged outbreak of KPC-3-producing Enterobacter cloacae and Klebsiella pneumoniae driven by multiple mechanisms of resistance transmission at a large academic burn center

Klebsiella pneumoniae carbapenemase (KPC)-producing Enterobacter cloacae have been recently recognized in the United States. Whole-genome sequencing (WGS) has become a useful tool for analysis of outbreaks and for determining transmission networks of multidrug-resistant organisms in healthcare settings, including carbapenem-resistant Enterobacteriaceae (CRE). We experienced a prolonged outbreak of CRE of E. cloacae and K. pneumoniae over a three-year period at a large academic burn center despite rigorous infection control measures. To understand the molecular mechanisms that sustained this outbreak, we investigated the CRE outbreak isolates using WGS. Twenty-two clinical isolates of CRE, including E. cloacae (N=15) and K. pneumoniae (N=7), were sequenced and analyzed genetically. WGS revealed that this outbreak, which seemed epidemiologically unlinked, was in fact genetically linked over a prolonged period. Multiple mechanisms were found to account for the ongoing outbreak of KPC-3-producing E. cloacae and K. pneumoniae . This outbreak was primarily maintained by a clonal expansion of E. cloacae ST114 with distribution of multiple resistance determinants. Plasmid and transposon analysis suggested that the majority of bla KPC-3 was transmitted via an identical Tn 4401 b element on part of a common plasmid. WGS analysis demonstrated complex transmission dynamics within the burn center at levels of strain and/or plasmid in association with transposon, highlighting the versatility of KPC-producing Enterobacteriaceae in their ability to utilize multiple modes to resistance-gene propagation

Genomic Analysis of Multidrug-Resistant Escherichia coli from North Carolina Community Hospitals: Ongoing Circulation of CTX-M-Producing ST131- H 30Rx and ST131- H 30R1 Strains

ABSTRACT Escherichia coli sequence type 131 (ST131) predominates globally among multidrug-resistant (MDR) E. coli strains. We used whole-genome sequencing (WGS) to investigate 63 MDR E. coli isolates from 7 North Carolina community hospitals (2010 to 2015). Of these, 39 (62%) represented ST131, including 37 (95%) from the ST131- H 30R subclone: 10 (27%) from its H 30R1 subset and 27 (69%) from its H 30Rx subset. ST131 core genomes differed by a median of 15 (range, 0 to 490) single-nucleotide variants (SNVs) overall versus only 7 within H 30R1 (range, 3 to 12 SNVs) and 11 within H 30Rx (range, 0 to 21). The four isolates with identical core genomes were all H 30Rx. Epidemiological and clinical characteristics did not vary significantly by strain type, but many patients with MDR E. coli or H 30Rx infection were critically ill and had poor outcomes. H 30Rx isolates characteristically exhibited fluoroquinolone resistance and CTX-M-15 production, had a high prevalence of trimethoprim-sulfamethoxazole resistance (89%), sul1 (89%), and dfrA17 (85%), and were enriched for specific virulence traits, and all qualified as extraintestinal pathogenic E. coli . The high overall prevalence of CTX-M-15 appeared to be possibly attributable to its association with the ST131- H 30Rx subclone and IncF[F2:A1:B−] plasmids. Some phylogenetically clustered non-ST131 MDR E. coli isolates also had distinctive serotypes/ fimH types, fluoroquinolone mutations, CTX-M variants, and IncF types. Thus, WGS analysis of our community hospital source MDR E. coli isolates suggested ongoing circulation and differentiation of E. coli ST131 subclones, with clonal segregation of CTX-M variants, other resistance genes, Inc-type plasmids, and virulence genes

Longitudinal Pooled Deep Sequencing of the Plasmodium vivax K12 Kelch Gene in Cambodia Reveals a Lack of Selection by Artemisinin

The emergence of artemisinin resistance among Plasmodium falciparum in the Greater Mekong subregion threatens malaria control interventions and is associated with multiple unique mutations in K13 (PF3D7_1343700). The aim of this study was to survey Cambodian Plasmodium vivax for mutations in the K13 ortholog (K12, PVX_083080) that might similarly confer artemisinin resistance. Extracted DNA from Cambodian isolates collected between 2009 and 2012 was pooled by province and year and submitted for next-generation sequencing. Single-nucleotide polymorphisms (SNPs) were identified using a pile-up approach that detected minority SNPs. Among the 14 pools, we found six unique SNPs, including three nonsynonymous SNPs, across six codons in K12. However, none of the SNPs were orthologous to artemisinin resistance–conferring mutations in PF3D7_1343700, and nonsynonymous changes did not persist through time within populations. These results suggest a lack of selection in the P. vivax population in Cambodia due to artemisinin drug pressure

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia

In Cambodia, where Plasmodium vivax and Plasmodium falciparum are coendemic and intense multimodal malaria-control interventions have reduced malaria incidence, P. vivax malaria has proven relatively resistant to such measures. We performed comparative genomic analyses of 150 P. vivax and P. falciparum isolates to determine whether different evolutionary strategies might underlie this species-specific resilience. Demographic modeling and tests of selection show that, in contrast to P. falciparum, P. vivax has experienced uninterrupted growth and positive selection at multiple loci encoding transcriptional regulators. In particular, a strong selective sweep involving an AP2 transcription factor suggests that P. vivax may use nuanced transcriptional approaches to population maintenance. Better understanding of P. vivax transcriptional regulation may lead to improved tools to achieve elimination

Partner-Drug Resistance and Population Substructuring of Artemisinin-Resistant Plasmodium falciparum in Cambodia

Plasmodium falciparum in western Cambodia has developed resistance to artemisinin and its partner drugs, causing frequent treatment failure. Understanding this evolution can inform the deployment of new therapies. We investigated the genetic architecture of 78 falciparum isolates using whole-genome sequencing, correlating results to in vivo and ex vivo drug resistance and exploring the relationship between population structure, demographic history, and partner drug resistance. Principle component analysis, network analysis and demographic inference identified a diverse central population with three clusters of clonally expanding parasite populations, each associated with specific K13 artemisinin resistance alleles and partner drug resistance profiles which were consistent with the sequential deployment of artemisinin combination therapies in the region. One cluster displayed ex vivo piperaquine resistance and mefloquine sensitivity with a high rate of in vivo failure of dihydroartemisinin-piperaquine. Another cluster displayed ex vivo mefloquine resistance and piperaquine sensitivity with high in vivo efficacy of dihydroartemisinin-piperaquine. The final cluster was clonal and displayed intermediate sensitivity to both drugs. Variations in recently described piperaquine resistance markers did not explain the difference in mean IC90 or clinical failures between the high and intermediate piperaquine resistance groups, suggesting additional loci may be involved in resistance. The results highlight an important role for partner drug resistance in shaping the P. falciparum genetic landscape in Southeast Asia and suggest that further work is needed to evaluate for other mutations that drive piperaquine resistance

Pooled Deep Sequencing of Plasmodium falciparum Isolates: An Efficient and Scalable Tool to Quantify Prevailing Malaria Drug-Resistance Genotypes

Molecular surveillance for drug-resistant malaria parasites requires reliable, timely, and scalable methods. These data may be efficiently produced by genotyping parasite populations using second-generation sequencing (SGS). We designed and validated a SGS protocol to quantify mutant allele frequencies in the Plasmodium falciparum genes dhfr and dhps in mixed isolates. We applied this new protocol to field isolates from children and compared it to standard genotyping using Sanger sequencing. The SGS protocol accurately quantified dhfr and dhps allele frequencies in a mixture of parasite strains. Using SGS of DNA that was extracted and then pooled from individual isolates, we estimated mutant allele frequencies that were closely correlated to those estimated by Sanger sequencing (correlations, >0.98). The SGS protocol obviated most molecular steps in conventional methods and is cost saving for parasite populations >50. This SGS genotyping method efficiently and reproducibly estimates parasite allele frequencies within populations of P. falciparum for molecular epidemiologic studies

Pooled Amplicon Deep Sequencing of Candidate Plasmodium falciparum Transmission-Blocking Vaccine Antigens

Polymorphisms within Plasmodium falciparum vaccine candidate antigens have the potential to compromise vaccine efficacy. Understanding the allele frequencies of polymorphisms in critical binding regions of antigens can help in the designing of strain-transcendent vaccines. Here, we adopt a pooled deep-sequencing approach, originally designed to study P. falciparum drug resistance mutations, to study the diversity of two leading transmission-blocking vaccine candidates, Pfs25 and Pfs48/45. We sequenced 329 P. falciparum field isolates from six different geographic regions. Pfs25 showed little diversity, with only one known polymorphism identified in the region associated with binding of transmission-blocking antibodies among our isolates. However, we identified four new mutations among eight non-synonymous mutations within the presumed antibody-binding region of Pfs48/45. Pooled deep sequencing provides a scalable and cost-effective approach for the targeted study of allele frequencies of P. falciparum candidate vaccine antigens

Variation in the Microbiota of Ixodes Ticks with Regard to Geography, Species, and Sex

Ixodes scapularis is the principal vector of Lyme disease on the East Coast and in the upper Midwest regions of the United States, yet the tick is also present in the Southeast, where Lyme disease is absent or rare. A closely related species, I. affinis, also carries the pathogen in the South but does not seem to transmit it to humans. In order to better understand the geographic diversity of the tick, we analyzed the microbiota of 104 adult I. scapularis and 13 adult I. affinis ticks captured in 19 locations in South Carolina, North Carolina, Virginia, Connecticut, and New York. Initially, ticks from 4 sites were analyzed by 454 pyrosequencing. Subsequently, ticks from these sites plus 15 others were analyzed by sequencing with an Illumina MiSeq machine. By both analyses, the microbiomes of female ticks were significantly less diverse than those of male ticks. The dissimilarity between tick microbiomes increased with distance between sites, and the state in which a tick was collected could be inferred from its microbiota. The genus Rickettsia was prominent in all locations. Borrelia was also present in most locations and was present at especially high levels in one site in western Virginia. In contrast, members of the family Enterobacteriaceae were very common in North Carolina I. scapularis ticks but uncommon in I. scapularis ticks from other sites and in North Carolina I. affinis ticks. These data suggest substantial variations in the Ixodes microbiota in association with geography, species, and sex