Gut Microbiome Alterations in COVID-19.

Since the outset of the coronavirus disease 2019 (COVID-19) pandemic, the gut microbiome in COVID-19 has garnered substantial interest, given its significant roles in human health and pathophysiology. Accumulating evidence is unveiling that the gut microbiome is broadly altered in COVID-19, including the bacterial microbiome, mycobiome, and virome. Overall, the gut microbial ecological network is significantly weakened and becomes sparse in patients with COVID-19, together with a decrease in gut microbiome diversity. Beyond the existence of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), the gut microbiome of patients with COVID-19 is also characterized by enrichment of opportunistic bacteria, fungi, and eukaryotic viruses, which are also associated with disease severity and presentation. Meanwhile, a multitude of symbiotic bacteria and bacteriophages are decreased in abundance in patients with COVID-19. Such gut microbiome features persist in a significant subset of patients with COVID-19 even after disease resolution, coinciding with 'long COVID' (also known as post-acute sequelae of COVID-19). The broadly-altered gut microbiome is largely a consequence of SARS-CoV-2infection and its downstream detrimental effects on the systemic host immunity and the gut milieu. The impaired host immunity and distorted gut microbial ecology, particularly loss of low-abundance beneficial bacteria and blooms of opportunistic fungi including Candida, may hinder the reassembly of the gut microbiome post COVID-19. Future investigation is necessary to fully understand the role of the gut microbiome in host immunity against SARS-CoV-2 infection, as well as the long-term effect of COVID-19 on the gut microbiome in relation to the host health after the pandemic.

Introduction

The ongoing global pandemic of coronavirus disease 2019 (COVID-19) is a disease caused by the RNA virus severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), primarily infecting the respiratory tract and resulting in various symptoms at various severity levels in patients after infection [1]. Around 5%–33% COVID-19 patients had gastrointestinal (GI) symptoms, including diarrhea, nausea, and vomiting [2], [3], [4]. Several studies have detected SARS-CoV-2 in stool samples and anal swabs [5], [6], suggesting that the digestive tract might be an extra-pulmonary site for SARS-CoV-2 infection. Although most cases of COVID-19 are mild, disease can be severe and result in hospitalization, respiratory failure, or death [1]. Such remarkable differences in individual’s presentations and symptoms of COVID-19 arise from the heterogeneous immune statuses and responses against SARS-CoV-2 infection [7], [8], [9]. The GI tract is the largest immune organ in humans, playing critical roles in combating infections of pathogens [10]. Living inside the gut of humans are trillions of microorganisms — bacteria, fungi, viruses, and other life forms that are collectively known as the microbiome — regulating host immunity [11]. As of now, accumulating evidence suggests that the gut microbiome ecology is broadly altered in patients with COVID-19 and that the gut microbiome configurations are associated with immune responses and disease presentations in COVID-19 [12], [13], [14], [15]. The SARS-CoV-2 infection course is crucial for the alterations in the ecology and dynamics of human gut microbiome, in both the short term and long term, which in return influence the human host’s health. Moreover, the presence of active SARS-CoV-2 virus in the gut and altered ecology of the gut microbiome may lead to an unfavorable gut milieu, which facilitates opportunistic bloom of certain fungi and pathogenic bacteria, further hindering the community assembly and function of the gut microbiome, as well as weakening the host immunity [9], [12], [13], [14], [15], [16]. Herein, we summarize the impact of COVID-19 on the human gut microbiome in association with disease phenotypes, from the perspective of the gut microbial ecology, including that of bacteria, fungi, and viruses.

The gut bacterial microbiome in COVID-19

Compositional changes of the gut bacterial microbiome

