Skip to main content

Comparison of SARS-CoV-2 spike RNA sequences in feces and nasopharynx indicates intestinal replication



Little is known of possible selection and replication of SARS-CoV-2 in the intestines and if viral load in feces is associated with severity of disease. Therefore, sequence variations of the spike region in strains collected from feces and nasopharynx (NPH) from the same patients were compared. It was also investigated whether viral load in feces related to severity of COVID-19 in hospitalized patients.


SARS-CoV-2 RNA was found in 88 (79%) fecal samples from 112 patients. The complete spike region could be sequenced in 15 fecal and 14 NPH samples. Fourteen Alpha-variants and one Beta-variant of SARS-CoV-2 were identified. The majority of the viral genetic variants (viral populations) in two fecal samples, but none in NPH, had a reversion of the H69/V70 amino acid deletion normally seen in the Alpha variants. Nine fecal samples contained up to nine minority variants, each which may constitute a separate viral population. Five NPH samples had one genetic variant each, and one NPH sample contained nine minority populations of SARS-CoV-2 spike genes.


The higher genomic diversity of SARS-CoV-2 in feces compared to NPH, and the reversion of the H69/V70 deletion in Alpha variants from feces indicate a selection of viral strains and replication of SARS-CoV-2 in the gastrointestinal tract.


Individuals hospitalized due to COVID-19 usually present primarily with respiratory symptoms such as shortness of breath, oxygen desaturation and cough. However, many other symptoms are commonly found, including fever, malaise, joint pain, headache and gastrointestinal symptoms. Quantitative reverse transcription PCR (RT-qPCR) of nasopharyngeal (NPH) swabs is the most widely-used method for detection of SARS-CoV-2. A meta-analysis by Wong et al. [1] found that a significant proportion of infected patients had positive SARS-CoV-2 RNA detected in their fecal specimens. Most studies on SARS-CoV-2 genetic divergences have focused on sequencing strains from NPH samples [2,3,4,5,6,7] and few studies have sequenced the SARS-CoV-2 genomes in stool samples from patients with COVID-19 [8,9,10,11,12]. A question remains if the viral RNA detected in stool samples is from replicating viruses in the intestines or from viruses having replicated in the respiratory tract. We therefore wanted to study SARS-CoV-2 genomes by next generation sequencing in both NPH and fecal samples to further understand the relevance of viral RNA in feces. Patients with moderate to severe disease were sampled to compare the viral load and viral sequences of the spike region in fecal and NPH samples. The spike region was selected since it has been shown to have the highest number of amino acid changes and one of the highest mutation frequencies of all genomic regions of SARS-CoV-2 [5] The aim was to investigate if there was a correlation between the viral load and severity of disease of the patient, and if there were differences between viral sequences obtained from the two sampling sites, which could indicate compartmentalized replication with specific selections of viral strains depending on replication site.


Factors associated with severity of disease in patients with COVID-19

In 88 of the 112 (79%) patients SARS-CoV-2 could be detected in fecal samples.

Clinical characteristics in relation to viral load are presented in Table 1. Advanced age, systemic inflammation (high CRP-levels), and lymphocytopenia were associated with a more severe clinical outcome (data not shown). Fifty-five patients had moderate disease, 33 patients had severe disease, and 24 patients, of whom 8 were intubated, had critical disease requiring intensive care, as defined by the WHO Clinical Progression scale [13].

Table 1 Demographic factors related to the viral load of SARS-CoV-2 in feces (N = 88)

As shown in Fig. 1, the Ct values of SARS-CoV-2 RNA in fecal and NPH samples did not differ significantly between severity-groups. The Ct values for SARS-CoV-2 in NPH were, however, significantly lower than those in the fecal samples (p < 0.001), corresponding to higher viral load in the NPH. Median [IQR] days of symptoms before sampling were significantly higher for fecal samples (11 [5]) compared with NPH samples (9.5 [5]), (p < 0.001).

Fig. 1
figure 1

