Next Article in Journal
Rethinking Keratoplasty for Patients with Acanthamoeba Keratitis: Early “Low Load Keratoplasty” in Contrast to Late Optical and Therapeutic Keratoplasty
Previous Article in Journal
Comparison of Power Ultrasound and NALC-NaOH Decontamination Methods for Stool Mycobacterial Culture: A Prospective Study
Previous Article in Special Issue
Ethanol-Producing Enterocloster bolteae Is Enriched in Chronic Hepatitis B-Associated Gut Dysbiosis: A Case–Control Culturomics Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Cultured Gut Microbiota Using MALDI-TOF MS in COVID-19 Patients from Serbia during the Predominance of the SARS-CoV-2 Omicron Variant

by
Aleksandra Patić
1,2,
Gordana Kovačević
1,3,*,
Vladimir Vuković
1,4,
Ivana Hrnjaković Cvjetković
1,2,
Mioljub Ristić
1,4,
Biljana Milosavljević
1,
Deana Medić
1,2,
Milan Djilas
1,
Jelena Radovanov
1,3,
Aleksandra Kovačević
5,
Tatjana Pustahija
1,4,
Dragana Balać
1,6 and
Vladimir Petrović
1,4
1
Institute of Public Health of Vojvodina, 21000 Novi Sad, Serbia
2
Department of Microbiology with Parasitology and Immunology, Faculty of Medicine, University of Novi Sad, 21000 Novi Sad, Serbia
3
Faculty of Medicine, University of Novi Sad, 21000 Novi Sad, Serbia
4
Department of Epidemiology, Faculty of Medicine, University of Novi Sad, 21000 Novi Sad, Serbia
5
Clinic of Nephrology and Clinical Immunology, University Clinical Center of Vojvodina, 21000 Novi Sad, Serbia
6
Department of Hygiene, Faculty of Medicine, University of Novi Sad, 21000 Novi Sad, Serbia
*
Author to whom correspondence should be addressed.
Microorganisms 2024, 12(9), 1800; https://doi.org/10.3390/microorganisms12091800
Submission received: 20 July 2024 / Revised: 8 August 2024 / Accepted: 13 August 2024 / Published: 30 August 2024
(This article belongs to the Special Issue The Role of the Gut Microbiome in Bacteria and Virus Infections)

Abstract

:
The currently dominant SARS-CoV-2 omicron variant, while causing mild respiratory symptoms, exhibits high transmissibility, drug resistance, and immune evasion. We investigated whether the presence of the SARS-CoV-2 affected the dynamics of fecal microbial composition isolated in culture in moderate COVID-19 patients. Blood, stool, and medical records were collected from 50 patients with confirmed SARS-CoV-2 infection. Two samples were taken per patient, at disease onset (within 5 days) and after symptom resolution (30–35 days). The part of the gut microbiota identifiable using MALDI-TOF MS was analyzed, and inflammatory cytokines and blood markers were measured in serum. The analysis identified 566 isolates at the species level, including 83 bacterial and 9 fungal species. Our findings indicate a change in the gut microbiota composition isolated in culture during the initial phase of infection, characterized by the proliferation of opportunistic bacteria such as Enterococcus spp. and Citrobacter spp., at the expense of beneficial commensal bacteria from the genus Bacillus and Lactobacillus. Additionally, the enrichment of fungal pathogens in fecal samples collected 30 days after the cessation of disease symptoms might suggest a prolonged disruption of the gut microbiota even after the resolution of COVID-19 symptoms. This study contributes to a growing body of evidence on the systemic effects of SARS-CoV-2 and highlights the importance of considering gastrointestinal involvement in the management and treatment of COVID-19.

1. Introduction

The term gut microbiota refers to the collection of resident microorganisms that inhabit the gastrointestinal tract (GIT), including bacteria, fungi, viruses, phages, and archaea [1,2]. Various factors, such as nutrition, genetics, environment, and lifestyle, influence the composition of gut microbiota [3,4]. Bacteria are essential components of this community due to their numerous functions, including food fermentation, protection against pathogens, immune response stimulation, and vitamin production [5,6]. Generally, bacteria inhabiting the gastrointestinal tract are categorized into six phyla: Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria, Fusobacteria, and Verrucomicrobia, with Firmicutes and Bacteroidetes being the dominant taxa [4]. The gut microbiome plays an important role in immune system regulation, and alterations in its composition have been observed in various infectious diseases [7,8]. Different bacterial and viral infections, including SARS-CoV-2, have previously been associated with persistent gastrointestinal symptoms and post-infectious irritable bowel syndrome [9]. A healthy gut microbiota helps protect against the colonization of harmful pathogens by competing for resources and producing antimicrobial substances [10]. This competitive exclusion helps maintain a balanced and protective environment in the gut. Understanding the intricate relationship between the gut microbiota and human health is an active area of research, and ongoing studies continue to uncover new insights into the role of these microorganisms in maintaining overall human health [11,12].
The pathogenesis of coronavirus disease (COVID-19) is a consequence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus’s tropism, which is determined by the distribution of the angiotensin-converting enzyme 2 (ACE2) receptor. In humans, the ACE2 receptor is present in the nose, lungs, kidneys, heart, and gastrointestinal tract. This distribution enables SARS-CoV-2 to infect not only human respiratory epithelial cells but also the cells of the gastrointestinal (GI) system [13,14]. In addition, several studies conducted in different regions have suggested that the gut microbiota may contribute to the pathogenesis of COVID-19 and the disease outcomes [2,15]. It has also been confirmed that SARS-CoV-2 infection is associated with an altered composition of the intestinal microbiota and correlated with inflammatory and immune responses [10,13]. However, the largest number of clinical and preclinical studies on the impact of SARS-CoV-2 on the gut microbiota have focused on earlier SARS-CoV-2 variants that caused severe forms of the disease [16]. Since the start of the COVID-19 pandemic, the World Health Organization (WHO) has designated several strains as variants of concern (VOCs). These include strains B.1.1.7 (alpha), B.1.351 (beta), P.1 (gamma), B.1.617.2 (delta), and B.1.1.529 (omicron) [17]. Compared to previous variants, omicron is characterized by increased transmissibility and milder forms of the disease, accompanied by a reduced rate of hospitalization and a shorter recovery time [18]. The virus has continued to evolve and accumulate mutations, particularly in the receptor-binding domain of the spike protein. This has led to the simultaneous rise of multiple omicron descendants, sharing common mutations that enhance the virus’s ability to evade neutralizing antibodies. Omicron subvariants currently circulating globally include JN.1, KP.2, KP.3, and KP.1 [19].
Given the global predominance of the omicron variant and its subvariants and the lack of research on the impact of SARS-CoV-2 on intestinal microflora in Serbia and the region, our study aims to explore the dynamics of microbial composition isolated in culture in SARS-CoV-2-positive patients with moderate clinical symptoms during their active disease and in the convalescence period. We seek to analyze the microbial composition isolated in culture at the onset of illness and one month after recovery and to determine whether any potential disruptions in the gut microbiota composition resolve after SARS-CoV-2 viral clearance and whether they are associated with disease severity in patients with COVID-19. Additionally, we investigated whether alterations in the gut microbiota composition and diversity correlate with biomarkers of inflammatory response in COVID-19 patients.

2. Materials and Methods

2.1. Sample Collection

This prospective cohort study included 50 adult patients with a laboratory-confirmed diagnosis of COVID-19, defined as a positive SARS-CoV-2 result from nose/throat swabs validated by either the RT-PCR method or a rapid antigen test. The research was conducted at the Institute of Public Health of Vojvodina, Novi Sad, Serbia, from July 2022 to August 2023. Within five days after a confirmed positive result, initial stool and blood samples were collected from study participants, followed by additional stool and blood samples 30–35 days after the cessation of disease symptoms. At the initial sampling, participants provided data on their socio-demographic characteristics, current illness symptoms, history of COVID-19, COVID-19 vaccination status, and comorbidities. The progression or regression of disease symptoms, along with details on therapy and treatment, was monitored every three days until symptom resolution.
The study included only individuals who had not used antibiotics and probiotics for at least a month prior to enrollment, were not undergoing chemotherapy or radiation therapy, were not immunocompromised (e.g., HIV/AIDS, transplant, cancer, etc.), and did not have a history of inflammatory bowel disease. Pregnant women and children under 18 years of age were not included in the research. An overview of the study procedure and general sampling process is presented in Figure 1.
Informed written consent was obtained from each participant. The study protocol was approved by the Ethics Committee of the Institute for Public Health of Vojvodina (Decision number: 01-991/1-1 of 27 June 2022).

2.2. Detection of SARS-CoV-2 RNA in Stool Samples

