Systematic review of wastewater surveillance of antimicrobial resistance in human populations

OBJECTIVES

We systematically reviewed studies using wastewater for AMR surveillance in human populations, to determine: (i) evidence of concordance between wastewater-human AMR prevalence estimates, and (ii) methodological approaches which optimised identifying such an association, and which could be recommended as standard. We used Lin's concordance correlation coefficient (CCC) to quantify concordance between AMR prevalence estimates in wastewater and human compartments (where CCC = 1 reflects perfect concordance), and logistic regression to identify study features (e.g. sampling methods) associated with high agreement studies (defined as >70% of within-study wastewater-human AMR prevalence comparisons within ±10%).

RESULTS

Of 8,867 records and 441 full-text methods reviewed, 33 studies were included. AMR prevalence data was extractable from 24 studies conducting phenotypic-only (n = 7), genotypic-only (n = 1) or combined (n = 16) AMR detection. Overall concordance of wastewater-human AMR prevalence estimates was reasonably high for both phenotypic (CCC = 0.85 [95% CI 0.8-0.89]) and genotypic approaches (CCC = 0.88 (95% CI 0.84-0.9)) despite diverse study designs, bacterial species investigated and phenotypic/genotypic targets. No significant relationships between methodological approaches and high agreement studies were identified using logistic regression; however, this was limited by inconsistent reporting of study features, significant heterogeneity in approaches and limited sample size. Based on a secondary, descriptive synthesis, studies conducting composite sampling of wastewater influent, longitudinal sampling >12 months, and time-/location-matched sampling of wastewater and human compartments generally had higher agreement.

CONCLUSION

Wastewater-based surveillance of AMR appears promising, with high overall concordance between wastewater and human AMR prevalence estimates in studies irrespective of heterogenous approaches. However, our review suggests future work would benefit from: time-/location-matched sampling of wastewater and human populations, composite sampling of influent, and sampling >12 months for longitudinal studies. Further research and clear and consistent reporting of study methods is required to identify optimal practice.

The 2021 WHO catalogue of Mycobacterium tuberculosis complex mutations associated with drug resistance: A genotypic analysis

Background

Molecular diagnostics are considered the most promising route to achieving rapid, universal drug susceptibility testing for Mycobacterium tuberculosiscomplex (MTBC). We aimed to generate a WHO endorsed catalogue of mutations to serve as a global standard for interpreting molecular information for drug resistance prediction.

Methods

A candidate gene approach was used to identify mutations as associated with resistance, or consistent with susceptibility, for 13 WHO endorsed anti-tuberculosis drugs. 38,215 MTBC isolates with paired whole-genome sequencing and phenotypic drug susceptibility testing data were amassed from 45 countries. For each mutation, a contingency table of binary phenotypes and presence or absence of the mutation computed positive predictive value, and Fisher's exact tests generated odds ratios and Benjamini-Hochberg corrected p-values. Mutations were graded as Associated with Resistance if present in at least 5 isolates, if the odds ratio was >1 with a statistically significant corrected p-value, and if the lower bound of the 95% confidence interval on the positive predictive value for phenotypic resistance was >25%. A series of expert rules were applied for final confidence grading of each mutation.

Findings

15,667 associations were computed for 13,211 unique mutations linked to one or more drugs. 1,149/15,667 (7·3%) mutations were classified as associated with phenotypic resistance and 107/15,667 (0·7%) were deemed consistent with susceptibility. For rifampicin, isoniazid, ethambutol, fluoroquinolones, and streptomycin, the mutations' pooled sensitivity was >80%. Specificity was over 95% for all drugs except ethionamide (91·4%), moxifloxacin (91·6%) and ethambutol (93·3%). Only two resistance mutations were classified for bedaquiline, delamanid, clofazimine, and linezolid as prevalence of phenotypic resistance was low for these drugs.

Interpretation

This first WHO endorsed catalogue of molecular targets for MTBC drug susceptibility testing provides a global standard for resistance interpretation. Its existence should encourage the implementation of molecular diagnostics by National Tuberculosis Programmes.

Funding

UNITAID, Wellcome, MRC, BMGF.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infectivity by Viral Load, S Gene Variants and Demographic Factors, and the Utility of Lateral Flow Devices to Prevent Transmission

BACKGROUND

How severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity varies with viral load is incompletely understood. Whether rapid point-of-care antigen lateral flow devices (LFDs) detect most potential transmission sources despite imperfect clinical sensitivity is unknown.

METHODS

We combined SARS-CoV-2 testing and contact tracing data from England between 1 September 2020 and 28 February 2021. We used multivariable logistic regression to investigate relationships between polymerase chain reaction (PCR)-confirmed infection in contacts of community-diagnosed cases and index case viral load, S gene target failure (proxy for B.1.1.7 infection), demographics, SARS-CoV-2 incidence, social deprivation, and contact event type. We used LFD performance to simulate the proportion of cases with a PCR-positive contact expected to be detected using 1 of 4 LFDs.

