Understanding severe coronavirus disease in humans from the analysis of clinical samples with RNA sequencing

Coronaviruses in humans have been of concern since the emergence of SARS-CoV, MERS-CoV and SARS-CoV-2 over the past two decades. Coronavirus disease in humans can range from asymptomatic to mild or severe where symptomatic disease is associated with fever, cough, and respiratory symptoms. Next generation sequencing and phylogeny studies can provide insight into viral evolution due to the nucleotide polymorphisms which arise due to the inherent error rates. Nanopore sequencing can facilitate these studies rapidly and are therefore a useful public health tool for genome surveillance. MERS-CoV emerged in Saudi Arabia in 2012 and is associated with sporadic outbreaks. To facilitate rapid genomic surveillance of MERS-CoV in Saudi Arabia, a PCR amplification sequencing method compatible with Nanopore was designed. This approach is useful as data derived from this methodology can be used for phylogeny and variant analysis which supports investigation into transmission events and viral evolution. Upon the emergence of SARS-CoV-2 in late 2019, this approach was then repurposed and utilised on a subset of patients from the UK. Alternative viral genome sequencing approaches were then employed to assess the tissue tropism of SARS-CoV-2 in fatal COVID-19 cases, where tropism was widespread, while inflammation was exclusive to pulmonary tissues. To complement the analysis conducted on the tissue of fatal COVID-19 patients, blood from patients at point of care were utilised for blood transcriptomics analysis. Both illumina and nanopore sequencing methodologies were employed to assess differences in transcript abundance in these patients. Transcriptomic profiles from COVID-19 patients were compared to profiles from Influenza patients and healthy controls, in addition to fatal and non-fatal COVID-19 cases. The key finding from this analysis was that immunoglobulin domains transcripts exhibited altered abundance when comparing COVID and Influenza patients, as well as between fatal and non-fatal COVID-19 cases. From this insight, it is hypothesised that an early adaptive immune response is associated with survival, or a delayed adaptive response is associated with fatality. As it is challenging to control variables from patient data, an hACE2 mice model was utilised to explore the host response against Influenza A virus and SARS-CoV-2 as independent and sequential infections using lung tissue. Transcriptomic analysis revealed a sustained cytokine and interferon response in coinfected mice. Transcripts changing in abundance in both human blood and mice lungs were compared to generate a subset of transcripts that may be essential to the response to SARS-CoV-2 infection. In summary, the results described within this thesis provide insights into the novel coronavirus SARS-CoV-2 and COVID-19 disease in humans. Additionally, the outputs of this thesis provide a foundation for further investigation and development of Nanopore sequencing methodologies as a prognostic tool

The Stereotypic Response of the Pulmonary Vasculature to Respiratory Viral Infections: Findings in Mouse Models of SARS-CoV-2, Influenza A and Gammaherpesvirus Infections

The respiratory system is the main target of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of coronavirus disease 19 (COVID-19) where acute respiratory distress syndrome is considered the leading cause of death. Changes in pulmonary blood vessels, among which an endothelialitis/endotheliitis has been particularly emphasized, have been suggested to play a central role in the development of acute lung injury. Similar vascular changes are also observed in animal models of COVID-19. The present study aimed to determine whether the latter are specific for SARS-CoV-2 infection, investigating the vascular response in the lungs of mice infected with SARS-CoV-2 and other respiratory viruses (influenza A and murine gammaherpesvirus) by in situ approaches (histology, immunohistology, morphometry) combined with RNA sequencing and bioinformatic analysis. Non-selective recruitment of monocytes and T and B cells from larger muscular veins and arteries was observed with all viruses, matched by a comparable transcriptional response. There was no evidence of endothelial cell infection in any of the models. Both the morphological investigation and the transcriptomics approach support the interpretation that the lung vasculature in mice mounts a stereotypic response to alveolar and respiratory epithelial damage. This may have implications for the treatment and management of respiratory disease in humans

Analysis of SARS-CoV-2 in Nasopharyngeal Samples from Patients with COVID-19 Illustrates Population Variation and Diverse Phenotypes, Placing the Growth Properties of Variants of Concern in Context with Other Lineages

