Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

16S rRNA gene sequencing of rectal swab in patients affected by COVID-19

  • Antonio Mazzarelli ,

    Contributed equally to this work with: Antonio Mazzarelli, Maria Letizia Giancola

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Maria Letizia Giancola ,

    Contributed equally to this work with: Antonio Mazzarelli, Maria Letizia Giancola

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Anna Farina,

    Roles Formal analysis

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Luisa Marchioni,

    Roles Methodology

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Martina Rueca,

    Roles Methodology

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Cesare Ernesto Maria Gruber,

    Roles Formal analysis

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Barbara Bartolini,

    Roles Funding acquisition, Supervision

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Tommaso Ascoli Bartoli,

    Roles Formal analysis

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Gaetano Maffongelli,

    Roles Supervision

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Maria Rosaria Capobianchi,

    Roles Supervision

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Giuseppe Ippolito,

    Roles Supervision

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Antonino Di Caro ,

    Roles Supervision, Writing – review & editing

    antonino.dicaro@inmi.it

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Emanuele Nicastri ,

    Roles Conceptualization, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation National Institute for Infectious Diseases, INMI “Lazzaro Spallanzani”, IRCCS, Rome, Italy

  • Valerio Pazienza ,

    Roles Conceptualization, Formal analysis, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation Division of Gastroenterology, Fondazione-IRCCS “Casa Sollievo della Sofferenza” Hospital, S. Giovanni Rotondo (FG), Italy

  • INMI COVID-19 study group

    Membership of the Collaborators Members of the National Institute for Infectious Diseases (INMI) COVID-19 study group is provided in the acknowledgments.

Abstract

COronaVIrus Disease-2019 (COVID-19) is a pandemic respiratory infection caused by a new betacoronavirus, the Severe Acute Respiratory Syndrome-CoronaVirus-2 (SARS-CoV-2). Few data are reported on the gut microbiota in COVID-19 patients. 16S rRNA gene sequencing was performed to reveal an altered composition of the gut microbiota in patients with COVID-19 pneumonia admitted in intensive care unit (ICU) (i-COVID19), or in infectious disease wards (w-COVID19) as compared to controls (CTRL). i-COVID19 patients showed a decrease of Chao1 index as compared to CTRL and w-COVID19 patients indicating that patients in ICU displayed a lower microbial richness while no change was observed as for Shannon Index. At the phylum level, an increase of Proteobacteria was detected in w-COVID19 patients as compared to CTRL. A decrease of Fusobacteria and Spirochetes has been found, with the latter decreased in i-COVID19 patients as compared to CTRL. Significant changes in gut microbial communities in patients with COVID-19 pneumonia with different disease severity compared to CTRL have been identified. Our preliminary data may provide valuable information and promising biomarkers for the diagnosis of the disease and, when validated in larger cohort, it could facilitate the stratification of patients based on the microbial signature.

1. Introduction

COronaVIrus Disease-2019 (COVID-19) is a pandemic respiratory infection caused by a new betacoronavirus infection, the Severe Acute Respiratory Syndrome-CoronaVirus-2 (SARS-CoV-2). It may progress rapidly to acute respiratory distress syndrome with remarkable morbidity and mortality [1]. However, SARS-CoV-2 can be detected in specimens from different sites and therefore it could potentially be transmitted in other ways than respiratory droplets [2]. In this regard, recent studies reported that SARS-CoV-2 RNA has been found in anal swabs, meaning that the virus could potentially be transmitted also through oral-fecal route [3]. These findings might suggest that other organs apart from lung might be additional sites for virus entry, repository and/or replication. Gastrointestinal (GI) symptoms, such as diarrhea (2%-10.1%), nausea and vomiting (1%-3.6%), are not very common at present in COVID-19 patients [4]. Nevertheless, an important proportion of patients were observed during the global pandemic showing atypical gastrointestinal symptoms [5].

Recently, several studies have demonstrated that respiratory infections are associated with changes of the gut microbiota composition [6, 7]. Typically, Bacteroidetes and Firmicutes prevail in the gut microbiota while potentially pathogenic species, such as some of those belonging to the phylum Proteobacteria, are present in a minor percentage [8, 9]. Current studies evaluate the relationship between the lung and the GI microbiota but this connection is not completely understood [10]. Patients with respiratory infections generally have gut dysfunction, which is related to a more severe clinical course of the disease. This phenomenon can also be observed in COVID-19 patients [5].

The gut microbiota has been shown to affect pulmonary health through cross-talk between the gut microbiota and the lungs, which is referred to as the “gut-lung axis” [11]. However, the “gut-lung axis” is supposed to be bidirectional: the gut microbiota, through microbial products and immune-modulators released upon recognition of commensals and pathogens by intestinal immune cells, can regulate lung immunity, influence the lung microbiota, and vice versa [12, 13]. Previous studies have shown that the modulation of the gut microbiota can reduce the severity of enteritis and ventilator-associated pneumonia by interacting with early replication of the viruses in the pulmonary epithelium, as in the case of influenza virus [14].