RESULTS

In total, 231 498/2 474 066 (9%) contacts of 1 064 004 index cases tested PCR-positive. PCR-positive results in contacts independently increased with higher case viral loads (lower cycle threshold [Ct] values), for example, 11.7% (95% confidence interval [CI] 11.5-12.0%) at Ct = 15 and 4.5% (95% CI 4.4-4.6%) at Ct = 30. B.1.1.7 infection increased PCR-positive results by ~50%, (eg, 1.55-fold, 95% CI 1.49-1.61, at Ct = 20). PCR-positive results were most common in household contacts (at Ct = 20.1, 8.7% [95% CI 8.6-8.9%]), followed by household visitors (7.1% [95% CI 6.8-7.3%]), contacts at events/activities (5.2% [95% CI 4.9-5.4%]), work/education (4.6% [95% CI 4.4-4.8%]), and least common after outdoor contact (2.9% [95% CI 2.3-3.8%]). Contacts of children were the least likely to test positive, particularly following contact outdoors or at work/education. The most and least sensitive LFDs would detect 89.5% (95% CI 89.4-89.6%) and 83.0% (95% CI 82.8-83.1%) of cases with PCR-positive contacts, respectively.

CONCLUSIONS

SARS-CoV-2 infectivity varies by case viral load, contact event type, and age. Those with high viral loads are the most infectious. B.1.1.7 increased transmission by ~50%. The best performing LFDs detect most infectious cases.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infectivity by Viral Load, S Gene Variants and Demographic Factors, and the Utility of Lateral Flow Devices to Prevent Transmission

BACKGROUND

How severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity varies with viral load is incompletely understood. Whether rapid point-of-care antigen lateral flow devices (LFDs) detect most potential transmission sources despite imperfect clinical sensitivity is unknown.

METHODS

We combined SARS-CoV-2 testing and contact tracing data from England between 1 September 2020 and 28 February 2021. We used multivariable logistic regression to investigate relationships between polymerase chain reaction (PCR)-confirmed infection in contacts of community-diagnosed cases and index case viral load, S gene target failure (proxy for B.1.1.7 infection), demographics, SARS-CoV-2 incidence, social deprivation, and contact event type. We used LFD performance to simulate the proportion of cases with a PCR-positive contact expected to be detected using 1 of 4 LFDs.

RESULTS

In total, 231 498/2 474 066 (9%) contacts of 1 064 004 index cases tested PCR-positive. PCR-positive results in contacts independently increased with higher case viral loads (lower cycle threshold [Ct] values), for example, 11.7% (95% confidence interval [CI] 11.5-12.0%) at Ct = 15 and 4.5% (95% CI 4.4-4.6%) at Ct = 30. B.1.1.7 infection increased PCR-positive results by ~50%, (eg, 1.55-fold, 95% CI 1.49-1.61, at Ct = 20). PCR-positive results were most common in household contacts (at Ct = 20.1, 8.7% [95% CI 8.6-8.9%]), followed by household visitors (7.1% [95% CI 6.8-7.3%]), contacts at events/activities (5.2% [95% CI 4.9-5.4%]), work/education (4.6% [95% CI 4.4-4.8%]), and least common after outdoor contact (2.9% [95% CI 2.3-3.8%]). Contacts of children were the least likely to test positive, particularly following contact outdoors or at work/education. The most and least sensitive LFDs would detect 89.5% (95% CI 89.4-89.6%) and 83.0% (95% CI 82.8-83.1%) of cases with PCR-positive contacts, respectively.

CONCLUSIONS

SARS-CoV-2 infectivity varies by case viral load, contact event type, and age. Those with high viral loads are the most infectious. B.1.1.7 increased transmission by ~50%. The best performing LFDs detect most infectious cases.

Time of Day of Vaccination Affects SARS-CoV-2 Antibody Responses in an Observational Study of Health Care Workers

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a global crisis with unprecedented challenges for public health. Vaccinations against SARS-CoV-2 have slowed the incidence of new infections and reduced disease severity. As the time of day of vaccination has been reported to influence host immune responses to multiple pathogens, we quantified the influence of SARS-CoV-2 vaccination time, vaccine type, participant age, sex, and days post-vaccination on anti-Spike antibody responses in health care workers. The magnitude of the anti-Spike antibody response is associated with the time of day of vaccination, vaccine type, participant age, sex, and days post-vaccination. These results may be relevant for optimising SARS-CoV-2 vaccine efficacy.

An end-to-end heterogeneous graph attention network for Mycobacterium tuberculosis drug-resistance prediction