New variants of SARS-CoV-2 are continuing to emerge and dominate the global sequence landscapes. Several variants have been labeled variants of concern (VOCs) because they may have a transmission advantage, increased risk of morbidity and/or mortality, or immune evasion upon a background of prior infection or vaccination. Placing the VOCs in context with the underlying variability of SARS-CoV-2 is essential in understanding virus evolution and selection pressures. Dominant genome sequences and the population genetics of SARS-CoV-2 in nasopharyngeal swabs from hospitalized patients were characterized. Nonsynonymous changes at a minor variant level were identified. These populations were generally preserved when isolates were amplified in cell culture. To place the Alpha, Beta, Delta, and Omicron VOCs in context, their growth was compared to clinical isolates of different lineages from earlier in the pandemic. The data indicated that the growth in cell culture of the Beta variant was more than that of the other variants in Vero E6 cells but not in hACE2-A549 cells. Looking at each time point, Beta grew more than the other VOCs in hACE2-A549 cells at 24 to 48 h postinfection. At 72 h postinfection there was no difference in the growth of any of the variants in either cell line. Overall, this work suggested that exploring the biology of SARS-CoV-2 is complicated by population dynamics and that these need to be considered with new variants. In the context of variation seen in other coronaviruses, the variants currently observed for SARS-CoV-2 are very similar in terms of their clinical spectrum of disease. IMPORTANCE SARS-CoV-2 is the causative agent of COVID-19. The virus has spread across the planet, causing a global pandemic. In common with other coronaviruses, SARS-CoV-2 genomes can become quite diverse as a consequence of replicating inside cells. This has given rise to multiple variants from the original virus that infected humans. These variants may have different properties and in the context of a widespread vaccination program may render vaccines less effective. Our research confirms the degree of genetic diversity of SARS-CoV-2 in patients. By comparing the growth of previous variants to the pattern seen with four variants of concern (VOCs) (Alpha, Beta, Delta, and Omicron), we show that, at least in cells, Beta variant growth exceeds that of Alpha, Delta, and Omicron VOCs at 24 to 48 h in both Vero E6 and hACE2-A549 cells, but by 72 h postinfection, the amount of virus is not different from that of the other VOCs

Tissue proteomic analysis identifies mechanisms and stages of immunopathology in fatal COVID-19

Funding: This work was funded by UK Research and Innovation (UKRI) (Coronavirus Disease [COVID-19] Rapid Response Initiative; MR/V028790/1 to C.D.L., D.A.D., and J.A.H.), LifeArc (through the University of Edinburgh STOPCOVID funding award, to K.D, D.A.D., C.D.L), The Chief Scientist Office (RARC-19 Funding Call, ‘Inflammation in Covid-19: Exploration of Critical Aspects of Pathogenesis; COV/EDI/20/10’ to D.A.D, C.D.L, C.D.R, J.K.B and D.J.H), and Medical Research Scotland (CVG-1722- 2020 to DAD, CDL, CDR, JKB, and DJH). C.D.L is funded by a Wellcome Trust Clinical Career Development Fellowship (206566/Z/17/Z). J.K.B. and C.D.R. are supported by the Medical Research Council (grant MC_PC_19059) as part of the ISARIC Coronavirus Clinical Characterisation Consortium (ISARIC-4C). C.D.R. is supported by an Edinburgh Clinical Academic Track (ECAT)/Wellcome Trust PhD Training Fellowship for Clinicians award (214178/Z/18/Z). J.A.H. is supported by the U.S. Food and Drug Administration (contract 75F40120C00085, Characterization of severe coronavirus infection in humans and model systems for medical countermeasure development and evaluation’). G.C.O is funded by an NRS Clinician award. N.N.G. is funded by a Pathological Society Award. A.R.A. is supported by a Cancer Research UK Clinician Scientist Fellowship award (A24867).Immunopathology occurs in the lung and spleen in fatal COVID-19, involving monocytes/macrophages and plasma cells. Anti-inflammatory therapy reduces mortality but additional therapeutic targets are required. We aimed to gain mechanistic insight into COVID-19 immunopathology by targeted proteomic analysis of pulmonary and splenic tissues. Lung parenchymal and splenic tissue was obtained from 13 post-mortem examinations of patients with fatal COVID-19. Control tissue was obtained from cancer resection samples (lung) and deceased organ donors (spleen). Protein was extracted from tissue by phenol extraction. Olink® multiplex immunoassay panels were used for protein detection and quantification. Proteins with increased abundance in the lung included MCP-3, antiviral TRIM21 and pro-thrombotic TYMP. OSM and EN-RAGE/S100A12 abundance was correlated, and associated with inflammation severity. Unsupervised clustering identified ‘early viral’ and ‘late inflammatory’ clusters with distinct protein abundance profiles, and differences in illness duration prior to death and presence of viral RNA. In the spleen, lymphocyte chemotactic factors and CD8A were decreased in abundance, and pro-apoptotic factors were increased. B-cell receptor signalling pathway components and macrophage colony stimulating factor (CSF-1) were also increased. Additional evidence for a sub-set of host factors (including DDX58, OSM, TYMP, IL-18, MCP-3 and CSF-1) was provided by overlap between (i) differential abundance in spleen and lung tissue, (ii) meta-analysis of existing datasets, and (iii) plasma proteomic data. This proteomic analysis of lung parenchymal and splenic tissue from fatal COVID-19 provides mechanistic insight into tissue anti-viral responses, inflammation and disease stages, macrophage involvement, pulmonary thrombosis, splenic B-cell activation and lymphocyte depletion.Publisher PDFPeer reviewe

