Nosocomial transmission of extensively drug resistant Acinetobacter baumannii strains in a tertiary level hospital

Abstract

Acinetobacter baumannii is an opportunistic infectious agent that affects primarily immunocompromised individuals. A. baumannii is highly prevalent in hospital settings being commonly associated with nosocomial transmission and drug resistance. Here, we report the identification and genetic characterization of A. baumannii strains among patients in a tertiary level hospital in Mexico. Whole genome sequencing analysis was performed to establish their genetic relationship and drug resistance mutations profile. Ten genetically different, extensively drug resistant strains were identified circulating among seven wards. The genetic profiles showed resistance primarily against aminoglycosides and beta-lactam antibiotics. Importantly, no mutants conferring resistance to colistin were observed. The results highlight the importance of implementing robust classification schemes for advanced genetic characterization of A. baumannii clinical isolates and simultaneous detection of drug resistance markers for adequate patient’s management in clinical settings.

About the Authors:

Paúl Alexis López-Durán

Roles Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliations Departamento de Bioinformática y Biotecnología Genómica, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Ciudad de México, México, Laboratorio de Investigación en Enfermedades Infecciosas; Unidad de Investigación Multidisciplinaria, Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México, México City, Edo. de México, México, Facultad de Ciencias de la Salud, Universidad Anáhuac, Campus Norte, Naucalpan de Juárez, Estado de México, México

Salvador Fonseca-Coronado

Roles Conceptualization, Investigation, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Laboratorio de Investigación en Enfermedades Infecciosas; Unidad de Investigación Multidisciplinaria, Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México, México City, Edo. de México, México

ORCID logo http://orcid.org/0000-0003-3853-6313

Lucila Maritza Lozano-Trenado

Roles Conceptualization, Data curation, Investigation, Methodology, Resources, Supervision, Validation

Affiliation: Escuela Militar de Graduados de Sanidad, Secretaría de la Defensa Nacional, Ciudad de México, México

Sergio Araujo-Betanzos

Roles Conceptualization, Data curation, Investigation, Resources, Supervision

Affiliation: Escuela Militar de Graduados de Sanidad, Secretaría de la Defensa Nacional, Ciudad de México, México

Deniria Alejandra Rugerio-Trujillo

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Supervision

Affiliation: Hospital Central Militar, Secretaría de la Defensa Nacional, Ciudad de México, México

Gilberto Vaughan

Roles Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Facultad de Ciencias de la Salud, Universidad Anáhuac, Campus Norte, Naucalpan de Juárez, Estado de México, México

Elsa Saldaña-Rivera

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Escuela Militar de Graduados de Sanidad, Secretaría de la Defensa Nacional, Ciudad de México, México

Introduction

Acinetobacter baumannii is a Gram negative coccobacillus that affects primarily immunocompromised individuals. Patients undergoing intensive care, emergency surgery, mechanical ventilation, tracheal intubation, central vein catheterization, urinary catheter, continuous renal replacement therapy, and other invasive procedures are at increased risk of infection [1–3].

The World Health Organization lists the carbapenem resistant A. baumannii on the pathogen critical priority list [4], making it an increasingly growing public health problem worldwide. Multidrug (MDR) resistance among A. baumannii strains is driven by rapid gene mutations or transfer of exogenous resistance genes by mobile genetic elements, such as plasmids, transposons or insertion sequences. Resistance to beta-lactam and aminoglycosides is also a rapidly growing public health problem globally. Likewise, circulation of carbapenem-resistant A. baumannii strains aggravates the antibiotic resistance scenery, making patient’s management challenging [5]. As a result, occurrence of transmission networks associated with MDR in clinical settings has become a growing health issue in different regions of the world [6]. MDR refers to strains that exhibit resistance to more than three or more antimicrobial drug classes, while extensively drug resistant (XDR) strains display resistance to all but two drug classes [7]. Pan-drug resistance (PDR) indicates resistance to all drug classes including carbapenems, colistin, and other polymyxins [8].