Antimicrobial resistance (AMR) poses a threat to global public health. To mitigate the impacts of AMR, it is important to identify the molecular mechanisms of AMR and thereby determine optimal therapy as early as possible. Conventional machine learning-based drug-resistance analyses assume genetic variations to be homogeneous, thus not distinguishing between coding and intergenic sequences. In this study, we represent genetic data from Mycobacterium tuberculosis as a graph, and then adopt a deep graph learning method—heterogeneous graph attention network (‘HGAT–AMR’)—to predict anti-tuberculosis (TB) drug resistance. The HGAT–AMR model is able to accommodate incomplete phenotypic profiles, as well as provide ‘attention scores’ of genes and single nucleotide polymorphisms (SNPs) both at a population level and for individual samples. These scores encode the inputs, which the model is ‘paying attention to’ in making its drug resistance predictions. The results show that the proposed model generated the best area under the receiver operating characteristic (AUROC) for isoniazid and rifampicin (98.53 and 99.10%), the best sensitivity for three first-line drugs (94.91% for isoniazid, 96.60% for ethambutol and 90.63% for pyrazinamide), and maintained performance when the data were associated with incomplete phenotypes (i.e. for those isolates for which phenotypic data for some drugs were missing). We also demonstrate that the model successfully identifies genes and SNPs associated with drug resistance, mitigating the impact of resistance profile while considering particular drug resistance, which is consistent with domain knowledge.

Antimicrobial resistance determinants are associated with Staphylococcus aureus bacteraemia and adaptation to the healthcare environment: a bacterial genome-wide association study

Staphylococcus aureus is a major bacterial pathogen in humans, and a dominant cause of severe bloodstream infections. Globally, antimicrobial resistance (AMR) in S. aureus remains challenging. While human risk factors for infection have been defined, contradictory evidence exists for the role of bacterial genomic variation in S. aureus disease. To investigate the contribution of bacterial lineage and genomic variation to the development of bloodstream infection, we undertook a genome-wide association study comparing bacteria from 1017 individuals with bacteraemia to 984 adults with asymptomatic S. aureus nasal carriage. Within 984 carriage isolates, we also compared healthcare-associated (HA) carriage with community-associated (CA) carriage. All major global lineages were represented in both bacteraemia and carriage, with no evidence for different infection rates. However, kmers tagging trimethoprim resistance-conferring mutation F99Y in dfrB were significantly associated with bacteraemia-vs-carriage (P=10-8.9-10-9.3). Pooling variation within genes, bacteraemia-vs-carriage was associated with the presence of mecA (HMP=10-5.3) as well as the presence of SCCmec (HMP=10-4.4). Among S. aureus carriers, no lineages were associated with HA-vs-CA carriage. However, we found a novel signal of HA-vs-CA carriage in the foldase protein prsA, where kmers representing conserved sequence allele were associated with CA carriage (P=10-7.1-10-19.4), while in gyrA, a ciprofloxacin resistance-conferring mutation, L84S, was associated with HA carriage (P=10-7.2). In an extensive study of S. aureus bacteraemia and nasal carriage in the UK, we found strong evidence that all S. aureus lineages are equally capable of causing bloodstream infection, and of being carried in the healthcare environment. Genomic variation in the foldase protein prsA is a novel genomic marker of healthcare origin in S. aureus but was not associated with bacteraemia. AMR determinants were associated with both bacteraemia and healthcare-associated carriage, suggesting that AMR increases the propensity not only to survive in healthcare environments, but also to cause invasive disease.

Anti-spike antibody response to natural SARS-CoV-2 infection in the general population

Understanding the trajectory, duration, and determinants of antibody responses after SARS-CoV-2 infection can inform subsequent protection and risk of reinfection, however large-scale representative studies are limited. Here we estimated antibody response after SARS-CoV-2 infection in the general population using representative data from 7,256 United Kingdom COVID-19 infection survey participants who had positive swab SARS-CoV-2 PCR tests from 26-April-2020 to 14-June-2021. A latent class model classified 24% of participants as 'non-responders' not developing anti-spike antibodies, who were older, had higher SARS-CoV-2 cycle threshold values during infection (i.e. lower viral burden), and less frequently reported any symptoms. Among those who seroconverted, using Bayesian linear mixed models, the estimated anti-spike IgG peak level was 7.3-fold higher than the level previously associated with 50% protection against reinfection, with higher peak levels in older participants and those of non-white ethnicity. The estimated anti-spike IgG half-life was 184 days, being longer in females and those of white ethnicity. We estimated antibody levels associated with protection against reinfection likely last 1.5-2 years on average, with levels associated with protection from severe infection present for several years. These estimates could inform planning for vaccination booster strategies.

Changes in paediatric respiratory infections at a UK teaching hospital 2016-2021; impact of the SARS-CoV-2 pandemic