Relationship between cycle threshold (Ct) values and severity of disease. Ct values, with a maximum of 45 cycles, were derived from PCR analyses for SARS-CoV-2 RNA in nasopharyngeal (blue triangles) and stool (red circles) samples acquired from 112 hospitalized patients with COVID-19 from March 2021 through May 2021. Severity of disease was categorized as Moderate (N = 55, non-High Flow Nasal Oxygen (HFNO)-dependent), Severe (N = 33, HFNO-dependent), and Critical (N = 24, Intensive Care-dependent) (13) Lines represent median

Next generation sequencing of SARS-CoV-2 in feces and NPH

The qPCR Ct values for SARS-CoV-2 were, on average, 33.4 ± 3.4 (mean ± SD) in the 88 fecal samples. Sequencing was attempted in strains from 29 qPCR-positive fecal samples with a Ct value ≤ 32 (n = 29; 33%) for SARS-CoV-2. The complete spike genomic region, known to more frequently mutate than other regions of the genome [5], was obtained for the strains in 15 of these fecal samples, and the NPH sample of corresponding patients, although one NPH sample was not recovered (patient 3). For another 9 fecal samples partial spike region sequences were obtained. On average, Ct values for SARS-CoV-2 in the fecal samples with a complete spike genomic region sequence were 29.8 ± 2.1 and the Ct values of SARS-CoV-2 in the 14 corresponding NPH samples were 24 ± 5. Median [IQR] days between collection of NPH- and fecal samples was 1 [4]. The Alpha variant of SARS-CoV-2 was identified in the samples from 14 of the patients, whereas patient 9 was infected with the Beta variant. During the study period Alpha was the dominant variant (79%) in Sweden [14].

As shown in Fig. 2, SARS-CoV-2 strains from 12 of the patients had identical consensus sequences of the spike region in the NPH- and fecal samples. The amino acid H69/V70 deletion, normally seen in Alpha variant, was not present in the fecal samples from patients 8 and 12, but was present in the NPH samples from these patients. The Alpha variant infecting four patients each had one nucleotide substitution in the consensus sequence of the spike protein gene in strains both from NPH and feces, one of these patients had an Alpha variant whose nucleotide substitution was non-synonymous [15], resulting in an amino acid change. Another Alpha variant from patient 1 had two nucleotide substitutions, both non-synonymous, compared to the Alpha reference strain accession number MZ344997 (Table 2). Strains from patient 9, the only patient infected with the Beta variant, had one nucleotide substitution that was non-synonymous, compared to the Beta variant reference strain accession number MW598419 (Table 2).

Fig. 2
figure 2

Neighbor joining tree based on 3822 nucleotides of SARS-CoV-2 spike protein gene of consensus sequences from fecal samples, NPH-samples and selected sequences from Genbank. Strains are named with Genbank accession number, country of origin and year of collection. F fecal sample (red); N NPH sample (blue)

Table 2 Changes in the SARS-CoV-2 spike protein gene in the feces and the nasopharynx in patients hospitalized with COVID-19 (N = 14)

Between one to nine different SARS-CoV-2 minority populations were found in strains from nine fecal and six NPH samples (Fig. 3 and Table 2). Minority populations were found in strains from samples collected from 11 of the 14 patients.

Fig. 3
figure 3

Neighbor joining tree based on 3822 nucleotides of SARS-CoV-2 spike protein gene with consensus sequences and variants from fecal samples (red color) and NPH-samples (blue color). The tree also includes selected sequences from GenBank, designated with accession number, country of origin and year of collection. Each variant has one change that differs from the consensus sequence. F fecal sample, N NPH sample. Variants are named with V1 variant 1, V2 variant 2 etc.

For five patients, minority populations were only found in strains from feces and not in NPH, while samples from two patients’ minority populations were only found in strains from NPH. The number of minority populations of strains found in all fecal samples were larger than those found in NPH (Table 2). The different variants in the fecal samples were found in up to 39% of the sequences. Most, 82.8%, of the nucleotide substitutions of the variants were non-synonymous, resulting in 4 amino acid changes in the N-terminal domain, 4 amino acid changes in the receptor binding domain and 12 amino acid changes in the S2 region of the spike protein.