Here, we describe the circulation of ten genetically distinguishable strains among patients in seven wards in a third level hospital. Timely identification of transmission networks of multidrug resistant strains is of the utmost importance to prevent the spread of highly virulent bacteria and appropriate patient’s management.

Methods

Bacterial isolates and susceptibility testing

Ten isolates were obtained from patients hospitalized at the “Hospital Central Militar” in Mexico City collected between May and November of 2018. A written informed consent was obtained from each patient. The study was approved by the Research Ethics Committee at Hospital Central Militar (approval no. 045/2017). Patients for whom sample isolates were identified as Acinetobacter spp. were included in the study. Demographic and clinical data is listed on Table 1.

[Figure omitted. See PDF.]

Table 1. Patients’ clinical characteristics.

https://doi.org/10.1371/journal.pone.0231829.t001

Bacterial identification and antibiotic minimum inhibitory concentration (MIC) testing were performed using the Vitek 2 Compact system (bioMérieux, Basingstoke, UK) following the manufacturer's instructions.

Isolation of DNA

A. baumannii genomic DNA was isolated and purified using the QIAamp DNA Mini Kit (Qiagen) following the manufacturer’s recommendations. The Quant-iT Pico Green dsDNA Reagent (Invitrogen, USA) was used to determine DNA concentration. Fluorescence was quantified using the Hidex Sense Microplate Reader. Final DNA concentration was adjusted at 0.2 ng/μl.

Whole genome sequencing (WGS)

Total DNA was used to perform sequencing on a MiSeq instrument (Illumina, San Diego, CA) using 2x300 Nextera XT library preparation kit (Illumina, San Diego, CA), following the manufacturer’s instructions. Average library sizes were evaluated using a High Sensitivity DNA kit on the Agilent 2100 Bioanalyzer according to the manufacturer’s instructions. The obtained files are available at the NCBI SRA database under BioProject ID PRJNA587517. Quality of sequencing reads was determined using FASTQC version 0.11.8 [9]. Reads of poor quality, adapter sequences removal as well as sizing was performed using Trimmomatic version 0.39 [10]. FASTQC was used to evaluate quality scores before trimming. Kmergenie version 1.7051 [11] was used to determine the best kmer for genome assembly. All sequences were assembled using SPAdes version 3.13.0 [12] and ABySS version 2.1.5 [13]. QUAST was used to derive assembly metrics including N50, L50 and number of contigs [14]. All sequences were annotated with PROKKA version 1.13.7 [15]. KmerFinder 3.1 was used to determine the species (https://cge.cbs.dtu.dk/services/KmerFinder/) [16].

Multilocus sequencing typing (MLST)

All loci (gltA, gyrB, gdhB, recA, cpn60, gpi, rpoD) included in the Oxford and Pasteur schemes (cpn60, fusA, gltA, pyrG, recA, rplB and rpoB) were derived from WGS originated from all ten isolates. Clustering was performed using UPGMA dendrograms as implemented in the Orange suite (https://orange.biolab.si).

Core genome multilocus sequencing typing (cgMLST)

Allele calling was performed using the chewBBACA algorithm [17], with the corresponding A. baumannii training set (https://github.com/mickaelsilva/prodigal_training_files). Additionally, ten A. baumannii strains (GenBank accession numbers CP014528.1, CP018254.1, CP035672.1, CP027611.1, CP021326.1, CP033869.1, CP040259.1, CP041035.1, CP046654.1, CP009257.1) were used as reference sequences. The cgMLST analysis included a total of 2390 alleles across the genome. Tree visualization was carried out with PHYLOViZ [18].

Characterization of drug resistant mutations

Characterization of resistance-mutations was performed using ResFinder v3.1.0 [19].

Results

Clinical data

Samples were collected between May and November 2018 from patients hospitalized for more than five days at seven different wards. Patients’ clinical information is shown in Table 1.

Genetic relatedness

All ten isolates were subjected to WGS and analyzed for genetic relatedness and drug resistance profile. The MLST analysis according to Pasteur scheme showed that all strains belonged to ST2 (Fig 1A). Additionally, two unrelated reference strains (CP021326.1 and CP035672.1) were also classified as ST2. All other reference strains belonged to other distinctive ST.

[Figure omitted. See PDF.]

Fig 1. Genetic analysis of Acinetobacter baumannii strains.

WGS were used to extract all loci required for Pasteur, Oxford schemes and cgMLST approach. The corresponding dendrograms were generated based on genetic distances among isolates and reference sequences. A) Pasteur scheme classified all isolates belonging to one distinctive sequence type (ST2). B) Oxford scheme grouped all outbreak sequences into two different sequence types (ST369 and ST1837). The unrelated sample derived from patient P006 was classified as ST208 along with reference strain CP021326.1. All other unrelated reference strains were classified as different ST. C) The cgMLST analysis showed a main branch with two distinctive clusters. Samples P004 and P008 formed a pair included in the small cluster while two pair of sequences (P010:P012 and P011:P007) were included in the major cluster along with samples P001, P003 and P009. Unrelates train P006 remained genetically distant from all other clinical samples.