Our early study showed that the gut bacterial microbiome of patients with COVID-19 was significantly altered compared with healthy controls, characterized by depletion of beneficial commensals and enrichment of opportunistic pathogens in the gut [15]. The feces of COVID-19 patients was enriched for opportunistic pathogens known to cause bacteremia, including Clostridium hathewayi, Actinomyces viscosus, and Bacteroides nordii (Figure 1; Table 1) [15], as a secondary infection/bloom post onset of COVID-19 due to disrupted gut microbial ecology and colonization resistance [17], [18]. The patients treated with antibiotics at hospitalization displayed a further depletion of bacterial species, particularly symbionts beneficial to host immunity including Faecalibacterium prausnitzii, Lachnospiraceae bacterium 5_1_63FAA, Eubacterium rectale, Ruminococcus obeum, and Dorea formicigenerans [15]. Such alterations in the bacterial microbiome ecology persisted over the disease course of COVID-19 and even after clearance of SARS-CoV-2 from the respiratory tract [15]. Consistently, another study also showed a similar pattern of gut microbiome dysbiosis in COVID-19 patients [19]. The abundance of butyrate-producing bacteria, such as Faecalibacterium prausnitzii, Clostridium butyricum, Clostridium leptum, and Eubacterium rectale, was significantly decreased in patients with COVID-19 compared to controls [19]. In contrast, the abundance of the common opportunistic pathogens Enterobacteriaceae and Enterococcus was significantly increased in patients with COVID-19 compared to controls [19]. At the genus level, the genera Streptococcus, Rothia, Veillonella, and Actinomyces (all opportunistic pathogens) were enriched in the feces of COVID-19 patients, whereas the genera Romboutsia, Faecalibacterium, and Fusicatenibacter were enriched in the feces of healthy controls [20]. An ecological network analysis revealed significant positive correlations across COVID-19-enriched genera [20], indicating co-expansion of opportunistic bacteria dominating the ecological network of the gut microbiome due to SARS-CoV-2 infection. A high baseline abundance of opportunistic bacteria Coprobacillus, Clostridium ramosum, and Clostridium hathewayi in patients’ feces at hospitalization was associated with a more severe COVID-19 disease course, whereas the anti-inflammatory bacterium Faecalibacterium prausnitzii showed an inverse correlation [15], suggesting baseline gut microbiome calibration of host immunity, thereby affecting disease response upon SARS-CoV-2 infection.

Figure 1

Alterations in the gut bacterial, fungal, and viral microbiome in patients with COVID-19

The gut bacterial microbiome in COVID-19 is characterized by decreased diversity and richness, and persistent bacterial microbiome dysbiosis even after disease resolution. The gut mycobiome in COVID-19 is characterized by increased fecal fungal load and increased beta-diversity (more heterogeneous), and it is unstable over time and also persistently altered after disease resolution. SARS-CoV-2 shows infectivity in the gut. Delayed SARS-CoV-2 viral shedding and persistent gut virome dysbiosis are both present after disease resolution. The gastrointestinal tract epithelial barrier is impaired in a subset of COVID-19 patients. The figure is created with BioRender.com.

Table 1

Gut microbes with altered abundance in COVID-19

Kingdom	Species	Change	Feature	Refs.
Bacteria	Actinomyces viscosus	↑	Opportunistic pathogens known to cause bacteremia	[15]
	Clostridium hathewayi	↑	Associated with human infection and bacteremia	[15]
	Bacteroides nordii	↑	Opportunistic pathogens to cause bacteremia	[15]
	Eubacterium ventriosum	↓	Butyrate-producing	[15]
	Dorea formicigenerans	↓	Butyrate-producing	[15]
	Faecalibacterium prausnitzii	↓	Butyrate-producing	[15], [19]
	Eubacterium rectale	↓	Butyrate-producing	[15], [19]
	Ruminococcus obeum	↓	Preventing increases in intestinal permeability	[15]
	Lachnospiraceae bacterium 5_1_63FAA	↓	Butyrate-producing	[15]
	Enterococcus faecium	↓	Opportunistic pathogens	[15]
	Clostridium ramosum	↓	Immune-regulatory characteristics	[15]
Fungi	Candida albicans	↑	Opportunistic infection and impairing microbiome assembly	[15]
	Candida auris	↑	Healthcare-associated invasive infection	[14]
	Aspergillus flavus	↑	Associated with respiratory symptoms	[14]
	Escheruchia virus	↑	Bacteriophage	[14]
