Hostname: page-component-78c5997874-fbnjt Total loading time: 0 Render date: 2024-11-05T05:11:29.922Z Has data issue: false hasContentIssue false

Representativeness of whole-genome sequencing approaches in England: the importance for understanding inequalities associated with SARS-CoV-2 infection

Published online by Cambridge University Press:  20 September 2023

Katherine A. Twohig
Affiliation:
COVID-19 National Epidemiology Cell, UKHSA, London, UK
Katie Harman
Affiliation:
COVID-19 Vaccines and Epidemiology Division, Public Health Programmes, Clinical and Public Health Group, UKHSA, London, UK
Asad Zaidi
Affiliation:
COVID-19 Vaccines and Epidemiology Division, Public Health Programmes, Clinical and Public Health Group, UKHSA, London, UK
Shirin Aliabadi
Affiliation:
COVID-19 National Epidemiology Cell, UKHSA, London, UK
Sophie G. Nash
Affiliation:
COVID-19 Vaccines and Epidemiology Division, Public Health Programmes, Clinical and Public Health Group, UKHSA, London, UK
Mary Sinnathamby
Affiliation:
COVID-19 Vaccines and Epidemiology Division, Public Health Programmes, Clinical and Public Health Group, UKHSA, London, UK
Ian Harrison
Affiliation:
Pathogen Genomics, Science Group, UKHSA, London, UK
Eileen Gallagher
Affiliation:
TARZET Division, Clinical and Emerging Infections Directorate, Clinical and Public Health Group, UKHSA, London, UK
Natalie Groves
Affiliation:
TARZET Division, Clinical and Emerging Infections Directorate, Clinical and Public Health Group, UKHSA, London, UK
Frank Schwach
Affiliation:
TARZET Division, Clinical and Emerging Infections Directorate, Clinical and Public Health Group, UKHSA, London, UK
Clare Pearson
Affiliation:
COVID-19 National Epidemiology Cell, UKHSA, London, UK
Alicia Thornton
Affiliation:
COVID-19 National Epidemiology Cell, UKHSA, London, UK
Richard Myers
Affiliation:
TARZET Division, Clinical and Emerging Infections Directorate, Clinical and Public Health Group, UKHSA, London, UK
Meera Chand
Affiliation:
TARZET Division, Clinical and Emerging Infections Directorate, Clinical and Public Health Group, UKHSA, London, UK
Simon Thelwall
Affiliation:
COVID-19 Vaccines and Epidemiology Division, Public Health Programmes, Clinical and Public Health Group, UKHSA, London, UK
Gavin Dabrera*
Affiliation:
COVID-19 Vaccines and Epidemiology Division, Public Health Programmes, Clinical and Public Health Group, UKHSA, London, UK
*
Corresponding author: Gavin Dabrera; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Whole-genome sequencing (WGS) information has played a crucial role in the SARS-CoV-2 (COVID-19) pandemic by providing evidence about variants to inform public health policy. The purpose of this study was to assess the representativeness of sequenced cases compared with all COVID-19 cases in England, between March 2020 and August 2021, by demographic and socio-economic characteristics, to evaluate the representativeness and utility of these data in epidemiological analyses. To achieve this, polymerase chain reaction (PCR)-confirmed COVID-19 cases were extracted from the national laboratory system and linked with WGS data. During the study period, over 10% of COVID-19 cases in England had WGS data available for epidemiological analysis. With sequencing capacity increasing throughout the period, sequencing representativeness compared to all reported COVID-19 cases increased over time, allowing for valuable epidemiological analyses using demographic and socio-economic characteristics, particularly during periods with emerging novel SARS-CoV-2 variants. This study demonstrates the comprehensiveness of England’s sequencing throughout the COVID-19 pandemic, rapidly detecting variants of concern, and enabling representative epidemiological analyses to inform policy.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© Crown Copyright - UK Health Security Agency, 2023. Published by Cambridge University Press

Introduction