https://doi.org/10.1371/journal.pone.0231829.g001

The Oxford scheme classified the outbreak strains into two different types. Samples P004, P007, P008 and P009 belonged to ST369 while samples P001, P003, P010, P011 and P012 were classified as ST1837 (Fig 1B). Additionally, samples P006 and reference strain CP021326.1 grouped together as ST208. Sample P006 was recovered from a patient originally diagnosed with A. baumannii infection in a hospital in Tijuana, Baja California, who was subsequently referred to the Hospital Central Militar in Mexico City. Thus, this isolate was expected, to a certain degree, to remain genetically distant from the main cluster identified in our hospital.

cgMLST analysis using WGS was able to discriminate strains from each other with better resolution. All samples grouped together in one branch exhibiting two distinctive clusters. Cluster 1 included two sample pairs (P010:P012 and P011:P007) besides samples P001, P003 and P009 (Fig 1C). Cluster 2 include a pair composed by samples P004 and P008. All samples were closely genetically related but different with the exception of sample P006. Sample P006 was not part of any of the clusters in the main branch. This was expected since this patient was infected by a different source in a different geographical region, and therefore unrelated to all other cases. Distance Matrix between all isolates corroborated the phylogenetic findings (S1 Fig).

Antibiotic susceptibility profiles

The MIC method displayed a high degree of drug resistance among these isolates. All ten strains were classified as XDR according to the scheme by Magiorakos et al. in 2012 [20]. All strains displayed resistance to penicillin’s, cephalosporins, carbapenems, fluoroquinolones, aminoglycosides, and sulfonamides. The complete susceptibility pattern is shown in Fig 2A.

[Figure omitted. See PDF.]

Fig 2. Antibiotic resistance profile.

A) All strains were subjected to antibiotic resistance testing using amikacin, ampicillin/sulbactam, cefepime, cefotaxime, ceftazidime, ceftriaxone, ciprofloxacin, gentamicin, levofloxacin, meropenem, tetracycline, tobramycin and trimethoprim/sulfamethoxazole. Two distinctive patterns were observed. All samples exhibited the same profile. Sample P010 had an incomplete antibiotic resistance testing; and therefore, could not be classified accordingly. B) Drug resistance mutation detection was inferred from WGS. Six distinctive patterns were identified.

https://doi.org/10.1371/journal.pone.0231829.g002

Likewise, the corresponding sequence analysis confirm the presence of drug resistance mutations confirming the XDR nature of these strains. The main antibiotic resistance genes identified among these isolates were associated with aminoglycosides (genes: aac(6')-Ib3, aac(6')-Il, aadA1, aph(3')-Ia, aph(3'')-Ib, aph(6)-Id, armA) and beta-lactam antibiotics (blaADC-25, blaOXA-237, blaOXA-66, blaOXA-72, blaTEM-171, blaTEM-1D). In addition, mutations associated with resistance to fluoroquinolones (aac(6')-Ib-cr), MLS—macrolide, lincosamide and streptogramin B (msr(E), mph(E)), phenicol (catB8), sulfonamides (sul1) and tetracycline (tet(B)) were also identified, although with less frequency (Fig 2B). Importantly, no mutations associated with colistin resistance was identified.