Objective To describe the impact of the SARS-CoV-2 pandemic on the incidence of paediatric viral respiratory tract infection in Oxfordshire, UK. Methods Data on paediatric Emergency Department (ED) attendances (0-15 years inclusive), respiratory virus testing, vital signs and mortality at Oxford University Hospitals were summarised using descriptive statistics. Results Between 1-March-2016 and 30-July-2021, 155,056 ED attendances occurred and 7,195 respiratory virus PCRs were performed. Detection of all pathogens was suppressed during the first national lockdown. Rhinovirus and adenovirus rates increased when schools reopened September-December 2020, then fell, before rising in March-May 2021. The usual winter RSV peak did not occur in 2020/21, with an inter-seasonal rise (32/1,000 attendances in 0-3 yr olds) in July 2021. Influenza remained suppressed throughout. A higher paediatric early warning score (PEWS) was seen for attendees with adenovirus during the pandemic compared to pre-pandemic (p = 0.04, Mann-Witney U test), no other differences in PEWS were seen. Conclusions SARS-CoV-2 caused major changes in the incidence of paediatric respiratory viral infection in Oxfordshire, with implications for clinical service demand, testing strategies, timing of palivizumab RSV prophylaxis, and highlighting the need to understand which public health interventions are most effective for preventing respiratory virus infections.

A genomic epidemiological study shows that prevalence of antimicrobial resistance in Enterobacterales is associated with the livestock host, as well as antimicrobial usage

Enterobacterales from livestock are potentially important reservoirs for antimicrobial resistance (AMR) to pass through the food chain to humans, thereby increasing the AMR burden and affecting our ability to tackle infections. In this study 168 isolates from four genera of the order Enterobacterales, primarily Escherichia coli, were purified from livestock (cattle, pigs and sheep) faeces from 14 farms in the United Kingdom. Their genomes were resolved using long- and short-read sequencing to analyse AMR genes and their genetic context, as well as to explore the relationship between AMR burden and on-farm antimicrobial usage (AMU), in the three months prior to sampling. Although E. coli isolates were genomically diverse, phylogenetic analysis using a core-genome SNP tree indicated pig isolates to generally be distinct from sheep isolates, with cattle isolates being intermediates. Approximately 28 % of isolates harboured AMR genes, with the greatest proportion detected in pigs, followed by cattle then sheep; pig isolates also harboured the highest number of AMR genes per isolate. Although 90 % of sequenced isolates harboured diverse plasmids, only 11 % of plasmids (n=58 out of 522) identified contained AMR genes, with 91 % of AMR plasmids being from pig, 9 % from cattle and none from sheep isolates; these results indicated that pigs were a principle reservoir of AMR genes harboured by plasmids and likely to be involved in their horizontal transfer. Significant associations were observed between AMU (mg kg⁻¹) and AMR. As both the total and the numbers of different antimicrobial classes used on-farm increased, the risk of multi-drug resistance (MDR) in isolates rose. However, even when AMU on pig farms was comparatively low, pig isolates had increased likelihood of being MDR; harbouring relatively more resistances than those from other livestock species. Therefore, our results indicate that AMR prevalence in livestock is not only influenced by recent AMU on-farm but also livestock-related factors, which can influence the AMR burden in these reservoirs and its plasmid mediated transmission.

Epidemiology of Mycobacterium abscessus in England: an observational study

BACKGROUND

Mycobacterium abscessus has emerged as a significant clinical concern following reports that it is readily transmissible in health-care settings between patients with cystic fibrosis. We linked routinely collected whole-genome sequencing and health-care usage data with the aim of investigating the extent to which such transmission explains acquisition in patients with and without cystic fibrosis in England.

METHODS

In this retrospective observational study, we analysed consecutive M abscessus whole-genome sequencing data from England (beginning of February, 2015, to Nov 14, 2019) to identify genomically similar isolates. Linkage to a national health-care usage database was used to investigate possible contacts between patients. Multivariable regression analysis was done to investigate factors associated with acquisition of a genomically clustered strain (genomic distance <25 single nucleotide polymorphisms [SNPs]).

FINDINGS

2297 isolates from 906 patients underwent whole-genome sequencing as part of the routine Public Health England diagnostic service. Of 14 genomic clusters containing isolates from ten or more patients, all but one contained patients with cystic fibrosis and patients without cystic fibrosis. Patients with cystic fibrosis were equally likely to have clustered isolates (258 [60%] of 431 patients) as those without cystic fibrosis (322 [63%] of 513 patients; p=0·38). High-density phylogenetic clusters were randomly distributed over a wide geographical area. Most isolates with a closest genetic neighbour consistent with potential transmission had no identifiable relevant epidemiological contacts. Having a clustered isolate was independently associated with increasing age (adjusted odds ratio 1·14 per 10 years, 95% CI 1·04-1·26), but not time spent as an hospital inpatient or outpatient. We identified two sibling pairs with cystic fibrosis with genetically highly divergent isolates and one pair with closely related isolates, and 25 uninfected presumed household contacts with cystic fibrosis.

INTERPRETATION

Previously identified widely disseminated dominant clones of M abscessus are not restricted to patients with cystic fibrosis and occur in other chronic respiratory diseases. Although our analysis showed a small number of cases where person-to-person transmission could not be excluded, it did not support this being a major mechanism for M abscessus dissemination at a national level in England. Overall, these data should reassure patients and clinicians that the risk of acquisition from other patients in health-care settings is relatively low and motivate future research efforts to focus on identifying routes of acquisition outside of the cystic fibrosis health-care-associated niche.

