key: cord-0829395-64rpg9jg
authors: Han, Wanda G. H.; Swart, Arno; Bonačić Marinović., Axel; Eggink, Dirk; Reimerink, Johan; Wijsman, Lisa A.; van der Veer, Bas; van den Brink, Sharon; van den Brandt, Anne-Marie; van Tol, Sophie; Godeke, Gert-Jan; Brouwer, Fion; Hoogerwerf, Marieke; Reukers, Daphne F. M.; Rots, Nynke; Reusken, Chantal; Meijer, Adam
title: SARS-CoV-2 RNA and antibody dynamics in a Dutch household study with dense sampling frame
date: 2022-05-13
journal: Sci Rep
DOI: 10.1038/s41598-022-11480-6
sha: df81b4ee15e6ab01b995774ec43c2d9747217238
doc_id: 829395
cord_uid: 64rpg9jg

This study investigated the dynamics of SARS-CoV-2 infection and diagnostics in 242 household members of different ages and with different symptom severity after SARS-CoV-2 exposure early in the pandemic (March–April 2020). Households with a SARS-CoV-2 confirmed positive case and at least one child in the Netherlands were followed for 6 weeks. Naso (NP)- and oropharyngeal (OP) swabs, oral fluid and feces specimens were analyzed for SARS-CoV-2 RNA and serum for SARS-CoV-2-specific antibodies. The dynamics of the presence of viral RNA and the serological response was modeled to determine the sampling time-frame and sample type with the highest sensitivity to confirm or reject a SARS-CoV-2 diagnosis. In children higher viral loads compared to adults were detected at symptom onset. Early in infection, higher viral loads were detected in NP and OP specimens, while RNA in especially feces were longer detectable. SARS-CoV-2-specific antibodies have 90% probability of detection from 7 days (total Ig) and 18 days (IgG) since symptom onset. For highest probability of detection in SARS-CoV-2 diagnostics early in infection, RT-PCR on NP and OP specimens are more sensitive than on oral fluid and feces. For SARS-CoV-2 diagnostics late after infection, RT-PCR on feces specimens and serology are more valuable.

www.nature.com/scientificreports/

Study protocol. A prospective cohort study was performed following households where one symptomatic household member was tested RT-PCR positive for SARS-CoV-2 in the period 24 March-6 April 2020 3 . In brief, persons 18 years and older testing positive for SARS-CoV-2 RT-PCR (i.e. the index case) who had at least one child in their household below the age of 18 could be included in this study (METC nr: NL13529.041.06).

The study was carried out according to National guidelines and regulations. The protocol was approved by the Medical-Ethical Review Committee (METC) of the University Medical Center Utrecht. Informed consent was obtained from all subjects and/or their legal guardian(s). A total of 242 participants from 55 complete households were included in this study (See Reukers et al. 3 for more details). The number of analyses performed per assay and specimen type at the various timepoints with the day of the first home visit (so the start of the study within the particular household) defined as day 1 are described in Table 1 and Table S1 . We adopted a different age classification than Reukers et al., foregoing the 'adolescent' category, and defining adults as individuals of 18 years of age or older. Individuals are considered SARS-CoV-2 infection positive when they tested positive in at least one RT-PCR or serological assay.

COVID-19 severity. The day of onset of possible COVID-19 associated symptoms, i.e. respiratory symptoms (including sore throat, cough, dyspnea or other respiratory difficulties, rhinorrhoea), fever, chills, headache, anosmia or ageusia, muscle pain, joint ache, diarrhoea, nausea, vomiting, loss of appetite or fatigue, as reported by the participant was defined as 1 dps. Severity, independent of being infected by SARS-CoV-2 or not, was defined based on self-reported clinical symptoms, consultation of a medical specialist and being hospitalized. Cases with any clinical symptoms other than dyspnea (used as proxy for possible pneumonia) and not reporting consultation of a medical specialist or being hospitalized were defined as mild cases. Moderate cases reported dyspnea and did not report consultation of a medical specialist or hospitalization. Severe cases reported dyspnea and consulted a medical specialist for their symptoms and/or reported having been admitted to the hospital 3, 4 . Severity was determined at visits 1, 2, 3, and overall, where the status at visit 1 included also reports on the 14 days before the visit.