Strains from three NPH samples also had minority populations differing from the corresponding majority population only by reversion of the H69/V70 deletion. Two of the other three strains from NPH had minority populations with only one nucleotide substitution (Table 2). The strain from the sixth NPH sample (patient 14) was an outlier with nine base substitutions in the minority population. This patient also differed clinically as the only patient in this study in need of a ventilator, of the sequenced samples, and who had detectable SARS-CoV-2 RNA in NPH for approximately 1.5 months. NPH samples had an average of 1.8 base substitutions per sample with a minority population (0.4 base substitutions when patient 14 was omitted) with frequencies up to 19%, where 72.7% were non-synonymous (50% non-synonymous when patient 14 was omitted).

Except for reversions to wild type, only one nucleotide substitution, C2042A found in 5% of the sequences in the fecal minority population of patient 9, infected with the Beta variant, was shared with all the sequences from the other patients’ samples, as this substitution is normally present in the SARS-CoV-2 Alpha variant.

As can be seen in Fig. 3, the different variants from each patient sample formed clades in the phylogenetic tree. Most variants obtained from both NPH and fecal samples were similar. There are, however, exceptions, as can be seen in the clades formed by variants in samples from patients 2, 4, 8, 11 and 12, due to the H69/V70 deletion. Other than the H69/V70 deletion, all mutations are nucleotide substitutions and a minority of these substitutions are reversions to wild type SARS-CoV-2 of mutations present in the Alpha or Beta variant.


This study showed a greater viral genomic diversity of SARS-CoV-2 in fecal compared to NPH samples from hospitalized patients with COVID-19. We found that two consensus strains and one minority strain out of 13 Alpha variants of SARS-CoV-2 had reverted the Alpha specific H69/V70 deletion [16] in the fecal samples, although it was present in NPH. In addition, there was a greater variety of virus strains in feces compared to NPH. Both the reversion of H69/V70 deletion in the SARS-CoV-2 genomes in the fecal viral populations as well as the emergence of new mutations in the minority populations are possibly explained by ongoing replication in both the respiratory tract and the gastrointestinal tract causing a divergence between the two populations. So far, studies have suggested intestinal viral replication of SARS-CoV-2 [17], with virus identified in feces, known expression of ACE2 in the intestinal tract and replication of the virus in intestinal organoids. However, little data is available regarding replication of SARS-CoV-2 in human intestine. We believe that our data adds to this knowledge, and indicates that intestinal replication and furthermore local intestinal selection of SARS-CoV-2 strains is occurring.

There have been conflicting reports on whether patients with PCR-positive rectal swabs or stool samples for SARS-CoV-2 have a higher risk for a more severe clinical course of COVID-19 [18, 19]. However, two other studies [20, 21] found that there was no increased risk for severe disease in patients with PCR-positive stool samples. In the present study, stool sample positivity rate did not correlate with severity of disease. The patients with severe and critical disease were sampled at a later stage of the disease course (p < 0.001) than the patients with moderate disease, which makes it difficult to evaluate a possible correlation without taking symptom duration into account. As also previously shown, advanced age, systemic inflammation, and lymphopenia were significantly linked to a more severe clinical course [22].

The clinical utility of detection of SARS-CoV-2 RNA in stool samples from patients hospitalized with COVID-19 appears to be low, as severe disease cannot be prognosticated from stool sample results and the viral load was lower in stool samples compared to NPH samples. The detection rate could probably be increased by collecting stool samples rather than swabs, as done in this study. Typing may, however, be of importance if prolonged shedding of the virus would induce selection of different strains in feces compared to NPH.

Relative to what has been published to date, our study has strength in numbers and percentage of patients with SARS-CoV-2 RNA in feces [8,9,10,11,12]. However, this study has several limitations. The patients were not sampled with paired fecal and NPH samples on the same day and different methods were used to quantify viral load in the two different sample types, which could have had impact on the positivity rates and comparisons of viral loads in NPH and feces. Also, the duration from onset of symptoms to sampling was considerably longer for patients with severe and critical disease, which could impact positivity rates and viral load comparisons between patients with different severity of COVID-19.

