Menu
September 22, 2019

PBHoover and CigarRoller: a method for confident haploid variant calling on Pacific Biosciences data and its application to heterogeneous population analysis

Motivation: Single Molecule Real-Time (SMRT) sequencing has important and underutilized advantages that amplification-based platforms lack. Lack of systematic error (e.g. GC-bias), complete de novo assembly (including large repetitive regions) without scaffolding, can be mentioned. SMRT sequencing, however suffers from high random error rate and low sequencing depth (older chemistries). Here, we introduce PBHoover, software that uses a heuristic calling algorithm in order to make base calls with high certainty in low coverage regions. This software is also capable of mixed population detection with high sensitivity. PBHoovertextquoterights CigarRoller attachment improves sequencing depth in low-coverage regions through CIGAR-string correction. Results: We tested both modules on 348 M.tuberculosis clinical isolates sequenced on C1 or C2 chemistries. On average, CigarRoller improved percentage of usable read count from 68.9% to 99.98% in C1 runs and from 50% to 99% in C2 runs. Using the greater depth provided by CigarRoller, PBHoover was able to make base and variant calls 99.95% concordant with Sanger calls (QV33). PBHoover also detected antibiotic-resistant subpopulations that went undetected by Sanger. Using C1 chemistry, subpopulations as small as 9% of the total colony can be detected by PBHoover. This provides the most sensitive amplification-free molecular method for heterogeneity analysis and is in line with phenotypic methodstextquoteright sensitivity. This sensitivity significantly improves with the greater depth and lower error rate of the newer chemistries. Availability and Implementation: Executables are freely available under GNU GPL v3+ at http://www.gitlab.com/LPCDRP/pbhoover and http://www.gitlab.com/LPCDRP/CigarRoller. PBHoover is also available on bioconda: https://anaconda.org/bioconda/pbhoover.


September 22, 2019

Potential survival and pathogenesis of a novel strain, Vibrio parahaemolyticus FORC_022, isolated from a soy sauce marinated crab by genome and transcriptome analyses.

Vibrio parahaemolyticus can cause gastrointestinal illness through consumption of seafood. Despite frequent food-borne outbreaks of V. parahaemolyticus, only 19 strains have subjected to complete whole-genome analysis. In this study, a novel strain of V. parahaemolyticus, designated FORC_022 (Food-borne pathogen Omics Research Center_022), was isolated from soy sauce marinated crabs, and its genome and transcriptome were analyzed to elucidate the pathogenic mechanisms. FORC_022 did not include major virulence factors of thermostable direct hemolysin (tdh) and TDH-related hemolysin (trh). However, FORC_022 showed high cytotoxicity and had several V. parahaemolyticus islands (VPaIs) and other virulence factors, such as various secretion systems (types I, II, III, IV, and VI), in comparative genome analysis with CDC_K4557 (the most similar strain) and RIMD2210633 (genome island marker strain). FORC_022 harbored additional virulence genes, including accessory cholera enterotoxin, zona occludens toxin, and tight adhesion (tad) locus, compared with CDC_K4557. In addition, O3 serotype specific gene and the marker gene of pandemic O3:K6 serotype (toxRS) were detected in FORC_022. The expressions levels of genes involved in adherence and carbohydrate transporter were high, whereas those of genes involved in motility, arginine biosynthesis, and proline metabolism were low after exposure to crabs. Moreover, the virulence factors of the type III secretion system, tad locus, and thermolabile hemolysin were overexpressed. Therefore, the risk of foodborne-illness may be high following consumption of FORC_022 contaminated crab. These results provided molecular information regarding the survival and pathogenesis of V. parahaemolyticus FORC_022 strain in contaminated crab and may have applications in food safety.


September 22, 2019

Comparative genomics of Escherichia coli sequence type 219 clones from the same patient: Evolution of the IncI1 blaCMY-carrying plasmid in vivo.