Virus/phage	Edafosvirus	↑	Unknown host from the environment	[13]
	Solumvirus	↑	Unknown host from the environment	[13]
	Bandra megavirus	↑	Unknown host from the environment	[13]
	Streptococcus virus	↑	Bacteriophage	[13]
	Paguma larvata circovirus	↑	Animal host	[13]
	Microcystis virus	↑	Bacteriophage	[13]
	Ralstonia phage	↑	Bacteriophage	[13]
	EBPR podovirus	↑	Bacteriophage	[13]
	Dishui lake phycodnavirus	↑	Algae host	[13]
	Homavirus	↑	Unknown host from the environment	[13]
	Cannes 8 virus	↑	Amoeba host	[13]
	Pigeonpox virus	↑	Columba host	[13]
	Sucra jujuba nucleopolyhedrovirus	↑	Geometrid host	[13]
	Phormidium phage	↑	Bacteriophage	[13]
	Lactococcus phage	↑	Bacteriophage	[13]
	Enterobacter phage	↑	Bacteriophage	[13]
	Columbid alphaherpesvirus	↑	Columba host	[13]
	Faustovirus	↑	Amoeba host	[13]
	Burkholderia virus	↓	Bacteriophage	[13]
	Lysinibacillus phage	↓	Bacteriophage	[13]
	Bdellovibrio phage	↓	Bacteriophage	[13]
	Citrobacter phage	↓	Bacteriophage	[13]
	Ranid herpesvirus	↓	Animal host	[13]
	Dunaliella viridis virus	↓	Algae host	[13]
	Marseillevirus	↓	Amoeba host	[13]
	Raccoonpox virus	↓	Animal host	[13]
	Kaumoebavirus	↓	Amoeba host	[13]
	Acinetobacter virus	↓	Bacteriophage	[13]
	Halomonas phage	↓	Bacteriophage	[13]
	Mycobacterium phage BGlluviae	↓	Bacteriophage	[13]
	Golden marseillevirus	↓	Amoeba host	[13]
	Silicibacter phage	↓	Bacteriophage	[13]
	Tokyovirus	↓	Amoeba host	[13]
	Sulfolobales virus YNP2	↓	Archaea host	[13]
	Nocardia phage	↓	Bacteriophage	[13]
	Alteromonas phage	↓	Bacteriophage	[13]
	Hyperthermophilic archaeal virus 2	↓	Archaea host	[13]
	Crassphage	↓	Bacteriophage	[13]
	Nodularia phage	↓	Bacteriophage	[13]
	Hamiltonella virus	↓	Bacteriophage	[13]
	Hudisavirus	↓	Human host	[13]
	Mesorhizobium phage	↓	Bacteriophage	[13]
	Uncultured crAssphage	↓	Bacteriophage	[13]
	Halocynthia phage	↓	Bacteriophage	[13]

Note: “↑”, abundance increased; “↓”, abundance decreased.

Evidence has been accumulating that a substantial number of COVID-19 patients experienced systemic and/or organ-specific afflictions during follow-up after disease resolution, including fatigue, muscle weakness, sleep difficulties, anxiety, depression, diarrhea, and poor glycemic controls [21], [22], [23], [24], a phenomenon known as ‘long COVID’. Interestingly, the GI tract is also affected in a long term in COVID-19, as demonstrated by a prolonged shedding of viral RNA in stool specimens up to 42 days and the presence of SARS-CoV-2 virus in the gut epithelium up to 90 days after disease resolution in some patients [25], [26]. Concordantly, long-lasting gut microbiome dysbiosis is also consistently observed in subjects recovered from COVID-19 [12], [15], [27], [28], implying that gut microbiome is closely linked to host health in a post-COVID-19 age.

Relationship between gut microbiome changes, infection, and host immunity