Transmission categories. We categorized the household transmission patterns in three groups. In the 'no transmission' category, SARS-CoV-2 infection was only detected in the index case. In the 'adult transmission' category, SARS-CoV-2 infection was detected in adults other than the index case only. In the 'family transmission' category, SARS-CoV-2 infection was detected in children and possibly also adults other than the index case. In order to assess the significance of differences in severities over transmission categories a generalized linear model (GLM) for count data with Poisson family and exponential link was set up with interactions between the transmission category and severity classes.

Molecular diagnostics. Nasopharyngeal swabs (NP) and oropharyngeal swabs (OP) were collected in gelatin-lactalbumin-yeast (GLY) viral transport medium (Mediaproducts BV, Groningen, The Netherlands), transported to the laboratory in a cooling box and stored at maximum three days at 4 °C until being processed for RT-PCR. Feces specimens were self-collected by the patient and send to the laboratory by regular mail, stored frozen at − 20 °C until being processed for RT-PCR. Oral fluid specimens were collected with a Oracol sponge (Malvern Medical Developments Ltd, U.K.), transported to the laboratory in a cooling box, processed for storage according to the manufacturer's instructions, and aliquots stored frozen at − 80 °C until being used for RT-PCR. Total nucleic acid was extracted from NP swabs, OP swabs, oral fluid or feces using MagNApure 96 (MP96) with total nucleic acid kit small volume (Roche). Of the feces specimens a 5% suspension was made in Minimal Essential Medium (MEM) with Hanks' salts and penicillin and streptomycin, vortexed for 15 secs and www.nature.com/scientificreports/ centrifuged for 1 min at 16,000 Relative Centrifugal Force. Two-hundred µl supernatant was mixed with 275 µl MP96 lysis buffer including equine arteritis virus (EAV) internal control and yeast tRNA stabilizer. Total nucleic acid was eluted in 50 µl Tris EDTA buffer. RT-qPCR was performed on 5 µl total nucleic acid using TaqMan® Fast Virus 1-Step Master Mix (Thermo Fisher) on Roche LC480 II thermal cycler with SARS-like beta coronavirus (Sarbeco) specific E-gene primers and probe and EAV as described previously 5, 6 . As no other Sarbeco viruses are currently detected in humans, a positive Sarbeco E-gene RT-qPCR is validly taken as positive for SARS-CoV-2. For modeling purposes no detection of SARS-CoV-2 RNA was given an artificial cycle threshold (Ct) value of 40.

Serological diagnostics. The Wantai SARS-CoV-2 total antibody ELISA (Beijing Wantai Biological Pharmacy Enterprise, Beijing, China; catalogue number WS1096) was performed according to the manufacturer's instructions 7 . This assay is a double-antigen sandwich ELISA using the recombinant receptor-binding domain of SARS-CoV-2 as antigen. Optical density (OD) is measured at 450 nm and the antibody OD ratio for each sample is calculated as the ratio of the OD of that sample to the reading of a calibrator (included in the kit). Sera with a ratio ≥ 1 were considered positive for antibodies against SARS-CoV-2. Sera were tested for the presence of IgG antibodies reactive with the SARS-CoV-2 S1 and SARS-CoV-2N antigens in a protein microarray, in duplicate twofold serial dilutions starting at 1:20, essentially as described previously 8 . For each antigen, a 4-parameter loglogistic calibration curve was generated. Antibody titers (EC50 value) were defined as the interpolated serum dilution that gave a fluorescence intensity of 50% of the corresponding calibration curve. Therefore, sera testing negative in the assay are given a titer of 10 (20 dvided by 2). Sera with a titer > 10 were considered positive for antibodies against SARS-CoV-2. Raw data were processed with the R 4.04 statistical software as described previously 9 .

Modeling RT-PCR. All available RT-PCR outcomes ( Table 1) were modelled by a Bayesian hierarchical model of the form Here, Ct i is the measured Ct value for sample i , and σ the overall variation. There is a part dependent on the days since onset of symptoms ( d i ) with coefficients β , and a constant part with coefficients α . Both the α and β parameters include several contributions stratified by categorical variables: id[i] is the person specific identifier for sample i enabling longitudinal modelling, agecat[i] is either 'child' or 'adult' , severity[i] is either 'asymptomatic' , 'mild' , 'moderate' or 'severe' , and test[i] is either 'pcr NP, 'pcr OP, 'pcr oral fluid' , 'pcr feces' , 'wantai' , 'microarray S1' or 'microarray N' .