Study participants received detailed instructions for the proper collection and transport of fecal samples. Immediately upon arrival at the laboratory, fecal samples were processed to create a 10% suspension by homogenizing 1 g of fecal material in 9 mL of sterile balanced salt solution (0.9% NaCl) while the remainder was sent for microbiological analysis. Approximately 150 µL of this prepared suspension served as the starting material for the extraction of viral nucleic acid. If fecal suspensions could not be immediately subjected to nucleic acid isolation, they were stored at −80 °C. Viral RNA was extracted using the Ribovirus extraction kit (Sacace, Biotechnologies, Como, Italy) following the manufacturer’s instructions. The isolated RNA was stored at −80 °C until amplification. Detection of SARS-CoV-2 RNA was performed using the commercially available GeneProof SARS-CoV-2 PCR kit (GeneProof a.s., Brno, Czech Republic).

2.3. Bacterial Identification and Characterization

Identification and characterization of bacterial species within fecal samples were performed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) technology, a soft ionization technique achieved by mixing a matrix with bacterial colonies on metal plates. Mass spectra were acquired using a Microflex BioTyper spectrometer (Bruker Daltonics, Billerica, MA, USA) equipped with a nitrogen laser and controlled by Flex Control software ver. 3.4, Build 135.10 (Bruker Daltonics). Spectra were generated by measuring the time of flight of ions through the tube to the detector. Spectral comparison with the data from the available defined database was automatically performed and permitted identification.
A total of 100 stool samples (comprising both the first and second samples from COVID-19 patients) were prepared for analysis. Stool specimens were inoculated onto several nutrient media, including Columbia Blood agar, Endo agar, Salmonella Shigella (SS) agar, Schaedler agar, and Sabouraud agar. Following the culture method of seeding samples on nutrient media and subsequent incubation, the most numerous isolated colonies from each nutrient medium were selected and inoculated onto new plates to obtain pure cultures. For the direct colony method, bacteria were applied as thin films to 96-spot, polished, stainless steel target plates using sterile toothpicks (Bruker Daltonik GmbH, Leipzig, Germany). The bacteria were then left to dry at room temperature for 1 min. Subsequently, 1 μL of the matrix solution, containing a saturated α-cyano-4-hydrocinnamic acid (Bruker Daltonik) in 50% acetonitrile (Sigma-Aldritch, Burlington, MA, USA) and 2.5% trifluoroacetic acid (Sigma-Aldritch), was applied to the samples and co-crystallized with them at room temperature for 10 min. Spectra in the mass range of 2 to 20 kDa were collected using the Auto Execute option by accumulating 240 laser shots acquired at 30–40% of maximum laser power.

2.4. Detection of SARS-CoV-2 IgG Antibodies and Inflammatory Biomarkers

All serum samples collected from project participants were tested to determine the presence of anti-SARS-CoV-2 IgG antibodies. The analysis was conducted on serum samples with volumes ranging from 0.5 to 1 mL utilizing the “LIAISON® SARS-CoV-2 TrimericS IgG” test developed by LIASON DiaSorin (DiaSorin SpA, Saluggia, VC, Italy). The tests were executed in a fully automated manner on the LIAISON® XL (DiaSorin SpA, Saluggia, VC, Italy) Analyzer. Values greater than or equal to 13.0 AU/mL were interpreted as positive.
An immunoturbidimetric method was used to determine human C-reactive protein (CRP) levels in serum. Values between 1 and 5 mg/L were considered normal. Levels of pro-inflammatory and anti-inflammatory cytokines were quantified using commercially available enzyme-linked immunosorbent assay (ELISA) kits. Specifically, immunoassays for human interleukin 1 alpha (IL-1α), human interleukin 6 (IL-6), and human interleukin (IL-10), manufactured by Elabscience Bionovation Inc. (Houston, TX, USA), were employed. Absorbance readings at 450 nm were performed using an automated ELISA reader for microtiter plates. The concentrations of the target analytes were determined based on the standard curve and expressed in pg/mL. Using CurveExpert Professional 2.7.3 software, the target concentrations of the samples were interpolated from the standard curve.
Since the manufacturer, Elabscience, did not specify the reference range for pro- and anti-inflammatory cytokines, reference values were taken from the relevant scientific literature. The reference range in healthy individuals is approximately 5.186 pg/mL for IL-6 [20]; 4.8–9.8 pg/mL for IL-10 [21]; and 0–5 pg/mL for IL-1α [22].

2.5. Statistical Analysis

All collected data were inserted into a specially created database and used in the analyses. Results were presented as means ± standard deviation (SD). The chi-square test was applied to compare the frequencies of different bacterial taxa. For data following a normal distribution, paired t-tests were utilized. Categorical data were analyzed using the McNemar test. Outcomes were compared, including changes in the presence of specific bacterial taxa at the onset of infection and after the resolution of symptoms. The Pearson correlation coefficient (r) was computed to assess the strength and direction of correlations. Statistical analyses and data visualization were performed using GraphPad Prism (8.0.1 software, GraphPad software Corporation, San Diego, CA, USA), Microsoft Power BI, and Microsoft Excel 2016 (Microsoft Office, Redmond, WA, USA, Professional Plus 2016). A p-value of <0.05 was considered statistically significant.

3. Results

3.1. Analysis of SARS-CoV-2 RNA in Stool Samples

A cohort of 50 patients with SARS-CoV-2 infection included 22 men and 28 women. The mean age of the cohort was 50.5 years. The clinical presentation varied, with the majority (64%, 32/50) presenting with a moderate illness, followed by 24% (12/50) with a moderately severe profile. A severe illness was noted in 8% (4/50) of the cases, while 4% (2/50) were asymptomatic. Within five days of symptom onset, a total of 66% (33/50) of patients excreted SARS-CoV-2 RNA in their first stool sample. Notably, only one patient (2%) continued to shed viral RNA 30–35 days after the cessation of disease symptoms. About 42 (84%) of the study participants had received at least two doses of the COVID-19 vaccine. Approximately 47 (94%) of the study participants had SARS-CoV-2 IgG antibodies at the beginning of the infection (first blood sample), while the second measurement after 30–35 days showed that all participants had SARS-CoV-2 IgG antibodies. Half of the study participants reported gastrointestinal (GI) complaints in addition to other symptoms. Table 1 compares detailed characteristics of individuals with and without SARS-CoV-2 RNA fecal shedding in their first stool sample, including data on their sex, age group, COVID-19 severity, duration of symptoms, gastrointestinal symptoms, prior infection status, and vaccination status. As demonstrated in Table 1, the presence of SARS-CoV-2 RNA was not associated with the explored demographic and clinical characteristics of the patients.

3.2. Identification and Analysis of Bacterial Taxa in Fecal Samples from COVID-19 Patients

A total of 100 fecal samples were analyzed using MALDI-TOF mass spectrometry. Of these, 50 were collected during the initial phase of COVID-19 (within 5 days of a laboratory-confirmed positive COVID-19 test result) and 50 were collected 30 to 35 days after the cessation of disease symptoms. Extensive analysis identified 566 isolates at the species level (Table S1 in the Supplementary Materials). Using various cultivation techniques, 83 bacterial species and nine fungal species were identified. In fecal samples collected during the initial phase of the disease, a total of 272 isolates were identified. Subsequent sampling after 30–35 days revealed an increase in species diversity, with 294 isolates identified (Figure 2A).
Among the 83 identified bacterial species, 13 (15.7%) were observed only once, 45 (54.2%) were detected in two to five different samples, 14 species (16.9%) were found in six to ten different samples, and 11 species (13.2%) were identified in more than ten distinct samples. At the phylum level, Firmicutes, Proteobacteria, and Actinobacteria were the most dominant, while Bacteroidetes were identified in a limited number of isolates. These phyla were present in both sample groups with no significant changes in relative abundance based on the χ2-test (Figure 2B). Enterococcus faecium was the most commonly identified bacterial species, with 71 isolates, followed by Escherichia coli, with 47 isolates, Hafnia alvei, with 40 isolates, Citrobacter braakii, with 18 isolates, and both Citrobacter freundii and Enterococcus faecalis, with 17 isolates each. The analysis of bacterial genera showed an abundance of Enterococcus, which was present in 21.3% of the first fecal samples and 15.6% of the second samples. Other genera identified in notable quantities in the first samples included Bacillus (9.1%), Escherichia (8.4%), Citrobacter (7.7%), Hafnia (7.3%), and Pseudomonas (6.9%). Additionally, in the first samples, fungi were found in 2.9% of isolates. In the second samples, collected 30 to 35 days after the cessation of disease symptoms, the genera Bacillus (11.3%), Pseudomonas (8.1%), Escherichia (8.1%), Hafnia (6.8%), and Citrobacter (6.1%), as well as fungi (6.1%), were found. These findings were visually represented using percentages in a stacked bar plot (Figure 2C).
To assess the diversity of gut community species at the onset of illness and 30–35 days after the cessation of disease symptoms, Shannon’s diversity index (Figure 2D) and Simpson’s index (Figure 2E) were applied. Although both indices demonstrated species abundance in the community, the observed differences were not statistically significant based on the t-test (Simpson, p = 0.906; Shannon, p = 0.620).

