Skip to main content
Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Microbiol Resour Announc. 2019 Nov; 8(48): e01269-19.
Published online 2019 Nov 27. doi: 10.1128/MRA.01269-19
PMCID: PMC6883113
PMID: 31776226

Complete Circular Genome Sequence of a Multidrug-Resistant Escherichia coli Strain from Cuba Obtained with Nanopore and Illumina Hybrid Assembly

David Rasko, Editor
David Rasko, University of Maryland School of Medicine;

The complete genome sequence of a multidrug-resistant Escherichia coli strain isolated from a healthy pig in Cuba was determined using short and long reads. This strain carried four plasmids, including a 42,683-kb IncX1 plasmid, which contains the third-generation cephalosporin resistance gene blaCTX-M-32 together with other disinfectant and antibiotic resistance genes.

ABSTRACT

The complete genome sequence of a multidrug-resistant Escherichia coli strain isolated from a healthy pig in Cuba was determined using short and long reads. This strain carried four plasmids, including a 42,683-kb IncX1 plasmid, which contains the third-generation cephalosporin resistance gene blaCTX-M-32 together with other disinfectant and antibiotic resistance genes.

ANNOUNCEMENT

Commensal Escherichia coli from animals represents a reservoir of acquired antibiotic resistance elements, which can be transferred to pathogenic E. coli strains (1). The multidrug-resistant (MDR) E. coli strain described here is part of a larger study screening pigs for the presence of third-generation cephalosporin-resistant E. coli in Cuba, where blaCTX-M-32 was found to be a predominant third-generation cephalosporin resistance gene. We completely sequenced one of these CTX-M-32-containing E. coli strains with both short and long reads to gain information on the structure of the mobile genetic elements and localization of the antibiotic resistance genes.

A rectal swab from a healthy pig was collected in a swine facility of the Matanzas municipality in July 2016 and was grown overnight at 37°C on a MacConkey (BioCen, Cuba) agar plate supplemented with cefotaxime (4 μg/ml) (Sigma-Aldrich, USA). Identification of the E. coli (strain PK6) isolate was performed by matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) mass spectrometry (Microflex LT, Germany). Genomic DNA was extracted from an overnight culture in agitated LB broth at 37°C using the DNeasy blood and tissue extraction kit (Qiagen, Germany) and purified using the AMPure XP kit (Beckman Coulter, USA). DNA was quantified using a Qubit 3.0 fluorometer (Invitrogen, USA), and 1.5 μg of DNA was used for sequencing.

Whole-genome sequencing was performed with both long and short reads to obtain a complete scaffold of the genome and an accurate sequence. Short reads were obtained using the NovaSeq 6000 S2 reagent kit (300 cycles) and an S2 flow cell on a NovaSeq 6000 (2 × 150-bp paired-end reads) system (Illumina, USA) at Eurofins Genomics GmbH, Germany, yielding 6,411,400 reads (1.9 Gbp of data and 385× coverage). Default parameters were used for all the following bioinformatic software. In order to remove adaptors and low-quality reads (quality value [QV], ≤20), short reads were trimmed using Trimmomatic v0.36 (illuminaclip:TruSeq3-PE.fa:2:30:10, LEADING:3, TRAILING:3, SLIDINGWINDOW:4:15, and MINLEN:36) (2). For long-read sequencing, the genomic DNA was sheared using Covaris g-TUBES to generate fragments of around 20 kb, and the resulting total DNA was used. The library was prepared using the ligation sequencing kit 1D SQK-LSK108 and native barcoding kit 1D EXP-NBD103 in an R9.4 SpotON flow cell with a MinION Mk1B device from Oxford Nanopore Technologies (ONT; United Kingdom), yielding 269,705 reads (1.5 Gbp of data, 309× coverage, mean read length of 6,664 bp, and N50 of 11,844 bp). The reads were base called and demultiplexed using Guppy basecaller (v2.3.7) and Guppy barcoder (v2.3.7), respectively (ONT). Hybrid de novo assembly of both short and long reads and circularization of the replicons were performed using the Unicycler v0.4.4 pipeline (3). The genome of E. coli PK6 was annotated using the NCBI Prokaryotic Genome Annotation Pipeline (4). The complete genome consisted of five circular contigs that are described in Table 1, with a sum of 4,789,128 bp and a G+C content of 50.8%, 4,582 coding sequences, 89 tRNAs, 22 rRNAs, 1 transfer-messenger RNA (tmRNA), and 10 noncoding RNAs (ncRNAs).

