Skip to main content


Front. Nutr., 13 May 2021
Sec. Nutrition and Microbes
Volume 8 - 2021 |

Gut Microbiota May Not Be Fully Restored in Recovered COVID-19 Patients After 3-Month Recovery

Yu Tian1 Kai-yi Sun1 Tian-qing Meng1 Zhen Ye1 Shi-meng Guo1 Zhi-ming Li1 Cheng-liang Xiong1 Ying Yin2* Hong-gang Li1* Li-quan Zhou1*
  • 1Institute of Reproductive Health, Center for Reproductive Medicine, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
  • 2School of Basic Medicine, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Coronavirus disease 2019 (COVID-19) has infected over 124 million people worldwide. In addition to the development of therapeutics and vaccines, the evaluation of the sequelae in recovered patients is also important. Recent studies have indicated that COVID-19 has the ability to infect intestinal tissues and to trigger alterations of the gut microbiota. However, whether these changes in gut microbiota persist into the recovery stage remains largely unknown. Here, we recruited seven healthy Chinese men and seven recovered COVID-19 male patients with an average of 3-months after discharge and analyzed their fecal samples by 16S rRNA sequencing analysis to identify the differences in gut microbiota. Our results suggested that the gut microbiota differed in male recovered patients compared with healthy controls, in which a significant difference in Chao index, Simpson index, and β-diversity was observed. And the relative abundance of several bacterial species differed clearly between two groups, characterized by enrichment of opportunistic pathogens and insufficiency of some anti-inflammatory bacteria in producing short chain fatty acids. The above findings provide preliminary clues supporting that the imbalanced gut microbiota may not be fully restored in recovered patients, highlighting the importance of continuous monitoring of gut health in people who have recovered from COVID-19.


The pandemic of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in more than 124 million confirmed cases according to Johns Hopkins University Data Archive (1), and the impact of COVID-19 on global public health could last for years to come. Therefore, it is of utmost importance to identify any potential sequelae and long-term effects given the large population of recovered patients worldwide. Physiologically, the infecting of cells with SARS-CoV-2 is mediated by the receptor angiotensin converting enzyme 2 (ACE2) and cofactor transmembrane serine protease 2 (TMPRSS2) (2). In addition to lungs, intestinal tract tissues also express a high level of ACE2 and TMPRSS2 (3). Clinical evidence indicated that gastrointestinal symptoms are common in COVID-19 patients (4) and SARS-CoV-2 RNA was detected in the fecal samples and gastrointestinal tissues (3, 5). The direct infection of the human intestinal tract and the active replication of SARS-CoV-2 were also established in the model of intestinal organoids (6). These results together support the impairment of intestines caused by SARS-CoV-2.

Notably, recent studies have demonstrated that COVID-19 inpatients have significant alterations in gut microbiota, consisting of an increase of opportunistic pathogens and loss of commensal bacteria (79). Several bacteria have been identified as biomarkers to distinguish healthy people and COVID-19 patients, and an abundance of specific microbes was found to be associated with disease severity, indicating the role of gut microbiota in the progression of COVID-19 (8, 10). The perspective that gut microbiota functioned in regulating host immunity and balancing inflammation has been accepted widely. This immunologic modulatory function is mainly executed by metabolites of microbes, such as short-chain fatty acids (SCFAs) (11). Changes in microbiota may cause abnormal intestinal immunity, which may lead to susceptibility to some diseases and involvement in their pathogenesis (12). The alternation of gut microbiota was also confirmed in other respiratory virus infections, which was recognized as the possible cause of intestinal inflammation and immune injury in influenza patients (13). Because of the hundreds of millions of recovered patients, the potential long-standing effect of infection on gut microbiota has attracted worldwide attention. However, current studies have centered mostly on the inpatients with COVID-19, whether these adverse effects continue after recovery and the characteristic of gut microbiota in the recovered patients are still unclear.

A better understanding of gut microbiota in recovered COVID-19 patients will be informative regarding the degree of recovery and further intervention and nursing they may need. In this study, the 16S rRNA sequencing analysis of the gut microbiota profile was performed in seven male recovered COVID-19 patients and seven healthy controls, with the goal of evaluating differences between the two groups. To the best of our knowledge, this is the first description of intestinal microbiota in recovered COVID-19 patients.

