SARS-CoV-2 genome sequence prevalence and growth rate update: 1 May 2024
Updated 8 August 2024
Applies to England
Please note that these briefings will be moving from monthly to quarterly publication until further notice.
Moving forward, the aim of UK genomic surveillance of SARS-CoV-2 will be to provide variant prevalence assessments at suitable time points to inform vaccine adjustment and selection and to contribute to the global surveillance picture to ensure significant changes such as non-Omicron variants are detected rapidly. In line with these changes, the samples that are available for sequencing will reduce in number, and these briefings will be moving to quarterly reporting until further notice.
Variant prevalence
Testing policy and sequencing should be considered when interpreting variant data. Information about surveillance systems for England are reported by UKHSA in the National Influenza and COVID-19 surveillance report.
The prevalence of lineages amongst UK sequences by Phylogenetic Assignment of Named Global Outbreak Lineages (Pangolin) designation is presented in Figure 1. Lineages are shown if there are more than or equal to 5,000 sequences since 29 January 2024 or if they represent more than or equal to 1% of sequences within a single week over the last 6 weeks. Lineages that do not meet these criteria are combined with their parent lineage (for example, BA.2.4 is combined with BA.2).
The lineages have been assigned using the accurate Ultrafast Sample placement on Existing tRee (UShER) mode and version 1.26 of the Pangolin data.
Figure 2 shows the relationship between lineages shown in Figure 1. The figure details the hierarchical relationship of the lineages from left to right, starting with B.1.1.529 through 14 sublineage levels, including aliases and the parent lineages for any recombinants included in the prevalence figures (for example, BJ.1 and BM.1.1.1 are the parent lineages for the XBB recombinant). Lineages that have been combined with a lineage ancestral to B.1.1.529 in Figure 1 are not included in Figure 2. The percentage of each lineage in the most recent week of data in Figure 1 is shown next to the relevant lineage labels in Figure 2. Where a sublineage has a percentage in Figure 2, it is excluded from the percentage given for any parent lineages. Any lineages that are not present in the most recent week of genomic data will not have a percentage in Figure 2. The colours in Figure 2 correspond to the lineage colours in Figure 1.
Figure 1. Prevalence of Pangolin lineages in the UK sequence data with a specimen date from week beginning 29 January 2024 to week beginning 8 April 2024, as of 26 April 2024
The total number of valid sequence results per week is shown by the black line. The ‘Other’ category in this plot contains all lineages that do not meet the relevant criteria after combining smaller sub-lineages. ‘Unassigned’ are sequences that could not be assigned a lineage by Pangolin. Lineages present in at least 2% of sequences in the most recent week are labelled to the right of the plot.
Figure 2. Sankey diagram showing the relationship between Pangolin lineages observed in UK sequence data since 29 January 2024. Data shown as of 26 April 2024. Proportions are given for lineages that are observed in sequences with a specimen date between 8 April 2024 and 14 April 2024
Lineage colours match those in Figure 1
Variant modelling
Two models are currently used to estimate the growth advantage of emerging lineages: a logistic regression generalised linear model (GLM) and a generalised additive model (GAM). Models are fit to a geographically stratified sample of Pillar 1 cases to ensure that relative growth rates are estimated in relation to a local set of co-circulating lineages. Tests associated with travel are excluded. A full description of methods can be found in the variant technical briefing series. In recent months there are fewer tests to model than earlier in the year (see Figure 1). This is due to reduced sampling effort and lower prevalence. Moreover, the proliferation of lineages to monitor means that sample sizes for specific lineages can be small. Uncertainty in our modelled relative growth rates is therefore increased, which is reflected in larger confidence intervals on the estimates.
We aim to select lineages and/or groups of lineages that are both specific enough to pick up on emerging signals but broad enough to maximise statistical power. Any lineage that has made up more than 1.5% of the total samples and with at least 50 sequenced cases within 6 weeks of the most recent specimen date is modelled separately. Lineages that do not meet these criteria will be added to a group with their closest parent lineage but will not be aggregated further than one. If the condition cannot be met for a particular lineage (for example, it does not have a close parent lineage at high enough prevalence) it will not be modelled. Lineages with a different high-level parent will never be aggregated together (for example, we will not aggregate BA.2 and BA.5 to B.1.1.529). Unassigned lineages are excluded from this analysis. Note that these aggregations will often be broader than in the prevalence plot presented in Figure 1. The relative growth rate of broad lineage classes (for example, parent lineages that include many distant child lineages) will be less informative than explicit modelling of specific sublineages. This is why we set a limit of one level when aggregating child lineages into parents. Methods of lineage collapsing for the growth rate analysis are therefore still being refined.
Growth rates were based on sequences sampled through Pillar 1 testing (primarily positive tests conducted in hospital) in England (Table 1). The sampling range for both the logistic regression GLM and GAM is from 11 November 2023 to 17 April 2024. The model fit for any lineage with a positive growth rate advantage (with 95% confidence intervals (CIs) that do not cross zero) are shown in Figure 3. The lineages that were estimated to have a positive growth rate advantage with reasonable uncertainty (CIs not below zero in the GAM) were JN.1.11.1 (38.92%, GAM), KP.2 (25.09%, GAM) and JN.1.7.1 (24.91%, GAM).
Table 1. Growth rate (GR) of English sequence lineages as of 17 April 2024†
Lineage* | Lineage Group Composition** | Pillar 1 Sample Size*** | Weekly growth rate advantage (GAM) | Estimated prevalence¥ (GAM) | Weekly growth rate advantage (GLM) |
---|---|---|---|---|---|
JN.1.11.1 (BA.2.86.1.1.11.1) | JN.1.11.1 (92.47%); KP.1 (7.53%) | 93 | 38.92% (95% CI: 30.93 to 46.92) | 8.99% (95% CI: 6.7 to 11.95) | 71.64% (95% CI: 18.93 to 124.34) |
KP.2 (BA.2.86.1.1.11.1.2) | KP.2 (100%) | 165 | 25.09% (95% CI: 1.25 to 48.92) | 14.79% (95% CI: 11.5 to 18.83) | 20.95% (95% CI: -3.74 to 45.65) |
JN.1.7.1 (BA.2.86.1.1.7.1) | JN.1.7.1 (100%) | 52 | 24.91% (95% CI: 19.11 to 30.71) | 4.26% (95% CI: 2.86 to 6.28) | 22.68% (95% CI: -35.74 to 81.1) |
XDK | XDK (75%); XDK.1 (25%) | 64 | 12.42% (95% CI: -4.63 to 29.48) | 4.29% (95% CI: 2.92 to 6.26) | 44.45% (95% CI: 3.53 to 85.38) |
JN.1.16 (BA.2.86.1.1.16) | JN.1.16 (100%) | 152 | 8.16% (95% CI: -5.06 to 21.38) | 10.87% (95% CI: 8.26 to 14.17) | 1.15% (95% CI: -24.94 to 27.24) |
JN.1.4 (BA.2.86.1.1.4) | JN.1.4 (94.55%); JN.1.4.3 (4.46%); JN.1.4.2 (0.5%); JN.1.4.1 (0.5%) | 202 | 2.77% (95% CI: -21.52 to 27.07) | 7.87% (95% CI: 5.86 to 10.48) | -3.39% (95% CI: -28.95 to 22.17) |
JN.1.18 (BA.2.86.1.1.18) | JN.1.18 (100%) | 116 | 2.39% (95% CI: -8.42 to 13.2) | 5.88% (95% CI: 4.25 to 8.07) | -7.05% (95% CI: -37.02 to 22.92) |
JN.1.7 (BA.2.86.1.1.7) | JN.1.7 (98.4%); JN.1.7.2 (1.6%) | 312 | -9.62% (95% CI: -22.56 to 3.32) | 15.42% (95% CI: 12.41 to 19.01) | 0.62% (95% CI: -17.43 to 18.67) |
JN.1 (BA.2.86.1.1) | JN.1 (78.95%); JN.1.1 (5.11%); JN.1.5 (4.49%); JN.1.8 (2.48%); JN.1.19 (2.48%)… | 646 | -29.74% (95% CI: -42.44 to -17.04) | 18.85% (95% CI: 15.86 to 22.25) | -21.81% (95% CI: -36.9 to -6.73) |
*Listed parent lineages include all sub-lineages, other than those explicitly modelled.
** The top 5 contributing lineages to the modelled group in the most recent 6 weeks (6 March 2024 to 17 April 2024). More than 5 sublineages are indicated by “…”
*** Sample size is for Pillar 1 samples in England in the most recent 6 weeks (6 March 2024 to 17 April 2024).
¥ Estimated prevalence for the 14 March 2024.
† Sampling range for both logistic regression and generalised additive models (GAM) is from 1 November 2023 to 17 April 2024.
CI = confidence intervals
Figure 3. Modelled prevalence of lineage groups with a growth rate advantage over other circulating lineages
The black line shows the central estimate and blue shaded regions the 95% confidence intervals. Points show the national level proportions, with the size being indicative of the sample size for that particular lineage. The grey portion of the ribbon denotes that this period of time is likely to be backfilled with more sequenced cases, making proportions unreliable.
Sources and acknowledgments
Data sources
Data used in this investigation is derived from the Cloud Infrastructure for Microbial Bioinformatics (CLIMB) data set, the UK Health Security Agency (UKHSA) genomic programme data set and the UKHSA Second Generation Surveillance System. International data has been gathered from the Global Initiative on Sharing All Influenza Data (GISAID).
Authors of this report
UKHSA Genomics Public Health Analysis Team
UKHSA Infectious Disease Modelling Team
UKHSA TARZET Technical Secretariat