Limited intestinal inflammation despite diarrhea, fecal viral RNA and SARS-CoV-2-specific IgA in patients with acute COVID-19

Patient demographics and clinical characteristics

We collected stool samples from 44 symptomatic COVID-19 patients hospitalized in New York City between April 15 and May 21 2020 (Table S1 and Date file S1). Samples were collected longitudinally during acute and convalescent phases of illness. The median time from symptom onset to sample collection of first stool sample was 16 days (range 2-66) and of second stool sample was 24.5 days (11-72) (Table S1 and Date file S2).

The majority of patients had at least one comorbidity (Table 1) and fifty percent of patients (n=22; 50%) presented to the hospital with moderate disease (see Methods for definitions of disease severity). During hospitalization 10 of the 22 (45.5%) patients who presented with moderate disease developed severe COVID-19 disease, while none of the patients who presented with mild symptoms had clinical worsening during hospitalization (Table 1 and Date file S1).

Seven (15.9%) patients died and 13 patients (29.5%) were admitted to the intensive care unit (ICU) during their admission, with 14 patients (31.8%) meeting the composite outcome of ICU admission or death (Table 2). The majority of patients received antibiotics (28 (63.6%)) and therapeutic anticoagulation (31 (70.5%)) during their admission (Table S2 and Date file S1).

Given the aim of the study we prioritized enrollment of patients with GI symptoms. Thirty-one (70.5%) of patients in the cohort presented with GI symptoms (Figure 1A). The most reported GI symptom was diarrhea (59.1%), followed by nausea (34.1%) and then vomiting (15.9%) (Figure 1B). Patients with GI symptoms were significantly younger than patients without GI symptoms (mean (SD) age was 52.9 (14.6) vs 63.1 (14.2) respectively; p=0.04) (Table 1). Distribution of race and ethnicity and of comorbidities was similar in patients who presented with or without GI symptoms (Table 1). Ten (32.2%) of patients with and 6 (46.2%) without GI symptoms presented with severe disease. Overall distribution of disease severity on admission was similar in patients presenting with or without GI symptoms (p=0.60) as was peak severity during admission (p=0.75) (Table 1). Similar rates of mortality were observed in patients presenting with GI symptoms (n=5; 16.1%) compared to those without GI symptoms (n=2; 15.4%) (p > 0.99); however, significantly lower proportion of patients with GI symptoms (n=6; 19.4%) were admitted to the ICU compared with those without GI symptoms (n=7; 53.8%%, p=0.03; Fisher’s exact test) (Table 2). As expected, more severe disease at presentation (p=0.09) and more severe disease during admission (p=0.02), were associated with higher frequencies of ICU admission or death (Table S3).

• Download figure
• Open in new tab

Fig. 1. Gastrointestinal (GI) symptoms and serologic parameters in hospitalized COVID-19 patients.

(A) Proportion of patients with GI manifestations in this cohort of hospitalized COVID-19 patients in whom stool samples obtained. (B) Serum concentrations of IL-1β, IL-6, IL-8, and TNFα and select laboratory values at the time of admission in patients with and without GI symptoms and (C) in patients with severe compared with non-severe COVID-19 disease. Each point represents an individual value for a patient, the box plot represents the median and the interquartile range and the p-values calculated using the Mann-Whitney test with significance defined as p<0.05.

When looking at inflammatory markers in serum, there were no significant differences observed between those with GI symptoms and those without, although those with GI symptoms tended to have lower serum IL-6 and IL-8 (Fig. 1B). No notable differences in serum inflammatory markers or laboratory values were noted in patients with GI symptoms or with severe disease (Fig. 1B, S1A and B).

Peak disease severity was associated with significantly higher serum IL-6, IL-8, and TNF-α consistent with prior reports [19, 20, 23] (Fig. 1C). Additionally, severe disease was associated with higher serum C-reactive protein (CRP), Lactate dehydrogenase (LDH), D-dimer, ferritin and procalcitonin (Fig. 1C and Fig. S1B). There were no significant differences in liver enzymes in those with severe disease compared to those with mild or moderatedisease (Fig. S1B).