Methods and Materials

Participants Recruiting and Sample Collection

Participants were recruited from the Center for Reproductive Medicine, Tongji Medical College, Huazhong University of Science and Technology. Inclusion criteria for recovered COVID-19 patients included documented recovery, age from 25 to 45 years, and having detailed medical records during hospitalization and discharge certificate. The diagnosis of COVID-19 was determined by the New Coronavirus Pneumonia Prevention and Control Program (7th edition) published by the National Health Commission of China (14). Exclusion criteria included asymptomatic cases, having received antibiotics or probiotics within 2-months, gastrointestinal diseases, and severe basic diseases. Healthy controls were individuals who were enrolled at regular physical checkups in the same hospital. Ethical approval for the study was obtained by the Ethics Committee for Clinical Research of Reproductive Medicine Center, Tongji Medical College, Huazhong University of Science and Technology. Seven male cases and seven male healthy controls were recruited in this study. All participants came from Hubei province and shared a similar moderate-fat dietary habit (mainly wheat flour, rice, pork, eggs, vegetables), in which fat provided about 30% energy. Meanwhile, participants were instructed to avoid any food and drugs affecting gastrointestinal function before sampling. The collection and storage of fecal samples adhered to the relevant standards. All fecal samples were collected immediately in sterile plastic tubes and stored at −80°C until analysis. Fecal SCFAs levels were detected using the Agilent 7,890B-7,000D GC-MS/MS platform with detailed methods described in the Supplementary Methods. Blood samples were collected from participants for analyzing serum antibody and cytokines with detailed methods described in the Supplementary Methods.

DNA Extraction and 16S rRNA Sequencing

Total bacterial genomic DNA was extracted from fecal samples using the E.Z.N.A.® soil DNA Kit (Omega Bio-tek, Norcross, GA, U.S.) according to the manufacturer's protocols. The final DNA concentrations were determined by NanoDrop 2,000 UV-vis spectrophotometer (Thermo Scientific, Wilmington, USA), Thereafter, the quality of genomic DNA was checked through 1.0% agarose gel electrophoresis. The V3–V4 region of the 16S rRNA gene was amplified using primers 338F (5'-ACT CCT ACG GGA GGC AGC AG-3') and 806R (5'-GGA CTA CHV GGG TWT CTA AT-3') on the GeneAmp 9,700 PCR System (Applied Biosystems, Foster City, CA, US). The PCR reactions were conducted according to the following program: 3 min of denaturation at 95°C, 27 cycles of 30 s at 95°C, 30 s for annealing at 55°C, and 45 s for elongation at 72°C, and a final extension at 72°C for 10 min. PCR reactions were performed in triplicate 20 μl mixture containing 4 μl of 5 × FastPfu Buffer, 2 μl of 2.5 mM dNTPs, 0.8 μl of each primer (5 μM), 0.4 μl of FastPfu Polymerase and 10 ng of template DNA. Amplification products were extracted from 2% agarose gels and purified by the DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, U.S.) and quantified by QuantiFluor ™ -ST (Promega, USA). Subsequently, the PCR amplification products were pooled in equimolar and paired-end sequenced (2 × 300) on an Illumina MiSeq platform (Illumina, San Diego, USA) at the Shanghai Majorbio Bio-pharm Technology Co., Ltd (Shanghai, China).

Gut Microbiota Data Processing and Statistical Analysis