For the id's a hierarchical model is built (i.e. a random effect), Whenever a Ct-value of 40 is encountered, we apply censoring in the model by changing the probability density function f (x; µ i , σ ) into its cumulative counterpart 1 − F(40; µ i , σ ) , thereby encoding that we have an unknown Ct value which would either indicate the absence of amplifiable RNA or presence of RNA but well below the detection limit of the used RT-PCR. All parameters are given weakly informative priors, and the posterior distributions are obtained using the JAGS software 10 , interfaced from R 11 . Bayesian credible intervals were obtained from the samples of the posterior as calculated by JAGS. Prediction intervals were calculated by drawing randomly from N µ j ,σ j , where the indicated means and standard deviations are samples from the posterior distributions. The posterior probability of being positive is modelled by F 40;μ j ,σ j . The modeling did not include the Ct values of the inclusion RT-PCR NP + OP which was performed in index cases just before start of the study, since these values are not known to us.

Modeling serology. The dynamics of serology cannot be assumed to be linear as is the case for Ct-values.

Rather, seronegative individuals have a titer (OD ratio for Wantai or EC50 for protein microarray) varying around a low value, and seropositive individuals have a titer varying around a high value. In the case of the ELISA-test and microarray-based assays used in the current study, we find that a cut-off value to distinguish seropositives and seronegatives works well, since the two components are well separated (Fig. S1 ). Using the cut-off values 1 for Wantai (according to manufacturer's instructions) and 10 for microarray 8 , we classify each measurement X i as either positive or negative. Using a Bernouilli distribution and logit link for the probability we model the outcomes as Parameter estimation proceeds analogous to the RT-PCR model. www.nature.com/scientificreports/ Assessing differences between factors. We assessed the difference between posterior estimates of parameters using the Region of Practical Equivalence (ROPE) [12] [13] [14] 

Household transmission SARS-CoV-2. To identify different transmission patterns, we visualized SARS-CoV-2 infection detection by the different assays and specimen types per participant and household in heatmaps. We identified the transmission pattern 'no transmission' in 16 households (Fig. 1A) , 'adult transmission' in 11 households (Fig. 1B) and 'family transmission' in 28 households (Fig. 1C) . Eight of the 28 households in the 'family transmission' category did not show transmission to adults. Symptom severity of COVID-19 index cases did not correlate with transmission of SARS-CoV-2 within a household (Fig. 2) . In the no transmission category, more than half (56.3%) of the index cases had mild symptoms, whilst only 12.5% had severe symptoms. In the adult-and family transmission category, there were 3 out of 10 (30%) and 8 out of 28 (28.6%) severe index cases respectively, but these results were not significant.

We investigated the SARS-CoV-2 infection dynamics in the participants of the study. Using an 'upset plot' 15 , Fig. 3 shows patterns of positive and negative results in the various molecular and serological assays. Of the 242 participants, 136 individuals were positive for SARS-CoV-2 infection by either molecular and/or serological diagnostics. Of these 136 SARS-CoV-2 infection positive individuals, 125 (91.9%) were symptomatic and all severe symptomatic individuals were RT-PCR positive. Only one severly ill individual was serology-negative (Fig. 3B ). Most individuals were found SARS-CoV-2 positive by multiple diagnostic assays and/or materials, but 19 individuals tested positive with only one assay type and/or material (only Wantai n = 5, RT-PCR NP n = 3, RT-PCR oral fluid n = 3, RT-PCR OP n = 2, RT-PCR feces n = 2, Inclusion PCR NP + OP n = 2, MA-S n = 1, MA-N n = 1) during the study period (Fig. 3A) . In most of the infected cases (81.6%, 111 of the 136) both SARS-CoV-2 RNA and SARS-CoV-2-specific antibodies were detected (Fig. 3B ). For 12 individuals only one positive test was found at one timepoint ( Fig. 1, red rectangles) .

Next, for the individuals with at least one RT-PCR and one serological result at visit 1, 2 and 3 (n = 198), we could analyze rough dynamics of the infection process (Fig. 4) . The median dps relative to visit 1 is indicated in Fig. 4 . Six common patterns in 173 individuals, ranked A-F based on frequency, could be identified. Within laboratory confirmed SARS-CoV-2 infected individuals, the common patterns B (n = 28) and C (n = 27) included individuals with a positive PCR and serological assay at visit 1. As can be expected, SARS-CoV-2 viral RNA was not detected anymore in these cases at the end of the study (4-6 weeks after inclusion), while SARS-CoV-2 antibodies remained present. Individuals with pattern E (n = 12) did not have detectable SARS-CoV-2 RNA at any visit, but did have detectable SARS-CoV-2 antibodies at visit 1, 2 and 3. This pattern is in line with an earlier onset of symptoms compared to the individuals with pattern B and C, thus these individuals were included in the study later in their infection process resulting in already diminished viral RNA and present antibodies at visit 1. Pattern D (n = 18) and F (n = 9) included individuals with a positive RT-PCR at visit 1 and developed antibodies after visit 1. Compared to pattern B and C, these individuals reported their onset of symptoms 2-4 days later, thus at study inclusion (visit 1) they were earlier in their infection process. Pattern A included individuals with negative RT-PCR and serology results at all visits. These rough patterns underline that there are optimal time windows in which detection of SARS-CoV-2 viral RNA or SARS-CoV-2-specific antibodies are most appropriate in diagnostics. To further investigate this we used a modelling approach.

Dynamics SARS-CoV-2 diagnostics. Bayesian modeling on all available RT-PCR data demonstrated a difference in Ct-values at symptoms onset between adults and children (Fig. 5A, B and S2-G1) and in Ct-values at symptom onset and in Ct-value increase per day between different specimen types ( Fig. 5C and S2-A-F) . The predicted Ct value (inversely correlated with viral load), was on average 2.6 Ct lower in children (Ct 27.5; all specimens) compared to adults (Ct 30.1; all specimens) at the day of symptom onset (intercept) (Fig. 5A, B and S2-G1). The decay over time in viral load (slope) was comparable between adults and children ( Fig. 5B and S2-G2 ). In line with this, there is a longer probability of SARS-CoV-2 RNA detection with increasing dps in children (99% detection until 13 dps), compared to adults (99% detection until 7.6 dps) (Fig. S3A) . When analyzing all ages, the predicted viral load seems slightly higher in NP and OP swabs (Ct 28.9 and 28.8) compared to feces (Ct 30.1) and oral fluid (Ct 30.8) at the day of symptom onset (intercept) (Fig. 5C) indicating higher sensitivity of SARS-CoV-2 RNA detection in NP and OP specimens compared to oral fluid (Fig. S2-B1 and D1) and possibly feces (Fig. S2-C1 and E1 ). In contrast, there seems te be a slower decay (slope) in viral load in oral fluid and feces specimens (0.25 and 0.22 Ct per day) compared to NP and OP specimens (0.35 and 0.36 Ct per day) (Fig. 5C) . The relevance of these findings is uncertain as there is partial overlap between the ROPE and HDI (Fig. S2-B1-E1 and B2-E2) . Overall, as a most likely estimate, the estimated viral load is lower in oral fuid compared to NP swabs (until 21.1 dps) and, OP swabs (until 19.0 dps) (Fig. S4-B and D) . Also, feces samples have lower viral load than NP (until 10.2 dps) and OP (until 9.7 dps) (Fig. S4-C and E) . However, there is a longer probability of SARS-CoV-2 RNA detection with increasing dps in feces (90% detection until 27 dps), compared to NP, OP and www.nature.com/scientificreports/ revealed that there were less individuals in the "Severe" category than in the "Mild" category (p = 0.1), and that households in "Family transmission" category was overrepresented in the "Severe" category (p = 0.03). There were no asymptomatic index cases. Maximum severity score of index is used. 5D and Figure S5 -C, E and F). Similar trends are shown for 50% and 10% detection probability (Figs. S5 and S6). We could not find a clear correlation in severity of symptoms and the dynamics of the SARS-CoV-2 RNA detection (Fig. S10) . Furthermore, we investigated the dynamics of SARS-CoV-2 antibody detection (Fig. 6 ). The Wantai assay (total Ig) demonstrated a higher sensitivity for detection of anti-S1 antibodies than the micro-array (IgG) as the probability for detection was earlier using Wantai upon onset of illness ( Fig. 6A and S7-A) . The dps at which 90% detection probability was reached for Wantai was 7.1 compared to 16.9 and 18 for Nucleoprotein (N)-and S1-protein microarray respectively (Fig. 6B) . The protein microarray for S1 and N had comparable sensitivity, in line with a previous study 8 . The probability of detecting N-specific IgG antibodies in children was delayed by 3.0 days (at 90% probability detection) versus adults (Fig. 6C,D) , while this was not the case for detection of S1-specific antibodies (Fig. 6C) . The relevance of this finding is uncertain as there is much overlap between the distributions in children and adults (Fig. S8) . The N-specific IgG antibody titers (at visit 3; convalescent phase) were not significantly lower in children compared to adults (Fig. S9) . Furthermore, we could not find a correlation in severity of symptoms and the dynamics of the SARS-CoV-2 infection detection by the serological assays (Fig. S11 ).