3.3. Analysis of Commensal and Opportunistic Bacterial Taxa Components

Considering pathogenic potential, host interaction, environmental conditions, host status, and genetic and phenotypic characteristics, we classified the identified bacterial species and genera into opportunistic and commensal categories. We examined the differences in their presence between the initial fecal samples collected within 5 days of a laboratory-confirmed positive COVID-19 test and those collected 30–35 days after the cessation of disease symptoms. Our results indicated that opportunistic taxa (including genera such as Enterococcus, Pseudomonas, Klebsiella, Citrobacter, Enterobacter, Proteus Moellerella, Yersinia and Serratia), along with fungal pathogens, were more prevalent than commensal species at the onset of the disease (60.3% vs. 39.7%). However, in samples collected 30–35 days after the cessation of disease symptoms, this difference was no longer observed (49.2% vs. 50.8%).
Further, we focused on the presence of bacterial genera to assess whether they changed significantly over time using McNemar’s test with a significance level of 0.05 (Table 2). Since we performed multiple comparisons (10 tests), Bonferroni correction was applied to adjust the significance threshold. The adjusted significance level was α′ = 0.005. Enterococcus faecium emerged as the most frequently detected species among COVID-19 patients, along with the genus Enterococcus. Despite this, McNemar’s test indicated no significant change in the number of Enterococcus detections from the beginning of the disease to 30–35 days after the cessation of disease symptoms, likely due to the high number of patients with persistent Enterococcus presence. Moreover, the presence of opportunistic Proteobacteria belonging to the order Eubacteriales (including genera such as Klebsiella, Citrobacter, Enterobacter, Proteus, Moellerella, Yersinia, and Serratia) decreased during the follow-up period. However, this change was not statistically significant after the Bonferroni correction. Additionally, while the initial analysis suggested a higher presence of fungi in samples collected 30–35 days after the cessation of disease symptoms, this finding did not remain statistically significant after the correction for multiple comparisons.
The total number of commensal species detected at the onset of the disease (39.7%) increased in samples collected 30–35 days after the cessation of disease symptoms (50.8%). However, there were notable shifts in their distribution. Members of the order Bacillales (including the genera Bacillus, Sporosarcina, Lysinibacillus, Psychrobacillus, Listeria, Paenibacillus, and Peribacillus) showed a depletion. Conversely, species from the order Lactobacillales (including Lactobacillus, Carnobacterium, Filifactor Limosilactobacillus, Streptococcus, Ligilactobacillus, and Tissierella) showed an enrichment in the second sample compared to the first. Initially, some p-values related to changes in the presence of species from Bacillales and Lactobacillales were below the traditional threshold of 0.05, suggesting potential significance. However, these did not meet the stricter criteria established by the Bonferroni correction for multiple testing, indicating that the observed differences were not statistically significant. Likewise, several commensal genera that were rarely identified or present in small numbers (such as Corynebacterium, Paenarthrobacter, Microbacterium, Micrococcus, Brevibacterium, Nocardiopsis, and Rhodococcus) exhibited increased presence after the cessation of disease symptoms. Despite these trends, none of these changes were statistically significant after the correction for multiple comparisons. In contrast, the number of bacteria from the genera Hafnia, Escherichia, and Raoultella remained relatively constant throughout the disease course.

3.4. Relationship between the Detection Rates of Opportunistic Bacteria and Demographic and Clinical Factors

Table 3 presents the detection rates of opportunistic bacteria in COVID-19 patients at two different time points: at the onset of the disease (stool sample 1) and 30 days after the cessation of disease symptoms (stool sample 2). The data are stratified by various demographic and clinical characteristics. A statistically significant difference was found in the percentage of opportunistic bacteria between paired samples, with the first stool samples showing a higher mean percentage (mean = 60.26%, SD = 27.10%) compared to the second stool samples (mean = 49.20%, SD = 27.94%) in patients with COVID-19 (p = 0.048).
A statistically significant difference was also observed between paired samples from individuals aged 50 to 59 years (p = 0.0357), as well as those who did not excrete SARS-CoV-2 through stool (p = 0.0332). Although males showed a decrease from a mean percentage of 49.68% to 36.13%, the difference found was not statistically significant (p = 0.062). The reduction in opportunistic taxa was not associated with comorbidities, body mass index (BMI), or symptom severity.
Regarding the duration of disease symptoms, two (4%) patients reported a duration of 1–5 days, 10 (20%) reported 6–9 days, 32 (64%) reported 10–14 days, and six (12%) reported more than 14 days. Further examination revealed no significant correlation between the percentage of opportunistic bacteria present in the samples and the duration of COVID-19 symptoms (Figure 3). The correlation coefficients were close to zero (r = −0.177 for sample I and r = −0.025 for sample II), and the p-values (p = 0.219 for sample I and p = 0.862 for sample II) indicated that any observed trends were likely due to random variation rather than a real effect.

3.5. Analysis of Inflammatory Biomarkers and SARS-CoV-2 IgG Antibody Levels in COVID-19 Patients

Our study further investigated serum levels of CRP, IL-1α, IL-6, and IL-10 as predictive markers for identifying patients at risk of worsening COVID-19. Figure 4 illustrates the levels of various cytokines and inflammatory markers in patients with COVID-19 at two different time points: within 5 days of the first positive test result (≤5 days) and 30–35 days after the cessation of disease symptoms (30–35 days).
CRP levels exhibited significantly higher values within 5 days of the first positive test compared to 30–35 days after the cessation of disease symptoms, indicating an elevated inflammatory response during the acute phase (Wilcoxon signed-rank test, p < 0.001). The variability in CRP levels was also more pronounced in this period, reflecting the dynamic inflammatory status early in the disease course. Additionally, the results indicate a significant increase in IgG antibody levels over time, with p-values less than 0.001 at both time points. In contrast, the IL-6, IL-1α, and IL-10 levels showed no statistically significant differences between the two time points. The median value for IL-6 was slightly above 5 pg/mL, with most values clustering around that level. The data spread showed some variability but remained relatively concentrated. The IL-1α levels also exhibited slightly higher variability within 5 days of the first positive test compared to 30–35 days after the cessation of disease symptoms, showing a higher peak at the earlier time point. Additionally, the IL-10 levels showed a similar distribution at both time points, with slightly higher levels within 5 days of the first positive test.

4. Discussion

