Post-acute COVID-19 associated with evidence of bystander T-cell activation and a recurring antibiotic-resistant bacterial pneumonia
Abstract
Here, we describe the case of a COVID-19 patient who developed recurring ventilator-associated pneumonia caused by Pseudomonas aeruginosa that acquired increasing levels of antimicrobial resistance (AMR) in response to treatment. Metagenomic analysis revealed the AMR genotype, while immunological analysis revealed massive and escalating levels of T-cell activation. These were both SARS-CoV-2 and P. aeruginosa specific, and bystander activated, which may have contributed to this patient’s persistent symptoms and radiological changes.
Introduction
The COVID-19 pandemic has brought with it the largest ever cohort of patients requiring mechanical ventilation. Complications associated with such severe viral infections are many-fold, and include increased susceptibility to secondary bacterial infections (Zhou et al., 2020; Langford et al., 2020), as well as post-acute COVID-19, where patients experience symptoms extending beyond 3 weeks from the onset of their first COVID-19 symptoms (Greenhalgh et al., 2020). The first report of secondary infections in COVID-19 patients was from Wuhan in March 2020, where 15% of hospitalized patients developed secondary infections, and of those who did not survive their SARS-CoV-2 infection, 50% had a secondary bacterial infection (Zhou et al., 2020). Since then many COVID-19 studies reporting secondary infections have been published, with a recent meta-analysis of 24 independent studies that included 3338 patients from five countries reporting that 14.3% of hospitalized COVID-19 patients developed a secondary bacterial infection, which is associated with significant morbidity, mortality and the financial costs associated with prolonged hospitilisation (Langford et al., 2020). The incidence of post-acute COVID-19 varies depending on the group of patients considered, with approximately 10% of patients who have tested positive for SARS-CoV-2 virus remaining unwell beyond 3 weeks (Greenhalgh et al., 2020). However, this can be as high as 74% when hospitalised patients are considered, where symptoms include breathlessness and excessive fatigue, with abnormal radiological features reported in 12% of this cohort (Arnold et al., 2020).
The DISCOVER study (DIagnostic and Severity markers of COVID-19 to Enable Rapid triage, REC: 20/YH/0121) was established in March 2020 to collect and analyse longitudinal samples from COVID-19 patients. One study participant, an otherwise healthy male between 45 and 55 years of age presented to hospital with Type-1 respiratory failure (hypoxaemia), 20 days after he tested positive for SARS-CoV-2 by RT-PCR. At the time of testing he was asymptomatic (tested as a house-hold contact of a health-care worker), and this represents Day one on the time-line presented in Figure 1. He became symptomatic for COVID-19 13 days after this, which is within the incubation period for SARS-CoV-2 as defined by the CDC of 14 days (COVID-19 (Coronavirus Disease), 2020), and his health declined over the following week. Upon admission to hospital a chest X-ray was taken, and he was admitted to the ICU where he was mechanically ventilated (Figure 1). An RT-PCR test on an endotracheal sample collected at this time did not detect SARS-CoV-2 suggesting that he had cleared the viral infection.
Once admitted to our intensive care unit he was assessed, as all our patients are daily, for signs of additional infection through a process led by a consultant medical microbiologist along with senior pharmacist. At these meetings, the clinical status of the patient, radiological features, oxygen requirements, inflammatory indices (WCC, CRP, and procalcitonin) and any culture results are discussed. We had no concerns at the time of admission of this patient to the ICU that he had any infection other than that of SARS-CoV-2. Computerised tomography (CT) scans of the patient were collected through his time in the ICU and the classic features of COVID-19, including ground-glass opacities, were evident throughout (Figure 3—figure supplement 1).
After a week in the ICU he was diagnosed with ventilator-associated pneumonia (VAP) based on the clinical, radiological, biochemical, and microbiological parameters described above. An antibiotic susceptible Pseudomonas aeruginosa strain was cultured from an aspirate collected from his endotracheal tube (ETT-1, Figure 1), (minimum inhibitory concentration (MIC) meropenem 0.5 mg/l; piperacillin-tazobactam (pip/taz) 8 mg/l). He was prescribed a seven-day course of pip/taz (4.5 g every 6 hr), and clinically recovered from this bacterial infection. Eight days after finishing this first course of antibiotics his VAP recurred (diagnosed as described above), and a pip/taz resistant P. aeruginosa was cultured from ETT-2 (minimum inhibitory concentration (MIC) >16 mg/l). He was prescribed a seven-day course of meropenem (1 g every 8 hr) and showed signs of clinical improvement. There were no further concerns about a bacterial infection during his ICU stay (i.e. he did not meet any of the clinical diagnostic criteria described above), although a third P. aeruginosa was cultured from a sample (ETT-3) collected two days after he completed his course of meropenem that was resistant to both pip/taz and meropenem (MICs > 16 mg/l and >8 mg/l respectively). His health continued to improve, and he was discharged to the respiratory ward a week later. The patient attended a follow-up clinic 3 weeks after hospital discharge where he reported on-going symptoms such as breathlessness and myalgia, and there was a significant drop in his oxygen saturations after mild exertion to 84%. Eight months following his positive test he is still reporting breathlessness and fatigue, common symptoms of the newly defined ‘Long COVID’ condition.
To understand the cellular and molecular dynamics of this post-acute COVID-19 case from the perspective of both the pathogen and patient’s immune response, longitudinal respiratory and blood samples were collected and analysed with a view to identifying early diagnostic biomarkers of infection onset and potential opportunities for immunotherapeutic intervention.
Results
Metagenomics was used to characterise in depth the composition and genomic features of the bacteria present in this patient’s lower respiratory tract from when he first developed VAP (ETT-1, Figure 1) to when he had recovered from his second VAP (ETT-3, Figure 1). ETT-2 collected during his second bout of VAP was unfortunately not available for sequencing. We extracted the entire genetic material from 500 μl of ETTs 1 and 3 with no bacterial enrichment or human DNA depletion steps, and to ensure data for the bacterial component of the samples was generated, these were sequenced on a PromethION (Oxford Nanopore technology). The bacterial DNA from ETT-1 was entirely that of P. aeruginosa corresponding to multi-locus sequence type ST253 (Jolley et al., 2018), a world-wide clone frequently associated with AMR epidemics (Treepong et al., 2018; Figure 2). There was a greater diversity of bacterial DNA in ETT-3, as the patient recovered from his SARS-CoV-2 and P. aeruginosa infections (Figure 2b), and we were able to construct a whole genome of the P. aeruginosa strain from within this, which again corresponded to multi-locus sequence type ST253. P. aeruginosa is a ubiquitous environmental bacterium, and while it can be found as a commensal in the respiratory tract of some individuals, this is typically only when they have some underlying chronic suppurative lung disease such as bronchiectasis. Prior to his COVID-19 diagnosis this relatively young patient was healthy with no prior medical issues, which suggests the P. aeruginosa infections were hospital or ICU acquired.
Resistance to beta-lactam antibiotics by P. aeruginosa can be multifactorial and includes the acquisition of single nucleotide polymorphisms (SNPs) in efflux and porin genes that affect the passage of the antibiotic into and out of the bacterial cell (Blair et al., 2015). Analysis of the sequence data from ETT-3 generated a whole genome for P. aeruginosa, again corresponding to ST253, that aligned with >99% identify to the genome from ETT-1. There were however two noteworthy SNPs in this later sample that explain the increased AMR of this isolate. The multi-drug efflux pump, MexAB, is transcriptionally regulated by the translated product of the mexR gene, MexR. We found a SNP in mexR that converts a key amino-acid (Arg91-Gln) in the DNA-binding domain of the MexR protein, that would likely result in the de-repression of this efflux system (Saito et al., 2003). A second SNP was found that introduced a premature stop codon (Tyr120-stop) in the gene encoding the outer membrane protein OprD, a protein with a well-established role in the entry of meropenem into the bacterial cell (Quale et al., 2006). The two P. aeruginosa strains could represent a parent strain that acquired resistance upon exposure to antibiotics, that wasn’t fully cleared and so re-emerged; or a new infection or colonisation event with a closely related AMR strain circulating within the ICU. Unfortunately, we are unable to tell which of these scenario is more likely from these data.
Having recovered from both a viral and a recurring bacterial infection, we also sought to analyse the kinetics of the patient’s innate and adaptive immune response across his time in the ICU (at days 23, 28, 38, 45, and 58 from testing positive for SARS-CoV2, Figure 1), as previous work has demonstrated striking immunological features associated with COVID-19 patients (Mathew et al., 2020). Robust activation of a broad range of immune cell subsets was revealed by flow cytometry, when compared to healthy controls. What was particularly striking was the level of activation and proliferation of CD4+, CD8+, and TCR-γδ T-cells, that appeared to wane between days 23–28 but were subsequently boosted after day 28, concomitant with the onset of the secondary bacterial infection (Figure 3a–c and Figure 3e–g, gating strategies in Figure 3—figure supplement 1). Similar levels of activation of T-cells has been reported for subsets of COVID-19 patients in other studies (Mathew et al., 2020). From day 28 a robust and steady increase in the activation and proliferation of conventional and TCR-γδ-T-cells was evident with approximately 20% and 40% of total CD4+ and CD8+ T-cells co-expressing the activation markers HLA-DR and CD38, respectively by day 58. A similar steady increase in activation levels, albeit at lower magnitudes, was observed for Natural Killer (NK) CD56dim and CD56bright cells (Figure 3d and h, gating strategies in Figure 3—figure supplement 1). Similar perturbations in the frequency and activation phenotypes of monocytes, blood monocyte-derived macrophages as well as neutrophils could be detected from day 28 onwards (Figure 3—figure supplement 2a–e). A robust IgG response to both SARS-CoV-2 and P. aeruginosa antigen was also detected in this patient (Figure 3—figure supplement 3).
Given the scale of the cellular response of this patient, we investigated whether the large expansions of CD4+ and CD8+ T-cells were due to a T-cell response targeting SARS-CoV-2, the secondary bacterial infection or to bystander T-cell activation (Sandalova et al., 2010; Rivino et al., 2015). We included the analysis of bystander activation, which is T-cell receptor-independent and cytokine-mediated, as it is known to occur during other types of acute viral infection (Sandalova et al., 2010; Rivino et al., 2015). The role of this type of bystander activated T-cells in the recovery of patients remains largely unclear, as these cells can participate in protective immunity towards the virus but can also contribute to tissue damage (Sandalova et al., 2010; Rivino et al., 2015; Kim and Shin, 2019; Maini et al., 2000). We performed a brief stimulation of the patient’s peripheral blood mononuclear cells (PBMCs), with or without overlapping peptides spanning the sequences of SARS-CoV-2 proteins (spike, membrane (M) and nucleoprotein (N)); of an immunodominant P. aeruginosa antigen (OprF); as well as an immunodominant human Cytomegalovirus protein (HCMV pp65) as an indication of non-T-cell receptor (TCR) driven bystander activation. This was followed by intracellular cytokine staining to detect production of IFN-γ and TNF-α by the specific T cells. Following the encounter with specific peptides we observed a robust CD4+ T-cell response to SARS-CoV-2 which decreased from day 23 to day 58, and a more modest SARS-CoV-2-specific CD8+ T-cell response. This is in line with published work, suggesting higher magnitudes of SARS-CoV-2-specific CD4+ versus CD8+ T-cells in severe COVID-19 patients (Peng et al., 2020), as well as consistent detection of virus-specific CD4+ T cells in recovered patients (Grifoni et al., 2020; Figure 4a–c and Figure 4—figure supplement 1a–b). CD4+ and CD8+ T-cell responses targeting P. aeruginosa peptides could also be detected, and increased over time concomitant with the onset of the recurring bacterial infection (Figure 4b–c).
The number of T-cells responding specifically to SARS-CoV-2 and P. aeruginosa represent only a fraction of those activated, which suggests that many of those detected in the patient’s blood may be bystander activated. In support of this we detected an increased magnitude of CMV-specific T-cell responses over time, particularly in the CD8+ compartment (Figure 3b–c). We therefore considered whether CMV-specific T-cells were activated and proliferating in the blood of the COVID-19 patient. To test this we performed a brief stimulation of PBMCs with or without HCMV pp65 peptide (CMV) or SARS-CoV-2 spike 1 one peptide pools (the predominant SARS-CoV-2 target in this patient, within those tested) or with a positive control, and assessed the expression of the activation and proliferation markers HLA-DR, CD38 and Ki67 by the cytokine-producing cells that were responding to the peptide pools. Peptide stimulations and negative controls were performed in duplicate wells (Figure 4—figure supplement 1a–b). Both CMV and spike-1-specific CD8+ T-cells showed increased co-expression of HLA-DR and CD38 and of Ki67, indicating that these cells were activated and proliferating ex vivo, while this was less evident for their CD4+ T-cell counterparts (Figure 4d–e and g–h). This observation is in line with previous reports of bystander T-cell activation for CMV-specific CD8+ T-cell populations (Sandalova et al., 2010; Rivino et al., 2015). We furthermore analysed the phenotype of CMV and spike 1-specific T-cells on the basis of expression of CCR7 and CD45RA, which identifies naïve (CCR7+ CD45RA+), central memory (TCM: CCR7+ CD45RA-), effector memory (TEM: CCR7- CD45RA-) and effector memory re-expressing RA T-cells (TEMRA: CCR7-CD45RA+). CMV and spike 1-specific CD4+ T cells were mainly contained respectively within the TEM and the TEM/TCM subsets, while CMV and spike 1-specific CD8+ T cells were mainly contained respectively within the TEMRA and the TEM subsets. The predominant TEMRA phenotype of SARS-CoV-2 specific CD8+ T-cells in convalescent COVID patients is in line with recent work (Neidleman et al., 2020). The presence of a large proportion of CD8+ TEMRA within the CMV and spike 1-specific CD8+ T cell populations also supports the notion that these cells are stimulated through cytokine receptors rather than through the TCR, as TEMRA cells are believed to arise following antigen withdrawal and cytokine-mediated stimulation (for e.g. IL-15), as TCR ligation downregulates CD45RA expression both in vitro and during responses in vivo (Henson et al., 2012).
In summary, the immune profiles show robust innate and adaptive immune activation in a COVID-19 patient experiencing a secondary bacterial infection, and a detectable T-cell and IgG response targeting both SARS-CoV-2 and P. aeruginosa. Furthermore, we show that CD8+ T-cells targeting viral antigens from SARS-CoV-2 and from the persistent virus HCMV display expression of markers of activation and proliferation and are predominantly TEMRA and TEM cells. Although it is a challenging task to demonstrate that these cells have not encountered their cognate viral antigens in vivo, the TEMRA phenotype of these cells makes it unlikely and suggests that these cells may be proliferating in a bystander way through the action of cytokines present in the inflammatory milieu.
Discussion
In this case report we demonstrate the potential of the application of multidisciplinary technologies to longitudinally-collected patient samples to define the complex dynamics of patient-pathogen-therapy interactions in real-time. Metagenomics directly applied to respiratory samples facilitated the identification of the SNPs responsible for the AMR phenotype of the later P. aeruginosa isolate. However, the most striking feature of this COVID-19 case was the escalating number of circulating activated T-cells more than two months after testing positive for SARS-CoV-2, and 6 weeks after clearing the viral infection. While we could attribute some of this to the recurring bacterial infection, the scale of the activation considered alongside our evidence of increased frequencies of T-cells specific for unrelated antigens, suggests there may be a significant amount of bystander activation.
Studies on large COVID-19 patient cohorts have shown that the period of the peak of T-cell activation in COVID-19 is prolonged compared to other acute viral infections or after vaccination with live attenuated viruses (Mathew et al., 2020), which may be due to the inability of the immune system to downregulate its response. Here, we speculate that due to SARS-CoV2 infection the patient displays a heightened immune system, which was further stimulated by the recurring P. aeruginosa infections, which led to bystander activation of T cells specific for antigens unrelated to either SARS-CoV2 or P. aeruginosa. Bystander T-cell activation could have played a critical role in the severity of illness and the longer-term complications associated with the development of post-acute COVID-19 experienced by this patient (Grifoni et al., 2020). Given the recent appreciation for the role of corticosteroids in reducing the risk of death following SARS-CoV-2 infection by 20% (Sterne et al., 2020), this case suggests that targeting its use to patients with immunological responses as described here, may also improve their long-term recovery.
Materials and methods
Patient recruitment
Request a detailed protocolThe patient was enrolled onto the DISCOVER study (Diagnostic and Severity markers of COVID-19 to Enable Rapid triage study), a single centre prospective study recruiting consecutive patients admitted with COVID-19, from 30.03.2020 until present (Ethics approval via South Yorkshire REC: 20/YH/0121, CRN approval no: 45469). Blood/serum samples from pre-pandemic healthy controls and asymptomatic healthy controls were obtained under the Bristol Biobank (NHS Research Ethics Committee approval ref 14/WA/1253).
Clinical microbiology
Request a detailed protocolAll ETT samples were processed in the Severn Infection Services laboratory as per standard operating procedures. Antibiotic susceptibility testing was performed by either disk diffusion (according to EUCAST version nine guidance EUCAST, 2019) or on a VITEK two machine (BioMeriux, France).
SARS-CoV-2 test
Request a detailed protocolSARS-CoV-2 test was performed by an in-house RT-PCR at the regional South West Public Health England Regional Virology laboratory, utilising a PHE approved assay at the time of testing.
DNA extractions
Request a detailed protocolDNA from the endotrachaeal aspirates were extracted using the CTAB/Phenol:Chloroform:Isoamyl alcohol and bead-beating approach of Griffiths et al., 2000 with modifications of DeAngelis et al., 2009. Phase lock gel tubes (ThermoFisher) and linear polyacrylamide (Sigma) were included to increase nucleic acid yields and total DNA was resuspended in 50 µl of DNase/RNase free water before storage at −20°C.
DNA sequencing
Request a detailed protocolSamples were prepared for sequencing using SQK-LSK109 kit (Oxford Nanopore) with 1 μg DNA starting input as per manufacture’s protocol. Briefly 1 μg of DNA was end repaired and a tailed using NEBNext Ultra II module E7546 (3.5 μl End Repair Buffer, 2 ul FFPE repair mix, 3.5 μl Ultra II end-prep reaction buffer and 3 μl of Ultra II end-prep enzyme mix to 1 μg DNA in a total of reaction volume of 60 μl). This was incubated at 20°C for 5 min followed by 65°C for 5 min. Clean-up was performed using AMPure XP beads (Beckman Coulter) in a 1X ratio. Adaptors were ligated by adding 5 μl adaptor mix (Oxford Nanopore) 25 μl ligation buffer (Oxford Nanopore) and 10 μl Quick T4 ligase (NEB Module MO202). Following a 20 min incubation at room temperature the adaptor ligated DNA was cleaned using AMPure beads in a 0.8 X ratio and washed using Long Fragment Buffer (Oxford Nanopore) before eluting in 25 μl of elution buffer (Oxford Nanopore). Final quantification by fluorometry (Qubit) was performed and 300 ng DNA prepared for sequencing according to the manufacturer’s instructions (Oxford Nanopore). Sequencing was performed on a PromethION R9.4.1 flow cell (FLO-PRO002) and run for 48 hr using live basecalling, files were outputted in Fast5 and Fastq format. Human-filtered sequencing data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession PRJEB40239.
Bioinformatics
Request a detailed protocolBioinformatics was orchestrated using reticulatus (https://github.com/SamStudio8/reticulatus/ Köster and Rahmann, 2012) with configuration: Flye (v2.6) metagenomic assembly (Kolmogorov et al., 2019), Racon polishing (v1.4.9 + GPU, two rounds) (Vaser et al., 2017) and Medaka polishing (v0.8.0 + GPU, one round). Contigs were assigned to a taxon by Kraken 2 (Wood et al., 2019). SNP differences between the P. aeruginosa contigs for ETT-1 and ETT-3 was determined using Mauve (Darling et al., 2004). Contigs were loaded into BRIG (v0.95) (Alikhan et al., 2011) as concentric rings and compared against the PA14 or JCSC1435 reference genomes using blastn (ncbi-blast 2.10.1+) (Alikhan et al., 2011). MLST types were identified by mapping contigs against P. aeruginosa or S. haemolyticus according to the schemes held in the pubmlst databases (Arnold et al., 2020).
PBMC isolation
Request a detailed protocolBlood samples were collected from the COVID-19 patient in EDTA vacutainer tubes and PBMCs isolated from peripheral blood by Ficoll gradient purification and cryopreserved. Healthy donor PBMCs were obtained from Bristol Biobank (REC: 14/WA/1253).
Synthetic peptides
Request a detailed protocol15-mer peptides overlapping by 10 amino acids and spanning the sequences of SARS-CoV-2 spike (Accession Number: NC_045512.2, Protein ID: YP_009724390.1) and HCMV pp65 (AD169 strain) were purchased from Mimotopes (Australia). The purity of the peptides was >80% (Spike) or >70% (pp65) and peptides were dissolved as described previously (Rivino et al., 2015). SARS-CoV-2 M and N Peptivator peptide pools were purchased from Miltenyi. An OprF (PA1777) peptide library comprising 20-mer peptides overlapping by 10 amino acids was synthesized by GL Biochem Ltd., Shanghai, China (Quigley et al., 2015).
Flow cytometry staining and PBMC stimulation
Request a detailed protocolPBMCs were thawed and either stained ex vivo or stimulated in AIMV 2% FCS with or without peptide pools from SARS-CoV2 spike, M, N, HCMV pp65 (all 1 μg/ml), OrpF PA (10 μg/ml) or with PMA/iono (PMA 10 ng/ml, Iono 100 ng/ml, Sigma Aldrich) for 5 hr at 37°C in the presence of brefeldin A (BD, 5 μg/ml). To assess degranulation, CD107a FITC antibody was added to the cells at the beginning of the stimulation. Cells were stained with a viability dye Zombie Aqua (Biolegend) for 10 min at room temperature and with antibodies targeting surface markers (20 min 4°C, diluted in PBS (HyClone) 1% BSA (Sigma Aldrich)). Cells were fixed for 45 min/overnight in eBioscience Foxp3/Transcription factor fixation/permeabilization buffer (Invitrogen) and intracellular staining was performed using eBioscience Foxp3/Transcription factor permeabilization buffer (Invitrogen) for Ki67 or intracellular cytokines (30 min on ice). Cells were acquired on a BD Fortessa X20 and data analysed using FlowJo software v10.7. A complete list of antibodies is included in Supplementary file 1.
Antibody quantification
Request a detailed protocola. Pseudomonas Aeruginosa: The PA antigen OprF was diluted in coating buffer, added to the wells of MaxiSorp flat-bottom 96-Well immunoplates (Thermo Fisher Scientific, USA) and incubated for 2 hr at 37°C. Immunoplates were blocked with 1% BSA/PBS for 1 hr at 37°C. Sera from each sample were diluted and added to the immunoplates in triplicate and incubated at 4°C overnight. Biotin mouse anti-human IgG (BD Pharmingen, USA) was added to the plate and incubated for 1 hr at RT. Streptavidin-HRP (R and D Systems, UK) was added and incubated for 30 mins at RT. ELISAs were developed using TMB substrate. Results were recorded using a Multiskan GO Microplate Spectrophotometer (Thermo Scientific, USA) and absorbance values (OD 450) were converted to ELISA units. OprF antigen was produced using the recombinant vector pSUMO-OprF as described previously (Quigley et al., 2015). b. ELISA methods were based on the protocols for measuring SARS-CoV-2 seroconversion (Amanat et al., 2020; Stadlbauer et al., 2020) with some modifications. A full-length trimeric stabilized version of the SARS-CoV-2 Spike protein was expressed in insect cells using the MultiBac baculovirus expression system and affinity purified, whereas the SARS-CoV-2 N protein was expressed in, and purified from, E. coli. High-binding plates (MaxiSorp, NUNC) were coated with either 10 mg/ml Spike or 20 mg/ml N protein in PBS and left overnight at 4°C followed by blocking for 1 hr (3% BSA/0.1% Tween/PBS). Dilution buffer (1% BSA/0.1% Tween/PBS) containing patient plasma samples collected at different time points, or serum controls, were added to wells either in duplicate (for patient, pre-pandemic (PP) or pooled serum standard (PS)) or in single wells (healthy controls, HC) at a final dilution of 1 in 450 (adjusting for dilution of plasma if obtained via PBMC collection) and left for 2 hr at RT. Goat anti-human IgG-HRP (Southern Biotech) was added at 1 in 25,000 in dilution buffer for 1 hr, followed by development with OPD solution (Sigma) as per the manufacturer’s instructions. Development was quenched after 30 min using 3M HCl and absorbance measured at 492 and 620 nm using a BMG FLUOstar Omega Spectrophotometer. Plots show mean and standard deviation where samples were run in duplicate.
Control samples: Serum from n = 6 adult pre-pandemic donors recruited in 2008 for a vaccine trial were obtained from the Bristol Biobank (NHS REC 14/WA/1253) and were used as known-negative controls on the same plate as the COVID19 case samples, in duplicate. In addition, we have included the OD values from the healthy control donor samples which were used as controls for the T cell assays; these donors were collected during the pandemic but reporting no exposure or symptoms prior to sampling (also obtained via the Biobank as explained in the manuscript methods section). These samples were run in single wells on different plates to the COVID19 case samples; however, all plates included a pooled serum (PS) standard control, in duplicate, to assess for inter-plate variability. The PS average ODs across six plates are therefore shown. The PS, HC and PP samples used in this experiment were heat inactivated at 56°C for 30 min, while the patient plasma were not (after testing negative for SARS-CoV-2). Despite the difference in sample type and handling, we are confident these do not impact on antibody levels or ELISA signal from our own in-house comparisons (data not shown) as well as data from others (Amanat et al., 2020).
Data availability
Human-filtered sequencing data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession PRJEB40239.
-
European Nucleotide ArchiveID PRJEB40239. Metagenomic analysis of respiratory samples from a COVID-19 ICU patient.
References
-
A serological assay to detect SARS-CoV-2 seroconversion in humansNature Medicine 26:1033–1036.https://doi.org/10.1038/s41591-020-0913-5
-
Molecular mechanisms of antibiotic resistanceNature Reviews Microbiology 13:42–51.https://doi.org/10.1038/nrmicro3380
-
Mauve: multiple alignment of conserved genomic sequence with rearrangementsGenome Research 14:1394–1403.https://doi.org/10.1101/gr.2289704
-
Rapid method for coextraction of DNA and RNA from natural environments for analysis of ribosomal DNA- and rRNA-based microbial community compositionApplied and Environmental Microbiology 66:5488–5491.https://doi.org/10.1128/AEM.66.12.5488-5491.2000
-
Properties of end-stage human T cells defined by CD45RA re-expressionCurrent Opinion in Immunology 24:476–481.https://doi.org/10.1016/j.coi.2012.04.001
-
The activation of bystander CD8+ T cells and their roles in viral infectionExperimental & Molecular Medicine 51:1–9.https://doi.org/10.1038/s12276-019-0316-1
-
Assembly of long, error-prone reads using repeat graphsNature Biotechnology 37:540–546.https://doi.org/10.1038/s41587-019-0072-8
-
Snakemake--a scalable bioinformatics workflow engineBioinformatics 28:2520–2522.https://doi.org/10.1093/bioinformatics/bts480
-
Bacterial co-infection and secondary infection in patients with COVID-19: a living rapid review and meta-analysisClinical Microbiology and Infection 26:1622–1629.https://doi.org/10.1016/j.cmi.2020.07.016
-
The role of Virus-Specific Cd8+ cells in liver damage and viral control during persistent hepatitis B virus infectionJournal of Experimental Medicine 191:1269–1280.https://doi.org/10.1084/jem.191.8.1269
-
Interplay of efflux system, ampC, and oprD expression in carbapenem resistance of Pseudomonas aeruginosa clinical isolatesAntimicrobial Agents and Chemotherapy 50:1633–1641.https://doi.org/10.1128/AAC.50.5.1633-1641.2006
-
Chronic infection by Mucoid Pseudomonas aeruginosa Associated with Dysregulation in T-Cell Immunity to Outer Membrane Porin FAmerican Journal of Respiratory and Critical Care Medicine 191:1250–1264.https://doi.org/10.1164/rccm.201411-1995OC
-
Virus-specific T lymphocytes home to the skin during natural dengue infectionScience Translational Medicine 7:278ra35.https://doi.org/10.1126/scitranslmed.aaa0526
-
Global emergence of the widespread Pseudomonas aeruginosa ST235 cloneClinical Microbiology and Infection 24:258–266.https://doi.org/10.1016/j.cmi.2017.06.018
-
Improved metagenomic analysis with Kraken 2Genome Biology 20:257.https://doi.org/10.1186/s13059-019-1891-0
Article and author information
Author details
Funding
Southmead Hospital Charity
- Fergus Hamilton
Wellcome Trust (212258/Z/18/Z)
- Ruth C Massey
Elizabeth Blackwell Institute
- Laura Rivino
UK Research and Innovation (MR/S019553/1)
- Rosemary J Boyton
UK Research and Innovation (MR/R02622X/1)
- Daniel M Altmann
Cystic Fibrosis Trust (CF Trust SRC 015)
- David K Butler
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Kapil Gupta and Imre Berger for kindly providing us the spike protein, Natalie Di Bartolo and Ashley Toye for kindly providing us the N protein, both used for the SARS-CoV-2 serology work. The authors wish to acknowledge the assistance of Dr. Andrew Herman and Helen Rice and the University of Bristol Faculty of Biomedical Sciences Flow Cytometry Facility. We would also like to thank Keith Jolley and the Bristol University UNCOVER team for helpful discussions during the execution of this work and preparation of the manuscript. This work was supported by donations to Southmead Hospital Charity (Registered Charity Number: 1055900), by the Wellcome Trust (reference number: 212258/Z/18/Z) and by the Elizabeth Blackwell Institute, University of Bristol, with funding from the University’s alumni and friends. DKB is supported by a Cystic Fibrosis Trust PhD studentship (CF Trust SRC 015). RJB and DMA are supported by UKRI (MR/S019553/1 and MR/R02622X/1)
Ethics
Human subjects: The patient was enrolled onto the DISCOVER study (Diagnostic and Severity markers of COVID-19 to Enable Rapid triage study), a single centre prospective study recruiting consecutive patients admitted with COVID-19, from 30.03.2020 until present (Ethics approval via South Yorkshire REC: 20/YH/0121, CRN approval no: 45469). Blood/serum samples from pre-pandemic healthy controls and asymptomatic healthy controls were obtained under the Bristol Biobank (NHS Research Ethics Committee approval ref 14/WA/1253).
Copyright
© 2020, Gregorova et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,072
- views
-
- 345
- downloads
-
- 25
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
Variant calling is fundamental in bacterial genomics, underpinning the identification of disease transmission clusters, the construction of phylogenetic trees, and antimicrobial resistance detection. This study presents a comprehensive benchmarking of variant calling accuracy in bacterial genomes using Oxford Nanopore Technologies (ONT) sequencing data. We evaluated three ONT basecalling models and both simplex (single-strand) and duplex (dual-strand) read types across 14 diverse bacterial species. Our findings reveal that deep learning-based variant callers, particularly Clair3 and DeepVariant, significantly outperform traditional methods and even exceed the accuracy of Illumina sequencing, especially when applied to ONT’s super-high accuracy model. ONT’s superior performance is attributed to its ability to overcome Illumina’s errors, which often arise from difficulties in aligning reads in repetitive and variant-dense genomic regions. Moreover, the use of high-performing variant callers with ONT’s super-high accuracy data mitigates ONT’s traditional errors in homopolymers. We also investigated the impact of read depth on variant calling, demonstrating that 10× depth of ONT super-accuracy data can achieve precision and recall comparable to, or better than, full-depth Illumina sequencing. These results underscore the potential of ONT sequencing, combined with advanced variant calling algorithms, to replace traditional short-read sequencing methods in bacterial genomics, particularly in resource-limited settings.
-
- Microbiology and Infectious Disease
The Nix-TB clinical trial evaluated a new 6 month regimen containing three oral drugs; bedaquiline (B), pretomanid (Pa), and linezolid (L) (BPaL regimen) for the treatment of tuberculosis (TB). This regimen achieved remarkable results as almost 90% of the multidrug-resistant or extensively drug-resistant TB participants were cured but many patients also developed severe adverse events (AEs). The AEs were associated with the long-term administration of the protein synthesis inhibitor linezolid. Spectinamide 1599 is also a protein synthesis inhibitor of Mycobacterium tuberculosis with an excellent safety profile, but it lacks oral bioavailability. Here, we propose to replace L in the BPaL regimen with spectinamide (S) administered via inhalation and we demonstrate that inhaled spectinamide 1599, combined with BPa ––BPaS regimen––has similar efficacy to that of the BPaL regimen while simultaneously avoiding the L-associated AEs. The BPaL and BPaS regimens were compared in the BALB/c and C3HeB/FeJ murine chronic TB efficacy models. After 4-weeks of treatment, both regimens promoted equivalent bactericidal effects in both TB murine models. However, treatment with BPaL resulted in significant weight loss and the complete blood count suggested the development of anemia. These effects were not similarly observed in mice treated with BPaS. BPaL and BPa, but not the BPaS treatment, also decreased myeloid to erythroid ratio suggesting the S in the BPaS regimen was able to recover this effect. Moreover, the BPaL also increased concentration of proinflammatory cytokines in bone marrow compared to mice receiving BPaS regimen. These combined data suggest that inhaled spectinamide 1599 combined with BPa is an effective TB regimen without L-associated AEs.