key: cord-1015791-re68t4yf authors: Majumdar, Ramanath; Vrana, Julie A.; Koepplin, Justin W.; Milosevic, Dragana; Roden, Anja C.; Garcia, Joaquin J.; Kipp, Benjamin R.; Moyer, Ann M. title: SARS-CoV-2 RNA Detection in Formalin-Fixed Paraffin-Embedded (FFPE) Tissue by Droplet Digital PCR (ddPCR) date: 2022-05-10 journal: Clin Chim Acta DOI: 10.1016/j.cca.2022.05.007 sha: ab4be615db6892ecc32593bd2724c2201d71917a doc_id: 1015791 cord_uid: re68t4yf BACKGROUND: SARS-CoV-2 is an RNA virus that causes primarily causes respiratory disease; however, infection of other tissue has been reported. Evaluation of SARS-CoV-2 in tissue specimens may increase understanding of SARS-CoV-2 pathobiology. MATERIALS AND METHODS: A qualitative test for detection of SARS-CoV-2 in formalin-fixed paraffin-embedded (FFPE) tissues was developed and validated using droplet digital PCR (ddPCR), which has a lower limit of detection than reverse transcription (RT)-qPCR. After extraction of total RNA from unstained FFPE tissue, SARS-CoV-2 nucleocapsid (N1, N2) target sequences were amplified and quantified, along with human RPP30 as a control using the Bio-Rad SARS-CoV-2 ddPCR kit. RESULTS: SARS-CoV-2 was detected in all 21 known positive samples and none of the 16 negative samples. As few as approximately 5 viral copies were reliably detected. Since January 2021, many tissue types have been clinically tested. Of the 195 clinical specimens, the positivity rate was 35% with placenta and fetal tissue showing the highest percentage of positive cases. CONCLUSION: This sensitive FFPE-based assay has broad clinical utility with applications as diverse as pregnancy loss and evaluation of liver transplant rejection. This assay will aid in understanding atypical presentations of COVID-19 as well as long-term sequelae. SARS-CoV-2 is a positive-sense, single-stranded RNA virus that causes life-threatening coronavirus disease 2019 in humans. SARS-CoV-2 infection can result in diverse, multiorgan pathology, the most significant being in the lungs, heart, kidney, central nervous system, liver, lymph nodes, bone marrow, vasculature, intestine, and placenta [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] . SARS-CoV-2 infection is usually identified by the detection of viral RNA using reverse transcriptase polymerase chain reaction (RT-qPCR) on nasopharyngeal or oropharyngeal swabs. Current or prior infection might also be recognized by a serum test for SARS-CoV-2 IgG antibodies. A subset of patients may present without the typical upper respiratory symptoms or with concurrent unusual features that potentially represent long-term sequelae of SARS-CoV-2 [11, 12] . Patients with current or previous SARS-CoV-2 infection who present with "unexplained" acute or acute on chronic respiratory failure or other organ-specific manifestations may undergo tissue biopsy. Furthermore, an autopsy may be performed for individuals with sudden unexplained death. Therefore, a formalin-fixed paraffin-embedded (FFPE)-tissue-based test was needed for evaluation of SARS-CoV-2 for both improved understanding of potential disease manifestations as well as for patient management. Various methods for the identification of SARS-CoV-2 infection in FFPE-tissue specimens have been published including in situ hybridization (ISH), immunohistochemistry (IHC), RT-qPCR, and electron microscopy [1, 6, 11, 13] . Moreover, due to various factors in the preanalytical, analytical, and postanalytical phase, virus detection in tissue biopsy samples presents specific diagnostic challenges. Therefore, potentially more sensitive techniques such as reverse transcription-droplet digital polymerase chain reaction (RT-ddPCR) might enhance detection of SARS-CoV-2 in FFPE tissue or other matrices [11] . In fact, ddPCR has been shown to have a significantly lower limit of detection of SARS-CoV-2 in nasopharyngeal swabs than reverse transcription RT-qPCR [14] . Consequently, detection of SARS-CoV-2 in deceased patients (autopsy tissue) or surgical specimens may confirm a suspected diagnosis among individuals with clinical and/or pathologic manifestations of COVID-19 and may increase understanding of SARS-CoV-2 pathobiology. The use of ddPCR for nucleotide target analysis has emerged in recent years as a robust and reliable methodology [15] [16] [17] . ddPCR can specifically measure single or multiple nucleotide targets in a single reaction [15, 16, 18] . With ddPCR, DNA or reverse transcribed cDNA is partitioned into thousands of droplets that are subsequently amplified, and fluorescently-labeled probe signals within each droplet are recorded as either positive or negative, depending on the presence or absence of a nucleotide target. With thousands of individual droplets, a Poisson calculation of target copies per droplet [−ln(1 − p)], where p is the fraction of positive droplets, is applied and the absolute concentration for each nucleotide target is obtained, providing a highly accurate estimation of nucleotide copy number that is crucial for clinical diagnostics [15, 18, 19] . Here the development and performance characteristics are described for a qualitative RT-ddPCR test for detection of SARS-CoV-2 in FFPE tissues, which can be implemented in most clinical molecular diagnostic laboratories. Thirty-seven tissue samples from 35 autopsy cases (20 positive, 15 negative), a cell block from a cell line infected with SARS-CoV-2, and a cell block from a cell line infected with influenza virus were utilized in this study for accuracy, precision, stability, linearity, and specificity studies. Multiple tissue types were used for test development, including lung, heart, brain, and liver tissue. Slides from autopsy cases of patients who tested positive for COVID-19 (n=20) by antemortem nasopharyngeal or oropharyngeal swab were reviewed by a pathologist. Blocks with histologic features suggestive of infection were tested by RNA in situ hybridization (ISH) and immunohistochemistry, as previously described, for comparison to the ddPCR assay [11] . Tissue samples from patients who underwent autopsy the year prior to the beginning of the COVID-19 pandemic were included as negative controls. These patients had lung pathology not due to SARS-CoV-2, such as diffuse alveolar damage or influenza virus infection (n=15). Note that although in the tables the positive samples are represented sequentially by lower numbers and the negative samples are represented sequentially by higher numbers, they were intermixed and the test developers were blinded as to which were positive and negative. They have been re-named and grouped for the purposes of publication. This assay was launched as a clinically orderable test in January 2021. All specimens analyzed clinically up to November 20, 2021 were included in a retrospective summary describing clinical applications of this test as well as positivity rates. This study was approved by the Mayo Clinic Institutional Review Board. An overview of the RT-ddPCR assay workflow is shown in Figure 1 . Total RNA, including viral RNA, was extracted from unstained FFPE tissue scrolls cut at 10 microns by either RNeasy® DSP FFPE Kit or miRNeasy FFPE Kit (Qiagen, Hilden, Germany). The RNA was quantified, ranging from 12-800 ng/μL, and stored at -80°C. The assay was performed according to the manufacturer protocol, as outlined in the "Instructions for Use" for the Bio-Rad SARS-CoV-2 ddPCR Test (Bio-Rad Laboratories, CA, USA), which had received Emergency Use Authorization (EUA) from the US Food and Drug Administration (FDA) for nasopharyngeal specimens at the time of development of this tissue-based assay. Briefly, 5.5 uL of RNA eluate (corresponding to a range of 66-800 ng of RNA) was added to 16.5 uL of ddPCR master mix containing 5.5 μL of 4x One Step-RT-ddPCR Supermix, 2.2 μL of reverse transcriptase, 1.1 μL of dithiothreitol, 1.1 μL of the 20x 2019-nCoV CDC ddPCR triplex probe. RNA samples were run both neat and at a 1:100 dilution. After the reaction was partitioned into droplets using the Bio-Rad AutoDG instrument, SARS-CoV-2 nucleocapsid (N1 and N2) target sequences, along with human ribonuclease P/MRP subunit 30 (RPP30) target as a control, were amplified using Veriti thermocyclers (Thermo Fisher Scientific, Inc., Waltham, MA) using the following protocol: RT steps: 25°C for 3 minutes and 50°C for 60 minutes, followed by the PCR steps: 95°C for 10 minutes, then 40 cycles of denaturation at 95°C for 30 seconds and annealing/extension at 55°C for 1 minute, and a final enzyme deactivation at 98°C for 10 minutes. Droplets were counted as negative or positive for one or more target by the QX200 Droplet Reader and signal data were analyzed using QuantaSoft Analysis Pro software version 1.0.596 (Bio-Rad Laboratories, Hercules, CA, USA). For each run, a low positive for SARS-CoV-2 and a negative sample were included for quality control. To evaluate the potential impact of the FFPE matrix and establish the limit of detection (LOD), samples with total SARS-CoV-2 RNA copies of 80, 40, 20, 10, and 5 were prepared in triplicate and samples with 0 copies were prepared in pentaplicate using Exact Diagnostics SARS-CoV-2 Standard (Exact Diagnostics LLC, Fort Worth, TX) mixed with extracted RNA from pooled SARS-CoV-2 negative FFPE tissues (523 ng total RNA in each reaction) as a diluent for testing with the RT-ddPCR assay. The LOD was determined based on the lowest signal that could be reproducibly detected and differentiated from background signal observed when no SARS-CoV-2 was present. All 37 known SARS-CoV-2 positive (n=21) and negative (n=16) tissues and cell blocks were utilized to determine accuracy. Intra-assay and inter-assay precision (reproducibility) were evaluated using pooled extracted RNA from two patient FFPE tissue specimens (one negative and one positive, pulmonary derived), one positive quality control (QC) sample derived from a SARS-CoV-2 embedded FFPE cell block, and one negative QC sample derived from an FFPE cell block of tissue infected with influenza A virus. For the intra-assay precision study, each of the four samples were run in triplicate together on the same run. To assess inter-assay precision, each sample was run on three different runs. Concordance of both the qualitative positive or negative result was evaluated, along with the quantitative value of N1 and N2 SARS-CoV-2 target. Extensive in silico analyses conducted by the kit manufacturer found no potential cross-reactivity of common respiratory tract pathogens (bacteria, fungi, and viruses including other coronavirus strains), and our laboratory previously identified no cross reactivity with other coronavirus strains when testing nasopharyngeal swab specimens (unpublished data). To further verify the analytical specificity in the FFPE matrix, two blocks of rhesus monkey kidney (RMK) cells infected with Influenza A and Parainfluenza were tested. QuantaSoft Analysis Pro software version 1.0.596 (Bio-Rad) was used for droplet cluster classification. The QX200 automated droplet reader counts every acceptable droplet and measures the fluorescence emissions from each droplet using 2 channels (FAM and HEX). Two droplet clusters and up to eight droplet clusters were observed for a normal (negative) and COVID-19 patient (positive) specimen, respectively, as illustrated in In some cases with strong RPP30 expression, a small number of N1+RPP30 and N2+RPP30 double positive droplets were observed without negative or single-positive droplets. The number of droplets generated per reaction ranged from 9000 to 18299 droplets. Any sample with <9000 droplets was rejected and repeated. SARS-CoV-2 N1 and N2 targets were detected for all 3 of the triplicate samples prepared with 5 virion copies ( Table 1 ). The three replicates had 7, 11, and 12 total positive N1 and N2 droplets, respectively. Two N1-positive droplets were observed in one replicate out of the 5 blank samples prepared with 0 copies of viral RNA, while the remainder had 0 positive N1 or N2 droplets. Additionally, the no template control (NTC) was run in triplicate with no N1 or N2 positive droplets. Therefore, the limit of blank was established at 2 positive droplets and the cut-offs for final calls were set at ≥5 N1 and N2 droplets = positive, 0-2 positive droplets = negative, and 3-4 positive droplets = indeterminate. Based on this experiment, 5 copies of viral RNA (equivalent to approximately 5 droplets of N1 and/or N2) can be reliably detected and is expected to be the approximate limit of detection. Additional linear regression analysis of the average observed copies of SARS-CoV-2 RNA detected by the RT-ddPCR assay as compared to the quantity expected in the sample based on the amount added demonstrates that the assay is linear and is shown in Figure 3 . RPP30 values may differ based on tissue type and size; therefore, a cut-off of 50 positive droplets was selected to ensure that RNA was successfully extracted. Total (100%) agreement of qualitative results (detected vs. not detected) was observed for 21 known-positive samples. Near complete agreement (15 of 16 samples; 93.8%) was observed for 16 known-negative samples ( Table 2) . One sample had a total of 4 positive droplets, which would be called indeterminate with the established cut-offs. This sample was heart tissue from autopsy of a patient known to be SARS-CoV-2 positive. As shown in Table 2 , the assay was saturated for several positive samples (N1 and N2 copies exceeded upper limits of approximately 20,000,000, but when diluted 1:100 with 1x TE buffer, these samples displayed the expected droplet patterns. Both intra-and inter-assay measurements were 100% concordant in qualitative positive or negative calls among all replicates for each sample. Although this is a qualitative assay, the coefficient of variation (CV) of the number of positive droplets detected and copy number were also evaluated. Across the 3 intra-assay and inter-assay measurements, the negative QC sample had 0 positive N1 and N2 droplets in each replicate, with the exception of 1 N2 droplet (corresponding to 2 copies of N2 in the reaction well) in the intra-assay study. This would have been considered negative in our qualitative assay. Similarly, the negative patient sample had one replicate in the inter-assay study with 1 N2 droplet (corresponding to 2 copies of N2 in the reaction well), and in the intra-assay study there was one replicate with 1 N1 droplet (corresponding to 1 copy of N1) and one replicate with 1 N2 droplet (corresponding to 2 N2 copies in the reaction). Again, each replicate would have been called negative in the qualitative assay. The positive QC sample had a mean of 16,575 copies of N1 per reaction detected, with 10% CV in the inter-assay study and 3% CV in the intra-assay study. This sample had a mean of 15,175 copies of N2 detected per reaction, with 10% CV in the inter-assay study and 3% CV in the intra-assay study. The positive patient sample had a mean of 156,938 copies of N1 detected per reaction, with 3% CV in the inter-assay study and 4% CV in the intra-assay study. This sample had a mean of 174,836 copies of N2 detected per reaction, with 10% CV in the interassay study and 8% CV in the intra-assay study. In this qualitative assay, each replicate of both the positive control and positive patient sample would have been called positive. No amplification of either of the nucleocapsid targets was detected (i.e., 0 droplet count for both N1 and N2) for any of the non-SARS-CoV-2 viruses tested. From January 22, 2021, which was the first clinical run using this test, until November 20, 2021, a total of 195 specimens have been run for clinical purposes (Table 3, Figure 4) . Additional testing has been performed on a research basis and not included in this number. Testing has been ordered both by physicians at our institution as well as through our reference laboratory practice. The overall SARS-CoV-2 positivity rate in clinical testing was 35% with placenta and fetal tissue showing the highest percentage of positive cases ( Table 3) . Lung tissue was most commonly tested. Among the 25 positive cases were 9 autopsies. Six of these patients were known to be COVID-positive prior to death, though three of those had been positive several weeks prior to death. The remaining three had no known history of a COVIDtest. The other positive cases were biopsy specimens, four of which were from patients with a remote history of COVID, while no clinical information was provided for the others. Among the negative lung specimens, 13 were from autopsy cases, two of which had a remote history of COVID -including one who had remained hospitalized for the 3 months since diagnosis. One patient had multiple negative COVID tests prior to death, and no COVID testing was known to have been performed on the remaining patients. The additional negative lung tissue specimens were biopsies and included two patients with a remote history of COVID while no clinical history was provided for the remainder. SARS-CoV-2 has been described to result in a variety of pathological manifestations affecting most organ systems. Because formalin fixation can impact the quality and integrity of nucleic acid, published methods for detection of SARS-CoV-2 are either not suitable for testing FFPE tissue specimens or lack sensitivity, and thus might not be helpful or feasible in routine clinical practice [1, 6, 7, 13, 20] . The aim of this study was to develop and validate an easy, reliable, and accurate assay using ddPCR technology for the detection of SARS-CoV-2 in FFPE tissue for routine clinical use. We have developed and validated a qualitative test protocol that directly detects the SARS-CoV-2 virus in human FFPE tissues. Of note, in this assay SARS-CoV-2 detection was possible and reproducible even at low viral load (i.e., 5 copies of viral RNA). Although many methods are available for detecting SARS-CoV-2 virus in tissues [6, 7, 13] , ddPCR accomplishes a level of sensitivity and specificity for SARS-CoV-2 detection unachievable by heavily used RT-qPCR because of inherent technical inconsistencies of this method. For example, PCR inhibitors have a large effect on amplification efficiency, and thus, qPCR assay performance [21, 22] . Concentrating target molecules in smaller reactions, ddPCR can reduce the bias introduced by PCR inhibitors often found in clinical specimens [21] [22] [23] . In addition, the calibration curves needed for qPCR are a major source of assay variability and are difficult to maintain. Consequently, one of the main advantages of ddPCR is the generation of thousands of droplets per sample and this method does not require replicates because it directly counts the number of target molecules rather than relying on a standard curve generated by reference standards or endogenous controls [15] [16] [17] [18] . Therefore, these properties make ddPCR a promising SARS-CoV-2 detection method for samples with low template abundance derived as may be observed in non-respiratory tissues [14] . Since the test became clinically available, it has been ordered for a wide variety of clinical and research purposes by physicians at our institution and across the country. Initially, common indications for clinical testing included evaluating brain or lung tissue, often in the setting of autopsy. Over time, fewer tests were ordered for these indications, and instead testing was used to evaluate lung tissue that showed morphologic features or sequela of an acute or acute on chronic lung injury [24] , skin rashes, gastrointestinal tract ulceration and/or inflammation, and other findings in the setting of current or recent COVID-19. In addition, testing has been used in the transplant setting to test organs from COVID-positive donors, as well as the donated organ after transplant in the setting of recipients with SARS-CoV-2 infection to better understand a patient's clinical course and to distinguish organ damage from SARS-CoV-2 from damage due to rejection. Testing to evaluate placental pathology has also increased over the time the test has been available. Initially testing was often performed to evaluate subtle lesions identified in the placenta, while after the delta variant became wide-spread in the United States, an increase in testing of placenta tissue or products of conception for the indication of fetal demise was noted. Concurrently, although this is not a quantitative test and the raw signal has not been formally evaluated, we noted an increase in the number of placenta samples with a saturated (highly positive) signal after the emergence of the delta variant. In summary, this work demonstrates the utility and sensitivity of tissue-based testing for detection of SARS-CoV-2 in various FFPE tissue types, autopsy to surgical biopsy samples, in a large group of patients and is a feasible method for use in clinical laboratories. The clinical utility is broad with the assay used for applications as diverse as pregnancy loss and transplant rejection. We showed that ddPCR provides both accuracy and sensitivity in detecting patients with low viral load. SARS-CoV-2 ddPCR testing will also aid in understanding atypical presentations of COVID-19, including in patients who present with "unexplained" symptoms or potential long-term sequelae. Funding: This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Tables Table 1 : Detection of SARS-CoV-2 targets at different dilutions tested by RT-ddPCR. Average and range of N1 copies, N2 copies, N1 positive droplets, N2 positive droplets, and total droplets for LOD experiments with known SARS-CoV-2 virion copies, as indicated. SARS-CoV2 vertical transmission with adverse effects on the newborn revealed through integrated immunohistochemical, electron microscopy and molecular analyses of Placenta Postmortem Findings in Italian Patients With COVID-19: A Descriptive Full Autopsy Study of Cases With and Without Comorbidities China Medical Treatment Expert Group for, Clinical Characteristics of Coronavirus Disease 2019 in China Bone Marrow and Peripheral Blood Findings in Patients Infected by SARS-CoV-2 SARS-CoV-2 infection of the placenta Comparison of RNA In Situ Hybridization and Immunohistochemistry Techniques for the Detection and Localization of SARS-CoV-2 in Human Tissues Central nervous system involvement by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Vertical transmission of coronavirus disease 2019: severe acute respiratory syndrome coronavirus 2 RNA on the fetal side of the placenta in pregnancies with coronavirus disease 2019-positive mothers and neonates at birth Pathologic features of COVID-19: A concise review Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Situ Hybridization, Immunohistochemistry, and Reverse Transcription Droplet Digital Polymerase Chain Reaction for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Testing in Tissue Cardiac involvement in the long-term implications of COVID-19 Detection methods for SARS-CoV-2 in tissue ddPCR: a more accurate tool for SARS-CoV-2 detection in low viral load specimens Basic Concepts and Validation of Digital PCR Measurements Evaluation of a droplet digital polymerase chain reaction format for DNA copy number quantification Fundamentals of multiplexing with digital PCR Application of droplet digital PCR to detect the pathogens of infectious diseases Detection of SARS-CoV-2 in formalin-fixed paraffin-embedded tissue sections using commercially available reagents Comparative study of sensitivity, linearity, and resistance to inhibition of digital and nondigital polymerase chain reaction and loop mediated isothermal amplification assays for quantification of human cytomegalovirus Droplet Digital PCR versus qPCR for gene expression analysis with low abundant targets: from variable nonsense to publication quality data Tolerance of droplet-digital PCR vs realtime quantitative PCR to inhibitory substances Late Complications of COVID-19: A Morphologic, Imaging, and Droplet Digital Polymerase Chain Reaction Study of Lung Tissue Negative 0 0 0 0 Negative PA- 22 Negative 0 0 0 0 Negative PA- 23 Negative