The present study included a cohort of 50 SARS-CoV-2 positive patients receiving outpatient care due to their stable health condition. The research sample was collected between July 2022 and May 2023, a period coinciding with the dominance of the SARS-CoV-2 omicron variant [23]. In Serbia, as well as globally, the omicron variant has led to a significant number of cases, although the incidence of severe illness and mortality has remained relatively low compared to previous viral strains.
The clinical manifestations among the study population were predominantly moderate, characterized by symptoms such as fever, sneezing and cough. Previous investigations into different SARS-CoV-2 variants of concern (VOCs) during earlier waves found that 40 to 85% of infected individuals excreted viral RNA in their feces [24]. In our study, approximately two-thirds (66%) of the omicron SARS-CoV-2 positive patients exhibited viral RNA shedding in their stool at the time of diagnosis, while only one patient (2%) continued to shed the virus 30–35 days after the cessation of disease symptoms. These findings confirm a well-documented observation that, despite omicron causing milder respiratory symptoms [25], the mechanisms underlying fecal virus excretion are likely similar to those of previous VOCs.
Although our study cohort encompassed individuals with mild to moderately severe symptoms, our findings elucidate that during the initial five days of omicron-associated COVID-19, there was a proliferation of opportunistic bacteria at the expense of commensal species, consistent with observations from previous waves [15,16,26]. It is noteworthy that geographical, socio-economic, and cultural factors play significant roles in modulating the composition of the microbial community [27]. Therefore, the predominant bacterial species observed in our study may exhibit certain variations compared to investigations conducted in other regions. Based on the Shannon index, we cannot conclude that there was a decrease in diversity in the fecal samples collected from patients at the onset of the disease and 30–35 days after the cessation of disease symptoms. However, it is evident that the relative abundance of certain taxa changed significantly without affecting overall diversity. Shannon’s and Simpson’s indices primarily assess the overall species richness and evenness within a community, which remained stable despite fluctuations in the relative abundance of individual taxa. The significant compositional changes observed indicate a dynamic response of the gut microbiota to the Omicron infection, likely involving specific taxa rather than a broad alteration in community diversity. Therefore, while the overall diversity metrics did not indicate significant differences, the underlying microbial community composition experienced notable alterations during the course of the disease. Among patients from our region, we observed a reduction in beneficial commensal bacteria belonging to the orders Bacillales (including genera such as Bacillus, Lisinibacillus, and Psychrobacillus) and Lactobacillales (including genera such as Lactobacillus, Leuconostoc, Carnobacterium, and Limosilactobacillus). Specifically, there was a substantial enrichment of the gut microbiota with the genus Enterococcus (E. faecalis, E. faecium) and a significant presence of Proteobacteria belonging to the order Enterobacterales, such as the Citrobacter genus (C. braakii and C. freundii).
Enrichment of intestinal microbiota with the species Enterococcus faecalis could be highlighted as a significant characteristic of our studied cohort. While we cannot assert that high proportions of enterococci are exclusive to COVID-19, it is worth noting that other studies have also reported an unexpectedly high frequency of enterococci in COVID-19-positive patients [28,29,30]. For instance, in some Italian patients with moderate/severe pneumonia, the gut microbiota was almost mono-dominated by Enterococcus spp., mostly E. faecium, E. hirae, E. faecalis, and E. villorum [31]. Additionally, our findings regarding the significant presence of Citrobacter freundii align with previous studies reporting an increase in Proteobacteria, particularly Citrobacter spp., in COVID-19 patients. This enrichment of Citrobacter freundii may contribute to the dysbiosis observed in COVID-19 patients, potentially impacting the immune response and disease progression [32]. Several mechanisms may underlie the increased prevalence of enterococci in COVID-19 patients. First, it is important to recognize that Enterococcus spp. are common inhabitants of the human gastrointestinal tract, capable of utilizing a wide range of substrates for growth and adapting to diverse environmental conditions [33]. The SARS-CoV-2 infection can lead to significant alterations in gut metabolism, including changes in nutrient availability and gut pH. Their exceptional adaptability may allow enterococci to gain a competitive advantage over other bacterial taxa under these altered metabolic conditions [34,35]. Moreover, SARS-CoV-2 infection is associated with considerable immune dysregulation, including altered cytokine profiles and impaired local immune responses, which can suppress the growth of competing bacteria and reduce mucosal immunity. This immune imbalance may facilitate the proliferation of Enterococcus spp. [36]. Therefore, it is reasonable to hypothesize that SARS-CoV-2-induced intestinal dysbiosis creates a niche that favors the proliferation of enterococci [31,36,37]. The synergy between SARS-CoV-2 and enterococci may have clinical implications. SARS-CoV-2 can infect enterocytes in the gut, potentially compromising the gut barrier and allowing the translocation of gut bacteria into the bloodstream [31,36,37]. Similarly, the presence of Citrobacter freundii, another opportunistic pathogen, could exacerbate gut dysbiosis and contribute to systemic inflammation, further complicating the COVID-19 management [32].
Additionally, the genus Enterococcus includes some of the most virulent bacterial strains implicated in human infections, such as Enterococcus faecium (known for its resistance to vancomycin) and Enterococcus faecalis (which shows resistance to linezolid) [36]. This resistance can pose a challenge during COVID-19 management, especially if secondary bacterial infections arise [31]. Furthermore, Enterococcus spp. possess various virulence factors that contribute to their pathogenicity, including the ability to form biofilms, produce toxins [33], and resist immune responses. These factors can complicate COVID-19, particularly in hospitalized patients or those with weakened immune systems.
Several studies, including our own, have identified a specific enrichment of opportunistic fungal pathogens following the resolution of COVID-19 symptoms [38]. Our findings indicated that fecal samples collected 30–35 days after the cessation of disease symptoms showed an increased abundance of the genera Candida, Geotrichum, Pichia, and Mucor. In contrast, other studies have primarily reported the prevalence of Mucor spp. and Pneumocystis jirovecii. According to one report, abnormalities in the gut mycobiota composition persisted after recovery in patients with mild and moderate COVID-19 symptoms [39]. There is increasing evidence that one of the major risk factors for invasive fungal infections is the use of corticosteroids. We hypothesize that the use of antibiotics in some participants could have contributed to creating favorable conditions for fungal proliferation.
Our study was conducted during a period when empirical antibiotic treatment for COVID-19 was common, although it is no longer recommended [40]. Despite the initial exclusion criteria, some participants were administered antibiotics after diagnosis and sample collection as part of their routine clinical care, which was beyond our control. It should be noted that the data on the impact of antibiotics on the microbiota of COVID-19 patients are conflicting. Several studies have reported that patients treated with antibiotics showed greater depletion of beneficial bacteria than patients who did not receive antibiotics [13,41]. Conversely, Vestad et al. (2022) found similar microbiota in antibiotic-treated and untreated COVID-19 patients [41]. Comparing our results with those of other studies is challenging due to variability in clinical conditions, disease duration, use of therapeutic drugs, treatment protocols, and general lifestyle and nutrition. In our study, the overall composition of the intestinal microbiota isolated in culture was more diverse 30–35 days after the cessation of disease symptoms compared to the initial phase of the disease, with a slightly higher proportion of commensal genera such as Bacillus and Lactobacillus. However, fungal pathogens persisted even after SARS-CoV-2 was cleared from the stool. The observed enrichment of opportunistic bacteria, such as Enterococcus, and the persistence of fungal pathogens post-recovery highlight the need for ongoing monitoring of gut microbiota in COVID-19 patients. This aligns with existing scientific evidence suggesting that gut dysbiosis can exacerbate COVID-19 symptoms, emphasizing the potential role of probiotics and dietary interventions as adjunct treatments [42]. Several studies support the beneficial effects of nutritional therapy, including the use of probiotics, prebiotics, and synbiotics combined with a healthy diet, on human health. These interventions may help mitigate SARS-CoV-2-related gut dysbiosis, as evidenced by recent findings [43].
Gut microbiota can both modulate and be modulated by cytokines produced in the gut, thereby directly and indirectly influencing immune responses [44]. Notably, the levels of IL-6 and IL-10 did not show statistically significant changes among the patients in our study despite an imbalance in the composition of intestinal microbiota favoring opportunistic bacteria at the onset of the disease. However, the observed variations in CRP and IL-1α levels during the early phase of infection highlight the dynamic changes in inflammatory responses that occurred during the course of the disease. CRP is a highly sensitive but nonspecific biomarker that indicates inflammation, tissue damage, and infection, and numerous studies have reported elevated CRP levels in patients with COVID-19 [16]. Although CRP is generally well correlated with IL-6, our study did not observe significant changes in IL-6 levels. The mean calculated IL-6 concentration was approximately 5 pg/mL, which is consistent with values reported by other researchers for patients with mild to moderate COVID-19 symptoms associated with the omicron variant [45]. SARS-CoV-2 infection of enterocytes may lead to subtle disturbances in the intestinal environment, potentially explaining the localized intestinal response observed in our study [46]. Additionally, factors such as effective immune regulation preventing systemic inflammation, antibiotic use, and individual predispositions likely contribute to maintaining a mild systemic immune response while allowing minor changes in gut microbiota. Moreover, our findings showed that approximately 94% of participants had detectable SARS-CoV-2 IgG antibodies at the beginning of the infection with varying levels of antibody titers. All participants showed a presence of IgG antibodies 30–35 days after symptom resolution. This widespread seroconversion indicates a sufficiently strong immune response, which may have contributed to the observed normalization of gut microbiota over time.
The results of our study highlight several important trends in the detection rates of opportunistic taxa in patients with COVID-19 over time. Overall, a reduction in opportunistic taxa was observed from the initial stage of the disease to 30–35 days after the cessation of disease symptoms across the entire population. This reduction might suggest a potential recovery or normalization of the gut microbiota as the symptoms of the disease subside. Furthermore, a decrease in opportunistic taxa over time was noted in patients taking vitamin D supplements. The significant impact of vitamin D supplementation could indicate its role in restoring a healthy microbiota. Recent studies have confirmed the anti-biofilm activity of vitamin D against tested Gram-negative strains [47]. However, this reduction in the number of opportunistic pathogens did not reach statistical significance in individuals with comorbidities. Chronic diseases often involve a deterioration of one or more physiological functions, which can modulate the course of COVID-19 [48].
This study has several limitations. We utilized MALDI-TOF MS for bacterial species identification, a method that relies on aligning the mass spectrum of the measured isolate with an established database. The absence of a comparative analysis between MALDI-TOF MS results and 16S rRNA gene sequencing represents a notable limitation. It should be noted that the 16S rRNA sequencing method, commonly used in gut microbiota studies, may not be sensitive enough to detect small microbiota alterations [49]. Therefore, future studies should be based on newer techniques such as shotgun metagenome sequencing. Anaerobic bacteria are more difficult to culture, so bacterial species that were identified by using MALDI-TOF MS represent only a part of the community composition of the gut microbiota. Furthermore, not all laboratories have sequencing capabilities due to the relatively high cost of the technology and the need for specialized bioinformatic tools for analysis. Therefore, in the absence of sequencing resources, culture and MALDI-TOF MS identification can provide valuable insight into the gut microbiota. We believe that this limitation did not substantially impact our findings due to MALDI-TOF MS’s demonstrated ability to accurately identify closely related species. Secondly, this study did not include a control group of uninfected individuals due to a high variability in their baseline characteristics, like lifestyle factors and health conditions (i.e., nutrition, comorbidities, immune responses, medications, etc.), which all might have an impact on the composition of microbiota isolated in culture in both health and disease, as demonstrated in previous research [50]. Instead, we used a longitudinal approach collecting samples from all participants at the time of disease onset and during convalescence. Third, a relatively small number of participants may limit the generalizability of our findings to the broader population from which the sample was taken. Future studies with larger cohorts are needed to validate these results. Furthermore, in future research, a combination of traditional diversity indices with principal component analysis (PCA) could reveal new patterns and trends in the microbial community structure, thereby providing a more comprehensive understanding of the microbiota dynamics.