Most known human fecal/orally transmitted viruses such as calici, rota, astro, adeno and enteroviruses are non-enveloped making them less susceptible to the harsh environment of the digestive system [23]. SARS-CoV-2 is an enveloped virus and should therefore not be able to remain infectious in the gastrointestinal tract. However, four different coronaviruses can infect the gastrointestinal tract in pigs [24] and the enveloped avian influenza virus reproduces in the intestines [25]. The possibility for infection in humans by enveloped viruses should therefore not be dismissed solely based on the viral structure. One possible explanation for how SARS-CoV-2 infects the gastrointestinal tract is by haematogenous spreading [17, 26]. However, ACE2, the viral receptor expressed in various different human tissues [27], is mainly expressed on the apical surfaces of the enterocytes, which suggests that the cells are exposed for infection mainly from the intestinal lumen [27,28,29]. Further, Zang et al. found that SARS-CoV-2 mainly infects the apical surface of human intestinal epithelial cells in vitro [30]. After a primary infection in the respiratory tract is established, virions may also be swallowed down together with saliva and mucus. This could act as a protective layer shielding the virions from the gastric acid as well as bile salts and pancreatic juices. This has been shown for both influenza virus A and B viruses, which when mixed with highly viscous artificial mucus remained infectious after several hours in environments simulating bile, pancreatic juices and gastric acids [31]. Amino acid alterations in receptor regions may change viral pH-resistance as is shown for avian influenza A [32] and viral tropism from intestinal to respiratory replication of the virus as is shown for porcine respiratory coronavirus [33]. The amino acids H69/V70 deletion was present in the spike protein of SARS-CoV-2 virions in NPH, as these patients were infected with the alpha variant, but in two patients this deletion had reverted in the majority of spike sequences detected in feces. The role if this specific deletion has been shown to increase the infectivity of the virus in Hek293T and Calu-3 cell lines [16] as well as be important for neutralizing antibodies [34]. Our data indicates selection of specific viruses in the intestinal tract, but the role of these amino acid residues in intestinal infection of SARS-CoV-2 requires further investigation. Also, the higher genomic diversity of SARS-CoV-2 spike gene in feces compared to NPH found in this study needs further understanding. This appears to occur regardless of the severity of the disease and may be important for intestinal replication of SARS-CoV-2. The selection of virus strains in the gastrointestinal tract during infection identified in this study may also have consequences for the spread and evolution of the virus and should be further investigated.


Next Generation Sequencing of SARS-CoV-2 RNA in fecal- and NPH samples showed a greater number of sequence variantion in feces as compared to corresponding NPH as well as a reversion to wild type of the H69/V70 deletion in two out of 13 Alpha variants in fecal samples but not in NPH samples. This can possibly be explained by ongoing replication in the intestinal- and the respiratory tract. An increased understanding of SARS-CoV-2 as an intestinal pathogen can have implications for public policy making and for finding effective cures.


Fecal samples were collected from 112 patients (74:38 [male: female]; 60 ± 14 [age ± SD] years) who were hospitalized with COVID-19 and admitted to the Department of Infectious Diseases, Sahlgrenska University Hospital, Gothenburg, Sweden from March through May 2021. The study was a sub-study of an ongoing prospective COVID-19 cohort study [35, 36]. All patients were hospitalized and had laboratory-confirmed COVID-19 determined by RT-qPCR for SARS-CoV-2 from nasopharyngeal swabs. The first day of symptoms was defined by the first day of at least one of the following symptoms prior to hospitalization: fever or chills, cough, shortness of breath or difficulty breathing, fatigue, muscle or body aches, headache, new loss of taste or smell, sore throat, congestion or runny nose, nausea or vomiting, or diarrhea [37].

Patients were divided into three groups depending on the severity of disease: moderate–hospitalized, but not high flow nasal oxygen (HFNO)-dependent; severe–HFNO-dependent; and Critical–Intensive care-dependent [13]. Clinical data and results from routine laboratory tests were retrieved from medical records and recorded in a case report form.

Feces was collected in tubes at the hospital ward and stored at 4 °C before they were transferred to the laboratory, the samples were stored at – 80 °C if not extracted upon arrival.

RNA extraction and quantification