We studied SARS-CoV-2 RNA and antibody kinetics in a household cohort during the early phase of the pandemic using an unusually dense sampling schedule allowing for high resolution analysis. For highest probability of detection of SARS-CoV-2 RNA by RT-PCR early in infection NP and OP are more suitable than oral fluid and feces. SARS-CoV-2-specific antibodies have a 90% probability of detection from 7 dps with the Wantai assay and 18 dps with the microarray S1 and N assay. This study has been performed in a naïve population during the early phase of the pandemic. Although SARS-CoV-2 vaccination, previous infection and the circulation of other SARS-CoV-2 variants may influence the dynamics of SARS-CoV-2 infection and thereby diagnostics, our study provides valuable reference insights into this subject.

Households present close-contact settings with high risk of SARS-CoV-2 transmission after introduction of SARS-CoV-2 in the household 3, 16 . We observed a positive correlation, although not significant, between transmission of SARS-CoV-2 between adults and children and the severity of disease in the household indexes. This is in line with studies that report that the severity of SARS-CoV-2 infection of the index case was associated with higher infectiousness 16, 17 . It should be noted that at the time of the study, SARS-CoV-2 testing in the Netherlands was limited to symptomatic healthcare workers and symptomatic vulnerable individuals. The index cases were thus mainly symptomatic healthcare workers 3 . As schools and daycare centres were closed during the study period, transmission outside the households among children was minimalized. Due to limited sample size, we could only categorize children as those of 17 years of age or younger for the Bayesian hierarchical modelling. It would be of interest to stratify the children in more groups, e.g. primary school age and adolescent age. www.nature.com/scientificreports/ Furthermore, we did not have data on compliance to the governmental strongly advised self-isolation protocol 3 .

