Introduction

Coronavirus Disease 2019 (COVID-19) vaccines have demonstrated protection against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)1,2,3. However, the efficacy against acquisition of infection is imperfect in that vaccines have an efficacy VES < 100%, particularly against variants of concern4,5,6. Breakthrough infections in vaccinated individuals have been documented4,5,6,7,8,9, but the transmission potential of these infections is poorly understood. It is conceivable that vaccinated persons who acquire the infection may be less infectious than unvaccinated persons who acquire the infection, as the vaccine-primed immune response may attenuate the natural history of infection, by reducing viral replication and accelerating viral clearance, leading to lower viral load and shorter duration of infection10,11. There is evidence that seems to support this hypothesis for SARS-CoV-2 infection12,13,14. Therefore, COVID-19 vaccines may be not only efficacious against acquisition of infection (VES), but also against transmission of infection10,15,16, thereby adding an additional efficacy for each vaccine, denoted as VEI and defined as the proportional reduction in infectiousness among those infected but vaccinated, compared to those infected but unvaccinated10.

Leveraging the national, federated databases that have captured all SARS-CoV-2 vaccinations and polymerase chain reaction (PCR) testing in Qatar since the start of the epidemic (“Methods”), we investigated the effect of vaccination on infectiousness by comparing SARS-CoV-2 real-time (quantitative) reverse transcription-PCR (RT-qPCR) cycle threshold (Ct) values of individuals infected and fully vaccinated with the values of those infected and unvaccinated. The RT-qPCR Ct value is a measure of the inverse of viral load and correlates strongly with culturable virus;17 thus, it can be used as a proxy of SARS-CoV-2 infectiousness17,18,19,20,21. We also investigated the effect of prior infection on infectiousness at reinfection by comparing the RT-qPCR Ct values among those reinfected with SARS-CoV-2 with the values among those with primary infection. For standardization of Ct values, the comparisons were conducted using only the RT-qPCR-confirmed infections diagnosed using the TaqPath COVID-19 Combo Kits (Thermo Fisher Scientific, USA22) which were used for >85% of all RT-qPCR tests in Qatar4,5,6,22.

The comparisons were implemented utilizing: (i) the national cohort of all 384,452 RT-qPCR-confirmed primary infections since epidemic onset (February 28, 2020) until the end of the study (July 11, 2021; Fig. 1a); (ii) the national cohort of all 1695 RT-qPCR-confirmed reinfections during the same period (Fig. 1b); (iii) the national cohort of all 898,648 individuals vaccinated with the BNT162b21 (Pfizer-BioNTech) vaccine since the first recorded vaccination in Qatar on December 16, 2020 until the end of the study (July 11, 2021; Fig. 1c); and (iv) the national cohort of all 468,872 vaccinated individuals with the mRNA-12732 (Moderna) vaccine during the same period (Fig. 1d).

Fig. 1: Flowchart illustrating the population selection process.
figure 1

Flowchart illustrating the selection of the cohorts of (a) primary infections in unvaccinated individuals, (b) reinfections in unvaccinated individuals, (c) BNT162b2-vaccine breakthrough infections, and (d) mRNA-1273-vaccine breakthrough infections.

Full size image

The BNT162b2 and mRNA-1273 vaccines have been the vaccines of choice in the national immunization campaign in Qatar4,5,6,8,9,23,24. In total, there has been 4777 breakthrough infections in those fully vaccinated with BNT162b2 (0.61% of those fully vaccinated; Fig. 1c) and 306 mRNA-1273 breakthrough infections in those fully vaccinated with mRNA-1273 (0.09% of those fully vaccinated; Fig. 1d). Mass immunization started with the BNT162b2 vaccine, and the mRNA-1273 vaccine was introduced only several weeks later. During this vaccination drive, Qatar experienced two back-to-back epidemic waves dominated by the Alpha25 (B.1.1.7) and Beta25 (B.1.351) variants, with the Beta wave contributing to most of the diagnosed infections during 2021 (Supplementary Fig. 1)4,5,6,8,9,26,27,28. While incidence since July 2021 has been dominated by the Delta25 (B.1.617.2) variant4,5,6,8,9,26,27,28, the number of Delta infections remained relatively low with no epidemic wave materializing as of the end date of this study, July 11, 2021.

