Keywords
SARS-CoV-2, COVID-19, gut microbiota
This article is included in the Pathogens gateway.
This article is included in the Emerging Diseases and Outbreaks gateway.
This article is included in the Cell & Molecular Biology gateway.
This article is included in the Coronavirus collection.
SARS-CoV-2, COVID-19, gut microbiota
Recent scientific articles and reviews1–3 discuss the relationship between gastrointestinal microbiota and COVID-19 disease. In particular, the prolonged presence of SARS-CoV-2 RNA in human faecal samples from COVID-19 patients has been reported4 and the potential role of orofecal transmission of SARS-CoV-2 has been examined in systematic reviews5,6 and open evidence briefs.7,8 Cases of SARS-CoV-2 detection in faecal samples from patient with typical symptoms of COVID-19 but negative to multiple SARS-CoV-2 real-time reverse transcription polymerase chain reaction (rRT-PCR) tests on oropharyngeal and nasopharyngeal swabs have been reported.9 SARS-CoV-2 faecal viral activity was depicted in association with gut microbiota composition in patients with COVID-19,10 and the replicating virus was detected in faeces.11 At the same time, Wölfel et al.12 observed high viral RNA concentration in stool samples, but reported isolation of infectious virus only from throat- and lung-derived samples, while Yao et al.13 had indication of viable SARS-CoV-2 particles in stool samples, denoting that the detailed biology of SARS-CoV-2 is not yet fully elucidated. Our experiments further explored the relationship between COVID-19 disease and SARS-CoV-2 infected faeces to provide data relevant for pandemic understanding and disease management. The results however did not correspond with current thinking of the epidemiology of SARS-CoV-2 and, therefore, we believe a quick sharing with the scientific community of our findings is imperative.
The experimental design included:
1) the inoculation of NutriSelect™ Plus nutrient broth at 37°C, fit for the growth of more fastidious bacteria, with a faecal sample (stool) from one subject positive to SARS-CoV-2 and from a healthy individual (here called sample A and sample B, respectively) following written informed consent.
2) The assessment of the presence of SARS-CoV-2 RNA in both samples, after seven days of culture, using the Luminex technology; with confirmation of the presence of SARS-CoV-2 RNA in sample A, and of its absence in sample B.
3) Inoculation of sample B with the supernatant of sample A, obtained after centrifugation (hereafter called sample B(A+)) and resuspension of the formed pellet (sample C).
4) Incubation of all samples (A, B, B(A+) and C) for 30 days under the same conditions in NutriSelect™ Plus nutrient broth at 37°C with measurement of the viral RNA load in each sample at days 1, 2, 3, 7, 14, 21, and 30, following the date of inoculation (day 0).
5) Antibiotic treatment on 18 aliquots derived from sample B(A+) at day 21, consisting in the addition of a specific antibiotic (each of the following: metronidazole, clindamycin, lincomycin, piperacillin+tazobactam, vancomycin, amoxicillin, ampicillin, cefixime, ceftriaxone, meropenem, rifaximin, azithromycin, erythromycin, gentamicin, ciprofloxacin, colistin, levofloxacin, and teicoplanin) to each aliquot. SARS-CoV-2 RNA load was measured by Luminex technology in each aliquot before and 3 days after antibiotic administration.
6) In parallel, additional analyses were performed to evaluate and monitor over time the bacterial growth and metabolic activity of all samples and all aliquots of sample B(A+), using SANIST Biotyper according to the method described by Cristoni et al.14
7) Purification and analysis of the peptides present in sample B(A+) at day 30.
The details of the procedures and protocols used are presented in the extended data, together with a schematic representation of the experimental design (Graphical abstract).35
The experimental design included a series of analyses (performed on all samples A, B, B(A+) and C) aimed at verifying: 1) the permanence/survival over time and the eventual multiplication of SARS-CoV-2 RNA in vitro; 2) the presence/synthesis of SARS-CoV-2 peptides in the cultures having confirmed SARS-CoV-2 RNA presence; 3) the effect of antibiotics administration in sample B(A+); 4) the concomitant presence of other metabolites; and 5) the characterisation of the bacterial samples, including the verification of the presence of eukaryotic cells.
The results presented and discussed here, carried out over a period of 30 days, confirmed the extra-corporal multiplication of SARS-CoV-2 RNA: viral load highly increased over time in sample B(A+), slightly increased in sample A, decreased in sample C while, as expected, sample B was found constantly negative (Figure 1).
In order to verify the reproducibility of our results, the whole experiment was repeated independently three times using the same infected and healthy samples (with the exception that the repetition experiments were stopped at day 14 instead of day 30). The results of the SARS-CoV-2 RNA load measurements in the repetitions are reported in Figure 2, where results are depicted as the average of the measurements of the three repetitions, together with the calculated standard deviations. The trend was confirmed, with the increase over time of measured viral RNA load in sample A and sample B(A+). Decrease in sample C and no detection in sample B were also confirmed, but they are not reported in Figure 2.
In addition, three new couples of faecal samples from different “infected donors” (i.e. sources of A) and “healthy recipients” (i.e. sources of B) have been recruited, and subject to the same experimental procedure. Samples were collected from anonymous donors, and no information (i.e. age, sex, blood serotype, severity of the disease, time of the collection, fatality, etc.) is available. All combinations of “infected donors” sources (A1, A2 and A3) and “healthy recipients” donors (B1, B2 and B3) were subject to the same experimental procedure. Although with certain differences, the observed trends are similar, confirming the increase over time of SARS-CoV-2 RNA load in samples of type A and in samples of type B(A+), independently from the sources of A and B.
Aliquots of sample B(A+) tested after three days of culture in the presence of the single different antibiotics belonging to different classes were analysed and the SARS-CoV-2 RNA load measured in each of them. SARS-CoV-2 RNA load was found to be influenced by the presence of antibiotics in different ways (Figure 4):
• SARS-CoV-2 RNA load was reduced to undetectable levels in the four aliquots treated with metronidazole, vancomycin, amoxicillin and azithromycin, respectively.
• SARS-CoV-2 RNA load decreased by 20% to 85% in the aliquots treated with piperallicin+tazobactam, ampicillin, cefixime, ceftriaxone, meropenem, gentamicin, ciprofloxacin and teicoplanin. For example, cefixime induced a decrease of viral RNA load of 85%, ciprofloxacin of 61% and teicoplanin of 56%.
• SARS-CoV-2 RNA load was not substantially affected by the presence of clindamycin, lincomycin, rifaximin, erythromycin, colistin and levofloxacin.
After 30 days of bacterial growth in culture, aliquots of samples B(A+) (from couple 0) were collected and tested for the presence of SARS-CoV-2-related peptides using mass spectrometry (details are described in extended data35). Several peptides found in the aliquot from sample B(A+) were assigned to SARS-CoV-2 proteins. As shown in Table 1, the sequence of some of the peptides (pep51 and pep121, matching on NSP3; pep199, matching on the spike protein; pep25 and pep68, matching on NS3 and N, respectively) have one or more amino acid (AA) changes (highlighted in red) with respect to the translations of CDS regions reported in the reference ‘Severe acute respiratory syndrome coronavirus 2 isolate Wuhan-Hu-1, complete genome sequence’ (GenBank LOCUS: NC_045512.2). No SARS-CoV-2-related peptide was identified in the aliquot of sample B.
The identified AA changes have been checked for their existence among the observed variations in SARS-CoV-2 sequenced isolates available in GISAID15 at time of writing. As shown in Table 2, all of them except NSP3:A274K in pep51 have been already reported in humans; the majority of them have been never reported in the country of origin of the samples (Italy), the remaining ones have been observed in samples sequenced in Italy, but after the time of collection of the infected sample A (February 2020). Some of the found peptides mapping on the SARS-CoV-2 spike protein are shown in Figure 5.
We have already described the detection of toxin-like peptides in plasma, urine and faecal samples from COVID-19 affected individuals (Cristoni et al.,16 under review). The evaluation on the potential release of toxic-like peptides in aliquots from sample B(A+) has been assessed by performing the same analyses. Toxin-like peptides have been observed, but their presence was completely reduced to negligible levels in the aliquots treated with metronidazole and vancomycin administration (data not shown). These results need to be carefully interpreted, taking into account the different antimicrobials kinetics.
Samples A and B(A+) were found to contain some bacterial genera particularly abundant and metabolically active during the whole experiment, as shown in Figure 6.
Aliquots of samples A, B, and B(A+), collected at different times, were analysed with both transmission electron microscope (TEM) and scanning electron microscope (SEM) to verify the presence of eukaryotic cells. More than 30 different preparations (including at 30 days of culture) have been observed: none was found to contain any structure resembling cells with nuclei, but only bacterial cells. Analyses of images on samples A and B(A+) revealed the presence of virus-like particles interacting with bacterial cells. Immune electron microscopy is ongoing for confirming that these particles are of SARS-CoV-2 origin (in preparation).
Our results indicate that the SARS-CoV-2 genome, in addition to its known interactions with eukaryotic cells, is additionally capable of replicating outside the human body, suggesting a possible ‘bacteriophage-like’ mode of action. It is not clear yet whether the SARS-CoV-2 genome could just be replicated by its RNA polymerase (which would correspond to a bacteriophage pseudo-lysogenic mechanism), or if the production of full-blown SARS-CoV-2 replicating particles within the bacteria occur (which would correspond to the typical lytic cycle of bacteriophages). In either case, according to our knowledge, this is a novel observation and has never described before for SARS-CoV-2.
The experiment here described was repeated three additional times using the same samples. In addition, independent replications were performed in a 3×3 design using different starting material. In all of them, very similar trends were observed and the increase of SARS-CoV-2 RNA load in sample A and sample B(A+) was confirmed in all experiments.
In the cases of replications with different starting material (i.e. faecal samples from different “infected donors” as sources of A, and “healthy recipients” as sources of B), the trends are similar and confirm that the experiment is reproducible. A whole-genome metagenomic sequencing on the samples is ongoing, aimed at characterising further which bacterial species are candidate target(s) of the observed behaviour of SARS-CoV-2.
Whereas our experimental design was intended to grow bacterial cells, the possibility that SARS-CoV-2 RNA increase could be due to replication in human cells present in the original faecal samples, was considered. The human cells most abundantly present in faecal samples are colonic epithelial cells (colonocytes). Loktionov,16 reported that cell exfoliation events from colonic epithelium are rare under normal conditions, while they dramatically increase in cases of uncontrolled growth of cells not under physiologic control (like in neoplasia), when cell removal by apoptosis does not function properly. In addition, Iyengar et al.17 reported that colonic epithelial cells terminally differentiated are devoid of proliferative activity. More recently, Nair et al.18 and Chandel et al.19 developed specific methodologies to recover viable colonocytes from stool. In our case, both sample A and B originated from adult individuals with no diagnosis of cancer. In addition, it is unlikely that human cells potentially present in samples A and B are able to:
• grow in a culture medium typically formulated for bacteria and not containing growth factors, serum, nor other important components for eukaryotic cell sustainment;
• survive in such a medium for 30 days, and in co-occurrence with an event of SARS-CoV-2 infection;
• multiply in the absence of specific CO2 concentration conditions.
Also, the possibility of interaction between SARS-CoV-2 and other eukaryotic organisms potentially present in the cultures, as e.g. parasitic nematodes and fungal cells, has been considered.
During the whole experiment, parasitic nematodes were not noted at visual inspections by eye. In addition, stool of sample B was independently analysed and certified to be free of known parasites and microbial pathogens (certification provided by the Italian diagnostic laboratory Biomolecular Diagnostic Srl). Parasitic nematodes are usually not able to survive outside the host and many intestinal roundworms (like those of genus Ascaris) release antimicrobial factors that interfere with bacterial growth,20 in contrast with the found high increase of metabolic activities of some bacterial genera. In the used medium, chemical elements relevant for (parasitic and not) nematodes (e.g. cholesterol and traces of metals) are missing.
The possibility of involvement of the mycobiome fraction present in the stool was considered. As highlighted by Chin et al.,21 multifaceted and multidisciplinary approaches are necessary to identify uncultivatable, low-abundance, permanent and transient fungal species residing in the gut, confirming that the human mycobiome is not yet fully characterised. Accordingly, while the ability of unknown fungi to grow in the used culture medium cannot be excluded, no significant metabolic activity of Candida albicans, most commonly found in the microbiota, was observed.
Finally, inspections of images from TEM and SEM on more than 30 different preparations did not reveal presence of eukaryotic cells (in preparation). If on the one hand, the possibility that a nematode or another unknown eukaryotic cell is able grow in the medium cannot be excluded, the used conditions make this possibility very unlikely. Anyhow, the ability of SARS-CoV-2 to interact either with nematodes or with fungal cells has never been observed before and would be a novel and surprising observation as well.
As indicated above, several peptides matching to SARS-CoV-2 proteins were found in the aliquot from sample B(A+). The identification of peptides with amino acid changes, compared to the translations of CDS regions of the reference SARS-CoV-2 genome, is intriguing but is compatible with the mechanism of viral replication in bacteria. RNA viruses such as SARS-CoV-2 inhabit the host as a population of variants called quasispecies, i.e. a group made of different variants that are genetically linked through mutation events, and contribute collectively to the characteristics of the whole (viral) population in the host.22 Recent studies highlighted the significant amount of intra-host genomic diversity in SARS-CoV-2 samples.23,24 In a ‘bacteriophage-like’ mode of action, as bacteria were grown for 30 days, it can’t be excluded that the observed amino acid changes represent viral quasispecies emerged through replication events in bacterial hosts. In relation to this, recent studies25–27 evidenced hypermutations occurrences in SARS-CoV-2 genomes and suggested APOBEC and ADAR deaminases as the possible responsible of these phenomena. The APOBEC family is related to bacterial, yeast, and plant deaminases all possessing highly conserved amino acid motifs responsible for coordination of zinc in the active site.28 As no sequencing was performed on the original infected stool sample, the presence of SARS-CoV-2 haplotypes in the initial SARS-CoV-2 population used to perform the experiments, therefore justifying the amino acid changes observed, cannot be excluded. However, all the amino acid changes found have been reported in sequences of SARS-CoV-2 found for the first time after the date of collection of the infected sample A (February 2020), and some of them have never been reported in Italy, the country of origin of the samples.
On the other hand, other mechanisms like those at the basis of diversity-generating retroelements (DGR) systems29 have that could contribute to SARS-CoV-2 hypermutation phenomena have recently been described in bacteria, and could therefore be responsible of the AA changes found.
These results can potentially provide new insights in the epidemiology of SARS-CoV-2. Considering the possible impact and implications that such relationship has on the manifestation, therapy and control of COVID-19 disease, some questions immediately arise like e.g.:
• Can this ‘bacteriophage-like’ behaviour of SARS-CoV-2 explain the long-term presence of SARS-CoV-2 observed in some recovered patients?30
• Can antibiotics and/or bacteriophage-based therapies play a role in the treatment of COVID-19 affected patients?31
• How would the (antecedent) administration of antibiotics to patients, influencing the microbiota population, impact the clinical course of the disease?32
• Can the involvement of bacteria in COVID-19 epidemiology help to explain clinical observations, like the elevated serum C-reactive protein, procalcitonin, D-dimer, and ferritin associated with poor outcomes in COVID-19?33
These questions are only examples of the plethora of questions to be addressed. Our results support the way to tackle COVID-19 pandemic proposed by Mushi,34 i.e. by using the One Health holistic approach. If individuals are considered not only human bodies, but as ‘holobionts’, i.e. discrete ecological units that need to be studied and treated as such, a deeper understanding of the role of the microbial community living in the human body is fundamental to tackle COVID-19 disease.
Faecal samples were collected and handled by CranioMed S.R.L. from anonymous donors who agreed to participate in this study by signing informed consent, as foreseen by Italian legislation. No personal information (i.e. age, sex, blood serotype, severity of the disease, time of the collection, fatality, etc.) were collected.
The study is compliant with the JRC Scientific Integrity and Research Ethics guidance.
The scientific output expressed does not imply a policy position of the European Commission. Neither the European Commission nor any person acting on behalf of the Commission is responsible for the use that might be made of this publication.
Zenodo: Underlying data for ‘Increase of SARS-CoV-2 RNA load in faecal samples prompts for rethinking of SARS-CoV-2 biology and COVID-19 epidemiology’. https://doi.org/10.5281/zenodo.4723550.35
The project contains the following underlying data:
Mass spectrometry raw data of the peptides.
NCBI Genome: Severe acute respiratory syndrome coronavirus 2 isolate Wuhan-Hu-1, complete genome. Accession number: NC_045512.2; https://www.ncbi.nlm.nih.gov/nuccore/NC_045512.2
NCBI Protein: Surface glycoprotein [Severe acute respiratory syndrome coronavirus 2]. Accession number: YP_009724390.1; https://www.ncbi.nlm.nih.gov/protein/1796318598
Zenodo: Extended data for ‘Increase of SARS-CoV-2 RNA load in faecal samples prompts for rethinking of SARS-CoV-2 biology and COVID-19 epidemiology’. https://doi.org/10.5281/zenodo.4723550.35
Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).
We gratefully acknowledge the Authors and the Originating laboratories where the clinical specimen or virus isolate was first obtained and the Submitting laboratories, where sequence data have been generated and submitted to GISAID.
The efforts and commitment put into this research work are dedicated to all EU citizens who suddenly left us because of COVID-19, and to their families.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Partly
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Microbiota, Probiotics, taxonomy, Microbiology, Molecular biology, Biotechnology.
Is the work clearly and accurately presented and does it cite the current literature?
Partly
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
References
1. Schenck LP, Surette MG, Bowdish DM: Composition and immunological significance of the upper respiratory tract microbiota.FEBS Lett. 2016; 590 (21): 3705-3720 PubMed Abstract | Publisher Full TextCompeting Interests: No competing interests were disclosed.
Reviewer Expertise: Virology and immunology
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 3 (revision) 01 Jul 21 |
read | |
Version 2 (revision) 23 Jun 21 |
read | read |
Version 1 11 May 21 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)