Raw FASTQ files were demultiplexed, quality-filtered by Trimmomatic, and merged by FLASH (15, 16) with the following criteria: (i) the reads were truncated at any site receiving an average quality score <20 over a 50 bp sliding window; (ii) primers were exactly matched allowing 2 nucleotide mismatching, and reads containing ambiguous bases were removed; (iii) sequences whose overlap longer than 10 bp were merged according to their overlap sequence. And subsequent high-quality reads were clustered into operational taxonomic units (OTUs) at 97% sequence similarity via USEARCH (Version 8.1). The taxonomy of each 16S rRNA gene sequence was classified using the Ribosomal Database Project (RDP) Classifier tool (Version 2.2) against the SILVA 119 16S rRNA database with a confidence threshold of 70%. OTUs were used for alpha diversity (Shannon, Simpson), richness (Chao) with a 97% threshold. Non-metric multidimensional analysis (NMDS) and principal coordinates analysis (PCoA) were performed based on the Bray-Curtis distance matrix calculated using OTU information from each sample. The differentially enriched microbes between two groups from phylum to genus were analyzed via the Wilcoxon rank-sum test. In addition, QIIME2- DADA2 flows were also run for reexamination, and the detailed process as well as results are available in Supplementary Materials. Statistical analyses were conducted in SPSS 26.0 (Chicago, Illinois, USA). Continuous variables were expressed as the mean ± SD and analyzed with the Student t-test. Categorical variables were expressed as percentages and compared by the chi-square test. P < 0.05 were accepted as indicating a significant difference.


Characteristics of Recruited Recovered COVID-19 Participants

We recruited seven male recovered COVID-19 patients, and their clinical characteristics are summarized in Table 1. The mean length of hospital stay was 19.4 days, and the time between discharge and sampling from 78 to 106 days, with an average time of 90 days. Two patients had complications during the course of the disease; of note, gastrointestinal symptoms were reported in one of the patients. According to the New Coronavirus Pneumonia Prevention and Control Program (7th edition), the degree of disease severity was classified as mild, moderate, or severe. The CoV-IgG of all patients was positive and CoV-IgM was found in only one of seven patients, which means that most patients have entered the recovery phase entirely. Moreover, the pro-inflammatory cytokines IL-6 and TNF-α showed no significant differences between recovered patients and healthy people (Supplementary Figure 1). The comparison of the baseline characteristics of COVID-19 group and healthy controls is shown in Supplementary Table 1. There was no significant difference among the groups.


Table 1. Clinical characteristics of recruited COVID-19 patients.

The Gut Microbiota Compositions of Recovered COVID-19 Patients and Controls

After merging and filtering, 756,359 high-quality sequences were acquired from 14 fecal samples by 16S rDNA gene sequencing with an average sequence number of 54,025 per sample. Subsequently, 397 OTUs were clustered at a 97% similarity level. The rarefaction curves became flat, and the Shannon index also approached a clear plateau, which together suggested that a near-complete diversity has been captured from each sample (Supplementary Figures 2A,B).

The compositions of gut microbiota of study subjects in the level of phylum and genus were shown in Figure 1A. Eight phyla were identified through bacteria microbiota analysis, and both the COVID-19 patients and healthy controls were dominated by four phyla, Firmicutes, Bacteroidota, Proteobacteria, Actinobacteriota, and Fusobacteriota. The Firmicutes were the predominant phylum, contributing 72.0 and 69.8% of the microbiota in healthy volunteers and COVID-19 patients, respectively. Actinobacteriota was the fourth most dominant phyla in healthy volunteers (9.7%), while a clear decrease in COVID-19 patients (1.3%). The second, third, and fifth most dominant phyla in healthy controls were Bacteroidota (13.4%), Proteobacteria (4.0%), and Fusobacteriota (0.8%); however, all of these phyla were elevated in the COVID-19 groups (Bacteroidota 16.1%, Proteobacteria 7.8%, Fusobacteriota 4.1%) (Figure 1A). The dominant genera of gut microbiota in both groups were Blautia, Bacteroides, Agathobacter, Faecalibacterium, and Escherichia-Shigella, with plenty of other sporadic genera (Figure 1B). The composition of the intestinal microbiota of each participant is shown in Supplementary Figure 3. Taken together, the gut microbiota of the two groups remains consistent in the overall microbiota composition but shows differences in proportions of some bacteria.


Figure 1. The relative abundance of gut microbiota at the Phylum (A) and Genus (B) levels in the healthy control group and the recovered COVID-19 patient group. The category “Other” covers all other phyla or genera with a low taxa abundance. RPs, recovered patients; HCs, healthy controls.

Comparison of Microbiota Between Recovered Patients and Healthy Controls