Discussion

Here we describe the molecular characterization and drug resistance profiles of ten XDR A. baumannii strains associated with nosocomial transmission in a tertiary hospital, the data presented here is not representative of the molecular epidemiology of A. baumannii in the country. As reported previously, it has been observed that A. baumannii is transmitted primarily by contact or invasive procedures, such as surgery, mechanical respirator or tracheostomy [6,21]. It has been also described that burned patients are particularly susceptible to an Acinetobacter infection during hospitalization [22]. All patients enrolled in this study underwent an invasive procedure and/or were immunocompromised.

The extensive drug resistance profile displayed by the analyzed strains showed a remarkable resistance to primarily aminoglycosides and β-lactam antibiotics. However, resistance to fluoroquinolones, phenicol, sulfonamides, tetracycline was also commonly observed among the isolates. Less frequent was the identification of mutants conferring resistance to macrolides/lincosamides and streptogramin. Although the P010 antibiogram is not complete, it can be inferred by genetic analysis that its profile is very similar to the other strains. This resistance pattern among XDR strains has been previously described in other countries [23,24]. Different studies have reported the circulation, and corresponding transmission of A. baumannii strains in other clinical settings in Mexico [1,25,26]. However, and to the best of our knowledge, this is the first report describing the circulation of XDR strains in the country.

A. baumannii infections that are resistant to carbapenems have been associated with high mortality rates. The most important mechanism in A. baumannii for carbapenem resistance is the production of carbapenem-hydrolyzing class D OXA-type β-lactamases [27]. Of note, most of the strains are carrying mutations associated with resistance in blaOXA-66 and blaOXA-72.

The blaOXA-66 is the gene that encode for the most frequent carbapenemase in the world, coexisting with others carbapenemases just like the blaOXA-23 and blaNDM-1 genes [28]. blaOXA-72 is the main allelic variant within the blaOXA24/40 group, and strains producing OXA-72 has previously been reported in hospitals and other sources in Eastern Europe, Southern Asia and South America [27]. Special attention must be point out on this element that represent a flexible way to transfer carbapenem resistance.

Regarding the phylogenetic classification, in our study Pasteur scheme showed no diversity in the strains being all assigned to the sequence type 2, it has largely showed that this scheme exhibits less discerning capacity among closely related strains [29]. According to the Oxford Scheme, two clusters of transmission were observed among patients sharing XDR A. baumannii strains (ST1837 y ST369), The ST183 and ST369 has been reported at other hospitals in Mexico, the regional Hospital Ignacio Zaragoza, Mexico City (n = 18), the Dr. José Eleuterio González University Hospital, Monterrey (n = 13) and the Hospital Central Monterrey (n = 1) [1,25,26]. About the foreign strain ST208, this has been reported in the Nuevo Leon University Hospital [26].

WGS analysis showed the circulation of closely genetically related strains. However, no clonal transmission was identified in this study. This could be explained by multiple sources of infections in this setting. Similarly, no common procedure or direct contact between cases was reported. The enrolling of new cases is ongoing, and the corresponding clinical investigation is expected to be continued. It is noteworthy that other conventional typing approaches for A. baumannii did not provide sufficient resolution to accurately characterized these isolates. Analyses for both, the Oxford and Pasteur typing schemes, were carried out using loci extracted from WGS. Neither typing approach resulted in sufficient discrimination among these isolates to carry out the clinical investigation.

Implementation of WGS as advanced molecular characterization of infectious agents in tertiary hospitals could improve both, identification and genetic characterization including resistance profiles, facilitating outbreak investigations and molecular surveillance. This information will likely provide insights on the molecular epidemiology and routes of transmission exploited by different infectious agents in clinical settings.