SARS-CoV-2 has demonstrated itself to be a rapidly mutating virus, which has highlighted the value of large-scale, accessible, and timely genomic surveillance [Reference Robishaw1]. The COVID-19 pandemic has been the first pandemic where genomic technology has been widely available at such scale [Reference van Dorp LJH, Richard and Balloux2]. In the United Kingdom (UK), the COVID-19 Genomics UK Consortium (COG-UK) was established in April 2020 to provide large-scale, rapid whole-genome sequencing (WGS) for SARS-CoV-2 [3]. The UK has been at the forefront of global sequencing throughout the pandemic, with over two million genomes sequenced by February 2022 [4].

WGS has significantly contributed to the understanding of genomic diversity and evolution of the SARS-CoV-2 virus, as well as clinical and epidemiological characteristics of the infection. Through the rapid identification of novel variants, sequencing has been crucial in providing evidence to inform the implementation of public health policies, such as those that were established to manage the response to the alpha (B.1.1.7) variant and other lineages in early 2021 [Reference Vöhringer5, Reference du Plessis6]. To maximise population-level epidemiological insights, WGS must be as representative of the population as possible. Representativeness is beneficial both to ensure analyses are unbiased and to aid in global sharing of sequences, as per the World Health Organization recommendations [7].

It has been evidenced throughout the pandemic that COVID-19 has had varying impacts on population sub-groups, with older age, black, Asian, and minority ethnicity (BAME), and residence in areas of greater socio-economic deprivation, indicating greater risk of infection, as well as more severe outcomes including death [Reference Brown8, 9]. The capability and capacity to monitor health inequalities associated with emerging SARS-CoV-2 variants are essential to inform public health policy, but they rely on sequencing of cases providing sufficient information on groups in relation to key characteristics. It is therefore critical to understand the representativeness of sequenced cases in relation to confirmed cases overall.

We evaluated representativeness by assessing the proportion of sequenced COVID-19 cases in England, including changes over time, and by key demographic characteristics including sex, age, geography, indices of deprivation (IMD), ethnicity, and travel status.

Methods

Polymerase chain reaction (PCR)-confirmed COVID-19 cases in England reported to the national laboratory system, Second Generation Surveillance System (SGSS) [Reference Clare10], with specimen dates between 1 March 2020 and 31 August 2021, were linked with sequencing data uploaded to the Cloud Infrastructure for Big Data Microbial Bioinformatics (CLIMB) [Reference Nicholls11]. This linkage was based on specimen identifiers assigned at the diagnostic sites and submitted to Public Health England (PHE) either through secure file transfer or uploaded to CLIMB [12]. Linkage between patient data and sequencing results was performed securely within the PHE environment. Records were only linked if the sequence passed quality assurance thresholds, deeming it suitable for genomic analysis.

Key attributes of the case and test results were extracted from SGSS, including sex, age, geography of residence, and IMD, where quintile 1 represents the most deprived and 5 represents the least deprived [13], ethnicity, and reporting pillar of the first positive test. The test pillar represented the laboratory and reporting pathway of the positive result. Pillar 1 (P1) includes tests undertaken by public health, National Health Service (NHS), and privately contracted laboratories, and some targeted testing such as people in hospital and workplace screening. Pillar 2 (P2) tests were generally community-based and were reported into SGSS through NHS Digital Platforms. Information about recent international travel was also assessed, which was defined as arrival from outside of the UK within 14 days before the positive test date. This was derived from the linkage of five sources: arrival forms from recent travellers, contact tracing information, travel information included on test request forms, reports from the international arrival testing programme and questionnaires submitted from regional health protection teams.

Overall, the proportion of sequenced cases was assessed by demographic and epidemiological characteristics, with a focus on three time periods: (A) March to July 2020, (B) August 2020 to April 2021, and (C) May to August 2021. These intervals were calculated based on the specimen collection dates and included key changes in the epidemiology, reporting, or sequencing capacity of COVID-19 in England. Period A reflected the initial epidemic, after the first sporadic cases, and the introduction of the P2 testing pathway. Period B started from the beginning of the second wave of cases in England and through the winter months when sequencing capacity was further increased, particularly for P2 [14]. Period C began in the spring of 2021 following the winter peak and included the transition to greater sequencing capacity being taken over by public health laboratories instead of academic sequencing partners. The proportions of cases in each period and category were calculated for both total cases and sequenced cases, and then, a ratio of these proportions was calculated. Chi-squared tests were conducted to compare whether the distributions of characteristics were similar between total cases and sequenced cases (P < 0.001).