5. Conclusions

In conclusion, this study provides valuable insights into the gut microbiota alterations, inflammatory responses, and fecal shedding patterns in patients infected with the omicron SARS-CoV-2 variant. To the best of our knowledge, this is the first study of its kind conducted in Serbia and the region. Despite the milder respiratory symptoms associated with omicron, we found that approximately two-thirds of patients exhibited viral RNA shedding in their stool at the time of diagnosis. This indicates that the mechanisms of fecal virus excretion remain consistent across different SARS-CoV-2 variants of concern (VOCs).
Our findings contribute to the understanding of the gut microbiota in COVID-19 patients by focusing on the identifiable species within the gut, as determined by culture and MALDI-TOF MS identification. It is crucial to highlight the importance of integrating different methodologies to achieve a comprehensive understanding of gut microbial communities. Access to cultured human microbiota offers detailed functional characterization of bacteria and facilitates the discovery of their biological activities during host–bacterial and inter-bacterial interactions in both health and disease.
Additionally, our findings highlight a significant shift in the gut microbiota composition isolated in culture, characterized by a proliferation of opportunistic bacteria such as Enterococcus spp. and Citrobacter spp. at the expense of beneficial commensal bacteria from the genera Bacillus and Lactobacillus. Moreover, the enrichment of fungal pathogens after the cessation of disease symptoms might suggest a prolonged microbiota disruption. The reduction in opportunistic bacterial taxa and stabilization of commensal species over time indicates a potential restoration of gut microbial balance as patients recover from the acute phase of the disease. A significant decrease in CRP levels further supports this recovery trend.
This research contributes to the growing body of evidence on the systemic effects of SARS-CoV-2, highlighting the importance of the gastrointestinal tract in COVID-19. Furthermore, these findings emphasize the need for targeted therapeutic strategies to restore a healthy gut microbiota in COVID-19 patients, potentially involving probiotics, prebiotics, and dietary interventions. Likewise, we believe that this study raises several important questions and provides a basis for conducting research on a larger cohort of patients in our region in order to fully understand the long-term effects of COVID-19 on gut health and the overall health of patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms12091800/s1, Table S1: Results of identification of bacterial and fungal isolates from fecal samples using MALDI-TOF MS.

Author Contributions

