Abstract
Background Previous metabarcoding studies based on 16S rRNA sequencing in patients with extrapulmonary sepsis have found early pulmonary dysbiosis associated with a poor prognosis. To further discern this association, here we aimed to better characterize the pulmonary bacterial communities in these patients by leveraging metagenomics and to evaluate if the presence of antibiotic resistance genes (ARGs) could explain the higher mortality of the patients.
Material and methods Metagenomic sequencing was performed using the Nextera XT Library Prep Kit and HiSeq 4000 (Illumina Inc.) on tracheal aspirate samples that were obtained within 24 h from diagnosis from patients with extrapulmonary sepsis admitted to the Intensive Care Unit (ICU). Analysis involved MetaSpades for contig assembly, Kraken2 and Metaphlan4 for taxonomic classification, and CARD and GTDB-tk for ARGs annotation and assignment to the bacterial species. The relationship between the presence of antibiotic resistance and ICU mortality was evaluated using the Wilcoxon test and logistic regression models adjusting for clinical and demographic variables.
Results In total, 127 different ARGs were detected circumscribed only to seven patients. The most common ARGs found were from the antibiotic groups of aminoglycosides and beta-lactams, both present in most of the patients. These ARGs were found, almost entirely linked to Klebsiella pneumoniae, Escherichia coli, and Pseudomonas aeruginosa. The results also show a significant enrichment of ARGs among patients who died while admitted in the ICU (57%, 95% confidence interval [CI]: 18-90%) compared to surviving patients (20%, 95% CI: 7-40%) (p=0.022). Analyses adjusting for clinical and demographic variables did not alter this result.
Conclusion Metagenomic sequencing has allowed an unprecedented characterization of the sepsis lung microbiome showing that antibiotic resistance is common among these patients. The results also suggest a relationship between the early accumulation of ARGs in the lung of patients with extrapulmonary sepsis who die while admitted in the ICU. Studies in independent samples will be needed to validate our findings.
Introduction
Sepsis is a severe and life-threatening complication of infection characterized by a systemic inflammatory response and organ dysfunction (Singer et al., 2016). Despite significant advances in medical care, sepsis leads the causes of mortality in intensive care units (ICUs) and remains a major cause of morbidity among survivors (Herrán-Monge et al., 2019). The pathogenesis of sepsis is complex and multifactorial, involving interactions between host immune responses and the microbiome (Akrami & Sweeney, 2018; van der Poll et al., 2021).
Rapid sepsis diagnosis and identification of the causative pathogen is crucial for the optimal management of the disease, including the selection of appropriate antimicrobial therapy (Kumar et al., 2006). However, the traditional microbiological culture-based approach for pathogen identification and characterization has several limitations, including the low sensitivity and specificity, the requirement of a pure culture isolation, and the inability to detect fastidious or non-cultivable organisms (Handelsman, 2004). As a result, in >30% of cases the aetiologic pathogen is not identified by these methods (Novosad et al., 2016), reflecting the limitations of current culture-based microbiologic diagnostics.
Consequently, antimicrobial therapies that are initiated early after sepsis diagnosis or ICU admission cannot rely on culture-based information (Strich et al., 2020). However, in this context, the rising trend of antimicrobial resistance (AMR) constitutes a main problem among patients admitted to ICUs, especially for those staying longer periods and under high antibiotic pressure. As a figure, about 0.7 million people die annually worldwide from multidrug-resistant (MDR) strains of microorganisms and it is estimated to increase to 10 million by 2050 (de Kraker et al., 2016).
Our previous studies in mechanically ventilated ICU patients with extrapulmonary sepsis revealed that pulmonary dysbiosis within 8 h from sepsis diagnosis was associated with ICU mortality (Guillen-Guio et al., 2020). These observations offer promising options that could be technologically adapted to rapidly identify those patients at higher risk of death (Rodríguez-Pérez et al., 2022). However, they relied on 16S rRNA sequencing experiments which impedes deep characterization of the bacterial communities that could help to further substantiate the association.
To further discern the association of lung dysbiosis with the ICU mortality among sepsis patients, here we used metagenomic sequencing for improved characterization of the bacterial lung communities and for adding support to the possible enrichment of AMRs among the deceased patients.
Materials and Methods
Ethics statement
The study was approved by the Research Ethics Committee of the Hospital Universitario Nuestra Señora de Candelaria (PI-30/14) and performed according to The Code of Ethics of the World Medical Association (Declaration of Helsinki). An informed consent was obtained from all patients or from their legal representatives.
Study patients, sample preparation, and shotgun sequencing
Samples were collected and stored at -20 °C until use from tracheal aspirates from 32 mechanically ventilated patients (25 survivors, 7 deceased) admitted to the medical-surgical ICU within the first 8 h from sepsis diagnosis (Singer et al., 2016) between January 2015 and January 2019 and followed up until discharge or death. These were critical patients with a diagnosis of acute abdomen and constitute a subset of those with sufficient DNA material from aspirates whose clinical characteristics have been described elsewhere (Guillen-Guio et al., 2020; Rodríguez-Pérez et al., 2022).
DNA was extracted using the QIAmp UCP Pathogen Mini Kit (Qiagen). Following quantification with a Qubit 3.0 fluorometer using a High Sensitivity DNA Analysis Kit (Thermo Fisher Scientific) and sample normalization, libraries were prepared with the Nextera XT Library Prep Kit (Illumina Inc.). Library sizes and concentrations were assessed on a TapeStation 4200 (Agilent Technologies) and the High Sensitivity DNA Analysis Kit. Sequencing was performed using a HiSeq 4000 System (Illumina Inc.) using 150 base paired-end reads along with 1% of a PhiX control V3 (Illumina Inc.).
Bioinformatic analysis
We designed and deployed a bioinformatics pipeline as a collection of BASH scripts to analyze the obtained metagenomic sequence data (https://github.com/genomicsITER/metagenomics_suppmaterial). In brief, the pipeline consisted in the following stages. Sequences underwent first a quality control step using Trimmomatic v0.39 (Bolger et al., 2014). Also, as a preprocessing step, a rapid taxonomic assignment of reads was performed with Kraken2 v2.1.3 (Wood et al., 2019) to identify the human reads. All the identified human reads were removed from the analyses using KrakenTools v1.2 (Lu et al., 2022). On the remaining reads, the subsequent analysis steps included using MetaPhlAn4 (Blanco-Míguez et al., 2023) to determine the general composition of each sample. In parallel, microbial reads were assembled de novo using Spades v3.15.5 (Bankevich et al., 2012) with a specific parametrization for metagenomic assembly. The contigs generated from the assembly step were then binned with Metabat2 (Kang et al., 2019) to produce contig sets or bins, each one belonging to a different organism. In order to identify the AMRs, all binned contigs were annotated for ARGs based on the Comprehensive Antibiotic Resistance Database (CARD) 3.2.7 (Alcock et al., 2019). The CARD hits with coverage below 85% were excluded from further analysis to prevent false positives. Additionally, we used GTDB-tk v2.3.0 to classify each bin based on the Genome Database Taxonomy (GTDB) v214.1 (Parks et al., 2022) to link the identified ARGs to their respective taxonomy at the species level. Finally, CheckM (Parks et al., 2015) was used to estimate the genome completeness and contamination.
Statistical analysis
The association between the presence of ARGs and sepsis mortality in the ICU was tested using logistic regression models using R v4.3.0 (R Foundation for Statistical Computing, 2021). Subsequent logistic regression-based sensitivity analyses were performed to investigate the effect of clinical and demographic variables, controlling for the level of diversity as measured by the Shannon diversity index.
Data availability
The datasets analysed in the study are available from the corresponding author on reasonable request.
Results
Sequencing data yield an average of more than 32 million reads per sample (range of 10 to 70 million reads) (see Supplementary Table S1). After human DNA sequences removal, between 0.25% and 12.8% of the reads per sample identified as non-human were kept for downstream analyses. At this stage, we observed that a higher proportion of non-human reads in the sample was associated with ICU mortality in the cases (Odds Ratio [OR]: 1.37; 95% Confidence Interval [CI]: 1.03-2.21; p=0.035).
After metagenomic assembly and genome classification of the non-human reads, 14 genomes with an estimated completeness of 80%-97% and low contamination were assembled across patient samples. All these genomes were classified as Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa, Acinetobacter baumannii, Achromobacter xylosoxidans, Staphylococcus aureus, and Haemophilus influenzae, which correspond to the most abundant taxa detected at read level by MetaPhlAn4 in the samples.
ARGs were detected only in nine patients (127 ARGs in total), four of whom died in the ICU (Figure 1). We observed that the capacity to detect AMRs was associated with the proportion of non-human reads that were recovered from the samples (OR: 11.02; 95% CI: 1.94-144.58; p=0.023). CARD ontology and bin classifications indicated that the most frequent AMR findings were linked to aminoglycosides and beta-lactam antibiotics. These were identified, almost entirely, in the assembled genomes from E. coli (48), K. pneumoniae (35), and P. aeruginosa (31) (Figure 2).
Finally, we assessed if the abundance of AMRs was associated with ICU mortality in this cohort of patients. We found that the number of AMRs was significantly higher (OR: 1.19; 95% CI: 1.03-1.40; p=0.022) in deceased patients (57%, 95% CI: 18-90%) compared to surviving patients (20%, 95% CI: 7-40%). These results were robust to model adjustments for relevant demographic and clinical covariates in the sensitivity analysis (Table 1).
Discussion
Dysbiosis of the lung microbiome has potential implications for disease pathogenesis and therapeutic interventions (Akrami & Sweeney, 2018; Mukherjee & Hanidziar, 2018). Our previous findings showed reduced lung bacterial diversity at genus and species levels associated with extrapulmonary sepsis mortality (Guillen-Guio et al., 2020; Rodríguez-Pérez et al., 2022). Here, a metagenomics-based analysis allowed us to further characterize the lung dysbiosis in these patients, beyond assessing bacterial composition and diversity, by focusing on genomic assemblies and the identification of AMR profiles. Our analysis revealed that the abundance of AMRs in the lungs within the first hours of sepsis diagnosis correlated with ICU mortality rates, providing a qualitative survey of bacterial co-infections likely contributing to the sepsis severity in these patients and making it the first study linking lung microbiome AMRs to ICU mortality in sepsis patients.
Predominant Gram-negative bacteria like E. coli, P. aeruginosa, and K. pneumonia harboring diverse resistance mechanisms have been previously identified as common sepsis-causing pathogens (Bloomfield et al., 2023; Mu et al., 2023). Some of these species have been linked to gut-to-lung bacterial migration and AMR acquisition and the worsening of sepsis (Wheatley et al., 2022). A significant proportion of the detected resistance determinants in these species, as well as in most of the AMR-positive samples from our study, belongs to the aminoglycoside, beta-lactam, and peptide antibiotic classes, which have been linked to a substantial impact on the efficacy of antimicrobial treatment (Fatsis-Kavalopoulos et al., 2022; Foudraine et al., 2021). Independent studies have found that poor outcomes in mechanically ventilated critically ill patients were associated with increased bacterial burden in lung microbiome (Dickson et al., 2020). This agrees with our observations where a higher proportion of non-human reads was associated with ICU mortality in our patient series. Taken together, our results emphasize the significance of understanding bacterial pathogens in sepsis severity and early response strategies.
This work has a series of limitations which include a small sample size from a single ICU, preventing optimal confounding effect assessments, and the study’s focus on bacterial reads due to the prevalence of bacterial sepsis, which restricts conclusions regarding non-bacterial pathogens (Faria et al., 2018). In addition, we did not have the possibility to assess an independent sample or access to patient samples before the onset of sepsis to serve as controls. One cannot either conclude that the detected bacteria were pathogens causing the sepsis or if they casually contributed to the patient severity. Moreover, the study relied on a limited proportion of reads non-aligning to the human genome, limiting our ability to assemble highly continuous bacterial genomes which could have consequences in the capacity to improve the detection of AMRs. Despite the advancements in high-throughput sequencing, metagenomic data and antibiotic resistance profiles of lung samples are underrepresented in literature compared to gut or oral microbiomes (Dickson et al., 2016; Ibironke et al., 2020). Within the last decade, metagenomic applications have emerged for precise pathogen identification and gene annotation, albeit challenges like standardized protocols and complex data interpretation remain (Chiang & Dekker, 2020). We suggest that future metagenomic applications could potentially enhance early-response strategies to sepsis, anticipating a growing relevance in clinical applications.
Conclusion
Sepsis remains a major cause of morbidity and mortality in the ICUs worldwide and early diagnosis and identification of causative pathogens are crucial for optimal patient management. Metagenomics by next-generation sequencing could overcome some of the limitations inherent to the traditional microbiological culture-based approaches and provide important clinical information. Here we applied metagenomics to 32 tracheal aspirate samples from ICU patients with non-pulmonary sepsis to identify the AMRs and to link them to bacterial species. The study results support the association between AMR abundance and patient mortality.
Data Availability
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Funding
This work was supported by Instituto de Salud Carlos III [CB06/06/1088, PI14/00844, PI17/00610, FI18/00230, PI19/00141, CD22/00138 and the IMPaCT-Data programme IMP/00019] and co-financed by the European Regional Development Funds, “A way of making Europe” from the European Union; Ministerio de Ciencia e Innovación [RTC-2017-6471-1, AEI/FEDER, UE]; Cabildo Insular de Tenerife [CGIEU0000219140 and A0000014697]; Fundación Canaria Instituto de Investigación Sanitaria de Canarias [PIFUN48/18 and PIFIISC21/37]; Fundación DISA [OA23/074]; Wellcome Trust [221680/Z/20/Z]; and by the agreements with Instituto Tecnológico y de Energías Renovables (ITER) [OA17/008 and OA23/043]. For the purpose of open access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.