Results

Between 01 March 2020 and 31 August 2021, there were 5,810,945 PCR-confirmed COVID-19 cases reported in England, of which 688,203 (11.8%) were linked with quality-assessed WGS results from CLIMB. Through all three time periods, P1 had a lower testing volume than P2, yet the proportion of linked sequences varied over time and by testing pillar (Figure 1). The highest proportion of P1 cases linked to sequences occurred when case numbers were very low at the beginning of the pandemic, with a slow increase to a further peak in June 2021. During periods with higher case numbers, P1 sequencing volumes increased, as shown by the similar proportions of sequenced cases in April 2020, January 2021, and August 2021.

Figure 1. (a) PCR-confirmed COVID-19 cases reported in England, by pillar. (b) PCR-confirmed COVID-19 cases linked with quality-assessed whole-genome sequencing results from CLIMB. (c) 7-day rolling per cent of cases that were sequenced.

At the start of the study period (Period A), when cases were less than 100 per day in March 2020, approximately 80% of cases were sequenced (Figure 1). P2 testing began in mid-April 2020, and sequencing for P2 tests started in June 2020. The proportion of P2 cases sequenced was highest when case numbers were low, particularly in June 2020, when the proportion occasionally exceeded 50%, and, in the spring of 2021, when the highest proportion throughout the study was observed at 73.5% in May 2021 (Figure 1).

From August 2020 (Period B), case numbers began to rise, particularly in P2, which had increased testing capacity. Although the absolute number of sequenced cases increased, they reflected a smaller proportion of the overall cases. P1 sequencing volumes peaked in the middle of January 2021 at 12.8%, whereas P2 sequencing peaked at 13.0% in the end of January 2021. The proportion sequenced then rose again in line with decreases in case numbers, with P1 reaching a peak of 43.3% sequenced in mid-June 2021, and P2 reaching 73.5% at the end of May. During Period C, P2 sequencing capacity greatly increased, and the proportion of cases sequenced remained high following this rise in capacity. In August 2021, an average of 24.3% of P1 cases were sequenced, compared with 16.8% for P2.

Overall, the ratio of the proportion of sequenced cases to the proportion of total cases shows that sequencing results were broadly representative of the underlying case populations (Table 1); however, there were some over-represented groups. For all time periods, there were a higher proportion of P1 cases that were sequenced, compared with the proportion of P1 among total cases.

Table 1. Demographic and epidemiological characteristic breakdowns of overall testing and sequencing, 01 March 2020 to 31 August 2021

a Ratio is defined as the proportion of total cases to the proportion of sequenced cases.

b P-values <0.001, obtained using chi-squared tests for the stratum.

Throughout all time periods, sequenced cases were broadly representative of total cases by age group and sex, with the ratio of total to sequenced cases being close to one throughout, aside from Period A, which saw slightly larger proportions of younger age groups sequenced.

There was strong statistical evidence that the proportion of cases that were sequenced varied by geography (P < 0.001). The proportion of cases that were sequenced was highest in the East of England (21.7%) and lowest in the West Midlands (4.1%) in Period A, with the ratio of total to sequenced cases being 2.1 and 0.4, respectively. In Period B, the proportion sequenced was highest in the North West (20.2%) and lowest in the South West (4.1%), with ratios of 1.3 and 0.7, respectively. In Period C, there was the least amount of geographic variation, with a range of ratios from 0.8 (lowest proportion sequenced in the North East with 5.3%) and 1.1 (highest sequenced in the North West with 18.7%).

By IMD quintile, the proportion of cases that were sequenced was highest in the most deprived quintile, IMD-1, (27.6%) and lowest in the least deprived IMD-5 (14.7%) in Period A. The decreasing trend in proportions sequenced across IMD quintiles remained through Periods B and C, with the least differences observed in Period C ranging from IMD-1 (22.5%) to IMD-5 (18.1%).

Those of Asian ethnicity were more highly represented in sequenced cases during Period A, with a ratio of total to sequenced cases of 1.4, followed by those of mixed ethnicity with a ratio of 1.2. There were no major disproportionalities in Period B, but again in period C those of Asian ethnicity were more highly represented with a ratio of 1.2.