Conceptualization, A.P., V.V., V.P., and G.K.; methodology, A.P., B.M., D.M., and I.H.C.; software, V.V.; validation, A.P., B.M., and D.M.; formal analysis, A.P., G.K., and V.V.; investigation, A.P., G.K., J.R., M.D., A.K., T.P., and D.B.; data curation, A.P., G.K., V.V., J.R., A.K. T.P., and D.B.; writing—original draft preparation, A.P., G.K., V.V.; writing—review and editing, V.P., M.R., and I.H.C.; visualization, J.R., M.D., T.P., and M.R.; supervision, V.P. and I.H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Provincial Secretariat for Higher Education and Scientific Research within the framework of the Short-term project of special interest for sustainable development in the AP Vojvodina, Republic of Serbia, grant number 142-451-2390/2022-01/01.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors would like to thank all study participants for their involvement in this research. Additionally, the authors would like to express their gratitude to their colleagues at the Center of Virology and the Center of Microbiology of the Institute of Public Health of Vojvodina, Novi Sad, Serbia for their administrative and technical support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sender, R.; Fuchs, S.; Milo, R. Revised Estimates for the Number of Human and Bacteria Cells in the Body. PLoS Biol. 2016, 14, e1002533. [Google Scholar] [CrossRef]
  2. Hussain, I.; Cher, G.L.Y.; Abid, M.A.; Abid, M.B. Role of Gut Microbiome in COVID-19: An Insight Into Pathogenesis and Therapeutic Potential. Front. Immunol. 2021, 12, 765965. [Google Scholar] [CrossRef]
  3. Afzaal, M.; Saeed, F.; Shah, Y.A.; Hussain, M.; Rabail, R.; Socol, C.T.; Hassoun, A.; Pateiro, M.; Lorenzo, J.M.; Rusu, A.V.; et al. Human gut microbiota in health and disease: Unveiling the relationship. Front. Microbiol. 2022, 26, 999001. [Google Scholar] [CrossRef]
  4. Syromyatnikov, M.; Nesterova, E.; Gladkikh, M.; Smirnova, Y.; Gryaznova, M.; Popov, V. Characteristics of the Gut Bacterial Composition in People of Different Nationalities and Religions. Microorganisms 2022, 18, 1866. [Google Scholar] [CrossRef]
  5. Hillman, E.T.; Lu, H.; Yao, T.; Nakatsu, C.H. Microbial Ecology along the Gastrointestinal Tract. Microbes Environ. 2017, 32, 300–313. [Google Scholar] [CrossRef]
  6. Hou, K.; Wu, Z.X.; Chen, X.Y.; Wang, J.Q.; Zhang, D.; Xiao, C.; Zhu, D.; Koya, J.B.; Wei, L.; Li, J.; et al. Microbiota in health and diseases. Signal Transduct. Target. Ther. 2022, 7, 135. [Google Scholar] [CrossRef]
  7. Wu, H.J.; Wu, E. The role of gut microbiota in immune homeostasis and autoimmunity. Gut Microbes 2012, 3, 4–14. [Google Scholar] [CrossRef]
  8. Wiertsema, S.P.; van Bergenhenegouwen, J.; Garssen, J.; Knippels, L.M.J. The Interplay between the Gut Microbiome and the Immune System in the Context of Infectious Diseases throughout Life and the Role of Nutrition in Optimizing Treatment Strategies. Nutrients 2021, 13, 886. [Google Scholar] [CrossRef]
  9. Ng, Q.X.; Yaow, C.Y.L.; Moo, J.R.; Koo, S.W.K.; Loo, E.X.; Siah, K.T.H. A systematic review of the association between environmental risk factors and the development of irritable bowel syndrome. J. Gastroenterol. Hepatol. 2024. [Google Scholar] [CrossRef]
  10. 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, 70–89. [Google Scholar] [CrossRef]
  11. Belkaid, Y.; Hand, T.W. Role of the microbiota in immunity and inflammation. Cell 2014, 157, 121–141. [Google Scholar] [CrossRef]
  12. Madhogaria, B.; Bhowmik, P.; Kundu, A. Correlation between human gut microbiome and diseases. Infect. Med. 2022, 1, 180–191. [Google Scholar] [CrossRef]
  13. 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, 944–955.e8. [Google Scholar] [CrossRef]
  14. Giobbe, G.G.; Bonfante, F.; Jones, B.C.; Gagliano, O.; Luni, C.; Zambaiti, E.; Perin, S.; Laterza, C.; Busslinger, G.; Stuart, H.; et al. SARS-CoV-2 infection and replication in human gastric organoids. Nat. Commun. 2021, 12, 6610. [Google Scholar] [CrossRef]
  15. Newsome, R.C.; Gauthier, J.; Hernandez, M.C.; Abraham, G.E.; Robinson, T.O.; Williams, H.B.; Sloan, M.; Owings, A.; Laird, H.; Christian, T.; et al. The gut microbiome of COVID-19 recovered patients returns to uninfected status in a minority-dominated United States cohort. Gut Microbes 2021, 1, 1–15. [Google Scholar] [CrossRef]
  16. Yeoh, Y.K.; Zuo, T.; Lui, G.C.; 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, 698–706. [Google Scholar] [CrossRef]
  17. WHO. World Health Organization Tracking SARS-CoV-2 Variants. Available online: https://www.who.int/en/activities/tracking-SARS-CoV-2-variants (accessed on 26 June 2024).
  18. Wills, C.P.; Perez, B.; Moore, J. Coronavirus Disease 2019: Past, Present, and Future. Emerg. Med. Clin. N. Am. 2024, 42, 415–442. [Google Scholar] [CrossRef]
  19. Infectious Diseases Society of America. “COVID-19 Variant Update”. Available online: https://www.idsociety.org/covid-19-real-time-learning-network/diagnostics/covid-19-variant-update/#/+/0/publishedDate_na_dt/desc/,GISAID (accessed on 26 June 2024).
  20. Said, E.A.; Al-Reesi, I.; Al-Shizawi, N.; Jaju, S.; Al-Balushi, M.S.; Koh, C.Y.; Al-Jabri, A.A.; Jeyaseelan, L. Defining IL-6 levels in healthy individuals: A meta-analysis. J. Med. Virol. 2021, 93, 3915–3924. [Google Scholar] [CrossRef] [PubMed]
  21. Sarris, A.H.; Kliche, K.O.; Pethambaram, P.; Preti, A.; Tucker, S.; Jackow, C.; Messina, O.; Pugh, W.; Hagemeister, F.B.; McLaughlin, P.; et al. Interleukin-10 levels are often elevated in serum of adults with Hodgkin’s disease and are associated with inferior failure-free survival. Ann. Oncol. 1999, 10, 433–440. [Google Scholar] [CrossRef] [PubMed]
  22. Alecu, M.; Geleriu, L.; Coman, G.; Gălăţescu, L. The interleukin-1, interleukin-2, interleukin-6 and tumour necrosis factor alpha serological levels in localised and systemic sclerosis. Rom. J. Intern. Med. 1998, 36, 251–259. [Google Scholar] [PubMed]
  23. OurWorldInData. SARS-CoV-2 Sequences by Variant, Serbia. 21 November 2022. Available online: https://ourworldindata.org/grapher/covid-variants-bar?country=~SRB (accessed on 26 June 2024).
  24. Natarajan, A.; Zlitni, S.; Brooks, E.F.; Vance, S.E.; Dahlen, A.; Hedlin, H.; Park, R.M.; Han, A.; Schmidtke, D.T.; Verma, R.; et al. Gastrointestinal symptoms and fecal shedding of SARS-CoV-2 RNA suggest prolonged gastrointestinal infection. Med 2022, 3, 371–387.e9. [Google Scholar] [CrossRef]
  25. Antia, A.; Alvarado, D.M.; Zeng, Q.; Casorla-Perez, L.A.; Davis, D.L.; Sonnek, N.M.; Ciorba, M.A.; Ding, S. SARS-CoV-2 Omicron BA.1 Variant Infection of Human Colon Epithelial Cells. Viruses 2024, 16, 634. [Google Scholar] [CrossRef] [PubMed]
  26. Moreno-Corona, N.C.; López-Ortega, O.; Pérez-Martínez, C.A.; Martínez-Castillo, M.; De Jesús-González, L.A.; León-Reyes, G.; León-Juárez, M. Dynamics of the Microbiota and Its Relationship with Post-COVID-19 Syndrome. Int. J. Mol. Sci. 2023, 24, 14822. [Google Scholar] [CrossRef] [PubMed]
  27. Hasan, N.; Yang, H. Factors affecting the composition of the gut microbiota, and its modulation. PeerJ 2019, 7, e7502. [Google Scholar] [CrossRef] [PubMed]
  28. Bardi, T.; Pintado, V.; Gomez-Rojo, M.; Escudero-Sanchez, R.; Azzam Lopez, A.; Diez-Remesal, Y.; Martinez Castro, N.; Ruiz-Garbajosa, P.; Pestaña, D. Nosocomial Infections Associated to COVID-19 in the Intensive Care Unit: Clinical Characteristics and Outcome. Eur. J. Clin. Microbiol. Infect. Dis. 2021, 40, 495–502. [Google Scholar] [CrossRef] [PubMed]
  29. Bartoletti, M.; Giannella, M.; Scudeller, L.; Tedeschi, S.; Rinaldi, M.; Bussini, L.; Fornaro, G.; Pascale, R.; Pancaldi, L.; Pasquini, Z.; et al. Development and Validation of a Prediction Model for Severe Respiratory Failure in Hospitalized Patients With SARS-CoV-2 Infection: A Multicenter Cohort Study (PREDI-CO Study). Clin. Microbiol. Infect. 2020, 26, 1545–1553. [Google Scholar] [CrossRef] [PubMed]
  30. Galeeva, J.S.; Fedorov, D.E.; Starikova, E.V.; Manolov, A.I.; Pavlenko, A.V.; Selezneva, O.V.; Klimina, K.M.; Veselovsky, V.A.; Morozov, M.D.; Yanushevich, O.O.; et al. Microbial Signatures in COVID-19: Distinguishing Mild and Severe Disease via Gut Microbiota. Biomedicines 2024, 12, 996. [Google Scholar] [CrossRef] [PubMed]
  31. Gaibani, P.; D’Amico, F.; Bartoletti, M.; Lombardo, D.; Rampelli, S.; Fornaro, G.; Coladonato, S.; Siniscalchi, A.; Re, M.C.; Viale, P.; et al. The Gut Microbiota of Critically Ill Patients with COVID-19. Front. Cell Infect. Microbiol. 2021, 11, 670424. [Google Scholar] [CrossRef] [PubMed]
  32. Tariq, F.N.; Shafiq, M.; Khawar, N.; Habib, G.; Gul, H.; Hayat, A.; Rehman, M.U.; Moussa, I.M.; Mahmoud, E.A.; Elansary, H.O. The functional repertoire of AmpR in the AmpC β-lactamase high expression and decreasing β-lactam and aminoglycosides resistance in ESBL Citrobacter freundii. Heliyon 2023, 9, e19486. [Google Scholar] [CrossRef]
  33. Chiopu, P.; Toc, D.A.; Colosi, I.A.; Costache, C.; Ruospo, G.; Berar, G.; Gălbău, S.G.; Ghilea, A.C.; Botan, A.; Pană, A.G.; et al. An Overview of the Factors Involved in Biofilm Production by the Enterococcus Genus. Int. J. Mol. Sci. 2023, 24, 11577. [Google Scholar] [CrossRef]
  34. Devaux, C.A.; Lagier, J.C.; Raoult, D. New insights into the physiopathology of COVID-19: SARS-CoV-2-associated gastrointestinal illness. Front. Med. 2021, 8, 640073. [Google Scholar] [CrossRef]
  35. Rafiqul Islam, S.M.; Foysal, M.J.; Hoque, M.N.; Mehedi, H.H.; Rob, M.A.; Salauddin, A.; Tanzina, A.Y.; Biswas, S.; Noyon, S.H.; Siddiki, A.Z.; et al. Dysbiosis of oral and gut microbiomes in SARS-CoV-2 infected patients in Bangladesh: Elucidating the role of opportunistic gut microbes. Front. Med. 2022, 9, 821777. [Google Scholar] [CrossRef]
  36. Toc, D.A.; Mihaila, R.M.; Botan, A.; Bobohalma, C.N.; Risteiu, G.A.; Simut-Cacuci, B.N.; Steorobelea, B.; Troanca, S.; Junie, L.M. Enterococcus and COVID-19: The Emergence of a Perfect Storm? Int. J. Transl. Med. 2022, 2, 220–229. [Google Scholar] [CrossRef]
  37. Giacobbe, D.R.; Labate, L.; Tutino, S.; Baldi, F.; Russo, C.; Robba, C.; Ball, L.; Dettori, S.; Marchese, A.; Dentone, C.; et al. Enterococcal bloodstream infections in critically ill patients with COVID-19: A case series. Ann. Med. 2021, 53, 1779–1786. [Google Scholar] [CrossRef]
  38. Naveen, K.V.; Saravanakumar, K.; Sathiyaseelan, A.; MubarakAli, D.; Wang, M.H. Human fungal infection, immune response, and clinical challenge—A perspective during COVID-19 pandemic. Appl. Biochem. Biotechnol. 2022, 194, 4244–4257. [Google Scholar] [CrossRef]
  39. Maeda, Y.; Motooka, D.; Kawasaki, T.; Oki, H.; Noda, Y.; Adachi, Y.; Niitsu, T.; Okamoto, S.; Tanaka, K.; Fukushima, K.; et al. Longitudinal alterations of the gut mycobiota and microbiota on COVID-19 severity. BMC Infect. Dis. 2022, 22, 572. [Google Scholar] [CrossRef]
  40. Granata, G.; Schiavone, F.; Pipitone, G.; Taglietti, F.; Petrosillo, N. Antibiotics Use in COVID-19 Patients: A Systematic Literature Review. J. Clin. Med. 2022, 11, 7207. [Google Scholar] [CrossRef]
  41. Vestad, B.; Ueland, T.; Lerum, T.V.; Dahl, T.B.; Holm, K.; Barratt-Due, A.; Kåsine, T.; Dyrhol-Riise, A.M.; Stiksrud, B.; Tonby, K.; et al. Respiratory Dysfunction Three Months after Severe COVID-19 Is Associated with Gut Microbiota Alterations. J. Intern. Med. 2022, 291, 801–812. [Google Scholar] [CrossRef]
  42. Dhar, D.; Mohanty, A. Gut microbiota and COVID-19- possible link and implications. Virus Res. 2020, 285, 198018. [Google Scholar] [CrossRef]
  43. Morais, A.H.A.; Passos, T.S.; Maciel, B.L.L.; da Silva-Maia, J.K. Can Probiotics and Diet Promote Beneficial Immune Modulation and Purine Control in Coronavirus Infection? Nutrients 2020, 12, 1737. [Google Scholar] [CrossRef] [PubMed]
  44. Maciel-Fiuza, M.F.; Muller, G.C.; Campos, D.M.S.; do Socorro Silva Costa, P.; Peruzzo, J.; Bonamigo, R.R.; Veit, T.; Vianna, F.S.L. Role of gut microbiota in infectious and inflammatory diseases. Front. Microbiol. 2023, 14, 1098386. [Google Scholar] [CrossRef] [PubMed]
  45. Liu, Y.; Liao, W.; Wan, L.; Xiang, T.; Zhang, W. Correlation between Relative Nasopharyngeal Virus RNA Load and Lymphocyte Count Disease Severity in Patients with COVID-19. Viral Immunol. 2021, 34, 330–335. [Google Scholar] [CrossRef] [PubMed]
  46. Ruff, W.E.; Greiling, T.M.; Kriegel, M.A. Host–microbiota interactions in immune-mediated diseases. Nat. Rev. Microbiol. 2020, 18, 521–538. [Google Scholar] [CrossRef] [PubMed]
  47. Lutfi, L.L.; Shaaban, M.I.; Elshaer, S.L. Vitamin D and vitamin K1 as novel inhibitors of biofilm in Gram-negative bacteria. BMC Microbiol. 2024, 24, 173. [Google Scholar] [CrossRef] [PubMed]
  48. Zsichla, L.; Müller, V. Risk Factors of Severe COVID-19: A Review of Host, Viral and Environmental Factors. Viruses 2023, 15, 175. [Google Scholar] [CrossRef] [PubMed]
  49. Ng, Q.X.; Yau, C.E.; Yaow, C.Y.L.; Chong, R.I.H.; Chong, N.Z.; Teoh, S.E.; Lim, Y.L.; Soh, A.Y.S.; Ng, W.K.; Thumboo, J. What Has Longitudinal ‘Omics’ Studies Taught Us about Irritable Bowel Syndrome? A Systematic Review. Metabolites 2023, 13, 484. [Google Scholar] [CrossRef]
  50. Duncanson, K.; Williams, G.; Hoedt, E.C.; Collins, C.E.; Keely, S.; Talley, N.J. Diet-microbiota associations in gastrointestinal research: A systematic review. Gut Microbes 2024, 16, 2350785. [Google Scholar] [CrossRef]