This study investigates the evolution of an Escherichia coli sequence type 219 clone in a patient with recurrent urinary tract infection, comparing isolate EC974 obtained prior to antibiotic treatment and isolate EC1515 recovered after exposure to several ß-lactam antibiotics (ceftriaxone, cefixime, and imipenem). EC974 had a smooth colony morphology, while EC1515 had a rough colony morphology on sheep blood agar. RAPD-PCR analysis suggested that both isolates belonged to the same clone. Antimicrobial susceptibility tests showed that EC1515 was more resistant to piperacillin/tazobactam, cefepime, cefpirome, and ertapenem than EC974. Comparative genomic analysis was used to investigate the genetic changes of EC974 and EC1515 within the host, and showed three plasmids with replicons IncI1, P0111, and IncFII in both isolates. P0111-type plasmids pEC974-2 and pEC1515-2, contained the antibiotic resistance genes aadA2, tetA, and drfA12. IncFII-type plasmids pEC974-3 and pEC1515-3 contained the antibiotic resistance genes blaTEM-1, aadA1, aadA22, sul3, and inuF. Interestingly, blaCMY-111 and blaCMY-4 were found in very similar IncI1 plasmids that also contained aadA22 and aac(3)-IId, from isolates EC974 (pEC974-1) and EC1515 (pEC1515-1), respectively. The results showed in vivo amino acid substitutions converting blaCMY-111 to blaCMY-4 (R221W and A238V substitutions). Conjugation experiments showed a high frequency of IncI1 and IncFII plasmid co-transference. Transconjugants and DH5a cells harboring blaCMY-4 or blaCMY-111 showed higher levels of resistance to ampicillin, amoxicillin, cefazolin, cefuroxime, cefotaxime, cefixime, and ceftazidime, but not piperacillin/tazobactam, cefpime, or ertapenem. All known genes (outer membrane proteins and extended-spectrum AmpC ß-lactamases) involved in ETP resistance in E. coli were identical between EC974 and EC1515. This is the first study to identify the evolution of an IncI1 plasmid within the host, and to characterize blaCMY-111 in E. coli.


September 22, 2019

The draft genomes of Elizabethkingia anophelis of equine origin are genetically similar to three isolates from human clinical specimens.

We report the isolation and characterization of two Elizabethkingia anophelis strains (OSUVM-1 and OSUVM-2) isolated from sources associated with horses in Oklahoma. Both strains appeared susceptible to fluoroquinolones and demonstrated high MICs to all cell wall active antimicrobials including vancomycin, along with aminoglycosides, fusidic acid, chloramphenicol, and tetracycline. Typical of the Elizabethkingia, both draft genomes contained multiple copies of ß-lactamase genes as well as genes predicted to function in antimicrobial efflux. Phylogenetic analysis of the draft genomes revealed that OSUVM-1 and OSUVM-2 differ by only 6 SNPs and are in a clade with 3 strains of Elizabethkingia anophelis that were responsible for human infections. These findings therefore raise the possibility that Elizabethkingia might have the potential to move between humans and animals in a manner similar to known zoonotic pathogens.


September 22, 2019

Integrating long-range connectivity information into de Bruijn graphs.

The de Bruijn graph is a simple and efficient data structure that is used in many areas of sequence analysis including genome assembly, read error correction and variant calling. The data structure has a single parameter k, is straightforward to implement and is tractable for large genomes with high sequencing depth. It also enables representation of multiple samples simultaneously to facilitate comparison. However, unlike the string graph, a de Bruijn graph does not retain long range information that is inherent in the read data. For this reason, applications that rely on de Bruijn graphs can produce sub-optimal results given their input data.We present a novel assembly graph data structure: the Linked de Bruijn Graph (LdBG). Constructed by adding annotations on top of a de Bruijn graph, it stores long range connectivity information through the graph. We show that with error-free data it is possible to losslessly store and recover sequence from a Linked de Bruijn graph. With assembly simulations we demonstrate that the LdBG data structure outperforms both our de Bruijn graph and the String Graph Assembler (SGA). Finally we apply the LdBG to Klebsiella pneumoniae short read data to make large (12 kbp) variant calls, which we validate using PacBio sequencing data, and to characterize the genomic context of drug-resistance genes.Linked de Bruijn Graphs and associated algorithms are implemented as part of McCortex, which is available under the MIT license at https://github.com/mcveanlab/mccortex.Supplementary data are available at Bioinformatics online.