SARS-CoV-2 RNA is frequently detected in stool and higher viral loads are associated with diarrhea and death

To detect the presence of SARS-CoV-2 RNA in stool, we performed real-time quantitative PCR (qPCR) targeting three regions of the SARS-CoV-2 genome (see Methods for details). A positive result was defined as two out of three SARS-CoV-2 primers with cycle thresholds (Ct) <40. Eighteen of the 44 patients (40.9%) had/showed a positive fecal PCR over this threshold during admission (21 of 62 samples tested (33.9%)).

To investigate the length of time SARS-CoV-2 RNA may be found in the stool, we stratified our qPCR results by the time from symptom onset. SARS-CoV-2 RNA was unlikely to be detected in stool samples more than 28 days after symptom onset (p = 0.052; Fig. 2A-B and S2A). Only 2/17 (11.8 %) tests were positive for viral RNA in the stool after 28 days from symptom onset as opposed to 19/45 (42.2%) of tests before this point (p=0.03, Fisher exact test; Fig. 2B and S2A). However, within the initial 28-day period there was no difference as to when stool samples tested positive for viral RNA (p = 0.24; Fig. 2A and S2A). Interestingly, we never detected viral genome in stool collected after a patient had a recorded negative nasopharyngeal PCR result (n=7 stool samples from 5 different patients collected after a negative nasal PCR; median 11 days between negative nasal PCR and stool sample collection in this subset).

• Download figure
• Open in new tab

Fig. 2. Detection and sequencing of SARS-CoV-2 virus genome in stool.

(A). The median fecal SARS-CoV-2 Ct of each sample, plotted according to the time of collection relative to either the onset of symptoms or the first positive nasal PCR. P value - Kruskal-Wallis. (B). The median fecal SARS-CoV-2 Ct of each sample, plotted according to the relative time of collection: day 28 as a cut off. Relative to either the onset of symptoms or the first positive nasal PCR. P value - Mann-Whitney. (C). The lowest fecal SARS-CoV-2 Ct measured for each donor shows higher peak viral loads in donors who did not survive. P value - Mann-Whitney. (D). The lowest fecal SARS-CoV-2 Ct value measured from each donor during admission, plotted according to the peak severity of COVID19 symptoms- P value - Kruskal-Wallis. (E). Fecal SARS-CoV-2 Ct values in samples from donors with and without diarrhea. P value - Mann-Whitney. (F). The relationship between the fecal SARS-CoV-2 Ct values and the Bristol Score of each sample. Spearman rank test. (G) Phylogenetic relationships of the SARS-CoV-2 isolates from stool (Stool, dark red) and other nasopharyngeal virus isolates from the Mount Sinai Health System (MSHS) in NYC (Nasopharyngeal, salmon), in a background of global strains (grey). Main phylogenetic clades are indicated. The scale bar represents divergence from the root (reference genome NC_045512) as number of single nucleotide variants (SNVs). The panel in the bottom right indicates the clade distribution of the MSHS nasopharyngeal samples from 403 SARS-CoV-2 genomes from February 29 to July 15. (H) Sequence alignment depicting clade and cluster defining (teal), and other (grey) SNVs present in the stool SARS-CoV-2 isolates from stool (S) or nasopharynx (N) of the indicated donors. The clade and amino acid substitutions are indicated when applicable.

To determine whether fecal viral genome loads were associated with COVID-19 symptoms or outcomes, we stratified our results and found that patients who did not survive had significantly higher fecal SARS-CoV-2 genome loads than survivors (p = 0.021; Fig. 2C). However, there was no association between the fecal Ct values and disease severity (p = 0.63; Fig. 2D), and it was notable that 3/6 (50%) of patients with mild COVID-19 disease had detectable virus in stool (Fig. 2D). We observed higher viral genome loads in samples from patients reporting diarrhea (p 0.049; Fig. 2E and S2B) and a significant correlation between viral genome load and Bristol stool scores of the samples (Spearman ρ = −0.29, p = 0.025; Fig. 2F).

