1. Introduction
Coronavirus disease 2019 (COVID-19) is caused by the infection of novel, highly contagious severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The virus quickly spread across the continents, and in March 2020, a pandemic was officially declared by the World Health Organization (WHO) [1]. In recent epidemiological reports, the burden of comorbidities, immunocompromised states, older age, and obesity have been strongly associated with the disease’s severity [2,3]. COVID-19 most commonly affects the pulmonary system, with the most frequent clinical presentation involving dyspnea, fever, cough, and sore throat [4]. Clinical manifestation is diverse and can range from the absence of symptoms to multi-organ failure and death [5]. In advanced cases, a patient’s clinical course can be complicated by pneumonia, acute respiratory distress syndrome, acute hypoxic respiratory failure, and death [6,7]. It is worth highlighting that SARS-CoV-2-induced intense cytokine and chemokine response (cytokine storm syndrome) plays an important role in the development of acute respiratory distress syndrome (ARDS) and multi-organ dysregulation [8]. Although respiratory symptoms are most predominantly reported amongst patients with COVID-19, extrapulmonary symptoms from gastrointestinal (GI), hepatobiliary, cardiovascular, renal, and neurological systems are also noted [5]. Moreover, following recovery, a subgroup of patients may develop autoinflammatory symptoms including multisystem inflammatory syndrome and Kawasaki-like disease in children [9,10].
Studies have provided direct evidence that the spike protein of SARS-CoV-2 serves as a crucial mode of entry into the host organism [11]. It binds to angiotensin-converting enzyme 2 (ACE2) and is cleaved by transmembrane cellular serine protease (TMPRSS2) [12]. This interaction increases the production and replication of infectious de novo virus molecules [13]. These receptors are extensively expressed in the lungs, heart, kidneys, and gut epithelial cells. This raises the possibility that ACE2 modulates intestinal immunity and the composition of gut microbiota by maintaining amino-acid transporters [14]. It is, therefore, crucial to consider the importance of the microbiome and the gut–lung axis in the context of SARS-CoV-2 infection. Emerging data revealed that COVID-19 infection cannot be detected merely in oral swabs but also in fecal specimens and anal/rectal swabs, suggesting crosstalk between the gut and the lungs [15,16,17]. Growing evidence showed that adequate interaction between the microbiota in the gut and the lungs plays a pivotal role in maintaining homeostasis. Accordingly, any alterations in the microbial community increase host vulnerability and disease outcomes and promote states of chronic low-grade inflammation [18]. Notably, our previous investigation has also confirmed the impact of dysbiosis on the inflammatory process accompanying the disease [19].
Considering the points raised, the aim of this study is to discuss the disturbances of gut microbiota based on a complete hypervariable sequence of the 16S rRNA gene in patients following three months of recovery from SARS-CoV-2 juxtaposed with healthy controls. Additionally, our goal is to verify obtained results by performing a meta-analysis with similar publicly available data.
2. Materials and Methods
2.1. Study Subject and Design
This study was conducted according to the Declaration of Helsinki and was approved by the local Ethics Committee of Poznan University of Medical Sciences (resolution no. 95/22, approved on 17 February 2022). In this prospective investigation, a total of 22 individuals—8 subjects with recovered COVID-19 infection and 14 healthy controls—were recruited. The diagnosis of SARS-CoV-2 was confirmed by a positive SARS-CoV-2 real-time reverse transcription-polymerase chain reaction (RT-PCR) throat swab. All participants provided written informed consent to participate in this research and agreed to the publication of the study results. Fecal sample collection was conducted between January 2021 and March 2021. Inclusion criteria for SARS-CoV-2 recovered subjects included age > 18 years and documented recovery history, whereas exclusion criteria encompassed asymptomatic patients. Both patients and healthy individuals have neither received probiotic or antibiotic treatment within six months nor had any severe additional diseases (including GI disorders, cancer, and diabetes).
2.2. Metaprofiling and Data Processing
2.2.1. DNA Extraction
Microbial genomic DNA was extracted from the stool samples using the QIAamp PowerFecal Pro DNA Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol. Fecal samples were collected and frozen immediately at −80 °C until DNA extraction. The isolated bacterial material was quantified and characterized using a spectrophotometer NanoDrop™ 2000 (Thermo Fisher Scientific, Inc., Waltham, MA, USA).
2.2.2. Next-Generation Sequencing (NGS)
Bacterial 16S ribosomal RNA (rRNA) gene and fungal internal transcribed spacer (ITS) libraries for metaprofiling analysis were prepared using QIAseq 16S/ITS Panels Reagent Kit Screening Panel (QIAGEN) according to the manufacturer’s protocol. The first step involved the amplification of microbial target DNA sequences in three multiplexed PCRs: (1) V1–V2, V4–V5, ITS; (2) V2–V3, V5–V7; and (3) V3–V4, V7–V9. The obtained amplification products were combined, mixed, and then cleaned using the magnetic QIAseq Beads (QIAGEN) and subjected to indexing with unique oligonucleotide sequences placed on a commercial plate QIAseq 16S/ITS 96-Index I (QIAGEN). Amplicon libraries were purified, and the following were checked using a High Sensitivity DNA Assay on an Agilent 2100 Bioanalyzer System (Agilent, Santa Clara, CA, USA) and pooled to one collective NGS library with a final concentration of 6 pM. Bacteriophage Φ-X174 (PhiX) control library was added at 3%. Sequencing was performed on an Illumina MiSeq platform using a MiSeq Reagent Kit v3 (600 cycles).
2.2.3. Bioinformatic Analysis
The 16S rRNA gene and ITS sequences were obtained and processed using a CLC Genomic Workbench 11.0 and CLC Microbial Genomics Module 1.2. (Qiagen Bioinformatics, Aarhus, Denmark). Reads were demultiplexed and trimmed. Such data have been grouped under operational taxonomic units (OTUs), relative to the SILVA v119 97% 16S rRNA database. Bioinformatics processing made it possible to obtain and export Excel files with a qualitative–quantitative list of microorganisms on different taxonomic levels: phylum, class, order, family, genus, and species. Alpha diversity was calculated using the OTUs richness, the Chao1, Shannon entropy, and Simpson index. Beta diversity was specified based on a Bray–Curtis dissimilarity matrix and with unweighted and weighted Unifrac distance.
2.3. Statistical Analysis
The normality of demographic data and microbiota percentage abundance levels at different taxonomic levels were assessed using the Shapiro–Wilk test. For nominal data, the chi-square test and t-test for independent groups were used. When the data did not follow the normal distribution pattern, the non-parametric Mann–Whitney test for median values was applied for the comparison. Furthermore, a correction concerning the multiple hypothesis testing was performed. The adjusted p-values were calculated according to the Benjamini–Hochberg procedure. To investigate the relationship between the particular microbial phyla, principal component analysis (PCA) was performed. Calculations were carried out using PQStat 1.8.2 (PQStat Software, Poznan, Poland) and Graph Pad Prism 8 (GraphPad Software, Inc., San Diego, USA) software, and all tests were considered significant at p < 0.05.
2.4. Search and Selection of Studies for Integrative Research
We conducted a systematic search of studies in PubMed, Embase, The Cochrane Library, and the Web of Science database including the keywords (‘Intestinal’ [MeSH] OR ‘Gut’ [MeSH]) AND (‘Bacteria’ [MeSH] OR ‘Microbiota’ [MeSH] OR ‘Microflora’ [MeSH] OR ‘Microbiome’ [MeSH] OR ‘Dysbiosis’ [MeSH]) AND (‘COVID-19’ [MeSH] OR ‘COVID19’ [MeSH] OR ‘SARS-CoV-19’ [MeSH] OR ‘SARS-CoV-19’ [MeSH]) AND (‘16S rRNA’ [MeSH] OR ‘NGS’ [MeSH]) with the following filters: publication date from 1 January 2020–25 April 2022, Humans, Adult. We also performed manual searches. The last search update was conducted on 25 April 2022. Including criteria were case–control studies with adult patients diagnosed with COVID-19 with nose–pharyngeal swab-positive PCR for SARS-CoV-2. Healthy individuals without a history of disease were used as controls. Assessments of microbiota composition were based on stool samples using NGS technology covering V1-V9 regions of 16S rRNA gene. Studies with publicly available raw 16S data for each case and control samples were used (mostly downloaded from online repositories, as Sequence Read Archive, SRA).
3. Results
A total of 22 subjects, 8 patients infected with SARS-CoV-2 (five females and three males) and 14 healthy individuals (eleven males and three females), were included in the study. The clinical characteristics of patients are presented in Table 1. The differences in sex, age, and BMI were not statistically significant between the analyzed groups (p > 0.05). With regard to supplements, two COVID-19 patients and two healthy controls were taking magnesium supplements, whereas two COVID-19 and one healthy subject were consuming vitamin B complex. The origin of the study group is from one region of the country. There were no subjects on a special diet such as vegan, vegetarian, gluten-free, or high- or low-protein.
In the NGS analysis, the obtained parameters indicate the high quality of the results. The flow cell cluster density was 911 ± 32. The Phred base-calling score GQX ≥ 30 of obtained sequencing data was an average of 78.14%, and the reads passing filter (PF) was 82.72%. Therefore, the mean PF reads were 43,870 ± 11,398 for patients and 33,603 ± 14,978 for controls (Figure 1A). The PF reads of all groups were efficiently classified as bacterial taxa in SILVA v119 97% 16S rRNA database.
In our study, fecal samples of COVID-19 patients demonstrated greater microbial taxonomic diversity compared to the samples collected from healthy individuals based on the deepest available analysis of all regions V1-V9 of the 16S rRNA gene. Our strategy, following the latest data from the literature [20], was aimed at obtaining model results because simplified analyses covering selected regions of this gene (e.g., V3-V4) may cause a quantitative and qualitative disturbance of the identified microbiota taxa. The total OTUs richness for the COVID-19 group (mean OTUs value 615) was statistically significantly higher (p < 0.001) compared to control subjects (mean OTUs value 420). Similarly, the Chao1 index was significantly greater (p < 0.01) for patients (Figure 1A). Moreover, Shannon entropy and Simpson’s index were calculated to present the alpha diversity of both groups (Figure 1A). Analysis of the beta diversity determined by the Bray–Curtis dissimilarity revealed that the bacterial microbiota of the COVID-19 patient group did not cluster apart from that of a control group. UniFrac distance calculation based on the gut microbiota composition and the phylogeny diversity (sequence distance) both for the absence (unweighted) and for the relative abundance (weighted) of microbiota indicated no significant clustering for both COVID-19 patients and healthy controls (Figure 1B).
We have analyzed the taxonomic composition of the gut microbiota in both the recovered COVID-19 patients and the control group at the phylum, class, order, family, genus, and species. A total of 17 phyla were detected (Table S1). It has been brought to our attention that the relative abundance of phylum Verrucomicrobia, Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria were the dominant bacteria in both groups (Figure 2).
The mean abundance of the phylum Verrucomicrobia was statistically greater (p = 0.004) and 30 times higher in SARS-CoV-2 patients compared to healthy controls (0.7378% versus 0.0241%, respectively) (Figure 3B). Nonetheless, the phyla Proteobacteria, Bacteroidetes, and Actinobacteria were found to have decreased, whereas the Firmicutes level was higher in COVID-19 individuals than in healthy controls (Figure 3B). Moreover, the Bacteroidetes/Firmicutes ratio was reduced significantly in both COVID-19 patients and healthy controls (a value below 2) (Figure 3A).
To investigate the correlations between the gut microbiota at the phylum level between COVID-19 patients and healthy controls, principal component analysis (PCA) in two dimensions (PC) was performed. Obtained results demonstrated that in the COVID-19 subjects, Bacteroidetes and Verrucomicrobia richness were strongly correlated with each other. Additionally, a similar association was observed for Cyanobacteria and Lentisphaerae. Interestingly, these relationships did not appear in the control group. Moreover, for the healthy subjects, the abundance of Cyanobacteria and SHA-109 were in a strong correlation, which was not found in COVID-19 patients (Figure 4).
Furthermore, a total of 37 classes of microorganisms were detected, and the mean abundance of the Verrucomicrobiae class in patients was 30 times higher than that of healthy individuals (0.7316% versus 0.0238%, relatively; Figure 5A; Table S2). Metaprofiling analysis detected a total of 66 orders of microorganisms. An increased abundance of Verrucomicrobiales, Neisseriales, and Cytophagales was observed over the course of COVID-19 recovery compared to controls (Figure 5B; Table S3). Moreover, a comparative analysis of the microbial community in both groups revealed the presence of a total of 112 families of microorganisms. Seven families, including the Verrucomicrobiaceae, Neisseriaceae, Leuconostocaceae, Phyllobacteriaceae, Lachnospiraceae, Clostridiaceae 1, and Uncultured bacterium-06, were initially significantly enriched in COVID-19 patients compared to healthy controls (p < 0.05) (Table S4). After the false discovery rate (FDR) correction, statistical significance remained for the family Verrucomicrobiaceae (Figure 5C).
Additionally, using NGS metaprofiling, we have identified a total of 232 genera of microorganisms. The abundance of Akkermansia, Bilophila, Adlercreutzia, and Neisseria was enriched in COVID-19 patients (p < 0.05) (Table S5; Figure 5D,E). Following SARS-CoV-2 recovery, over 50% of bacteria in patients belonged to four genera: Bacterioides, Incertae Sedis-04, Ruminococcus, and Blautia. Comparatively, when healthy individuals were included, other four genera were dominant: Bacterioides, Incertae Sedis-04, Faecalibacterium, and Ruminococcus, which constituted 43.8% of the intestinal microbiota. Furthermore, a total of 561 species of microorganisms was detected (Table S6; Figure 6A).
Up to 95.5% of microorganisms that make up the intestinal microbiome of COVID-19 patients have not been cultured in vitro and/or are unknown. Importantly, almost all initial significant differences after the Mann–Whitney test in the quantitative microbiota composition of patients and controls (except for Parabacteroides distasonis and Bilophila wadsworthia) were observed among the previously uncharacterized and uncultured in vitro microorganisms (average abundance 0.0315% versus 0.0024%, p = 0.0289 and 0.0211% versus 0.0044%, p = 0.0404, Table S6; Figure 6B). It is worth noting that in patients, of the known species, Escherichia coli dominated, which on average accounts for 3.42% of the microbiota composition, whereas in the healthy controls, it encompassed only 0.23%. However, this difference is not statistically significant, similar to the other functionally important bacterial species (Table S6; Figure 6C).
Our electronic search identified 18 potentially eligible entries (Figure 7). Among these results, only one was compliant with our methodology (NGS analysis of V1–V9 regions of the 16S rRNA gene) [21]. The rest of the research was excluded using the narrowed-down NGS methodology. Moreover, not all of them were concerned with stool samples but also rectal swabs [22], nasopharyngeal swabs [23], or saliva samples [24]. Finally, after a detailed text analysis of the included manuscript, the determinant of its rejection was the criterion for selecting the studied groups. The subject of this work was not the gut microbiota of COVID-19 patients but the impact of the psychological stress of frontline healthcare workers fighting against COVID-19 on intestinal dysbiosis [21].
4. Discussion
In the current study, we compared the gut microbiota compositions of subjects recovered from SARS-CoV-2 and healthy controls. Interestingly, our data have revealed that fecal samples of COVID-19 patients showed statistically greater microbial taxonomic diversity compared to the samples collected from healthy subjects.
The human microbiota consists of numerous microorganisms, including bacteria, viruses, fungi, and archaea, which play an imperative role in the preservation of intestinal microbiota [25]. The healthy human intestinal flora has a greater abundance of phyla Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria; however, the colon often encompasses bacterial families of Bacteroidaceae, Lachnospiraceae, Prevotellaceae, Riknellaceae, and Ruminococcaceae [26]. Studies have shown that host diet, genetics, mode of birth, age, environmental exposure, and medication consumption contribute to the function and composition of the microbiota [27]. Therefore, its qualitative and quantitative arrangement is dynamic and changes under the influence of the aforementioned factors. Due to its essential role, it is plausible to assume that the microbial community plays a fundamental role in the maintenance of regulatory pathways and the coordination of innate and adaptive responses [28]. It is worth highlighting that about 70–80% of the body’s immune cells are located in the (GI) tract [29]. Notably, studies have reported the importance of crosstalk between the gut–lung axis and the intestinal microbiota–mucosal immune system [30].
Moreover, growing evidence suggests that lung immunity can be affected by alterations in the taxonomic composition and reduced diversity of the microbiota. In the literature, dysbiosis is defined as an imbalance in the gut microbiota composition or a change in the bacterial distribution or metabolic activities within the bowel [31]. These alterations have been proven to disrupt the intestinal mucosal barrier integrity, perpetuating systemic inflammation. It has been reported that increased gut permeability may lead to bacterial translocation from distant organs into the intestines via the lymphatic and circulatory systems [32]. This mechanism is believed to be associated with the pathogenesis of sepsis and acute respiratory distress syndrome [33]. An interesting study presented direct evidence that SARS-CoV-2 is translocated from the pulmonary system into the enterocytes, such that it influences the microbial community composition [34]. This disruption of the corresponding intestinal microbiome explains why over 60% of patients infected with SARS-CoV-2 present symptoms such as nausea, vomiting, and diarrhea [35].
Previous reports have highlighted the importance of the gut microbiota–lung axis and related respiratory infections [36]. Consistently, an increased abundance of opportunistic pathogens, including phylum Actinomyces, Erysipelatoclostridium, Streptococcus, Rothia, and Veillonella, in COVID-19 patients was reported [37]. In our study, the fecal metagenomic analysis revealed a statistically greater abundance of the genera Akkermansia, Bilophila, Adlercreutzia, and Neisseria, indicating that the dysbiosis persisted during SARS-CoV-2 infection and was not fully restored after 3 three months’ recovery. Moreover, the depletion of potentially beneficial bacterial families Ruminococcaceae and Lachnospiraceae was noted, which is in line with other studies [37]. Additionally, a previous report revealed an increased abundance of opportunistic bacteria Neisseria and Rothia in oral and intestinal microbiota during the first few days after COVID-19 symptom onset [38]. Our metaprofiling analysis has shown a statistically significant increase in Neisseriaceae, which appear vastly in the upper airways. These correlations may imply the bidirectional relationship between lungs and gut microbiota. It is worth noting that the fecal microbiota of subjects infected with SARS-CoV-2 was characterized by the reduction of beneficial bacteria and enriched by potentially harmful commensals (Bacteroides nordii and Clostridium hathewayi). Our findings are consistent with these results; however, the increased abundance of Bacteroidetes and Firmicutes did not reach statistical significance. More importantly, the severity of COVID-19 infection has been vastly related to gut dysbiosis [39]. In our study, the class Clostridia in recovered patients was significantly more abundant than in the control group. In line with that, a pilot study reported that the baseline abundance of Clostridium hathewayi, Clostridium ramosum, and Coprobacillus positively correlates with disease severity, whereas the reduced abundance of Faecalibacterium prausnitzii inversely correlates with SARS-CoV-2 infectivity [37]. Notably, increased levels of blood markers and inflammatory cytokines such as C-reactive protein, gamma-glutamyltransferase, aspartate aminotransferase, lactate dehydrogenase, IL-10, and TNF-alfa have been reported. Interestingly, in a cohort study, stool samples of COVID-19 patients revealed a lower relative abundance of Eubacterium rectale, Faecalibacterium prausnitzii, and multiple Bifidobacterium species [40]. Our findings are consistent with a previous study conducted in SARS-CoV-2 patients, in which a greater abundance of Akkermansia muciniphila and Bacteroides dorei was noted [40]. Moreover, these microorganisms positively correlated with C-X-C motif ligand 8 (CXCL8), interleukin-6 (IL-6), and IL-1β [40]. In a recent cross-sectional study, authors reported significantly greater levels also of IL-6 [41]. Moreover, subjects with COVID-19 had an increased abundance of potentially harmful pathogens including Actinomyces, Veillonella, Streptococcus, and Rothia. It is worth emphasizing that studies have shown that the prevalence of Bacteroidetes stercoris, Parabacteroides merdae, Alistipes onderdonkii, and Lachnospiraceae bacterium correlates with decreased COVID-19 infectivity [42]. In light of reports, the depletion of B. adolescentis, F. prausnitzii, E. rectale, R. (Blautia) obeum, and D. formicigenerans markedly increases microbial-induced immune dysregulation [40]. As such, NF-κB suppression by B. adolescentis has been proven to reduce the expression of pro-inflammatory cytokines. Moreover, in patients with severe COVID-19, markedly elevated levels of inflammatory markers such as IL-1, IL-6, and TNF-alfa have been detected [43]. Furthermore, clinical observations point out that subjects infected with SARS-CoV-2 have a significantly reduced abundance of short-chain-fatty-acid-producing bacteria (short-chain fatty acids—SCFAs), including Eubacterium hallii, Fusicatenibacter, and Faecalibacterium praustnitzii [41,44]. These findings are in line with our study, in which fecal samples collected from the recovered COVID-19 subjects showed a decreased abundance of the butyrate-producing bacteria of the family Lachnospiraceae and Ruminococcaceae. Emerging data revealed that lactate/sugar-utilizing bacteria including Anaerostipes species along with Eubacterium hallii can convert acetate and lactate into butyrate [44]. Similarly, our COVID-19-recovered patients had reduced fecal abundance of Anaerostipes species compared to healthy controls. Furthermore, apart from the aforementioned microorganisms, phylum Actinobacteria, including Bifidobacterium species, yields SCFA lactate and acetate during carbohydrate fermentation [45]. Moreover, our metaprofiling analysis revealed an increased abundance of phylum Verrucomicrobia and the genus Akkermansia. Apart from acetate and propionate production, this microorganism stimulates mucus degradation [46]. Some studies have also reported the importance of Akkermansia muciniphila in intestinal wound healing, induction of the immune response, and production of cytokines, IL-10, and SCFAs [47,48]. These molecules have been reported to regulate the intestinal immune system and mucosal permeability by preventing circulation and translocation of enteric endotoxins and their metabolites. Studies have shown that butyrate, one of the SCFAs in the intestinal epithelial cells, increases the formation and release of IL-18, which is crucial in maintaining homeostasis during inflammatory processes [49,50]. Mounting evidence suggests that at the molecular level, butyrate enhances the intestinal barrier by favoring the redistribution of tight-junction proteins occludin and zonula occludent-1 (ZO-1) [51]. Moreover, this key molecule increases oxygen formation and, therefore, stabilizes the hypoxia-inducible transcription factor (HIF) engaged in antimicrobial peptide production [52]. Furthermore, data have shown that consuming a fiber-rich diet increases the relative abundance of families Rikenellaceae, Porphyromonadaceae, and Lachnospiracea. Additionally, recent data have revealed that in subjects infected with SARS-CoV-2, proline and arginine metabolism was upregulated, indicating intestinal barrier dysfunction. These amino acids play a fundamental role in numerous chemical processes, such as the urea cycle, protein synthesis, and nitric oxide production [53]. Interestingly, it was also reported that one of the mechanisms of the pathogenesis of gut dysbiosis is linked to angiotensin-converting enzyme 2 expression. Recent investigations in the murine colon have shown that Bacteroides species such as Bacteroides massiliensis, Bacteroides thetaiotaomicron, Bacteroides dorei, and Bacteroides ovatus downregulate ACE2 receptor expression, which inversely correlates with COVID-19 severity [37]. Moreover, receptor downregulation markedly reduces antimicrobial peptide secretion as well as enteric tryptophan (Try) absorption [54]. A recent report demonstrated that plasma levels of this essential amino acid have been significantly reduced in ACE2-depleted mice [55]. Tryptophan is actively absorbed from the intestines and serves as the main precursor for the synthesis of tryptamine and serotonin neurotransmitters [56]. Emerging evidence showed that its deficiency plays a key role in the development of anxiety, depression, fatigue, muscle weakness, sleep, and mood disturbances. Moreover, Ninomiya et al. revealed that low levels of tryptophan underlie the pathogenesis of muscle atrophy [57]. Furthermore, this propensity may imply that Try depletion causes weakening of the diaphragm muscle and, therefore, dyspnea.
Furthermore, our study has revealed an enrichment of pathobionts in the gut microbiota of COVID-19 subjects. These include E. coli, which is known for its pro-inflammatory properties, and B. wadsworthia, which utilizes the activity of the sulfated amino acid taurine in the production of hydrogen sulfide and the rapid catalase reaction, thereby contributing to the pathogenesis of various diseases. Additionally, an increased abundance of other bacteria with pathogenic potentials, such as Parabacteroides distasonis and Enterococcus sp. CGLBL186, were noted. Moreover, our metaprofiling analysis showed an increased level of E. hallii, which contributes to intestinal propionate formation and glucose utilization [58]. On the other hand, the depletion of toxin-producing pathogenic Citrobacter freundii was reported. Similarly, a reduced abundance of B. fragilis was noted, which may be responsible for the development of endocarditis, meningitis, septic arthritis, and osteomyelitis. Interestingly, it has been brought to our attention that the beneficial commensals, i.e., Bifidobacterium ruminantium and Faecalibacterium prausnitzii, were decreased in the COVID-19 group (Figure 6C). It is plausible to assume that a reduced Firmicutes/Bacteroidetes ratio in SARS-CoV-2 subjects may indicate ongoing inflammatory processes and may be used as an inflammatory marker [59] (Figure 3).
The research we have carried out has several important strengths. Our metaprofiling analysis, including all hypervariable regions (V1–V9) of the 16S rRNA marker gene, gives a complete, deep picture of the microbial community composition in samples collected from both studied groups. In contrast to the commonly used 16S rRNA taxonomic marker narrowed down to the V3–V4 regions, our research allowed us to identify all the possible families, genera, and species of bacteria with equal affinity. Earlier studies have shown that analyses based on all variable regions of the 16S rRNA gene can deliver much more data than the selected regions of this gene [20]. Moreover, it is worth emphasizing that our studied cohorts are homogeneous and chosen according to strict criteria (Table 1). The patient’s group included only subjects with a mild course of COVID-19 disease, without antibiotic and probiotic intake in the last six months. At the same time, the collected controls of healthy individuals were checked during internal control analysis to confirm group homogeneity. For this purpose, the control group we divided into two subgroups, and the metadata of microbiota composition at all taxonomic levels we compared. No significant differences were found.
It is acknowledged that this pilot study has certain limitations. Firstly, the small sample size of enrolled individuals, particularly in the SARS-CoV-2 recovered patients. Thus, to verify our findings, we conducted integrative research with other similar studies applying the same methodology pipeline in the comparative analysis of COVID-19 patients versus healthy controls microbiota. Surprisingly, we found only one manuscript (from two months ago) that included a full 16S rRNA gene analysis but, nevertheless, in a different clinical aspect [21] which confirms our approach and results are pioneering, and the meta-analysis was not able to be performed. Therefore, studies including a larger population are required for further validation.
Secondly, stool collection was taken at only a one-time point at 3 months following the COVID-19 symptom onset. Thus, future exploratory studies should ideally collect fecal material at SARS-CoV-2 onset and during the disease course, followed by post-infection recovery. Thirdly, presented microbiome alterations not only may be caused by COVID-19 but also may be associated with co-infections with different pathogens or other diseases.
5. Conclusions
In conclusion, our findings highlight the importance of microbial community in the course of SARS-CoV-2 infection. Moreover, the gut microbiota has the propensity to drive inflammatory responses by influencing the release of pro- and anti-inflammatory cytokines, including antiviral responses at distal mucosal sites, including the lungs. Additionally, current research provides strong evidence for the gut–lung crosstalk in the pathogenesis of COVID-19 and possible detection of the viral particles in the feces following recovery. It is plausible to assume that the extent of intestinal dysbiosis and increased gut permeability may contribute to COVID-19 infection severity. Considering the data published to date on SARS-CoV-2, we may conclude that adjunctive therapies focusing on a correction of intestinal dysbiosis may be a crucial therapeutic approach, as previously shown for other disorders, such as inflammatory bowel diseases [19,48]. This long-range protection conferred by the microbiota has been attributed to metabolites secreted by specific bacterial species. Hence, metabolomics may also be the direction of further research into the influence of the microbiota on COVID-19 infection.
Considering the gut microbiota’s ability to disturb the immune response, our findings suggest the importance of the enteric microbiota in the course of SARS-CoV-2 infection. Mounting evidence reveals that increased intestinal permeability and enhanced release of bacterial toxins and metabolites contribute to the severity of COVID-19 infection. Moreover, this pilot study shows that the composition of the microbial community may not be fully restored in individuals with SARS-CoV-2 following a 3-month recovery.
D.M.-W.—conceptualization, writing—review and editing; A.W.—writing—original draft; J.Z. and M.S.-Z.—bioinformatic and statistical analysis and data interpretation; J.Z., O.Z.-B. and M.S.-Z.—molecular analysis and graphical presentation; M.G., D.M. and A.B.—patients’ recruitment and sample collection, R.S. and A.D.—supervision. All authors have read and agreed to the published version of the manuscript.
This study was conducted in accordance with the Declaration of Helsinki and was approved by the local Ethics Committee of Poznan University of Medical Sciences (resolution no. 95/22 approved on 17 February 2022).
Informed consent was obtained from all subjects involved in the study.
The datasets generated for this study were deposited in the Short Read Archive database (SRA) (
Authors declare no conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Figure 1. Alpha and beta diversity of the gut microbiota between COVID-19 patients and healthy controls. (A) Richness with operational taxonomic units (OTUs) number and within-sample diversity calculated with the Chao1, Shannon entropy, and Simpson index (median, quartile Q1–Q3, min–max). Asterisks denote statistically significant differences given by the Mann–Whitney test (** p < 0.01, *** p < 0.001) (B). Comparison between sample diversities by principal coordinates analysis (PCoA) with a Bray–Curtis dissimilarity matrix and with unweighted and weighted Unifrac distance. Red and blue represent COVID-19 patients and healthy control data, respectively.
Figure 2. Individual (A,B) and group (C) profiles of gut microbial composition at the phylum and genus level in the COVID-19 patient and healthy control groups. Charts for group profiles show mean relative abundance.
Figure 3. Bacteroidetes/Firmicutes ratio (A) and the main phylum-level abundance (B) in the gut microbiota among COVID-19 patients and healthy controls (median, quartile Q1–Q3, min–max, and mean values with SD, respectively). ** p < 0.01 calculated using the non-parametric Mann–Whitney test based on medians.
Figure 4. Projection of the microbial composition at phylum level on the two-dimensions biplot for COVID-19 patients (A) and healthy controls (B). PC1—principal component 1, PC2—principal component 2.
Figure 5. Differences in microbial composition between COVID-19 patients and healthy controls at class (A), order (B), family (C), and genus (D,E) level (mean values with SD). These data only concern initial statistically significant results calculated using the non-parametric Mann–Whitney test based on medians, and asterisks mark statistical significance obtained after multiple hypothesis testing corrections according to the Benjamini–Hochberg procedure * p < 0.05.
Figure 6. Species abundance distribution and difference in microbial composition between COVID-19 patients and healthy controls. (A) Venn diagram for identified in the study gut microbiota species. Red and blue represent COVID-19 patients and healthy controls data, respectively. The intersection shows the species (543) shared with both groups. (B) Heatmap of abundance differences in studied group at species level (mean values). These data only concern initial statistically significant results calculated using the non-parametric Mann–Whitney. (C) Observed multiple differences in the amount of bacteria abundance (↓—reduction, ↑—increase) between COVID-19 patients’ microbiota compared to healthy controls.
Figure 7. Flowchart of literature search. N—number of entries (research manuscripts). V1–V9—hypervariable regions of the 16S rRNA gene.
Demographic and clinical data of the study cohort.
Parameter | Total |
COVID-19 Patients |
Healthy Controls |
p-Value |
---|---|---|---|---|
Gender | 0.0541 | |||
Female, n (%) | 8 (36.36) | 5 (62.5) | 3 (21.43) | |
Male, n (%) | 14 (63.64) | 3 (37.5) | 11 (78.57) | |
Age, (years) mean ± SD |
36.23 ± 12.51 |
34.88 ± |
37 ± 110.91 |
0.6680 |
BMI (kg/m2), mean ± SD |
25,64 ± 3.525 |
24.86 ± 3.784 |
26,09 ± 3.431 |
0.4468 |
Smoker, n (%) | 0 (0.00) | 0 (0.00) | 0 (0.00) | NA |
Probiotics taking, n (%) | 0 (0.00) | 0 (0.00) | 0 (0.00) | NA |
Supplements taking, n (%) | 7 (31.82) | 4 (50.00) | 3 (21.43) | 0.1663 |
Antibiotics, last 6 months, n (%) | 0 (0.00) | 0 (0.00) | 0 (0.00) | NA |
COVID-19 symptoms | ||||
Cough, n (%) | 3 (13.64) | 3 (37.50) | - | NA |
Fever, n (%) | 5 (22.73) | 5 (62.50) | - | NA |
Dyspnea, n (%) | 1 (4.55) | 1 (12.50) | - | NA |
Headache, n (%) | 3 (13.64) | 3 (37.50) | - | NA |
Sinus pain | 2 (9.09) | 2 (25.00) | - | NA |
Myalgia, n (%) | 3 (13.64) | 3 (37.50) | - | NA |
Chest pain, n (%) | 2 (9.09) | 2 (25.00) | - | NA |
Fatigue, n (%) | 6 (27.27) | 6 (75.00) | - | NA |
Taste loss, n (%) | 5 (22.73) | 5 (62.50) | - | NA |
Sense of smell loss, n (%) | 4 (18.18) | 4 (50.00) | - | NA |
Symptom duration (days), mean ± SD (min.–max.) | - | 7.125 ± 3.40 (4–14) | - | NA |
SD—standard deviation; NA—not analyzed; p-values calculated by the non-parametric Mann–Whitney test (age, BMI) and Chi-square test (the remaining).
Supplementary Materials
The following are available online at
References
1. WHO Director-General’s Opening Remarks at the Media Briefing on COVID-19—11 March 2020. Available online: https://www.who.int/director-general/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-COVID-19---11-march-2020 (accessed on 16 December 2021).
2. Sharma, A.; Bhatt, N.S.; St Martin, A.; Abid, M.B.; Bloomquist, J.; Chemaly, R.F.; Dandoy, C.; Gauthier, J.; Gowda, L.; Perales, M.-A. et al. Clinical Characteristics and Outcomes of COVID-19 in Haematopoietic Stem-Cell Transplantation Recipients: An Observational Cohort Study. Lancet Haematol.; 2021; 8, pp. e185-e193. [DOI: https://dx.doi.org/10.1016/S2352-3026(20)30429-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33482113]
3. Jordan, R.E.; Adab, P.; Cheng, K.K. Covid-19: Risk Factors for Severe Disease and Death. BMJ; 2020; 368, m1198. [DOI: https://dx.doi.org/10.1136/bmj.m1198] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32217618]
4. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X. et al. Clinical Features of Patients Infected with 2019 Novel Coronavirus in Wuhan, China. Lancet; 2020; 395, pp. 497-506. [DOI: https://dx.doi.org/10.1016/S0140-6736(20)30183-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31986264]
5. Johnson, K.D.; Harris, C.; Cain, J.K.; Hummer, C.; Goyal, H.; Perisetti, A. Pulmonary and Extra-Pulmonary Clinical Manifestations of COVID-19. Front. Med.; 2020; 7, 526. [DOI: https://dx.doi.org/10.3389/fmed.2020.00526]
6. Chen, G.; Wu, D.; Guo, W.; Cao, Y.; Huang, D.; Wang, H.; Wang, T.; Zhang, X.; Chen, H.; Yu, H. et al. Clinical and Immunological Features of Severe and Moderate Coronavirus Disease 2019. J. Clin. Investig.; 2020; 130, pp. 2620-2629. [DOI: https://dx.doi.org/10.1172/JCI137244]
7. Mao, R.; Qiu, Y.; He, J.-S.; Tan, J.-Y.; Li, X.-H.; Liang, J.; Shen, J.; Zhu, L.-R.; Chen, Y.; Iacucci, M. et al. Manifestations and Prognosis of Gastrointestinal and Liver Involvement in Patients with COVID-19: A Systematic Review and Meta-Analysis. Lancet Gastroenterol. Hepatol.; 2020; 5, pp. 667-678. [DOI: https://dx.doi.org/10.1016/S2468-1253(20)30126-6]
8. Ye, Q.; Wang, B.; Mao, J. The Pathogenesis and Treatment of the ‘Cytokine Storm’ in COVID-19. J. Infect.; 2020; 80, pp. 607-613. [DOI: https://dx.doi.org/10.1016/j.jinf.2020.03.037]
9. Abid, M.A.; Nunley, L.; Abid, M.B. Could Coronavirus Disease 2019 (COVID-19) Render Natural Immunity to Re-Infections? A Spotlight on the Therapeutic Pipeline. Front. Immunol.; 2020; 11, 1294. [DOI: https://dx.doi.org/10.3389/fimmu.2020.01294]
10. Viner, R.M.; Whittaker, E. Kawasaki-like Disease: Emerging Complication during the COVID-19 Pandemic. Lancet; 2020; 395, pp. 1741-1743. [DOI: https://dx.doi.org/10.1016/S0140-6736(20)31129-6]
11. Li, W.; Moore, M.J.; Vasilieva, N.; Sui, J.; Wong, S.K.; Berne, M.A.; Somasundaran, M.; Sullivan, J.L.; Luzuriaga, K.; Greenough, T.C. et al. Angiotensin-Converting Enzyme 2 Is a Functional Receptor for the SARS Coronavirus. Nature; 2003; 426, pp. 450-454. [DOI: https://dx.doi.org/10.1038/nature02145]
12. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.-H.; Nitsche, A. et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell; 2020; 181, pp. 271-280.e8. [DOI: https://dx.doi.org/10.1016/j.cell.2020.02.052] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32142651]
13. Yu, Z.; Yang, Z.; Wang, Y.; Zhou, F.; Li, S.; Li, C.; Li, L.; Zhang, W.; Li, X. Recent Advance of ACE2 and Microbiota Dysfunction in COVID-19 Pathogenesis. Heliyon; 2021; 7, e07548. [DOI: https://dx.doi.org/10.1016/j.heliyon.2021.e07548] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34296023]
14. Vuille-dit-Bille, R.N.; Camargo, S.M.; Emmenegger, L.; Sasse, T.; Kummer, E.; Jando, J.; Hamie, Q.M.; Meier, C.F.; Hunziker, S.; Forras-Kaufmann, Z. et al. Human Intestine Luminal ACE2 and Amino Acid Transporter Expression Increased by ACE-Inhibitors. Amino Acids; 2015; 47, pp. 693-705. [DOI: https://dx.doi.org/10.1007/s00726-014-1889-6]
15. Xu, Y.; Li, X.; Zhu, B.; Liang, H.; Fang, C.; Gong, Y.; Guo, Q.; Sun, X.; Zhao, D.; Shen, J. et al. Characteristics of Pediatric SARS-CoV-2 Infection and Potential Evidence for Persistent Fecal Viral Shedding. Nat. Med.; 2020; 26, pp. 502-505. [DOI: https://dx.doi.org/10.1038/s41591-020-0817-4]
16. Zhang, W.; Du, R.-H.; Li, B.; Zheng, X.-S.; Yang, X.-L.; Hu, B.; Wang, Y.-Y.; Xiao, G.-F.; Yan, B.; Shi, Z.-L. et al. Molecular and Serological Investigation of 2019-NCoV Infected Patients: Implication of Multiple Shedding Routes. Emerg. Microbes Infect.; 2020; 9, pp. 386-389. [DOI: https://dx.doi.org/10.1080/22221751.2020.1729071]
17. Tang, A.; Tong, Z.; Wang, H.; Dai, Y.; Li, K.; Liu, J.; Wu, W.; Yuan, C.; Yu, M.; Li, P. et al. Detection of Novel Coronavirus by RT-PCR in Stool Specimen from Asymptomatic Child, China. Emerg. Infect. Dis.; 2020; 26, pp. 1337-1339. [DOI: https://dx.doi.org/10.3201/eid2606.200301]
18. Noce, A.; Marrone, G.; Di Daniele, F.; Ottaviani, E.; Wilson Jones, G.; Bernini, R.; Romani, A.; Rovella, V. Impact of Gut Microbiota Composition on Onset and Progression of Chronic Non-Communicable Diseases. Nutrients; 2019; 11, 1073. [DOI: https://dx.doi.org/10.3390/nu11051073]
19. Zakerska-Banaszak, O.; Tomczak, H.; Gabryel, M.; Baturo, A.; Wolko, L.; Michalak, M.; Malinska, N.; Mankowska-Wierzbicka, D.; Eder, P.; Dobrowolska, A. et al. Dysbiosis of Gut Microbiota in Polish Patients with Ulcerative Colitis: A Pilot Study. Sci. Rep.; 2021; 11, 2166. [DOI: https://dx.doi.org/10.1038/s41598-021-81628-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33495479]
20. Sirichoat, A.; Sankuntaw, N.; Engchanil, C.; Buppasiri, P.; Faksri, K.; Namwat, W.; Chantratita, W.; Lulitanond, V. Comparison of Different Hypervariable Regions of 16S RRNA for Taxonomic Profiling of Vaginal Microbiota Using Next-Generation Sequencing. Arch. Microbiol.; 2021; 203, pp. 1159-1166. [DOI: https://dx.doi.org/10.1007/s00203-020-02114-4]
21. Gao, F.; Guo, R.; Ma, Q.; Li, Y.; Wang, W.; Fan, Y.; Ju, Y.; Zhao, B.; Gao, Y.; Qian, L. et al. Stressful Events Induce Long-Term Gut Microbiota Dysbiosis and Associated Post-Traumatic Stress Symptoms in Healthcare Workers Fighting against COVID-19. J. Affect. Disord.; 2022; 303, pp. 187-195. [DOI: https://dx.doi.org/10.1016/j.jad.2022.02.024]
22. Mazzarelli, A.; Giancola, M.L.; Farina, A.; Marchioni, L.; Rueca, M.; Gruber, C.E.M.; Bartolini, B.; Ascoli Bartoli, T.; Maffongelli, G.; Capobianchi, M.R. et al. 16S RRNA Gene Sequencing of Rectal Swab in Patients Affected by COVID-19. PLoS ONE; 2021; 16, e0247041. [DOI: https://dx.doi.org/10.1371/journal.pone.0247041] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33596245]
23. Nardelli, C.; Gentile, I.; Setaro, M.; Di Domenico, C.; Pinchera, B.; Buonomo, A.R.; Zappulo, E.; Scotto, R.; Scaglione, G.L.; Castaldo, G. et al. Nasopharyngeal Microbiome Signature in COVID-19 Positive Patients: Can We Definitively Get a Role to Fusobacterium Periodonticum?. Front. Cell Infect. Microbiol.; 2021; 11, 625581. [DOI: https://dx.doi.org/10.3389/fcimb.2021.625581] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33659220]
24. Uehara, O.; Abiko, Y.; Nagasawa, T.; Morikawa, T.; Hiraki, D.; Harada, F.; Kawano, Y.; Toraya, S.; Matsuoka, H.; Paudel, D. et al. Alterations in the Oral Microbiome of Individuals with a Healthy Oral Environment Following COVID-19 Vaccination. BMC Oral Health; 2022; 22, 50. [DOI: https://dx.doi.org/10.1186/s12903-022-02093-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35241064]
25. Matijašić, M.; Meštrović, T.; Čipčić Paljetak, H.; Perić, M.; Barešić, A.; Verbanac, D. Gut Microbiota beyond Bacteria—Mycobiome, Virome, Archaeome, and Eukaryotic Parasites in IBD. IJMS; 2020; 21, 2668. [DOI: https://dx.doi.org/10.3390/ijms21082668] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32290414]
26. Hall, A.B.; Tolonen, A.C.; Xavier, R.J. Human Genetic Variation and the Gut Microbiome in Disease. Nat. Rev. Genet.; 2017; 18, pp. 690-699. [DOI: https://dx.doi.org/10.1038/nrg.2017.63] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28824167]
27. Cryan, J.F.; O’Riordan, K.J.; Cowan, C.S.M.; Sandhu, K.V.; Bastiaanssen, T.F.S.; Boehme, M.; Codagnone, M.G.; Cussotto, S.; Fulling, C.; Golubeva, A.V. et al. The Microbiota-Gut-Brain Axis. Physiol. Rev.; 2019; 99, pp. 1877-2013. [DOI: https://dx.doi.org/10.1152/physrev.00018.2018]
28. Belkaid, Y.; Hand, T.W. Role of the Microbiota in Immunity and Inflammation. Cell; 2014; 157, pp. 121-141. [DOI: https://dx.doi.org/10.1016/j.cell.2014.03.011]
29. Yoo, B.B.; Mazmanian, S.K. The Enteric Network: Interactions between the Immune and Nervous Systems of the Gut. Immunity; 2017; 46, pp. 910-926. [DOI: https://dx.doi.org/10.1016/j.immuni.2017.05.011]
30. Budden, K.F.; Gellatly, S.L.; Wood, D.L.A.; Cooper, M.A.; Morrison, M.; Hugenholtz, P.; Hansbro, P.M. Emerging Pathogenic Links between Microbiota and the Gut–Lung Axis. Nat. Rev. Microbiol.; 2017; 15, pp. 55-63. [DOI: https://dx.doi.org/10.1038/nrmicro.2016.142] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27694885]
31. DeGruttola, A.K.; Low, D.; Mizoguchi, A.; Mizoguchi, E. Current Understanding of Dysbiosis in Disease in Human and Animal Models. Inflamm. Bowel Dis.; 2016; 22, pp. 1137-1150. [DOI: https://dx.doi.org/10.1097/MIB.0000000000000750]
32. Cardinale, V.; Capurso, G.; Ianiro, G.; Gasbarrini, A.; Arcidiacono, P.G.; Alvaro, D. Intestinal Permeability Changes with Bacterial Translocation as Key Events Modulating Systemic Host Immune Response to SARS-CoV-2: A Working Hypothesis. Dig. Liver Dis.; 2020; 52, pp. 1383-1389. [DOI: https://dx.doi.org/10.1016/j.dld.2020.09.009]
33. Dickson, R.P.; Singer, B.H.; Newstead, M.W.; Falkowski, N.R.; Erb-Downward, J.R.; Standiford, T.J.; Huffnagle, G.B. Enrichment of the Lung Microbiome with Gut Bacteria in Sepsis and the Acute Respiratory Distress Syndrome. Nat. Microbiol.; 2016; 1, 16113. [DOI: https://dx.doi.org/10.1038/nmicrobiol.2016.113] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27670109]
34. Burcelin, R.; Garidou, L.; Pomié, C. Immuno-Microbiota Cross and Talk: The New Paradigm of Metabolic Diseases. Semin. Immunol.; 2012; 24, pp. 67-74. [DOI: https://dx.doi.org/10.1016/j.smim.2011.11.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22265028]
35. Jin, X.; Lian, J.-S.; Hu, J.-H.; Gao, J.; Zheng, L.; Zhang, Y.-M.; Hao, S.-R.; Jia, H.-Y.; Cai, H.; Zhang, X.-L. et al. Epidemiological, Clinical and Virological Characteristics of 74 Cases of Coronavirus-Infected Disease 2019 (COVID-19) with Gastrointestinal Symptoms. Gut; 2020; 69, pp. 1002-1009. [DOI: https://dx.doi.org/10.1136/gutjnl-2020-320926] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32213556]
36. Hanada, S.; Pirzadeh, M.; Carver, K.Y.; Deng, J.C. Respiratory Viral Infection-Induced Microbiome Alterations and Secondary Bacterial Pneumonia. Front. Immunol.; 2018; 9, 2640. [DOI: https://dx.doi.org/10.3389/fimmu.2018.02640] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30505304]
37. Zuo, T.; Zhang, F.; Lui, G.C.Y.; Yeoh, Y.K.; Li, A.Y.L.; Zhan, H.; Wan, Y.; Chung, A.C.K.; Cheung, C.P.; Chen, N. et al. Alterations in Gut Microbiota of Patients With COVID-19 During Time of Hospitalization. Gastroenterology; 2020; 159, pp. 944-955.e8. [DOI: https://dx.doi.org/10.1053/j.gastro.2020.05.048]
38. Xu, R.; Lu, R.; Zhang, T.; Wu, Q.; Cai, W.; Han, X.; Wan, Z.; Jin, X.; Zhang, Z.; Zhang, C. Temporal Association between Human Upper Respiratory and Gut Bacterial Microbiomes during the Course of COVID-19 in Adults. Commun. Biol.; 2021; 4, 240. [DOI: https://dx.doi.org/10.1038/s42003-021-01796-w]
39. Liu, T.F.D.; Philippou, E.; Kolokotroni, O.; Siakallis, G.; Rahima, K.; Constantinou, C. Gut and Airway Microbiota and Their Role in COVID-19 Infection and Pathogenesis: A Scoping Review. Infection; 2021; 50, pp. 815-847. [DOI: https://dx.doi.org/10.1007/s15010-021-01715-5]
40. Yeoh, Y.K.; Zuo, T.; Lui, G.C.-Y.; Zhang, F.; Liu, Q.; Li, A.Y.; Chung, A.C.; Cheung, C.P.; Tso, E.Y.; Fung, K.S. et al. Gut Microbiota Composition Reflects Disease Severity and Dysfunctional Immune Responses in Patients with COVID-19. Gut; 2021; 70, pp. 698-706. [DOI: https://dx.doi.org/10.1136/gutjnl-2020-323020]
41. Gu, S.; Chen, Y.; Wu, Z.; Chen, Y.; Gao, H.; Lv, L.; Guo, F.; Zhang, X.; Luo, R.; Huang, C. et al. Alterations of the Gut Microbiota in Patients With Coronavirus Disease 2019 or H1N1 Influenza. Clin. Infect. Dis.; 2020; 71, pp. 2669-2678. [DOI: https://dx.doi.org/10.1093/cid/ciaa709]
42. Zuo, T.; Liu, Q.; Zhang, F.; Lui, G.C.-Y.; Tso, E.Y.; Yeoh, Y.K.; Chen, Z.; Boon, S.S.; Chan, F.K.; Chan, P.K. et al. Depicting SARS-CoV-2 Faecal Viral Activity in Association with Gut Microbiota Composition in Patients with COVID-19. Gut; 2020; 70, pp. 276-284. [DOI: https://dx.doi.org/10.1136/gutjnl-2020-322294] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32690600]
43. Han, H.; Ma, Q.; Li, C.; Liu, R.; Zhao, L.; Wang, W.; Zhang, P.; Liu, X.; Gao, G.; Liu, F. et al. Profiling Serum Cytokines in COVID-19 Patients Reveals IL-6 and IL-10 Are Disease Severity Predictors. Emerg. Microbes Infect.; 2020; 9, pp. 1123-1130. [DOI: https://dx.doi.org/10.1080/22221751.2020.1770129] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32475230]
44. Louis, P.; Flint, H.J. Diversity, Metabolism and Microbial Ecology of Butyrate-Producing Bacteria from the Human Large Intestine. FEMS Microbiol. Lett.; 2009; 294, pp. 1-8. [DOI: https://dx.doi.org/10.1111/j.1574-6968.2009.01514.x]
45. Rivière, A.; Selak, M.; Lantin, D.; Leroy, F.; De Vuyst, L. Bifidobacteria and Butyrate-Producing Colon Bacteria: Importance and Strategies for Their Stimulation in the Human Gut. Front. Microbiol.; 2016; 7, 979. [DOI: https://dx.doi.org/10.3389/fmicb.2016.00979] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27446020]
46. Derrien, M.; Vaughan, E.E.; Plugge, C.M.; de Vos, W.M. Akkermansia Muciniphila Gen. Nov., Sp. Nov., a Human Intestinal Mucin-Degrading Bacterium. Int. J. Syst. Evol. Microbiol.; 2004; 54, pp. 1469-1476. [DOI: https://dx.doi.org/10.1099/ijs.0.02873-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15388697]
47. Naito, Y.; Uchiyama, K.; Takagi, T. A Next-Generation Beneficial Microbe: Akkermansia muciniphila. J. Clin. Biochem. Nutr.; 2018; 63, pp. 33-35. [DOI: https://dx.doi.org/10.3164/jcbn.18-57]
48. Mańkowska-Wierzbicka, D.; Stelmach-Mardas, M.; Gabryel, M.; Tomczak, H.; Skrzypczak-Zielińska, M.; Zakerska-Banaszak, O.; Sowińska, A.; Mahadea, D.; Baturo, A.; Wolko, Ł. et al. The Effectiveness of Multi-Session FMT Treatment in Active Ulcerative Colitis Patients: A Pilot Study. Biomedicines; 2020; 8, 268. [DOI: https://dx.doi.org/10.3390/biomedicines8080268]
49. Pascoal, L.B.; Rodrigues, P.B.; Genaro, L.M.; dos Santos Pereira Gomes, A.B.; Toledo-Teixeira, D.A.; Parise, P.L.; Bispo-Dos-Santos, K.; Simeoni, C.L.; Guimarães, P.V.; Buscaratti, L.I. et al. Microbiota-Derived Short-Chain Fatty Acids Do Not Interfere with SARS-CoV-2 Infection of Human Colonic Samples. Gut Microbes; 2021; 13, pp. 1-9. [DOI: https://dx.doi.org/10.1080/19490976.2021.1874740]
50. Singh, N.; Gurav, A.; Sivaprakasam, S.; Brady, E.; Padia, R.; Shi, H.; Thangaraju, M.; Prasad, P.D.; Manicassamy, S.; Munn, D.H. et al. Activation of Gpr109a, Receptor for Niacin and the Commensal Metabolite Butyrate, Suppresses Colonic Inflammation and Carcinogenesis. Immunity; 2014; 40, pp. 128-139. [DOI: https://dx.doi.org/10.1016/j.immuni.2013.12.007]
51. Peng, L.; Li, Z.-R.; Green, R.S.; Holzman, I.R.; Lin, J. Butyrate Enhances the Intestinal Barrier by Facilitating Tight Junction Assembly via Activation of AMP-Activated Protein Kinase in Caco-2 Cell Monolayers. J. Nutr.; 2009; 139, pp. 1619-1625. [DOI: https://dx.doi.org/10.3945/jn.109.104638]
52. Saeedi, B.J.; Kao, D.J.; Kitzenberg, D.A.; Dobrinskikh, E.; Schwisow, K.D.; Masterson, J.C.; Kendrick, A.A.; Kelly, C.J.; Bayless, A.J.; Kominsky, D.J. et al. HIF-Dependent Regulation of Claudin-1 Is Central to Intestinal Epithelial Tight Junction Integrity. Mol. Biol. Cell; 2015; 26, pp. 2252-2262. [DOI: https://dx.doi.org/10.1091/mbc.E14-07-1194] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25904334]
53. Gambardella, A.; Licata, G.; Calabrese, G.; De Rosa, A.; Pagliuca, F.; Alfano, R.; Argenziano, G. CMV Infection: A Clinical Challenge in Biological Therapy? The Case of Asymptomatic Patients with Persistent Positive Immunoglobulin M Anti-CMV Treated with Secukinumab. Psoriasis Targets Ther.; 2020; 10, pp. 57-60. [DOI: https://dx.doi.org/10.2147/PTT.S284701] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33282718]
54. Liévin-Le Moal, V.; Servin, A.L. The Front Line of Enteric Host Defense against Unwelcome Intrusion of Harmful Microorganisms: Mucins, Antimicrobial Peptides, and Microbiota. Clin. Microbiol. Rev.; 2006; 19, pp. 315-337. [DOI: https://dx.doi.org/10.1128/CMR.19.2.315-337.2006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16614252]
55. Perlot, T.; Penninger, J.M. ACE2—From the Renin–Angiotensin System to Gut Microbiota and Malnutrition. Microbes Infect.; 2013; 15, pp. 866-873. [DOI: https://dx.doi.org/10.1016/j.micinf.2013.08.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23962453]
56. Heine, W.; Radke, M.; Wutzke, K.-D. The Significance of Tryptophan in Human Nutrition. Amino Acids; 1995; 9, pp. 91-205. [DOI: https://dx.doi.org/10.1007/BF00805951]
57. Ninomiya, S.; Nakamura, N.; Nakamura, H.; Mizutani, T.; Kaneda, Y.; Yamaguchi, K.; Matsumoto, T.; Kitagawa, J.; Kanemura, N.; Shiraki, M. et al. Low Levels of Serum Tryptophan Underlie Skeletal Muscle Atrophy. Nutrients; 2020; 12, 978. [DOI: https://dx.doi.org/10.3390/nu12040978]
58. Pickard, J.M.; Zeng, M.Y.; Caruso, R.; Núñez, G. Gut Microbiota: Role in Pathogen Colonization, Immune Responses, and Inflammatory Disease. Immunol. Rev.; 2017; 279, pp. 70-89. [DOI: https://dx.doi.org/10.1111/imr.12567]
59. Di Pierro, F. Gut Microbiota Parameters Potentially Useful in Clinical Perspective. Microorganisms; 2021; 9, 2402. [DOI: https://dx.doi.org/10.3390/microorganisms9112402]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
It is crucial to consider the importance of the microbiome and the gut–lung axis in the context of SARS-CoV-2 infection. This pilot study examined the fecal microbial composition of patients with COVID-19 following a 3-month recovery. Using for the first time metagenomic analysis based on all hypervariable regions (V1-V9) of the 16S rRNA gene, we have identified 561 microbial species; however, 17 were specific only for the COVID-19 group (n = 8). The patients’ cohorts revealed significantly greater alpha diversity of the gut microbiota compared to healthy controls (n = 14). This finding has been demonstrated by operational taxonomic units (OTUs) richness (p < 0.001) and Chao1 index (p < 0.01). The abundance of the phylum Verrucomicrobia was 30 times higher in COVID-19 patients compared to healthy subjects. Accordingly, this disproportion was also noted at other taxonomic levels: in the class Verrucomicrobiae, the family Verrucomicrobiaceae, and the genus Akkermansia. Elevated pathobionts such as Escherichia coli, Bilophila wadsworthia, and Parabacteroides distasonis were found in COVID-19 patients. Considering the gut microbiota’s ability to disturb the immune response, our findings suggest the importance of the enteric microbiota in the course of SARS-CoV-2 infection. This pilot study shows that the composition of the microbial community may not be fully restored in individuals with SARS-CoV-2 following a 3-month recovery.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details






1 Department of Gastroenterology, Dietetics and Internal Diseases, Poznan University of Medical Sciences, 60-355 Poznan, Poland
2 Institute of Human Genetics, Polish Academy of Sciences, 60-479 Poznan, Poland