Supporting information

[Figure omitted. See PDF.]

S1 Fig.

https://doi.org/10.1371/journal.pone.0231829.s001

(TIFF)

Citation: López-Durán PA, Fonseca-Coronado S, Lozano-Trenado LM, Araujo-Betanzos S, Rugerio-Trujillo DA, Vaughan G, et al. (2020) Nosocomial transmission of extensively drug resistant Acinetobacter baumannii strains in a tertiary level hospital. PLoS ONE 15(4): e0231829. https://doi.org/10.1371/journal.pone.0231829

References

1. Alcántar-Curiel MD, Rosales-Reyes R, Jarillo-Quijada MD, Gayosso-Vázquez C, Fernández-Vázquez JL, Toledano-Tableros JE, et al. Carbapenem-Resistant Acinetobacter baumannii in Three Tertiary Care Hospitals in Mexico: Virulence Profiles, Innate Immune Response and Clonal Dissemination. Front Microbiol 2019;10:2116. pmid:31616391

2. Kumburu HH, Sonda T, van Zwetselaar M, Leekitcharoenphon P, Lukjancenko O, Mmbaga BT, et al. Using WGS to identify antibiotic resistance genes and predict antimicrobial resistance phenotypes in MDR Acinetobacter baumannii in Tanzania. J Antimicrob Chemother 2019;74:1484–93. pmid:30843063

3. Zhou H, Yao Y, Zhu B, Ren D, Yang Q, Fu Y, et al. Risk factors for acquisition and mortality of multidrug-resistant Acinetobacter baumannii bacteremia: A retrospective study from a Chinese hospital. Medicine (Baltimore) 2019;98:e14937. https://doi.org/10.1097/MD.0000000000014937.

4. Mulani MS, Kamble EE, Kumkar SN, Tawre MS, Pardesi KR. Emerging Strategies to Combat ESKAPE Pathogens in the Era of Antimicrobial Resistance: A Review. Front Microbiol 2019;10:539. pmid:30988669

5. Beigverdi R, Sattari-Maraji A, Emaneini M, Jabalameli F. Status of carbapenem-resistant Acinetobacter baumannii harboring carbapenemase: First systematic review and meta-analysis from Iran. Infect Genet Evol 2019;73:433–443. pmid:31176030

6. Takoi H, Fujita K, Hyodo H, Matsumoto M, Otani S, Gorai M, et al. Acinetobacter baumannii can be transferred from contaminated nitrile examination gloves to polypropylene plastic surfaces. Am J Infect Control 2019;0. https://doi.org/10.1016/j.ajic.2019.04.009.

7. Leski TA, Vora GJ, Barrows BR, Pimentel G, House BL, Nicklasson M, et al. Molecular Characterization of Multidrug Resistant Hospital Isolates Using the Antimicrobial Resistance Determinant Microarray. PLoS ONE 2013;8:e69507. pmid:23936031

8. Lee C-M, Lim H-K, Liu C-P, Tseng H-K. Treatment of pan-drug resistant Acinetobacter baumannii. Scand J Infect Dis 2005;37:195–9. pmid:15849052

9. Brown J, Pirrung M, McCue LA. FQC Dashboard: integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics 2017;33:3137–9. pmid:28605449

10. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014;30:2114–20. pmid:24695404

11. Chikhi R, Medvedev P. Informed and automated k-mer size selection for genome assembly. Bioinformatics 2014;30:31–7. pmid:23732276

12. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. J Comput Biol 2012;19:455–77. pmid:22506599

13. Jackman SD, Vandervalk BP, Mohamadi H, Chu J, Yeo S, Hammond SA, et al. ABySS 2.0: resource-efficient assembly of large genomes using a Bloom filter. Genome Res 2017;27:768–77. pmid:28232478

14. Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 2013;29:1072–5. pmid:23422339

15. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinforma Oxf Engl 2014;30:2068–9. https://doi.org/10.1093/bioinformatics/btu153.

