- 1Department of Cardiovascular Surgery, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 2Genoxor Medical Science and Technology Inc., Zhejiang, China
Respiratory infections are complicated biological processes associated with an unbalanced microbial community and a wide range of pathogens. To date, robust approaches are still required for distinguishing the pathogenic microorganisms from the colonizing ones in the clinical specimens with complex infection. In this study, we retrospectively analyzed the data of conventional culture testing and metagenomic next-generation sequencing (mNGS) of the sputum samples collected from 50 pulmonary infected patients after cardiac surgery from December 2020 and June 2021 in Ruijin Hospital. Taxonomic classification of the sputum metagenomes showed that the numbers of species belonging to bacteria, fungi, and viruses were 682, 58, and 21, respectively. The full spectrum of microorganisms present in the sputum microbiome covered all the species identified by culture, including 12 bacterial species and two fungal species. Based on species-level microbiome profiling, a reference catalog of microbial abundance detection limits was constructed to assess the pathogenic risks of individual microorganisms in the specimens. The proposed screening procedure detected 64 bacterial pathogens, 10 fungal pathogens, and three viruses. In particular, certain opportunistic pathogenic strains can be distinguished from the colonizing ones in the individual specimens. Strain-level identification and phylogenetic analysis were further performed to decipher molecular epidemiological characteristics of four opportunistic etiologic agents, including Klebsiella pneumoniae, Corynebacterium striatum, Staphylococcus aureus, and Candida albicans. Our findings provide a novel metagenomic insight into precision diagnosis for clinically relevant microbes, especially for opportunistic pathogens in the clinical setting.
Introduction
As is well known, postoperative pulmonary infections are characterized by cough, phlegm, shortness of breath, chest pain, and temperature above 38 degrees (Conde and Lawrence, 2008). The etiological agents causing respiratory tract infections (RTIs) are usually related to a wide variety of pathogenic microbes (Cappelletty, 1998; Chowdhary et al., 2016; Kutter et al., 2018). For instance, a study on 426 patients with suspected lower RTIs has revealed that the most prevalent bacterial species isolated from the specimens (i.e. sputum, endotracheal secretion, and bronchial washing) were Pseudomonas aeruginosa and Haemophilus influenzae from Gram-negative organisms, and Streptococcus pneumonia and Staphylococcus aureus from Gram-positive (Khan et al., 2015). Several studies on respiratory infections have reported that Schizophyllum commune belonging to filamentous basidiomycetes is an emerging fungal pathogen, and its culture and species identification need a couple of weeks under specific culture media (Chowdhary et al., 2013; Chowdhary et al., 2016). Although traditional culture methods are still the gold standard for clinically microbiological testing, the new omics approaches particularly metagenomics have facilitated culture-free, faster detection for different types of pathogens (e.g. bacteria, viruses, and fungi) in complicated infectious diseases (Deurenberg et al., 2017; Ferone et al., 2020). Besides, many members of the human respiratory tract microbiota are colonizing opportunistic pathogens that are normally harmless and, when the microbial homeostasis becomes disrupted, can cause infections (Price et al., 2017).
Taxonomic profiling via metagenomic next-generation sequencing (mNGS) techniques has enabled higher sensitivity and resolution than does conventional culture approaches for better understanding of microbial compositions present in the sputum specimens (Human Microbiome Project Consortium, 2012; Finch et al., 2015; Yan et al., 2016; Ditz et al., 2020; Dicker et al., 2021b). According to 16S rRNA gene amplicon sequencing of sputum samples and clinically measured phenotypes of patients with bronchiectasis, a prospective observational study has revealed that the risk of exacerbation and long-term outcomes are associated with the reduced diversity of sputum microbiome, particularly dominated by bacterial genera Pseudomonas, Enterobacteriaceae, and Stenotrophomonas (Dicker et al., 2021b). Another study has also demonstrated that the dominance of Proteobacteria in the sputum microbiome is significantly associated with the neutrophil activation pathway and increased mortality in chronic obstructive pulmonary disease (Dicker et al., 2021a). Besides, Yan et al., have reported that the “near-complete” bacterial genomes could be reconstructed for the highest-abundant pathogens through sputum massive metagenomic sequencing (Yan et al., 2016).
To date, how to distinguish microbial colonization and/or infection is still a vital issue for pathogen identification under hospital settings. To address the above issue, the sputum microbiota of a cohort of postoperative patients with pulmonary infection were investigated based on the data of mNGS and culture testing. Metagenomic analysis was performed to reveal the taxonomic diversity of the sputum microbial community. Based on species-level microbiome profiling, a reference catalog of abundance detection limits was constructed to assess the pathogenic risk of individual microorganisms, which was further compared with the species detected by culture. Subsequently, strain-level identification and phylogenetic analysis were carried out to characterize the clinically relevant pathogens. The scheme of microbial abundance detection limits proposed in this manuscript may provide a useful reference to metagenomic surveillance for opportunistic pathogenic microbes.
Materials and Methods
Patients and Samples
This retrospective analysis was performed for the patients with suspected lung infections after cardiac surgery in Ruijin Hospital. In this study, the postoperative patients were enrolled based on abnormal imaging changes represented by chest X-ray, including chest infiltration, consolidation, cavity, and pleural effusion. Pulmonary infection was determined if at least two of the following three options were satisfied: 1) fever greater than 38℃, 2) the number of white blood cells beyond the routine detection limits, 3) occurrence of purulent airway secretion. Sputum specimens of the patients were collected from December 2020 and June 2021 (Table S1).
Specimen Culture
According to the cultivation procedure for bacteria and fungi in the diagnostic laboratory of Ruijin Hospital, the routine isolation media were used, including blood agar, eosin methylene-blue (EMB) agar, chocolate agar, Mueller-Hinton agar, sabouraud dextrose agar, and cooked meat medium (Oxoid, UK). Chocolate agar and EMB agar plates were incubated in 5-10% CO2 at 37°C for 48-72 hours; blood agar, Mueller-Hinton agar, and cooked meat medium plates were incubated at 37°C for 24-48 hours. Sabouraud dextrose agar plates suitable for the isolation of Candida albicans were incubated at 28°C and 35°C for 24-48 hours, respectively. Species identification of isolated strains was then carried out using the VITEK-2 Compact Instrument (bioMérieux, France) (Melhem et al., 2014).
Experiments of mNGS
For the preprocessing of specimens, sputum samples were homogenized by the dithiothreitol (DTT) treatment as proposed by Terranova et al. (Terranova et al., 2018). Approximately 200 μL sputum per sample was resuspended by 600 μL of DTT solution (1% concentration) and 20 μL of proteinase K (20 mg mL-1), followed by incubation at 56°C for 20 mins. Cells were removed by centrifugation to reduce host-background nucleic acid. Total DNA was extracted using HostZERO Microbial DNA Kit (Zymo Research, USA) following the manufacturer’s instructions. DNA extraction yield was quantified using a Quant-iT dsDNA HS Assay Kit and Qubit 3.0 Fluorometer (Thermo Scientific, USA). Enzymatic shearing was employed for the fragmentation (~200 bp) of DNA molecules and the libraries were then constructed using the Nextera XT DNA Library Preparation Kit (Illumina, USA). A certain number of PCR cycles (5-7) were applied to produce 1 µg of the metagenomic library. The quality of the libraries was assessed by a 2100 Bioanalyzer using the High Sensitivity DNA Assay (Agilent Technologies, USA). Metagenome shotgun sequencing in a single-end 75-bp mode was performed using the NextSeq 500/550 High Output Kit (92 cycles) on an Illumina NextSeq 550 sequencer. The samples as the No-Template Control (NTC) were sequenced simultaneously to assess contaminations during the wet-lab experiments.
Bioinformatics Analysis of Species-Level Abundance Profiling
Raw sequencing data were first subjected to a quality control process for trimming adapter sequences, removing low-quality tails and reads by Trimmomatic v0.36 (Bolger et al., 2014). Next, the reads mapping to the human reference genome GRCh37 were excluded using the short-read alignment tool Bowtie v2.2.6 (Langmead and Salzberg, 2012). To remove identical duplicated reads considered as technical artifacts by PCR (Martin et al., 2018), de-duplication was performed using an in-house Perl script. Taxonomic classification of microbial reads was conducted using Kraken v2.0.9-beta (Wood et al., 2019) and a custom k-mer database constructed using 51,543 genomes of ~27,000 species from the NCBI assembly databases (Kitts et al., 2016). To produce species-level abundance estimates, the number of reads in the Kraken classification report was further estimated by the Bayesian algorithm implemented by Bracken v2.2 (Lu et al., 2017). The microbial taxa bearing putative contamination during the experiments were removed from the final profiles based on the information of NTC samples in each batch: the ratio of read count (sputum sample)/read count (NTC) less than 10 for a specified species. Using a similar way to TPM (transcripts per million) described by Li et al. (Li et al., 2009), the estimates of percentage relative abundance of each species were computed for bacteria, viruses (excluding phages), and fungi, respectively.
To investigate pathogenic signals in the sputum microbiome, we proposed a rule of thumb to evaluate detection limits of pathogenic risks for individual organisms through species-level abundance profiling. Briefly, the values of the mean and standard deviation of relative abundances were calculated according to the estimates of each species across the samples. In this study, limits of detection (LOD) for prompting pathogenic potential of a specified species were estimated through the filtering criteria below: read count > 100; relative abundance > 0.05; presence in at least three samples; the mean plus the standard deviation of the abundance values as the LOD of the pathogenic risk.
Strain-Level Phylogenetic Analysis
Strain-level profiling analysis was conducted based on single-nucleotide variants (SNVs) using the StrainPhlAn v3.0 package (Truong et al., 2017). Briefly, metagenomic reads per sample were first aligned to the marker gene database of MetaPhlAn v3.0 (Beghini et al., 2021). The sample-specific consensus sequences from the markers were reconstructed for each species’ most abundant strain and the species-specific markers were extracted from MetaPhlAn for identifying the homologous sequences in the reference isolate genomes using blastn (Altschul et al., 1990). Based on sequence alignment, maximum-likelihood phylogenetic trees were reconstructed using RAxML (Stamatakis, 2014) and then visualized using MEGA10 (Kumar et al., 2018). Genome sequences of the isolated strains were retrieved from the NCBI Assembly database and were used for phylogenetic reconstruction together with the strains detected in the metagenomic samples. The program MetaMLST v1.2.2 was employed to predict the known and novel Sequence Type (ST) of the most abundant strains in metagenomic samples (Zolfo et al., 2017). Since it is complicated to resolve the phylogeny of the diploid fungi, an alternative pipeline was used to build the phylogenetic tree according to multilocus sequence typing (MLST) analysis of seven loci (Hirakawa et al., 2015b). Seven MLST allele sequences were reconstructed from the metagenomes using MetaMLST, and the related gene sequences from the isolate genomes were processed using FastMLST (Guerrero-Araya et al., 2021). Multiple sequence alignment was performed with MUSCLE (Edgar, 2004), and the phylogeny was inferred with RAxML.
Statistical Analysis of Community Diversity
The analysis of α-diversity and two sample T-test were conducted using the package fossil v0.4.0 (Chao, 1987) based on the Chao1 index representing community richness. To estimate the community compositional variation of the samples between groups, β-diversity was measured by calculating the Bray-Curtis distance matrix using Vegan v2.5-7 (Dixon, 2003). Differences in the rank dissimilarities within and between groups were inferred using the non-parametric ANOSIM test of 1,000 permutations. All statistical analyses and visualizations were performed in R v4.1.0 (R Core Team, 2021).
Results
Characteristics of Patients and Samples
In this retrospective study, we analyzed the data of microbiological testing on sputum samples from postoperative patients that were diagnosed with pulmonary infection. During the 7-month study, the specimens were collected from 50 adults aged from 34 to 81 years with a median of 64 years. Approximately three quarters (n = 36; 72%) were male. Participant characteristics are listed in Table S1, including body mass index (BMI), dates of both cardiac surgery and specimen sampling, patient outcomes, National Early Warning Score (NEWS) (Smith et al., 2013), and Sequential Organ Failure Assessment (SOFA) scores upon admission (Vincent et al., 1998). The median time of sample collection was day 3 after the cardiac surgery. The samples were then subject to conventional culture and mNGS testing for species identification. DNA mNGS testing yielded an average of 22.9 million reads for each specimen (Table S2).
Species Identification and Pathogenic Risk
As shown in Table S1, the culture testing reported 12 bacterial species and two fungal species across all the samples. It resulted in one species detected in each of the 16 sputum samples, two species in five samples, and three species in one sample (Table 1). None of the microorganisms was detected in the remaining tested samples (n = 28; 56%) by the culture method. On the contrary, a substantial number of microorganisms were identified by mNGS testing. The numbers of microorganisms assigned to the domains bacteria, fungi, and viruses were 682, 58, 21, respectively. Particularly, the numbers of bacterial and fungal species detected by mNGS testing were much higher than those by culture in the individual samples (Table 1). The average number of bacterial species in each sample was 90 with a range from 4 to 219. The average number of fungal species in each sample was 5 with a range from 1 to 17. Besides, the viral species were identified in more than half (n = 31; 62%) of the samples through mNGS; whereas the routine culture approach used in our study can support isolation and identification for bacteria and fungi, but not for viruses.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/855839/fcimb-12-855839-HTML/image_m/fcimb-12-855839-t001.jpg)
Table 1 Summary of the number of microorganisms detected by culture and mNGS testing for the individual sputum specimens.
To assess the pathogenic risk of individual microorganisms, we proposed the abundance LOD for screening candidate pathogens based on mNGS assays and species-level abundance profiles of bacteria, fungi, and viruses, respectively (Table S3). Most of the bacterial and fungal species detected by the culture approach were present in the list of the pathogens that were screened out in the individual samples (Table 2; the full list shown in Table S1). For instance, both culture and mNGS testing for patient #2 identified B. cenocepacia and C. albicans, which were predominant pathogens in the bacterial (86.3% abundance) and fungal community (98.7%), respectively. For patient #44, the mNGS approach identified five pathogens C. striatum (56.2%), S. aureus (39.2%), C. albicans (93.1%), A. flavus (5.5%), and Human alphaherpesvirus 1 (99.8%). Of these, S. aureus and C. albicans were also detected by culture. According to the LOD (30.8%) of K. pneumoniae, this bacterium was assessed as a pathogen in the four patients: 85.4% for #1, 87.1% for #4, 97.7% for #16, and 72.1% for #40. The sputum culture was positive for K. pneumoniae in these patients except for #4. The mNGS approach for patient #4 identified K. pneumoniae (87.1%) and C. albicans (97.8%) as pathogens; whereas, none of them was detected by culture. C. albicans, an opportunistic pathogenic yeast, was the most prevalent fungal pathogen detected in the sputum metagenomes of 15 patients (Table S1). In addition, human alphaherpesvirus 1 was the most frequently detected pathogen that was identified in six patients. Generally, according to our proposed mNGS screening criteria, the predicted microbial targets bearing pathogenic risk not only well matched the ones detected by culture, but also provided additional and prioritized pathogens missed by culture. It should provide complementary evidence for microbiological diagnosis in the clinical setting.
Community Structure of the Sputum Microbiota
To explore organismal composition in the sputum microbiota, the mNGS dataset containing all sequenced metagenomes was analyzed. We focused on the diversity of bacterial communities bearing high species richness compared to the fungal and viral communities (Table 1). The species accumulation curve tended to be smooth, indicating that the number of bacterial species would vary slowly with the increase in sample size (Figure 1A). The profile of percentage taxonomic abundances across the sputum metagenomes was displayed in Figure 1B. According to the average values of abundances, the top five abundant species were K. pneumonia (7.3%), S. maltophilia (5.6%), S. epidermidis (4.9%), S. infantis (4.8%), and A. baumannii (3.6%) (Figure 1B). It was observed that certain species that were remarkably dominating in the community of the individual samples were estimated to be etiologic agents mentioned above, such as K. pneumonia in #1/#4/#16/#40, A. baumannii in #25/#36, N. flavescens in #9/#17/#46, and S. anginosus in #13/#19/#21/#24/#26 (Table S1). Next, we compared the community structural dynamics of the sputum microbiota between two clinical outcomes: Recovery and Worse. The α-diversity analysis of community richness indicated that no significant difference (p = 0.125) was found between the two groups (Figure 1C). Based on the β-diversity ANOSIM test (R = 0.173, p = 0.008), it was apparent that the mean of ranked dissimilarities between groups was lower than that within the group Worse (Figure 1D). Meanwhile, the low R-value suggested a weak effect of clinical outcomes on the microbial communities. Both α- and β-diversity analyses indicated that the bacterial community diversity was less associated with clinical outcomes of postoperative infected patients, probably due to diverse pathogens dominating the individual samples.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/855839/fcimb-12-855839-HTML/image_m/fcimb-12-855839-g001.jpg)
Figure 1 Bacterial community structure and diversity of the sputum microbiota from the post-surgery patients with pulmonary infection. (A) An accumulation curve of the microorganisms in the sputum samples. The boxplot denotes the species richness according to the observed species numbers for a given number of samples collected from 100 random permutations of all the samples. (B) Percent stacked barplot for species-level taxonomic profiling of the sputum bacterial microbial communities. Labels stand for the top 12 abundant species. The samples are grouped by clinical outcomes Recover and Worse. (C) Chao1 index of the sputum microbiota from the two groups Recover and Worse. (D) The ANOSIM boxplot ranks between and within the Recover and Worse groups.
Strain-Level Characteristics of Pathogens
Strain identification is of significant importance for better understanding microbial pathogenicity and transmission of clinically relevant microorganisms (Nayfach et al., 2016). According to the pathogen list assessed by mNGS testing (Table S1), strain-level profile analysis was performed for four opportunistic pathogenic species K. pneumoniae, C. striatum, S. aureus, and C. albicans. StrainPhlAn analysis displays strain-level phylogenetic trees visualizing the relatedness between the metagenome-recovered strains and clinical bacterial isolates with publicly genome sequences (Figure 2). K. pneumoniae strains were reconstructed from the metagenomic samples of three patients #1/#4/#16 (Figure 2A). The analysis of metaMLST showed that the strain from patient #4 was K. pneumoniae of ST29 and the other two strains were represented by novel STs. Consistently, a K. pneumoniae strain derived from #4 was most closely related to a multidrug-resistant ST29 strain (CP024489/INF249) isolated from urinary tract infection (Gorrie et al., 2018) (Figure 2A). The strain from #16 was closely related to the K. pneumoniae strain NTUH-K2044 (AP006725) causing liver abscess and meningitis (Wu et al., 2009). In addition, two C. striatum strains from #27/#44 were phylogenetically placed together (Figure 2B), and their closest relatives were the two strains (VCOZ00000000 and VCOY00000000) of C. striatum isolated from influenza patients with secondary lower respiratory tract bacterial infections (Qin et al., 2020). Besides, the S. aureus strain from #44 was neighboring to the strain 55/2053, which caused the global pandemic of Panton-Valentine leukocidin-producing S. aureus in Europe and the United States in the 1950s (Zuo et al., 2021) (Figure 2C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/855839/fcimb-12-855839-HTML/image_m/fcimb-12-855839-g002.jpg)
Figure 2 Maximum-likelihood phylogenetic trees of three bacterial pathogens. Using StrainPhlAn trees are built based on the sample-specific consensus sequences of marker genes from K. pneumoniae (A), C. striatum (B), and S. aureus (C), respectively. The strains detected in the metagenomic samples are labeled in red. The isolate genomes are labeled by the GenBank accessions, strain identifiers, isolation country, and year. Bootstrap support values greater than 80% are displayed.
Figure 3 shows the phylogeny of three metagenomic strains of C. albicans from #13/#20/#29 together with the isolated strains of C. albicans, which is a well-known human fungal pathogen responsible for painful mucosal infections (Kim and Sudbery, 2011; Gulati and Nobile, 2016). Phylogenetic relationships of these fungal strains were investigated by the StrainPhlAn-based tree and the MLST-based tree, respectively. Both trees show that the C. albicans strain from #20 was closely related to the strains CHN1, Ca529L, A123, and A67.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/855839/fcimb-12-855839-HTML/image_m/fcimb-12-855839-g003.jpg)
Figure 3 Maximum-likelihood phylogenetic trees of the fungal pathogen C. albicans. Trees are built according to the StrainPhlAn analysis (A) and MLST analysis (B), respectively. The strain of C. albicans in the metagenome sample #13 is excluded from the MLST tree due to its low sequencing depth and failure to reconstruct the MLST loci-sequences.
Discussion
Rapid species identification followed by assessment of microbial colonization and pathogenesis associated with infections is always a significant issue for public health (Stressmann et al., 2017; Martin and Bachman, 2018). In recent years, culture-free high-throughput sequencing methods have successively uncovered many microbial colonizers in the nasopharynx, oropharynx, airway, trachea, and lung of healthy humans, such as Streptococcus, Haemophilus, Corynebacterium, Staphylococcus, Klebsiella, et al. (Charlson et al., 2011; Martin and Bachman, 2018; Paling et al., 2020). Microbiota balance of specific body sites, including the respiratory tract, is thought to be closely related to host health and disease (Man et al., 2017). Disruption of the microbiota homeostasis underlying the state of diseases enables changing the proportions of certain opportunistic pathogenic strains, thereby facilitating the proliferation of a complex community dominated by these strains and further causing infectious diseases (Haldar et al., 2020; Dicker et al., 2021b). Inference on the risk of colonization and infection by opportunistic pathogens has important implications for the diagnostics of these clinically relevant strains. Thus, we explored the screening procedure of mNGS-based species abundance profiles from the patients with pulmonary infection, further proposing LOD to assess the pathogenic risk in the sputum microbiome in this manuscript.
Microbiome-wide estimation on pathogenic risk for diverse microorganisms is beneficial for distinguishing the pathogenic strains from the colonizing strains in the clinical specimens. In this retrospective analysis, the final pathogen list generated by mNGS-based screening procedures encompassed nearly all the species identified by culture and some additional pathogens associated with infection. Of these pathogens detected in the sputum microbiome, many are human colonizers and also opportunistic etiologic agents, e.g. K. pneumoniae (Chung, 2016), S. marcescens (Khanna et al., 2013), S. oralis (Adams et al., 2017), S. mitis (Adams et al., 2017), C. striatum (Alibi et al., 2017), and S. aureus (Feng et al., 2008). Gram-positive bacterium S. pneumoniae is one of the prevalent commensal residents of the nasopharynx and also an opportunistic pathogen that can lead to pneumonia, bacteremia, meningitis, and otitis media (Shak et al., 2013). Herein, S. pneumoniae was present in more than half (n = 28) of the samples (Table S3). Based on the LOD of 3.6% for this bacterium, it was identified as a pathogen in the three patients: 15.7% abundance in #6; 7.4% in #33; 8.7% in #42. S. marcescens, a Gram-negative opportunistic nosocomial pathogen, is commonly involved in catheter-associated bacteremia, urinary tract infections, and wound infections of hospitalized patients (Khanna et al., 2013). S. marcescens was detected as a singular pathogen (56.8%) in #10, which was consistent with species identification by culture. C. albicans (36.3%) is normally a harmless commensal yeast and also a prevalent opportunistic pathogen that can cause symptomatic infections of mucosal membranes (Kim and Sudbery, 2011). This common fungus was present in about half (n = 23) of the samples (Table S3), and it was further assessed as a high-risk pathogen dominating the fungal community of 30% of the samples. These results suggest that the proposed LOD value of pathogenic risk may be a promising clinical indicator for infections caused by opportunistic etiological agents. However, it still needs a bigger sample group to provide robust references of microbial abundances for statistical estimation on LOD of pathogenic risk.
On the other hand, mining metagenomic data can provide more accurate taxonomic classification even at the single-strain resolution, which can be further used to interpret microbial pathogenesis, antimicrobial resistance, and transmission of the newly detected strain (Nayfach et al., 2016; Truong et al., 2017). For example, the sputum culture was negative for K. pneumoniae in patient #4, whereas K. pneumoniae was identified in the corresponding specimen by mNGS. Interestingly, strain profiling of the sputum metagenome revealed that this highly abundant K. pneumoniae strain was phylogenetically related to the multidrug-resistant and virulent isolate INF249 belonging to K. pneumoniae of ST29 and KL30 (Gorrie et al., 2018) (Figure 2A). Additionally, the K. pneumoniae strain detected in patient #16 was related to NTUH-K2044, a hypervirulent strain bearing the K1 capsule and ST23 (Wu et al., 2009). Besides, for the fungal typing of C. albicans, we identified three strains that were distantly related to each other in the phylogenetic tree (Figure 2D). The strain of C. albicans isolated from patient #13 was most closely related to C. albicans strain P57072, which is isolated from a bloodstream infection (Hirakawa et al., 2015). Since opportunistic etiological agents normally include both disease-causing and commensal strains (Paczosa and Mecsas, 2016; Braz et al., 2020), strain-level phylogeny is an alternative to providing accurate detection of whether it is a pathogenic strain or not.
In summary, we utilized mNGS-based species abundance profiling data to propose a practical screening strategy for the assessment of microbial pathogenic risk in sputum of pulmonary infected patients. The strategy enables reducing the disturbance from colonizers in a microbial community; meanwhile, it could also distinguish pathogenic strains from the colonizing strains in specific body sites. Our work provides a novel metagenomic insight into precision diagnosis for clinically relevant microbes, especially for opportunistic pathogens in the clinical setting.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository and accession number can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA765792.
Ethics Statement
The studies involving human participants were reviewed and approved by The Ethics Committee of Ruijin Hospital, Shanghai Jiaotong University School of Medicine Review Board (2021-69). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Conceived and designed the study: JC and JG. Performed the experiments: JC, XL, and QY. Analyzed the data: XC, LS, and KQ. Wrote and revised the manuscript: JC, LS, and JG. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (81671832).
Conflict of Interest
The authors LS and XC are employed by Genoxor Medical Science and Technology Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.855839/full#supplementary-material
References
Adams, H. M., Joyce, L. R., Guan, Z., Akins, R. L., Palmer, K. L. (2017). Streptococcus Mitis and S. Oralis Lack a Requirement for CdsA, the Enzyme Required for Synthesis of Major Membrane Phospholipids in Bacteria. Antimicrobial. Agents Chemother. 61 (5), e02552–e02516. doi: 10.1128/AAC.02552-16
Alibi, S., Ferjani, A., Boukadida, J., Cano, M. E., Fernández-Martínez, M., Martínez-Martínez, L., et al. (2017). Occurrence of Corynebacterium Striatum as an Emerging Antibiotic-Resistant Nosocomial Pathogen in a Tunisian Hospital. Sci. Rep. 7 (1), 9704. doi: 10.1038/s41598-017-10081-y
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., Lipman, D. J. (1990). Basic Local Alignment Search Tool. J. Mol. Biol. 215 (3), 403–410. doi: 10.1016/s0022-2836(05)80360-2
Beghini, F., McIver, L. J., Blanco-Míguez, A., Dubois, L., Asnicar, F., Maharjan, S., et al. (2021). Integrating Taxonomic, Functional, and Strain-Level Profiling of Diverse Microbial Communities With Biobakery 3. eLife 10, e65088. doi: 10.7554/eLife.65088
Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 30 (15), 2114–2120. doi: 10.1093/bioinformatics/btu170
Braz, V. S., Melchior, K., Moreira, C. G. (2020). Escherichia Coli as a Multifaceted Pathogenic and Versatile Bacterium. Front. Cell Infect. Microbiol. 10 (793). doi: 10.3389/fcimb.2020.548492
Cappelletty, D. (1998). Microbiology of Bacterial Respiratory Infections. Pediatr. Infect. Dis. J. 17 (8 Suppl), S55–S61. doi: 10.1097/00006454-199808001-00002
Chao, A. (1987). Estimating the Population Size for Capture-Recapture Data With Unequal Catchability. Biometrics 43 (4), 783–791. doi: 10.2307/2531532
Charlson, E. S., Bittinger, K., Haas, A. R., Fitzgerald, A. S., Frank, I., Yadav, A., et al. (2011). Topographical Continuity of Bacterial Populations in the Healthy Human Respiratory Tract. Am. J. Respir. Crit. Care Med. 184 (8), 957–963. doi: 10.1164/rccm.201104-0655OC
Chowdhary, A., Agarwal, K., Meis, J. F. (2016). Filamentous Fungi in Respiratory Infections. What Lies Beyond Aspergillosis and Mucormycosis? PloS Pathog. 12 (4), e1005491–e1005491. doi: 10.1371/journal.ppat.1005491
Chowdhary, A., Randhawa, H. S., Gaur, S. N., Agarwal, K., Kathuria, S., Roy, P., et al. (2013). Schizophyllum Commune as an Emerging Fungal Pathogen: A Review and Report of Two Cases. Mycoses 56 (1), 1–10. doi: 10.1111/j.1439-0507.2012.02190.x
Chung, P. Y. (2016). The Emerging Problems of Klebsiella Pneumoniae Infections: Carbapenem Resistance and Biofilm Formation. FEMS Microbiol. Lett. 363 (20). doi: 10.1093/femsle/fnw219
Deurenberg, R. H., Bathoorn, E., Chlebowicz, M. A., Couto, N., Ferdous, M., García-Cobos, S., et al. (2017). Application of Next Generation Sequencing in Clinical Microbiology and Infection Prevention. J. Biotechnol. 243, 16–24. doi: 10.1016/j.jbiotec.2016.12.022
Dicker, A. J., Huang, J. T. J., Lonergan, M., Keir, H. R., Fong, C. J., Tan, B., et al. (2021a). The Sputum Microbiome, Airway Inflammation, and Mortality in Chronic Obstructive Pulmonary Disease. J. Allergy Clin. Immunol. 147 (1), 158–167. doi: 10.1016/j.jaci.2020.02.040
Dicker, A. J., Lonergan, M., Keir, H. R., Smith, A. H., Pollock, J., Finch, S., et al. (2021b). The Sputum Microbiome and Clinical Outcomes in Patients With Bronchiectasis: A Prospective Observational Study. Lancet Respir. Med. 9 (8), 885–896. doi: 10.1016/s2213-2600(20)30557-9
Ditz, B., Christenson, S., Rossen, J., Brightling, C., Kerstjens, H. A. M., van den Berge, M., et al. (2020). Sputum Microbiome Profiling in COPD: Beyond Singular Pathogen Detection. Thorax 75 (4), 338–344. doi: 10.1136/thoraxjnl-2019-214168
Dixon, P. (2003). VEGAN, a Package of R Functions for Community Ecology. J. Vegetation. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Edgar, R. C. (2004). MUSCLE: Multiple Sequence Alignment With High Accuracy and High Throughput. Nucleic Acids Res. 32 (5), 1792–1797. doi: 10.1093/nar/gkh340
Feng, Y., Chen, C.-J., Su, L.-H., Hu, S., Yu, J., Chiu, C.-H. (2008). Evolution and Pathogenesis of Staphylococcus Aureus : Lessons Learned From Genotyping and Comparative Genomics. FEMS Microbiol. Rev. 32 (1), 23–37. doi: 10.1111/j.1574-6976.2007.00086.x
Ferone, M., Gowen, A., Fanning, S., Scannell, A. G. M. (2020). Microbial Detection and Identification Methods: Bench Top Assays to Omics Approaches. Compr. Rev. Food Sci. Food Saf. 19 (6), 3106–3129. doi: 10.1111/1541-4337.12618
Finch, S., McDonnell, M. J., Abo-Leyah, H., Aliberti, S., Chalmers, J. D. (2015). A Comprehensive Analysis of the Impact of Pseudomonas Aeruginosa Colonization on Prognosis in Adult Bronchiectasis. Ann. Am. Thorac. Soc. 12 (11), 1602–1611. doi: 10.1513/AnnalsATS.201506-333OC
Gorrie, C. L., Mirceta, M., Wick, R. R., Judd, L. M., Wyres, K. L., Thomson, N. R., et al. (2018). Antimicrobial-Resistant Klebsiella Pneumoniae Carriage and Infection in Specialized Geriatric Care Wards Linked to Acquisition in the Referring Hospital. Clin. Infect. Dis. an Off. Publ. Infect. Dis. Soc. America 67 (2), 161–170. doi: 10.1093/cid/ciy027
Guerrero-Araya, E., Muñoz, M., Rodríguez, C., Paredes-Sabja, D. (2021). FastMLST: A Multi-Core Tool for Multilocus Sequence Typing of Draft Genome Assemblies. Bioinf. Biol. Insights 15, 11779322211059238. doi: 10.1177/11779322211059238
Gulati, M., Nobile, C. J. (2016). Candida Albicans Biofilms: Development, Regulation, and Molecular Mechanisms. Microbes Infect. 18, 310–321. doi: 10.1016/j.micinf.2016.01.002
Haldar, K., George, L., Wang, Z., Mistry, V., Ramsheh, M. Y., Free, R. C., et al. (2020). The Sputum Microbiome is Distinct Between COPD and Health, Independent of Smoking History. Respir. Res. 21 (1), 183. doi: 10.1186/s12931-020-01448-3
Hirakawa, M. P., Martinez, D. A., Sakthikumar, S., Anderson, M. Z., Berlin, A., Gujja, S., et al. (2015). Genetic and Phenotypic Intra-Species Variation in Candida Albicans. Genome Res. 25 (3), 413–425. doi: 10.1101/gr.174623.114
Human Microbiome Project Consortium (2012). Structure, Function and Diversity of the Healthy Human Microbiome. Nature 486 (7402), 207–214. doi: 10.1038/nature11234
Khanna, A., Khanna, M., Aggarwal, A. (2013). Serratia Marcescens- a Rare Opportunistic Nosocomial Pathogen and Measures to Limit its Spread in Hospitalized Patients. J. Clin. Diagn. Res. 7 (2), 243–246. doi: 10.7860/JCDR/2013/5010.2737
Khan, S., Priti, S., Ankit, S. (2015). Bacteria Etiological Agents Causing Lower Respiratory Tract Infections and Their Resistance Patterns. Iran BioMed. J. 19 (4), 240–246. doi: 10.7508/ibj.2015.04.008
Kim, J., Sudbery, P. (2011). Candida Albicans, a Major Human Fungal Pathogen. J. Microbiol. (Seoul. Korea). 49 (2), 171–177. doi: 10.1007/s12275-011-1064-7
Kitts, P. A., Church, D. M., Thibaud-Nissen, F., Choi, J., Hem, V., Sapojnikov, V., et al. (2016). Assembly: A Resource for Assembled Genomes at NCBI. Nucleic Acids Res. 44, D73–D80. doi: 10.1093/nar/gkv1226
Kumar, S., Stecher, G., Li, M., Knyaz, C., Tamura, K. (2018). MEGA X: Molecular Evolutionary Genetics Analysis Across Computing Platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Kutter, J. S., Spronken, M. I., Fraaij, P. L., Fouchier, R. A. M., Herfst, S. (2018). Transmission Routes of Respiratory Viruses Among Humans. Curr. Opin. Virol. 28, 142–151. doi: 10.1016/j.coviro.2018.01.001
Langmead, B., Salzberg, S. L. (2012). Fast Gapped-Read Alignment With Bowtie 2. Nat. Methods 9 (4), 357–359. doi: 10.1038/nmeth.1923
Li, B., Ruotti, V., Stewart, R. M., Thomson, J. A., Dewey, C. N. (2009). RNA-Seq Gene Expression Estimation With Read Mapping Uncertainty. Bioinformatics 26 (4), 493–500. doi: 10.1093/bioinformatics/btp692
Lu, J., Breitwieser, F. P., Thielen, P., Salzberg, S. L. (2017). Bracken: Estimating Species Abundance in Metagenomics Data. PeerJ. Comput. Sci. 3, e104. doi: 10.7717/peerj-cs.104
Man, W. H., de Steenhuijsen Piters, W. A. A., Bogaert, D. (2017). The Microbiota of the Respiratory Tract: Gatekeeper to Respiratory Health. Nat. Rev. Microbiol. 15 (5), 259–270. doi: 10.1038/nrmicro.2017.14
Martin, R. M., Bachman, M. A. (2018). Colonization, Infection, and the Accessory Genome of Klebsiella Pneumoniae. Front. Cell Infect. Microbiol. 8. doi: 10.3389/fcimb.2018.00004
Martin, T. C., Visconti, A., Spector, T. D., Falchi, M. (2018). Conducting Metagenomic Studies in Microbiology and Clinical Research. Appl. Microbiol. Biotechnol. 102 (20), 8629–8646. doi: 10.1007/s00253-018-9209-9
Melhem, M. S. C., Bertoletti, A., Lucca, H. R. L., Silva, R. B. O., Meneghin, F. A., Szeszs, M. W. (2014). Use of the VITEK 2 System to Identify and Test the Antifungal Susceptibility of Clinically Relevant Yeast Species. Braz. J. Microbiol. 44 (4), 1257–1266. doi: 10.1590/S1517-83822014005000018
Nayfach, S., Rodriguez-Mueller, B., Garud, N., Pollard, K. S. (2016). An Integrated Metagenomics Pipeline for Strain Profiling Reveals Novel Patterns of Bacterial Transmission and Biogeography. Genome Res. 26 (11), 1612–1625. doi: 10.1101/gr.201863.115
Paczosa, M. K., Mecsas, J. (2016). Klebsiella Pneumoniae: Going on the Offense With a Strong Defense. Microbiol. Mol. Biol. Rev. 80 (3), 629–661. doi: 10.1128/MMBR.00078-15
Paling, F. P., Hazard, D., Bonten, M. J. M., Goossens, H., Jafri, H. S., Malhotra-Kumar, S., et al. (2020). Association of Staphylococcus Aureus Colonization and Pneumonia in the Intensive Care Unit. JAMA Network. Open 3 (9), e2012741–e2012741. doi: 10.1001/jamanetworkopen.2020.12741
Price, L. B., Hungate, B. A., Koch, B. J., Davis, G. S., Liu, C. M. (2017). Colonizing Opportunistic Pathogens (COPs): The Beasts in All of Us. PloS Pathog. 13 (8), e1006369. doi: 10.1371/journal.ppat.1006369
Qin, T., Geng, T., Zhou, H., Han, Y., Ren, H., Qiu, Z., et al. (2020). Super-Dominant Pathobiontic Bacteria in the Nasopharyngeal Microbiota as Causative Agents of Secondary Bacterial Infection in Influenza Patients. Emerging. Microbes Infect. 9, 605–615. doi: 10.1080/22221751.2020.1737578
R Core Team (2021) R: A Language and Environment for Statistical Computing R Foundation for Statistical Computing (Vienna, Austria). Available at: https://www.R-project.org/ (Accessed November 1, 2021).
Shak, J. R., Vidal, J. E., Klugman, K. P. (2013). Influence of Bacterial Interactions on Pneumococcal Colonization of the Nasopharynx. Trends Microbiol. 21 (3), 129–135. doi: 10.1016/j.tim.2012.11.005
Smith, G. B., Prytherch, D. R., Meredith, P., Schmidt, P. E., Featherstone, P. I. (2013). The Ability of the National Early Warning Score (NEWS) to Discriminate Patients at Risk of Early Cardiac Arrest, Unanticipated Intensive Care Unit Admission, and Death. Resuscitation 84 (4), 465–470. doi: 10.1016/j.resuscitation.2012.12.016
Stamatakis, A. (2014). RAxML Version 8: A Tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033
Stressmann, F. A., Couve-Deacon, E., Chainier, D., Chauhan, A., Wessel, A., Durand-Fontanier, S., et al. (2017). Comparative Analysis of Bacterial Community Composition and Structure in Clinically Symptomatic and Asymptomatic Central Venous Catheters. mSphere 2 (5), e00146–e00117. doi: 10.1128/mSphere.00146-17
Terranova, L., Oriano, M., Teri, A., Ruggiero, L., Tafuro, C., Marchisio, P., et al. (2018). How to Process Sputum Samples and Extract Bacterial DNA for Microbiota Analysis. Int. J. Mol. Sci. 19 (10), 3256. doi: 10.3390/ijms19103256
Truong, D. T., Tett, A., Pasolli, E., Huttenhower, C., Segata, N. (2017). Microbial Strain-Level Population Structure and Genetic Diversity From Metagenomes. Genome Res. 27, 626–638. doi: 10.1101/gr.216242.116
Vincent, J. L., de Mendonça, A., Cantraine, F., Moreno, R., Takala, J., Suter, P. M., et al. (1998). Use of the SOFA Score to Assess the Incidence of Organ Dysfunction/Failure in Intensive Care Units: Results of a Multicenter, Prospective Study. Working Group on "Sepsis-Related Problems" of the European Society of Intensive Care Medicine. Crit. Care Med. 26 (11), 1793–1800. doi: 10.1097/00003246-199811000-00016
Wood, D. E., Lu, J., Langmead, B. (2019). Improved Metagenomic Analysis With Kraken 2. Genome Biol. 20 (1), 257. doi: 10.1186/s13059-019-1891-0
Wu, K.-M., Li, L.-H., Yan, J.-J., Tsao, N., Liao, T.-L., Tsai, H.-C., et al. (2009). Genome Sequencing and Comparative Analysis of Klebsiella Pneumoniae NTUH-K2044, a Strain Causing Liver Abscess and Meningitis. J. Bacteriol. 191, 4492–4501. doi: 10.1128/JB.00315-09
Yan, Q., Cui, S., Chen, C., Li, S., Sha, S., Wan, X., et al. (2016). Metagenomic Analysis of Sputum Microbiome as a Tool Toward Culture-Independent Pathogen Detection of Patients With Ventilator-Associated Pneumonia. Am. J. Respir. Crit. Care Med. 194 (5), 636–639. doi: 10.1164/rccm.201601-0034LE
Zolfo, M., Tett, A., Jousson, O., Donati, C., Segata, N. (2017). MetaMLST: Multi-Locus Strain-Level Bacterial Typing From Metagenomic Samples. Nucleic Acids Res. 45, e7. doi: 10.1093/nar/gkw837
Keywords: pulmonary infection, sputum microbiome, metagenomic next-generation sequencing, pathogen risk, limit of detection, opportunistic pathogen, strain profiling
Citation: Chen J, Sun L, Liu X, Yu Q, Qin K, Cao X and Gu J (2022) Metagenomic Assessment of the Pathogenic Risk of Microorganisms in Sputum of Postoperative Patients With Pulmonary Infection. Front. Cell. Infect. Microbiol. 12:855839. doi: 10.3389/fcimb.2022.855839
Received: 16 January 2022; Accepted: 11 February 2022;
Published: 03 March 2022.
Edited by:
Costas C. Papagiannitsis, University of Thessaly, GreeceReviewed by:
Ibrahim Bitar, Charles University, CzechiaAna Elena Pérez-Cobas, Ramón y Cajal Institute for Health Research, Spain
Copyright © 2022 Chen, Sun, Liu, Yu, Qin, Cao and Gu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianwei Gu, gjw11428@rjh.com.cn
†These authors have contributed equally to this work