September 22, 2019

Sequencing of pT5282-CTXM, p13190-KPC and p30860-NR, and comparative genomics analysis of IncX8 plasmids.

This study proposes a replicon-based scheme for typing IncX plasmids into nine separately clustering subgroups, including IncX1a, IncX1ß and IncX2-8. The complete nucleotide sequences of three IncX8 plasmids, namely pT5282-CTXM and p30860-NR from Enterobacter cloacae and p13190-KPC from Klebsiella pneumoniae, were determined and were compared with two other previously sequenced IncX8 plasmids (pCAV1043-58 and pCAV1741-16). These five plasmids possessed conserved IncX8 backbones with limited genetic variation with respect to gene content and organisation, and each of them carried one or three accessory modules that harboured resistance markers and metabolic gene clusters as well as transposons, insertion sequence (IS)-based transposition units and miniature inverted repeat transposable elements (MITEs), indicating that the relatively small IncX8 backbones were able to integrate various foreign genetic contents. The resistance genes blaCTX-M-3 and blaTEM-1 (ß-lactam resistance), blaKPC-2 (carbapenem resistance) and ?blaTEM-1, and tet(A) (tetracycline resistance) and mph(E) (macrolide resistance) were found in pT5282-CTXM, p13190-KPC and pCAV1741-16, respectively, whilst p30860-NR and pCAV1043-58 carried no resistance genes. The data presented here provide an insight into the diversification and evolution history of IncX8 plasmids. Copyright © 2018 Elsevier B.V. and International Society of Chemotherapy. All rights reserved.


September 22, 2019

Identification and characterization of conjugative plasmids that encode ciprofloxacin resistance in Salmonella.

This study aimed to characterize novel conjugative plasmids that encode transferrable ciprofloxacin resistance in Salmonella In this study, 157 non-duplicated Salmonella isolates were recovered from food products, 55 out of which were found to be resistant to ciprofloxacin. Interestingly, 37 out of the 55 (67%) CipRSalmonella isolates did not harbor any mutations in the Quinolone resistance determine regions (QRDR). Interestingly, six Salmonella isolates were shown to carry two novel types of conjugative plasmids that could transfer ciprofloxacin resistance phenotype to E. coli J53 (AziR). The first type belonged to the ~110kb IncFIB type conjugative plasmid carrying qnrB-bearing and aac(6′)-Ib-cr-bearing mobile elements. Transfer of the plasmid between E. coli or Salmonella could confer CIP MIC to 1 to 2µg/ml. The second type of conjugative plasmid belonged to ~240kb IncH1/IncF plasmids carrying a single PMQR gene, qnrS Importantly, this type of conjugative ciprofloxacin resistance plasmids could be detected in clinical isolates of Salmonella Dissemination of these conjugative plasmids that confer ciprofloxaicn resistance poses serious public health impact and Salmonella infection control. Copyright © 2018 American Society for Microbiology.


September 22, 2019

Emergence of an XDR and carbapenemase-producing hypervirulent Klebsiella pneumoniae strain in Taiwan.

Carbapenemase-producing Klebsiella pneumoniae causes high mortality owing to the limited therapeutic options available. Here, we investigated an emergent carbapenem-resistant K. pneumoniae strain with hypervirulence found among KPC-2-producing strains in Taiwan.KPC-producing K. pneumoniae strains were collected consecutively from clinical specimens at the Taipei Veterans General Hospital between January 2012 and December 2014. Capsular types and the presence of rmpA/rmpA2 were analysed, and PFGE and MLST performed using these strains. The strain positive for rmpA/rmpA2 was tested in an in vivo mouse lethality study to verify its virulence and subjected to WGS to delineate its genomic features.A total of 62 KPC-2-producing K. pneumoniae strains were identified; all of these belonged to ST11 and capsular genotype K47. One strain isolated from a fatal case with intra-abdominal abscess (TVGHCRE225) harboured rmpA and rmpA2 genes. This strain was resistant to tigecycline and colistin, in addition to carbapenems, and did not belong to the major cluster in PFGE. TVGHCRE225 exhibited high in vivo virulence in the mouse lethality experiment. WGS showed that TVGHCRE225 acquired a novel hybrid virulence plasmid harbouring a set of virulence genes (iroBCDN, iucABCD, rmpA and rmpA2, and iutA) compared with the classic ST11 KPC-2-producing strain.We identified an XDR ST11 KPC-2-producing K. pneumoniae strain carrying a hybrid virulent plasmid in Taiwan. Active surveillance focusing on carbapenem-resistant hypervirulent K. pneumoniae strains is necessary, as the threat to human health is imminent.