RNA was extracted from the feces using the AllPrep PowerViral DNA/RNA Kit (Qiagen, Germany). Approximately 0.25 g of feces was added to 600 µl of the lysis buffer supplied in the kit and was homogenized at 3000 RPM in Minilys (Bertin instruments, France) for 4 min. After homogenization, the manufacturer’s protocol was followed and nucleic acids were eluted in 100 µl RNase free water. All samples were analyzed for SARS-CoV-2 RNA with qPCR. Nucleic acids in samples with cycle threshold (Ct) values ≤ 32 were also extracted with NucliSENS easyMAG (bioMérieux, France) as sequencing quality was improved when sufficient SARS-CoV-2 RNA could be extracted. Approximately 0.25 g of feces was transferred to a 10 ml tube with 2 ml PBS-buffer and glass beads. The material was then homogenized by vortexing until the solids in the feces had dissolved. The tube was centrifuged for 5 min at 3000 g, 250 µl of the supernatant was then extracted using the manufacturers off board lysis protocol. The elution volume was 110 µl.

Real-time quantitative PCR was performed using the Applied Biosystems 7300 System (ThermoFisher Scientific, Massachusetts, USA) as described previously [38] with the modification that the mastermix contained 40 U RNaseOUT (ThermoFisher Scientific, Massachusetts, USA) and one dilution (1/106) of 2 µg/ml plasmid (Eurofins genomics, Luxembourg) was used as a positive control.

Quantification of SARS-CoV-2 in nasopharyngeal samples had been determined by using either the cobas© SARS-CoV-2 test on the cobas© 6800 system (Roche Diagnostics, Switzerland) or by using the Xpert Xpress SARS-CoV-2/Flu/RSV kit on the GeneXpert IV system (Cepheid, California, USA). Nasopharyngeal samples that underwent sequencing were extracted by using the MagNA Pure 96 DNA and Viral NA Small Volume Kit on the MagNA Pure 96 System (Roche Diagnostics, Switzerland).

Next generation sequencing

SARS-CoV-2 RNA, extracted from the fecal samples with a Ct value below or equal to 32 (n = 29) were sequenced, using Ion Torrent (ThermoFisher Scientific, Massachusetts, USA) as described previously [39].

SARS-CoV-2 RNA in nasopharyngeal samples was sequenced, using either Ion Torrent [39] or Illumina (California, USA). When using Illumina, the Illumina COVIDSeq test kit was used in accordance with the protocol from the manufacturer (llumina, California, USA). Briefly, cDNA was constructed from RNA extracted with MagnaPure 96 by using random hexamer primers. Samples were amplified with COVIDSeq Primer Pool 1 HT and COVIDSeq Primer Pool 2 HT primers (llumina, California, USA). Amplified cDNA was fragmented, tagged with adaptor sequences, amplified and pooled using the automated Tecan Fluent 780 system (Tecan, Männedorf, Switzerland). The pooled libraries were quality checked, quantified and sequenced using the NextSeq 500 system with Mid Output kit V2 (300 c) or the Novaseq 6000 system with S2 reagent kit (300 c).

Assembly of consensus sequences and the determination of SARS-CoV-2 lineages in samples sequenced with Ion Torrent were made with IRMAreport plug-in (ThermoFisher Scientific, Massachusetts, USA) and assigned to variants, using Pangolin methods [39]. For strains in samples sequenced with Illumina, assembly and determination of SARS-CoV-2 lineages were made with the Arctic Networks nCoV-2019 novel coronavirus bioinformatics protocol running on a nextflow pipeline developed by Connor Lab [40]. CLC Genomic Workbench 12 (Qiagen, Hilden, Germany) was used for variant calling, minimum frequency was set to 5% and minimum depth was set to 100 reads. GenBank accession numbers for the consensus sequences of the spike region are presented in Additional file 1: Table S1.

Phylogenetic analysis

The complete spike gene, 3822 nucleotides, obtained from strains in fecal and NPH samples from 14 patients were aligned with 39 sequences of the corresponding region obtained from GenBank. The sequences from GenBank were selected to represent different strains of Alpha, Beta or Delta variants from various countries. The evolutionary distance was calculated using Maximum Composite Likelihood method by using a gamma distribution with alpha 0.03 in the program Mega6 and phylogenetic trees were constructed using neighbor-joining (NJ) method [41].

Statistical analysis

Continuous variables were compared using one-way ANOVA-test and comparisons of proportions were made using Pearson’s Chi-squared or Fischer’s exact tests as appropriate. P values < 0.05 were considered statistically significant. All statistical analyses were done with the SPSS software package version (IBM, Armonk, New York, US).

