Environmental drivers of SARS-CoV-2 lineage B.1.1.7 transmission in England, October to December 2020 ===================================================================================================== * Thomas P. Smith * Ilaria Dorigatti * Swapnil Mishra * Erik Volz * Patrick G. T. Walker * Manon Ragonnet-Cronin * Michael Tristem * William D. Pearse ## Abstract Previous work has shown that environment affects SARS-CoV-2 transmission, but it is unclear whether emerging strains show similar responses. Here we show that lineage B.1.1.7 spread with greater transmission in colder and more densely populated parts of England. We also find evidence of B.1.1.7’s transmission advantage at warmer temperatures versus other strains, implying that spring conditions may facilitate B.1.1.7’s invasion in Europe and across the Northern hemisphere, undermining the effectiveness of public health interventions. Keywords * SARS-CoV-2 * lineage B.1.1.7 * transmission * environment * population density * climate ## Environmental drivers of SARS-CoV-2 transmission A year into the COVID-19 pandemic, quantifying the factors driving SARS-CoV-2 transmission is still key for the optimal implementation of control strategies. SARS-CoV-2 transmission is influenced by local environmental conditions [1], and we have previously shown that colder temperatures and higher population densities explain spatial and temporal variation in transmission intensity [2]. A critical insight from previous work has been that human behaviour drives SARS-CoV-2 transmission, and that environment only plays a marginal effect when effective non-pharmaceutical interventions (NPIs) are in place [1,2]. While several studies have investigated the impact of environmental factors on SARS-CoV-2 transmission during the early stages of the pandemic, numerous new variants have recently emerged (in for example, the UK [3], Brazil [4], and South Africa [5]) with greater transmissibility than strains present during the early stages of the pandemic [6,7]. Any effect of environment must therefore be re-assessed in the context of a landscape of emerging viral lineages. ### Impacts of environment on SARS-CoV-2 transmission in England The UK SARS-CoV-2 ‘Variant of Concern’ (VOC; also known as lineage B.1.1.7 [3]) rapidly spread through the UK population in late 2020, showing greater transmission intensity (as measured by the reproduction number, *R*) than other circulating strains [7]. The wealth of data collected for this strain makes the VOC a strong test-case for investigating environmental effects in the context of emerging lineages. To assess the impact of environment on VOC transmission, we extracted weekly *R* estimates calculated for fine-grained regions of England (Lower Tier Local Authorities; LTLAs) [8], and regressed these against VOC frequency, temperature, and population density for each LTLA. Population density and daily temperature data were averaged across each LTLA, using the same methods as Smith et al. [2]. Daily temperatures were combined into weekly mean temperatures, corresponding to each of the weekly *R* estimates. The weekly VOC frequencies were calculated for each LTLA as the percentage of true positive rate (TPR)-adjusted S-gene target negative test to positive test results, as in Volz et al. [7]. We assessed the impact of the environment by investigating the association between weekly *R* estimates and the predictors (*i*.*e*., spatial variation in *R*), for sequential weeks. To account for the VOC’s emergence in the South-East (which is among the warmest and most densely populated parts of England), and its subsequent spread into colder, less densely populated parts of the country, we incorporated interactions between VOC frequency and each of temperature and population density. Without including this interaction, any effect of temperature or population density on *R* might be compounded or masked by spatial autocorrelation with VOC frequency. Following Volz et al. [7], we measured the effects of the weekly predictors on *R* in the subsequent week, to account for the generation time of SARS-CoV-2. As observed by Volz et al. [7], we find a strong positive effect of VOC frequency on transmission intensity, particularly during the initial emerging phase when VOC frequency was less consistent among regions (Table 1). While the UK’s national lockdown initially decreased SARS-CoV-2 transmission intensity, as the VOC spread, *R* increased across regions (Figure 1). We find that the impact of the environment on transmission is mediated by VOC frequency, and it is only at higher VOC frequencies that the effects of temperature and population density are pronounced (Table 2 and Figure 1). Furthermore, we find a greater environmental effect after the UK moved from full lockdown to a (less strict) tiered system, confirming our previous observations that the effect of temperature becomes pronounced when NPIs are relaxed [2]. Thus, the effects of environmental drivers, and even those of the increased transmissibility of the VOC, are secondary to differences in human behaviour driven by differences in NPIs. As the VOC has now spread to dominate the country’s viral population, we now expect environmental factors to become the dominant driver of spatial variation in *R*, with colder areas facing higher transmission intensities than warmer regions, unless NPIs or the accumulation of immunity (either naturally acquired or vaccine derived) sufficiently reduce transmission. View this table: [Table 1.](http://medrxiv.org/content/early/2021/03/12/2021.03.09.21253242/T1) Table 1. VOC frequency is the strongest driver of transmission during early stages of the VOC sweep (week 46, starting 9 November 2020). Linear regression model with VOC frequency, temperature, and population density as predictors of Rt, r2 = 48%, F5,293 = 53.5,p < 0.001. Predictors were scaled to have mean = 0 and SD = 1 and thus coefficients are measures of the relative importance of each variable. Frequency of the variant is a comparably strong driver of R. * = p < 0.05. View this table: [Table 2.](http://medrxiv.org/content/early/2021/03/12/2021.03.09.21253242/T2) Table 2. Both temperature and population density significantly interact with VOC frequency during late stages of the VOC sweep (week 50, starting 7 December 2020). r2 = 55%, F5,298 = 71.6, p < 0.001. Predictors were scaled to have mean = 0 and SD = 1. Significant interaction effects between variant frequency and environmental parameters are observed. Significant interactions show that at high VOC frequencies, there is a negative effect of temperature and a positive effect of population density on Rt. * = p < 0.05. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/12/2021.03.09.21253242/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/03/12/2021.03.09.21253242/F1) Figure 1: Transmission dynamics during emergence and spread of SARS-CoV-2 lineage B.1.1.7. (**A**.) Median (bold) and interquartile range (box) of the transmission intensity estimates (R) for weeks 43-50 of 2020, with point estimates for each LTLA overlaid. Points are coloured by the frequency of the VOC in positive COVID-19 tests (%), boxplots are coloured by the median VOC frequency across all LTLAs. The first date of each week is given within parentheses. As the VOC emerged, the UK entered a national lockdown, leading to an initial decrease in R across regions. However, as the VOC spread through space to become the dominant lineage, R increased across regions despite the national lockdown. When the UK re-entered a tiered NPI system on the 2nd of December 2020, average R further increased across regions. Panels **B** and **C** show the observed (points) and expected (background) effects of temperature and VOC frequency on R during early (week 46) and late (week 50) phases of the variant sweep, respectively (legend shared across both plots). ### The role of natural immunity To investigate the potential role of immunity in our results, we perform a sensitivity analysis and re-run our regression analyses on *R* estimates corrected for the relative attack rate (AR) estimates. This accounts for the potential slowdown in *R* due to the accumulation of natural immunity across the population through time. We obtain qualitatively identical results in the AR-corrected analyses, but generally lower correlation (additive decrease in *r*2 of 16% in weeks 45 and 46) as the VOC emerged, versus in the later stages (average decrease in *r*2 of only 3%, Supplementary Table S1). Note that the timeframe of our analyses was prior to mass vaccination of the UK public and even in the latest timepoint (week 50), median attack rate across regions was less than 10% (Supplementary Figure S1). ### Environmental response of lineage B.1.1.7 transmission We additionally investigated whether the VOC may respond differently to the environment as compared to non-VOC strains. Volz et al. [7] estimated the transmission intensity of VOC vs non-VOC strains and here we regressed the ratios of VOC *R* to non-VOC *R* against temperature, to understand under which environmental conditions the VOC may have an advantage. This was performed at a larger spatial scale than our LTLA analyses (English NHS sustainability and transformation partnerships – STP regions). We found a small, negative correlation between temperature and the *R* ratio (Supplementary Table S2), *i*.*e*., at warmer temperatures the VOC was *even more* transmissible than non-VOC strains, than in colder conditions. This may be expected, as differences in fitness between strains are likely to be magnified in harsher environments where fitness is reduced overall [9]. This effect may have played into the dynamics of the spread of the VOC in the South East of England by enhancing its spread even more in regions where previous strains’ transmission intensities were lower. ## Discussion Whereas, previously, a central question has been whether winter conditions would enhance SARS-CoV-2 transmission, in the context of these new strains it is important to ask whether warmer summer conditions could slow transmission. Here we have shown that temperature can drive variation in transmission intensity between regions and between variants, but this effect is likely secondary to the overall increased transmission intensity of the new variant. Our results show that variant B.1.1.7 is more transmissible in colder areas. Using the model fitted to the week 50 data, our results suggest that when the variant accounts for 100% of SARS-CoV-2 cases in the population, *R* increases by approximately 0.12 per °C of temperature decrease. However, we also find that in warmer conditions, the difference in transmissibility between VOC and non-VOC strains is greater. Our results suggest that for each °C of temperature increase, the ratio of VOC to non-VOC *R* may increase by approximately 0.12-0.22 (Supplementary Table S2). While we caution against extrapolating from winter to summer conditions, our results highlight that there is no reason to suppose that summer weather alone will slow down the invasion dynamics of B.1.1.7 and significantly reduce the transmission intensity of SARS-CoV-2. Thus, it is imperative to quantify and continue to monitor the impact of interventions (*e*.*g*., vaccines and NPIs) to inform policy. In that regard, the speed with which the VOC has spread through the UK, while concerning from a public health perspective, provides the perfect opportunity to parameterise models of its responses. Importantly in the context of new strains, transmission of respiratory viruses in tropical climates (such as Brazil, where a new variant has arisen [4]) are often more sensitive to rainfall than temperature [10]. Better characterisation of past, and new potential relationships between environment and transmission, as new variants emerge and spread internationally, is a research priority. ## Supporting information Supplementary Material [[supplements/253242_file02.pdf]](pending:yes) ## Data Availability No new data was collected for this study. The data used is archived in a github repository alongside the analysis code. [https://github.com/smithtp/covid19-variant-N501Y](https://github.com/smithtp/covid19-variant-N501Y) ## Footnotes * **Conflict of interests:** The authors declare no conflicts of interest. * **Funding Statement:** This work was funded by a Natural Environment Research Council grant NE/V009710/1 and by the Imperial College COVID-19 Research Fund. * We emphasise that this pre-print manuscript has not been peer-reviewed, and thus its findings should be considered, at best, preliminary. * Received March 9, 2021. * Revision received March 9, 2021. * Accepted March 12, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1.Poirier C, Luo W, Majumder MS, Liu D, Mandl KD, Mooring TA, et al. The Role of Environmental Factors on Transmission Rates of the COVID-19 Outbreak: An Initial Assessment in Two Spatial Scales. Sci Rep 2020;10:17002. 2. 2.Smith TP, Flaxman S, Gallinat AS, Kinosian SP, Stemkovski M, Watson OJ, et al. Environment influences SARS-CoV-2 transmission in the absence of non-pharmaceutical interventions. medRxiv [Preprint]. 2020 [cited 2021 Mar 9]. Available from: [https://www.medrxiv.org/content/10.1101/2020.09.12.20193250v1](https://www.medrxiv.org/content/10.1101/2020.09.12.20193250v1) doi: 10.1101/2020.09.12.20193250 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wOS4xMi4yMDE5MzI1MHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMy4wOS4yMTI1MzI0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. 3.Rambaut A, Loman N, Pybus O, Barclay W, Barrett J, Carabelli A, et al. Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations. Virological. 2020 [cited 2021 Mar 9]. Available from: [https://virological.org/t/563](https://virological.org/t/563) 4. 4.Voloch CM, da Silva RF Jr., de Almeida LGP, Cardoso CC, Brustolini OK, Gerber AL, et al. Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil. medRxiv [preprint]. 2020 [cited 2021 Mar 9]. Available from: [https://www.medrxiv.org/content/10.1101/2020.12.23.20248598v1](https://www.medrxiv.org/content/10.1101/2020.12.23.20248598v1) doi: 10.1101/2020.12.23.20248598 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMi4yMy4yMDI0ODU5OHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMy4wOS4yMTI1MzI0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. 5.Tegally H, Wilkinson E, Giovanetti M, Iranzadeh A, Fonseca V, Giandhari J, et al. Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa. medRxiv [preprint]. 2020 [cited 2021 Mar 9]. Available from: [https://www.medrxiv.org/content/10.1101/2020.12.21.20248640v1](https://www.medrxiv.org/content/10.1101/2020.12.21.20248640v1) doi: 10.1101/2020.12.21.20248640 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMi4yMS4yMDI0ODY0MHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMy4wOS4yMTI1MzI0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 6. 6.Faria NR, Mellan TA, Whittaker C, Claro IM, Candido DS, Mishra S, et al. Genomics and epidemiology of a novel SARS-CoV-2 lineage in Manaus, Brazil. medRxiv [preprint]. 2021 [cited 2021 Mar 9]. Available from: [https://www.medrxiv.org/content/10.1101/2021.02.26.21252554v1](https://www.medrxiv.org/content/10.1101/2021.02.26.21252554v1) doi: 10.1101/2021.02.26.21252554 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wMi4yNi4yMTI1MjU1NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMy4wOS4yMTI1MzI0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 7. 7.Volz E, Mishra S, Meera Chand, Barrett JC, Johnson R, Geidelberg L, et al. Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data. medRxiv [preprint]. 2021 [cited 2021 Mar 9]. Available from: [https://www.medrxiv.org/content/10.1101/2020.12.30.20249034v2](https://www.medrxiv.org/content/10.1101/2020.12.30.20249034v2) doi: 10.1101/2020.12.30.20249034 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMi4zMC4yMDI0OTAzNHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMy4wOS4yMTI1MzI0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 8. 8.Mishra S, Scott J, Zhu H, Ferguson NM, Bhatt S, Flaxman S et al. A COVID-19 Model for Local Authorities of the United Kingdom. medRxiv [preprint]. 2021 [cited 2021 Mar 9]. Available from: [https://www.medrxiv.org/content/10.1101/2020.11.24.20236661v1](https://www.medrxiv.org/content/10.1101/2020.11.24.20236661v1) doi: 10.1101/2020.11.24.20236661 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMS4yNC4yMDIzNjY2MXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMy4wOS4yMTI1MzI0Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. 9.Soberon J, Peterson AT. Interpretation of Models of Fundamental Ecological Niches and Species’ Distributional Areas. Biodiversity Informatics 2005;2:1–10. 10. 10.Pica N, Bouvier NM. Environmental factors affecting the transmission of respiratory viruses Curr Opin Virol. 2012;2(1):90–95.