16. Hasman H, Saputra D, Sicheritz-Ponten T, Lund O, Svendsen CA, Frimodt-Moller N, et al. Rapid Whole-Genome Sequencing for Detection and Characterization of Microorganisms Directly from Clinical Samples. J Clin Microbiol 2014;52:139–46. pmid:24172157

17. Silva M, Machado MP, Silva DN, Rossi M, Moran-Gilad J, Santos S, et al. chewBBACA: A complete suite for gene-by-gene schema creation and strain identification. Microb Genomics 2018;4. https://doi.org/10.1099/mgen.0.000166.

18. Nascimento M, Sousa A, Ramirez M, Francisco AP, Carriço JA, Vaz C. PHYLOViZ 2.0: providing scalable data integration and visualization for multiple phylogenetic inference methods. Bioinformatics 2017;33:128–9. pmid:27605102

19. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother 2012;67:2640–4. pmid:22782487

20. Magiorakos A-P, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect 2012;18:268–81. pmid:21793988

21. Maragakis LL, Perl TM. Acinetobacter baumannii: Epidemiology, Antimicrobial Resistance, and Treatment Options. Clin Infect Dis 2008;46:1254–63. pmid:18444865

22. Lima WG, Silva Alves GC, Sanches C, Antunes Fernandes SO, de Paiva MC. Carbapenem-resistant Acinetobacter baumannii in patients with burn injury: A systematic review and meta-analysis. Burns 2019;45:1495–508. pmid:31351820

23. Pinzón JO, Mantilla JR, Valenzuela EM, Fernández F, Arturo Álvarez C, Osorio EDEJ. Caracterización molecular de aislamientos de Acinetobacter baumannii provenientes de la unidad de quemados de un hospital de tercer nivel de Bogotá. Infectio 2006;10:71–8.

24. Saavedra SY, Nuñez JC, Pulido IY, Gonzalez EB, Valenzuela EM, Reguero MT, et al. Characterisation of carbapenem-resistant Acinetobacter calcoaceticus–A. baumannii complex isolates in a third-level hospital in Bogotá, Colombia. Int J Antimicrob Agents 2008;31:389–91. pmid:18308519

25. Bocanegra-Ibarias P, Peña-López C, Camacho-Ortiz A, Llaca-Díaz J, Silva-Sánchez J, Barrios H, et al. Genetic characterisation of drug resistance and clonal dynamics of Acinetobacter baumannii in a hospital setting in Mexico. Int J Antimicrob Agents 2015;45:309–13. pmid:25561030

26. Gonzalez-Villoria AM, Tamayo-Legorreta E, Garza-Ramos U, Barrios H, Sanchez-Pérez A, Rodríguez-Medina N, et al. A Multicenter Study in Mexico Finds Acinetobacter baumannii Clinical Isolates Belonging to Clonal Complexes 636 B (113 B) and 92 B Harboring OXA-72, OXA-239, and OXA-469. Antimicrob Agents Chemother 2016;60:2587–8. pmid:26833167

27. Jia H, Sun Q, Ruan Z, Xie X. Characterization of a small plasmid carrying the carbapenem resistance gene blaOXA-72 from community-acquired Acinetobacter baumannii sequence type 880 in China. Infection and Drug Resistance. 2019;Volume 12:1545–1553. pmid:31239730

28. Lv W, Zhang X, Hou M, Han D, Li Y, Xiong W. Draft genome sequence of an OXA-23, OXA-66, ADC-25 and TEM-1D co-producing Acinetobacter baumannii ST195 isolated from a patient with neonatal pneumonia in China. Journal of Global Antimicrobial Resistance. 2019;16:1–3. pmid:30445210

29. Gaiarsa S, Batisti Biffignandi G, Esposito E, Castelli M, Jolley K, Brisse S et al. Comparative Analysis of the Two Acinetobacter baumannii Multilocus Sequence Typing (MLST) Schemes. Frontiers in Microbiology. 2019;10.