In a six-month follow-up study on the gut microbiome of patients with COVID-19, significant decreases in the richness (Chao1 index) of gut microbiome were observed across the acute, convalescence, and post-convalescence phases of COVID-19 [27]. In addition, COVID-19 patients had a significantly reduced gut bacterial diversity [20], [29]. Microbial diversity is a critical determinant of microbial ecosystem stability [30]. Stable ecosystems provide colonization resistance to opportunistic pathogens [31]. Therefore, the reduction in gut microbiota diversity and richness may somewhat contribute to the expansion of opportunistic bacteria and have long-term impact in patients with COVID-19 [32]. Concordantly, a critical proportion of patients with COVID-19 also experienced persistent symptoms following disease resolution and hospital discharge, known as ‘long COVID-19’ [21], [24]. Patients with lower post-convalescence bacterial microbiome richness had higher levels of COVID-19 severity (worse pulmonary functions) and blood C-reactive protein (CRP) during the acute phase [27], suggesting a relationship between gut dysbiosis and hyper-inflammatory response in COVID-19. A more recent study also showed that the gut microbiome ecology was stratified well with COVID-19 severity, as demonstrated in the principal component analysis (PCA) visualization that the gut microbiome communities followed a continuum along the mild, moderate, severe, and critical gradients of COVID-19 severity [12]. Moreover, the gut microbiota composition was correlated with plasma concentrations of inflammatory cytokines and blood parameters, such as CRP, lactate dehydrogenase, aspartate aminotransferase, and gamma-glutamyl transferase [12]. These data together suggest that SARS-CoV-2 infection may cause immune-pathophysiological changes in the human host, including the gut, resulting in gradual changes in the gut microbial ecology in relation to illness severity. In favor of this hypothesis, a recent proof-of-principle study in a mouse model of COVID-19 demonstrated that SARS-CoV-2 infection elicited immune/infection-related gene expression in the gut epithelial cells, leading to a change in the gut milieu where the microbiota were affected [33]. Following that, we found that the SARS-CoV-2 activity in the gut might be a prominent factor in shaping the gut microbiome composition [34]. Patients with high SARS-CoV-2 infectivity in the gut displayed a high abundance of the bacterial species Collinsella aerofaciens, Collinsella tanakaei, Streptococcus infantis, and Morganella morganii (Figure 1; Table 1), as well as a high functional capacity for nucleotide de novo biosynthesis, amino acid biosynthesis, and glycolysis [34]. However, patients with low-to-none SARS-CoV-2 infectivity in the gut displayed a high abundance of short-chain fatty acid (SCFA)-producing bacteria, Alistipes onderdonkii, Parabacteroides merdae, Bacteroides stercoris, and Lachnospiraceae bacterium 1_1_57FAA [34]. Among them, Alistipes onderdonkii was a bacterial species, the abundance of which also showed a negative correlation with COVID-19 severity [15]. Interestingly, Alistipes species are indole-positive, involved in the serotonin precursor tryptophan metabolism and in maintaining gut immune homeostasis [35], [36]. This is later validated in animals that tryptophan metabolism in the gut was altered as a result of SARS-CoV-2 infection [28]. In addition, the hyper-inflammatory response of COVID-19 patients was associated with disrupted gut permeability and microbial translocation [16], [37]. The amount of fecal calprotectin, a marker of intestinal inflammation as a consequence of translocation of granulocytes and monocytes/macrophages into the gut lumen, was elevated in the feces of patients with COVID-19 [38], indicating immune dysfunction of the gut and altered gut niche in COVID-19 patients. Taken together, the compositional alterations in the gut microbiome of COVID-19 patients are likely the result of host immune responses and altered gut milieu during SARS-CoV-2 infection.

To validate that SARS-CoV-2 infection is the driving force of the alterations in the gut microbiome ecology, Sokol et al. used a non-human primate model (rhesus macaques and cynomolgus macaques) challenged with SARS-CoV-2 and subsequently analyzed the impact of SARS-CoV-2 infection on dynamic changes of the gut microbiome [28]. Strikingly, the gut microbiome gradually changed from day 0 until day 13 post infection, at which time point the gut microbiome was most different from the baseline microbiome in terms of fecal microbial community structure [28]. This result indicates that SARS-CoV-2 infection virtually induces alterations in the gut microbiome ecology in COVID-19. To delineate the dynamic changes caused by SARS-CoV-2 infection, the authors compared the composition at each time point post infection with time points before infection. Consistent with the findings in humans [15], [20], the abundance of opportunistic bacteria from the Proteobacteria phylum was increased, whilst the abundance of beneficial members from the Firmicutes phylum (especially those from the Ruminococcaceae and Lachnospiraceae families) was decreased after SARS-CoV-2 infection [28]. Although some alterations in the gut microbiome were resolved at later time points, certain perturbations persisted even after disease resolution [28], analogous to the observations in humans [15], [20]. This finding further addresses that SARS-CoV-2 infection may have long-lasting impact on the gut microbiome ecology. Ecological network analysis of the bacterial–bacterial interactions in the gut microbiome of macaques before vs. after SARS-CoV-2 infection revealed a sparse, atrophied bacterial microbiome ecological network after SARS-CoV-2 infection compared to a dense, interconnected network before SARS-CoV-2 infection [28]. The gut microbiome ecological network reflects the complex interplay of microbial communities [39]. In a steady state, the gut microbiome exhibits a dense, intricate microbial ecological network, whereas under gut inflammation conditions, such as inflammatory bowel disease (IBD) and Clostridioides difficile infection (CDI), it manifests a significant sparse one [11], [40], [41]. The significantly weakened ecological microbial network after SARS-CoV-2 infection both in humans and macaques implies that SARS-CoV-2 infection may induce host inflammatory responses resulting in disrupted gut microbiome ecology.