FUNDING

The National Institute for Health Research, Health Data Research UK, The Wellcome Trust, The Medical Research Council, and Public Health England.

Quantitative SARS-CoV-2 anti-spike responses to Pfizer-BioNTech and Oxford-AstraZeneca vaccines by previous infection status

OBJECTIVES

We investigated determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) anti-spike IgG responses in healthcare workers (HCWs) following one or two doses of Pfizer-BioNTech or Oxford-AstraZeneca vaccines.

METHODS

HCWs participating in regular SARS-CoV-2 PCR and antibody testing were invited for serological testing prior to first and second vaccination, and 4 weeks post-vaccination if receiving a 12-week dosing interval. Quantitative post-vaccination anti-spike antibody responses were measured using the Abbott SARS-CoV-2 IgG II Quant assay (detection threshold: ≥50 AU/mL). We used multivariable logistic regression to identify predictors of seropositivity and generalized additive models to track antibody responses over time.

RESULTS

3570/3610 HCWs (98.9%) were seropositive >14 days post first vaccination and prior to second vaccination: 2706/2720 (99.5%) were seropositive after the Pfizer-BioNTech and 864/890 (97.1%) following the Oxford-AstraZeneca vaccines. Previously infected and younger HCWs were more likely to test seropositive post first vaccination, with no evidence of differences by sex or ethnicity. All 470 HCWs tested >14 days after the second vaccination were seropositive. Quantitative antibody responses were higher after previous infection: median (IQR) >21 days post first Pfizer-BioNTech 14 604 (7644-22 291) AU/mL versus 1028 (564-1985) AU/mL without prior infection (p < 0.001). Oxford-AstraZeneca vaccine recipients had lower readings post first dose than Pfizer-BioNTech recipients, with and without previous infection, 10 095 (5354-17 096) and 435 (203-962) AU/mL respectively (both p < 0.001 versus Pfizer-BioNTech). Antibody responses >21 days post second Pfizer vaccination in those not previously infected, 10 058 (6408-15 582) AU/mL, were similar to those after prior infection followed by one vaccine dose.

CONCLUSIONS

SARS-CoV-2 vaccination leads to detectable anti-spike antibodies in nearly all adult HCWs. Whether differences in response impact vaccine efficacy needs further study.

Epidemiological data and genome sequencing reveals that nosocomial transmission of SARS-CoV-2 is underestimated and mostly mediated by a small number of highly infectious individuals

OBJECTIVES

Despite robust efforts, patients and staff acquire SARS-CoV-2 infection in hospitals. We investigated whether whole-genome sequencing enhanced the epidemiological investigation of healthcare-associated SARS-CoV-2 acquisition.

METHODS

From 17-November-2020 to 5-January-2021, 803 inpatients and 329 staff were diagnosed with SARS-CoV-2 infection at four Oxfordshire hospitals. We classified cases using epidemiological definitions, looked for a potential source for each nosocomial infection, and evaluated genomic evidence supporting transmission.

RESULTS

Using national epidemiological definitions, 109/803(14%) inpatient infections were classified as definite/probable nosocomial, 615(77%) as community-acquired and 79(10%) as indeterminate. There was strong epidemiological evidence to support definite/probable cases as nosocomial. Many indeterminate cases were likely infected in hospital: 53/79(67%) had a prior-negative PCR and 75(95%) contact with a potential source. 89/615(11% of all 803 patients) with apparent community-onset had a recent hospital exposure. Within 764 samples sequenced 607 genomic clusters were identified (>1 SNP distinct). Only 43/607(7%) clusters contained evidence of onward transmission (subsequent cases within ≤ 1 SNP). 20/21 epidemiologically-identified outbreaks contained multiple genomic introductions. Most (80%) nosocomial acquisition occurred in rapid super-spreading events in settings with a mix of COVID-19 and non-COVID-19 patients.

CONCLUSIONS

Current surveillance definitions underestimate nosocomial acquisition. Most nosocomial transmission occurs from a relatively limited number of highly infectious individuals.

Whole-genome sequencing for surveillance of tuberculosis drug resistance and determination of resistance level in China

OBJECTIVES

Phenotypic drug susceptibility testing for prediction of tuberculosis (TB) drug resistance is slow and unreliable, limiting individualized therapy and monitoring of national TB data. Our study evaluated whole-genome sequencing (WGS) for its predictive accuracy, use in TB drug-resistance surveillance and ability to quantify the effects of resistance-associated mutations on MICs of anti-TB drugs.

METHODS

We used WGS to measure the susceptibility of 4880 isolates to ten anti-TB drugs; for pyrazinamide, we used BACTEC MGIT 960. We determined the accuracy of WGS by comparing the prevalence of drug resistance, measured by WGS, with the true prevalence, determined by phenotypic susceptibility testing. We used the Student-Newman-Keuls test to confirm MIC differences of mutations.

