Genomic evolution of the globally disseminated multidrug-resistant Klebsiella pneumoniae clonal group 147

The rapid emergence of multidrug-resistant Klebsiella pneumoniae is being driven largely by the spread of specific clonal groups (CGs). Of these, CG147 includes 7-gene multilocus sequence typing (MLST) sequence types (STs) ST147, ST273 and ST392. CG147 has caused nosocomial outbreaks across the world, but its global population dynamics remain unknown. Here, we report a pandrug-resistant ST147 clinical isolate from India (strain DJ) and define the evolution and global emergence of CG147. Antimicrobial-susceptibility testing following European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines and genome sequencing (Illumina and Oxford Nanopore Technologies, Unicycler assembly) were performed on strain DJ. Additionally, we collated 217 publicly available CG147 genomes [National Center for Biotechnology Information (NCBI), May 2019]. CG147 evolution was inferred within a temporal phylogenetic framework (beast) based on a recombination-free sequence alignment (Roary/Gubbins). Comparative genomic analyses focused on resistance and virulence genes and other genetic elements (BIGSdb, Kleborate, PlasmidFinder, phaster, ICEfinder and CRISPRCasFinder). Strain DJ had a pandrug-resistance phenotype. Its genome comprised the chromosome, seven plasmids and one linear phage-plasmid. Four carbapenemase genes were detected: bla NDM-5 and two copies of bla OXA-181 in the chromosome, and a second copy of bla NDM-5 on an 84 kb IncFII plasmid. CG147 genomes carried a mean of 13 acquired resistance genes or mutations; 63 % carried a carbapenemase gene and 83 % harboured bla CTX-M. All CG147 genomes presented GyrA and ParC mutations and a common subtype I-E CRISPR-Cas system. ST392 and ST273 emerged in 2005 and 1995, respectively. ST147, the most represented phylogenetic branch, was itself divided into two main clades with distinct capsular loci: KL64 (74 %, DJ included, emerged in 1994 and disseminated worldwide, with carbapenemases varying among world regions) and KL10 (20 %, emerged in 2002, predominantly found in Asian countries, associated with carbapenemases NDM and OXA-48-like). Furthermore, subclades within ST147-KL64 differed at the yersiniabactin locus, OmpK35/K36 mutations, plasmid replicons and prophages. The absence of IncF plasmids in some subclades was associated with a possible activity of a CRISPR-Cas system. K. pneumoniae CG147 comprises pandrug-resistant or extensively resistant isolates, and carries multiple and diverse resistance genes and mobile genetic elements, including chromosomal bla NDM-5. Its emergence is being driven by the spread of several phylogenetic clades marked by their own genomic features and specific temporo–spatial dynamics. These findings highlight the need for precision surveillance strategies to limit the spread of particularly concerning CG147 subsets.

Hackflex: low-cost, high-throughput, Illumina Nextera Flex library construction

We developed a low-cost method for the production of Illumina-compatible sequencing libraries that allows up to 14 times more libraries for high-throughput Illumina sequencing to be generated for the same cost. We call this new method Hackflex. The quality of library preparation was tested by constructing libraries from Escherichia coli MG1655 genomic DNA using either Hackflex, standard Nextera Flex (recently renamed as Illumina DNA Prep) or a variation of standard Nextera Flex in which the bead-linked transposase is diluted prior to use. In order to test the library quality for genomes with a higher and a lower G+C content, library construction methods were also tested on Pseudomonas aeruginosa PAO1 and Staphylococcus aureus ATCC 25923, respectively. We demonstrated that Hackflex can produce high-quality libraries and yields a highly uniform coverage, equivalent to the standard Nextera Flex kit. We show that strongly size-selected libraries produce sufficient yield and complexity to support de novo microbial genome assembly, and that assemblies of the large-insert libraries can be much more contiguous than standard libraries without strong size selection. We introduce a new set of sample barcodes that are distinct from standard Illumina barcodes, enabling Hackflex samples to be multiplexed with samples barcoded using standard Illumina kits. Using Hackflex, we were able to achieve a per-sample reagent cost for library prep of A$7.22 (Australian dollars) (US $5.60; UK £3.87, £1=A$1.87), which is 9.87 times lower than the standard Nextera Flex protocol at advertised retail price. An additional simple modification and further simplification of the protocol by omitting the wash step enables a further price reduction to reach an overall 14-fold cost saving. This method will allow researchers to construct more libraries within a given budget, thereby yielding more data and facilitating research programmes where sequencing large numbers of libraries is beneficial.