Enrichment of SARS-CoV-2 sequence from nasopharyngeal swabs whilst identifying the nasal microbiome

Simultaneously characterising the genomic information of coronaviruses and the underlying nasal microbiome from a single clinical sample would help characterise infection and disease. Metatranscriptomic approaches can be used to sequence SARS-CoV-2 (and other coronaviruses) and identify mRNAs associated with active transcription in the nasal microbiome. However, given the large sequence background, unenriched metatranscriptomic approaches often do not sequence SARS-CoV-2 to sufficient read and coverage depth to obtain a consensus genome, especially with moderate and low viral loads from clinical samples. In this study, various enrichment methods were assessed to detect SARS-CoV-2, identify lineages and define the nasal microbiome. The methods were underpinned by Oxford Nanopore long-read sequencing and variations of sequence independent single primer amplification (SISPA). The utility of the method(s) was also validated on samples from patients infected seasonal coronaviruses. The feasibility of profiling the nasal microbiome using these enrichment methods was explored. The findings shed light on the performance of different enrichment strategies and their applicability in characterising the composition of the nasal microbiome

Sequential Infection with Influenza A Virus Followed by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Leads to More Severe Disease and Encephalitis in a Mouse Model of COVID-19

COVID-19 is a spectrum of clinical symptoms in humans caused by infection with SARS-CoV-2. The coalescence of SARS-CoV-2 with seasonal respiratory viruses, particularly influenza viruses, is a global health concern. To understand this, transgenic mice expressing the human ACE2 receptor (K18-hACE2) were infected with influenza A virus (IAV) followed by SARS-CoV-2 and the host response and effect on virus biology was compared to K18-hACE2 mice infected with IAV or SARS-CoV-2 alone. The sequentially infected mice showed reduced SARS-CoV-2 RNA synthesis, yet exhibited more rapid weight loss, more severe lung damage and a prolongation of the innate response compared to the singly infected or control mice. Sequential infection also exacerbated the extrapulmonary encephalitic manifestations associated with SARS-CoV-2 infection. Conversely, prior infection with a commercially available, multivalent live-attenuated influenza vaccine (Fluenz Tetra) elicited the same reduction in SARS-CoV-2 RNA synthesis, albeit without the associated increase in disease severity. This suggests that the innate immune response stimulated by IAV inhibits SARS-CoV-2. Interestingly, infection with an attenuated, apathogenic influenza vaccine does not result in an aberrant immune response and enhanced disease severity. Taken together, the data suggest coinfection (‘twinfection’) is deleterious and mitigation steps should be instituted as part of the comprehensive public health and management strategy of COVID-19

Emerging variants of canine enteric coronavirus associated with seasonal outbreaks of severe canine gastroenteric disease

