Identifying Risk Factors for Secondary Infection Post-SARS-CoV-2 Infection in Patients With Severe and Critical COVID-19

Emerging evidence has unveiled the secondary infection as one of the mortal causes of post-SARS-CoV-2 infection, but the factors related to secondary bacterial or fungi infection remains largely unexplored. We here systematically investigated the factors that might contribute to secondary infection. By clinical examination index analysis of patients, combined with the integrative analysis with RNA-seq analysis in the peripheral blood mononuclear cell isolated shortly from initial infection, this study showed that the antibiotic catabolic process and myeloid cell homeostasis were activated while the T-cell response were relatively repressed in those with the risk of secondary infection. Further monitoring analysis of immune cell and liver injury analysis showed that the risk of secondary infection was accompanied by severe lymphocytopenia at the intermediate and late stages and liver injury at the early stages of SARS-CoV-2. Moreover, the metagenomics analysis of bronchoalveolar lavage fluid and the microbial culture analysis, to some extent, showed that the severe pneumonia-related bacteria have already existed in the initial infection.

Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a kind of beta coronavirus that infects the host cell via its interaction with the angiotensin-converting enzyme–related carboxypeptidase (ACE2), has become the greatest scourge for humanity since the 1918 influenza (1, 2). Unlike other viruses caused severe acute respiratory syndrome, as SARS-CoV-2’s targeted receptor ACE2 is mainly expressed in the epithelial cells of the cardiopulmonary, SARS-CoV-2 infection could lead to dysfunction of alveolar epithelial cells and acute inflammation response, exerting severe acute respiratory syndrome (3). Compared to the SARS and the Middle East Respiratory Syndrome Coronavirus (MERS), the virulence of the SARS-CoV-2 is relatively low, only a small portion of infected people developed into severe and critical syndromes (4, 5). When taking the extremely highly spreading ability of this virus into consideration, in particular the discovery of the new SARS-CoV-2 variant “ (VoC) 202012/01” in the UK and “B.1.617” in India, a greater number of secondary infections and deaths is inevitably accompanied by the enhanced version of the virus (6, 7). Thus, a better understanding of the causes of severe/critical secondary infections and death is urgently needed.

Many efforts have been devoted to unveiling the life-threatening causes of COVID-19, one of the major causes is the secondary infection in those severe or critical patients (8–13). Secondary infection, also known as co-infection, belongs to one of the leading causes of virus-related mortalities, especially the respiratory infections such as viral pneumonia (14, 15). Generally, there are many opportunistic pathogens existing in human body, such as Mycoplasma pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Haemophilus influenzae, Staphylococcus aureus, Streptococcus pneumoniae and other resistant Enterobacteriaceae (16–18), with the changes of the immune system and micro-environment when fighting against the viral infection, these opportunists look to expand rapidly and cause severe secondary infection post-SARS-CoV-2-infection (19, 20).

In SARS-CoV-2, many secondary infection events have been observed, including bacterial and fungal infection, especially in those severe/critical cases (9–13). The overall secondary infection rate post-SARS-CoV-2 infection is about 3.2% ~ 15%, and the mortality of secondary infection was much higher in those severe and critical COVID-19 cases (8, 9, 21). Still, what caused the secondary pathogens infection remains largely unexplored, and a better understanding of the potential mechanisms underlying secondary infection risk could help us predict the risk of severe disease and take action before it happens.

In this study, we systematically determined the factors contributed to elevated risk factors of secondary SARS-CoV-2 infection in single-center cases from Shanghai Public Health Clinical Center Affiliated to Fudan University. We performed RNA-seq analysis in 43 peripheral blood mononuclear cells (PBMCs) and 5 bronchoalveolar lavage (BAL) samples from COVID-19 patients at initial diagnosis stage. Also, we investigated the correlation between lymphocytopenia, multi-organ damage, coagulation and secondary bacterial infection of SARS-CoV-2. This study will provide prospective evidence of a rational and scientific therapeutic response to the SARS-CoV-2 pandemic and secondary co-infection.

Materials and Methods

Subjects and Data Collection

This study was retrospectively conducted in Shanghai Public Health Clinical Center, designated as a designated hospital for COVID-19 treatment in East China, from February to May 2020. Among all the patients, their symptoms and groups were identified strictly in accordance with the novel Coronavirus infection diagnosis and treatment plan (Trial version 7) published by Chinese General Office of the National Health Commission. A total of 458 patients were included in this study. Secondary infection was diagnosed based on the patients’ clinical, radiological and laboratory data. A total of 43 patients including 21 mild/moderate and 22 severe/critical were studied within five days post-symptoms. All the clinical data were recorded in a computerized database in the hospital laboratory medicine department.

Ethics Statement