Transcriptomic analysis of Pseudomonas ogarae F113 reveals the antagonistic roles of AmrZ and FleQ during rhizosphere adaption

Rhizosphere colonization by bacteria involves molecular and cellular mechanisms, such as motility and chemotaxis, biofilm formation, metabolic versatility, or biosynthesis of secondary metabolites, among others. Nonetheless, there is limited knowledge concerning the main regulatory factors that drive the rhizosphere colonization process. Here we show the importance of the AmrZ and FleQ transcription factors for adaption in the plant growth-promoting rhizobacterium (PGPR) and rhizosphere colonization model Pseudomonas ogarae F113. RNA-Seq analyses of P. ogarae F113 grown in liquid cultures either in exponential and stationary growth phase, and rhizosphere conditions, revealed that rhizosphere is a key driver of global changes in gene expression in this bacterium. Regarding the genetic background, this work has revealed that a mutation in fleQ causes considerably more alterations in the gene expression profile of this bacterium than a mutation in amrZ under rhizosphere conditions. The functional analysis has revealed that in P. ogarae F113, the transcription factors AmrZ and FleQ regulate genes involved in diverse bacterial functions. Notably, in the rhizosphere, these transcription factors antagonistically regulate genes related to motility, biofilm formation, nitrogen, sulfur, and amino acid metabolism, transport, signalling, and secretion, especially the type VI secretion systems. These results define the regulon of two important bifunctional transcriptional regulators in pseudomonads during the process of rhizosphere colonization.

High quality genome assembly of the amitochondriate eukaryote Monocercomonoides exilis

Monocercomonoides exilis is considered the first known eukaryote to completely lack mitochondria. This conclusion is based primarily on a genomic and transcriptomic study which failed to identify any mitochondrial hallmark proteins. However, the available genome assembly has limited contiguity and around 1.5 % of the genome sequence is represented by unknown bases. To improve the contiguity, we re-sequenced the genome and transcriptome of M. exilis using Oxford Nanopore Technology (ONT). The resulting draft genome is assembled in 101 contigs with an N50 value of 1.38 Mbp, almost 20 times higher than the previously published assembly. Using a newly generated ONT transcriptome, we further improve the gene prediction and add high quality untranslated region (UTR) annotations, in which we identify two putative polyadenylation signals present in the 3′UTR regions and characterise the Kozak sequence in the 5′UTR regions. All these improvements are reflected by higher BUSCO genome completeness values. Regardless of an overall more complete genome assembly without missing bases and a better gene prediction, we still failed to identify any mitochondrial hallmark genes, thus further supporting the hypothesis on the absence of mitochondrion.

Genome placement of alpha-haemolysin cluster is associated with alpha-haemolysin sequence variation, adhesin and iron acquisition factor profile of Escherichia coli

Since the discovery of haemolysis, many studies focused on a deeper understanding of this phenotype in Escherichia coli and its association with other virulence genes, diseases and pathogenic attributes/functions in the host. Our virulence-associated factor profiling and genome-wide association analysis of genomes of haemolytic and nonhaemolytic E. coli unveiled high prevalence of adhesins, iron acquisition genes and toxins in haemolytic bacteria. In the case of fimbriae with high prevalence, we analysed sequence variation of FimH, EcpD and CsgA, and showed that different adhesin variants were present in the analysed groups, indicating altered adhesive capabilities of haemolytic and nonhaemolytic E. coli . Analysis of over 1000 haemolytic E. coli genomes revealed that they are pathotypically, genetically and antigenically diverse, but their adhesin and iron acquisition repertoire is associated with genome placement of hlyCABD cluster. Haemolytic E. coli with chromosome-encoded alpha-haemolysin had high frequency of P, S, Auf fimbriae and multiple iron acquisition systems such as aerobactin, yersiniabactin, salmochelin, Fec, Sit, Bfd and hemin uptake systems. Haemolytic E. coli with plasmid-encoded alpha-haemolysin had similar adhesin profile to nonpathogenic E. coli, with high prevalence of Stg, Yra, Ygi, Ycb, Ybg, Ycf, Sfm, F9 fimbriae, Paa, Lda, intimin and type 3 secretion system encoding genes. Analysis of HlyCABD sequence variation revealed presence of variants associated with genome placement and pathotype.