Breach of compliance could have had an influence on household-transmission. Seroconversion rates in mild to severe symptomatic SARS-CoV-2 RT-PCR positive cases have been reported in the range of 93-100% after 3-4 weeks [18] [19] [20] . During the course of this study, a vast majority of RT-PCR confirmed SARS-CoV-2 infected participants (91.0%, 111/122) developed SARS-CoV-2 specific antibodies (Fig. 3B) . In 6 RT-PCR positive cases, serology data was missing. In 11 of 122 (9.0%) RT-PCR confirmed SARS-CoV-2 infected individuals there was no seroconversion (Fig. 3B) . These individuals may have experienced a relatively mild infection, or in the cases with only one positive RT-PCR test (including the available results of the extra sampling between visit 1 and 2, n = 8) and no evidence of antibody response, the SARS-CoV-2 diagnosis of these individuals is disputable. Four of the 8 RT-PCR-negative cases that did have SARS-CoV-2-specific antibodies In the current study we used Bayesian hierarchical models to determine the sampling time-frame and sample type with the highest sensitivity to confirm or reject a SARS-CoV-2 diagnosis. Due to the close-contact settings with high probability of (pre/a-symptomatic) transmission of SARS-CoV-2 after introduction of the virus in households and the longitudinal dense sampling performed in this study, the analysis yielded results with Figure 6 . Dynamics of SARS-CoV-2 infection diagnosis by different serological assays. All available serology outcomes (Table 1) were incuded in these analyses. Probability of SARS-CoV-2-specific antibody detection (A) and average dps from when the Wantai (Spike-specific IgM and IgG), microarray S1 (Spike-specific IgG) and microarray N (Nucleoprotein-specific IgG) assays have at least 90% detection probability (B) for all ages. (C) Probability of SARS-CoV-2-specific antibody detection by Wantai, microarray S1 and microarray N in adults and children (D) Average dps from when microarray N has at least 90% detection probability for adults and children. The shadows in (C) and (D) indicate the 95% Bayesian confidence interval. www.nature.com/scientificreports/ relatively narrow credibility intervals which support our findings. For highest probability of a SARS-CoV-2 diagnosis (lower Ct), RT-PCR on NP and OP specimens were more suitable than feces and oral fluid until 10 dps and 20 dps, respectively (Fig. S4C,E and B,D) . Our study confirms that NP and OP (or combined NP + OP) were the preferred sample type for RT-PCR-based SARS-CoV-2 diagnostics within 1 week upon onset of symptoms [22] [23] [24] .