Angiotensin-converting enzyme 2 (ACE2) is the main receptor of SARS-CoV [14] and of SARS-CoV-2 [15]. This receptor is highly expressed in both the respiratory tract and GI, so it is possible to consider that SARS-CoV-2 uses ACE2 receptor to get into both body districts [16, 17].

All these virus characteristics raise a remarkable possibility that the pulmonary disease caused by SARS-CoV-2 may influence the gut microbiota [5].

In this pilot study, we performed the 16S RNA sequencing of fecal samples from COVID-19 in patients admitted to the National Institute for Infectious Disease (INMI) L. Spallanzani in Rome, Italy, between April 15, 2020 and May 31, 2020.

2. Methods

2.1. Study design

All the subjects agreed to participate according to the ethical guidelines of the 2013 Declaration of Helsinki signing an informed consent under the Ethical committee of the National Institute for Infectious Diseases Lazzaro Spallanzani—IRCCS approval number (n. 9/2020; n. 3 23.12.2019) and followed the same pre-analytical and analytical procedures, including fecal samples collection and storage. Data were analyzed anonymously.

From April 15 2020 and May 31 2020, rectal swabs were collected from patients hospitalized to INMI Spallanzani in Rome with confirmed or suspected SARS-CoV-2 infections. Fifteen out of 23 patients were affected by COVID-19 while 8 patients were COVID-19 negative.

All patients had pneumonia and were classified in three groups: a) COVID-19 patients with nose-pharyngeal swab positive for SARS-CoV-2, >18 years of age admitted in infectious disease wards at the time of the rectal swab execution (w-COVID19); b) COVID-19 patients with nose-pharyngeal swab positive for SARS-CoV-2, >18 years of age admitted in ICU (i-COVID19); c) ‘controls’: patients admitted in the same time period with nose-pharyngeal swab negative for SARS-CoV-2 infections, hospitalized in ICU and/or in infectious disease ward (CTRL). We performed the rectal swab one or two days after the hospitalization.

Patient data, including laboratory test results and clinical manifestations, were obtained from Laboratory Information Systems (LIS) and clinical records.

Rectal swabs from COVID-19 patients were processed in the laboratory within 4 hours after collection or stored at -80°C until analysis.

2.2. DNA extraction