September 22, 2019

Chromosomally encoded mcr-5 in colistin non-susceptible Pseudomonas aeruginosa.

Whole genome sequencing (WGS) of historical Pseudomonas aeruginosa clinical isolates identified a chromosomal copy of mcr-5 within a Tn3-like transposon in P. aeruginosa MRSN 12280. The isolate was non-susceptible to colistin by broth microdilution and genome analysis revealed no mutations known to confer colistin resistance. To the best of our knowledge, this is the first report of mcr in colistin non-susceptible P. aeruginosa.


September 22, 2019

Characterization of LE3 and LE4, the only lytic phages known to infect the spirochete Leptospira.

Leptospira is a phylogenetically unique group of bacteria, and includes the causative agents of leptospirosis, the most globally prevalent zoonosis. Bacteriophages in Leptospira are largely unexplored. To date, a genomic sequence is available for only one temperate leptophage called LE1. Here, we sequenced and analysed the first genomes of the lytic phages LE3 and LE4 that can infect the saprophyte Leptospira biflexa using the lipopolysaccharide O-antigen as receptor. Bioinformatics analysis showed that the 48-kb LE3 and LE4 genomes are similar and contain 62% genes whose function cannot be predicted. Mass spectrometry led to the identification of 21 and 23 phage proteins in LE3 and LE4, respectively. However we did not identify significant similarities with other phage genomes. A search for prophages close to LE4 in the Leptospira genomes allowed for the identification of a related plasmid in L. interrogans and a prophage-like region in the draft genome of a clinical isolate of L. mayottensis. Long-read whole genome sequencing of the L. mayottensis revealed that the genome contained a LE4 phage-like circular plasmid. Further isolation and genomic comparison of leptophages should reveal their role in the genetic evolution of Leptospira.


September 22, 2019

Linking genotype and phenotype in an economically viable propionic acid biosynthesis process