Canine enteric coronavirus (CECoV) variants have an emerging role in severe outbreaks of canine gastroenteritis. Here we used syndromic health data from a sentinel network of UK veterinary practices to identify an outbreak of severe canine gastroenteritis. Affected dogs frequently presented with vomiting, diarrhoea and inappetence. Data from sentinel diagnostic laboratories showed similar seasonal increases in CECoV diagnosis. Membrane glycoprotein (M) gene sequence analysis implied wide geographical circulation of a new CECoV variant. Whole genome sequencing suggested the main circulating 2022 variant was most closely related to one previously identified in 2020 with additional spike gene recombination; all variants were unrelated to CECoV-like viruses recently associated with human respiratory disease. Identifying factors that drive population-level evolution, and its implications for host protection and virulence, will be important to understand the emerging role of CECoV variants in canine and human health, and may act as a model for coronavirus population adaptation more widely

The P323L Substitution in the SARS-CoV-2 Polymerase (NSP12) Confers a Selective Advantage During Infection

BACKGROUND: The mutational landscape of SARS-CoV-2 varies at the dominant viral genome sequence and minor genomic variant population. During the COVID-19 pandemic, an early substitution in the genome was the D614G change in the spike protein, associated with an increase in transmissibility. Genomes with D614G are accompanied by a P323L substitution in the viral polymerase (NSP12). However, P323L is not thought to be under strong selective pressure. RESULTS: Investigation of P323L/D614G substitutions in the population shows rapid emergence during the containment phase and early surge phase during the first wave. These substitutions emerge from minor genomic variants which become dominant viral genome sequence. This is investigated in vivo and in vitro using SARS-CoV-2 with P323 and D614 in the dominant genome sequence and L323 and G614 in the minor variant population. During infection, there is rapid selection of L323 into the dominant viral genome sequence but not G614. Reverse genetics is used to create two viruses (either P323 or L323) with the same genetic background. L323 shows greater abundance of viral RNA and proteins and a smaller plaque morphology than P323. CONCLUSIONS: These data suggest that P323L is an important contribution in the emergence of variants with transmission advantages. Sequence analysis of viral populations suggests it may be possible to predict the emergence of a new variant based on tracking the frequency of minor variant genomes. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions

Emerging Variants of Canine Enteric Coronavirus Associated with Outbreaks of Gastroenteric Disease.

A 2022 canine gastroenteritis outbreak in the United Kingdom was associated with circulation of a new canine enteric coronavirus closely related to a 2020 variant with an additional spike gene recombination. The variants are unrelated to canine enteric coronavirus-like viruses associated with human disease but represent a model for coronavirus population adaptation

Evaluating the Immune Response in Treatment-Naive Hospitalised Patients With Influenza and COVID-19

The worldwide COVID-19 pandemic has claimed millions of lives and has had a profound effect on global life. Understanding the body’s immune response to SARS-CoV-2 infection is crucial in improving patient management and prognosis. In this study we compared influenza and SARS-CoV-2 infected patient cohorts to identify distinct blood transcript abundances and cellular composition to better understand the natural immune response associated with COVID-19, compared to another viral infection being influenza, and identify a prognostic signature of COVID-19 patient outcome. Clinical characteristics and peripheral blood were acquired upon hospital admission from two well characterised cohorts, a cohort of 88 patients infected with influenza and a cohort of 80 patients infected with SARS-CoV-2 during the first wave of the pandemic and prior to availability of COVID-19 treatments and vaccines. Gene transcript abundances, enriched pathways and cellular composition were compared between cohorts using RNA-seq. A genetic signature between COVID-19 survivors and non-survivors was assessed as a prognostic predictor of COVID-19 outcome. Contrasting immune responses were detected with an innate response elevated in influenza and an adaptive response elevated in COVID-19. Additionally ribosomal, mitochondrial oxidative stress and interferon signalling pathways differentiated the cohorts. An adaptive immune response was associated with COVID-19 survival, while an inflammatory response predicted death. A prognostic transcript signature, associated with circulating immunoglobulins, nucleosome assembly, cytokine production and T cell activation, was able to stratify COVID-19 patients likely to survive or die. This study provides a unique insight into the immune responses of treatment naïve patients with influenza or COVID-19. The comparison of immune response between COVID-19 survivors and non-survivors enables prognostication of COVID-19 patients and may suggest potential therapeutic strategies to improve survival