RESULTS

Resistance to isoniazid, rifampin and ethambutol was highly accurately predicted with at least 92.92% (95% confidence interval [CI], 88.19-97.65) sensitivity, resistance to pyrazinamide with 50.52% (95% CI, 40.57-60.47) sensitivity, and resistance to six second-line drugs with 85.05% (95% CI, 80.27-89.83) to 96.01% (95% CI, 93.89-98.13) sensitivity. The rpoB S450L, katG S315T and gyrA D94G mutations always confer high-level resistance, while rpoB L430P, rpoB L452P, fabG1 C-15T and embB G406S often confer low-level resistance or sub-epidemiological cutoff (ECOFF) MIC elevation.

CONCLUSION

WGS can predict phenotypic susceptibility with high accuracy and could be a valuable tool for drug-resistance surveillance and allow the detection of drug-resistance level; It can be an important approach in TB drug-resistance surveillance and for determining therapeutic schemes.

Ten-year longitudinal molecular epidemiology study of Escherichia coli and Klebsiella species bloodstream infections in Oxfordshire, UK

Background

The incidence of Gram-negative bloodstream infections (BSIs), predominantly caused by Escherichia coli and Klebsiella species, continues to increase; however, the causes of this are unclear and effective interventions are therefore hard to design.

Methods

In this study, we sequenced 3468 unselected isolates over a decade in Oxfordshire (UK) and linked this data to routinely collected electronic healthcare records and mandatory surveillance reports. We annotated genomes for clinically relevant genes, contrasting the distribution of these within and between species, and compared incidence trends over time using stacked negative binomial regression.

Results

We demonstrate that the observed increases in E. coli incidence were not driven by the success of one or more sequence types (STs); instead, four STs continue to dominate a stable population structure, with no evidence of adaptation to hospital/community settings. Conversely in Klebsiella spp., most infections are caused by sporadic STs with the exception of a local drug-resistant outbreak strain (ST490). Virulence elements are highly structured by ST in E. coli but not Klebsiella spp. where they occur in a diverse spectrum of STs and equally across healthcare and community settings. Most clinically hypervirulent (i.e. community-onset) Klebsiella BSIs have no known acquired virulence loci. Finally, we demonstrate a diverse but largely genus-restricted mobilome with close associations between antimicrobial resistance (AMR) genes and insertion sequences but not typically specific plasmid replicon types, consistent with the dissemination of AMR genes being highly contingent on smaller mobile genetic elements (MGEs).

Conclusions

Our large genomic study highlights distinct differences in the molecular epidemiology of E. coli and Klebsiella BSIs and suggests that no single specific pathogen genetic factors (e.g. AMR/virulence genes/sequence type) are likely contributing to the increasing incidence of BSI overall, that association with AMR genes in E. coli is a contributor to the increasing number of E. coli BSIs, and that more attention should be given to AMR gene associations with non-plasmid MGEs to try and understand horizontal gene transfer networks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00947-2.

Antibody responses to SARS-CoV-2 vaccines in 45,965 adults from the general population of the United Kingdom

We report that in a cohort of 45,965 adults, who were receiving either the ChAdOx1 or the BNT162b2 SARS-CoV-2 vaccines, in those who had no prior infection with SARS-CoV-2, seroconversion rates and quantitative antibody levels after a single dose were lower in older individuals, especially in those aged >60 years. Two vaccine doses achieved high responses across all ages. Antibody levels increased more slowly and to lower levels with a single dose of ChAdOx1 compared with a single dose of BNT162b2, but waned following a single dose of BNT162b2 in older individuals. In descriptive latent class models, we identified four responder subgroups, including a 'low responder' group that more commonly consisted of people aged >75 years, males and individuals with long-term health conditions. Given our findings, we propose that available vaccines should be prioritized for those not previously infected and that second doses should be prioritized for individuals aged >60 years. Further data are needed to better understand the extent to which quantitative antibody responses are associated with vaccine-mediated protection.

Human and Porcine Transmission of Clostridioides difficile Ribotype 078, Europe

Genomic analysis of a diverse collection of Clostridioides difficile ribotype 078 isolates from Ireland and 9 countries in Europe provided evidence for complex regional and international patterns of dissemination that are not restricted to humans. These isolates are associated with C. difficile colonization and clinical illness in humans and pigs.

The Duration, Dynamics, and Determinants of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Antibody Responses in Individual Healthcare Workers

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) immunoglobulin G (IgG) antibody measurements can be used to estimate the proportion of a population exposed or infected and may be informative about the risk of future infection. Previous estimates of the duration of antibody responses vary.

METHODS

We present 6 months of data from a longitudinal seroprevalence study of 3276 UK healthcare workers (HCWs). Serial measurements of SARS-CoV-2 anti-nucleocapsid and anti-spike IgG were obtained. Interval censored survival analysis was used to investigate the duration of detectable responses. Additionally, Bayesian mixed linear models were used to investigate anti-nucleocapsid waning.