A systematic review concluded that of the alternative specimens to NP and OP swabs, oral fluid (saliva) has an estimated sensitivity of 83.9% (95% CI 77.4-88.8) and specificity of 96.4% (95% CI 89.5-98.8) compared to reference NP and OP swabs in nucleic acid assays 24 . The sensitivity and specificity of feces specimens seems much lower, although limited data is available 24 . For SARS-CoV-2 diagnostics late in infection or in past infections, RT-PCR on feces and oral fluid specimens are more valuable than NP and OP specimens, since viral RNA in especially feces remains present over a longer time compared to NP and OP swab specimens. This is in line with findings of other studies that indicate that SARS-CoV-2 RNA can be detected up to 126 days in feces compared to 83 days in respiratory specimens and that beyond 10 dps, feces sampling may be preferred 23, 25, 26 .

Although SARS-CoV-2 RNA can remain present in respiratory and feces specimens for a long time, the duration of presence of viable virus is relatively short-lived 27, 28 . However, this pattern might change depending on the characteristics of emerging virus variants 29 and the presence of pre-existing antibodies due to experienced infection in the past or vaccination 30 . Alternatively, from a week after symptom onset testing for the presence of SARS-CoV-2 Spike-specific total Ig using Wantai ELISA can confirm a recent or past SARS-CoV-2 infection. SARS-CoV-2 diagnostics using protein microarray detecting SARS-CoV-2 Spike-and Nucleoprotein-specific IgG antibodies, is useful 2 weeks after infection or symptom onset. The infection dynamics of SARS-CoV-2 may be influenced by characteristics of the tested population, such as age and the severity of COVID-19. We, however, could not find a clear correlation in severity of symptoms or age and the RNA and SARS-CoV-2-specific antibody kinetics, although our study may be underpowered to detect these differences (data not shown). Children in general report milder symptoms compared to adults (Table S1) 3 . It is known that with age, the expression of ACE-2 increases in nasal epithelium 31 . Since SARS-CoV-2 uses the ACE-2 receptor for host entry, a lower expression of ACE-2 in children relative to adults might explain the lower susceptibility and milder infection course in children. However, our findings and other studies show that viral loads in children are similar or higher than viral loads in adults 32, 33 . In our study children displayed lower Ct values (higher viral loads) at the day of symptom onset compared to adults, while the decay in viral load was comparable (Fig. 5 ). This suggests that if children become infected with SARS-CoV-2, they can carry high loads of virus for a longer time compared to adults. Therefore, children are potentially longer infectious than adults after symptom onset. Whether this observation holds for new Variants of Concerns (VOCs) e.g. delta or omicron, warrants further investigation. No clear differences between adults and children were found in the dynamics of SARS-CoV-2 serology, yet the detection of N-specific antibodies seems slightly delayed in children compared to adults (Fig. 6C,D and S8) . A study showed a reduced breadth of anti-SARS-CoV-2-specific antibodies, predominantly generating IgG antibodies specific for the S protein but not the N protein in children compared to adults 34 . Whether this has consequences for the development of immunity to SARS-CoV-2 is not yet clear.

A drawback of the curent study design is that infections originating from outside the household cannot be excluded. This could be tested by submitting the samples to molecular typing, and comparing the sequences to community variants, as performed in 35 . Another improvement that could be implemented in future househould transmission studies is a denser temporal sampling frame. Having more time-points would allow for capturing the initial rising phase of antibody levels and viral load, and also have more precise estimates of the waning of antibodies and decline of viral loads. This in turn would provide valuable information on transmissibility and detectability of SARS-CoV-2.