Studies have shown that intestinal microbiota can affect viral replication and systemic pathogenesis [42], [43], [44], [45]. Depletion of the intestinal microbiota in mice by antibiotics rendered the mice less susceptible to poliovirus disease and supported minimal viral replication in the intestines of mice [43]. Exposure of poliovirus to bacteria enhanced host–cell association and infection, since poliovirus binds to lipopolysaccharide [43]. The pathology of reovirus (an unrelated enteric virus) infection was also more severe in the presence of intestinal microbes [43]. In addition, antibiotics prevented persistent murine norovirus infection, which was reversed by replenishment of the bacterial microbiota [45]. In parallel, enteric bacteria were also found to promote human and mouse norovirus infection of B cells [44]. These studies together suggest that gut microbes influence virus infection and that viruses may exploit intestinal microbes for replication and transmission. That being said, the gut microbiome may play a role in SARS-CoV-2 susceptibility and infectivity, which remains to be verified in future studies.

COVID-19 is essentially a lung disease, and it has been established that gut can affect lung through the gut–lung axis [46], [47]. Beyond the local immune regulation by the gut microbiota, the far-reaching immune impact of gut microbiota is also well recognized, especially on the pulmonary immune system [48]. SCFAs, a group of prototypic metabolites produced by gut bacteria, translocate across the intestinal barrier, reach the systemic circulation, and modulate the lung immune response [49], [50]. They are mainly produced by bacterial degradation and fermentation of dietary fibers, acting as signaling molecules in the lungs on resident antigen-presenting cells to attenuate the inflammatory and allergic responses [49], [51], [52]. Decreases in the abundance of SCFA-producing bacteria observed in the gut microbiota of patients with COVID-19 [12], [15], [34] may represent one of the critical mechanisms contributing to the gut–lung crosstalk and thereby disease severity in COVID-19.

Angiotensin-converting enzyme 2 and the gut microbiome

Studies have provided direct evidence that angiotensin-converting enzyme 2 (ACE2) is the binding site of SARS-CoV-2 for host entry [53], [54]. ACE2 is highly expressed in the respiratory tract and the intestines, especially in nasal epithelial cells and colonocytes of humans [55]. ACE2 has also been demonstrated to regulate amino acid transport, expression of antimicrobial peptides, microbial ecology, and inflammation in the gut [56]. These lines of evidence underscore an interplay between ACE2 expression, SARS-CoV-2 infection, and the gut microbiome in the host. Bacterial species from the Bacteroidetes phylum were shown to down-regulate ACE2 expression, while species from the Firmicutes phylum displayed variable effects in modulating ACE2 expression in the murine colon [57]. Interestingly, our study in the gut microbiome of COVID-19 patients showed that the fecal abundance of the Bacteroidetes species, Alistipes onderdonkii and Bacteroides ovatus, was inversely correlated with COVID-19 severity, and the abundance of 4 Bacteroidetes species, Bacteroides dorei, Bacteroides thetaiotaomicron, Bacteroides massiliensis, and Bacteroides ovatus, showed inverse correlation with the fecal viral load of SARS-CoV-2 [15], [34]. Amongst these Bacteroidetes species, Bacteroides dorei was previously shown to inhibit colonic ACE2 expression [57] and to calibrate host immune response [58], [59]. Intriguingly, subjects with pre-existing chronic diseases (such as diabetes mellitus, hypertension, obesity, and coronary artery disease) were characterized by a low abundance of Bacteroides species and had the highest COVID-19 mortality and morbidity [60], [61], [62]. Collectively, these data imply a sophisticated quaternary relationship between SARS-CoV-2, gut microbiome, ACE2 expression, and host immunity, underlying the varying anti-SARS-CoV-2 immune responses and thereby disease severity in the host.

Functionality changes of the gut microbiome in COVID-19