To investigate whether an infection has a long-term effect on the gut microbiota, we evaluated the richness and diversity of recovered patients and healthy controls using the α-diversity including the Shannon index, Chao index, and Simpson index (Figures 2AC). Recovered COVID-19 patients showed a decreasing trend in microbial diversity and richness, there was a statistically significant difference among the groups in light of richness estimator, Chao index (p = 0.0298), and diversity estimator Simpson index (p = 0.015). But there was no significant difference in another diversity estimator, the Shannon index (p = 0.055). Other α-diversity estimators are summarized in Supplementary Table 2. Meanwhile, the number of OTUs in the recovered patient group and health group were 291 and 346, and 240 OTUs overlapped between the two groups (Figure 2D). β-diversity was calculated through NMDS and PCoA in the level of OTUs. NMDS indicated that a clear separation between recovered patients and healthy controls (R = 0.02449, p = 0.0250) (Figure 2E). PCoA also produced a similar result (Supplementary Figure 4). Taken together, our results indicated a role of SARS-CoV-2 infection in affecting the fecal microbiota structure.


Figure 2. Differences in the diversity and richness of fecal microbiota between recovered COVID-19 patients and healthy controls. The Shannon (A), Chao (B), and Simpson index (C) were calculated at a 3% distance. (D) A Venn diagram indicates the overlapping and unique OTUs among the groups. (E) Non-metric multidimensional analysis (NMDS) was conducted at the OTU level to show the different distribution of recovered patients and controls.

Differences in Specific Microbes of Recovered COVID-19 Patients and Controls

Many microbes differ substantially between the two groups, characterized by the decrease in beneficial bacteria and increase in opportunistic pathogens (Supplementary Table 3). There was no significant difference in the phylum level. At the class level, recovered COVID-19 patients had a significant decrease in Coriobacteriia (p = 0.030), with an increase in Acidimicrobiia (0.031) compared with controls, partly caused by the abundant change of order Coriobacteriales and Microtrichales (Figure 3A). Additionally, at the order level, the relative abundance of Eubacteriales (0.011) and Micrococcales (0.021) increased in the recovered patients, while Oscillospirales (0.030) decreased compared with the controls (Figure 3B). At the family level, the relative abundance of Ruminococcaceae (p = 0.010) and the dramatic reduction in Coriobacteriaceae (p = 0.037) in recovered patients, with an elevation in the abundance of Eubacteriaceae (p = 0.011), Micrococcaceae (p = 0.015), Microtrichaceae (0.031) (Figure 3C). At the Genus level, recovered COVID-19 patients showed a clear decrease in Faecalibacterium (0.021), Eubacterium hallii group (0.004), Collinsella (0.037), Erysipelotrichaceae UCG-003 (0.001), NK4A214 group (0.031), and a concomitant increase in Flavonifractor (0.041), Eubacterium (0.011), Rothia (0.030), Candidatus Microthrix (0.030), especially the Erysipelatoclostridium (0.001) (Figure 3D). This is partly similar to previous results on COVID-19 inpatients (7). On the one hand, a remarkable depletion was observed in commensals, particularly in the anti-inflammatory bacteria, such as Faecalibacterium and the Eubacterium_hallii_group. On the other hand, opportunistic pathogens such as Rothia and Erysipelatoclostridium showed a higher relative abundance than controls.


Figure 3. The taxonomic differences of recovered COVID-19 patients and controls in the level of Class (A), Order (B), Family (C), and Genus (D).

To identify potential bacterial taxon biomarkers between recovered patient and healthy controls, linear discriminant analysis (LDA) effect size (LEfSe) was used to find species that differed significantly in abundance between the groups. A cladogram was used to show microbiota structures at the phylum and genus levels by an LDA score >3.5 (Figure 4). The gut microbiome of the COVID-19 group was mainly dominated by Acidimicrobiia, Microtrichales, Candidatus Microthrix, Microtrichaceae, Rothia, Micrococcales, Erysipelatoclostridium, and Micrococcaceae. These species may act as biomarkers to distinguish and analyze the outcomes of recovered patients on gut health.


Figure 4. Identifying bacterial differences by the LEfSe. (A) Taxonomic cladogram shows the taxa that were considered statistically significant between two groups from the Phylum to the Genus level. (B) the LDA score of differentially enriched taxa in healthy controls and COVID-19 group (LDA threshold value >3.5).