Samples were treated with 500 μl of Lysis Buffer ATL (QIAGEN, Hilden, Germany) at 56°C for 10 min with 20 μl Proteinase K (Darmstadt, Germany) before DNA extraction. Microbial DNA was extracted from 500 μl of sample using Qiasymphony automatic extractor (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol.

2.3. S sequencing and analysis

DNA samples were generated from PCR amplicons targeting the hypervariable regions V2, V4, V8 and V3-6, 7–9 of the 16S gene and libraries were processed using the Ion 16S metagenomics Kit. Ion Xpress Plus Fragment Library kit was used for libraries obtainment. Sequencing was performed on Ion 530 chip by Ion S5 sequencer (Ion Torrent-ThermoFisher Scientific).

The sequencing run has generated in total 16x10^6 reads with the 77% of high quality reads (21% low quality, 2% test fragments). Finally we obtained 5.5x10^5 reads per sample (reads length were 244 bp mean, 260 bp median and 289 bp mode) and all analyzed specimens showed a suitable library’s profile. The analysis was performed by 16S Metagenomics GAIA 2.0 software and DESeq2 package software. Sequence data generated as FASTQ files, were analyzed using the 16S Metagenomics GAIA 2.0 software which performs the quality control of the reads/pairs (i.e., trimming, clipping and adapter removal steps) through FastQC and BBDuk. The reads/pairs are mapped with BWA-MEM against the 16S databases (GAIA based on NCBI). Differential expression analysis using DESeq2 package to test for differential analysis by use of negative binomial generalized linear models was used. Only changes with FDR below 0.05 were considered significant. The percent similarity used to determine species and genus calls was 93% at genus, 97% at species. PCoA analysis was obtained with GAIA software based on Bray-Curtis dissimilarities.

2.4. Statistical analysis

To evaluate if any clinical or laboratory variables are significantly different between patient cohorts, one-way ANOVA test and Kruskal-Wallis rank sum test were performed in R environment (www.cran.r-project.org) using aov and kruskal.test functions, respectively. Venn diagrams were obtained with Venny 2.1.0.

3. Results

3.1. Study population

The study population included 23 hospitalized inpatients; 15 out of 23 were patients with confirmed SARS-CoV-2 infection, 9 w-COVID19 and 6 i-COVID19, 8 were CTRL (three hospitalized in ICU and five in floor). Clinical data of the study patients are shown in Table 1. Overall, thirteen patients (56%) were male; median age was 67 (IQR 44–83). All patients presented pneumonia (for our CTRL six out of eight had bacteria pneumonia while 2 CTRL had non-COVID-19 viral pneumonia) and none of them had diarrhea when the rectal swab was performed (one or two days after the hospitalization). ANOVA and Kruskal-Wallis analysis in Table 1, showed that the p-values of the variables between our patient cohorts was not significant except for ferritin level that is significantly lower in ICU patients.

thumbnail
Table 1. Clinical and serologic data of the 23 enrolled patients.

https://doi.org/10.1371/journal.pone.0247041.t001

Eleven patients (48%) were receiving antibiotic therapy one or at most two days before the rectal swab was collected: 5 w-COVID19, 3 i-COVID19 and 3 CTRL.

Nineteen patients (83%) had one or more comorbidities, mainly cardiovascular and brain disorders. A concomitant infection (5 patients) and/or colonization (5 patients, although harbor a potentially pathological bacterium in their intestines namely Enterococcus faecium and Enterococcus faecalis, without developing the disease and consequently without the need of therapies or further tests) was found in seven patients (30%), 5 of them were hospitalized in ICU.

3.2. Microbial richness and diversity indices in COVID-19 patients as compared to CTRL

To understand the gut microbiota alterations between w-COVID19 and i-COVID19 patients with CTRL, as a first step we evaluated the richness and Shannon indices among the different groups. As shown in Fig 1A, Chao1 index was significantly decreased in i-COVID19 as compared to w-COVID19 (p = 0.02) and CTRL (p = 0.006). The same trend was also observed for Shannon index without reaching the statistical significance (Fig 1B).

thumbnail
Fig 1.

Box-plots of Chao1 index of species richness (A) and Shannon index of species diversity (B) in w-COVID19, i-COVID19 patients and CTRL. Triangles indicate the medians and Q1 and Q3 are reported.

https://doi.org/10.1371/journal.pone.0247041.g001

Principal coordinates analysis (PCoA) was performed to cluster the microbial communities at the Family operational taxonomic unit (OTU) level based on Bray-Curtis distances. Fig 2 displayed distinct patterns among the three groups CTRL (red) w-COVID19 (green) and i-COVID19 (blue).

thumbnail
Fig 2. Principal coordinates plot (PCoA) based on Bray-Curtis distances at family level showing a clustering pattern among samples obtained from controls (red), w-COVID19 (green) and i-COVID19 (blue).

https://doi.org/10.1371/journal.pone.0247041.g002

3.3. Microbiota profiles of w-COVID19 and i-COVID19 patients as compared to CTRL

At the Phylum level, Proteobacteria were significantly increased in w-COVID19 patients as compared to CTRL (17.1% vs 11.3% respectively FDR = 0.03) while Spirochaetes and Fusobacteria were decreased (0% vs 0.08% FDR = 0.00 and 0.02% vs 0.04% FDR = 0.00 respectively). When comparing w-COVID19 and i-COVID19 patients’ microbiota no significant changes were observed at the Phylum level (Fig 3 and S1S3 Tables).

thumbnail
Fig 3. Microbiota composition of w-COVID19, i-COVID19 patients and CTRL, at the phylum level.

The mean value of all the detected taxa is represented.

https://doi.org/10.1371/journal.pone.0247041.g003

At the Family level, a number of potential pathogenic bacteria such as Peptostreptococcaceae, Enterobacteriaceae, Staphylococcaceae, Vibrionaceae, Aerococcaceae, Dermabacteraceae, Actinobacteria and others were increased in w-COVID19 as compared to CTRL while Nitrospiraceae, Propionibacteriaceae, Aeromonadaceae, Moraxellaceae, Mycoplasmataceae were significantly reduced together with others reported in Fig 3 and S4 Table. When considering the i-COVID19 as compared to CTRL, in addition to some bacteria in common with w-COVID19 patients (i.e. Staphylococcaceae, Aerococcaceae, Dermabacteraceae, Actinobacteria and so on Fig 4, and Fig 8A and S5 Table) Erysipelotrichaceae, Microbacteriaceae, Mycobacteriaceae, Pseudonocardiaceae, Brevibacteriaceae, and others reported in S5 Table were significantly increased while Carnobacteriaceae, Coriobacteriaceae and Mycoplasmataceae were significantly reduced.

thumbnail
Fig 4. Microbiota composition of w-COVID19, i-COVID19 patients and CTRL, at the family level.

The mean value of all the detected taxa is represented.

https://doi.org/10.1371/journal.pone.0247041.g004

Nevertheless Staphylococcaceae, Microbacteriaceae, Micrococcaceae, Pseudonocardiaceae, Erysipelotrichales and others reported in S6 Table were significantly higher in i-COVID-19 as compared to w-COVID19. Carnobacteriaceae, Pectobacteriaceae, Moritellaceae, Selenomonadaceae, Micromonosporaceae, Coriobacteriaceae and few others were significantly decreased in i-COVID19 as compared to w-COVID19. Individual microbiota profiles are provided in Figs 5 and 6 at the Phylum and Family level.

thumbnail
Fig 5. Taxonomy bar plot showing the individual microbiota profile at phylum levels.

https://doi.org/10.1371/journal.pone.0247041.g005

thumbnail
Fig 6. Taxonomy bar plot showing the individual microbiota profile at family levels.

https://doi.org/10.1371/journal.pone.0247041.g006

Moreover, unsupervised hierarchical analysis at the family level (Fig 7) revealed a characteristic microbial signature in CTRL segregated from that one of COVID-19 positive patients. Strikingly, a distinct profile can be distinguished between i-COVID19 and w-COVID19 with the latter being closer to CTRL. Microbiota analysis at lower taxonomic levels (genera and species, reported in S7S12 Tables).

thumbnail
Fig 7. Heatmap of one-way hierarchical clustering of differentially abundant families among the three cohorts.

A dual-color code counts for species up- (red) and down-represented (blue), respectively.

https://doi.org/10.1371/journal.pone.0247041.g007

3.4. Differences in microbial populations of w-COVID19 and i-COVID19 patients in comparison to CTRL

At lower taxonomic levels, many differences with CTRL group emerged in both i- and w-COVID19 patients. VENN diagrams showed the number of families (Fig 8A) genera (Fig 8B) and species (Fig 8C) shared or distinctive of the two different groups as compared to CTRL used as reference. Although w-COVID19 and i-COVID19 patients share a number of increased and decreased bacteria, a distinctive bacteria profile can be also observed when compared to CTRL (S13 Table).

thumbnail
Fig 8.

Venn diagrams showing the number of distinct and shared families (A), genera (B) and species (C) up and decreased between subjects grouped by w-COVID19, i-COVID19 patients as compared to CTRL used as referenced.

https://doi.org/10.1371/journal.pone.0247041.g008

4. Discussion

It is nowadays well recognized that virus infections can alter the host’s microbiota at different sites [18, 19], however is less clear whether changes of microbiota have direct or indirect effects i.e. limiting or promoting viral infections. Microbiota’s products such as short chain fatty acids, metabolites or bacteriocine may directly interact with viral particles to alter infectivity or responses to therapy [20, 21].

In our pilot study we demonstrated that SARS-CoV-2 infection is associated with major changes in gut microbiota profile of the patients. The main findings are the reduction of microbial richness in i-COVID19 as compared to CTRL and w-COVID19 indicating that patients in ICU displayed a lower microbial richness as measured by Chao1 index. For the Shannon index the same trend was also observed, but without reaching statistical significance. Our results are in line with those recently obtained by Zuo et al, in which enrichment of opportunistic pathogens and loss of beneficial bacteria was observed [22].

Moreover, ANOVA and Kruskal-Wallis analysis reveal that the variables between our patient cohorts was not noteworthy except for ferritin level that is significantly lower in ICU patients.

Ferritin is a marker of inflammation and the high levels of ferritin detected in i-COVID19 patients in comparison to w-COVID19, may be associate with a greater severity of the disease and adverse outcomes. Normally, ferritin is able to activate macrophages that when stimulated begin to secrete cytokines that at low concentrations, help to protect the body from viruses and bacteria. On the other hand, high levels of ferritin activate more macrophages that produce the so-called "cytokine storm" which can be lethal for the body [23].

As for gut microbiota, antibiotics use can obviously determine a further loss of heterogeneity and composition, leading to down regulation of beneficial symbionts and exacerbation of gut dysbiosis, and for this reason the avoidance of unnecessary antibiotics use in the treatment of viral pneumonia is strongly suggested, as antibiotics can eliminate beneficial bacteria and weaken the gut barrier [24]. Although an increase of pro-inflammatory and potential pathogenic bacteria such as Peptostreptococcaceae, Enterobacteriaceae, Staphylococcaceae, Vibrionaceae, Aerococcaceae, Dermabacteraceae, Actinobacteria [2527], is confirmed in w-COVID19 and i-COVID19 patients, with some of them found in both groups, hierarchical analysis shows a distinct profile between i- and w-COVID19 with the latter being closer to CTRL (Fig 4). Notably, a profound dysbiosis (Fig 5) was observed in one ward patient (90-year-old patient with diabetes, meningioma and osteoporosis in association with an increase of C-Reactive Protein (CRP) and lymphocytes), with a significant increase in Proteobacteria and a relevant reduction in Bacteriodetes, reflecting an important inflammatory state. Growing evidence has shown that perturbation of the gut microbial community may fuel blooms of otherwise low abundance and harmful bacteria which can further exacerbate the intestinal inflammation. Indeed, dysbiosis in the distal gut is often characterized by a decrease in the prevalence of strict anaerobes and an increased relative abundance of facultative anaerobic bacteria.

This could also contribute to the lower severity of symptoms of w-COVID19 as compared to patients admitted to ICU. This evidence is confirmed by our three CTRL patients hospitalzed in ICU that had a nasopharyngeal PCR negative for SARS-CoV2. The latter, showed an increase in Erysipelotrichaceae which are involved in inflammation-related disorders of the GI tract [28]. These data are in agreement with a previous study [29] where Erysipelotrichaceae were found associated with COVID-19 severity. Noteworthy a strong decrease of Faecalibacterium (which is associated to Crohn’s disease, colonrectal tumor, diabetes and non alcoholic steatohepatitis) [30] in i-COVID19 patients was detected together with a reduction of Ruminococcaceae, Clostridiaceae which are involved in Short Chain Fatty Acids (SCFAs) production among which butyrate presents potent antinflammatory properties [31]. Among the species (S9 Table), it is worth of note the decrease of Bacteroides dorei, Bacteroides thetaiotaomicron in w-COVID19 as compared to CTRL which are known to down-regulate ACE2 expression in the murine gut [19].

There are several limitations to this study. First, this is a pilot study conducted in a single center in urban area in central Italy with a limited number of enrolled patients. Our preliminary observations on the likely impact of SARS-CoV-2 infection on gut microbiota need to be confirmed in larger comparative trial including paucysimptomatic or asymptomatic COVID-19 patients a part from those admitted with severe or critical disease. Moreover, another limit of our work was related to our ICU patients; as reported in several papers in literature, the long stay of ICU patients can change the microbiota gut composition [32]. Definitely more patients should be included to evaluate the microbiota composition in non-COVID-19 ICU patients in comparison with w-COVID19 and i-COVID19. However, among our CTRL, three subjects were non-COVID-19 ICU patients and rectal swabs were performed one or two days after the hospitalization. Furthermore, the antibiotic treatment for ICU patients could have affected the microbiota profile, but the antibiotic therapy was administered a few days after the hospitalization in ICU, so could be assumed that this effect may not be significant. As for the link between clinical data and microbiota profiles, we are aware that age, gender and co-morbidities are factors that strongly influence the microbiota profiles, however in our study population these demographic features were not statistically significant different among the three groups.

The use of rectal swabs for gut microbiota analysis instead of standard fecal samples is not fully endorsed. Currently, feces or mucosal biopsy specimens are the biological samples most commonly used for standard 16S analysis [33]. However, as shown in several studies, stool and rectal swab are highly similar, indicating that these sampling methods could be used interchangeably to assess the community structure of the distal GI tract [34, 35]. Finally, no data on the molecular detection of SARS-CoV-2 in rectal swabs are presented in our patients. Current knowledge on whether fecal transmissibility (either orally, through fomites, or by aspiration of fecal contaminated droplets) is likely to be an important mode of COVID-19 transmission, is still limited [36]. Although out of our objective, this is an interesting research topic, particularly in health care facilities with incontinent residents.

In conclusion, significant changes in gut microbial communities in patients with COVID-19 pneumonia with different disease severity compared to CTRL have been identified. Specific microbial signatures in COVID-19 patients and roles of gut microbiota in different phase of disease and hospital setting are needed to be investigated and validated in larger cohorts.

Supporting information

S1 Table. Deseq2_phylum level_w-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s001

(XLSX)

S2 Table. Deseq2_phylum level_i-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s002

(XLSX)

S3 Table. Deseq2_phylum level_i-COVID19 vs w-COVID19 used as reference.

https://doi.org/10.1371/journal.pone.0247041.s003

(XLSX)

S4 Table. Deseq2_family level_w-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s004

(XLSX)

S5 Table. Deseq2_family level_i-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s005

(XLSX)

S6 Table. Deseq2_family level_i-COVID19 vs w-COVID19 used as reference.

https://doi.org/10.1371/journal.pone.0247041.s006

(XLSX)

S7 Table. Deseq2_genus level_w-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s007

(XLSX)

S8 Table. Deseq2_genus level_i-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s008

(XLSX)

S9 Table. Deseq2_genus level_i-COVID19 vs w-COVID19 used as reference.

https://doi.org/10.1371/journal.pone.0247041.s009

(XLSX)

S10 Table. Deseq2_species level_w-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s010

(XLSX)

S11 Table. Deseq2_species level_i-COVID19 vs CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s011

(XLSX)

S12 Table. Deseq2_specie level_i-COVID19 vs w-COVID19 used as reference.

https://doi.org/10.1371/journal.pone.0247041.s012

(XLSX)

S13 Table. Shared and distinct bacteria in w-COVID19 and i-COVID19 versus CTRL used as reference.

https://doi.org/10.1371/journal.pone.0247041.s013

(XLSX)

Acknowledgments

The authors gratefully acknowledge the Collaborators Members of the National Institute for Infectious Diseases (INMI) COVID-19 study group: Maria Alessandra Abbonizio, Amina Abdeddaim, Chiara Agrati, Fabrizio Albarello, Gioia Amadei, Alessandra Amendola, Mario Antonini, Andrea Antinori, Tommaso Ascoli Bartoli, Francesco Baldini, Raffaella Barbaro, Barbara Bartolini, Rita Bellagamba, Martina Benigni, Nazario Bevilacqua, Gianlugi Biava, Michele Bibas, Licia Bordi, Veronica Bordoni, Evangelo Boumis, Marta Branca, Donatella Busso, Marta Camici, Paolo Campioni, Maria Rosaria Capobianchi, Alessandro Capone, Cinzia Caporale, Emanuela Caraffa, Ilaria Caravella, Fabrizio Carletti, Concetta Castilletti, Adriana Cataldo, Stefano Cerilli, Carlotta Cerva, Roberta Chiappini, Pierangelo Chinello, Carmine Ciaralli, Stefania Cicalini, Francesca Colavita, Angela Corpolongo, Massimo Cristofaro, Salvatore Curiale, Alessandra D’Abramo, Cristina Dantimi, Alessia De Angelis, Giada De Angelis, Maria Grazia De Palo, Federico De Zottis, Virginia Di Bari, Rachele Di Lorenzo, Federica Di Stefano, Gianpiero D’Offizi, Davide Donno, Francesca Faraglia, Federica Ferraro, Lorena Fiorentini, Andrea Frustaci, Matteo Fusetti, Vincenzo Galati, Roberta Gagliardini, Paola Gallì, Gabriele Garotto, Saba Gebremeskel Tekle, Maria Letizia Giancola, Filippo Giansante, Emanuela Giombini, Guido Granata, Maria Cristina Greci, Elisabetta Grilli, Susanna Grisetti, Gina Gualano, Fabio Iacomi, Giuseppina Iannicelli, Giuseppe Ippolito, Eleonora Lalle, Simone Lanini, Daniele Lapa, Luciana Lepore, Raffaella Libertone, Raffaella Lionetti, Giuseppina Liuzzi, Laura Loiacono, Andrea Lucia, Franco Lufrani, Manuela Macchione, Gaetano Maffongelli, Alessandra Marani, Luisa Marchioni, Andrea Mariano, Maria Cristina Marini, Micaela Maritti, Alessandra Mastrobattista, Giulia Matusali, Valentina Mazzotta, Paola Mencarini, Silvia Meschi, Francesco Messina, Annalisa Mondi, Marzia Montalbano, Chiara Montaldo, Silvia Mosti, Silvia Murachelli, Maria Musso, Emanuele Nicastri, Pasquale Noto, Roberto Noto, Alessandra Oliva, Sandrine Ottou, Claudia Palazzolo, Emanuele Pallini, Fabrizio Palmieri, Carlo Pareo, Virgilio Passeri, Federico Pelliccioni, Antonella Petrecchia, Ada Petrone, Nicola Petrosillo, Elisa Pianura, Carmela Pinnetti, Maria Pisciotta, Silvia Pittalis, Agostina Pontarelli, Costanza Proietti, Vincenzo Puro, Paolo Migliorisi Ramazzini, Alessia Rianda, Gabriele Rinonapoli, Silvia Rosati, Martina Rueca, Alessandra Sacchi, Alessandro Sampaolesi, Francesco Sanasi, Carmen Santagata, Alessandra Scarabello, Silvana Scarcia, Vincenzo Schininà, Paola Scognamiglio, Laura Scorzolini, Giulia Stazi, Fabrizio Taglietti, Chiara Taibi, Roberto Tonnarini, Simone Topino, Francesco Vaia, Francesco Vairo, Maria Beatrice Valli, Alessandra Vergori, Laura Vincenzi, Ubaldo Visco-Comandini, Serena Vita, Pietro Vittozzi, and Mauro Zaccarelli.

References

  1. 1. Yi Ye, Philip N.P. Lagniton, Sen Ye, Enqin Li, and Ren-He Xu. COVID-19: what has been learned and to be learned about the novel coronavirus disease. Int J Biol Sci. 2020; 16(10): 1753–1766. pmid:32226295
  2. 2. Wang Wenling, Xu Yanli, Gao Ruqin, Lu Roujian, Han Kai, Wu Guizhen, et al. Detection of SARS-CoV-2 in Different Types of Clinical Specimens. JAMA 2020 Mar 11;323(18):1843–1844. pmid:32159775
  3. 3. Zhang Wei, Du Rong-Hui, Li Bei, Zheng Xiao-Shuang, Yang Xing-Lou, Hu Ben, et al. Molecular and serological investigation of 2019-nCoV infected patients: implication of multiple shedding routes. Emerg Microbes Infect. 2020;9(1):386–389. pmid:32065057
  4. 4. Chen Nanshan, Zhou Min, Dong Xuan, Qu Jieming, Gong Fengyun, Han Yang, et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet. 2020;395(10223):507–513. pmid:32007143
  5. 5. Qin Yan Gao, Ying Xuan Chen, Jing Yuan Fang. 2019 Novel Coronavirus Infection and Gastrointestinal Tract. J Dig Dis. 2020 Mar;21(3):125–126. pmid:32096611
  6. 6. Helen T Groves, Sophie L Higham, Miriam F Moffatt, Michael J Cox, John S Tregoning. Respiratory Viral Infection Alters the Gut Microbiota by Inducing Inappetence. mBio. 2020 Feb 18;11(1):e03236–19. pmid:32071269
  7. 7. Zhang Dapeng, Li Sha, Wang Ning, Tan Hor-Yue, Zhang Zhimin, Feng Yibin. The cross-talk between gut microbiota and lungs in common lung diseases. Front Microbiol. 2020 Feb 25;11:301. pmid:32158441
  8. 8. Fontana A, Panebianco C, Picchianti-Diamanti A, Laganà B, Cavalieri D, Potenza A, et al. Gut Microbiota Profiles Differ among Individuals Depending on Their Region of Origin: An Italian Pilot Study. Int J Environ Res Public Health. 2019 Oct 23;16(21):4065. pmid:31652705
  9. 9. Iebba V., Totino V., Gagliardi A., Santangelo F., Cacciotti F., Trancassini M., et al. Eubiosis and dysbiosis: The two sides of the microbiota. New Microbiol. 2016, 39, 1–12 pmid:26922981
  10. 10. Zuo T, Zhang F, Lui GCY, Yeoh YK, Li AYL, Zhan H, et al. Alterations in Gut Microbiota of Patients With COVID-19 During Time of Hospitalization. Gastroenterology. 2020 May 19:S0016-5085(20)34701-6. pmid:32442562
  11. 11. Marslnd BJ, Trompette A, Gollwitzer ES. The Gut-Lung Axis in Respiratory Disease. Ann Am Thorac Soc. 2015 Nov;12 Suppl 2:S150–6. pmid:26595731
  12. 12. Samuelson DR, Welsh DA, Shellito JE. Regulation of lung immunity and host defense by the intestinal microbiota. Front Microbiol. 2015 Oct 7;6:1085. pmid:26500629
  13. 13. Dumas Alexi, Bernard Lucie, Poquet Yannick, Geanncarlo Lugo-Villarino, Olivier Neyrolles. The Role of the Lung Microbiota and the Gut-Lung Axis in Respiratory Infectious Diseases. Cell Microbiol 2018 Dec;20(12):e12966. pmid:30329198
  14. 14. Konrad C Bradley Katja Finsterbusch, Schnepf Daniel, Crotta Stefania, Llorian Miriam, Davidson Sophia, et al. Microbiota-driven tonic interferon signals in lung stromal cells protect from influenza virus infection. Cell Rep. 2019;28(1):245–256.e4. pmid:31269444
  15. 15. W. L, Moore M. J., Vasilieva N., Sui J., Wong S. K., Berne M.A., et al., Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus, Nature 426 (2003) 450–454. pmid:14647384
  16. 16. Lu Roujian, Zhao Xiang, Li Juan, Niu Peihua, Yang Bo, Wu Honglong, et al. Genomic characterization and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet. 2020;395(10224):565–574. pmid:32007145
  17. 17. Liang Weichng, Feng Zhijie, Rao Shitao, Xiao Cuicui, Xue Xingyang, Lin Zexiao, et al. Diarrhoea May Be Underestimated: A Missing Link in 2019 Novel Coronavirus. Gut. 2020 Jun;69(6):1141–1143. pmid:32102928
  18. 18. John g, Barnaby Edward Young, Ong Sharon. COVID-19 in Gastroenterology: A Clinical Perspective. Gut. 2020 Jun;69(6):1144–1145. pmid:32198152
  19. 19. Christopher M Robinson, Julie K Pfeiffer. Viruses and the Microbiota. Annu Rev Virol. 2014;1:55–69. pmid:25821837
  20. 20. Zuo Tao, Zhang Fen, Grace C Y Lui, Yun Kit Yeoh, Amy Y L Li, Hui Zhan, et al. Alterations in Gut Microbiota of Patients With COVID-19 During Time of Hospitalization. Gastroenterology. 2020 May 19;S0016-5085(20)34701-6. pmid:32442562
  21. 21. Sharon K Kuss 1, Gavin T Best, Chris A Etheredge, Andrea J Pruijssers, Johnna M Frierson, Lora V Hooper, et al. Intestinal microbiota promote enteric virus replication and systemic pathogenesis. Science. 2011 Oct 14;334(6053):249–52. pmid:21998395
  22. 22. Zuo Tao, Liu Qin, Zhang Fen, Grace Chung-Yan Lui Eugene YK Tso, Yun Kit Yeoh, et al. Depicting SARS-CoV-2 faecal viral activity in association with gut microbiota composition in patients with COVID-19. Gut. 2020 Jul 20;gutjnl-2020-322294. pmid:32690600
  23. 23. Corona Shoenfeld Y. (COVID-19) time musings: Our involvement in COVID-19 pathogenesis, diagnosis, treatment and vaccine planning. Autoimmun Rev. 2020 Jun;19(6):102538. pmid:32268212
  24. 24. Lange Kathleen, Buerger Martin, Stallmach Andreas, Bruns Tony. Effects of Antibiotics on Gut Microbiota. Dig Dis. 2016;34(3):260–8. pmid:27028893
  25. 25. Numberger Daniela, LarsGanzert Luca Zoccarato, Kristin Mühldorfer Sascha Sauer, Hans-Peter Grossart, et al. Characterization of bacterial communities in wastewater with enhanced taxonomic resolution by full-length 16S rRNA sequencing. Sci Rep. 2019 Jul 4;9(1):9673. pmid:31273307
  26. 26. Kurtis F Budden, Shaan L Gellatly, David L A Wood, Matthew A Cooper, Mark Morrison, Philip Hugenholtz, et al. Emerging pathogenic links between microbiota and the gut-lung axis. Nat Rev Microbiol. 2017 Jan;15(1):55–63. pmid:27694885
  27. 27. Joseph M Pickard, Melody Y Zeng, Roberta Caruso, Gabriel Núñez. Gut microbiota: Role in pathogen colonization, immune responses, and inflammatory disease. Immunol Rev. 2017 Sep;279(1):70–89. pmid:28856738
  28. 28. Wang Jian, Li Fengqi, Sun Rui, Gao Xiang, Wei Haiming, Li Lan-Juan, et al. Bacterial colonization dampens influenza-mediated acute lung injury via induction of M2 alveolar macrophages. Nat Commun. 2013;4:2106. pmid:23820884
  29. 29. Nadeem O Kaakoush. Insights Into the Role of Erysipelotrichaceae in the Human Host. Front Cell Infect Microbiol. 2015 Nov 20;5:84. pmid:26636046
  30. 30. Miquel S, R Martín, Rossi O, Bermúdez-Humarán L G, Chatel J M, Sokol H, et al. Faecalibacterium prausnitzii and human intestinal health. Curr Opin Microbiol. 2013 Jun;16(3):255–61. pmid:23831042
  31. 31. Esquivel-Elizondo S, Ilhan Z E, Garcia-Peña E I, Krajmalnik-Brown R. Insights into Butyrate Production in a Controlled Fermentation System via Gene Predictions. mSystems. 2017 Jul 18;2(4):e00051–17. pmid:28761933
  32. 32. Ravi Anuradha, Fenella D Halstead Amy Bamford, Casey Anna, Thomson Nicholas M., Willem van Schaik, et al. Loss of microbial diversity and pathogen domination of the gut microbiota in critically ill patients. Microb Genom. 2019 Sep;5(9):e000293. pmid:31526447
  33. 33. Mottawea Walid, Butcher James, Li Jennifer, Abujamel Turki, Manoogian Juliana, Mack David, et al. The mucosal-luminal interface: an ideal sample to study the mucosa-associated microbiota and the intestinal microbial biogeography. Pediatr Res. 2019 May;85(6):895–903. pmid:30758325
  34. 34. Bassis CM, Moore NM, Lolans K, Seekatz AM, Weinstein RA, Young VB, et al.; CDC Prevention Epicenters Program. Comparison of stool versus rectal swab samples and storage conditions on bacterial community profiles. BMC Microbiol. 2017 Mar 31;17(1):78. pmid:28359329
  35. 35. Reyman Marta, Marlies A. van Houten, Kayleigh Arp, Sanders Elisabeth A. M. & Debby Bogaert. Rectal swabs are a reliable proxy for faecal samples in infant gut microbiota research based on 16S-rRNA sequencing. Sci Rep. 2019 Nov 5;9(1):16072. pmid:31690798
  36. 36. Amirian ES. Potential fecal transmission of SARS-CoV-2: Current evidence and implications for public health. Int J Infect Dis. 2020; 95:363–370. ). pmid:32335340