Beyond the compositional changes, the functionality of the gut microbiome is also changed in COVID-19. The fecal microbiome of patients with high SARS-CoV-2 infectivity in the gut was enriched for functional pathways involved in amino acid biosynthesis (L-lysine biosynthesis II and superpathway of L-serine and glycine biogenesis), carbohydrate metabolism (glycolysis II from fructose-6-phosphate), and nucleotide metabolism (superpathway of adenosine nucleotide de novo biosynthesis II, superpathway of adenosine nucleotide de novo biosynthesis I, superpathway of guanosine deoxyribonucleotide de novo biosynthesis II, purine ribonucleoside degradation, adenosine deoxyribonucleotide de novo biosynthesis II, and guanosine nucleotide de novo biosynthesis II) compared to samples with low-to-none SARS-CoV-2 infectivity [34]. The augmented functionality of amino acid and nucleotide biosynthesis and carbohydrate metabolism in the gut microbiome suggests an enhanced production of building blocks for macromolecules and energy extraction in bacterial cells. In the macaque model of COVID-19, a decrease in the concentration of SCFAs and alterations in the concentration of several bile acids and tryptophan metabolites in the feces of infected animals were observed, as revealed by targeted quantitative metabolomics [28]. The decreases in the concentration of SCFAs were in agreement with the microbiome compositional changes in COVID-19 patients, where the abundance of SCFA-producing bacteria was significantly decreased compared with controls [15], [34]. Reductions in SCFA production were also observed in animals during influenza infection, which contributed to further gut microbial dysbiosis and pulmonary pneumococcal superinfection [63]. Moreover, SARS-CoV-2 infection led to an impairment of the fecal bile acid pool, where the primary-to-secondary bile acid ratio (as a function of bile acid transformation by the gut microbiota) was changed and positively correlated with serum levels of chemokines such as C-X-C motif chemokine ligand 13 (CXCL13) [28]. Overall, an increase in the fecal bile acid levels was seen in the macaques infected with SARS-CoV-2, suggesting that the infection leads to accelerated transit and/or impaired bile acid reabsorption in the ileum [28]. Incidentally, the amount of tryptamine, a tryptophan metabolite of microbiota known to accelerate bowel transit [64], [65], was also increased in the feces of SARS-CoV-2-infected macaques [28]. Two end-products (quinolinic and picolinic acids) of tryptophan metabolism by the host indoleamine 2,3-dioxygenase (IDO) pathway was detected in the feces, indicating the presence of intestinal inflammation [64]. Collectively, the functional alterations in the gut microbiome of COVID-19 patients are likely the result of both compositional alterations of the gut microbiome and host immune responses against SARS-CoV-2 infection, which are intertwined with host pathophysiology.

The gut mycobiome in COVID-19

The human gut also harbors a large number of fungi, known as the gut mycobiome. The gut fungi have been demonstrated to be causally implicated in microbiome assembly, ecology, and immune development [66], [67]. Our study in patients with COVID-19 also showed alterations in the gut mycobiome, characterized by enrichment of Candida albicans and highly heterogeneous mycobiome configurations [14]. The abundance of opportunistic fungal pathogens, Candida albicans, Candida auris, and Aspergillus flavus was increased in the feces of COVID-19 patients during the disease course (Figure 1; Table 1) [14]. Fungal pathogens associated with pneumonia and respiratory symptoms, Aspergillus flavus and Aspergillus niger, were detected in the fecal samples from a subset of patients with COVID-19, even after disease resolution [14]. Unstable gut mycobiomes and prolonged dysbiosis persisted in a significant proportion (∼ 30%) of COVID-19 patients [14]. Another study investigated the gut mycobiota in both COVID-19 and H1N1-infected patients and found increased fungal load and enrichment of fungi, including Candida species, in both groups of patients [68]. Presence of Aspergillus niger was positively correlated with diarrhea, while the abundance of Penicillium citrinum was inversely correlated with blood levels of CRP [68]. Aspergillus infections were recently reported in respiratory tract secretions and tracheal aspirates in patients with COVID-19 [69], [70]. Aspergillus is a genus of ubiquitous fungi that cause a variety of pulmonary and respiratory symptoms [71]. Aspergillus may harness the host with immune dysfunction and affect the clinical features and disease course [71]. Cough was found to be more frequent in COVID-19 subjects with Aspergillus infections than those who were not infected [72]. COVID-19 patients who had Aspergillus flavus presence in feces also presented with cough during hospitalization, suggestive of a link of gut mycobiome in the gut–lung axis. These data suggest a gut mycobiome dysbiosis in COVID-19 and its relationship with a systemic dysregulation of host immunity.