Primary infection was defined as the first RT-qPCR-positive test for a given individual. Reinfection was defined as the first RT-qPCR-positive test that occurred ≥90 days after the primary infection29,30,31,32,33. Breakthrough infection in a vaccinated individual was defined as an RT-qPCR-positive test 14 or more days after the individual received the second vaccine dose, conditional on this RT-qPCR-positive test being the first ever positive for this individual.

Results

Study populations

Table 1 shows the Ct value of the RT-qPCR-positive tests stratified by reason for testing for all SARS-CoV-2 infections diagnosed in Qatar using the TaqPath COVID-19 Combo Kits platform from February 28, 2020 to July 11, 2021. Among symptomatic infections, the mean Ct value was 23.1 (95% CI: 23.0–23.1). Among asymptomatic infections, the mean Ct value was 25.2 (95% CI: 25.1–25.2).

Table 1 Characteristics of RT-qPCR Ct values in the study population.
Full size table

Figure 1 shows the process for identifying eligible primary infections, reinfections, BNT162b2 breakthrough infections, and mRNA-1273 breakthrough infections. Figure 2 schematizes the six pairwise comparisons conducted between these cohorts of infections, after exact matching to the first match in a 1:1 ratio by sex, 10-year age group, reason for RT-qPCR testing, and calendar week of the RT-qPCR test, to control for differences in biology by sex and age, as well as risk of exposure to SARS-CoV-2 in Qatar34,35,36,37,38 and variant exposure4,5,6,8,9,26,27,28 (“Methods”). Only the first RT-qPCR-positive test for each individual for each infection category was used for analysis.

Fig. 2: Schematic diagram showing the process of formulating the six pairwise comparisons between study cohorts.
figure 2

Schematic diagram showing the process of formulating the six pairwise comparisons between the cohorts of primary infections in unvaccinated individuals, reinfections in unvaccinated individuals, BNT162b2-vaccine breakthrough infections, and mRNA-1273-vaccine breakthrough infections, after 1:1 matching by sex, 10-year age group, reason for RT-qPCR testing, and calendar week of the RT-qPCR test.

Full size image

Supplementary Tables 1 and 2 show demographic characteristics of the study populations (age, sex, and nationality) in the six pairwise comparisons. Scatter plots of the distribution of RT-qPCR Ct values in each comparison were generated for all RT-qPCR-confirmed infections, regardless of the reason for the RT-qPCR testing (Supplementary Fig. 2), for only the randomly diagnosed (asymptomatic) infections (Fig. 3; “Methods”), and for only the symptomatic infections (Supplementary Fig. 3; “Methods”).

Fig. 3: RT-qPCR Ct values in the randomly diagnosed (asymptomatic) SARS-CoV-2 infections.
figure 3

Distribution of these Ct values (blue circles) in the six pairwise comparisons between primary infections in unvaccinated individuals, reinfections in unvaccinated individuals, BNT162b2-vaccine breakthrough infections, and mRNA-1273-vaccine breakthrough infections, a–f. A randomly diagnosed infection was defined as an RT-qPCR-positive test conducted with no prior reason to suspect infection and no reported presence of symptoms compatible with a respiratory tract infection. That is, the RT-qPCR test was conducted as part of a survey (random testing campaigns), for routine healthcare testing, for pre-travel requirement, or at port of entry upon arrival in Qatar. a includes, in each comparison group, n = 1584 biologically independent samples, b includes n = 158 biologically independent samples, c includes n = 140 biologically independent samples, d includes n = 987 biologically independent samples, e includes n = 421 biologically independent samples, and f includes n = 60 biologically independent samples, each over 1 experiment. Boxplots center lines indicate the median Ct values, box limits indicate the 25% and 75% quartiles, and whiskers indicate maximum and minimum observations within 1.5 of interquartile range. Paired t tests were used to compare the difference in means between study groups, with no adjustment for multiple comparisons. Two-sided p-values are reported.