Figure 1. Overview of the study procedure and general sampling process.
Figure 1. Overview of the study procedure and general sampling process.
Microorganisms 12 01800 g001
Figure 2. Fecal microbial composition isolated in culture and diversity in samples from SARS-CoV-2 positive patients. Fecal samples were taken at two different time points: within 5 days of the first positive test result for COVID-19 (marked as ≤5 days) and 30–35 days after the cessation of disease symptoms (marked as 30–35 days). (A): A total of 272 microbial isolates were identified at the onset of infection (≤5 days), while 294 isolates were identified 30–35 days after the cessation of disease symptoms. Stacked bar plots summarize the relative abundance of microbial taxa at the phylum level (B) and genus level (C) in samples taken at the two time points. Gut microbiota diversity was assessed using Shannon’s index (D) and Simpson’s index (E). Both indices showed no statistically significant differences in diversity between the two groups (p > 0.05; t-test).
Figure 2. Fecal microbial composition isolated in culture and diversity in samples from SARS-CoV-2 positive patients. Fecal samples were taken at two different time points: within 5 days of the first positive test result for COVID-19 (marked as ≤5 days) and 30–35 days after the cessation of disease symptoms (marked as 30–35 days). (A): A total of 272 microbial isolates were identified at the onset of infection (≤5 days), while 294 isolates were identified 30–35 days after the cessation of disease symptoms. Stacked bar plots summarize the relative abundance of microbial taxa at the phylum level (B) and genus level (C) in samples taken at the two time points. Gut microbiota diversity was assessed using Shannon’s index (D) and Simpson’s index (E). Both indices showed no statistically significant differences in diversity between the two groups (p > 0.05; t-test).
Microorganisms 12 01800 g002
Figure 3. Correlation between the percentage of opportunistic bacteria and the COVID-19 symptom duration. The percentage of opportunistic bacteria in stool samples from COVID-19 patients was measured at two different time points: (A) at the onset of the disease (Sample 1) and (B) 30–35 days after the cessation of disease symptoms (Sample 2).
Figure 3. Correlation between the percentage of opportunistic bacteria and the COVID-19 symptom duration. The percentage of opportunistic bacteria in stool samples from COVID-19 patients was measured at two different time points: (A) at the onset of the disease (Sample 1) and (B) 30–35 days after the cessation of disease symptoms (Sample 2).
Microorganisms 12 01800 g003
Figure 4. Levels of inflammatory biomarkers and SARS-CoV-2 IgG antibodies in COVID-19 positive patients. The levels of various biomarkers were measured in serum samples at two different time points: ≤5 days and 30–35 days. (A) The bar plot shows CRP levels; (B) The dot plot compares IL-6 levels; (C) The dot plot for IL-1α; (D) The dot plot for IL-10 levels. (E) The dot plot for SARS-CoV-2 IgG. The plots include individual data points, mean values (central lines), and variability indicators (mean 95% CI). Statistically significant differences between the two groups were observed for the levels of CRP and SARS-CoV-2 IgG according to the Wilcoxon signed-rank test, *** p < 0.001; ns: not significant.
Figure 4. Levels of inflammatory biomarkers and SARS-CoV-2 IgG antibodies in COVID-19 positive patients. The levels of various biomarkers were measured in serum samples at two different time points: ≤5 days and 30–35 days. (A) The bar plot shows CRP levels; (B) The dot plot compares IL-6 levels; (C) The dot plot for IL-1α; (D) The dot plot for IL-10 levels. (E) The dot plot for SARS-CoV-2 IgG. The plots include individual data points, mean values (central lines), and variability indicators (mean 95% CI). Statistically significant differences between the two groups were observed for the levels of CRP and SARS-CoV-2 IgG according to the Wilcoxon signed-rank test, *** p < 0.001; ns: not significant.
Microorganisms 12 01800 g004
Table 1. Demographic and clinical characteristics of the study cohort in relation to SARS-CoV-2 fecal shedding in samples collected at the onset of infection.
Table 1. Demographic and clinical characteristics of the study cohort in relation to SARS-CoV-2 fecal shedding in samples collected at the onset of infection.
ParametersNSARS-CoV-2 Fecal Shedding (n = 33; 100%) 1Without
SARS-CoV-2
Fecal Shedding (n = 17; 100%) 1
p-Value
Gender
Male2212 (36.4%)10 (58.8%)0.13
Female2821 (63.6%)7 (41.2%)
Age group
18–2922 (6.1%)00.386
30–3995 (15.2%)4 (23.5%)
40–491210 (30.3%)2 (11.8%)
50–59158 (24.2%)7 (41.2%)
60+128 (24.2%)4 (23.5%)
COVID-19 severity
Asymptomatic21 (3.0%)1 (5.9%)0.684
Mild illness3223 (69.7%)9 (52.9%)
Moderate illness127 (21.2%)5 (29.4%)
Severe illness42 (6.1%)2 (11.8%)
COVID-19 duration of symptoms (days)
1–522 (6.1%)00.538
6–9106 (18.2%)4 (23.5%)
10–133220 (60.6%)12 (70.6%)
≥1465 (15.2%)1 (5.9%)
GI symptoms
yes2517 (51.5%)8 (47.1%)0.765
no2516 (48.5%)9 (52.9%)
COVID-19 vaccine status (at least two doses received prior to infection)
yes4227 (81.8%)15 (88.2%)0.558
no86 (18.2%)2 (11.8%)
Prior COVID-19 infection
yes158 (24.2%)7 (41.2%)0.216
no3525 (75.8%)10 (58.8%)
1 Statistical analysis was performed on samples collected at the onset of infection (≤5 days), given that only one patient continued to shed SARS-CoV-2 in the sample collected 30–35 days after the cessation of disease symptoms.
Table 2. Statistical analysis overview of the presence and absence of various bacterial and fungal taxa over time.
Table 2. Statistical analysis overview of the presence and absence of various bacterial and fungal taxa over time.
Bacterial and Fungal Taxa
Order/Genus
Bacterial Taxa (Presence or Absence)Mc
Nemar
Chi-Squared
p-Value
5 Days/30–35 Days
+/+
n (%)
+/−
n (%)
−/+
n (%)
−/−
n (%)
Genus: Bacillus 1, Sporosarcina, Lysinibacillus, Psychrobacillus, Listeria, Paenibacillus, Peribacillus8
(16%)
17
(34%)
6
(12%)
19
(38%)
4.3480.037
Genus: Lactobacillus, Carnobacterium, Filifactor Limosilactobacillus, Streptococcus, Ligilactobacillus, Tissierella6
(12%)
5
(10%)
15
(30%)
24
(48%)
4.0500.044
Genus: Enterococcus31
(62%)
11
(22%)
7
(14%)
1
(2%)
0.5000.479
Genus: Pseudomonas2
(2%)
7
(14%)
2
(4%)
39
(78%)
1.7780.182
Genus: Klebsiella, Citrobacter, Enterobacter, Proteus Moellerella, Yersinia, Serratia8
(14%)
9
(18%)
5
(10%)
28
(56%)
0.6430.422
Genus: Escherichia14
(26%)
3
(6%)
3
(4%)
30
(60%)
0.1670.683
Genus: Hafnia8
(16%)
5
(10%)
10
(20%)
27
(54%)
1.0670.302
Genus: Raoultella6
(12%)
2
(4%)
0
(0.0%)
42
(84%)
1.1250.289
Genus: Corynebacterium, Micrococcus, Paenarthrobacter, Microbacterium, Brevibacterium6
(12%)
4
(8%)
1
( 2%)
39
(78%)
0.1000.752
Fungi: Candida, Geotrichum, Pichia, Mucor, Rhodotorula, Aspergillus, Cryptococcus4
(8%)
2
(4%)
10
(20%)
34
(68%)
4.0830.043
1 The genus Bacillus was analyzed excluding Bacillus cereus, which is considered an opportunistic pathogen. The table analyzes four categories indicating the presence or absence of bacterial and fungal taxa at two different time points: at the first measurement (within 5 days of onset) and at the second measurement (30–35 days after the cessation of disease symptoms). A mark of “+/+” indicates the taxon was observed in both samples; “+/−” indicates the taxon was found only in the first sample; “−/+” indicates the taxon was found only in the second sample; “−/−” indicates the taxon was not found in either sample. Results are presented along with McNemar’s chi-square test and Bonferroni corrected significance, considering a p-value below 0.005 as statistically significant
Table 3. Statistical analysis of the detection rates of opportunistic bacteria stratified by demographic and clinical characteristics.
Table 3. Statistical analysis of the detection rates of opportunistic bacteria stratified by demographic and clinical characteristics.
ParametersN (%)Detection Rate of
Opportunistic Taxa (%);
≤5 Days
Detection Rate of
Opportunistic Taxa (%);
30–35 Day
p-Value
Mean %SDMean %SD
Total5060.2627.1049.2027.940.048 *
Gender
Male 22 (44%)62.6719.3948.6228.340.062
Female28 (56%)58.3632.1349.6628.130.292
Age group
18–292 (4%)26.679.4354.1729.460.500
30–399 (18%)47.1429.2849.6026.860.865
40–4912 (24%)53.1020.2749.6033.950.752
50–5915 (30%)69.9224.9851.1119.770.036 *
60+12 (24%)70.7728.0345.2834.600.063
GI symptoms a
Yes25 (50%)57.7128.6250.6725.170.393
No25 (50%)62.8025.8247.7330.910.054
COVID-19 severity
Asymptomatic2 (4%)65.0021.2125.0035.360.500
Mild illness32 (64%)55.8625.6649.4229.460.364
Moderate illness12 (24%)75.5124.5456.1925.270.094
Severe illness4 (8%)47.3238.1738.5716.860.652
COVID-19-duration of symptoms (days)
1–52 (4%)85.7120.2035.0021.210.333
6–910 (20%)68.2722.2260.2523.320.378
10–1432 (64%)55.0228.4647.8128.180.328
≥14 6 (12%)66.3524.4342.9435.650.213
Comorbidities b
Yes25 (50%)60.8628.2348.5930.120.158
No25 (50%)59.6526.4949.8126.180.182
COVID-19 vaccine status (at least two doses)
Yes42 (84%)60.6326.3948.9228.230.056
No8 (16%)58.2832.5150.6928.130.621
Prior COVID-19 Infection
Yes15 (30%)53.9526.6347.1724.920.436
No35 (70%)62.9627.2350.0729.440.071
SARS-CoV-2 RNA fecal shedding c
Yes18 (36%)54.8430.8849.0726.820.588
No32 (64%)63.3024.7349.2728.970.033 *
Vitamin D Supplements
Yes17 (34%)59.3527.6260.7828.240.883
No33 (66%)60.7327.2543.2326.240.011 *
a—Gastrointestinal (GI) symptoms included: nausea, vomiting, diarrhea; b—Comorbidities (hypertension, diabetes, chronic lung disease, chronic heart disease, chronic kidney disease, or autoimmune disease); c—SARS-CoV-2 RNA fecal shedding observed in the sample collected within the initial 5 days of the onset of COVID-19. Statistical significance was assessed by * paired t-test.
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.