We sequenced the genome of SARS-CoV-2 genomes from stool samples of three patients using short read sequencing [24]. Each of these isolates belonged to clade 20C, the clade also most commonly isolated from nasopharyngeal samples in New York City during this period of time. The mutation profiles of these samples suggest that there were no unusual substitutions particularly associated with GI infection (Fig. 2G and H). From two of these patients we also obtained SARS-CoV-2 genomes isolated from the nasopharynx. These strains were identical to those isolated from stool of the same patient (Fig. 2H).

Fecal inflammatory cytokines are elevated in COVID-19 cases but are unrelated to gastrointestinal symptoms

The detection of SARS-CoV-2 genome in the stool of multiple donors and the inverse association with fecal consistency led us to investigate if viral infection was causing intestinal inflammation associated with diarrhea among this cohort. Cytokines can be measured in fecal samples, and this method has been used as a sensitive and non-invasive way to monitor intestinal immune responses in inflammatory disease and during enteric infection [25-31]. In particular, fecal levels of IL-6, IL-8 and IL-1β are elevated in the context of acute bacterial or viral gastroenteritis and ulcerative colitis [28, 29]. We measured concentrations of eight cytokines and calprotectin in the stool samples from the COVID-19 cohort and in a group of stool samples collected from healthy donors recruited before SARS-CoV-2 was endemic in New York City (see Methods for details).

We first examined differences in fecal cytokines between COVID-19 patients and uninfected controls. Levels of fecal IL-8 were significantly elevated while levels of IL-10 were significantly lower in COVID-19 patients compared to uninfected controls (Fig. 3A), while the remaining tested cytokines were not significantly different. Although we observed elevated levels of IL-1β and TNFα in some COVID-19 patients, this was not consistent across the cohort (Fig. 3A). Next, we analyzed fecal cytokines in COVID-19 patients stratified by disease severity. IL-23 was the only cytokine that was significantly different between groups with higher a concentration in patients with severe COVID-19 disease (Fig. 3B). The remainder of fecal cytokines analyzed (IL-1β, IL-1ra, IL-6, IL-8, IL-17, TNFα) were not significantly different in COVID-19 patients compared to controls nor associated with severity of disease (Fig. 3A and B). We found no difference in the concentration of any fecal cytokine between those patients reporting diarrhea compared to those without diarrhea (Fig. S3A). Fecal concentrations of IL-1β, IL-6, IL-8and TNFα were not associated with concentrations of these cytokines measured in serum at the time of admission, although it is important to note that these samples were collected on different study days (median 16 days, range 2-66 days between samples; Fig. S3B). Fecal calprotectin, a marker of intestinal inflammation, was not associated with COVID-19 disease severity (p=0.12; Fig. 3C) nor stool viral genome load (p=0.22; Fig. 3D).

• Download figure
• Open in new tab

Fig. 3. Fecal cytokine levels in COVID-19 patients.

(A). Concentrations of cytokines in fecal samples from control donors and COVID-19 patients. Where multiple samples from the same donor were collected, the median of these samples is plotted. (B). Concentrations of fecal cytokines in COVID-19 patients with severe disease compared to those with non-severe disease. (C) Concentrations of fecal calprotectin in COVID-19 patients with severe disease compared to those with non-severe disease. (D) Concentrations of fecal calprotectin are unrelated to fecal viral load. Boxplots show median +/- IQR and P values are calculated by Mann-Whitney test.

The gut microbiome of COVID-19 patients is altered by antibiotic exposure but unrelated to disease severity or GI symptoms