Full size image

Differences in RT-qPCR Ct values in all confirmed infections

In the comparisons including all RT-qPCR-confirmed infections, regardless of the reason for the RT-qPCR testing, the mean RT-qPCR Ct value was higher in all cohorts of breakthrough infections compared to the cohort of primary infections in unvaccinated individuals (Table 2 and Supplementary Fig. 2). The Ct value was 1.0 (95% confidence interval (CI): 0.7–1.2) cycles higher for BNT162b2 breakthrough infections, 3.5 (95% CI: 2.5–4.5) cycles higher for mRNA-1273 breakthrough infections, and 3.8 (95% CI: 3.4–4.2) cycles higher for reinfections in unvaccinated individuals. The Ct value was 3.1 (95% CI: 2.0–4.2) cycles higher in mRNA-1273 breakthrough infections than in BNT162b2 breakthrough infections. Compared to reinfections in unvaccinated individuals, the Ct value was 2.1 (95% CI: 0.33–3.8) cycles lower in mRNA-1273 breakthrough infections and 1.7 (95% CI: 1.1–2.3) cycles lower in BNT162b2 breakthrough infections. All differences in Ct values were statistically significant with p values ≤ 0.02.

Table 2 RT-qPCR Ct values in the six pairwise comparisons including all infections.
Full size table

To account for the possibility of prolonged primary infections39,40,41, a comparison was restricted to reinfections in which the RT-qPCR-positive test occurred between 90 and 180 days after the primary infection. Compared to primary infections in unvaccinated individuals, the Ct value was 6.2 (95% CI: 5.4–6.9) cycles higher for the reinfections in unvaccinated individuals (p values ≤ 0.001).

The Ct values across the different infection categories were analyzed using matched pairwise comparisons, rather than using an overall regression model, for assessment of pairwise differences, improved statistical precision of the differences, and better control of confounding factors, considering the relatively small number of breakthrough infections in some of the categories. For additional validation of the results, all RT-qPCR-confirmed infections were also analyzed using an overall regression model (Table 3). The results of this analysis were similar and consistent with the results of the main analysis using matched pairwise comparisons.

Table 3 Linear regression analyses for the RT-qPCR Ct value across the different categories of infections, adjusting for sex, age, nationality, reason for RT-qPCR testing, and calendar week of RT-qPCR testa.
Full size table

Differences in RT-qPCR Ct values in randomly diagnosed infections

In the comparisons including only the randomly diagnosed (asymptomatic) infections (“Methods”), which are perhaps most representative for differences between these cohorts of infections, given the random diagnosis, the mean RT-qPCR Ct value was also higher in all cohorts of breakthrough infections compared to the cohort of primary infections in unvaccinated individuals (Table 4 and Fig. 3). The Ct value was 1.3 (95% CI: 0.9–1.8) cycles higher for BNT162b2 breakthrough infections, 3.2 (95% CI: 1.9–4.5) cycles higher for mRNA-1273 breakthrough infections, and 4.0 (95% CI: 3.5–4.5) cycles higher for reinfections in unvaccinated individuals. The Ct value was 2.2 (95% CI: 0.9–3.6) cycles higher in mRNA-1273 breakthrough infections than in BNT162b2 breakthrough infections. Compared to reinfections in unvaccinated individuals, the Ct value was 2.0 (95% CI: 1.2–2.8) cycles lower in BNT162b2 breakthrough infections and 1.7 (95% CI: −0.3 to 3.7) cycles lower in mRNA-1273 breakthrough infections. All differences in Ct values were statistically significant with p values ≤ 0.003 except for the comparison between reinfections and mRNA-1273 breakthrough infections; likely a consequence of the small number of breakthrough infections that have been documented among those vaccinated with the mRNA-1273 vaccine in Qatar.