Travel data were not available for Period A, and the proportion of all cases that were travel-related during Periods B and C was 0.5% and 2.3% of all cases, respectively (Table 1). For the two periods when travel data were available, there were a higher proportion of travel-related cases that were sequenced compared with the proportion of all cases that were travel-related, with ratios of 3.6 and 1.4, respectively.

Discussion

More than 10% of all PCR-confirmed COVID-19 cases in England between March 2020 and August 2021 had WGS data available for national epidemiological analysis. As a novel pathogen with initially few cases, most of the early SARS-CoV-2 samples were sequenced to gain early insights into the virus. The proportion of sequenced cases was inversely related to the number of cases per day, with the highest proportions sequenced during periods of lower prevalence, likely reflecting sequencing capacity within the laboratory. As the pandemic progressed and sequencing capacity increased, sequenced cases more closely reflected the demographics of total cases and comprised more P2 (community) samples, which captured population-level assessment of the emergence and spread of new variants, providing useful information on variant prevalence. There were some sequencing disproportionalities by geographic location, likely a result of variation in sampling methodology, and operational or logistical considerations. However, the overall disproportionality between regions decreased over time from periods A to C. Evidence of targeted sequencing was also observed, including an emphasis on international travellers and P1 testing.

Sequencing was broadly representative of total cases when broken down by age group. However, variations were observed in ethnicity and deprivation, which were key health inequalities highlighted in the pandemic [9]. In particular, cases of Asian ethnicity were more likely to be sequenced than people of other ethnicities in Periods A and C and those of mixed ethnicity were more highly represented in Period A. Higher representativeness of different ethnic groups during certain periods was important in developing insights for minority at-risk population groups, including related inequalities, as demonstrated by analysis during the emergence of the delta variant [Reference Leeman, Campos-Matos and Dabrera15]. Overall, the moderate over-representation of some ethnic groups in the study findings informs our understanding of the inequitable impact of COVID-19 in specific communities, but further explanatory work, such as analysis of hospitalisation data, would be needed to reduce any impact of selection bias and more comprehensively assess disproportionate burden.

Throughout the pandemic, there has been a higher burden of disease in more deprived residential areas [Reference Leeman, Campos-Matos and Dabrera15]. Early sequencing coverage of cases from more deprived geographies in this study provides further insight into these trends, particularly in relation to the emergence of new variants. Overall, the representativeness of sequencing improved over time. During Period A, while the proportion of cases with sequencing data were broadly similar, those residing in more deprived areas (quintile 1) were more likely to be sequenced, which persisted during Period B despite the significant increase in sequencing capacity. By Period C, proportions sequenced aligned to be generally representative of the total cases in each IMD quintile. A higher proportion of sequenced cases from persons of Asian ethnicity and those from more deprived areas suggest that ethnicity and deprivation may have intersected with other characteristics that affected sample selection, such as more severe infections being hospitalised or targeted geographic sequencing related to emerging variants [Reference Sze16]. However, due to incomplete information on why people sought testing and which route they may have accessed for their initial positive test (P1 or P2), we cannot fully assess the relationship between severity and sample selection.

Some of the present study findings, particularly around the testing pillar, reflect the establishment and expansion of WGS. The testing pillar determined the journey of a specimen and affected the probability of being selected for sequencing. Examples of this include Period A, when a disproportionately greater number of cases in East of England were sequenced. This may have been an artefact reflecting the geographical bases of laboratories and the increasing sequencing capacity being established at this time. This was also observed in Period B, when hospitalisation rates were high in the northern regions and the sequencing through P1 was over-represented [17]. Later, fluctuations in geographical representativeness may also reflect periods when there was surge testing in certain geographical areas, increasing case ascertainment [Reference Iacobucci18].