A common feature of intestinal inflammation is disruption to the structure and diversity of the gut microbiota (dysbiosis) and there have been reports of gut microbiota changes in COVID-19 patients [7]. We evaluated the fecal microbiome of the COVID-19 patient cohort by both metagenomic shotgun and 16S rRNA gene amplicon sequencing. Analysis of both datasets revealed no features of composition, structure or diversity that associated with disease severity (Fig. 4A). Furthermore, we found no consistent changes in gut microbiome composition or richness in samples from patients reporting diarrhea, although the Shannon entropy, a metric of community richness and evenness, was slightly lower in samples with higher Bristol scores (p=0.034; Fig. 4B). Patients currently or recently treated with antibiotics, particularly vancomycin and/or ceftriaxone showed significant compositional disruption and reduced diversity as measured by both 16S (p=0.03) and metagenomic (p=0.0079) sequencing (Fig. 4C). These data further support a relative lack of severe inflammation, which would be expected to skew the composition of the microbiota, and suggests that antibiotic therapy, common among hospitalized COVID-19 patients, has a greater impact on the gut microbiome than SARS-CoV-2 infection itself, regardless of GI symptoms or clinical course.

• Download figure
• Open in new tab

Fig. 4. Gut microbiome and fecal SARS-CoV-2 genome sequencing of COVID-19 patients.

(A) PCoA (Bray-Curtis dissimilarity) of 16S rRNA gene amplicon and shotgun metagenomic sequencing of fecal samples from donors stratified by COVID-19 severity. (B) Alpha diversity (Shannon entropy) is not significantly altered in mild, moderate or severe COVID-19. (C) PCoA (Bray-Curtis dissimilarity) of 16S rRNA gene amplicon and shotgun metagenomic sequencing of fecal samples from donors with and without diarrhea. (D) Gut microbiome alpha diversity (Shannon entropy) in patients with and without diarrhea and as relates to the Bristol score of each sample. (E) PCoA (Bray-Curtis dissimilarity) of 16S rRNA gene amplicon and shotgun metagenomic sequencing of fecal samples from patients receiving or not receiving antibiotic therapy. (F). Gut microbiome alpha diversity (Shannon entropy) is reduced in hospitalized COVID-19 patients recently or currently receiving antibiotic therapy (Abx). Boxplots show median +/- IQR.

SARS-CoV-2 specific IgA is found in the stool of patients with severe COVID-19

A major component of the immune response against SARS-CoV-2 is believed to be the development of neutralizing antibodies against the viral Spike (S) protein, specifically the receptor binding domain (RBD) [32]. IgA is the major mucosal isotype and is known to provide long-term protection from enteric viruses, including rotavirus [33]. SARS-CoV-2-specific IgA has been detected in serum [32] and is found secreted in bronchoalveolar lavage fluid, nasal washes and in breast milk of convalescent donors [34-36]. We sought to determine if individuals with COVID-19 mounted an analogous gut mucosal IgA response against SARS-CoV-2, and if this was related to fecal viral load or other aspects of their clinical course.

We modified a direct enzyme linked immunosorbent assays (ELISA) assay developed to detect the presence of anti-RBD antibodies in the serum [32, 37] to measure fecal RBD-specific IgA (see Methods for details). Stool samples were treated with 0.5% Triton-X 100 for 1 hour to inactivate potentially infectious virus. Importantly, we found Triton-X 100 did not interfere with the detection of RBD-specific IgA (Fig. S4A).

The level of fecal RBD-reactive IgA was correlated with serum RBD-specific levels (R²=0.14, p=0.026, f-test; Fig. 5A). However, titers in feces were lower than in serum (Fig. 5A and S4B) and in contrast to serum we observed significant background signal at higher fecal concentrations, even among samples from healthy uninfected controls, which we attribute to non-specific binding of proteins from the high-complexity fecal samples (Fig. 5B). Across the cohort, there was no significant difference in the level of RBD-specific IgA detected in stool in healthy controls and COVID-19 patients (p=0.48, Mann-Whitney test; Fig. 5C), nor in the total amount of fecal IgA (p=0.76, Mann-Whitney test; Fig. S4C). Five out of forty-four COVID-19 patients (11%) demonstrated substantial fecal IgA reactivity to the SARS-CoV-2 RBD above the background level (AUC >3.5; Fig. 5B and 5C). The 5 patients with detectable levels of fecal RBD-specific IgA were all among those with severe COVID-19 (Fig. 5D) and notably, all these subjects were obese (median BMI of 40.7 vs 28.8 for those without detectable levels, p=0.007; Mann-Whitney test).

• Download figure
• Open in new tab

Fig. 5. SARS-CoV-2-specific IgA in stool of COVID-19 patients.

(A). Fecal and serum anti-SARS-CoV-2 RBD IgA are correlated. (B,C). Relative titers of anti-SARS-CoV-2 S1 RBD IgA in feces of control and COVID infected individuals. (D). Relative titers of anti-SARS-CoV-2 S1 RBD IgA in feces of COVID infected individuals with varying disease severity. (E). Anti-SARS-CoV-2 S1 RBD IgA and viral load are not correlated. Each point shows data from a single sample. (F). Anti-SARS-CoV-2 S1 RBD IgA and anti-S1 IgA in feces are correlated. Each point shows data from a single sample.

Of those patients with detectable fecal RBD-specific IgA (AUC > 3.5), 4/5 (80%) reported diarrhea, but across the whole cohort fecal RBD-specific IgA was not associated with diarrhea, Bristol stool score or the level of viral RNA in feces (p=0.89, Wilcox test; Fig. 5E Fig. S4D and S4E). We sampled COVID-19 patients longitudinally and found an individual’s fecal IgA reactivity to RBD was relatively stable over the sampling period – up to 30 days for some donors (Fig. S4F), similar to what is observed in systemic circulation[38]. Additionally, a subset of samples was tested for fecal spike protein IgA reactivity which was strongly correlated with fecal RBD reactivity (Fig. 5F). As with fecal RBD specific IgA, we did not observe significantly different spike-specific fecal IgA between COVID-19 and control donors (Fig. S4H).

Patient recruitment

We daily screened all consecutive patients admitted to the Mount Sinai Hospital for COVID-19 infection from April 15 2020 to May 21 2020. The study was reviewed and approved by the institutional review board at our medical center (Protocol #: HS# 16-00512/ GCO# 16-0583). We enrolled patients who were able to give informed consent and were able to provide a stool sample. We excluded patients who were not able to provide consent at the time of SARS-CoV-2 testing (sedation, endotracheal intubation prior or at the time of Emergency Room arrival, dementia or any other condition that prevented the subject to give fully informed consent). As our focus was to evaluate the role of gastrointestinal (GI) symptoms in the course of disease, we actively enrolled all patients with documented diarrhea, nausea and / or vomiting on presentation. We also enrolled patients without GI symptoms presenting at the same time to have a comparison group. We enrolled a total of 66 patients. Three patients were asymptomatic from SARS-CoV-2 infection and were excluded. Of the remaining 63, we collected stool samples from 44 and only these patients are considered in these analyses.

Clinical data collection; questionnaire and chart review

A questionnaire regarding symptoms onset and type of symptoms experienced was administered to all subjects. Detailed questions regarding gastrointestinal symptoms were included in the questionnaire (Supplemental Item 1). Chart review was used to collect subjects’ demographic characteristics (gender, age, race and ethnicity), comorbidities and medical history, COVID-19 clinical presentation, course of disease, laboratory data and COVID-19 treatment. Subject medical history was reviewed, and relevant comorbidities (hypertension, diabetes mellitus, obesity, chronic kidney disease, underlying chronic lung or cardiac diseases, inflammatory bowel disease, chronic HIV-infection, autoimmune diseases, current or prior history of cancer, current or prior history of organ transplant, concomitant pregnancy) were recorded.

COVID-19 disease severity was defined based on our institutional algorithm. Mild disease was defined as having a peripheral oxygen saturation (O2 Sat) ≥ 94% and no evidence of pneumonia on chest radiographic imaging (CXR); moderate disease was defined as either evidence of pneumonia on CXR or O2 Sat < 94%; severe disease was defined as evidence of pneumonia on CXR and the need for use of supplemental oxygen in the form of non-rebreather mask, high flow nasal cannula or mechanical ventilation. Severe disease was further subclassified as with or without evidence of end organ damage defined as: need for vasopressors support, creatinine clearance < 30 mL/min, alanine aminotransferase (ALT) > 5 times the upper limit of normal. Disease severity was recorded at the time of admission. Subjects were classified based on their symptoms at presentation in patients with and without GI symptoms. Qualifying GI symptoms were diarrhea, nausea or vomiting (at least one episode of self-reported symptom per day prior to presentation). Date of overall symptoms and of GI symptoms were also recorded. Clinical outcomes were collected as maximum level of oxygenation required during hospitalization (room air, nasal cannula, non-rebreather mask, high flow nasal cannula or mechanical ventilation), need for Intensive Care Unit (ICU) admission and death or discharge from the hospital.

Laboratory data were also collected via chart review, we included results for: SARS-CoV-2 nasal PCR, SARS-CoV-2 antibodies, peripheral WBC and lymphocyte percentage, Aspartate Aminotransferase (AST), Alanine Aminotransferase (ALT), Alkaline Phosphatase (ALKP), C-Reactive Protein, Procalcitonin, Lactic Acid Dehydrogenase (LDH), Ferritin, D-Dimer and cytokines that were implemented as part of clinical care IL-6, IL-8, TNFα and IL-1β. Serum cytokines were measured using the Ella platform as previously described[23]. Laboratory values were collected at baseline (at the time of admission) and at the time of discharge or death. Peak values for each lab during the admission was also collected.

COVID-19 treatment data was also collected, including use of hydroxychloroquine, azithromycin, systemic anticoagulation, steroids, convalescent plasma, and remdesivir. Information regarding subject’s participation on double blinded clinical trials was recorded. Use of antibiotics and class of antibiotics was also captured. All deidentified clinical data is provided in Data File S1.

Fecal samples from 22 healthy donors were previously collected from donors recruited in New York City between 2015 and 2016 as part of an unrelated study reviewed by the Mount Sinai Institutional Review Board (HS# 11-01669) and were stored at –80C until analysis. The mean (SD) age in the control group was 52.1 (11.3) compared to 55.9 (15.1) of the COVID-19 cohort (p=0.3, t-test) and 9 healthy donors were male (40.9%) compared with 47.7% of the COVID-19 cohort.

Stool sample collection and processing

Stool donations were collected into sterile containers (Covidien Precision™ Stool Collector) and kept at 4°C for no more than 24 hours before processing. Working under aseptic conditions in a BSL2+ laboratory, the sample was subdivided into aliquots for subsequent analyses. For the measurement of antibodies and cytokines, stool aliquots of 150-600 mg were frozen at −80°C. Samples were thawed on ice and homogenized to a concentration of 500 mg/mL in 1x PBS supplemented with protease inhibitor cocktail (cOmplete EDTA-free, Roche (1 tablet in 10mL) and phenylmethylsulfonyl fluoride (0.3 mM); to prevent protein degradation) and Triton X-100 (0.5% final concentration; to inactivate SARS-CoV-2 virus [44]). Samples were homogenized by vigorous shaking in a Mini-Beadbeater (Biospec) without beads for 60 seconds, kept on ice for one hour to ensure virus inactivation, centrifuged at 12,000 x g for 20 minutes at 4°C and the supernatants collected. For measurement of fecal calprotectin, stool aliquots of 150 - 600 mg were placed into 1 mL PBS containing protease inhibitor cocktail (cOmplete EDTA-free, Roche; 1 tablet in 10mL), placed at 60°C for 1 hour to inactivate SARS-CoV-2 before freezing at −80°C. For SARS-CoV-2 qPCR, stool aliquots (89.2 mg ± 77.8 (SD)) were immediately placed into 0.5 mL of Trizol and homogenized in a Mini-Beadbeater for 2 minutes. The mass of stool used for extraction did not affect positivity rate (p=0.3, Mann-Whitney test). The samples were stored at −80°C until RNA extraction. For culture independent gut microbiome sequencing, stool aliquots were immediately added to 1.3 mL of extraction buffer (0.55 mL 25:24:1 Phenol:Chloroform:Isoamylalcohol, 0.28 mL 2x STE Buffer (20mM Tris, 2mM EDTA, 2mM NaCl, pH 8), 0.27 mL QIAquick PM buffer (Qiagen) and 0.2 mL 20% SDS) and bead beat for 2 minutes with ∼400ul of 0.1mm silica beads. After centrifugation (8000rpm for 1 minute), the aqueous phase was removed and stored at −80°C for additional purification as described below.

Blood sample collection and processing

Research study blood samples were collected with the next clinically indicated blood draw after enrollment and were processed within four hours of collection. Serum was collected in a Tiger-top tube (serum separator tube) and stored directly at −80 degrees after being spun.

SARS-CoV-2 PCR from fecal samples

Measurement of fecal cytokines, calprotectin and total IgA

Fecal cytokines were assayed by sandwich ELISA (DuoSet® ELISA kits, R&D systems (IL-1β, IL-1ra, IL-6, IL-8, IL-10, IL-17 and TNFα or ThermoFisher Human ELISA kit (IL-23)) following the manufacturer’s instructions. For analysis of IL-23, fecal supernatants were diluted 1:2 in the kit-provided Sample Diluent A. For IL-1β and IL-1ra, fecal supernatants were diluted 1:5 in PBS with 1% BSA and for the remaining cytokines fecal supernatants were diluted 1:2 in PBS with 1% BSA. For measurement of fecal calprotectin, frozen heat-inactivated stool aliquots were further diluted in extraction buffer (BULHLMANN, EK-CAL2) and fecal Calprotectin levels were determined using the fCAL ELISA kit (BULHLMANN, EK-CAL2) according to manufacturer’s instructions. Total IgA was measured by ELISA using anti-IgA capture (1μg/mL; Southern Cat #2050-01) followed by anti-human IgA-HRP (1:8000, Southern Cat #2053-05) as previously described[46].

To comply with local biosafety protocols, it was necessary to treat samples to inactivate any potentially infectious virus present in the feces. We found that heating samples to 60°C for 60 minutes, a commonly used method for inactivation [32], led to degradation of recombinant cytokines spiked into stool samples. Some cytokines, including IL-1β, were particularly sensitive (Fig. S5). In contrast, 0.75% Triton X-100 included in the homogenization buffer (to achieve a final concentration of 0.5% when fecal samples were suspended 1:2 w:v), which is known to inactivate SARS-CoV-2 [44], and did not interfere with these assays. This pretreatment was used for all samples where proteins were measured (Fig. S5).

ELISA for fecal SARS-CoV-2-specific IgA

The protocol was adapted from a published method [32, 37]. Polystyrene 96 half-area plates (Corning# 3690) were coated with 60μL of diluted SARS-CoV-2 RBD or S1 protein (2μg/mL in 1x PBS, produced in the laboratory of Dr. Florian Krammer) and incubated at 4°C overnight. Plates were washed with PBS plus 0.05% Tween20 (PBS-T) and blocked with PBS-T with 3% (w:v) milk powder (AmericanBio) for 2-3 hours. Fecal extract (pre-diluted 1:2 w:v, see above) or serum was mixed with equal volume PBS-T with 1% (w:v) milk powder and serially diluted two-fold for eight iterations. 30μL of the analyte dilution series was added to the plate for 2 hours. Plates were washed, and 60μL of anti-IgA-horseradish peroxidase (HRP) (Southern, 1:8000) or anti-mouse IgA (Millipore Sigma, 1:8000) was added. After 2 hours, ELISA plates were again washed and developed using TMB substrate (Biolegend). Plates were quenched with 2N sulfuric acid and read at 450nm. A titrated serum sample from a seropositive donor was included on every assay plate and used to normalize data between plates, correcting for minor batch variation.

16S and metagenomic microbiome sequencing

DNA extraction buffer (4.3 mM Tris, pH 8; 0.4 mM EDTA, 43 mM NaCl, 3% v/w SDS, 42% v/v phenol/chloroform/IAA and 21% v/v QIAquick PM buffer) containing 400 μl 0.1 mm Zirconia/Silica Beads was added to fecal samples and the samples were bead-beat on a BioSpec Mini-Beadbeater-96 for 2 minutes. After centrifugation (5 min, 4000 rpm), the ∼400 μl of aqueous layer was mixed with 650 μl QIAquick PM buffer DNA was purified using the QIAquick 96 PCR Purification Kit (Qiagen) according to manufactures protocol. DNA concentration was quantified using the Quant-iT dsDNA Assay Kit, broad range (Life Technologies) and normalized to 2 ng/μl on a Beckman handling robot. Amplicon preparation and sequencing was performed as previously described [47]. Briefly, bacterial 16S rDNA PCR including no template controls were setup in a separate PCR workstation using dual-indexed primers. PCR reactions contained 1 µM for each primer, 4 ng DNA, and Phusion Flash High-Fidelity PCR Master Mix (Thermo Fisher Scientific). Reactions were held at 98°C for 30 s, proceeding to 25 cycles at 98°C for 10 s, 45°C for 30 s, and 72°C for 30 s and a final extension of 2 min at 72°C. Amplicons were evaluated by gel electrophoresis. The sequencing library was prepared by combining equivolume amounts of each amplicon, size-selected and concentrated using AMPure XP beads (0.8X, Beckman). Library concentration was quantified by Qubit and qPCR, mixed with 15% PhiX, diluted to 4 pM and subjected to paired-end sequencing (Reagent Kit V2, 2×150bp) on an Illumina MiSeq sequencer. Resulting fastq files were analyzed using QIIME2 version 2019.10[48] and the DADA2[49] denoised-paired plugin with a truncation length of 150bp and 145bp for the forward and reverse read, respectively. Amplicon Sequence Variants (ASVs) were classified using the Scikit-Learn plugin [50] using the Naive Bayes classifiers trained on the silva-132-99-515-806-nb-classifier. Resulting ASV tables were filtered using a minimum depth of 5000, a minimal ASV frequency of 10 and minimal sample frequency of 2 and core metrics were calculated using the Qiime2 core metrics plugin.

Metagenomic sequencing was performed as previously described [51]. Briefly, DNA (extracted as for 16S rRNA amplicon sequencing) was sonicated and Illumina sequencing libraries generated using the NEBNext Ultra II DNA Library Prep Kit. Ligation products of 500– 600 base pairs were purified using SPRIselect beads (Beckman Coulter) and enrichment PCR performed. Samples were pooled in equal proportions and size-selected using AMPure XP beads (Beckman Coulter) before sequencing with an Illumina HiSeq (paired-end 150 bp). One sample yielding less than 100,000 reads was excluded from further analysis. Remaining samples were sequenced to an average depth of 1.52 × 10⁶ reads (SD 3.73 × 10⁵) per sample. Reads were trimmed with Trimmomatic [52] and taxonomic assignments were generated with MetaPhlAn2^[53]. Data is deposited to NCBI under BioProject PRJNA660883.