Table 4 RT-qPCR Ct values including only the randomly diagnosed (asymptomatic) infections in the six pairwise comparisons.
Full size table

Differences in RT-qPCR Ct values in symptomatic infections

In the comparisons including only the symptomatic infections (“Methods”), the mean RT-qPCR Ct value was also higher in all cohorts of breakthrough infections compared to the cohort of primary infections in unvaccinated individuals, but the difference was smaller for BNT162b2 breakthrough infections (Table 5 and Supplementary Fig. 3). The Ct value was 0.2 (95% CI: −0.2 to 0.6) cycles higher for BNT162b2 breakthrough infections, 4.9 (95% CI: 2.3–7.4) cycles higher for mRNA-1273 breakthrough infections, and 3.8 (95% CI: 2.9–4.7) cycles higher for reinfections in unvaccinated individuals. The Ct value was 5.3 (95% CI: 2.6–8.1) cycles higher in mRNA-1273 breakthrough infections than in BNT162b2 breakthrough infections. Compared to reinfections in unvaccinated individuals, the Ct value was 1.8 (95% CI: 0.6–3.0) cycles lower in BNT162b2 breakthrough infections and 6.4 (95% CI: −0.04 to 12.9) cycles lower in mRNA-1273 breakthrough infections. The differences in Ct values were generally statistically significant with p values ≤ 0.05 except for the comparison between primary infections and BNT162b2 breakthrough infections. Notably, very few symptomatic mRNA-1273 breakthrough infections have been documented in Qatar. This may have resulted in the borderline significance in the comparison between mRNA-1273 breakthrough infections and reinfections in unvaccinated individuals (Table 5).

Table 5 RT-qPCR Ct values including only the symptomatic infections in the six pairwise comparisons.
Full size table

Discussion

Breakthrough infections in those vaccinated or who had a prior infection have higher RT-qPCR Ct values than primary infections in unvaccinated persons. While these breakthrough infections are not uncommon, the results indicate that they have lower viral load and are less likely to be infectious than primary infections. While some of these breakthrough infections could lead to secondary transmissions, and indeed some of them did have high viral loads (Fig. 3 and Supplementary Figs. 2 and 3), the risk of onward transmission is reduced, compared to primary infections. Thus, they are of less public health concern.

A consequence of these findings is that the public health benefits of vaccination may be underestimated. In addition to the conventional vaccine efficacy against acquisition of infection (VES), that is assessed in randomized clinical trials1,2,3, there is an additional “breakthrough” efficacy against transmission (VEI) that augments the benefits of VES, at least for the BNT162b2 and mRNA-1273 vaccines investigated in this study. The existence of this additional VEI efficacy implies that the reproduction number (R0) after vaccination is lower than current estimates10, and that SARS-CoV-2 incidence may decline faster with vaccine scale-up than previously thought10.

One finding of this study is that there appears to be a hierarchy in infectiousness of SARS-CoV-2 infections, where primary infections in unvaccinated persons are most infectious, followed by BNT162b2 breakthrough infections, mRNA-1273 breakthrough infections, and finally reinfections in unvaccinated persons. Strikingly, this hierarchy is the mirror image of the hierarchy observed in the efficacy against acquisition of infection. In our earlier studies on these national cohorts in Qatar, we found that those vaccinated with BNT162b2 had (relatively) the lowest protection against acquisition of infection (at 75% against the Beta variant4,6,23,42 that dominated incidence since the onset of vaccination4,5,6,8,9,23,24,26,27,28). Meanwhile, the protection against acquisition of infection was considerably higher among those vaccinated with mRNA-12735,6,42 and those with a prior infection29,30,31,43. This may indicate that both VES and VEI are essentially inter-related manifestations of the strength of the vaccine-induced (or natural-infection-induced) immune response. When the immune response is strong against acquisition of infection, it also appears strong at reducing viral replication upon acquisition of the virus, leading to lower viral load and faster infection clearance. Thus, less secondary transmission occurs in breakthrough infections.