Levels of Fecal SCFAs in Recovered Patients and Healthy Controls

Family Ruminococcaceae, genus Faecalibacterium and Eubacterium_hallii_group showed a significant decrease in recovered patients and are SCFAs-producing bacteria (1719). Decreased production of SCFAs has been found to be a consequence in mouse models of influenza A virus infection, which was capable of affecting immune responses in the lungs and modulating disease outcomes (20). Therefore, we examined the fecal concentrations of certain SCFAs, including acetic acid (AA), propionic acid (PA), isobutyric acid (IBA), butyric acid (BA), isovaleric acid (IVA), valeric acid (VA), and hexanoic acid (HA) among the two groups (Figures 5AG). Intriguingly, there was no significant difference in the level of any SCFA. This result indicated that the SCFA levels of recovered patients were normal, though the upstream bacteria varied. We further analyzed the relationship between SCFA levels and bacterial genera using correlation heatmap based on the spearman correlation coefficient (Figure 5H). Intriguingly, there was no significant association between altered microbes and SCFAs levels. However, other SCFA-producing bacteria, such as Alistipes, Dialister, Butyricimonas, and Phascolarctobacterium, demonstrate remarkable correlation with at least one SCFA, the relative abunance of which was not decreased according to analysis mentioned above.


Figure 5. The levels of SCFAs in the fecal samples of recovered patients and healthy controls. Examination of fecal samples of the two groups on levels of (A) Acetic acid (AA), (B) Propionic acid (PA), (C) Isobutyric acid (IBA), (D) Butyric acid (BA), (E) Isovaleric acid (IVA), (F) Valeric acid (VA), and (G) Hexanoic acid (HA). (H) Correlation between relative abundance of bacteria genera and levels of seven SCFAs. *p < 0.05; **p < 0.01.

Taken together, the drop in abundance of several anti-inflammatory bacteria that can secrete SCFAs does not impact SCFA levels, probably because of supplementation from other SCFA-producing bacteria, and other potential impacts of reduced anti-inflammatory bacteria should be paid attention to.


Our study may be the first to describe the gut microbiota of people recovered from COVID-19. Compared with healthy controls, there was a clear difference in the relative abundance of gut microbiota between recovered male patients and controls. A previous report revealed that several opportunistic pathogens at the phylum level showed a clear elevation in COVID-19 inpatients, including Streptococcus, Rothia, Veillonella, Erysipelatoclostridium, and Actinomyces (7). In our study, the genera Rothia and Erysipelatoclostridium in recovered patients was also significantly more abundant than in the control group, indicating a long-term effect on gut microbiota. Rothia belongs to the oral-related microbes, which were correlated with Th17-induced lung inflammation and pneumonia of immunocompromised patients (21, 22). The elevation of Rothia in gut may suggest a possible shift of microbes from the mouth and respiratory tract to the intestines. However, increases in Streptococcus, Actinomyces, and Veillonella were observed but did not approach statistical significance in our study. These differences may suggest potential dynamic changes with the progression of COVID-19 and reflect the distinction of different stages in this disease.

We also analyzed our data by QIIME2 flow. Although there were no obvious differences in α-diversity estimators, Venn diagrams, and NMDS showed more significant separation between the two groups (Supplementary Figure 5). Additionally, the differences in relative abundance of specific microbes obtained by QIIME2 were similar to Figure 3, especially in opportunistic pathogens and SCFAs-producing bacteria (Supplementary Figures 6A–D). Besides, Rothia, Micrococcales, Erysipelatoclostridium, and Micrococcaceae were also identified using LEfSe in QIIME2 analysis (Supplementary Figures 6E,F). Taken together, results obtained from the two flows are generally consistent.

Underlying mechanisms of bacteria changes should be considered. Direct infection of virus and subsequent pathological changes including inflammation and hypoxia could be the most obvious causes of the abnormal gut microbiota (23). Besides, the use of antiviral medications and antibiotics during hospitalization also has a severe impact on the delicate balance in bacteria community. Intriguingly, ACE2, the receptor of SARS-CoV-2 also plays a role in maintaining the intestinal microbiome homeostasis, and the infection of SARS-CoV-2 may down-regulate the availability of ACE2 and induce dysfunction of the microbiota (24). The duration of these adverse effects is a key factor in evaluating the intestinal health of recovered patients. Zuo et al. demonstrated that altered gut microbiota continued after clearance of the virus when patients were discharged (8). And our study further suggests that changes in gut microbiota were not fully restored after a 3-month recovery.