Genome-wide analysis of fitness-factors in uropathogenic Escherichia coli during growth in laboratory media and during urinary tract infections

Uropathogenic Escherichia coli (UPEC) UTI89 is a well-characterized strain, which has mainly been used to study UPEC virulence during urinary tract infection (UTI). However, little is known on UTI89 key fitness-factors during growth in lab media and during UTI. Here, we used a transposon-insertion-sequencing approach (TraDIS) to reveal the UTI89 essential-genes for in vitro growth and fitness-gene-sets for growth in Luria broth (LB) and EZ-MOPS medium without glucose, as well as for human bacteriuria and mouse cystitis. A total of 293 essential genes for growth were identified and the set of fitness-genes was shown to differ depending on the growth media. A modified, previously validated UTI murine model, with administration of glucose prior to infection was applied. Selected fitness-genes for growth in urine and mouse-bladder colonization were validated using deletion-mutants. Novel fitness-genes, such as tusA, corA and rfaG; involved in sulphur-acquisition, magnesium-uptake, and LPS-biosynthesis, were proved to be important during UTI. Moreover, rfaG was confirmed as relevant in both niches, and therefore it may represent a target for novel UTI-treatment/prevention strategies.

Accuracy of an amplicon-sequencing nanopore approach to identify variants in tuberculosis drug-resistance-associated genes

A rapid and accurate diagnostic assay represents an important means to detect Mycobacterium tuberculosis , identify drug-resistant strains and ensure treatment success. Currently employed techniques to diagnose drug-resistant tuberculosis include slow phenotypic tests or more rapid molecular assays that evaluate a limited range of drugs. Whole-genome-sequencing-based approaches can detect known drug-resistance-conferring mutations and novel variations; however, the dependence on growing samples in culture, and the associated delays in achieving results, represents a significant limitation. As an alternative, targeted sequencing strategies can be directly performed on clinical samples at high throughput. This study proposes a targeted sequencing assay to rapidly detect drug-resistant strains of M. tuberculosis using the Nanopore MinION sequencing platform. We designed a single-tube assay that targets nine genes associated with drug resistance to seven drugs and two phylogenetic-determining regions to determine strain lineage and tested it in nine clinical isolates and six sputa. The study’s main aim is to calibrate MinNION variant calling to detect drug-resistance-associated mutations with different frequencies to match the accuracy of Illumina (the current gold-standard sequencing technology) from both culture and sputum samples. After calibrating Nanopore MinION variant calling, we demonstrated 100% agreement between Illumina WGS and our MinION set up to detect known drug resistance and phylogenetic variants in our dataset. Importantly, other variants in the amplicons are also detected, decreasing the recall. We identify minority variants and insertions/deletions as crucial bioinformatics challenges to fully reproduce Illumina WGS results.

Worldwide distribution and environmental origin of the Adelaide imipenemase (AIM-1), a potent carbapenemase in Pseudomonas aeruginosa

Carbapenems are potent broad-spectrum β-lactam antibiotics reserved for the treatment of serious infections caused by multidrug-resistant bacteria such as Pseudomonas aeruginosa . The surge in P. aeruginosa resistant to carbapenems is an urgent threat, as very few treatment options remain. Resistance to carbapenems is predominantly due to the presence of carbapenemase enzymes. The assessment of 147 P . aeruginosa isolates revealed that 32 isolates were carbapenem non-wild-type. These isolates were screened for carbapenem resistance genes using PCR. One isolate from wastewater contained the Adelaide imipenemase gene (bla AIM-1) and was compared phenotypically with a highly carbapenem-resistant clinical isolate containing the bla AIM-1 gene. A further investigation of wastewater samples from various local healthcare and non-healthcare sources as well as river water, using probe-based qPCR, revealed the presence of the bla AIM-1 gene in all the samples analysed. The widespread occurrence of bla AIM-1 throughout Adelaide hinted at the possibility of more generally extensive spread of this gene than originally thought. A blast search revealed the presence of the bla AIM-1 gene in Asia, North America and Europe. To elucidate the identity of the organism(s) carrying the bla AIM-1 gene, shotgun metagenomic sequencing was conducted on three wastewater samples from different locations. Comparison of these nucleotide sequences with a whole-genome sequence of a P. aeruginosa isolate revealed that, unlike the genetic environment and arrangement in P. aeruginosa , the bla AIM-1 gene was not carried as part of any mobile genetic elements. A phylogenetic tree constructed with the deduced amino acid sequences of AIM-1 suggested that the potential origin of the bla AIM-1 gene in P. aeruginosa might be the non-pathogenic environmental organism, Pseudoxanthomonas mexicana .