While it is well-established that the Ct value correlates inversely with SARS-CoV-2 infectiousness17,18,19,20,21, the mathematical form of the relationship between Ct value and transmission risk, or the mathematical form of the relationship between viral load and transmission risk, are inadequately understood17,18,19,20,21,44,45,46,47. If the transmission risk has a linear decreasing relationship with Ct value, as suggested by a nationwide study from Denmark18, the observed differences in Ct value by vaccination and prior infection status may not necessarily translate into large differences in transmission risk. However, if the transmission risk has a linear or a non-linear power-law relationship with viral load, as is the case for the rigorously established transmission risk of HIV infection48,49,50, the observed differences in Ct value would translate into large differences in transmission risk.

This study has limitations. The number of documented mRNA-1273 breakthrough infections was small, thereby limiting the statistical precision of the comparisons involving these infections and leading to estimates with wider 95% CIs; perhaps making them also more prone to bias. The matching further reduced the numbers of infections used in analysis and resulted in different subsets of each category of breakthrough infections being used in the different comparisons. With the relatively small number of breakthrough infections, there was not adequate statistical power to investigate effect modifications by age or other demographic factors. With evidence for waning of vaccine immunity over time after the second dose8,51,52,53,54,55, some of the differences observed for those vaccinated with BNT162b2 versus mRNA-1273 may be explained by the shorter duration between the second dose and the breakthrough infection for those mRNA-1273 vaccinated, rather than differences in the biological immunity induced by each of these two vaccines. Future studies could investigate whether RT-qPCR Ct values are lower with increasing time since vaccination or secondary infection.

For symptomatic cases, the date of symptom onset was not available and thus an analysis of the duration between date of symptom onset and date of RT-qPCR testing, for the different categories of infections, was not possible. Only a small proportion of these documented infections was sequenced or RT-qPCR genotyped in Qatar;4,5,6,8,9,26,27,28 thus, we were unable to implement the above analyses for each variant separately. With the dominance of the Beta variant among cases included in this study (Supplementary Fig. 1), the results are most representative for this variant and may not be representative for the Alpha, Delta or other variants.

The study was implemented on documented infections, but other infections may have occurred and gone undocumented. It is possible that breakthrough infections in those vaccinated or who had a prior infection are less likely to be documented, perhaps because of minimal/mild or no symptoms. However, with the high rate of PCR testing in Qatar, the majority of infections are identified not because of testing symptomatic cases, but because of testing for other reasons, such as random testing campaigns, contact tracing, individual request, routine healthcare testing, pre-travel, and at ports of entry8. The results of our study were also consistent when we included only the randomly diagnosed infections.

Imperfect assay sensitivity and specificity of RT-qPCR testing may have affected infection ascertainment, but RT-qPCR testing was performed using a validated commercial platform that has been used globally and has essentially 100% sensitivity and specificity22 (“Methods”). Unlike blinded, randomized clinical trials, the investigated observational cohorts were neither blinded nor randomized. Our cohorts predominantly included working-age adults; therefore, results may not necessarily be generalizable to other population groups, such as children or the elderly. Matching was done for sex, age, reason for the RT-qPCR testing, and calendar week of the RT-qPCR test, but could not be done for other factors, such as comorbidities, as these were not available to study investigators. But inclusion of additional factors in the matching would have considerably reduced the sample sizes—breakthrough infections in those vaccinated or who had a prior infection are relatively uncommon. It is noteworthy that matching by age and sex may have served as a proxy for matching by co-morbidity, as comorbidities are associated with older age and may differ between females and males.

In conclusion, prior immunity, whether due to vaccination or prior infection, is associated with lower SARS-CoV-2 viral load upon infection. While breakthrough infections have been observed globally, they appear less infectious than primary infections; thus, they constitute a lesser public health concern. The public health benefits of vaccination are underestimated, as both the BNT162b2 and mRNA-1273 vaccines seem to protect not only against acquisition of infection, but also against transmission of infection. These findings justify optimism and stress the urgency to scale-up vaccination globally in order to robustly control infection transmission and the extent of the pandemic.