Propionic acid (PA) is used as a food preservative and increasingly, as a precursor for the synthesis of monomers. PA is produced mainly through hydrocarboxylation of ethylene, also known as the `oxo-process’; however, Propionibacterium species are promising biological PA producers natively producing PA as their main fermentation product. However, for fermentation to be competitive, a PA yield of at least 0.6 g/g is required.


September 22, 2019

Citrobacter freundii fitness during bloodstream infection.

Sepsis resulting from microbial colonization of the bloodstream is a serious health concern associated with high mortality rates. The objective of this study was to define the physiologic requirements of Citrobacter freundii in the bloodstream as a model for bacteremia caused by opportunistic Gram-negative pathogens. A genetic screen in a murine host identified 177 genes that contributed significantly to fitness, the majority of which were broadly classified as having metabolic or cellular maintenance functions. Among the pathways examined, the Tat protein secretion system conferred the single largest fitness contribution during competition infections and a putative Tat-secreted protein, SufI, was also identified as a fitness factor. Additional work was focused on identifying relevant metabolic pathways for bacteria in the bloodstream environment. Mutations that eliminated the use of glucose or mannitol as carbon sources in vitro resulted in loss of fitness in the murine model and similar results were obtained upon disruption of the cysteine biosynthetic pathway. Finally, the conservation of identified fitness factors was compared within a cohort of Citrobacter bloodstream isolates and between Citrobacter and Serratia marcescens, the results of which suggest the presence of conserved strategies for bacterial survival and replication in the bloodstream environment.


September 22, 2019

Novel enterobacter lineage as leading cause of nosocomial outbreak involving carbapenemase-producing strains.

We investigated unusual carbapenemase-producing Enterobacter cloacae complex isolates (n = 8) in the novel sequence type (ST) 873, which caused nosocomial infections in 2 hospitals in France. Whole-genome sequence typing showed the 1-year persistence of the epidemic strain, which harbored a blaVIM-4 ST1-IncHI2 plasmid, in 1 health institution and 2 closely related strains harboring blaCTX-M-15 in the other. These isolates formed a new subgroup in the E. hormaechei metacluster, according to their hsp60 sequences and phylogenomic analysis. The average nucleotide identities, specific biochemical properties, and pangenomic and functional investigations of isolates suggested isolates of a novel species that had acquired genes associated with adhesion and mobility. The emergence of this novel Enterobacter phylogenetic lineage within hospitals should be closely monitored because of its ability to persist and spread.


September 22, 2019

The complete methylome of an entomopathogenic bacterium reveals the existence of loci with unmethylated adenines.

DNA methylation can serve to control diverse phenomena in eukaryotes and prokaryotes, including gene regulation leading to cell differentiation. In bacteria, DNA methylomes (i.e., methylation state of each base of the whole genome) have been described for several species, but methylome profile variation during the lifecycle has rarely been studied, and only in a few model organisms. Moreover, major phenotypic changes have been reported in several bacterial strains with a deregulated methyltransferase, but the corresponding methylome has rarely been described. Here we report the first methylome description of an entomopathogenic bacterium, Photorhabdus luminescens. Eight motifs displaying a high rate of methylation (>94%) were identified. The methylome was strikingly stable over course of growth, but also in a subpopulation responsible for a critical step in the bacterium’s lifecycle: successful survival and proliferation in insects. The rare unmethylated GATC motifs were preferentially located in putative promoter regions, and most of them were methylated after Dam methyltransferase overexpression, suggesting that DNA methylation is involved in gene regulation. Our findings bring key insight into bacterial methylomes and encourage further research to decipher the role of loci protected from DNA methylation in gene regulation.


September 22, 2019

Characterization of a novel blaKLUC variant with reduced ß-lactam resistance from an IncA/C group plasmid in a clinical Klebsiella pneumoniae isolate.

Similar to other CTX-M family enzymes, KLUC is a recently identified and emerging determinant of cefotaxime resistance that has been recovered from at least three Enterobacteriaceae species, including Kluyvera cryocrescens, Escherichia coli, and Enterobacter cloacae. Whether this extended-spectrum ß-lactamase (ESBL) has been disseminated among commonly isolated Enterobacteriaceae is worthy of further investigation. In this study, we screened 739 nosocomial Enterobacteriaceae isolates (240 Klebsiella pneumoniae and 499 E. coli strains) and found that one K. pneumoniae and four E. coli isolates harbored the blaKLUC gene. Three blaKLUC determinants isolated from E. coli were entirely identical to a blaKLUC-3 gene previously recovered in the same hospital. PFGE of four blaKLUC-harboring E. coli strains showed that prevalence of these determinants was most likely mediated by horizontal gene transfer but not clonal dissemination. However, the variant isolated from K. pneumoniae belonged to a novel member of the KLUC enzyme group. This newly identified enzyme (KLUC-5) has an amino acid substitution compared with previously identified KLUC-1 (G18S) and KLUC-3 (G240D). Antimicrobial susceptibility tests showed that KLUC-5 significantly reduced resistance activity to almost all the selected antimicrobials compared to previously identified KLUC-3. Site-directed mutagenesis showed that blaKLUC-5-D240G and blaKLUC-5-S18G significantly enhanced the MIC against its best substrate. Conjugation and S1-PFGE indicated that blaKLUC-5 was located on a transferable plasmid, which was further decoded by single-molecule, real-time sequencing. Comparative genome analysis showed that its backbone exhibited genetic homology to the IncA/C incompatibility group plasmids. A transposable element, ISEcp1, was detected 256-bp upstream of the blaKLUC-5 gene; this location was inconsistent with the previously identified blaKLUC-1 but congruent with the variants recovered from E. coli in the same hospital. These data provide evidence of the increasingly emerging KLUC group of ESBLs in China.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.