ABSTRACT
Rationale In the upper respiratory tract replicating (culturable) SARS-CoV-2 is recoverable for ∼ 4 to 8 days after symptom onset, however, there is paucity of data about the frequency or duration of replicating virus in the lower respiratory tract (i.e. the human lung).
Objectives We undertook lung tissue sampling (needle biopsy), shortly after death, in 42 mechanically ventilated decedents during the Beta and Delta waves. An independent group of 18 ambulatory patents served as a control group.
Methods Lung biopsy cores from decedents underwent viral culture, histopathological analysis, electron microscopy, transcriptomic profiling and immunohistochemistry.
Results 38% (16/42) of mechanically ventilated decedents had culturable virus in the lung for a median of 15 days (persisting for up to 4 weeks) after symptom onset. Lung viral culture positivity was not associated with comorbidities or steroid use. Delta but not Beta variant lung culture positivity was associated with accelerated death and secondary bacterial infection (p<0.05). Nasopharyngeal culture was negative in 23.1% (6/26) of decedents despite lung culture positivity. This, hitherto, undescribed bio-phenotype of lung-specific persisting viral replication was associated with an enhanced transcriptomic pulmonary pro-inflammatory response but with concurrent viral culture positivity.
Conclusions Concurrent, rather than sequential active viral replication continues to drive a heightened pro-inflammatory response in the human lung beyond the second week of illness and was associated with variant-specific increased mortality and morbidity. These findings have potential implications for the design of interventional strategies and clinical management of patients with severe COVID-19 disease.
Scientific Knowledge on the Subject Investigations to understand SARS-CoV-2 viral shedding (determined by PCR or antigen testing) have extensively focused on samples from the upper respiratory tract. The widely accepted view is that acute severe SARS-CoV-2 infection is characterised by a viral replicative phase in the first week of symptomatic illness followed by a pro-inflammatory immunopathologic phase peaking in the second and third weeks of illness. However, it remains unclear whether detection of SARS-CoV-2 beyond 2 weeks after symptom onset in published studies represent active replication competent virus because it may represent residual genomic or antigenic material in the tissue.
What This Study Adds to the Field We have identified a, hitherto, undescribed bio-phenotype of acute severe COVID-19 characterised by persisting viral replication in the lung for up to 4 weeks after symptom onset. ∼40% of acute severe COVID-19 intensive care unit (ICU) decedents (n=42) had nasopharyngeal swab culture positivity at ∼2 weeks post-symptom onset versus only ∼5% in a group of ambulatory control patients (n=18). There was compartment-specific (nasopharynx versus lung) discordance. The phenotype of lung-specific persisting viral replication was associated with variant-specific accelerated death, an exaggerated inflammatory response, and attenuated T-cell immunity in the lung (based on histopathological and transcriptomic studies). This challenges the traditional view that viral replication occurs during the first 5 to 10 days of illness, which is followed by an effector or hyperinflammatory phase. This is the first study, to our knowledge, to systematically culture virus from the human lung and map out its related clinical determinants, and which describes the human lung transcriptomic profile of culture-positive versus culture-negative patients with severe COVID-19 disease.
Introduction
Coronavirus disease-19 (COVID-19) caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has been the foremost killer globally over the last 3 years. Case fatality risk in hospitalised patients, and particularly in mechanically ventilated patients, during the Beta and Delta waves was particularly high [∼50%-70%; (1)]. Even with the Omicron-related variants, case fatality risk remains significant in the elderly and immunocompromised persons, and in several countries including the UK, Italy, France, Brazil, and prominently in China where there is now an ongoing epidemic of severe COVID-19 disease (2–10). Better therapeutic interventions are needed. However, despite considerable research, the pathogenesis of severe COVID-19, relative to viral kinetics, remains incompletely understood.
SARS-CoV-2 detection (ascertained through PCR positivity or antigen detection) can persist for several weeks from symptom onset (11). Post-mortem studies have shown persistence of SARS-CoV-2 in tissues detected by PCR and immunohistochemistry for up to several weeks after symptom onset (12, 13). However, detection of SARS-CoV-2 in these studies may not represent replication competent virus (detectable only by viral culture) but residual genomic or antigenic material in the tissues. Shedding of replicating virus confirmed through serial viral culture (i.e. in vitro replication in human cell lines) from the upper respiratory tract (URT) has been shown to persist for only ∼2 to 8 days after symptom onset (11, 14–23). These findings have been confirmed in human lung challenge studies with viable pathogen where virus was cultured from the URT until a median of 4 days (and a maximum of 10 days) from symptom onset (24). However, hardly anything is known about the compartment-specific duration of actively replicating virus in the lower respiratory tract (LRT), particularly in acute severely ill hospitalised patients undergoing mechanical ventilation. We hypothesised that there is compartment-specific uncoupling of viral replication in severe COVID-19 i.e. replicating virus can persist in the LRT beyond 10 days from symptoms onset, independent of its persistence in the URT, and this persistence may be associated with an altered pulmonary immunity.
Methods
Patients
The decedents (n=42) were recruited from Chris Hani Baragwanath Academic Hospital, Johannesburg, South Africa (n=18; Beta group) and Groote Schuur Hospital, Cape Town, South Africa (n=24; Delta group). Figure 1A outlines an overview of the study plan. Ambulatory controls (n=18) were recruited at diagnosis (baseline; ∼5 days from symptom onset), 7 days and 14 days post diagnosis. Minimally invasive tissue samples (MITS) and nasopharyngeal swabs from decedents (n=42) in the Beta and Delta waves (Figure 1B) were taken immediately after death. In addition, heart, liver, kidney, and adipose tissue samples were also taken from the Delta variant decedent cohort only. Ethical approval was obtained from the Human Research Ethics Committee (HREC) of the University of Cape Town (HREC approval number 866/2020) and University of Witwatersrand (HREC approval number M200313). Biosafety approvals were obtained from the Faculty Biosafety Committee of the University of Cape Town (IBC008-2021).
Nasopharyngeal (NP) swabs from ambulatory COVID-19 controls were obtained approximately 5 days after symptom onset (diagnosis), and then at 12- and 19-days post symptom onset. Minimally invasive tissue samples (MITS) and NP swabs were retrieved from decedents shortly after death.
Viral culture
To establish the in vitro viral culture model, a SARS-CoV-2 viral stock was used to infect the human lung carcinoma cell line, H1299 ACE2, in a BSL3 laboratory and infection was confirmed by light microscopy (as assessed by cytopathic effects of the virus on the cell line) and confocal microscopy (Figure S1A and B). Serial dilutions of the viral stock were used to establish the limit of detection of the PCR assay at 1×101 copies/ml (Figure S1C). Viral culture was performed on the nasopharyngeal swab and lung biopsy samples as indicated in the study overview (Figure 1) and detailed in the online supplement. Viral culture result reproducibility was good (see online supplement).
Multiplex PCR to detect secondary bacterial infections
The lung biopsy cores, stored in universal transport medium, were briefly homogenised and 200µl of the supernatant was applied to the BioFireÒ FilmArrayÒ Pneumonia panel (Bioméieux, South Africa). The panel was run using protocol BAL v3.3 according to the manufacturer’s instructions, thus generating RT-PCR readouts for 33 bacterial and viral pathogens. Bronchopneumonia was defined as histological evidence of a neutrophilic alveolar infiltration together with the detection of bacterial genomic material in the biopsy cores.
Immunohistochemistry
Immunohistochemical staining was performed using the Roche Ventana Automated platform (Ventana XT autostainer) as indicated by the manufacturer. Tissue sections were prepared, stained, and viewed using standard techniques (25). Antibodies included anti-CD3 (2GV6), anti-CD4 (SP35), anti-CD8 (SP57) and anti-CD68 (KP-1) (Roche USA).
Haematoxylin & eosin (H&E) staining and transmission electron microscopy (TEM)
H&E staining and TEM were performed according to standard procedures (25). H&E-stained slides were viewed using an Olympus BX43 microscope. TEM tissue sections were viewed using a Carl Zeiss EM109 microscope.
SARS-CoV-2 whole genome sequencing
Total SARS-CoV-2 RNA was extracted from lung biopsy samples and whole genome sequencing was performed. The generated reads were analysed with the Exatype (https://exatype.com) software to identify minor and major variants. The assembled consensus sequences were analysed using Nextclade Web (https://clades.nextstrain.org) for further quality control and clade assignment.
RNAseq
Total RNA was extracted from lung biopsy samples from the Delta group, sequenced and mapped consecutively to the human and COVID reference genomes using the Spliced Transcripts Alignment to a Reference (STAR) software [version 2.7.7a, (26)]. A differential expression (DE) analysis was performed on the generated raw read count file with the edgeR (Version 3.38.4) R package (27). The DE results were ranked by fold change and the gseGO function, from the clusterProfiler R clusterProfiler [Version 4.0, (28)] R package was used to perform a gene set enrichment analysis (GSEA) for the Gene Ontology Biological Process pathways. Pathways with an FDR <0.05 were considered significant.
Confocal microscopy
The H1299 ACE2 cells were plated, infected with SARS-CoV-2 and allowed to adhere to coverslips slides overnight at 37°C. The next day the cells were stained with or without anti-SARS-CoV-2 S1 spike protein (ThermoFisher, USA) and the slides were mounted in Mowiol (Calbiochem, USA) containing n-propyl gallate (Sigma-Aldrich, Germany) as an anti-fading agent. Confocal microscopy was performed with a Zeiss Axiovert 200M LSM510 Meta NLO Confocal Microscope.
Sample size calculation and statistical analysis
We hypothesised that we would detect lung culture positivity at 14 days post-symptom onset in ∼33% of decedents. A sample size of ∼40 participants would allow us to ascertain that level of positivity with a 15% margin of error using 95% confidence and 80% power (OpenEpi, Version 3, opensource calculator).
The Fisher Exact test was employed for categorical variables and for continuous variables, Mann-Whitney test was used for non-parametrically distributed data between the culture-negative and culture-positive groups (Stata version 17 or GraphPad, Version 9.4.1). A p-value of < 0.05 was considered significant for all statistical analyses.
The multivariable analysis was performed in R by fitting a binomial Generalized Linear Model (GLM) to assess the association between steroid use and the presence of secondary bacterial infection on culture status. The tidymodels (version, 1.0.0) R package was used to perform predictive modelling using the glm binomial classification algorithm. To account for the small sample size, 1000 bootstraps were performed for each analysis.
Results
Demographics and clinical characteristics of the decedents
The demographics of patients enrolled in the study are shown in Table S1. The median age of the patients was 53 years with 48% being males (20/42). 40.5% (17/42) had a secondary bacterial infection and 11% (4/38) had bacterial bronchopneumonia (microbiologically and histopathological confirmed). The median time from onset of symptoms to death, ICU admission to death and high flow oxygen admission to death was 17 (IQR; 9-22), 5 (2–12) and 11 (6–15) days, respectively.
SARS-CoV-2 replicating persistence in the human lung of mechanically ventilated decedents
We first ascertained the frequency and duration of replicating virus in lung tissue (which to our knowledge has not been previously undertaken). Culturable virus in the lung was present in 38.1% (16/42; Figure 2A) of mechanically ventilated ICU decedents, at a median of 15 days (and up ∼4 weeks; Figure S2) from symptom onset to sampling/death (Figure 3A). As expected, 56% (10/18) of a prospectively recruited control group of ambulatory patients had culturable virus, using nasopharyngeal swab samples, at day 5 from symptom onset (Figure 2B). In the same group of patients after 12- and 19-days after symptom onset, only 5.5% (1/18) and 0% (0/18), respectively, had culturable virus from their nasopharyngeal swab (Figure 2B). By contrast, 38% of nasopharyngeal swabs from the mechanically ventilated ICU descendants had culturable virus (Figure 2B), at a median of 13 days from symptom onset to sampling/death (Figure S3B). Additionally, SARS-CoV-2 could be detected by PCR in multiple organs in lung culture-positive decedents in the Delta cohort (biopsies other than the lung were not performed in the Beta cohort) suggesting widespread multi-organ viral dissemination (Figure 2C). SARS-CoV-2 was also detected in adipose tissue of culture-positive decedents (hitherto undescribed). Clinical characteristics, such as, age and comorbidities were similar in the lung culture-positive versus culture-negative groups (Table S1). We found no association between viral genetic variant and the phenotype of replicating viral persistence (although this might have been a factor of the limited sample size; Table S5).
(A) Proportion of lung biopsy samples that were culture-positive from the decedents. (B) Proportion of ambulatory patient and decedent NP swab samples that were culture-positive. (C) PCR positivity of organs of lung culture-positive decedents from the Delta cohort. NP= nasopharyngeal.
(A) The days from symptom onset to death for the culture-negative (-ve; green) and culture-positive (+ve; red) groups for both groups combined, and for the Beta (B) and Delta (C) groups alone. Proportion of samples/participants with a secondary bacterial infection in culture-negative and culture-positive decedents overall i.e. the combined groups (Beta and Delta) (D), Beta group only (E), and Delta group only (F). (G) PCR cycle threshold (Ct) value at the time of death or at admission could not discriminate or predict lung culture status. The nasopharyngeal swab PCR Ct values at admission or death were missing for some participants because they were either diagnostically confirmed by antigen testing or the Ct value was not recorded. Due to the nature of the pandemic and the burden of the disease on the healthcare infrastructure at the time, Ct values at peak periods were not recorded. We have conducted sensitivity and imputation analyses indicating that these missing data points are redundant. The nasopharyngeal swab PCR Ct values at admission or death were missing for some participants because they were either diagnostically confirmed by antigen testing or the Ct value was not recorded. Due to the nature of the pandemic and the burden of the disease on the healthcare infrastructure at the time, Ct values at peak periods were not recorded. We have conducted sensitivity and imputation analyses indicating that these missing data points are redundant.
Time to death in the Delta and Beta groups and predictors of lung culture positivity
Next, we evaluated variant-specific relationships to clinical outcomes. Mechanically ventilated patients who were SARS-CoV-2 lung culture-positive in the Delta, but not the Beta group, had accelerated death (i.e. shorter duration from symptom onset to death; Figure 3C versus 3B; p=0.004), and a higher proportion of lung-specific secondary bacterial infection (Figure 3F versus 3E; p=0.032) compared to culture-negative decedents. Similarly, to the lung culture data, the nasopharyngeal swab culture-positive Delta, but not the Beta group, had accelerated death (Figure S3D versus S3C; p=0.026).
The bacterial species identified from the lung biopsies of both the Beta and Delta groups included Streptococcus, Staphylococcus, Haemophilus, Acinetobacter, Proteus spp, Escherichia, Klebsiella, Enterobacter and Serratia (Table S4). Overall, both groups were infected with one or more bacteria that were sensitive or resistant to β-lactams and/or carbapenems (Table S4). Key clinical and demographic characteristics such as differences in co-morbidities (age, obesity, diabetes, HIV positivity etc; Table S1) associated as drivers of severe COVID-19 disease and poor prognosis, could not explain these observations, despite the lower population-level vaccination and pre-existing COVID-19 exposure rates in the Beta cohort. Steroid usage was similar in the culture-positive and culture-negative groups (Table S1) and there was no significant (p>0.05) association between steroid use and the presence of secondary bacterial infection in a multivariable analysis.
Next, we interrogated whether nasopharyngeal PCR characteristics (Ct value), either at admission or close to death, could identify the phenotype of lung replicating viral persistence. However, nasopharyngeal Ct neither at admission, nor at the time of death was associated with lung culture positivity (Figure 3G). This suggests that the kinetics of viral replication was different in the upper and the lower respiratory tract.
Lung immunity and histology of the culture-negative versus culture-positive groups
We then ascertained whether the phenotype of replicating viral persistence was associated with attenuated or modulated lung immunity in the Delta decedents (transcriptomic and flow cytometric studies were only carried out in Cape Town, i.e the Delta decedents, due to location-specific availability of assays and limited Beta group biopsy cores that had been used for unrelated studies). Immunohistochemical staining indicated that there was significantly less infiltration of CD3+ T-cells, specifically CD8+ T-cells in the alveoli and interstitium of the SARS-CoV-2 culture-positive compared with the culture-negative individuals in the Delta decedents (Figure 4A and B).
(A) More CD3+ and CD8+ T-cells infiltrate into the alveoli and interstitial space of the lung culture-negative versus culture-positive group in the Delta decedents as assessed by immunohistochemistry. (B) Representative images (immunohistochemistry) showing increased T-cell infiltration into interstitial space (blue arrow) in the lung culture-negative versus the culture-positives in the Delta cohort. Histopathology findings (C) and representative images (D) associated with diffuse alveolar damage and microvascular thrombosis in the Delta decedents. The black arrows indicate key histopathological features.
The typical histological features of severe COVID-19 (e.g. diffuse alveolar damage and microvascular thrombosis) were similar in the SARS-CoV-2 culture-positive and the culture-negative phenotype suggesting that these events occurred in the early rather than the persistent viral replication phase (Figure 4C, D and S4; Table S2, S3). Interestingly, we observed that some features of leucocyte hyperactivation (i.e., hemophagocytic syndrome) were more common in the SARS-CoV-2 culture-negative versus the culture-positive group, potentially in keeping with an aberrant immune response characterised by a lack of immune regulation, as outlined above (Figure 4C and 4D; p=0.013).
The transcriptional analysis of post-mortem lung tissue identified 630 up- and 885 downregulated genes at an uncorrected p-value ≤ 0.05 when comparing the culture-positive versus negative-groups (Figure S5). After adjustment for multiple testing, the same comparison identified a total of 11 up- and 4 down-regulated genes (FDR<0.05; specific genes discussed further in the online supplement; Figure S6).
The GSEA performed using the full list of differentially expressed genes ranked by fold-change, identified activated pathways that were associated with a proinflammatory response related to cytokine signalling, neutrophil and monocyte chemotaxis/recruitment, and viral entry/defence, all of which are implicated in COVID-19-related hypercytokinaemia (29) (Figure 5, S5, S6 and S7, Table S6 and S7A). Significantly repressed pathways were generally associated with body homeostasis (Figure 5A and B). There was also in tandem upregulation of Th1 and Th17 signalling pathways (Table S7B) but to a substantially lesser extent than that of innate cellular and signalling pathways (IL-1, IL-6 and neutrophil-related; Table S7A). These features may be consistent with an aberrant immune response including a lack of activation of regulatory and immune-suppressive pathways. T-cell exhaustion consistent with upregulation of PD-1, CTLA-4 and LAG (Table S7C) known to be associated with severe COVID-19, was not observed.
Dot plot illustrating the significantly activated and suppressed pathways along with the gene count and ratio for each pathway (A), enrichment map illustrating the significantly activated and suppressed pathways along with the gene count and ratio for each pathway (B) and cnetplot similarly illustrating the overlap of genes and their fold changes for selected activated pathways (C).
The differential expression (DE) results also revealed that a number of SARS-CoV-2 genes were significantly upregulated (FDR<0.01) in the culture-positive versus the culture-negative group including nucleocapsid phosphoprotein (log2 FC=8.4) and ORF3a (log2 FC=5.5) while the surface/spike glycoprotein encoding gene had a log2 FC of 5.3 and an FDR of 0.067 (Table S6). A visual inspection of the mapped SARS-CoV-2 reads revealed that those that mapped to the 5’ end of the genes were spliced with a portion mapping to the 5’ leader sequence of the genome. This suggests the reads originated from sub-genomic mRNA (sgRNA) rather than genomic RNA which is consistent with the active viral replication observed in the culture-positive group.
Finally, we evaluated whether any of the DE genes could act as biomarkers discriminating between lung culture-positive and negative-individuals. Logistic regression predictive modelling revealed that GREM1 and FGFBP1 were associated with a sensitivity and specificity above 90% (further details in the online supplement). Future studies are warranted to determine if these lung-based biomarkers can predict patient culture status in blood samples.
Discussion
The widely accepted view in severe acute COVID-19 is that resolution of the initial viral replication phase in the first week after symptom onset is followed by an effector or hyperinflammatory phase in the second and third week of illness which is characterised by diffuse alveolar damage, thrombo-inflammation, and endotheliopathy (30). Indeed, the Infectious Disease Society of America (IDSA) recommends the use of remdesivir for only 5 days in patients with severe illness and not in mechanically ventilated patients (31). However, our results, based on post-mortem lung biopsies obtained using minimally invasive tissue sampling methods (MITS) shortly after death indicated that, in contradistinction to the URT where replication often ceases within ∼8 days from symptom onset, in the human lung virus is culturable in ∼40% of mechanically ventilated patients until death (median of 15 days and up to 4 weeks after symptom onset; see Figure 1 for the study overview). The culture-positive group in the Delta cohort had accelerated death and a higher proportion of secondary bacterial infection in the lung compared to the culture-negative group. Nasopharyngeal SARS-CoV-2 viral load (based on Ct value) neither at admission nor at death, was predictive of lung culture-positivity. SARS-CoV-2 culture-positivity in the lung of decedents was associated with attenuated pulmonary T-cell immunity and an exaggerated pro-inflammatory phenotype. Importantly, this was concurrent with, rather than sequential to the viral replication or viral culture-positive phase.
These findings challenge the traditional paradigm of an initial viral replicative phase in the first week of severe illness sequentially followed by an effector or inflammatory phase (30). Our data suggest that in ∼40% of ventilated patients, viral replication persisted until death (i.e. 3rd and 4th week of illness and a median of 15 days after symptom onset) compared to ∼2 to ∼8 days in the URT as outlined in several studies including a live virus human challenge study (11, 14–24). One outlier study reported culturing virus from the URT for up to 3 weeks after symptom onset (32). However, a large proportion of patients were immunocompromised, samples at diagnosis and follow-up were combined (skewing the results), a high proportion of participants were healthcare workers (re-infection may have been a confounder), and as the authors suggested a limitation was that the Vero cell line used was overtly permissive to infection compared to the human lung carcinoma cell line, H1299 ACE2, which is a biologically representative cell line (and one that we used). Another recent study showed that infectious virus production peaked in the human lung within 2 days, but this model used ex vivo agarose infused devascularised and explanted human lung slices, which are not representative of what is occurring in freshly harvested human lung (33). The culture-based findings in the afore-mentioned studies must be explicitly distinguished from studies that detected residual free viral genomic RNA (but not replicating virus) embedded in the respiratory tract tissue of patients that had severe disease for an extended period of time (13, 34, 35). Indeed, SARS-CoV-2 RNAs have been detected in patient tissue many months after recovery from acute infection (36–38). It was initially suggested that sgRNA (sub-genomic RNA; small strands of reversely transcribed RNA) could be used as a proxy to infer viral replication. However, several recent studies have indicated that it has poor predictive value as a proxy for viral replication (39, 40). Indeed, Stein et al (13) detected sgRNA in multiple post-mortem organ biopsies, including the brain, several months after symptom onset. Thus, the data presented in this manuscript is the first to do so conclusively and comprehensively using viral culture from lung tissue (the gold standard to detect viral replication) beyond two weeks after symptom onset.
We demonstrated active viral replication in the lungs of acutely ill ventilated patients for up to ∼4 weeks after symptom onset. This challenges the current practice of using antivirals like remdesivir for only 5 days and suggests that a longer duration of treatment may be required in critically ill patients. Furthermore, antivirals like remdesivir are not recommended for use by IDSA in mechanically ventilated patients (conditional recommendation) as they felt that such patients (often in the third week of their illness) are no longer in the viral replicative phase, and published controlled trial data showed no mortality benefit of remdesivir in such patients (31, 41). However, these studies demonstrated a group effect, and the analyses did not adjust for disease severity or the time from symptom onset to death in mechanically ventilated patients (42). Our data suggest that a significant number of patients may likely benefit from antivirals during mechanical ventilation. Indeed, several observational studies have found a survival benefit using remdesivir in mechanically ventilated patients, but this requires further clarification in appropriate trials (42–44). In a multivariate analysis we found no association between steroid usage and lung viral culture positivity - in fact, steroid usage was lower in the viral culture-positive group and often culture positivity persisted beyond the 10 days of steroid usage.
The transcriptomic data suggested that, in a significant number of patients, the hyperinflammatory and viral replication phase occur concurrently in the 3rd and 4th week of illness, in contradistinction to the widely held view that these are sequential phases. Antiviral and selective proinflammatory responses were over-represented in the SARS-CoV-2 culture-positive compared with the culture-negative decedents, and we did not detect attenuated type 1 interferon responses at the site of disease compared with other reports (45–49). Three prior studies (one that enrolled 5 COVID-19 patients) evaluated transcriptomic lung responses in patients with severe COVID-19 versus healthy controls (45, 46, 50). These first level studies logically attempted to address the significance of transcriptomic changes specific to COVID-19 by using healthy controls or non-diseased parts of the lung from lung cancer patients. However, we specifically sought to compare culture-positive versus culture-negative groups (hitherto not undertaken) to dissect out pathways that facilitate permissiveness to ongoing viral replication.
We identified two lung-based biomarkers (GREM1 and FGFBP1) that could predict culture-positivity. Although these are lung-specific biomarkers, this preliminary analysis in a limited number of samples suggests that in the future, RT-PCR of tracheal aspirates or blood (if they are concordant with lung findings), could potentially serve as biomarkers to identify and direct appropriate treatment protocols to culture-positive persons but further investigation is needed.
There are several limitations to our findings. Firstly, our findings are relevant to acute severe COVID-19 ARDS/pneumonia requiring mechanical ventilation and may not be applicable to milder forms of disease seen in hospitalised patients or chronic infection seen in immunocompromised patients. Secondly, we only studied patients with the Beta and Delta variants as these were the predominant variant at the time of the study. However, Omicron has also been associated with severe disease in several settings including the surge of severe COVID-19 unfolding in China. Thirdly, we did not study a control group comprising severe ARDS due to other causes because our express aim was to investigate the presence and duration of viral replication in the LRT in severe COVID-19 disease. Fourth, the sample size limited our ability to make conclusions about several aspects. However, the highly resource intensive and demanding nature of the study limited our ability to recruit higher numbers of participants. Fifthly, it could be suggested that there may have been sampling error and variability of the viral culture assay. However, the reproducibility of the viral culture technique using 6 samples across 2 separate runs had a low standard error, which was indicative of high reproducibility. Finally, the transcriptional signature and flow cytometric findings may have been affected by post-death sampling, but several detailed studies have shown (51) that most protein and RNA species are preserved and stable for several hours after death. Given that biopsies for the transcriptional studies were taken ∼2 hours after death, we feel they are broadly representative of the picture at the time of death.
In summary, our data suggest that in COVID-19 disease there is considerable heterogeneity in the frequency and duration of viral replication in the upper versus the lower respiratory tract (i.e. lungs) beyond the 2nd week of illness, and that in a significant proportion of seriously ill patients, persisting viral replication occurs concurrently and may drive an exaggerated proinflammatory response (higher than in culture-negative persons), rather than sequentially as it is widely believed. These findings have potential implications for the use of antiviral therapy in seriously ill patients with COVID-19 and suggest that better biomarkers are needed to identify patient phenotypes and subsets that might benefit from concurrent anti-inflammatory and antiviral therapy.
Data Availability
Individual participant data will be made available to researchers who provide a protocol that is approved by their respective human research ethics committee. All protocols will be reviewed and approved by the MITS consortium trial steering committee up to five years following publication. A data sharing agreement (DTA) will need to be concluded between the representatives of the requesting institution and the University of Cape Town Lung Institute. Data sharing requests should be directed to keertan.dheda@uct.ac.za
Contributions
K.D, MT, S.J, A.P, A.E, M.D, M.N and S.M conceived and designed experiments. K.D, A.E, S.O, L.K, T.P, A.S, I.J, J.P, D.T, D.F, M.M, M.N and S.M arranged medical ethical approval, recruitment of study participants and collection of study material. MT, SJ. A.P, L.W, A.K, S.M, K.P, R.R, R.K, R.M, C.M, J.B, R.L, M.D, C vdM, A.R, M.H, V.B, S. M, C.W, and R.J performed the experiments. M.T, S.J, A.P and M.D set up experimental assays. A.S provided the cell line. K.D, M.T, S.J, A.P, L.W, A.K, L.S, S.M, K.P, R.R, R.L, M.H, V.B, S.M, R.J, and N.M analysed and interpreted data. K.D, M.T, L.S, and S.M wrote the manuscript with input from all listed authors.
Declaration of interests
The authors have no competing interests.
Acknowledgments
The authors would like to thank Arnold-Day C, Crowther M, Fernandes N, and Mitchell L from the Division of Critical Care, Department of Anaesthesia and Perioperative Medicine, University of Cape Town, South Africa for identifying and referring potential patients into the study.
The study was funded by the Bill & Melinda Gates Foundation (grant number INV-017282) and the South African Medical Research Council (grant number SHIP NCD 96756) with partial support from the Department of Science and Technology and National Research Foundation: South African Research Chair Initiative in Vaccine Preventable Diseases. The KD lab acknowledges funding from the SA MRC (RFA-EMU-02-2017), EDCTP (TMA-2015SF-1043, TMA-1051-TESAIII, TMA-CDF2015), UK Medical Research Council (MR/S03563X/1), NIH (CRDF-OISE-16-62105) and the Wellcome Trust (MR/S027777/1). This work was co-funded by The Wellcome Centre for Infectious Diseases Research in Africa is supported by core funding from the Wellcome Trust (230135/Z/16/Z) and the European Union’s Horizon Europe Research and Innovation Actions (101046041) for genomic surveillance.
The authors would like to thank the families of the deceased who gave us permission to conduct the study, which may advance our understanding of SARS-CoV-2 pathogenesis to inform future clinical management of respiratory pathogen pandemics.
Footnotes
The manuscript was amended according to reviewers comments