RESULTS

Anti-spike IgG levels remained stably detected after a positive result, for example, in 94% (95% credibility interval [CrI] 91-96%) of HCWs at 180 days. Anti-nucleocapsid IgG levels rose to a peak at 24 (95% CrI 19-31) days post first polymerase chain reaction (PCR)-positive test, before beginning to fall. Considering 452 anti-nucleocapsid seropositive HCWs over a median of 121 days from their maximum positive IgG titer, the mean estimated antibody half-life was 85 (95% CrI 81-90) days. Higher maximum observed anti-nucleocapsid titers were associated with longer estimated antibody half-lives. Increasing age, Asian ethnicity, and prior self-reported symptoms were independently associated with higher maximum anti-nucleocapsid levels and increasing age and a positive PCR test undertaken for symptoms with longer anti-nucleocapsid half-lives.

CONCLUSIONS

SARS-CoV-2 anti-nucleocapsid antibodies wane within months and fall faster in younger adults and those without symptoms. However, anti-spike IgG remains stably detected. Ongoing longitudinal studies are required to track the long-term duration of antibody levels and their association with immunity to SARS-CoV-2 reinfection.

The Duration, Dynamics, and Determinants of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Antibody Responses in Individual Healthcare Workers

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) immunoglobulin G (IgG) antibody measurements can be used to estimate the proportion of a population exposed or infected and may be informative about the risk of future infection. Previous estimates of the duration of antibody responses vary.

METHODS

We present 6 months of data from a longitudinal seroprevalence study of 3276 UK healthcare workers (HCWs). Serial measurements of SARS-CoV-2 anti-nucleocapsid and anti-spike IgG were obtained. Interval censored survival analysis was used to investigate the duration of detectable responses. Additionally, Bayesian mixed linear models were used to investigate anti-nucleocapsid waning.

RESULTS

Anti-spike IgG levels remained stably detected after a positive result, for example, in 94% (95% credibility interval [CrI] 91-96%) of HCWs at 180 days. Anti-nucleocapsid IgG levels rose to a peak at 24 (95% CrI 19-31) days post first polymerase chain reaction (PCR)-positive test, before beginning to fall. Considering 452 anti-nucleocapsid seropositive HCWs over a median of 121 days from their maximum positive IgG titer, the mean estimated antibody half-life was 85 (95% CrI 81-90) days. Higher maximum observed anti-nucleocapsid titers were associated with longer estimated antibody half-lives. Increasing age, Asian ethnicity, and prior self-reported symptoms were independently associated with higher maximum anti-nucleocapsid levels and increasing age and a positive PCR test undertaken for symptoms with longer anti-nucleocapsid half-lives.

CONCLUSIONS

SARS-CoV-2 anti-nucleocapsid antibodies wane within months and fall faster in younger adults and those without symptoms. However, anti-spike IgG remains stably detected. Ongoing longitudinal studies are required to track the long-term duration of antibody levels and their association with immunity to SARS-CoV-2 reinfection.

Genomic network analysis of environmental and livestock F-type plasmid populations

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum β-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR may reflect their suitability for rapid niche adaptation.

Preferences for Medical Consultations from Online Providers: Evidence from a Discrete Choice Experiment in the United Kingdom

BACKGROUND

In the UK, consultations for prescription medicines are available via private providers such as online pharmacies. However, these providers may have lower thresholds for prescribing certain drugs. This is a particular concern for antibiotics, given the increasing burden of antimicrobial resistance. Public preferences for consultations with online providers are unknown, hence the impact of increased availability of online consultations on antibiotic use and population health is unclear.

OBJECTIVE

To conduct a discrete choice experiment survey to understand UK public preferences for seeking online consultations, and the factors that influence these preferences, in the context of having symptoms for which antibiotics may be appropriate.

METHODS

In a survey conducted between July and August 2018, general population respondents completed 16 questions in which they chose a primary care consultation via either their local medical centre or an online provider. Consultations were described in terms of five attributes, including cost and similarity to traditional 'face-to-face' appointments. Choices were modelled using regression analysis.

RESULTS

Respondents (n = 734) placed a high value on having a consultation via their local medical centre rather than an online provider, and a low value on consultations by phone or video. However, respondents characterised as 'busy young professionals' showed a lower strength of preference for traditional consultations, with a higher concern for convenience.

CONCLUSION

Before COVID-19, the UK public had limited appetite for consultations with online providers, or for consultations that were not face-to-face. Nevertheless, prescriptions from online providers should be monitored going forward, particularly for antibiotics, and in key patient groups.

Population-level faecal metagenomic profiling as a tool to predict antimicrobial resistance in Enterobacterales isolates causing invasive infections: An exploratory study across Cambodia, Kenya, and the UK

BACKGROUND