Availability of data and materials

GenBank accession numbers for the consensus sequences of the spike region are presented in Additional file 1: Table S1.


  1. Wong MC, Huang J, Lai C, Ng R, Chan FKL, Chan PKS. Detection of SARS-CoV-2 RNA in fecal specimens of patients with confirmed COVID-19: a meta-analysis. J Infect. 2020;81(2):e31–8.

    Article  CAS  Google Scholar 

  2. Czech-Sioli M, Günther T, Robitaille A, Roggenkamp H, Büttner H, Indenbirken D, et al. Integration of sequencing and epidemiological data for surveillance of SARS-CoV-2 infections in a tertiary-care hospital. Clin Infect Dis. 2022.

    Article  PubMed  Google Scholar 

  3. Faridl M, Mellyani K, Khoirunnisa K, Septiani P, Giri-Rachman EA, Nugrahapraja H, et al. RNA-seq analysis of nasopharyngeal swabs from asymptomatic and mildly symptomatic COVID-19 patients. Int J Infect Dis. 2022;122:449–60.

    Article  CAS  Google Scholar 

  4. Gaston DC, Miller HB, Fissel JA, Jacobs E, Gough E, Wu J, et al. Evaluation of metagenomic and targeted next-generation sequencing workflows for detection of respiratory pathogens from bronchoalveolar lavage fluid specimens. J Clin Microbiol. 2022.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Troyano-Hernáez P, Reinosa R, Holguín Á. Evolution of SARS-CoV-2 in Spain during the first two years of the pandemic: circulating variants, amino acid conservation, and genetic variability in structural, non-structural, and accessory proteins. Int J Mol Sci. 2022;23(12):6394.

    Article  Google Scholar 

  6. Van Elslande J, Kerckhofs F, Cuypers L, Wollants E, Potter B, Vankeerberghen A, et al. Two separate clusters of SARS-CoV-2 delta variant infections in a group of 41 students travelling from India: an illustration of the need for rigorous testing and quarantine. Viruses. 2022;14(6):1198.

    Article  Google Scholar 

  7. Wegrzyn RD, Appiah GD, Morfino R, Milford SR, Walker AT, Ernst ET, et al. Early detection of SARS-CoV-2 variants using traveler-based genomic surveillance at four US airports, September 2021-January 2022. Clin Infect Dis. 2022.

    Article  PubMed  Google Scholar 

  8. Papoutsis A, Borody T, Dolai S, Daniels J, Steinberg S, Barrows B, et al. Detection of SARS-CoV-2 from patient fecal samples by whole genome sequencing. Gut Pathog. 2021;13(1):7.

    Article  CAS  Google Scholar 

  9. Li Y, Chen YI, Ni H, Yi BO, Zhang D, Zhang J, et al. Discovery and full genome characterization of SARS-CoV-2 in stool specimen from a recovered patient. China Pol J Microbiol. 2021;70(3):401–4.

    Article  Google Scholar 

  10. Mohanty MC, Taur PD, Sawant UP, Yadav RM, Potdar V. Prolonged fecal shedding of SARS-CoV-2 in asymptomatic children with inborn errors of immunity. J Clin Immunol. 2021;41(8):1748–53.

    Article  CAS  Google Scholar 

  11. Xu Y, Kang L, Shen Z, Li X, Wu W, Ma W, et al. Dynamics of severe acute respiratory syndrome coronavirus 2 genome variants in the feces during convalescence. J Genet Genom. 2020;47(10):610–7.

    Article  Google Scholar 

  12. Zhang Y, Chen C, Song Y, Zhu S, Wang D, Zhang H, et al. Excretion of SARS-CoV-2 through faecal specimens. Emerg Microbes Infect. 2020;9(1):2501–8.

    Article  CAS  Google Scholar 

  13. Characterisation WHOWGotC. Management of CI. A minimal common outcome measure set for COVID-19 clinical research. Lancet Infect Dis. 2020;20(8):e192–7.

    Article  Google Scholar 

  14. Folkhälsomyndigheten. Statistik om SARS-CoV-2 virusvarianter av särskild betydelse 2022. Last updated 29 July 2022. Accessed 10 Apr 2022.

  15. Ohta T. Synonymous and nonsynonymous substitutions in mammalian genes and the nearly neutral theory. J Mol Evol. 1995;40(1):56–63.

    Article  CAS  Google Scholar 

  16. Meng B, Kemp SA, Papa G, Datir R, Ferreira I, Marelli S, et al. Recurrent emergence of SARS-CoV-2 spike deletion H69/V70 and its role in the alpha variant B117. Cell Rep. 2021;35(13):109292.

    Article  CAS  Google Scholar 

  17. Guo M, Tao W, Flavell RA, Zhu S. Potential intestinal infection and faecal-oral transmission of SARS-CoV-2. Nat Rev Gastroenterol Hepatol. 2021;18(4):269–83.

    Article  CAS  Google Scholar 

  18. Lin W, Xie Z, Li Y, Li L, Wen C, Cao Y, et al. Association between detectable SARS-COV-2 RNA in anal swabs and disease severity in patients with coronavirus disease 2019. J Med Virol. 2021;93(2):794–802.

    Article  CAS  Google Scholar 

  19. Chen GQ, Luo WT, Zhao CH, Li CN, Hong ZS, Xia JY, et al. Comparison of clinical characteristics between fecal/perianal swab nucleic acid-positive and -negative patients with COVID-19. J Infect Dev Ctries. 2020;14(8):847–52.

    Article  CAS  Google Scholar 

  20. Chen Y, Chen L, Deng Q, Zhang G, Wu K, Ni L, et al. The presence of SARS-CoV-2 RNA in the feces of COVID-19 patients. J Med Virol. 2020;92(7):833–40.

    Article  CAS  Google Scholar 

  21. Zheng S, Fan J, Yu F, Feng B, Lou B, Zou Q, et al. Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China, January-March 2020: retrospective cohort study. BMJ. 2020;369:m1443.

    Article  Google Scholar 

  22. Chen G, Wu D, Guo W, Cao Y, Huang D, Wang H, et al. Clinical and immunological features of severe and moderate coronavirus disease 2019. J Clin Invest. 2020;130(5):2620–9.

    Article  CAS  Google Scholar 

  23. Bushman FD, McCormick K, Sherrill-Mix S. Virus structures constrain transmission modes. Nat Microbiol. 2019;4(11):1778–80.

    Article  CAS  Google Scholar 

  24. Liu Q, Wang HY. Porcine enteric coronaviruses: an updated overview of the pathogenesis, prevalence, and diagnosis. Vet Res Commun. 2021;45(2–3):75–86.

    Article  Google Scholar 

  25. Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y. Evolution and ecology of influenza a viruses. Microbiol Rev. 1992;56(1):152–79.

    Article  CAS  Google Scholar 

  26. Lin L, Lu L, Cao W, Li T. Hypothesis for potential pathogenesis of SARS-CoV-2 infection-a review of immune changes in patients with viral pneumonia. Emerg Microbes Infect. 2020;9(1):727–32.

    Article  CAS  Google Scholar 

  27. Hikmet F, Mear L, Edvinsson A, Micke P, Uhlen M, Lindskog C. The protein expression profile of ACE2 in human tissues. Mol Syst Biol. 2020;16(7):e9610.

    Article  CAS  Google Scholar 

  28. Hamming I, Timens W, Bulthuis ML, Lely AT, Navis G, van Goor H. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J Pathol. 2004;203(2):631–7.

    Article  CAS  Google Scholar 

  29. Suarez-Farinas M, Tokuyama M, Wei G, Huang R, Livanos A, Jha D, et al. Intestinal inflammation modulates the expression of ACE2 and TMPRSS2 and potentially overlaps with the pathogenesis of SARS-CoV-2-related disease. Gastroenterology. 2021;160(1):287-301 e20.

    Article  CAS  Google Scholar 

  30. Zang R, Gomez Castro MF, McCune BT, Zeng Q, Rothlauf PW, Sonnek NM, et al. TMPRSS2 and TMPRSS4 promote SARS-CoV-2 infection of human small intestinal enterocytes. Sci Immunol. 2020;5(47):eabc3582.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Hirose R, Nakaya T, Naito Y, Daidoji T, Watanabe Y, Yasuda H, et al. Mechanism of human influenza virus RNA persistence and virion survival in feces: mucus protects virions from acid and digestive juices. J Infect Dis. 2017;216(1):105–9.

    Article  CAS  Google Scholar 

  32. Fujimoto Y, Ito H, Ono E, Kawaoka Y, Ito T. The low-pH resistance of neuraminidase is essential for the replication of influenza a virus in duck intestine following infection via the oral route. J Virol. 2016;90(8):4127–32.

    Article  CAS  Google Scholar 

  33. Krempl C, Schultze B, Laude H, Herrler G. Point mutations in the S protein connect the sialic acid binding activity with the enteropathogenicity of transmissible gastroenteritis coronavirus. J Virol. 1997;71(4):3285–7.

    Article  CAS  Google Scholar 

  34. Xie X, Liu Y, Liu J, Zhang X, Zou J, Fontes-Garfias CR, et al. Neutralization of SARS-CoV-2 spike 69/70 deletion, E484K and N501Y variants by BNT162b2 vaccine-elicited sera. Nat Med. 2021;27(4):620–1.

    Article  CAS  Google Scholar 

  35. Marklund E, Leach S, Axelsson H, Nystrom K, Norder H, Bemark M, et al. Serum-IgG responses to SARS-CoV-2 after mild and severe COVID-19 infection and analysis of IgG non-responders. PLoS ONE. 2020;15(10):e0241104.

    Article  CAS  Google Scholar 

  36. Yilmaz A, Marklund E, Andersson M, Nilsson S, Andersson LM, Lindh M, et al. Upper respiratory tract levels of severe acute respiratory syndrome coronavirus 2 RNA and duration of viral RNA shedding do not differ between patients with mild and severe/critical coronavirus disease 2019. J Infect Dis. 2021;223(1):15–8.

    Article  CAS  Google Scholar 

  37. Prevention CfDCa. Symptoms of Covid-19 CDC; 2021. Updated 220322. []

  38. Saguti F, Magnil E, Enache L, Churqui MP, Johansson A, Lumley D, et al. Surveillance of wastewater revealed peaks of SARS-CoV-2 preceding those of hospitalized patients with COVID-19. Water Res. 2021;189:116620.

    Article  CAS  Google Scholar 

  39. Ringlander J, Olausson J, Nystrom K, Harnqvist T, Jakobsson HE, Lindh M. Recurrent and persistent infection with SARS-CoV-2—epidemiological data and case reports from Western Sweden, 2020. Infect Dis. 2021;53(12):900–7.

    Article  CAS  Google Scholar 

  40. ConnorLab. ncov2019-artic-nf: GitHub repository; 2021. []

  41. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  Google Scholar 