In summary, our study allowed for a high resolution analysis of the sensitivity of molecular and serologybased detection of recent SARS-CoV-2 infections due to the unusually dense sampling strategy in a confined setting. For highest probability of SARS-CoV-2 diagnostics early in infection, PCR on NP and OP specimens are in favor over oral fluid and feces. For SARS-CoV-2 diagnostics late in infection or in past infection, RT-PCR on feces specimens and serology are more valuable. Children seem to carry higher loads of virus for a prolonged time in comparison to adults. The data presented here strengthen the evidence-basis for SARS-CoV-2 testing strategies. 

Timeline of WHO's response to COVID-19

Minister Tamara van Ark. Kamerbrief met nieuwe aanvullende maatregelen om de COVID 19 uitbraak te bestrijden. Medische Zorg en Sport OCeW, Basis-en Voortgezet Onderwijs en Media

High infection secondary attack rates of SARS-CoV-2 in Dutch households revealed by dense sampling

Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR

Diagnosis of human metapneumovirus and rhinovirus in patients with respiratory tract infections by an internally controlled multiplex real-time RNA PCR

Evaluation of nine commercial SARS-CoV-2 immunoassays

Accurate serology for SARS-CoV-2 and common human coronaviruses using a multiplex approach

Profiling of humoral immune responses to influenza viruses by using protein microarray

JAGS: A program for analysis of bayesian graphical models using gibbs sampling

R: A Language and Environment for Statistical Computing. Foundation for Statistical Computing

The Bayesian new statistics: Hypothesis testing, estimation, meta-analysis, and power analysis from a Bayesian perspective

Statistical rethinking: A Bayesian course with examples in R and STAN

Comparison of Bayesian predictive methods for model selection

UpSet: Visualization of Intersecting Sets

Household transmission of SARS-CoV-2: A systematic review and meta-analysis

The household secondary attack rate of SARS-CoV-2: A rapid review

Dynamics and correlation among viral positivity, seroconversion, and disease severity in COVID-19: A retrospective study

Antibody responses to SARS-CoV-2 in patients with COVID-19

Humoral response and PCR positivity in patients with COVID-19 in the New York City region, USA: An observational study

A case of transient existence of SARS-CoV-2 RNA in the respiratory tract with the absence of anti-SARS-CoV-2 antibody response

Diagnosis of SARS-Cov-2 infection by RT-PCR using specimens other than naso-and oropharyngeal swabs: A systematic review and meta-analysis

Towards a sensitive and accurate interpretation of molecular testing for SARS-CoV-2: A rapid review of 264 studies

Performance of Saliva, Oropharyngeal Swabs, and Nasal Swabs for SARS-CoV-2 Molecular Detection: a Systematic Review and Meta-analysis

At what times during infection is SARS-CoV-2 detectable and no longer detectable using RT-PCR-based tests? A systematic review of individual participant data

SARS-CoV-2, SARS-CoV, and MERS-CoV viral load dynamics, duration of viral shedding, and infectiousness: A systematic review and meta-analysis

Duration and key determinants of infectious virus shedding in hospitalized patients with coronavirus disease-2019 (COVID-19)

Virological assessment of hospitalized patients with COVID-2019

Viral dynamics and duration of PCR positivity of the SARS-CoV-2 Omicron variant

Viral dynamics of SARS-CoV-2 variants in vaccinated and unvaccinated persons

Nasal gene expression of angiotensin-converting enzyme 2 in children and adults

Age-related differences in nasopharyngeal severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) levels in patients with mild to moderate coronavirus disease 2019 (COVID-19)

Culture-competent SARS-CoV-2 in nasopharynx of symptomatic neonates, children, and adolescents

Distinct antibody responses to SARS-CoV-2 in children and adults across the COVID-19 clinical spectrum

Virologic features of severe acute respiratory syndrome coronavirus 2 infection in children

We thank the Public Health Service Utrecht for assistance in the recruitment of households. We thank Alper 

This study was funded by the Dutch Ministry of Health, Welfare, and Sport (VWS).

The authors declare no competing interests.