Methods

Data sources and study design

Analyses were conducted using the centralized, integrated, and standardized national severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) databases compiled at Hamad Medical Corporation (HMC), the main public healthcare provider and the nationally designated provider for all COVID-19 healthcare needs. Through a nationwide digital health information platform, these databases are complete, and have captured all SARS-CoV-2-related data as well as related demographic details with no missing information since the start of the epidemic, including all records of PCR testing, antibody testing, COVID-19 hospitalizations, vaccinations, infection severity classification per World Health Organization (WHO) guidelines56 (performed by trained medical personnel through individual chart reviews), and COVID-19 deaths, also assessed per WHO guidelines57.

Every real-time reverse transcription-PCR (RT-qPCR) test conducted in Qatar, regardless of location (outpatient clinic, drive-thru, or hospital, etc.), is classified on the basis of symptoms and the reason for testing (clinical symptoms, contact tracing, random testing campaigns (surveys), individual requests, routine healthcare testing, pre-travel, and port of entry). Qatar has unique demographics by sex and nationality, since expatriates from over 150 countries comprise 89% of the population34,58.

The BNT162b2 and mRNA-1273 vaccines have been the vaccines of choice in the national immunization campaign in Qatar4,5,6,8,9,23,24. Nearly all individuals in the cohorts of this study were vaccinated (free of charge) in Qatar, rather than elsewhere. In rare situations where an individual received vaccination outside Qatar, that individual’s vaccination details were still recorded in the health system at the port of entry upon return to Qatar, given the national requirements and to benefit from privileges associated with vaccination, such as quarantine exemption24.

Leveraging the national databases, effects of vaccination and of prior infection on SARS-CoV-2 infectiousness were investigated by comparing the RT-qPCR cycle threshold (Ct) values in matched cohorts of primary infections in unvaccinated individuals, reinfections in unvaccinated individuals, BNT162b2 breakthrough infections, and mRNA-1273 breakthrough infections. Only the first RT-qPCR-positive test for each individual for each infection category was used in the analysis. These types of infections are defined in the main text.

In total, six types of pairwise comparisons were conducted on the cohorts of this study, after exact matching to the first match in a 1:1 ratio by sex, 10-year age group, reason for RT-qPCR testing, and calendar week of the RT-qPCR test, to control for differences in biology by sex and age, as well as exposure risk34,35,36,37,38 and variant exposure4,5,6,8,9,26,27,28. It is noteworthy that the first SARS-CoV-2 epidemic wave in Qatar occurred before introduction of any variant of concern and peaked in late May, 202034,35. The second wave was triggered by introduction and expansion of the Alpha25 (B.1.1.7) variant and peaked in early March, 20214,5,6,8,9,26,27,28. The third wave was dominated by the Beta25 (B.1.351) variant, and peaked in the first week of April, 20214,5,6,8,9,26,27,28. The Delta25 (B.1.617.2) variant has been introduced only recently in Qatar, and as of July 11, 2021, it remains at low incidence4,5,6,8,9,26,27,28. There is no evidence that any other variant of concern is or has been responsible for appreciable community transmission in Qatar4,5,6,8,9,26,27,28.

Comparisons across the cohorts of infection were implemented for all RT-qPCR-confirmed infections, for only the symptomatic infections defined as RT-qPCR-positive tests conducted because of clinical suspicion due to symptoms compatible with a respiratory tract infection, and for only the randomly diagnosed (asymptomatic) infections defined as RT-qPCR-positive tests conducted with no prior reason to suspect infection and no reported symptoms compatible with a respiratory tract infection. The latter was strictly defined as an RT-qPCR-positive test conducted as part of a survey (random testing campaigns), for routine healthcare testing, as a pre-travel requirement, or at a port of entry upon arrival in the country4,5,6,34.