This study was approved by the Ethics Committee of Shanghai Public Health Clinical Center (#2020-Y025-01), and informed written consents from all human subjects were signed according to the Declaration of Helsinki. This study was conducted according to the guideline of Novel Coronavirus Laboratory Safety from the General Office of the National Health Commission letter in 2020 (Document No. 70, 2nd edition).

PBMCs Isolation and BAL Collection

The PBMCs were separated by the Ficoll density gradient centrifugation using the whole peripheral blood samples from remaining samples for clinical examination as previous processing (22). About 4 mL of fresh collected peripheral blood samples in anticoagulant tubes were resuspended with equal volume phosphate buffer saline (PBS), the diluted samples were then gently transferred into a 15 mL tube contains 4mL of Ficoll-paque PLUS (GE Healthcare Life Sciences). The tubes were then centrifugated for 30 min at 400 g, and the PBMCs in the middle layer of the tubes were collected and transferred to a new 15 mL tube (BD Falcon, USA). Then PBMCs were washed twice with PBS and lysed in 1 mL Trizol reagent (Thermo Fisher Scientific, USA). Pulmonary bronchoalveolar lavage fluid was collected and centrifuged for 10 min at 600 g. Then pellets were washed once with PBS and lysed in 1 mL Trizol reagent.

RNA-seq and RNA-seq Analysis

RNAs were extracted by the Trizol reagent according to the manufacturer’s instructions. The mRNA-seq library was constructed according to the manufacturer’s instructions of the VAHTSTM mRNA-seq V3 Library Prep Kit for Illumina^®. The libraries were sequenced in the Illumina’s Nova-seq 6000 platform.

The RNA-seq analysis was performed as previously described (23). All reads were mapped to the Homo sapiens (human) genome assembly GRCh38 (hg38) by the STAR software (24). The read counts were extracted by HTseq algorithm (25). Differential expressed genes were analyzed by DESeq2 (26), and genes with two-fold change and adjusted P-value < 0.01 were kept.

For an integrative analysis, the expression matrix of PBMC cells in patient with Flu, COVID-19 and healthy donors were obtained from the GEO database under the accession of GSE150728 (27) and GSE149689 (28). The datasets were processed with the Seurat suite (29) and the Harmony algorithm (30).

Clinical Laboratory Examination and Data Collection

All laboratory tests were conducted in the department of laboratory medicine in the Shanghai Public Health Clinical Center. The blood routine tests, including white blood cell count (WBC), lymphocyte count (LPC), platelet count (PLT), Hypersensitive C-reactive protein (hs-CRP), were analyzed by Blood cell analyzer SYSMEX XN-A100 (Hisense Meikang medical electronics, Shanghai Co., Ltd). The CD3⁺ cell counts, CD4⁺ cell counts, CD8⁺ cell counts, CD19⁺ cell counts, CD16⁺CD56⁺ cell counts, and CD4⁺/CD8⁺ percentage were detected by FACS Canto II (BD Biosciences). Lactate (LACT), lactate dehydrogenase (LDH), Cholinesterase (CHE), Prealbumin (PA), Albumin (ALB) were detected by the Biochemical Immune Automation Analysis Workstation (ARCHITECT 3600J, Abbott Laboratories Co., LTD).

Microbiology Laboratory Testing

The sputum and bronchoalveolar lavage were used for microbiology laboratory tests. The collected samples were cultured in Columbia Blood Agar plate and chocolate Agar plate (Shanghai Comarja Company) at 35°C, 5% CO₂, and smear and gram staining were performed simultaneously. The Microflex MALDI-TOF MS mass spectrometry (Brock, Germany) and automatic microbiological analysis system (Biomerieux VITEK231, France) were used for the authentication of the strains of bacteria.

Quantification and Statistical Analysis

All data in this study were expressed as Mean ± SD as indicated. For parametric analysis, the F test was used to determine the equality of variances between the groups compared; statistical significance across two groups was tested by Student’s t-test, one-way analysis of variance (ANOVA) were used to determine statistically significant differences between multiple groups. Values of P < 0.05 were considered to be statistically significant.

Results

Demographics, Disease Course, and Outcome of Patients With COVID-19

We studied 458 patients from Shanghai Public Health Clinical Center during January to October 2020. A total of 21 mild/moderate and 22 severe/critical COVID-19 patients were randomly included in this study (Table 1). The average age of enrolled patients of the mild/moderate was 59.52 ± 14.36 years, 65.47 ± 13.45 for the severe/critical. Length of hospital stay was 35.23 ± 12.74 VS 64.28 ± 22.63. Patients with mild/moderate disease did not resort to respiratory equipment during treatment. In the severe/critical patients, 6 (27.27%) cases were diagnosed with secondary infection. Four patients of severe/critical disease were deceased, all of whom had secondary infection.

TABLE 1

Table 1 Demographics of patients.

Abnormal Anti-Viral Response in Severe and Critical SARS-CoV-2

To explore the immune response post-SARS-CoV-2 infection, we collected the PBMCs from 21 mild/moderate and 22 severe/critical COVID-19 patients within five days post-symptoms and performed RNA-seq experiments with these collected samples. Gene expression pattern was analyzed to determine the range associated with the disease. Among these 22 severe/critical cases, six developed secondary infection post-SARS-CoV-2 infection. The principal component analysis showing the mild/moderate and severe/critical could be separated according to the expression patterns, indicating the mild/moderate and severe/critical COVID-19 have distinct gene expression signatures (Figure 1A). To disclose the distinct transcriptional regulation programs in mild/moderate and severe/critical cases, we analyzed the differentially expressed genes in PBMC cells between mild/moderate and severe/critical samples. A total of 523 genes were significantly highly expressed in mild/moderate and 1,341 genes were significantly highly expressed in severe/critical COVID-19 (Figure 1B and Table S1).

FIGURE 1

Figure 1 RNA-seq analysis showing the transcriptional programs in blood cells related to the severity of COVID-19. (A) PCA analyses of the transcriptome of PBMC cells in mild/moderate and severe/critical COVID-19. (B) Heatmap showing the differentially expressed genes in PBMC cells between mild/moderate and severe/critical COVID-19. (C, D) Functional annotation of genes correlated with the risk of secondary infection during the initial infection.

Functional annotation of these genes demonstrated that genes related to T cell activation, regulation of the developmental process, and chemotaxis were enriched in those genes highly expressed in mild/moderate cases, indicating an active lymphocytic response to viral infections in these mild/moderate COVID-19 patients (Figure 1C and Table S2). Unlike those mild/moderate cases, we noticed that the severe/critical COVID-19 significantly enriched at genes related to neutrophil degranulation, myeloid leukocyte mediated immunity, and exocytosis (Figure 1D and Table S2), implying a different response to viral infection in these severe/critical patients as compared to those mild/moderate cases.

Distinct Expression Signatures of the Risk of Secondary Infection

As the secondary infection is one of the leading causes of death, we aim to find some predictive factors associated the potential secondary infection after the initial SARS-CoV-2 infection for the 22 severe/critical cases. We noticed that the transcriptome pattern of cases with the risk of secondary infection was significantly different from those without these risks in the samples of initial infection (Figure 2A), highlighting the potential risk had occurred during the initial infection of different individuals. To explore the risk of secondary infection at the initial SARS-CoV-2 infection, we obtained the differentially expressed genes in those cases with or without secondary infection risk. A total of 291 risk-related genes and 11 non-risk related genes were identified (Figure 2B and Table S3).

FIGURE 2

Figure 2 Identifying the key transcriptional programs that might enhance the risk of secondary infection at the initial infection of SARS-COV-2. (A) PCA analyses of the severe cases with or without the risk of secondary infection. (B) Heatmap showing the differentially expressed genes in PBMC cells between risk and non-risk of secondary infection. (C, D) Functional annotation of genes correlated with risk of secondary infection during the initial infection. (E, F) scRNA-seq analyses showing the enriched of these risk genes in distinct immune cells from PBMC.

Functional annotation analysis of these genes showed that genes associated with T cell activation and lymphocyte co-stimulation were enriched within the non-risk related genes (Figure 2C), implying the functional T cell response might protect the severe/critical COVID-19 from secondary infection. Notably, genes involved in the regulation of myeloid cell homeostasis, oxygen transport, erythrocyte differentiation, and antibiotic catabolic process were correlated with an enhanced risk of secondary infection (Figure 2D and Table S4). Also, the activation of the antibiotic catabolic process was observed in these patients with secondary infection risk. These findings suggest that further monitoring of antibiotic metabolism might be required.

To further explore the cells in the PBMCs contributed to these transcriptional regulatory programs, we retrieved two independent RNA-seq datasets in PBMCs isolated from COVID-19, Flu, and healthy donors (Figure S1). And then, we analyzed the enrichment of these genes related to erythrocyte homeostasis, hydrogen peroxide catabolic process, myeloid cell homeostasis, tetrapyrrole metabolic process, myeloid cell development, antibiotic catabolic process, and heme metabolic process in distinct subtypes of cells in PBMCs (Figures 2E, F). We noticed that the majority of these processes were enriched in the erythrocyte cells and macrophage cells, indicating the secondary infection risk related to dysregulation of the transcriptional programs was mainly occurring within the erythrocyte cells and macrophage cells.

Secondary Infection Risk Was Accompanied by Lymphocytopenia

After the initial infection, the immune response was dynamically changed at different stages. We next traced the changes of immune cells in the initial infection stage, intermediate or later stage, and the potential convalescent stages and analyzed the distinct patterns in those patients with secondary infection risks. We found that CD4⁺, CD8⁺, CD19⁺, CD3⁺, and CD16⁺CD56⁺ cells in severe and critical cases were significantly lower than those in mild or moderate cases. Moreover, the levels of these cells were even significantly lower in cases with secondary infection risk as compared to those non-secondary risk severe and critical cases at the intermediate or later stage of COVID-19 infection (Figures 3A–F), indicating the secondary infection risk were accomplished with dysfunction of lymphoid cell response, which could also be supported by the genomic response of patients from RNA-seq analyses in Figure 2B.

FIGURE 3

Figure 3 Elevation of lymphocytopenia enhanced the risk of secondary infection. The levels of CD4+ (A), CD8+ (B), CD3+ (C), CD19+ (D), CD16+CD56+ (E) and CD4+/CD8+ (F) ratio in mild/moderate COVID-19, severe/critical COVID-19 without risk of secondary infection and severe/critical COVID-19 with the risk of secondary infection were shown. Each group of samples was further divided into three groups according to the time post-symptoms, including day 0 to day 9, day 10 to day 30, and post day 30. The barplot patients with the risk of secondary infection at day 0 to day 9 were not shown due to no available data. The chosen level of statistically significance were *P < 0.01 and **P < 0.05, respectively.

Liver Injury Is Associated With the Risk of Secondary Infection

Next, we explored the biochemical features in those secondary infection risk and non-risk cases. We found that liver injury-related factors were highly correlated with secondary infection. The levels of lactate, lactate dehydrogenase was significantly higher in severe/critical cases as compared to those mild/moderate cases, and the levels of these two factors were even higher in the cases with secondary infection risk at the initial infection stage (Day 0 to Day 9) (Figures 4A, B). Meanwhile, the levels of cholinesterase, pre-albumin, and albumin were significantly decreased in severe and critical COVID-19, and the levels of these three factors were even lower in the severe and critical cases with secondary infection risk at both the initial infection stage and intermediate or later stage (Figures 4C–F). Factors associated with kidney injury were not correlated with secondary infection risk (Data not shown).

FIGURE 4

Figure 4 Elevation of Liver injury enhanced the risk of secondary infection. The levels of Lactate (A), Lactate dehydrogenase (B), Cholinesterase (C), Prealbumin (D), Albumin (E), and CO₂ (F) in mild/moderate COVID-19, severe/critical COVID-19 without risk of secondary infection and severe/critical COVID-19 with the risk of secondary infection were shown. Each group of samples was further divided into three groups according to the time post-symptoms, including day 0 to day 9, day 10 to day 30, and post day 30. The chosen level of statistically significance were *P < 0.01 and **P < 0.05, respectively.

Initial Microbial Environment Affected the Risk of Secondary Bacterial or Fungal Infection

Next, we explored the correlation between the initial microbial environment and the secondary infection risks. We performed RNA-seq analyses in five bronchoalveolar lavage fluid (BLF) samples isolated shortly after the initial infection according to a clinical report from the Microbiology laboratory, including three developed into a secondary bacterial infection, and analyzed the metagenomic features in these samples (Table S5). Many bacteria related to the risk of severe pneumonia were found during the initial infection, including Klebsiella pneumonia, Stenotrophomonas maltophilia, Acinetobacter pittii PHEA-2, Acinetobacter nosocomialis and Ralstonia pickettii were found in the initial infected samples of COVID-19 patients who developed the secondary bacterial infection, despite these bacteria were not the major bacteria in its environment (Figure 5).

FIGURE 5

Figure 5 Metagenomic analysis of BALF samples in the initial infection. The top 10 enriched bacteria or fungi were plotted. RNA-seq was performed with the BALF sample in patients with severe/critical COVID-19. All of these patients did not undergo secondary infection yet. All pneumonia-related bacteria were highlighted in red color. Three cases with secondary infection and two cases without secondary infection were included in this analysis.

Further analysis about bacteria culture results could support these findings (Table 2). Of these SARS-CoV-2 patients with co-pathogens infection, the pulmonary infection was the main infection. A few co-pathogens infections occurred in the period after admission, the majority of cases in this study showed sputum/BLF culture-negative for pathogenic bacteria in the initial infection, but the bacteria infection became culture-positive post-SARS-CoV-2 infection for quite a long time. Through the microbial culture results, it is speculated that, to some extent, the secondary infection in some patients may be their own endogenous infection, because our sequencing results in cases 4, 7, 8 were matched with their positive results in the later culture. The dysregulated immune micro-environment, micro-ecology and liver injury provided the chance for expanding the conditional pathogenic bacteria and led to a severe secondary infection.

TABLE 2

Table 2 General microbiological test information of patients with secondary infection post-SARS-COV-2 infection.

Discussion

In the face of the COVID-19, currently rampaging around the world and continues to spread, the secondary infection has been reported to be an adverse prognostic marker for the severe illness and mortality of SARS-CoV-2 (8, 31). It was documented than the mortality rate of secondary bacterial infections among COVID-19 patients admitted to the intensive care unit (ICU) was as high as 95% (32). Identifying potential predictive markers for secondary infection could help a better precaution, but the risk factors but the risk factors contributing to secondary infection are unclear. There has previously documented that the mechanism of the virus‐bacteria co-infection mainly involves the lack of effective immune response and pathogenic synergism by different factors (33, 34). We here reported that genes related to myeloid cell homeostasis, antibiotic catabolic process, and erythrocyte homeostasis pathways were dysregulated in severe and critical SARS-CoV-2 with the potential risk of secondary infection. Furthermore, we demonstrated that the liver injury-related factors at the early stage of infections, lymphocytopenia risk-related factors at the intermediate and late stages, and severe/critical SARS-CoV-2 with secondary infection risks are more abnormal compared to severe/critical cases with non-secondary infection risks.

Noteworthy, viral-caused pneumonia, such as influenza virus and coronavirus, usually led to an exhausted and fragile immune system due to the decreased immune cells when it fights against the virus (35, 36). Thus, viral infection suppresses certain changes in the body’s immune system, and secondary infection can be a life-threatening factor, especially in virus-induced pneumonia (37). In accordance with previous observation, we observed genes correlated with T cell activation were less activated in those severe/critical COVID-19 compared to the mild/critical patients (3, 38, 39). Moreover, the T cell activation and lymphocyte activation programs in the severe/critical cases with secondary infection risks were less activated than those in non-secondary-infection-risk groups of severe/critical COVID-19. Thus, the initial immunities of changed situation were different, and the T cell immunity was abnormal in severe/critical cases, especially those with potential secondary infection risks (22). These observations could be further verified by the flow cytometry analyses of the subset of lymphocytes, and we found that the CD4⁺ and CD8⁺ T cells were extremely low in patients with risk of secondary infections. In our previous work, we found that the frequencies of total lymphocytes and T cells decreased significantly in acutely infected COVID-19 patients compared to healthy controls (40).

In the transcriptomic analysis in PBMC samples, we found that the antibiotic catabolic process was activated in patients with risks of secondary infection. Since antibiotics are the major therapeutic approaches to bacterial infection, the elevated levels of antibiotics metabolism might lower the levels of antibiotics in the circulation system, which might increase the risk of bacterial and fungal co-infection (41, 42). Still, these speculations should be further validated by monitoring the plasma concentration of the antibiotics post-therapy.

The liver is generally considered as an organ involved in the immune response (43). The liver injury might partially lead to immune system failure in response to viral and bacterial infections (44). We noted that the liver injury-related markers, including decreased pre-albumin, albumin, and cholinesterase levels, were observed at the initial stages, and the intermediate stage of SARS-CoV-2 infection in severe and critical COVID-19. Moreover, the albumin and cholinesterase levels were even lower in those cases with the risk of secondary infection, highlighting a link between liver injury at the initial infection stage and the secondary infections. Many studies in COVID-19 have reported liver injury in severe and critical cases (9, 45). We here noticed that the liver injury happened at the initial stages of COVID-19, and thus it might be used as an indicator for the risk of secondary infections in severe and critical cases.

Similar to other viral infections, emerging evidence suggests that an increasing number of COVID-19 patients are being diagnosed with bacterial co-infections during hospitalization (32, 46). Studies about secondary infection to influenza infection have demonstrated that neuraminidase (NA) of influenza virus can activate transforming growth factor-β (TGF-β) in the host, resulting in increased expression of adhesion molecules on the surface of host cells. It increases bacterial adherence to host lung cells, which causes secondary bacterial infections and predisposes the host to secondary bacterial pneumonia (47). The source and specific nature of these pathogens are yet to be fully explored, but there is some evidence suggesting that multidrug‐resistant bacteria are thought to be responsible for the development of these infections. In one study conducted in Iran, the incidence of A. baumannii and S. aureus were the most common in the ICU (32). However, given lack of the data about SARS-CoV-2 secondary infection, we investigated the bacterial population in the pulmonary microenvironment yielded the bacterial class of the initial infection, and by comparison with the results of microbial culture after infection, it was found that to a certain extent, the culture-positive bacteria in the infection of some patients may be their own endogenous bacteria, of course, our speculation was based on the evaluation of patients with secondary infection. Further work is required to investigate whether there are other factors associated co-infected with COVID-19.

In summary, increasing co-infection is a non-negligible factor in COVID-19, our study identified some abnormalities in the early stage of SARS-CoV-2 infection that might contribute to secondary bacterial and fungal infections. We emphasized the concern of bacterial infections in COVID-19 patients, this could help us better to prevent worse outcomes due to secondary infection. Also, we provided some predictors for the pre-treatment of antibiotics in those risks of secondary infection cases. Still, there are some limitations to our study, with a small sample size, we only evaluated the peripheral blood immune response and clinical manifestations of patients in different groups, other studies with an expanded sample size might be conducted to address these issues further.

Data AvailabilityStatement

The data presented in the study are deposited in the GEO database repository, accession number GSE184401.

Ethics Statement

The studies involving human participants were reviewed and approved by The Ethics Committee of Shanghai Public Health Clinical Center Affiliated to Fudan University (#2020-Y025-01). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

ZZ, MQG, and JC designed the experiments. MG, JG, and MQG performed the majority of experiments. MG and JC performed the bioinformatics analyses. TZ, XJ, and JF collected the clinical data. QYW, XL performed the laboratory detection. MQG wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Natural Science Foundation of China (NSFC, 8210061470), the Emergency Project of Shanghai Science and Technology Commission (No. 20411950502), the Shanghai Sailing Program (No. 21YF1438800) and Project of Shanghai Public Health Clinical Center (No. KY-GW-2020-15).

Conflict of Interest

The 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.

Acknowledgments

We also would like to thank Dr. Yun Tan in National Research Center for Translational Medicine at Shanghai for his guidance and assistance in bioinformatics in this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.715023/full#supplementary-material

References

1. Zhang X, Tan Y, Ling Y, Lu G, Liu F, Yi Z, et al. Viral and Host Factors Related to the Clinical Outcome of COVID-19. Nature (2020) 583(7816):437–40. doi: 10.1038/s41586-020-2355-0

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Long QX, Liu BZ, Deng HJ, Wu GC, Deng K, Chen YK, et al. Antibody Responses to SARS-CoV-2 in Patients With COVID-19. Nat Med (2020) 26(6):845–8. doi: 10.1038/s41591-020-0897-1

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Shi Y, Wang Y, Shao C, Huang J, Gan J, Huang X, et al. COVID-19 Infection: The Perspectives on Immune Responses. Cell Death Differ (2020) 27(5):1451–4. doi: 10.1038/s41418-020-0530-3

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Petersen E, Koopmans M, Go U, Hamer DH, Petrosillo N, Castelli F, et al. Comparing SARS-CoV-2 With SARS-CoV and Influenza Pandemics. Lancet Infect Dis (2020) 20(9):e238–e44. doi: 10.1016/S1473-3099(20)30484-9

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Suleyman G, Fadel RA, Malette KM, Hammond C, Abdulla H, Entz A, et al. Clinical Characteristics and Morbidity Associated With Coronavirus Disease 2019 in a Series of Patients in Metropolitan Detroit. JAMA Netw Open (2020) 3(6):e2012270. doi: 10.1001/jamanetworkopen.2020.12270

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Mahase E. Covid-19: What Have We Learnt About the New Variant in the UK? BMJ (Clinical Res ed) (2020) 371:m4944. doi: 10.1136/bmj.m4944

CrossRef Full Text | Google Scholar

7. Yadav PD, Sapkal GN, Abraham P, Ella R, Deshpande G, Patil DY, et al. Neutralization of Variant Under Investigation B.1.617 With Sera of BBV152 Vaccinees. Clin Infect Dis (2021) ciab411. doi: 10.1093/cid/ciab411

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Langford BJ, So M, Raybardhan S, Leung V, Westwood D, MacFadden DR, et al. Bacterial Co-Infection and Secondary Infection in Patients With COVID-19: A Living Rapid Review and Meta-Analysis. Clin Microbiol Infect (2020) 26(12):1622–9. doi: 10.1016/j.cmi.2020.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhang H, Zhang Y, Wu J, Li Y, Zhou X, Li X, et al. Risks and Features of Secondary Infections in Severe and Critical Ill COVID-19 Patients. Emerg Microbes Infect (2020) 9(1):1958–64. doi: 10.1080/22221751.2020.1812437

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Mirzaei R, Goodarzi P, Asadi M, Soltani A, Aljanabi HAA, Jeda AS, et al. Bacterial Co-Infections With SARS-CoV-2. IUBMB Life (2020) 72(10):2097–111. doi: 10.1002/iub.2356

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lansbury L, Lim B, Baskaran V, Lim WS. Co-Infections in People With COVID-19: A Systematic Review and Meta-Analysis. J Infect (2020) 81(2):266–75. doi: 10.1016/j.jinf.2020.05.046

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Rawson TM, Moore LSP, Zhu N, Ranganathan N, Skolimowska K, Gilchrist M, et al. Bacterial and Fungal Co-Infection in Individuals With Coronavirus: A Rapid Review to Support COVID-19 Antimicrobial Prescribing. Clin Infect Dis (2020) 71(9):2459–68. doi: 10.1093/cid/ciaa530

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Bengoechea JA, Bamford CG. SARS-CoV-2, Bacterial Co-Infections, and AMR: The Deadly Trio in COVID-19? EMBO Mol Med (2020) 12(7):e12560. doi: 10.15252/emmm.202012560

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Morens DM, Taubenberger JK, Fauci AS. Predominant Role of Bacterial Pneumonia as a Cause of Death in Pandemic Influenza: Implications for Pandemic Influenza Preparedness. J Infect Dis (2008) 198(7):962–70. doi: 10.1086/591708

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Rice TW, Rubinson L, Uyeki TM, Vaughn FL, John BB, Miller RR 3rd, et al. Critical Illness From 2009 Pandemic Influenza A Virus and Bacterial Coinfection in the United States. Crit Care Med (2012) 40(5):1487–98. doi: 10.1097/CCM.0b013e3182416f23

PubMed Abstract | CrossRef Full Text | Google Scholar

16. MacIntyre CR, Chughtai AA, Barnes M, Ridda I, Seale H, Toms R, et al. The Role of Pneumonia and Secondary Bacterial Infection in Fatal and Serious Outcomes of Pandemic Influenza a(H1N1)Pdm09. BMC Infect Dis (2018) 18(1):637. doi: 10.1186/s12879-018-3548-0

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Song SH, Chen TL, Deng LP, Zhang YX, Mo PZ, Gao SC, et al. Clinical Characteristics of Four Cancer Patients With SARS-CoV-2 Infection in Wuhan, China. Infect Dis Poverty (2020) 9(1):82. doi: 10.1186/s40249-020-00707-1

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Stawicki SP, Jeanmonod R, Miller AC, Paladino L, Gaieski DF, Yaffee AQ, et al. The 2019-2020 Novel Coronavirus (Severe Acute Respiratory Syndrome Coronavirus 2) Pandemic: A Joint American College of Academic International Medicine-World Academic Council of Emergency Medicine Multidisciplinary COVID-19 Working Group Consensus Paper. J Global Infect Dis (2020) 12(2):47–93. doi: 10.4103/jgid.jgid_86_20

CrossRef Full Text | Google Scholar

19. Brundage JF. Interactions Between Influenza and Bacterial Respiratory Pathogens: Implications for Pandemic Preparedness. Lancet Infect Dis (2006) 6(5):303–12. doi: 10.1016/S1473-3099(06)70466-2

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Randolph AG, Xu R, Novak T, Newhams MM, Bubeck Wardenburg J, Weiss SL, et al. Vancomycin Monotherapy May Be Insufficient to Treat Methicillin-Resistant Staphylococcus Aureus Coinfection in Children With Influenza-Related Critical Illness. Clin Infect Dis (2019) 68(3):365–72. doi: 10.1093/cid/ciy495

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Hughes S, Troise O, Donaldson H, Mughal N, Moore LSP. Bacterial and Fungal Coinfection Among Hospitalized Patients With COVID-19: A Retrospective Cohort Study in a UK Secondary-Care Setting. Clin Microbiol Infect (2020) 26(10):1395–9. doi: 10.1016/j.cmi.2020.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Li J, Guo M, Tian X, Wang X, Yang X, Wu P, et al. Virus-Host Interactome and Proteomic Survey Reveal Potential Virulence Factors Influencing SARS-CoV-2 Pathogenesis. Med (2020) 2(1):99–112. doi: 10.1016/j.medj.2020.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Overmyer KA, Shishkova E, Miller IJ, Balnis J, Bernstein MN, Peters-Clarke TM, et al. Large-Scale Multi-Omic Analysis of COVID-19 Severity. Cell Syst (2021) 12(1):23–40.e7. doi: 10.1016/j.cels.2020.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics (2013) 29(1):15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Anders S, Pyl PT, Huber W. HTSeq–a Python Framework to Work With High-Throughput Sequencing Data. Bioinformatics (2015) 31(2):166–9. doi: 10.1093/bioinformatics/btu638

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Love MI, Huber W, Anders S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol (2014) 15(12):550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wilk AJ, Rustagi A, Zhao NQ, Roque J, Martinez-Colon GJ, McKechnie JL, et al. A Single-Cell Atlas of the Peripheral Immune Response in Patients With Severe COVID-19. Nat Med (2020) 26(7):1070–6. doi: 10.1038/s41591-020-0944-y

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lee JS, Park S, Jeong HW, Ahn JY, Choi SJ, Lee H, et al. Immunophenotyping of COVID-19 and Influenza Highlights the Role of Type I Interferons in Development of Severe COVID-19. Sci Immunol (2020) 5(49):eabd1554. doi: 10.1126/sciimmunol.abd1554

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating Single-Cell Transcriptomic Data Across Different Conditions, Technologies, and Species. Nat Biotechnol (2018) 36(5):411–20. doi: 10.1038/nbt.4096

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. Fast, Sensitive and Accurate Integration of Single-Cell Data With Harmony. Nat Methods (2019) 16(12):1289–96. doi: 10.1038/s41592-019-0619-0

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Townsend L, Hughes G, Kerr C, Kelly M, O’Connor R, Sweeney E, et al. Bacterial Pneumonia Coinfection and Antimicrobial Therapy Duration in SARS-CoV-2 (COVID-19) Infection. JAC Antimicrob Resist (2020) 2(3):dlaa071. doi: 10.1093/jacamr/dlaa071

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Sharifipour E, Shams S, Esmkhani M, Khodadadi J, Fotouhi-Ardakani R, Koohpaei A, et al. Evaluation of Bacterial Co-Infections of the Respiratory Tract in COVID-19 Patients Admitted to ICU. BMC Infect Dis (2020) 20(1):646. doi: 10.1186/s12879-020-05374-z

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Paget C, Trottein F. Mechanisms of Bacterial Superinfection Post-Influenza: A Role for Unconventional T Cells. Front Immunol (2019) 10:336. doi: 10.3389/fimmu.2019.00336

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Mirzaei R, Goodarzi P, Asadi M, Soltani A, Aljanabi HAA, Jeda AS, et al. Bacterial Co-Infections With SARS-CoV-2. Proc Natl Acad Sci U S A (2020) 72(10):2097–111. doi: 10.1002/iub.2356

CrossRef Full Text | Google Scholar

35. Cadena AM, Hopkins FF, Maiello P, Carey AF, Wong EA, Martin CJ, et al. Concurrent Infection With Mycobacterium Tuberculosis Confers Robust Protection Against Secondary Infection in Macaques. PLoS Pathog (2018) 14(10):e1007305. doi: 10.1371/journal.ppat.1007305

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Martin-Loeches I, van Someren Greve F, Schultz MJ. Bacterial Pneumonia as an Influenza Complication. Curr Opin Infect Dis (2017) 30(2):201–7. doi: 10.1097/QCO.0000000000000347

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Hanada S, Pirzadeh M, Carver KY, Deng JC. Respiratory Viral Infection-Induced Microbiome Alterations and Secondary Bacterial Pneumonia. Front Immunol (2018) 9:2640. doi: 10.3389/fimmu.2018.02640

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wang F, Hou H, Luo Y, Tang G, Wu S, Huang M, et al. The Laboratory Tests and Host Immunity of COVID-19 Patients With Different Severity of Illness. JCI Insight (2020) 5(10) :e137799. doi: 10.1172/jci.insight.137799

CrossRef Full Text | Google Scholar

39. Lucas C, Wong P, Klein J, Castro TBR, Silva J, Sundaram M, et al. Longitudinal Analyses Reveal Immunological Misfiring in Severe COVID-19. Nature (2020) 584(7821):463–9. doi: 10.1038/s41586-020-2588-y

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Gao M, Liu Y, Guo M, Wang Q, Wang Y, Fan J, et al. Regulatory CD4(+) and CD8(+) T Cells are Negatively Correlated With CD4(+)/CD8(+) T Cell Ratios in Patients Acutely Infected With SARS-CoV-2. J Leuk Biol (2021) 109(1):91–7. doi: 10.1002/JLB.5COVA0720-421RR

CrossRef Full Text | Google Scholar

41. Manohar P, Loh B, Athira S, Nachimuthu R, Hua X, Welburn SC, et al. Secondary Bacterial Infections During Pulmonary Viral Disease: Phage Therapeutics as Alternatives to Antibiotics? Front Microbiol (2020) 11:1434. doi: 10.3389/fmicb.2020.01434

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Jamieson AM, Pasman L, Yu S, Gamradt P, Homer RJ, Decker T, et al. Role of Tissue Protection in Lethal Respiratory Viral-Bacterial Coinfection. Science (2013) 340(6137):1230–4. doi: 10.1126/science.1233632

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kubes P, Jenne C. Immune Responses in the Liver. Annu Rev Immunol (2018) 36:247–77. doi: 10.1146/annurev-immunol-051116-052415

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wu Z, Han M, Chen T, Yan W, Ning Q. Acute Liver Failure: Mechanisms of Immune-Mediated Liver Injury. Liver Int (2010) 30(6):782–94. doi: 10.1111/j.1478-3231.2010.02262.x

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Xu L, Liu J, Lu M, Yang D, Zheng X. Liver Injury During Highly Pathogenic Human Coronavirus Infections. Liver Int (2020) 40(5):998–1004. doi: 10.1111/liv.14435

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Garcia-Vidal C, Sanjuan G, Moreno-Garcia E, Puerta-Alcalde P, Garcia-Pouton N, Chumbita M, et al. Incidence of Co-Infections and Superinfections in Hospitalized Patients With COVID-19: A Retrospective Cohort Study. Clin Microbiol Infect (2021) 27(1):83–8. doi: 10.1016/j.cmi.2020.07.041

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Li N, Ren A, Wang X, Fan X, Zhao Y, Gao GF, et al. Influenza Viral Neuraminidase Primes Bacterial Coinfection Through TGF-Beta-Mediated Expression of Host Cell Receptors. Proc Natl Acad Sci U S A (2015) 112(1):238–43. doi: 10.1073/pnas.1414422112

PubMed Abstract | CrossRef Full Text | Google Scholar