Overall, such fungal bloom in the gut of patients with COVID-19 is likely a result of SARS-CoV-2 infection. Secondary fungal infection or co-infection in patients with COVID-19 during the pandemic was frequently observed [73], [74], [75]. Candida and Aspergillus lineages were amongst the specific opportunistic fungal pathogens enriched in patients with COVID-19 during the disease course, particularly Candida albicans [14], [68]. Candida albicans has been shown to impair gut microbiome assembly in both humans and mice, including gut microbiome reassembly after disruption by antibiotics and inflammation [76], [77]. A recent multi-kingdom microbiome study in preterm infants to elucidate the ecological drivers of gut microbiota assembly and dynamics found that between-kingdom interactions have a key role in community dynamics and that the single fungal species, Candida albicans, inhibits multiple dominant genera of gut bacteria (including Klebsiella and Escherichia) [78]. Our prior fecal microbiota transplantation (FMT) study in CDI also demonstrated that presence of Candida albicans in donors or recipients impairs colonization of donor bacteria into recipients, therefore nullifying FMT efficacy in clearing CDI, in both humans and mice [76]. Surprisingly, such inhibitory effect of a single fungus on the restoration of the gut bacterial microbiome can be extended to other fungi, including Aspergillus penicillioides and Penicillium brocae [76]. These studies suggest a crucial role of gut fungi in the gut microbiome ecology, revealing the centrality of simple microbial–microbial interactions in shaping host-associated microbiota. Moreover, gut colonization by Candida albicans can aggravate inflammation in the gut and non-gut tissues [79], [80]. Therefore, the opportunistic expansion of certain fungi in COVID-19 patients potentially has a deleterious role on gut microbiome assembly, where a persistent gut microbiome dysbiosis is consistently seen even after disease resolution and hospital discharge. The long-term effect of gut fungi on the gut microbiome and host health remains to be further investigated.

The gut virome in COVID-19

In addition to the bacteria and fungi, the human gut also harbors an immense diversity of viruses collectively known as the gut virome [81], [82]. Virome consists of both RNA and DNA viruses that chronically infect their eukaryotic (humans, animals, and plants) and prokaryotic hosts (bacteria) [81]. The gut virome serves to modulate the ecology of the co-resident gut bacterial microbiota as well as the immunity of the mammalian host [82].

By shotgun RNA sequencing of the fecal RNA virome, an active presence of SARS-CoV-2 was found in 47% of patients with COVID-19, even in the absence of GI symptoms and after respiratory clearance of SARS-CoV-2 [34]. Meanwhile, COVID-19 patients also had underrepresentation of pepper mild mottle virus (RNA virus), which may originate from diet [13]. By shotgun DNA sequencing of the fecal DNA virome, 19 virus species were identified to be enriched in COVID-19 patients, whereas 26 virus species were enriched in non-COVID-19 controls [13]. Among them, a majority (18 out of 26 virus species) of the DNA viruses enriched in the feces of non-COVID-19 controls were prokaryotic viruses, particularly bacteriophages (16 of 18) (Figure 1; Table 1) [13]. In contrast, more eukaryotic viruses (11 out of 19 virus species) were enriched in the feces of COVID-19 patients [13], which may be a result of SARS-CoV-2 infection. The eukaryotic viruses may harness the immune dysfunction of the host after SARS-CoV-2 infection to expand [81]. The gut virome in COVIID-19 showed more stress-, inflammation-, and virulence-associated gene coding capacities [13]. At patient baseline, the fecal abundance of the RNA virus, pepper chlorotic spot virus, and multiple bacteriophage species was inversely correlated with COVID-19 severity [13]. The abundance of these viruses was also inversely associated with blood levels of pro-inflammatory proteins, white cell count, and neutrophil count [13], indicating that gut resident viruses might tune host immune response to SARS-CoV-2 infection [81], [83]. These data highlight that the gut virome may contribute to immunological and physiological changes in the host during COVID-19. Administration of the antiviral medication lopinavir-ritonavir was associated with the decreased abundance of Listeria phage in COVID-19, suggesting that use of antivirals may tune host bacteriophage–bacteria ecology in the gut, also likely a result of its role in modulating host immune defense against SARS-CoV-2.