Another vital question is whether these alterations in gut microbiota have a substantial impact on human health. Gut microbiota is a vital determinant of normal gut function and immunity, and potential instability may contribute to multiple diseases. SCFAs are capable of regulating the inflammatory responses by binding the G-protein-coupled receptor 43 and suppressing the activation of nuclear factor kappa-B (25, 26), which has been established as a mechanistic link between the anti-inflammation effect and the gut microbiota. Abnormality of relative abundance of some SCFAs-producing bacteria were reported in COVID-19 inpatients (7, 8). We also observed that recovered patients had a remarkable decrease in SCFA-producing bacteria, including the family Ruminococcaceae and genus Faecalibacterium and the Eubacterium hallii group, but the levels of SCFAs in recovered patients were normal. The supplementation from other SCFA-producing bacteria could be the vital point. Additionally, Sokol et al. (27) analyzed the dynamic changes of the gut microbiota in macaques over the entire course of SARS-CoV-2 infection. SCFA production declined during infection and returned to a normal level after recovery. Similar shifts may occur in COVID-19 patients. Notably, Faecalibacterium also functions in preventing inflammatory bowel disease by secreting other anti-inflammatory productions such as salicylic acid and Microbial Anti-Inflammatory Molecule (28). And the enrichment of opportunistic pathogens indicates an increased risk of chronic inflammation in the gut (18, 29). Thus, persistent monitoring of gut health in recovered patients is strongly recommended.

We acknowledge some limitations in this study. First, this is a very small sample size, and objects center on male patients, which was due to the difficulty in recruiting patients recovered from COVID-19 for just 3-months. The results provided preliminary evidence, which should be interpreted with great caution. Multicenter studies involving more cases are invited to verify our results. Second, this is a cross-sectional study in which only recovered patients were evaluated. The microbiota configuration in other stages, particularly before SARS-CoV-2 infection and during hospitalization, is helpful to reveal the thorough effect of COVID-19 on gut microbiota. Third, whether imbalanced gut microbiota causes functional impairments in recovered patients is unknown due to a lack of clinical studies. Therefore, it must be emphasized that the adverse effects that COVID-19 has induced via gut microbiota should not be overstated for now (30).

To conclude, we reported the differences of gut microbiota between recovered COVID-19 patients and controls, which propel our knowledge regarding the persistence of implications on gut microbiota mediated by COVID-19. Our findings partly support the previous studies on inpatients with COVID-19 and also indicated a possible change in gut microbiota of COVID-19 patients in different stages. The remarkable decrease of anti-inflammatory bacteria may suggest a potential risk of chronic intestinal inflammation disorders for recovered patients; thus, interventions to restore the microbiota ecology should be considered. Longstanding retrospective studies are warranted to reveal the thorough impact of COVID-19 on human health via gut microbiota and to help develop appropriate nutritional interventions for recovered patients.

Data Availability Statement

The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at:, PRJNA683685.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee for Clinical Research of Reproductive Medicine Center, Tongji Medical College, Huazhong University of Science and Technology. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

L-qZ, HL, and C-lX designed the project. YY, T-qM, ZY, S-mG, Z-mL, K-yS, and YT performed the experiments. YT and KS analyzed the data and wrote the manuscript. LZ, H-gL, and YY revised the manuscript. All authors contributed to the article and approved the submitted version.


This work was supported by the National Key R&D Program of China [2018YFC1004502 and 2018YFC1004001], HUST COVID-19 Rapid Response Call (2020kfyXGYJ057), and the National Natural Science Foundation of China [NSFC 31771661] and program for HUST Academic Frontier Youth Team.

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.


The authors would like to express their thanks to Doctor Hui Zhou, Xing-jie Zhang, and Zhen-Yu Zhong for their assistance on this study. Thanks to Shanghai Majorbio Bio-pharm Technology Co., Ltd for providing assistance on bioinformatics analysis and Wuhan Metware Biotechnology Co, Ltd for assistance on metabolic profiling.