Overall, national co-ordination of sequencing sample selection was successful in ensuring that it was representative of COVID-19 cases in England. However, there may have been some bias introduced when different sampling strategies were employed. During the early stages of the pandemic, despite high proportions of cases being sequenced, access to PCR testing was limited and therefore confirmed cases would not have been representative of the true burden of infection in the population. Another data limitation is that pillar assignment reflects the testing route of an individual’s earliest positive specimen; therefore, in instances where someone tested positive across both pillars and their subsequent test was sequenced, the pillar assignment in this analysis would not match the sequenced result. Lastly, the representativeness of overall cases might have been affected by the use of lateral flow tests (LFTs), which were distributed to local authorities in November 2020 [19] and were available to the general public from April 2021 [20]. Despite guidance that positive LFTs should be followed by confirmatory PCR testing through the majority of the study period [21], this was not always followed and WGS is not possible without PCR samples. Data on the number of tests sent for sequencing and subsequent failure rate were not available, and differences in sequencing attempts by demographic group are not known.

This study is greatly strengthened by access to high-quality data on COVID-19 cases and consequent sequencing, allowing for an in-depth nationwide analysis. This analysis demonstrates that sequenced cases of COVID-19 between March 2020 and August 2021 in England were generally representative of all cases by key socio-economic characteristics, such as the most deprived who were disproportionately affected by the pandemic, providing important utility for these sequencing data. Future work will be needed to explore sequencing trends in later time periods and to consider the effects on representativeness from changes to SARS-CoV-2 testing approaches in England [Reference Halford22], particularly in relation to deprivation and ethnicity. This study has provided evidence to support the use of sequencing data for rigorous epidemiological analyses, including, but not limited to, assessing variant severity, vaccine effectiveness, and household transmission, which have been crucial in understanding health inequalities and informing the national COVID-19 pandemic response [Reference Allen23Reference Flannagan27].

Data availability statement

The individual-level nature of the data used risk individuals being identified, or being able to self-identify, if the data are released publicly. Requests for access to these non-publicly available data should be directed to UKHSA.

Acknowledgements

The authors acknowledge Ross Harris for his expert statistical advice.

Author contribution

Data curation: A.Z., S.G.N., I.H., A.T., C.P., E.G., F.S., K.A.T., M.S., N.G., S.A.; Writing – review & editing: A.Z., S.G.N., I.H., A.T., C.P., E.G., F.S., M.C., G.D., K.A.T., M.S., N.G., R.M., S.A., S.T., K.H.; Resources: M.C., G.D., R.M.; Conceptualization: G.D., K.A.T., M.S.; Supervision: G.D., S.T.; Formal analysis: K.A.T., K.H.; Writing – original draft: K.A.T.; Validation: K.H.

Financial support

This work was performed as part of UKHSA’s responsibility to monitor COVID-19.

Competing interest

G.D. declares that his employer’s predecessor organisation, Public Health England, received funding from GlaxoSmithKline for a previous research project related to influenza antiviral treatment. This preceded and had no relation to COVID-19, and G.D. had no role in and received no funding from the project. All other authors report no potential conflicts.

Ethical standard

UKHSA has legal permission, provided by Regulation 3 of the Health Service (Control of Patient Information) Regulations 2002 to process confidential patient information under Sections 3(i) (a) to (c), 3(i) (d) (i) and (ii), and 3(iii) as part of its outbreak response activities. This study falls within the research activities approved by the UKHSA Research Ethics and Governance Group.

Footnotes

K.A.T. and K.H. are joint first authors.

References