TABLE 1

Characteristics of the genome of E. coli strain PK6 from Cuba

Genetic elementSequence nameAssembly size (bp)G+C content (%)Multiplicity (k)a Incompatibility group (plasmids)b Antimicrobial resistance mechanismsc
ChromosomePK6CUB-RH_CHR4,726,52650.81.00NAQnrB19, TEM-1B, Tet(B), GyrA (S83L), ParC (S80I)
Plasmid 1pRHEcCUB-142,68348.51.30IncX1AadA1, AadA2, CmlA1, DfrA12, QacH2, Sul3 (In640); CTX-M-32
Plasmid 2pRHEcCUB-212,60245.52.94NI
Plasmid 3pRHEcCUB-34,59350.413.53NI
Plasmid 4pRHEcCUB-42,72447.74.14NI
aMultiplicity was calculated using Unicylcler to distinguish between the single-copy contig (k = 1) and the multiple contigs (k > 1) (3).
bNA, not applicable; NI, not identified using PlasmidFinder (v2.0) (6).
cAntimicrobial resistance mechanisms were AadA1 and AadA2, streptomycin/spectinomycin adenylyltransferases; CmlA1, chloramphenicol efflux protein; CTX-M-32, extended-spectrum β-lactamase for cephalosporin resistance; DfrA12, trimethoprim-resistant dihydrofolate reductase; GyrA (S83L) and ParC(S80I), serine-to-leucine substitution at position 83 in topoisomerase II (GyrA) and serine-to-isoleucine substitution at position 80 in topoisomerase IV (ParC) for high-level resistance to fluoroquinolones; QacH2, quaternary ammonium compound exporter protein; QnrB19, DNA gyrase protection protein for low-level resistance to fluoroquinolones; Sul3, sulfonamide-resistant dihydropteroate synthase; TEM-1B, β-lactamase; Tet(B), tetracycline efflux protein; and In640, integron containing the qacH2, aadA1, aadA2, cmlA1, and dfrA12 genes as determined using the INTEGRALL database.

Strain PK6 belongs to sequence type 1695 (ST1695), as determined by in silico analysis using MLST (v2.0) based on the Achtman scheme for E. coli (5). Only one plasmid, namely, pRHEcCUB-1, could be assigned to an incompatibility group (IncX1) using PlasmidFinder (v2.0) (6).

Antibiotic resistance genes were detected in silico with ResFinder (v3.1) (7) on the chromosome as well as on the 42,683-bp IncX1 plasmid pRHEcCUB-1. Plasmid pRHEcCUB-1 also contains a gene for resistance to a quaternary ammonium compound disinfectant (Table 1).

This report provides information on the complete and circularized genome of an MDR and extended-spectrum β-lactamase-producing E. coli isolate from Cuba. The sequence can serve as a baseline for future molecular epidemiological studies and for surveillance of antibiotic resistance in bacteria from humans and animals in Cuba.

ACKNOWLEDGMENTS

This study was financed by internal funds of the Institute of Veterinary Bacteriology of Bern and by the National Centre of Animal and Plant Health of Cuba. Rosa Elena Hernández-Fillor received a grant (number 2018.0714) from the Federal Commission for Scholarships for Foreign Students (FCS) of Switzerland.

REFERENCES

1. European Food Safety Authority and European Centre for Disease Prevention and Control. 2019. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2017. EFSA J 17:5598. doi: 10.2903/j.efsa.2019.5598. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
2. Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. doi: 10.1093/bioinformatics/btu170. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
3. Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol 13:e1005595. doi: 10.1371/journal.pcbi.1005595. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
4. Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J. 2016. NCBI Prokaryotic Genome Annotation Pipeline. Nucleic Acids Res 44:6614–6624. doi: 10.1093/nar/gkw569. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
5. Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-Pontén T, Ussery DW, Aarestrup FM, Lund O. 2012. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol 50:1355–1361. doi: 10.1128/JCM.06094-11. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
6. Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, Villa L, Møller Aarestrup FM, Hasman H. 2014. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 58:3895–3903. doi: 10.1128/AAC.02412-14. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
7. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. 2012. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother 67:2640–2644. doi: 10.1093/jac/dks261. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

Articles from Microbiology Resource Announcements are provided here courtesy of American Society for Microbiology (ASM)

-