Supplementary Material

The Supplementary Material for this article can be found online at:


1. Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. (2020) 20:533–4. doi: 10.1016/s1473-309930120-1

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell. (2020) 181:271–80.e8.

PubMed Abstract | Google Scholar

3. Xiao F, Tang M, Zheng X, Liu Y, Li X, Shan H. Evidence for gastrointestinal infection of SARS-CoV-2. Gastroenterology. (2020) 158:1831–3 e3. doi: 10.1053/j.gastro.2020.02.055

CrossRef Full Text | Google Scholar

4. Lin L, Jiang X, Zhang Z, Huang S, Zhang Z, Fang Z, et al. Gastrointestinal symptoms of 95 cases with SARS-CoV-2 infection. Gut. (2020) 69:997–1001. doi: 10.1136/gutjnl-2020-321013

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Xie C, Jiang L, Huang G, Pu H, Gong B, Lin H, et al. Comparison of different samples for 2019 novel coronavirus detection by nucleic acid amplification tests. Int J Infect Dis. (2020) 93:264–7. doi: 10.1016/j.ijid.2020.02.050

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Zhou J, Li C, Liu X, Chiu MC, Zhao X, Wang D, et al. Infection of bat and human intestinal organoids by SARS-CoV-2. Nat Med. (2020) 26:1077–83. doi: 10.1038/s41591-020-0912-6

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gu S, Chen Y, Wu Z, Chen Y, Gao H, Lv L, et al. Alterations of the gut microbiota in patients with COVID-19 or H1N1 influenza. Clin Infect Dis. (2020) 71:2669–78. doi: 10.1093/cid/ciaa709

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Zuo T, Zhang F, Lui GCY, Yeoh YK, Li AYL, Zhan H, et al. Alterations in gut microbiota of patients with COVID-19 during time of hospitalization. Gastroenterology. (2020) 159:944–55.e8. doi: 10.1053/j.gastro.2020.05.048

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zuo T, Liu Q, Zhang F, Lui GC, Tso EY, Yeoh YK, et al. Depicting SARS-CoV-2 faecal viral activity in association with gut microbiota composition in patients with COVID-19. Gut. (2020) 70:276–84. doi: 10.1136/gutjnl-2020-322294

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Gou W, Fu Y, Yue L, Chen G-D, Cai X, Shuai M, et al. Gut microbiota may underlie the predisposition of healthy individuals to COVID-19. medRxiv. (2020). doi: 10.1101/2020.04.22.20076091

CrossRef Full Text | Google Scholar

11. Rooks MG, Garrett WS. Gut microbiota, metabolites and host immunity. Nat Rev Immunol. (2016) 16:341–52. doi: 10.1038/nri.2016.42

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Budden KF, Gellatly SL, Wood DLA, Cooper MA, Morrison M, Hugenholtz P, et al. Emerging pathogenic links between microbiota and the gut–lung axis. Nat Rev Microbiol. (2016) 15:55–63. doi: 10.1038/nrmicro.2016.142

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang J, Li F, Wei H, Lian ZX, Sun R, Tian Z. Respiratory influenza virus infection induces intestinal immune injury via microbiota-mediated Th17 cell-dependent inflammation. J Exp Med. (2014) 211:2397–410. doi: 10.1084/jem.20140625

PubMed Abstract | CrossRef Full Text | Google Scholar

14. National Health Commission of the People's Republic of China. Guidelines for COVID-19 Diagnosis and Treatment (Trial version 7). (2020). Available online at:

Google Scholar

15. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. (2014) 30:2114–20. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. doi: 10.1093/bioinformatics/btr507

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Koh A, De Vadder F, Kovatcheva-Datchary P, Backhed F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell. (2016) 165:1332–45. doi: 10.1016/j.cell.2016.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ferreira-Halder CV, Faria AVS, Andrade SS. Action and function of Faecalibacterium prausnitzii in health and disease. Best Pract Res Clin Gastroenterol. (2017) 31:643–8. doi: 10.1016/j.bpg.2017.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Louis P, Young P, Holtrop G, Flint HJ. Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA-transferase gene. Environ Microbiol. (2010) 12:304–14. doi: 10.1111/j.1462-2920.2009.02066.x

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Sencio V, Barthelemy A, Tavares LP, Machado MG, Soulard D, Cuinat C, et al. Gut dysbiosis during influenza contributes to pulmonary pneumococcal superinfection through altered short-chain fatty acid production. Cell Rep. (2020) 30:2934–47.e6. doi: 10.1016/j.celrep.2020.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Segal LN, Clemente JC, Tsay JC, Koralov SB, Keller BC, Wu BG, et al. Enrichment of the lung microbiome with oral taxa is associated with lung inflammation of a Th17 phenotype. Nat Microbiol. (2016) 1:16031. doi: 10.1038/nmicrobiol.2016.31

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Maraki S, Papadakis IS. Rothia mucilaginosa pneumonia: a literature review. Infect Dis. (2015) 47:125–9. doi: 10.3109/00365548.2014.980843

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Trottein F, Sokol H. Potential causes and consequences of gastrointestinal disorders during a SARS-CoV-2 infection. Cell Rep. (2020) 32:107915. doi: 10.1016/j.celrep.2020.107915

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Hashimoto T, Perlot T, Rehman A, Trichereau J, Ishiguro H, Paolino M, et al. ACE2 links amino acid malnutrition to microbial ecology and intestinal inflammation. Nature. (2012) 487:477–81. doi: 10.1038/nature11228

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Maslowski KM, Vieira AT, Ng A, Kranich J, Sierro F, Yu D, et al. Regulation of inflammatory responses by gut microbiota and chemoattractant receptor GPR43. Nature. (2009) 461:1282–6. doi: 10.1038/nature08530

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Luhrs H, Gerke T, Muller JG, Melcher R, Schauber J, Boxberge F, et al. Butyrate inhibits NF-kappaB activation in lamina propria macrophages of patients with ulcerative colitis. Scand J Gastroenterol. (2002) 37:458–66. doi: 10.1080/003655202317316105

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sokol H, Contreras V, Maisonnasse P, Desmons A, Delache B, Sencio V, et al. SARS-CoV-2 infection in nonhuman primates alters the composition and functional activity of the gut microbiota. Gut Microbes. (2021) 13:1–19. doi: 10.1080/19490976.2021.1893113

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Miquel S, Leclerc M, Martin R, Chain F, Lenoir M, Raguideau S, et al. Identification of metabolic signatures linked to anti-inflammatory effects of Faecalibacterium prausnitzii. mBio. (2015) 6:e00300–15. doi: 10.1128/mBio.00300-15

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Hiippala K, Jouhten H, Ronkainen A, Hartikainen A, Kainulainen V, Jalanka J, et al. The potential of gut commensals in reinforcing intestinal barrier function and alleviating inflammation. Nutrients. (2018) 10:988. doi: 10.3390/nu10080988

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Walter J, Armet AM, Finlay BB, Shanahan F. Establishing or exaggerating causality for the gut microbiome: lessons from human microbiota-associated rodents. Cell. (2020) 180:221–32. doi: 10.1016/j.cell.2019.12.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: SARS-CoV-2, recovered COVID-19 patient, gut microbiota, 16S sequence, short chain fatty acids

Citation: Tian Y, Sun K-y, Meng T-q, Ye Z, Guo S-m, Li Z-m, Xiong C-l, Yin Y, Li H-g and Zhou L-q (2021) Gut Microbiota May Not Be Fully Restored in Recovered COVID-19 Patients After 3-Month Recovery. Front. Nutr. 8:638825. doi: 10.3389/fnut.2021.638825

Received: 07 December 2020; Accepted: 29 March 2021;
Published: 13 May 2021.

Edited by:

Michael Gänzle, University of Alberta, Canada

Reviewed by:

Maria De Angelis, University of Bari Aldo Moro, Italy
Senem Kamiloglu, Uludag University, Turkey

Copyright © 2021 Tian, Sun, Meng, Ye, Guo, Li, Xiong, Yin, Li and Zhou. 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: Li-quan Zhou,; Hong-gang Li,; Ying Yin,

These authors have contributed equally to this work