The SARS-CoV-2 B.1.351/V2 variant is outgrowing the B.1.1.7/V1 variant in French regions in April 2021 ====================================================================================================== * Bénédicte Roquebert * Sabine Trombert-Paolantoni * Stéphanie Haim-Boukobza * Emmanuel Lecorche * Laura Verdurme * Vincent Foulongne * Mircea T. Sofonea * Samuel Alizon ## Abstract SARS-CoV-2 variants threaten our ability to control COVID-19 epidemics. We analyzed 36,590 variant-specific RT-PCR tests performed on samples collected between April 12 and May 7, 2021 in France to compare variant spread. Contrarily to January to March 2021, we found that the V2 variant had a significant transmission advantage over V1 in some regions (15.1 to 16.1% in Île-de-France and 16.1 to 18.8% in Hauts-de-France). This shift in transmission advantage is consistent with the immune evasion abilities of V2 and the high levels of immunization in these regions. Keywords * COVID-19 * virus * variant * RT-PCR * epidemiology * statistical modelling * immune evasion ‘Variants of concern’ (VOC) are SARS-CoV-2 phenotypically distinct lineages that are associated with major epidemiological and/or clinical shifts. Variant V1 from lineage B.1.1.7 is currently causing the majority of infections in Europe and North America [1], whereas variant V2 from lineage B.1.351 is dominant in South-Africa [2], and variant V3 from lineage P.1 dominates in Brazil and South America [3]. The outcome of the (indirect) competition between variants is an open question. In France, the early introduction of the V2 variant in some regions makes it particularly important to monitor the spread of different variants [4]. Since January 2021, the national guideline is to test all positive samples with an additional RT-PCR to detect implicated variants [5,6]. Since April 2021, the guideline for variant-specific RT-PCRs is to target the N501Y and E484K mutations. ## Data We use the ID™ SARS-CoV-2/N501Y/E484K Quadruplex assay from ID Solution to test N=53,687 SARS-CoV-2 positive samples collected between Monday April 12 and Friday May 7 2021 in 13 French regions, mainly from the Ile-de-France region (47.9%). Some of the samples (7.95%) originate from hospitals (mostly hospitalized patients). We only analyzed data from individuals aged from 5 to 80 years old to minimize the weight of preschool children and aged care facilities in the analysis. 17.3% of the tests could not be interpreted, because of Ct value too high to keep an equal sensitivity for the N501Y and E484K targets, and were removed from the analysis. Therefore, to avoid biasing the variant screening, all tests with cycle thresholds (Ct) values above 30 were ignored (28%). Overall, we analyze 36,590 tests, i.e. 68.2% of the all the tests performed (Table 1). Raw proportions are shown in Figure 1. View this table: [Table 1.](http://medrxiv.org/content/early/2021/05/17/2021.05.12.21257130/T1) Table 1. Characteristics of the N=36,590 samples analyzed via variant-specific RT-PCR. The interquartile range (IQR) shown is the 95% confidence interval. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/17/2021.05.12.21257130/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/05/17/2021.05.12.21257130/F1) Figure 1: Raw daily frequencies of variant-specific test results in French regions. The number of tests performed is indicated in each panel. Regions with few tests can exhibit strong variations in frequencies. Wild type corresponds to N501Y-/E484K-, B.1.351/V2 to N501Y+/E484K+, B.1.525 to N501Y-/E484K+, and B.1.1.7 to N501Y+/E484-. See supplementary figure S1 for the raw numbers instead of the proportion. ![Figure S1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/17/2021.05.12.21257130/F3.medium.gif) [Figure S1:](http://medrxiv.org/content/early/2021/05/17/2021.05.12.21257130/F3) Figure S1: Raw daily occurrences of variant-specific test results in French regions. The number of tests performed is indicated in each panel. Wild type corresponds to N501Y-/E484K-, B.1.351/V2 to N501Y+/E484K+, B.1.525 to N501Y-/E484K+, and B.1.1.7 to N501Y+/E484-. ### Sequencing profiles The specificity of the variant-specific RT-PCR we used is limited since they only target 2 mutations. To gain additional insights, we sequenced all the samples collected on the same day (March 30) for which the Ct was lower than 28 using Twist Libraries and Illumina sequencing (currently being uploaded to GISAID). These 478 samples showed a majority of viruses from the B.1.1.7 / V1, lineage (79.1%). The other prevalent lineages were B.1.351 / V2 (7.9%), B.1.525 (4.4%), and B.1.214 (2.3%). Other lineages represented less than 2% of the samples (lineage P.1 / V3 being below 1%). Therefore, in the following, samples with only the N501Y mutation are treated as V1, samples with both mutations as V2, samples with only the E484K mutation as B.1.525, and samples with no mutation as wild types (B.1.214 being rare). ### Lineage spreading in France We used a multinomial log-linear model (the multinom function from the nnet R package) to analyze the lineage variable (V1 being the variant of reference). The explanatory variables were the patient age, individual location (in a hospital or not), and the interaction between the sampling region and date. The multinomial model revealed differences between lineages (Table 1). In terms of age, we find that older patients tend to be more infected by V1 (our reference), than by V2 or B.1.525. In hospital settings, we find an over-representation of V2 compared to V1. When analyzing region-specific temporal trends, we found that the odds of being infected by a wild type or a B.1.525 virus were either identical or lower than the one of being infected by V1. Conversely, in Île-de-France, and to a lesser extent in Hauts-de-France and Nouvelle-Aquitaine, we find that the risk of being infected by V2 instead of V1 significantly increases with time. View this table: [Table 1.](http://medrxiv.org/content/early/2021/05/17/2021.05.12.21257130/T2) Table 1. Odds ratios regarding the risk of lineage detection estimated using a multinomial log-linear model. The reference lineage in the analysis is B.1.1.7 (N501Y+/E484K-), i.e. the variant of concern V1. Only factor values with a least one significant effect are shown. P-values are obtained using a 2-tailed z test (symbols indicate p-values lower than 10−3 for \***|, than 10−2 for **, and than 0.05 for *). ### Transmission advantage Trends from the multinomial model should be treated with caution because of auto-correlation issues. Therefore, to investigate the temporal trends, we used the method described in [5] and fitted a logistic growth model to the fitted values of a generalized linear model (GLM) with 4 factors (sampling date, sampling department, and individual age) on the data sampled outside hospitals. For simplicity, we tested the transmission advantage of V2 compared to V1 and neglected the other strains in the analysis. We performed the analysis only in the Île-de-France, Hauts-de-France, and Nouvelle-Aquitaine regions. In the absence of adequate data to infer a serial interval in France, we use the one from [7]. We found a transmission advantage of 15.8% (95% confidence interval: [15.1,16.1]%) in Île-de-France and 17.1% (95%CI: [16.1,18.8]%) in Hauts-de-France (Figure 2). In Nouvelle-Aquitaine, the logistic growth was not significant. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/17/2021.05.12.21257130/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/05/17/2021.05.12.21257130/F2) Figure 2: Estimated proportion of V2 with respect to V1 (other lineages are excluded) in Île-de-France (blue) and Hauts-de-France (red). The triangles indicate the GLM-fitted values and the line is the output of the logistic growth model estimation. The top figures indicate the estimated transmission advantage of V2 with respect to V1 and the 95%-confidence interval. ## Discussion When analyzing variant-specific testing spanning from January to March 2021, we found that V1’s transmission advantage to the wild type was larger than that of V2 [8]. During April 2021, in at least two French regions, this trend appears to have shifted with V2 spreading more rapidly than V1. The V2 variant has known immune evasion properties [9,10]. Therefore, Île-de-France being one of the French regions the most impacted to date by the epidemic, it is sensible that we detect the shift in transmission advantage there because the proportion of individuals with natural immunity being high. Vaccination can also be expected to favor immune escape mutants [11] but the coverage is more homogeneous among French regions. These results call for more detailed analyses regarding the link between the transmission advantage of the V2 variant and the proportion of the population immunized in different French regions. Finally, note that France has entered a third national lock-down on April 3rd, which means that most of the tests analyzed here have been performed in a declining epidemic [12]. However, it is unlikely that the lock-down *per se* would affect differently the transmission of the variants. Given the high COVID-19 incidence in France and the uncertainties regarding the lock-down lifting [12], these results call for particular care to avoid another epidemic wave in the country. ## Data Availability Data will be made available upon manuscript publication. ## Ethics This study has been approved by the Institutional Review Board of the CHU of Montpellier and is registered at [ClinicalTrials.gov](http://ClinicalTrials.gov) with identifier [NCT04738331](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT04738331&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F17%2F2021.05.12.21257130.atom). ## Conflict of interest ‘None declared’. ## Funding Statement We thank the CNRs, the IRD, the ANR and the Région Occitanie (PHYEPI project) for funding. ## Acknowledgements We thank the ETE modelling team for discussion, as well as the CNRS, the IRD, the ANR, and the Région Occitanie for funding (PHYEPI grant). We thank the Infectious diseases department in Cerba laboratory, technical and IT teams for their technical assistance. * Received May 12, 2021. * Revision received May 12, 2021. * Accepted May 17, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Davies NG, Abbott S, Barnard RC, Jarvis CI, Kucharski AJ, Munday JD, et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science. 2021;eabg3055 2. 2.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. 2020;2020.12.21.20248640 3. 3.Faria NR, Mellan TA, Whittaker C, Claro IM, Candido D da S, Mishra S, et al. Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil. Science. 2021; doi: 10.1126/science.abh2644 [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiJzY2llbmNlLmFiaDI2NDR2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA1LzE3LzIwMjEuMDUuMTIuMjEyNTcxMzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. 4.The SARS-CoV-2 variant with lineage B.1.351 clusters investigation team. Linked transmission chains of imported SARS-CoV-2 variant B.1.351 across mainland France, January 2021. Eurosurveillance. 2021;26(13):2100333. 5. 5.Haim-Boukobza S, Roquebert B, Trombert-Paolantoni S, Lecorche E, Verdurme L, Foulongne V, et al. Detection of Rapid SARS-CoV-2 Variant Spread, France, January 26--February 16, 2021. Emerging Infectious Diseases. 2021;27(5):1496–9. 6. 6.Gaymard A, Bosetti P, Feri A, Destras G, Enouf V, Andronico 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. Eurosurveillance. 2021;26(9):2100133. 7. 7.Nishiura H, Linton NM, Akhmetzhanov AR. Serial interval of novel coronavirus (COVID-19) infections. International Journal of Infectious Diseases. 2020;93:284–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32145466&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F17%2F2021.05.12.21257130.atom) 8. 8.Roquebert B, Haim-Boukobza S, Trombert-Paolantoni S, Lecorche E, Verdurme L, Foulongne V, et al. SARS-CoV-2 variants of concern are associated with lower RT-PCR amplification cycles between January and March 2021 in France. medRxiv. 2021;2021.03.19.21253971. 9. 9.Zhou D, Dejnirattisai W, Supasa P, Liu C, Mentzer AJ, Ginn HM, et al. Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera. Cell. 2021;184(9):2348-2361.e6. 10. 10.Li Q, Nie J, Wu J, Zhang L, Ding R, Wang H, et al. SARS-CoV-2 501Y.V2 variants lack higher infectivity but do have immune escape. Cell. 2021;184(9):2362-2371.e9. 11. 11.Hacisuleyman E, Hale C, Saito Y, Blachere NE, Bergh M, Conlon EG, et al. Vaccine Breakthrough Infections with SARS-CoV-2 Variants. New England Journal of Medicine. 2021; doi: 10.1056/NEJMoa2105000 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2105000&link_type=DOI) 12. 12.Sofonea MT, Boennec C, Michalakis Y, Alizon S. Two waves and a high tide: the COVID-19 epidemic in France. Anaesthesia Critical Care and Pain Medicine. 2021; in press.