Antimicrobial resistance (AMR) in Enterobacterales is a global health threat. Capacity for individual-level surveillance remains limited in many countries, whilst population-level surveillance approaches could inform empiric antibiotic treatment guidelines.

METHODS

In this exploratory study, a novel approach to population-level prediction of AMR in Enterobacterales clinical isolates using metagenomic (Illumina) profiling of pooled DNA extracts from human faecal samples was developed and tested. Taxonomic and AMR gene profiles were used to derive taxonomy-adjusted population-level AMR metrics. Bayesian modelling, and model comparison based on cross-validation, were used to evaluate the capacity of each metric to predict the number of resistant Enterobacterales invasive infections at a population-level, using available bloodstream/cerebrospinal fluid infection data.

FINDINGS

Population metagenomes comprised samples from 177, 157, and 156 individuals in Kenya, the UK, and Cambodia, respectively, collected between September 2014 and April 2016. Clinical data from independent populations included 910, 3356 and 197 bacterial isolates from blood/cerebrospinal fluid infections in Kenya, the UK and Cambodia, respectively (samples collected between January 2010 and May 2017). Enterobacterales were common colonisers and pathogens, and faecal taxonomic/AMR gene distributions and proportions of antimicrobial-resistant Enterobacterales infections differed by setting. A model including terms reflecting the metagenomic abundance of the commonest clinical Enterobacterales species, and of AMR genes known to either increase the minimum inhibitory concentration (MIC) or confer clinically-relevant resistance, had a higher predictive performance in determining population-level resistance in clinical Enterobacterales isolates compared to models considering only AMR gene information, only taxonomic information, or an intercept-only baseline model (difference in expected log predictive density compared to best model, estimated using leave-one-out cross-validation: intercept-only model = -223 [95% credible interval (CI): -330,-116]; model considering only AMR gene information = -186 [95% CI: -281,-91]; model considering only taxonomic information = -151 [95% CI: -232,-69]).

INTERPRETATION

Whilst our findings are exploratory and require validation, intermittent metagenomics of pooled samples could represent an effective approach for AMR surveillance and to predict population-level AMR in clinical isolates, complementary to ongoing development of laboratory infrastructures processing individual samples.

Niche and local geography shape the pangenome of wastewater- and livestock-associated Enterobacteriaceae

Escherichia coli and other Enterobacteriaceae are diverse species with "open" pangenomes, where genes move intra- and interspecies via horizontal gene transfer. However, most analyses focus on clinical isolates. The pangenome dynamics of natural populations remain understudied, despite their suggested role as reservoirs for antimicrobial resistance (AMR) genes. Here, we analyze near-complete genomes for 827 Enterobacteriaceae (553 Escherichia and 274 non-Escherichia spp.) with 2292 circularized plasmids in total, collected from 19 locations (livestock farms and wastewater treatment works in the United Kingdom) within a 30-km radius at three time points over a year. We find different dynamics for chromosomal and plasmid-borne genes. Plasmids have a higher burden of AMR genes and insertion sequences, and AMR-gene-carrying plasmids show evidence of being under stronger selective pressure. Environmental niche and local geography both play a role in shaping plasmid dynamics. Our results highlight the importance of local strategies for controlling the spread of AMR.

Occurrence and characterization of Escherichia coli ST410 co-harbouring blaNDM-5, blaCMY-42 and blaTEM-190 in a dog from the UK

BACKGROUND/OBJECTIVES

Carbapenemase-producing Enterobacteriaceae (CPE) are a public health threat, and have been found in humans, animals and the environment. Carbapenems are not authorized for use in EU or UK companion animals, and the prevalence of carbapenem-resistant Gram-negative bacilli (CRGNB) in this population is unknown.

METHODS

We investigated CRGNB isolated from animal specimens received by one diagnostic laboratory from 34 UK veterinary practices (September 2015-December 2016). Any Gram-negative isolates from clinical specimens showing reduced susceptibility to fluoroquinolones and/or aminoglycosides and/or cephalosporins were investigated phenotypically and genotypically for carbapenemases. A complete genome assembly (Illumina/Nanopore) was generated for the single isolate identified to investigate the genetic context for carbapenem resistance.

RESULTS

One ST410 Escherichia coli isolate [(CARB35); 1/191, 0.5%], cultured from a wound in a springer spaniel, harboured a known carbapenem resistance gene (blaNDM-5). The gene was located in the chromosome on an integrated 100 kb IncF plasmid, also harbouring other drug resistance genes (mrx, sul1, ant1 and dfrA). The isolate also contained blaCMY-42 and blaTEM-190 on two separate plasmids (IncI1 and IncFII, respectively) that showed homology with other publicly available plasmid sequences from Italy and Myanmar.

CONCLUSIONS

Even though the use of carbapenems in companion animals is restricted, the concurrent presence of blaCMY-42 and other antimicrobial resistance genes could lead to co-selection of carbapenemase genes in this population. Further studies investigating the selection and flow of plasmids carrying important resistance genes amongst humans and companion animals are needed.