Among the COVID-19-enriched viruses, Escherichia phage and Enterobacter phage were prominent [13]. Expansion of these phages has been causally implicated in gut inflammation and host interferon response in mice and humans [41], [84]. In addition, the abundance of their host bacteria was also increased in the gut after SARS-CoV-2 infection [19], [28]. The co-expansion of Escherichia phage and Escherichia was also reported in gut inflammation, and the bloom of Escherichia phage is potentially triggered by lysis of its bacterial host Escherichia under inflammatory conditions [41], [85]. Gut inflammation per se is able to boost bacteriophage transfer between bacteria [86]. Therefore, the alterations in the ecology of gut virome, particularly in the bacteriophage community, are at least partly caused by the alterations of the bacterial microbiome under the influence of SARS-CoV-2 infection and the subsequent immune dysfunction. Similarly, the gut virome dysbiosis persisted along with the dysbiosis of the gut bacterial microbiome, even after disease resolution of COVID-19 [13]. A strong correlation between the composition of virome and bacterial microbiome in COVID-19 patients was observed [33]. Combined ecological network analysis of the virome and bacterial microbiome in COVID-19 revealed that three bacterial species, Faecalibacterium prausnitzii, Bacteroides vulgatus, and Ruminococcus gnavus (the abundance of these bacterial species was also associated with COVID-19 and/or disease severity [12], [15], [19], [28], [34]), and Microviridae bacteriophages constitute central network nodes [33]. These bacterial and viral species may be keystone species that play prominent roles in mediating microbial–microbial interactions in the gut microbial ecology.

Concluding remarks and perspectives

SARS-CoV-2 infection leads to complicated immunologic and pathophysiologic responses in the host. Along with the phenotypic changes in the host, the gut microbiome is broadly altered in COVID-19, including the bacterial microbiome, mycobiome, and virome. Moreover, subsequent blooms of opportunistic bacteria, fungi, and viruses under circumstances of SARS-CoV-2 infection and quiescent/overt gut inflammation in COVID-19 pose further threats to host health and gut microbiome restoration. Such expansions in certain microbial species and decreases in microbiome diversity in conjunction with the impaired host immunity may hinder reassembly of the gut microbiome post COVID-19. Consequently, the altered gut microbiome ecology persists even after disease resolution. Overall, the intricate microbiome ecological network in a steady state is significantly weakened in COVID-19, shifting to one predominated by COVID-19-enriched microbes.

It is well-known that confounding factors such as treatment and diet can significantly affect the gut microbiome composition. However, due to the acute nature of COVID-19, controlling for these confounding factors or including treatment-naïve COVID-19 patients seems infeasible. Therefore, some of the differences between the microbiomes of COVID-19 patients and controls, and of those between disease stages (i.e., mild vs. severe COVID-19 cases), could be attributed to treatment regimens and/or diet. Albeit, we observed consistent microbiome changes across studies, including decreases in the abundance of Eubaterium and SCFA-producing bacteria [12], [15], [19], [33], [34]. In addition, we observed that SARS-CoV-2 infection predominated over medications and diet in affecting the gut virome alterations in patients with COVID-19 [13]. These results together suggest that SARS-CoV-2 infection might be a crucial contributor to the gut microbiome dysbiosis in patients with COVID-19. Although studies have demonstrated that the infection of SARS-CoV-2 would lead to the altered gut microbiome, the causal relationships among the baseline gut microbiome (before infection) that regulates ACE2 expression and host immune status, infectivity/severity of SARS-CoV-2, and altered gut microbiome after infection are complicated. Potential causal loops may also exist, for example, compositions of the gut microbiome (baseline) favor the infection of SARS-CoV-2, and subsequent infection of SARS-CoV-2 induces the change of gut microbial ecology. Little is known about the relative contribution of the baseline status of the gut microbiome to the later-on infection and the dynamics of the altered gut microbiome. Beyond that, it is also paramount to further understand how the gut microbiome regulates host immunity against SARS-CoV-2 infection, therefore disease severity, as well as the long-term impact of COVID-19 on the gut microbiome reassembly in relation to host health after the pandemic.

CRediT author statement

Tao Zuo: Conceptualization, Methodology, Writing - original draft. Xiaojian Wu: Writing - review & editing, Supervision. Weiping Wen: Writing - review & editing, Supervision. Ping Lan: Conceptualization, Supervision, Project administration. All authors have read and approved the final manuscript.

Competing interests

The authors have declared no conflict of interest.