Abstract
Secondary bacterial and fungal infections are associated with respiratory viral infections and invasive mechanical ventilation. In Coronavirus disease 2019 (COVID-19), lung injury by SARS-CoV-2 and impaired immune response can provide a favorable environment for microorganism growth and colonization in hospitalized individuals. Recent studies suggest that secondary bacterial pneumonia is a risk factor associated with COVID-19. In Brazil, knowledge about microbiota present in COVID-19 patients is incipient. This work describes the microbiota of 21 COVID-19 patients admitted to intensive care units from two Brazilian centers. We identified respiratory, nosocomial and bacterial pathogens as prevalent microorganisms. Other bacterial opportunistic and commensal species are also represented. Virulence factors of these pathogenic species, metabolic pathways used to evade and modulate immunological processes and the interconnection between bacterial presence and virulence in COVID-19 progression are discussed.
Article Summary Line We identified respiratory, nosocomial and bacterial pathogens as prevalent microorganisms in 21 Brazilian COVID-19 patients admitted to Intensive Care Units. Pathogen virulence factors and immune response evasion metabolic pathways are correlated to COVID-19 severity.
Introduction
Severe COVID-19 cases are characterized by systemic hyperinflammation and immune dysfunction, which may lead to lung damage. Infection begins with the internalization of the viral complex through the interaction between SARS-CoV-2 spike protein and the angiotensin-converting enzyme 2 (ACE2) in the host cell. ACE2 expression reduces after SARS-CoV-2-infection, resulting in an increase in angiotensin II which lead to pathophysiological effects that include vasoconstriction, increased inflammation and fibrosis (1). Lung edema, endothelial and epithelial injuries are accompanied by an influx of neutrophils into the interstitium and bronchoalveolar space, resulting in impairment of arterial oxygenation (2). Infection of alveolar macrophages by SARS-CoV and neutrophil infiltration are important triggers of the cytokine storm, which is extensively associated with COVID-19 severity and mortality (3).
In the critical stage of COVID-19, cytokine storm is continuous, contributing to vascular congestion, complement cascade activation and over disseminated intravascular coagulation (4). Hypoxia is aggravated and invasive mechanical ventilation may be required for life maintenance. Under these conditions, variation in oxygen tension, alveolar ventilation, deposition of inhaled particles, blood flow and concentration of inflammatory cells are factors that directly affect lung microbial growth conditions (5, 6). In addition, irregular innate immune response in acute lung injury, mainly during the process of neutrophil recruitment, may favor bacterial infections, which are also regulated by virulence factors (7).
In Influenza and Severe Acute Respiratory Syndrome (SARS), secondary bacterial infection has been reported during intensive care unit admission and during use of invasive mechanical ventilation. Data in literature pointed out that occurrence of co-infection or secondary bacterial pneumonia was observed in 11% to 35% of individuals with laboratory confirmed influenza (8), and that Streptococcus pneumoniae and Staphylococcus aureus were present in 35% and 28% of patients, respectively. In Canada, Chlamydophila pneumoniae or Mycoplasma pneumoniae were abundant in up to 30% of SARS patients (9). Bacterial infection was also identified in COVID-19 cases (10 - 12). Studies conducted in China, the USA, Spain and Thailand, showed that secondary bacterial infection was present in 14% of patients, with occurrence of Streptococcus pneumoniae, Klebsiella pneumoniae and Haemophilus influenzae within 1–4 days of COVID-19 disease onset (13,14). According to Rawson et al. (2020)(15), the extensive use of broad-spectrum antibiotics can predispose COVID-19 patients to acquisition of bacterial infections and increase multidrug-resistance.
Colonizing species in patients infected with SARS-CoV-2, prevalence of bacterial infection and mechanism of co-or secondary-infection remain poorly understood. In Brazil, an investigation of bacterial microbiota in SARS-CoV-2 was conducted in only one patient and concise results showed the predominance of Lautropia, Prevotella and Haemophilus genera (16). In-depth work was not performed.
This study aims to determine the microbiota of 21 COVID-19 patients admitted to intensive care units from two Brazilian centers, focusing on prevalence of bacterial pathogens, and to correlate the microbiota with the immunological disorder characteristic of Coronavirus disease 2019.
Material and Methods
Study cohort and data collection
Twenty-one COVID-19 patients admitted to Intensive Care Units (ICU) from Hospital Naval Marcílio Dias (Rio de Janeiro) and Hospital Eduardo de Menezes (Belo Horizonte) were included in this study. Both hospitals are located in cities in the Southeast region of Brazil. All patients had positive RT-PCR tests for SARS-CoV-2 (Table 1). Samples identification (IDs) were modified to preserve the patient’s integrity. The study cohort included 12 males and 9 females and no age restrictions were applied (ranging from 37 to 89 years of age).The present study was approved by the Ethics Commission from Hospital Naval Marcílio Dias and Hospital Eduardo de Menezes (protocol number 32382820.3.0000.5256 and 31462820.3.0000.5149). Clinical and demographic characteristics of the 21 patients were obtained from medical records and summarized in Table1.
Summary of the clinical and demographic aspects of the 21 COVID-19 patients analyzed.
DNA extraction and sequencing
DNA was extracted using standard manufacturer’s protocols for the Qiagen QIAamp DNA Microbiome kit (Qiagen, Germany). Metagenomic libraries were constructed using the Nextera DNA Flex Library Preparation Kit (Illumina, USA). Sequencing was performed in a NextSeq 500 System with a NextSeq 500/550 High Output Kit v2.5 (300 Cycles) (Illumina, USA).
Bioinformatics analysis
Sequencing files were submitted to MG-RAST (version 4.0.3) (17) for taxonomical and functional inferences. Sequences were trimmed using default parameters and human sequences removed by screening against Homo sapiens. Species were identified from MG-RAST using RefSeq and functional abundances were obtained using SEED subsystem database. For both analyses, the Best Hit Classification criteria and alignment cutoff applied included an e-value > 10-5, a minimum identity of 60%, and a minimum alignment of 15. Low abundance species (≤ 50 absolute counts of the sum of all samples) were removed to avoid false-positives and minimize the error of taxonomical predictions (18). The number of species was resampled so that all samples had the same library size. Prevalent species were detected using the microbiome R package (19), considering only those that occurred in 75% of all patients. Phyloseq R package was used to estimate the Shannon diversity index (20).
Results and Discussion
The study included only adult patients, the majority of whom were above 65 years of age (42.8%). Males represent 57% of the patients sampled. The period between symptom onset and sample collection varied from 4 to 34 days. All individuals were admitted to intensive care units and invasive ventilatory support was necessary in 95.2% of the cases. No previous pulmonary impairment was reported, except for one individual. Most patients had at least one of the main comorbidities associated with disease severity. Diabetes represented 50% of the cases studied, followed by obesity and hypertension (33% in both) and cardiovascular disease (23.8%) (Table 1). Tracheal lavage samples with SARS-CoV-2 infectivity had higher abundances of bacterial respiratory and nosocomial species (Table 2). To check the relationship between microbial diversity and the proportion of microbial pathogens, we estimated diversity using the Shannon index and calculated the relative abundance of pathogens (Figure 1). The proportion of microbial pathogens was based on the relative abundance of individual pathogens divided by the total number of species in the microbiome. First, we observed a negative correlation (−0.9733005; p ≤ 0.05) between these two parameters, indicating that the prevalence of pathogens in all samples is inverse to diversity, causing a decrease in species richness. A high proportion of pathogens can be explained by the potential imbalance between microbial immigration and elimination, and regional growth conditions, as proposed by Dickson and collaborators (2015)(21).
Abundance of the fifty most prevalent species obtained from COVID-19 patients using MG-RAST inference.
Pearson correlation between the Shannon diversity index and the relative abundance of pathogenic species identified in the 21 COVID-19 patients.
Dickson and collaborators (2014, 2015) (6, 21) report that, while in healthy lungs the immigration and elimination processes are more determinant for microbial composition, in advanced lung diseases, pulmonary pathogenesis is a consequence of respiratory dysbiosis. Hyperventilation accelerates the influx of air-borne microbes and microbial elimination is reduced by bronchoconstriction and impaired mucociliary clearance. This condition provides a nutrient-rich growth medium and decreased oxygen concentration, influencing bacterial reproduction. In addition, these authors propose that the response to inflammation and endothelial and epithelial injury increases vascular permeability, promoting the escape of nutrient-rich fluid into the alveolar compartment, restoring the nutrient supply and selectively favoring the growth of specific lung pathogens. Our results suggest that COVID-19 reduces microbial diversity and increases the proportion of pathogens, similar to other acute and chronic lung diseases (22).
Almost 75% of the microbiota present in patients was composed of bacteria whose diversity was inferred from the fifty most prevalent species (Figure 2 and Table 2). Most respiratory pathogens were composed of Acinetobacter baumannii (corresponding to 16.5% of all bacterial species), Streptococcus pneumoniae (8.8%), Pseudomonas aeruginosa (3.6%) and Staphylococcus aureus (3.4%). Nosocomial pathogens mainly included fungi Candida albicans and Candida tropicalis, representing 11.2% and 8.1% of Eukarya domain, as well as the bacteria Enterococcus faecalis (1.5%) and Escherichia coli (0.5%).
Relative abundance of the 50 most prevalent species identified in the 21 COVID-19 patients. The characteristics of the individuals analyzed, such as age range, comorbidities and time (days) from onset of symptoms until sample collection were represented.
Furthermore, the patients analyzed in our work showed a diversity of other pathogenic microorganisms, such as the uncommon Listeria monocytogenes, Streptococcus (S. agalactiae and S. pyogenes), sporulating Bacillus (B. anthracis and B. thuringiensis) and Acidovorax spp., although less frequently (Figure 2). Curiously, the opportunistic microorganisms Delftia acidovorans and Coprobacillus spp. were found with an abundance of almost 6% and 2%, respectively. Delftia acidovorans is usually not pathogenic, but catheter-related bacteremia and cases of pneumonia with lung cavity formation has been observed (23, 24).
As expected, an elevated occurrence of the bacterial respiratory pathogen Streptococcus pneumoniae was found independent of the comorbidity, while in patients without comorbidity the occurrence of respiratory and nosocomial pathogens was drastically reduced, with exception of Candida albicans and Candida tropicalis (Figure 2 and Table 2). In contrast, Acinetobacter baumannii was significantly increased in patients with higher SARS-CoV2 infection period (RJ2 and RJ3 samples). No significant difference or direct association was observed between pathogen abundance and comorbidity type. However, the microbiota here identified reinforce the severity of infection under weakened health conditions.
Additionally, a large diversity of the oropharynx commensal species was identified, such as Rothia mucilaginosa, Parvimonas micra, Actinomyces odontolyticus, and Streptococcus anginosus (Table 2). R. mucilaginosa has been associated to pneumonia in patients with chronic obstructive pulmonary disease (COPD) (25). Furthermore, an association of co-infection by A. odontolyticus and P. micra or Streptococcus spp. to lung abscess and acute respiratory failure was reported (26, 27). Although respiratory disease caused by these commensals is a rare event, their presence in individuals with compromised immune responses should not be neglected.
The bacterial pathogens identified in this study, particularly Acinetobacter baumannii, Streptococcus pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, may promote virulence and evade the host immune response by biofilm production, induction of hemolysins, pneumolysin, phospholipases, iron acquisition factors, cytokines, adhesins, and complement system resistance, besides other virulence factors (28 - 30).
In pulmonary and catheter-associated infections, biofilm formation can be favored by microbial interaction in endotracheal tubes and mechanical-ventilation apparatus, resulting in a diverse microbial complex and enhanced antimicrobial resistance (31,32). Studies demonstrated that A. baumannii, Klebsiella pneumoniae, and Enterococcus faecium are also frequently recovered from endotracheal tubes (33-35). In addition, in clinical isolates with persistent infections, P. aeruginosa, S. aureus and C. albicans are the predominant biofilm-producers and species with more biofilm formation capacity are frequently observed to be multidrug-resistant (36-38). For some of these microorganisms, the increase of biofilm, production of bacterial virulence factors and antibiotic resistance are related to enrichment of macrophage secretory products in culture and oxygen-limiting conditions (39, 40). In S. aureus, for example, biofilm is related to accumulation of activated macrophages exhibiting anti-inflammatory and pro-fibrotic properties (41). Moreover, it was observed that, when in biofilm, bacterial species are more resistant to the immune response, evading neutrophil mediated phagocytosis (42).
Hemolysin, an important virulence factor of S. aureus, E. coli, and Enterococcus faecalis (43, 44), is a protein that causes erythrocyte lyse and disrupts the cell membrane, contributing to lung injury and pneumonia. In Streptococcus pneumoniae, when pulmonary surfactant is deficient, there is an expression of hemolysin, which is correlated with lung epithelial cell injury and induction of interleukin (IL)–8 (45).
Similar to hemolysin, pneumolysin is a pore-forming toxin, and most of its isoforms exhibit hemolytic activity. This protein is crucial for virulence and chronic bacteraemia of Streptococcus pneumoniae, and together with polysaccharide capsule proteins and adhesins reduce the S. pneumoniae entrapment in the pulmonary mucus, permitting adherence and inflammatory responses and activating complement and apoptosis (46).
For Candida species identified in this work, phospholipase and biofilm are described as the major virulence factors (47, 48). Phospholipases possess an important role in lung diseases, promoting change in membrane composition, stimulating chemokines and cytokines, and altering gas exchange and as a lung surfactant (49).
Hemolysin, phospholipase and biofilm are some of the many pathogenicity factors in Acinetobacter baumannii. Successful fitness of this bacteria, however, is a consequence of a large virulence repertoire, also composed by outer membrane proteins, secretion systems, surface adhesins, glycoconjugates and iron-chelating activities, which can explain its multi-resistance to antibiotics in nosocomial infections (28).
Biofilm-based infections are significantly less susceptible to antimicrobial agents and their treatment is extremely difficult. According to Sanchez et al., 2013 (37), the investigation of biofilm forming capacity in patients with diverse infective sources revealed that Staphylococcus aureus, Acinetobacter baumannii, Pseudomonas aeruginosa, Klebsiella pneumoniae, and Escherichia coli are strong biofilm producers, at levels greater than or equal to Staphylococcus epidermidis, a positive biofilm producing strain. The greatest number of biofilm producing strains corresponded to P. aeruginosa and S. aureus species.
Interestingly, the metabolic profile obtained in our study correlated to virulence factors of the species described, particularly biofilm (Tables 3-6). Within Central Carbohydrate Metabolism (Table 3), the TCA cycle was the predominant pathway (18%), with abundance of dihydrolipoamide dehydrogenase enzyme (Table 4). According to Lu and authors (2019) (50), alteration in the metabolites from the TCA cycle, amino acid metabolism, and glycerolipid pathways were observed during biofilm production. Furthermore, characterization of dihydrolipoamide dehydrogenase from Streptococcus pneumoniae showed that this enzyme plays an important role in pneumococcal infection, being necessary for the survival and capsular polysaccharide production of pneumococci within the host (51).
Functional pathway of Central Carbohydrate Metabolism, related to virulence factors reported for the species identified in this work.
Functional profile of TCA_Cycle pathway, related to virulence factors reported for the species identified in this work.
Besides enzymes associated with biofilm formation, we identified adhesins, permeases for antibiotic export, resistance and toxins, and superantigens belong to the Virulence, Disease and Defense functional category (Table 5). The most abundant adhesins are fibronectin/fibrinogen-binding proteins, although high numbers of sortases, related to attachment of specific proteins to the cell wall, were also found. Toxins and superantigens are composed of the streptolysin S locus, which produces the beta-haemolytic phenotype. Streptolysin S causes impaired phagocytic clearance and promotes epithelial cell cytotoxicity, enhancing the effects of M protein and streptolysin O (52). In the Membrane Transport category, protein secretion system Types I, V, and VI are highlighted, with abundance of exoprotein involved in heme utilization, cytolysin activator and adhesin (Table 6).
Functional profile of Virulence, Disease and Defense pathway, related to virulence factors reported for the species identified in this work.
Functional profile of Membrane Transport pathway, related to virulence factors reported for the species identified in this work.
The scenario highlights the high occurrence of potentially virulent pathogens, which can contribute to severity of lung infection in Coronavirus Disease 2019 patients admitted to intensive care units. Moreover, the prevalence of respiratory pathogens observed in COVID-19 patients analyzed can be orchestrated mainly by the complex and distinct immune events in response to lung damage, as well as by the ecological model proposed by Dickson and collaborators (2014, 2015) (26, 25), which indicate that several lung diseases can alter the growth of local microbiota, leading to an increase in bacterial abundance. A cellular model to explain the bacterial species described in our study correlated to processes that aggravate COVID-19 is presented (Figure 3).
Model for an interplay of the bacteria found in COVID-19 patients, virulence factors and free-iron in immune response and disease severity. In severe lung diseases, hyperventilation contributes to microbial influx and reduces the elimination of microorganisms. Tecidual injury promotes a rich environment for growth of specific lung pathogens into the alveolar compartment. Bacterial pathogens evade the host immune response by biofilm production, induction of hemolysins, iron acquisition factors, adhesins, and other virulence factors. The SARS-CoV-2 and bacterial pathogens promote erythrocyte lysis, causing heme(Hb)/iron(Fe) liberation. Heme uptake is required for bacterial host colonization and increased production of virulence factors. Bacterial virulence factors, as well as free-heme, activates the inflammasome, contributing to neutrophil inflammation, cytokine elevation and pathogenesis processes. This figure was generated using the BioRender (available at https://biorender.com/).
Impairment in neutrophil function, macrophage depletion and excessive inflammation has been reported as an important factor in the progression of lower respiratory tract pathologies (53). It was demonstrated that alveolar macrophage depletion could facilitate the bacterial infection by the establishment of a niche for secondary Streptococcus pneumoniae infection, altering cellular innate immunity and resulting in lethal pneumonia (54). Additionally, neutrophilia plays an important role in bacterial pulmonary diseases, including those caused by Streptococcus pneumoniae, Klebsiella pneumoniae, Haemophilus influenzae, and Staphylococcus aureus (55). Interestingly, Sodhi and collaborators (2019) (56) demonstrated that impaired neutrophil and innate immune response might be mediated by variations in the activity of angiotensin-converting enzyme 2 (ACE2). In bacterial lung infection, when pulmonary active ACE2 levels vary dynamically, neutrophil influx increases, compromising the host-defense capability and heightened inflammatory response, with subsequent elevation of infection severity and mortality
Similar to Coronavirus Disease 2019, patients with chronic obstructive pulmonary disease (COPD) show vascular permeabilization, platelet activation, thrombosis, adhesion molecule expression, leukocyte recruitment, and complement activation, which are also observed in injury-derived free heme (57). Su and authors (2020)(58) observed increased heme and heme oxygenase (HO)-1 in the blood system of severe COVID-19 patients. Moreover, an additional role of hemoglobinopathy, hypoxia and hyperferritinemia has been proposed in worsening COVID-19 infection (59).
Increase of cell-free hemoglobin (Hb) and HO-1 expression in lung tissue and bronchoalveolar lavage fluid has been correlated to COPD and acute respiratory distress syndrome (ARDS) severity (60). Simon and authors (2009)(61) demonstrated an ACE-like activity of Hb in the conversion of angiotensin I to its active metabolites induced by ferrous- and ferryl-Hb, with a possible effect on vasoconstriction. Moreover, intratracheal administration of heme in mice led to alveolar-capillary barrier dysfunction and increased alveolar permeability, contributing to acute lung and ARDS (62). Curiously, heme uptake by pathogens plays an important role during bacterial infection. These microorganisms induce erythrocyte hemolysis caused by pore-forming toxins, such as hemolysin and phospholipases. Free-heme and iron acquired from erythrocyte disruption is required for invasion, growth and successful pathogen colonization in the host (63). In addition, Dutra and collaborators (2014)(64) showed that heme promotes the processing of caspase-1, inflammasome activation and IL-1β secretion by macrophages that participate in the immune response induced by hemolysis, contributing to inflammation and pathogenesis. Strategies for interfering in inflammasome activation by bacterial pathogens have evolved (65), and together with SARS-Cov2, may contribute to a hyperactivated inflammatory response (66).
Interdependence between iron acquisition and virulence factors for bacterial pathogenesis have been demonstrated. In Streptococcus pneumoniae, iron is a critical factor for regulation of several enzymes. During pneumonia and bacteremia, the transcription of hemolysin, toxins, cell wall hydrolases and biofilm formation are increased in a positive Fe-dependent regulation (67). For Pseudomonas aeruginosa, high iron concentrations are correlated to stimulation of aggregation, adhesion and biofilm formation, in a positive-modulation manner (68). The acquisition of heme-iron bound to erythrocyte hemoglobin is the preferred source of iron during S. aureus infection initiation, as well as for its proliferation. In the presence of heme, regulatory systems induced exotoxins such as cytolysins, proteases, and other virulence factors, promote Staphylococcal cellular adherence to host cells and resistance to neutrophil killing. Furthermore, iron and oxygen conditions have been shown to influence the TCA cycle activity of Staphylococcus aureus (69 - 71). Free-iron uptake by ferritin and haemolytic activity is employed by Candida albicans. Interestingly, it has been speculated that co-habitation of C. albicans with both commensal and pathogenic bacteria, in a synergistic interaction, enables C. albicans to obtain iron more easily (72).
In conclusion, the presence of respiratory pathogens, nosocomial bacteria and opportunistic microorganisms during colonization of the host, secrete an extensive diversity of virulence factors. These molecules favor hemolysis, cytokine activation and inflammation. The presence of free-iron exacerbates bacterial pathogenesis in cyclic connection with the worsening of the immune response, favoring the resistance of bacterial species. These findings provide intriguing insights into SARS-CoV2-bacterial interactions contributing to COVID-19 severity.
Data Availability
NGS data generated in our study is publicly available in SRA-NCBI (www.ncbi.nlm.nih.gov/sra), Bioproject accession PRJNA683652.
Data availability
NGS data generated in our study is publicly available in SRA-NCBI (www.ncbi.nlm.nih.gov/sra), Bioproject accession PRJNA683652.
Author Bio
F.M.C is currently a postdoctoral fellow in the Laboratory of Bioinformatics (LABINFO) of the National Laboratory of Scientific Computing (LNCC). Research interests comprise metagenomic (clinical and environmental) and comparative genomics.
Acknowledgments
This work was developed under the frameworks of Corona-ômica-RJ (FAPERJ = E-26/210.179/2020) and Rede Corona-ômica BR MCTI/FINEP (FINEP = 01.20.0029.000462/20, CNPq = 404096/2020-4). A.T.R.V. is supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico - CNPq (303170/2017-4) and FAPERJ (E-26/202.903/20). R.S.A. is supported by CNPQ 312688/2017-2, 439119/2018-9; FAPERJ 202.922/2018 .C.C.C is supported by FAPERJ (E-26/202.791/2019) and S.P.C.B is supported by FAPERJ (E-26/010.000168/20).
Footnotes
R
References