Download references


Open access funding provided by University of Gothenburg. This study was supported by grants from the Swedish research council (2021-06545), Swedish Research Council for Sustainable Development, Formas (2020-02710) and the Sci-lifelab National COVID-19 Research Program, financed by the Knut and Alice Wallenberg Foundation (2020-0182 and 2020.0241) and by the Swedish Government, under an agreement between the Swedish government and the county councils (ALF agreement ALFGBG-965885 and 723,471).

Author information

Authors and Affiliations



MG, LA, HN and KN conceptualized and designed the study. TB conducted sampling and clinical data collection. AK and HW analysed the samples. AK, KN and HN interpreted sequenced data. TB and AK wrote the original draft. HN and KN substantially revised the manuscript. TB, AK, KN, HW, MG, LA and HN reviewed the manuscript. MG, HN and LA acquired funding. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kristina Nyström.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Swedish Ethical Review Authority (Dnr: 2020-01771 and 2021-01322) and patients were included after informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. GenBank accession numbers for the consensus sequences of the fecal- and nasopharyngeal samples.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Beck-Friis, T., Kärmander, A., Nyström, K. et al. Comparison of SARS-CoV-2 spike RNA sequences in feces and nasopharynx indicates intestinal replication. Gut Pathog 14, 35 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • SARS-CoV-2
  • COVID-19
  • Gastrointestinal tract
  • Feces
  • Genomic structural variation
  • Mutation