Increased household secondary attacks rates with Variant of Concern SARS-CoV-2 index cases ========================================================================================== * Sarah A. Buchan * Semra Tibebu * Nick Daneman * Michael Whelan * Thuva Vanniyasingam * Michelle Murti * Kevin A. Brown ## Abstract **IMPORTANCE** Higher secondary attack rates related to variant of concern (VOC) index cases have been reported, but have not been explored within households, which continue to be an important source of coronavirus disease 2019 (COVID-19) transmission **OBJECTIVE** To compare secondary attack rates in households with VOC versus non-VOC index cases. **DESIGN** A retrospective cohort study of household index cases reported from February 7 – 27, 2021. A propensity-score matched cohort was derived to calculate adjusted estimates. **SETTING** Ontario, Canada. **PARTICIPANTS** A population-based cohort of all private households with index cases. We excluded cases in congregate settings, as well as households with one individual or with >1 case with the same earliest symptom onset date. **EXPOSURE** VOC status, defined as either individuals confirmed as B.1.1.7 using whole genome sequencing or those that screened positive for the N501Y mutation using real-time PCR. **MAIN OUTCOME AND MEASURE** Household secondary attack rate, defined as the number of household secondary cases that occurred 1-14 days after the index case divided by the total number of household secondary contacts. **RESULTS** We included 1,259 index VOC and non-VOC cases in the propensity score-matched analysis. The secondary attack rate for VOC index cases in this matched cohort was 1.31 times higher than non-VOC index cases (RR=1.31, 95%CI 1.14-1.49), similar to the unadjusted estimate. In stratified analyses, the higher secondary attack rate for VOC compared to non-VOC index cases was accentuated for asymptomatic index cases (RR=1.91, 95% CI 0.96-3.80) and presymptomatic cases (RR=3.41, 95%CI 1.13-10.26) **CONCLUSIONS AND RELEVANCE** This study provides strong evidence of increased transmissibility in households due to VOCs and suggests that asymptomatic and pre-symptomatic transmission may be of particular importance for VOCs. Our study suggests that more aggressive public health measures will be needed to control VOCs and that ongoing research is needed to understand mechanisms of VOC transmissibility to curb their associated morbidity and mortality. ## Introduction The prevalence of variants having the N501Y mutation have rapidly increased globally,1 including in Ontario, Canada, where this prevalence increased dramatically in February, 2021. These patterns of rapid strain replacement suggest increased transmissibility of variants with the N501Y mutation, which is present across variants of concern (VOC), including B.1.1.7, B.1.351, and P.1 lineages.2 However the exact degree of increased transmissibility, and specific settings when increased transmissibility occurs, remains unclear. 1, 3,4,5 Higher secondary attack rates related to VOC index cases have been reported, 6 but have not been explored within households, which continue to be an important source of coronavirus disease 2019 (COVID-19) transmission.7 The objective of this study was to compare secondary attack rates in households with VOC versus non-VOC index cases in Ontario, Canada. ## Methods We identified individuals with laboratory-confirmed COVID-19 reported in the Public Health Case and Contact Management Solution (CCM), Ontario’s COVID-19 reporting system, and included households with index cases reported from February 7 to 27, 2021. VOC cases included either individuals confirmed as B.1.1.7 using whole genome sequencing or those that screened positive for the N501Y mutation using real-time PCR. B.1.1.7 accounted for 94% of confirmed VOC cases reported in Ontario until February 27, 2021.8 All PCR-positive specimens in Ontario with cycle threshold ≤35 underwent screening for the N501Y mutation during the study period using real-time PCR. Non-VOC cases were those that screened negative for the N501Y mutation. We grouped cases living in the same household based on residential address.9 Index cases were defined as the first case in the household based on symptom onset date (or specimen collection date, if symptom onset date was not available) and secondary cases were included if they occurred 1-14 days after the index case. We excluded cases in congregate settings, as well as households with one individual or with >1 case with the same earliest symptom onset date. We used reported household size to calculate secondary attack rates by dividing the number of secondary cases by the total number of household secondary contacts (i.e., household size minus one). Cases with and without symptoms were classified based on symptom information and onset date reported in CCM. Poisson regression was performed for the unadjusted and adjusted analyses, overall and by strata. The models were specified with the count of household secondary cases as the outcome, the logarithm of household size as the offset and VOC status as a binary exposure covariate. Clustering was accounted for with a random intercept for household to account for known overdispersion. Our unadjusted risk ratios (RR) included the entire cohort of households with a VOC or non-VOC index case in Ontario. Adjusted RRs were based on a propensity score matched analysis. VOC index cases were 1:1 matched to non-VOC index cases based on gender and age group, and matched on the logit of the propensity score, using a caliper width of 0.2 times the standard deviation.10 The propensity score was based on a logistic regression model of VOC status as a function of five covariates: reported date, time between symptom onset and testing, association with a reported outbreak, as well as the neighbourhood proportion of visible minority residents (non-white and non-Indigenous population) and household crowding as determined using 2016 Canadian Census data.11 Regression estimates were reported using RRs and 95% confidence intervals (CI). Statistical analysis was performed in R version 4.0.4.12 This study was approved by Public Health Ontario’s Research Ethics Board. ## Results We identified 5,617 index cases and 3,397 secondary cases across the study period. Amongst index cases, 1,318 were classified as VOC (151 B.1.1.7 and 1,167 N501Y) and 4,299 were classified as non-VOC. The overall secondary attack rate was higher for VOC index cases (25.9%) compared to non-VOC (20.5%, p<0.01) with consistently higher secondary attack rates for VOCs across individual characteristics of the index cases (Table 1). The secondary attack rate of VOC index cases was 1.28 times higher than that of non-VOC index cases (RR=1.28, 95% CI 1.16-1.42); a similar trend of increased secondary attack rate was observed across subgroups (Figure 1). View this table: [Table 1.](http://medrxiv.org/content/early/2021/04/05/2021.03.31.21254502/T1) Table 1. Characteristics of index cases by VOC status and household secondary attack rate ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/05/2021.03.31.21254502/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/04/05/2021.03.31.21254502/F1) Figure 1. Risk ratio comparing household secondary attack rate associated with VOC vs. non-VOC index cases by characteristic of index case, unadjusted (black) and adjusted by 1:1 propensity-score matching (green) estimates Note: Unadjusted risk ratio accounts for household clustering and includes the full cohort. Adjusted risk ratio refers to the propensity score-matched cohort. We included 1,259 index VOC and non-VOC cases in the propensity score-matched analysis. The secondary attack rate for VOC index cases in this matched cohort was 1.31 times higher than non-VOC index cases (RR=1.31, 95%CI 1.14-1.49), similar to the unadjusted estimate. In stratified analyses, the higher secondary attack rate for VOC compared to non-VOC index cases was accentuated for asymptomatic index cases (RR=1.91, 95% CI 0.96-3.80) and presymptomatic cases (RR=3.41, 95%CI 1.13-10.26) ## Discussion In our cohort, we estimated that the household secondary attack rate of VOC index cases was 31% higher than non-VOC index cases, providing evidence of increased transmissibility. This is consistent with previous VOC (i.e., B.1.1.7) secondary attack rate estimates from the United Kingdom (relative secondary attack rate (32%, 12.9% vs 9.7% among all contacts);6 however, households are an important contributor to COVID-19 and provide a valuable setting in which to examine transmission.7 Our estimates of increased transmission from asymptomatic and pre-symptomatic VOC index cases has not previously been reported and suggests an increased importance of prevention measures when individuals are not aware of their infection. The biological mechanism responsible for increased transmissibility has not been identified, but hypotheses include higher viral loads (i.e. increased transmission potential per contact), and prolonged viral shedding (i.e., longer infectious period).3 The N501Y mutation has also been associated with enhanced binding affinity of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) to angiotensin-converting enzyme 2 (ACE2) receptors.13,14 Limitations of this study include potential misclassification of secondary cases as index cases and small sample sizes in some subgroups. We may have underestimated secondary attack rates as we only captured diagnosed secondary cases and we lacked testing data on all household contacts; however, we do not believe this would be differential by VOC status of the index case. Ontario implemented more stringent measures for close contacts of all cases (not just VOC cases) in early February in response to VOC introductions, including increased frequency of testing during quarantine and outbreaks.15 If household contacts were more likely to test once VOC status was known, a similar increase in the risk ratio would be expected across symptom status, suggesting an independent effect of asymptomatic and presymptomatic VOC cases on transmission. This study provides strong evidence of increased transmissibility in households due to VOCs and suggests that asymptomatic and pre-symptomatic transmission may be of particular importance for VOCs. Our study suggests that more aggressive public health measures will be needed to control VOCs. While measures effective for persons with unknown disease status such as physical distancing and masking may continue to be highly effective, measures focused on symptomatic individuals, such as public health contact tracing, may be increasingly ineffective unless extremely rapid. Ongoing research is needed to understand mechanisms of VOC transmissibility to curb their associated morbidity and mortality. ## Data Availability Data sharing requests should be directed to Public Health Ontario. ## Footnotes * sarah.buchan{at}oahpp.ca * semra.tibebu{at}oahpp.ca * nick.daneman{at}sunnybrook.ca * michael.whelan{at}oahpp.ca * thuva.vanniyasignam{at}oahpp.ca * michelle.murti{at}oahpp.ca * kevin.brown{at}oahpp.ca * Received March 31, 2021. * Revision received March 31, 2021. * Accepted April 5, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.European Centre for Disease Prevention and Control. Risk related to the spread of new SARS-CoV-2 variants of concern in the EU/EEA-first update. 2021. [https://www.ecdc.europa.eu/sites/default/files/documents/COVID-19-risk-related-to-spread-of-new-SARS-CoV-2-variants-EU-EEA-first-update.pdf](https://www.ecdc.europa.eu/sites/default/files/documents/COVID-19-risk-related-to-spread-of-new-SARS-CoV-2-variants-EU-EEA-first-update.pdf). Accessed March 10, 2021. 2. 2.Martin DP, Weaver S, Tegally H, et al. The emergence and ongoing convergent evolution of the N501Y lineages coincides with a major global shift in the SARS-CoV-2 selective landscape [published online March 10, 2021]. medRxiv. doi:10.1101/2021.02.23.21252268 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wMi4yMy4yMTI1MjI2OHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMDUvMjAyMS4wMy4zMS4yMTI1NDUwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. 3.Davies NG, Abbott S, Barnard SC, et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England [published online March 3, 2021]. Science. doi:10.1126/science.abg3055 [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiJzY2llbmNlLmFiZzMwNTV2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA0LzA1LzIwMjEuMDMuMzEuMjEyNTQ1MDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. 4.Gaymard A, Bosetti P, Feri A, et al. Early assessment of diffusion and possible expansion of SARS-CoV-2 Lineage 20I/501Y.V1 (B.1.1.7, variant of concern 202012/01) in France, January to March 2021. Euro Surveill. 2021;26(9):pii=2100133. doi:10.2807/1560-7917.ES.2021.26.9.2100133 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2021.26.9.2100133&link_type=DOI) 5. 5.Volz E, Mishra S, Chand M, et al. Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data. 2020. [https://virological.org/t/transmission-of-sars-cov-2-lineage-b-1-1-7-in-england-insights-from-linking-epidemiological-and-genetic-data/576](https://virological.org/t/transmission-of-sars-cov-2-lineage-b-1-1-7-in-england-insights-from-linking-epidemiological-and-genetic-data/576). Accessed March 10, 2021 6. 6.Public Health England. Investigation of novel SARS-CoV-2 variant. Variant of Concern 202012/01: Technical briefing 5. 2021. [https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment\_data/file/](https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/) 959426/Variant\_of\_Concern\_VOC\_202012\_01\_Technical\_Briefing\_5.pdf. Accessed March 10, 2021. 7. 7.Madewell ZJ, Yang Y, Longini IM, Halloran ME, Dean NE. Household Transmission of SARS-CoV-2: A Systematic Review and Meta-analysis. JAMA Netw Open. 2020;3(12):e2031756. doi:10.1001/jamanetworkopen.2020.31756 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2020.31756&link_type=DOI) 8. 8.Ontario Agency for Health Protection and Promotion (Public Health Ontario). SARS-CoV-2 (COVID-19 virus) variant of concern (VoC) surveillance. Toronto, ON: Queen’s Printer for Ontario. 2021. [https://www.publichealthontario.ca/en/laboratory-services/test-information-index/covid-19-voc](https://www.publichealthontario.ca/en/laboratory-services/test-information-index/covid-19-voc). Accessed March 28, 2021. 9. 9.Paul LA, Daneman N, Brown KA, et al. Characteristics associated with household transmission of SARS-CoV-2 in Ontario, Canada: A cohort study [published online March 5, 2021]. Clin Infect Dis. doi:10.1093/cid/ciab186 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciab186&link_type=DOI) 10. 10.Austin PC. Optimal caliper widths for propensity-score matching when estimating differences in means and differences in proportions in observational studies. Pharm Stat. 2011;10(2):150–161. doi:10.1002/pst.433 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pst.433&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20925139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F05%2F2021.03.31.21254502.atom) 11. 11.Statistics Canada. Dictionary, Census of Population, 2016. 2017. [https://www12.statcan.gc.ca/census-recensement/2016/ref/dict/az1-eng.cfm?topic=az1](https://www12.statcan.gc.ca/census-recensement/2016/ref/dict/az1-eng.cfm?topic=az1) Accessed March 10, 2021. 12. 12.R: A language and environment for statistical computing. Version 4.0.4. 2021. R Foundation for Statistical Computing, Vienna, Austria. [http://www.R-project.org/](http://www.R-project.org/) 13. 13.Gu H, Chen Q, Yang G, et al. Adaptation of SARS-CoV-2 in BALB/c mice for testing vaccine efficacy. Science. 2020;369(6511):1603–1607. doi:10.1126/science.abc4730 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUxMS8xNjAzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMDUvMjAyMS4wMy4zMS4yMTI1NDUwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. 14.Barnes CO, Jette CA, Abernathy ME, et al. SARS-CoV-2 neutralizing antibody structures inform therapeutic strategies. Nature. 2020;588, 682–687. doi:10.1038/s41586-020-2852-1 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2852-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33045718&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F05%2F2021.03.31.21254502.atom) 15. 15.Ministry of Health. COVID-19 Variant of Concern: Case, Contact and Outbreak Management Interim Guidance. 2021. [https://www.health.gov.on.ca/en/pro/programs/publichealth/coronavirus/docs/VOC\_guidance.pdf](https://www.health.gov.on.ca/en/pro/programs/publichealth/coronavirus/docs/VOC_guidance.pdf) Accessed March 25, 2021.