Context-aware genomic surveillance reveals hidden transmission of a carbapenemase-producing Klebsiella pneumoniae

Genomic surveillance can inform effective public health responses to pathogen outbreaks. However, integration of non-local data is rarely done. We investigate two large hospital outbreaks of a carbapenemase-carrying Klebsiella pneumoniae strain in Germany and show the value of contextual data. By screening about 10 000 genomes, over 400 000 metagenomes and two culture collections using in silico and in vitro methods, we identify a total of 415 closely related genomes reported in 28 studies. We identify the relationship between the two outbreaks through time-dated phylogeny, including their respective origin. One of the outbreaks presents extensive hidden transmission, with descendant isolates only identified in other studies. We then leverage the genome collection from this meta-analysis to identify genes under positive selection. We thereby identify an inner membrane transporter (ynjC) with a putative role in colistin resistance. Contextual data from other sources can thus enhance local genomic surveillance at multiple levels and should be integrated by default when available.

The bacterial biome of ticks and their wildlife hosts at the urban–wildland interface

Advances in sequencing technologies have revealed the complex and diverse microbial communities present in ticks (Ixodida). As obligate blood-feeding arthropods, ticks are responsible for a number of infectious diseases that can affect humans, livestock, domestic animals and wildlife. While cases of human tick-borne diseases continue to increase in the northern hemisphere, there has been relatively little recognition of zoonotic tick-borne pathogens in Australia. Over the past 5 years, studies using high-throughput sequencing technologies have shown that Australian ticks harbour unique and diverse bacterial communities. In the present study, free-ranging wildlife (n=203), representing ten mammal species, were sampled from urban and peri-urban areas in New South Wales (NSW), Queensland (QLD) and Western Australia (WA). Bacterial metabarcoding targeting the 16S rRNA locus was used to characterize the microbiomes of three sample types collected from wildlife: blood, ticks and tissue samples. Further sequence information was obtained for selected taxa of interest. Six tick species were identified from wildlife: Amblyomma triguttatum, Ixodes antechini, Ixodes australiensis, Ixodes holocyclus, Ixodes tasmani and Ixodes trichosuri. Bacterial 16S rRNA metabarcoding was performed on 536 samples and 65 controls, generating over 100 million sequences. Alpha diversity was significantly different between the three sample types, with tissue samples displaying the highest alpha diversity (P<0.001). Proteobacteria was the most abundant taxon identified across all sample types (37.3 %). Beta diversity analysis and ordination revealed little overlap between the three sample types (P<0.001). Taxa of interest included Anaplasmataceae , Bartonella , Borrelia , Coxiellaceae , Francisella , Midichloria , Mycoplasma and Rickettsia . Anaplasmataceae bacteria were detected in 17.7% (95/536) of samples and included Anaplasma , Ehrlichia and Neoehrlichia species. In samples from NSW, ‘Ca. Neoehrlichia australis’, ‘Ca. Neoehrlichia arcana’, Neoehrlichia sp. and Ehrlichia sp. were identified. A putative novel Ehrlichia sp. was identified from WA and Anaplasma platys was identified from QLD. Nine rodent tissue samples were positive for a novel Borrelia sp. that formed a phylogenetically distinct clade separate from the Lyme Borrelia and relapsing fever groups. This novel clade included recently identified rodent-associated Borrelia genotypes, which were described from Spain and North America. Bartonella was identified in 12.9% (69/536) of samples. Over half of these positive samples were obtained from black rats (Rattus rattus), and the dominant bacterial species identified were Bartonella coopersplainsensis and Bartonella queenslandensis . The results from the present study show the value of using unbiased high-throughput sequencing applied to samples collected from wildlife. In addition to understanding the sylvatic cycle of known vector-associated pathogens, surveillance work is important to ensure preparedness for potential zoonotic spillover events.