Share and Cite

MDPI and ACS Style

Patić, A.; Kovačević, G.; Vuković, V.; Hrnjaković Cvjetković, I.; Ristić, M.; Milosavljević, B.; Medić, D.; Djilas, M.; Radovanov, J.; Kovačević, A.; et al. Analysis of Cultured Gut Microbiota Using MALDI-TOF MS in COVID-19 Patients from Serbia during the Predominance of the SARS-CoV-2 Omicron Variant. Microorganisms 2024, 12, 1800. https://doi.org/10.3390/microorganisms12091800

AMA Style

Patić A, Kovačević G, Vuković V, Hrnjaković Cvjetković I, Ristić M, Milosavljević B, Medić D, Djilas M, Radovanov J, Kovačević A, et al. Analysis of Cultured Gut Microbiota Using MALDI-TOF MS in COVID-19 Patients from Serbia during the Predominance of the SARS-CoV-2 Omicron Variant. Microorganisms. 2024; 12(9):1800. https://doi.org/10.3390/microorganisms12091800

Chicago/Turabian Style

Patić, Aleksandra, Gordana Kovačević, Vladimir Vuković, Ivana Hrnjaković Cvjetković, Mioljub Ristić, Biljana Milosavljević, Deana Medić, Milan Djilas, Jelena Radovanov, Aleksandra Kovačević, and et al. 2024. "Analysis of Cultured Gut Microbiota Using MALDI-TOF MS in COVID-19 Patients from Serbia during the Predominance of the SARS-CoV-2 Omicron Variant" Microorganisms 12, no. 9: 1800. https://doi.org/10.3390/microorganisms12091800

APA Style

Patić, A., Kovačević, G., Vuković, V., Hrnjaković Cvjetković, I., Ristić, M., Milosavljević, B., Medić, D., Djilas, M., Radovanov, J., Kovačević, A., Pustahija, T., Balać, D., & Petrović, V. (2024). Analysis of Cultured Gut Microbiota Using MALDI-TOF MS in COVID-19 Patients from Serbia during the Predominance of the SARS-CoV-2 Omicron Variant. Microorganisms, 12(9), 1800. https://doi.org/10.3390/microorganisms12091800

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop