Detection of SARS-CoV-2 in fecal samples with different pretreatment methods and PCR kits

Background Gastrointestinal symptoms are common in COVID-19 patients and SARS-CoV-2 RNA has been detected in the patients’ feces, which could lead to fecal–oral transmission. Therefore, fecal sample testing with real-time RT-PCR is highly recommended as a routine test for SARS-CoV-2 infection. However, varying rates of detection in fecal sample have been reported. The aim of this study was to provide insights into the detection rates of SARS-CoV-2 in COVID-19 patients’ fecal sample by using four real-time RT-PCR kits and two pretreatment methods (inactive and non-inactive). Results The detection rate of Trizol pretreatment group was slightly higher than that of Phosphate Buffered Saline (PBS) groups, showing that pretreatment and inactivation by Trizol had no influence to SARS-CoV-2 nucleic acid test (NAT) results. 39.29% detection rate in fecal sample by DAAN was obtained, while Bio-germ was 40.48%, Sansure 34.52%, and GeneoDx 33.33%. The former three kits had no significant difference. The DAAN kit detection rates of ORF1ab and N gene were nearly equal and Ct value distribution was more scattered, while the Bio-germ kit distribution was more clustered. The positive rate of SARS-COV-2 in fecal samples correlated with the severity of the disease, specifically, severe cases were less likely to be identified than asymptomatic infection in the DAAN group (adjusted OR 0.05, 95%CI = 0.00 ~ 0.91). Conclusions Trizol should be of choice as a valid and safe method for pretreatment of fecal samples of SARS-CoV-2. All real-time RT-PCR kits assessed in this study can be used for routine detection of SARS-CoV-2 in fecal samples. While DAAN, with high NAT positive rate, could be the best out of the 4 kits used in this study. SARS-CoV-2 positive rate in fecal sample was related to the severity of illness. Supplementary Information The online version contains supplementary material available at 10.1186/s12866-021-02118-0.


Introduction
The Coronavirus Disease-19 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in an ongoing pandemic. Globally, as of 27 September 2020, there have been 32,730,945 confirmed cases of COVID-19, including 991,224 deaths [1]. Unlike other human coronaviruses, SARS-CoV-2 infects the lower respiratory tract and could lead to the severe respiratory symptoms. While the majority of infections lead to mild clinical manifestations [2,3], such as cough, fever and malaise, some patients deteriorate to acute respiratory distress syndrome, multi-organ failure, shock, and blood clots, with high risk for a fatal outcome [2,4,5]. As SARS-CoV-2 invades cells through binding to angiotensin-converting enzyme 2 (ACE2) [6], which is expressed in various tissues apart from the lungs [7], COVID-19 could cause multiple organ involvements including gastrointestinal infection [8,9], and fecal samples have tested positive for the nucleic acid of SARS-CoV-2, however, this could be simultaneous with respiratory tract samples testing negative [10]. Real-time reverse transcriptase-polymerase chain reaction (realtime RT-PCR) assay is one of the most sensitive and specific assay for viral nucleic acid test (NAT) and is widely recommended for the detection of SARS-CoV-2 [11,12]. Different real-time RT-PCR positive rates in fecal samples have been reported by several studies [13][14][15][16], possibly due to different RNA isolation methods, or real-time RT-PCR kits. Coronaviruses belong to the Coronaviridae family in the Nidovirales order, with positivesense RNA that can express its replication and transcription complex like RNA-dependent RNA polymerase (RdRp) from open reading frame 1ab (ORF1ab), a large open reading frame [17]. Coronavirus structural proteins, including envelope protein (E), nucleocapsid protein (N), and spike protein (S), are expressed by generating subgenomic messenger RNAs [18]. The ORF1ab, N, RdRP, S and E genes are the main targets of the viral RT-PCR detection protocols recommended by institutes worldwide [19].
The detection of SARS-CoV-2 in feces has recently aroused the concern about the potential fecal-oral or fecal-respiratory transmission route [16,[20][21][22][23]. Therefore, more precise detection of the virus in the feces is essential for understanding and preventing viral transmission. Nevertheless, fecal samples are far more complicated than nasopharyngeal/oropharyngeal swabs, sputum, and other respiratory samples, therefore it is critical to design an appropriate methodology to implement optimal sampling and RNA extraction procedures, for biases introduced from this process could influence detection.
In our study, fecal samples were collected from confirmed COVID-19 cases following the Guidelines supported by the National Health Commission of China (Version 7) [24]. Fecal samples were pretreated with Phosphate Buffered Saline (PBS) and Trizol, respectively. Few studies have assessed different kits for the NAT of SARS-CoV-2 in feces, so we screened four real-time RT-PCR kits in this study. Moreover, the demographical and clinical characteristics, severity of illness, fecal sample category, sampling interval and general information of kits used for viral NAT were compared between those cases. Our study provides evidence for the fecal positive rate of patients diagnosed with COVID-19 and effective selection of pretreatment methods and SARS-CoV-2 detection kits for fecal samples.

Description of the population
For the 90 cases enrolled, median age was 46 years (interquartile range [IQR 35-59], range 7-83 years), with a balance gender ratio (M/F ratio, 42/48). Most patients were of moderate illness (56.67%), followed by severe illness (20.00%). Asymptomatic infection, mild, and critical patients were less than 10% each. Most fecal samples were collected 1 to 2 weeks after illness onset. The majority of fecal samples were categorized as formed stool, 6 cases of diarrhea were also collected. The patients covered a wide variety of careers, including peasant, unemployed or housewife (Table 1).

Pretreatment with Trizol had no influence on NAT
Two pretreatment methods were used to suspend 38 fecal samples. The double positive rates of PBS, Trizol groups were 55.26% (21/38) and 60.53% (23/38), respectively, when using Bio-germ Kit. The double positive rates of GeneoDx were 39.47% (15/38) and 44.74% (17/ 38), respectively. Paired chi-square test was used for analysis, and it was found that there was no statistical difference in the positive rate of the same kit using different pretreatment methods (P > 0.05), and there was no statistical difference between the two kits. However, the detection rate of Trizol group (60.53, 44.74%) was slightly higher than that of PBS group.
The Ct values of ORF1ab and N genes detected by the 4 kits were different to some extent (Fig. 1a). The positive numbers of ORF1ab and N gene detected with DAAN kit were nearly equal, i.e., 35 and 33, respectively. Followed by, Bio-germ, 36 and 40, respectively. The double positive rate of these two kits was also the highest (39.29, 40.48%), and the Ct value distribution of Biogerm was more clusterred, while that of DAAN was more scattered. The Ct values of ORF1ab gene ranged from 28 to 37.1, and N gene from 25.  Fig. 1b. Interestingly, the positive rate of Biogerm was still higher (34/84) than that of GeneoDx (22/ 66) after enlarging the sample size.
The positive rate of SARS-CoV-2 correlates with the severity of the disease When using kits from DAAN and Bio-germ, the positive rates of different clinical types were statistically significant (Fisher's Exact test, two-tailed, P < 0.05) (Additional file 1: Supplementary Table 1). Univariate logistic regression models were used to test the correlation between individual covariates and the outcome of a positive test, and multivariate logistic regression was used to detect variables that were associated with a positive test for SARS-CoV-2 after multiple imputations of missing values. In univariable analysis, there was no significant association between NAT positive rate and severity of illness. Moreover, in multivariable analysis, after adjusting for age, gender and fecal samples category, it was found that fecal samples from severe cases were less likely to be identified than asymptomatic infection when using the DAAN kit (adjusted odds ratio (OR) 0.05, 95%confidence interval (CI) = 0.00-0.91) ( Table 3). While a similar correlation was not observed in the Biogerm kit group (Additional file 1: Supplementary Table 2). This may have been due to the small sample size and the small number of grids in some clinical classifications, which makes it hard to calculate the risk coefficient or the large confidence interval.

There was no correlation between sampling interval and positive results
Logistic regression models were used to test the correlation between sampling interval and the outcome of a positive test. The sampling interval was defined as the days from illness onset to the sampling time (≤ 7 d, 7-14 d and > 14 d). We attempted to analyze NAT results from all the kits and adjusted for age, gender and fecal samples category, but no obvious differences were found in our results (Additional file 1: Supplementary Tables 3-6).

Discussion
With the rapid spread of the SARS-CoV-2, there have been a large number of research studies on the detection of the viral NAT, such as the positive detection rate in different clinical samples [16,25], comparison of clinical performance of different real-time RT-PCR kits [26], etc. But few articles mentioned the process of fecal samples pretreatment and the impact of different pretreatment methods on NAT. The composition of fecal samples is more complex than that of respiratory samples and there are many impurities, hence proper pretreatment must be performed before viral RNA isolation and purification. At present, the method of vortexing with NS or PBS to prepare a feces suspension followed by centrifuging to create a supernatant for RNA extraction is typically adopted. However, the infectivity, transmission ability, and survivability of this novel virus remain unclear. Vortexing and centrifuging are high-risk biohazardous operations which produce laboratory aerosols. Trizol is a common RNA extraction reagent, comprised of phenol and guanidine isothiocyanate, can lyse cells, ensure RNA integrity, and has a strong protein denaturation effect. Some studies have proven that Trizol could inactive high-titer viruses, including SARS-CoV-2 [27,28]. Therefore, Trizol can not only extract RNA effectively and ensure the integrity of RNA, but also inactivate virus, reducing infectivity and ensure biosafety. In our study, PBS and Trizol pretreatments were compared. The detection rate of the Trizol group was slightly higher than that of PBS group, though there was no statistical difference. So we concluded that fecal samples treated with Trizol can not only inactive the virus and reduce the risk of the experimental procedure, but also had no influence on the downstream NAT. Trizol is a monophasic solution which solubilizes biological material and denatures protein simultaneously [29]. Trizol can be recommended to treat fecal samples to help improve the sensitivity with low limit of detection, and maintain viral RNA integrity, and samples could be inactivated in Trizol, thus reduce infectivity and the chance of exposing healthcare workers [30]. However, further research is required before Trizol could be used as a routine reagent in laboratory protocols related to SARS-CoV-2.
The positive rates of SARS-CoV-2 in fecal samples of COVID-19 patients testing by real-time RT-PCR ranged from 25 to 82% [10,16,22,[31][32][33][34][35]. An astonishing study reported that 39 (53.42%) out of the 73 hospitalized COVID-19 patients' fecal samples were tested SARS-CoV-2 positive and remained positive for 1-12 days, 17  (23.29%) of whom persisted even after the respiratory samples turned negative for the viral detection [22]. In our study, all the detection rates in feces were more than 30% (DAAN 39.29%, Sansure 34.52%, Bio-germ 40.48%, GeneoDx 33.33%). The difference of detection rates of SARS-CoV-2 in fecal samples may result from different sample sizes, different kits and the various and unstandardized sampling time. Accumulating evidence has supported the potential for feces associated transmission of SARS-CoV-2 [8,10,16,22]. By evaluating the positive detection rate of the four kits, the distribution of Ct values, and the positive numbers of N and ORF1ab genes, we recommend the use of DAAN for NAT of fecal samples. Other research groups have found co-expression of ACE2 and transmembrane serine protease 2 (TMPRSS2) in the enterocytes, progenitor, and stem-like epithelial cells of the lower gastrointestinal tract, especially in the small intestine [36]. It is well established that the invasion of SARS-CoV-2 depends on the interaction of spike protein with ACE2 and TMPRSS2 [6,37,38]. In addition, staining of ACE2 and SARS-CoV-2 was simultaneously observed in gastrointestinal epithelium from those who tested positive for SARS-CoV-2 in fecal samples [22]. Surprisingly, a recent report showed that the toilet bowl, sink, and door handle in the room where the COVID-19 patient had resided were contaminated by SARS-CoV-2 [39]. This revealed that the SARS-CoV-2 can transmit through feces as well as respiratory droplets or contact transmission. Therefore, we believe, prevention of fecaloral transmission should be considered, and all COVID-19 patients' fecal samples should be as a routine test for the SARS-CoV-2. Once it is positive, the special regulations and nursing strategies must be used to prevent spread of the virus. Further research is needed to understand to what extent SARS-CoV-2 is transmitted via fecal-oral route. This will help decide whether to test the virus in the feces of COVID-19 patients.

Fecal sample collection
A total of 90 fecal samples of confirmed COVID-19 cases in Sichuan Province, China were collected and stored in sterile containers at − 80°C.

Pretreatment for the fecal samples and RNA extraction
To study the influence of different pretreatment methods on SARS-CoV-2 NAT results, 38 fecal samples were selected randomly to suspend in two ways. To be specific, 200 mg of stool was suspended in a 15 mL tube containing 2 mL PBS, 2 mL Trizol (Trizol @ Reagent, invitrogen, USA), respectively. The mixture was stirred gently and mixed. Keep still for 10 min and then supernatant were transferred to a new tube. Depending on the NAT results of the 38 randomly selected fecal samples, the remaining 52 samples were suspended in Trizol. 200 μL supernatant were used for RNA isolation by the NP968 Nucleic Acid Extraction System (Xi'an Tianlong Science & Technology Co., LTD, Xi'an, China). Extracted RNA was stored at − 20°C until use.
Real-time RT-PCR assays for the detection of SARS-CoV-2 Group 1 To find the best pretreatment method, the viral RNA extracted from 38 fecal samples suspend in two ways were tested with two kits (Shanghai Bio-germ Medical Co., Ltd., and Shanghai GeneoDx Biotech Co., Ltd), respectively.

Group 2
To study positive rate of SARS-CoV-2 in fecal samples and evaluate real time RT-PCR kits, the viral RNA extracted from 84 fecal samples suspend in Trizol were tested with four kits (Shanghai Bio-germ Medical Co., Ltd., and Shanghai GeneoDx Biotech Co., Ltd. DAAN Gene Co., Ltd. of Sun Yat-sen University and Sansure Biotech Inc.), respectively. As 6 stool samples did not contain enough material for RNA extraction, only 84 samples were tested in group 2. Due to a supply shortage, only 66 fecal samples in group 2 were tested with the GeneoDx kit. As previously described [40], only when the two target genes (ORF1ab and N) are simultaneously positive, can a positive detection of SARS-CoV-2 be reported.

Characteristics of the selected kits
Here, we provided a comparison of four readily available COVID-19 real-time RT-PCR kits from different  Table 4). Two of these kits have been proven to detect low concentration viral RNA: Sansure can achieve 1 copy/reaction, and Bio-germ can achieve 10 Copies/reaction [41]. The 4 commercial kits have been approved by the China National Medical Products Administration (NMPA), and 2 have received CE (Conformité Européenne) marking (DAAN and Sansure). All real-time RT-PCR kits assessed in this study have been used for routine diagnostics of SARS-CoV-2 by experienced molecular diagnostic laboratories.

Data analysis
Statistical analyses were performed by using SPSS26.0 software. All statistical tests were two-sided, and significant differences were considered at P < 0.05. Continuous variables were evaluated using the median and interquartile range (IQR) values. Chi-square or Fisher's exact tests were utilized to compare the proportions of the categorical variables. Chi-square test was used to compare inter-group differences, and logistic regression analysis was performed to analyze the risk factors for SARS-CoV-2 prevalence. Crude and adjusted risk ratios and 95%CIs for NAT positive rate were calculated using modified Poisson regression. Risk ratios were adjusted for age (as a continuous variable), gender, and fecal samples category.

Supplementary Information
The