All records of RT-qPCR testing in Qatar were examined in this study, but only samples of matched cohorts were included in the analysis. Individuals with a record of a SARS-CoV-2 antibody-positive test before the first RT-qPCR-positive test were excluded from analysis of those with primary infections. Individuals with a record of vaccination before the reinfection diagnosis were excluded from the analysis of those with reinfection. Only breakthrough infections in fully vaccinated individuals were included in the analysis. Being fully vaccinated was defined as completion of ≥14 days after the second dose. Only individuals who had their first and second doses with the same vaccine were included in the analysis.

Reporting of the study followed the STROBE guidelines (Supplementary Table 3).

Laboratory methods

Nasopharyngeal and/or oropharyngeal swabs were collected for PCR testing and placed in Universal Transport Medium (UTM). Aliquots of UTM were extracted on a QIAsymphony platform (QIAGEN, USA) and tested with RT-qPCR using TaqPath COVID-19 Combo Kits (100% sensitivity and specificity;22 Thermo Fisher Scientific, USA) on an ABI 7500 FAST (Thermo Fisher, USA); extracted using a custom protocol59 on a Hamilton Microlab STAR (Hamilton, USA) and tested using AccuPower SARS-CoV-2 Real-Time RT-PCR Kits (100% sensitivity and specificity;60 Bioneer, Korea) on an ABI 7500 FAST; or loaded directly into a Roche cobas 6800 system and assayed with a cobas SARS-CoV-2 Test (95% sensitivity, 100% specificity61; Roche, Switzerland). The first assay targets the viral S, N, and ORF1ab regions. The second targets the viral RdRp and E-gene regions, and the third targets the ORF1ab and E-gene regions.

For standardization of Ct values, only the RT-qPCR-confirmed infections diagnosed using the TaqPath COVID-19 Combo Kits platform (Thermo Fisher Scientific, USA22) were included in the analysis. This platform was used for >85% of all RT-qPCR tests in Qatar and reports individual Ct values for each of the N, ORF1ab, and S genes4,5,6,22. The average of these three Ct values was included in the analysis. The correlation between each pair of these Ct values across all RT-qPCR-positive tests was very strong with a Pearson correlation coefficient ≥0.976. In the case of a gene “target failure”, and specifically an S-gene “target failure” that was the defining characteristic of the Alpha cases62,63,64, the average was determined using only the two Ct values of the N and ORF1ab genes.

Antibodies against SARS-CoV-2 in serological samples were detected using a Roche Elecsys Anti-SARS-CoV-2 assay (99.5% sensitivity65, 99.8% specificity;65,66 Roche, Switzerland), an electrochemiluminescence immunoassay that uses a recombinant protein representing the nucleocapsid (N) antigen for antibody binding. Results were interpreted according to the manufacturer’s instructions (reactive: optical density (proxy for antibody titer67) cutoff index ≥1.0 vs. non-reactive: optical density cutoff index <1.0).

All testing was conducted at the HMC Central Laboratory or Sidra Medicine Laboratory, following standardized protocols.

Statistical analysis

Socio-demographic characteristics of study samples were described using frequency distributions and measures of central tendency. Differences in proportions across categorical variables between study groups were evaluated using Chi-square tests. The distributions of the RT-qPCR Ct values were illustrated using scatter plots and boxplots, and summarized using measures of central tendency and dispersion. Mean differences in the RT-qPCR Ct values between study groups and associated 95% CIs were calculated using paired t-tests. Two-sided p value of ≤0.05 indicated a significant association.

To provide an additional validation of the results, associations with a higher RT-qPCR Ct value across study groups and matching factors were explored using analysis of variance tests and univariable linear regressions. Beta coefficients (βs), 95% CIs, and two-sided p-values were reported. All covariates were included in the multivariable analysis to estimate the adjusted βs and associated 95% CIs and p values. Covariates with p value ≤ 0.05 in the multivariable model were considered associated with the RT-qPCR Ct value.

Statistical analyses were conducted in STATA/SE version 17.068.

Ethical approval

The study was approved by the Hamad Medical Corporation (MRC-05-011) and Weill Cornell Medicine—Qatar (20-00017) Institutional Review Boards with waiver of informed consent.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.