Robishaw, JD, et al. (2021) Genomic surveillance to combat COVID-19: Challenges and opportunities. Lancet Microbe 2(9), e481e484.CrossRefGoogle ScholarPubMed
van Dorp LJH, CJ, Richard, D and Balloux, F (2021) COVID-19, the first pandemic in the post-genomic era. Current Opinion in Virology 50, 4048.CrossRefGoogle Scholar
(COG-UK) C-GUC (2020) An integrated national scale SARS-CoV-2 genomic surveillance network. Lancet Microbe 1(3), e99e100.CrossRefGoogle Scholar
(UKHSA) UHSA. UK completes over 2 million SARS-CoV-2 whole genome sequences. GOV.UK (2022).Google Scholar
Vöhringer, HS, et al. (2021) Genomic reconstruction of the SARS-CoV-2 epidemic in England. Nature 600(7889), 506511.CrossRefGoogle ScholarPubMed
du Plessis, L, et al. (2021) Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK. Science 371(6530), 708712.CrossRefGoogle ScholarPubMed
WHO. Genomic sequencing of SARS-CoV-2: a guide to implementation for maximum impact on public health (2021).Google Scholar
Brown, AE, et al. (2021) Epidemiology of confirmed COVID-19 deaths in adults, England, March-December 2020. Emerging Infectious Diseases 27(5), 14681471.CrossRefGoogle ScholarPubMed
PHE. Disparities in the risk and outcomes of COVID-19 (2020).Google Scholar
Clare, T, et al. (2021) Timeliness and completeness of laboratory-based surveillance of COVID-19 cases in England. Public Health 194, 163166.CrossRefGoogle ScholarPubMed
Nicholls, SM, et al. (2021) CLIMB-COVID: Continuous integration supporting decentralised sequencing for SARS-CoV-2 genomic surveillance. Genome Biology 22(1), 196.CrossRefGoogle ScholarPubMed
(COG-UK) C-GUC. How do we manage data about people’s SARS-CoV-2 samples? (2021).Google Scholar
Ministry of Housing CLG. English indices of deprivation 2019 (2019).Google Scholar
(COG-UK) C-GUC. £12.2 million boost for SARS-CoV-2 real-time genomic surveillance (2020).Google Scholar
Leeman, D, Campos-Matos, I and Dabrera, G (2022) Inequalities associated with emergence of Delta SARS-CoV-2 variant of concern (B.1.617.2) in England: Awareness for future variants. Public Health 205, e14e15.CrossRefGoogle ScholarPubMed
Sze, S, et al. (2020) Ethnicity and clinical outcomes in COVID-19: A systematic review and meta-analysis. EClinicalMedicine 29, 100630.CrossRefGoogle ScholarPubMed
PHE. Weekly national Influenza and COVID19 surveillance report: Week 41 report (up to week 40 data) (2020).Google Scholar
Iacobucci, G (2021) Covid-19: Surge testing for the South Africa variant begins in England. BMJ 372, n305.CrossRefGoogle ScholarPubMed
DHSC. More rapid COVID-19 tests to be rolled out across England (2020).Google Scholar
DHSC. Twice weekly rapid testing to be available to everyone in England (2021).Google Scholar
UKHSA. COVID-19: Reintroducing confirmatory PCR testing (2021).Google Scholar
Halford, F, et al. (2022) Variation in reported SARS-CoV-2 cases after testing policy changes. The Lancet Infectious Diseases 22(10), 14181419.CrossRefGoogle ScholarPubMed
Allen, H, et al. (2022) Household transmission of COVID-19 cases associated with SARS-CoV-2 delta variant (B.1.617.2): National case-control study. Lancet Regional Health Europe 12, 100252.CrossRefGoogle ScholarPubMed
Twohig, KA, et al. (2022) Hospital admission and emergency care attendance risk for SARS-CoV-2 delta (B.1.617.2) compared with alpha (B.1.1.7) variants of concern: A cohort study. The Lancet Infectious Diseases 22(1), 3542.CrossRefGoogle ScholarPubMed
Lopez Bernal, J, et al. (2021) Effectiveness of Covid-19 vaccines against the B.1.617.2 (Delta) variant. The New England Journal of Medicine 385(7), 585594.CrossRefGoogle ScholarPubMed
Dabrera, G, et al. (2020 ) Assessment of mortality and hospital admissions associated with confirmed infection with SARS-CoV-2 alpha variant: A matched cohort and time-to-event analysis, England, October to December 2020. Euro Surveillance 27(20), 2100377.Google Scholar
Flannagan, J, et al. (2022) Increased residential clustering of COVID-19 cases associated with SARS-CoV-2 variant of concern B.1.1.7. Epidemiology 33(4), 465469.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. (a) PCR-confirmed COVID-19 cases reported in England, by pillar. (b) PCR-confirmed COVID-19 cases linked with quality-assessed whole-genome sequencing results from CLIMB. (c) 7-day rolling per cent of cases that were sequenced.

Figure 1

Table 1. Demographic and epidemiological characteristic breakdowns of overall testing and sequencing, 01 March 2020 to 31 August 2021