Skip to main content

SARS-CoV-2 vaccine ChAdOx1 nCoV-19 infection of human cell lines reveals low levels of viral backbone gene transcription alongside very high levels of SARS-CoV-2 S glycoprotein gene transcription

Abstract

Background

ChAdOx1 nCoV-19 is a recombinant adenovirus vaccine against SARS-CoV-2 that has passed phase III clinical trials and is now in use across the globe. Although replication-defective in normal cells, 28 kbp of adenovirus genes is delivered to the cell nucleus alongside the SARS-CoV-2 S glycoprotein gene.

Methods

We used direct RNA sequencing to analyse transcript expression from the ChAdOx1 nCoV-19 genome in human MRC-5 and A549 cell lines that are non-permissive for vector replication alongside the replication permissive cell line, HEK293. In addition, we used quantitative proteomics to study over time the proteome and phosphoproteome of A549 and MRC5 cells infected with the ChAdOx1 nCoV-19 vaccine.

Results

The expected SARS-CoV-2 S coding transcript dominated in all cell lines. We also detected rare S transcripts with aberrant splice patterns or polyadenylation site usage. Adenovirus vector transcripts were almost absent in MRC-5 cells, but in A549 cells, there was a broader repertoire of adenoviral gene expression at very low levels. Proteomically, in addition to S glycoprotein, we detected multiple adenovirus proteins in A549 cells compared to just one in MRC5 cells.

Conclusions

Overall, the ChAdOx1 nCoV-19 vaccine’s transcriptomic and proteomic repertoire in cell culture is as expected. The combined transcriptomic and proteomics approaches provide a detailed insight into the behaviour of this important class of vaccine using state-of-the-art techniques and illustrate the potential of this technique to inform future viral vaccine vector design.

Background

Since the emergence of SARS-CoV-2 in late 2019, there has been a global effort to develop effective vaccines and at least four different adenovirus-vectored vaccines are at various stages of development or deployment [1,2,3,4,5]. Adenoviruses are non-enveloped viruses containing a linear double-stranded DNA genome of approximately 36 kbp that efficiently deliver their DNA genome to the nuclei of host cells for viral genome replication. Adenoviruses of mammals are in the genus Mastadenovirus, and these are grouped into one of 7 species (A–G) [6]. Adenovirus-based vectors typically have two regions of the virus genome removed, known as E1 and E3 [7]. The E1 region contains early genes required to trigger a transcription cascade enabling viral replication; E1-deleted vectors therefore need to be grown in E1 trans-complementing cell lines such as HEK293 cells [8]. HEK293 cells have a 4-kbp region of human adenovirus type 5 (HuAd5) integrated into the cellular genome that provides the E1 genes in trans enabling efficient virus vector replication and recombinant virus production. The E3 region is comprised of genes encoding proteins that primarily act to subvert the immune response to adenovirus infection and are thus not needed for replication in cell culture and potentially undesirable from a vaccine platform perspective. Usually, the transgene to be expressed is inserted into the virus genome in place of the E1 region under the control of a highly active promoter.

Adenovirus-based vaccines are grown in HEK293 cells (or equivalents) to very high titres, purified and then administered to individuals. The virus will attach to a host cell, enter and deliver the recombinant genome to the nucleus where the desired gene of interest is expressed. Typically, the lack of E1-encoded proteins should ensure that there is no productive viral replication. However, it is noteworthy that removal of the E1 and E3 regions still leaves a number of other viral genes and promoters, including the E2 and E4 regions/promoters as well as the major late promoter, that normally drives expression of transcripts covering the major structural proteins of adenovirus (e.g. Fig. 1a). Moreover, background expression of at least some of these adenovirus genes has been previously identified in studies with human adenovirus-based vectors [9,10,11,12]. Expression of these viral vector genes may lead to the production of adenoviral antigens which in turn would bring the cell to the attention of the adaptive immune response [13,14,15]. More directly, pre-existing antibodies can be problematic, for example by neutralising the incoming vaccine vector. This has prompted the development of adenovirus vectors based on strains of adenovirus for which there is little or no pre-existing immunity in the general human population [16].

Fig. 1
figure 1

Transcription map of detected transcripts in non-permissive and permissive human cells. a A simplified schematic of the S glycoprotein expression cassette as inserted into the ChAdOx1 nCoV-19 vector backbone in place of the E1 region of ChAdOx1. The nucleotide positions of major features in the expression cassette are indicated. Black transcripts originate from the top strand of the dsDNA viral genome (illustrated by two black lines with circles at the 5′ ends) and blue transcripts originate from the complementary strand. Each arrow depicts the expected transcription start and polyadenylation site but omits splicing events within each transcript that give rise to mature mRNA. Most ChAdOx1 structural proteins are transcribed starting at the same promoter (known as the major late promoter) ending at one of 5 different polyadenylation sites (known as L1–L5). However, the 22K gene can also be expressed from its own promoter as shown. Typically, every member of this late transcript group has three obligate exons at the beginning known as the tripartite leader. In addition, an occasional 4th exon is included containing an ORF known as the i-leader protein. In each subsequent part (b to e), the drawing from IGV viewer illustrates the structure of the dominant transcripts expressed from the 35,539 bp ChAdOx1 nCoV-19 genome that codes for each indicated ORF as detected by dRNAseq in MRC5 cells (b), A549 cells (c) and 293 cells (d). For the E4 region, we have marked with an asterisk those E4 ORFs that are in fact derived from HuAd5. The rectangles indicate exons joined by lines indicating introns. In each case, the dominant transcript for the S glycoprotein of SARS-CoV-2 is indicated in red. Transcripts that map to the top strand of the virus genome are in black and those that map to the bottom strand (and are reversed in orientation) are in blue. e The structure of dominant transcripts derived from the 4323 bp of the HuAd5 E1 region integrated into the genome of HEK293 cells

The ChAdOx1 vector virus is derived from chimpanzee adenovirus Y25 and is deleted for E1 and E3 genes. The seroprevalence of clinically relevant neutralising antibodies against Y25 in humans has been assessed in random adult serum samples and was found to be absent in UK adults (n = 100) compared to 9% of Gambian adults (n = 57) [17]. In addition, parts of the E4 region of the Y25 virus have been replaced with the equivalent E4 regions from HuAd5 to facilitate growth in HEK293 cells [17]. We have previously published data using this vaccine platform in earlier studies [18,19,20,21,22], and we have more recently reported on the ChAdOx1 nCoV-19 vaccine’s effectiveness in a phase 3 clinical trial [23]. In this study, we wanted to determine in detail exactly which adenovirus genes on the ChAdOx1 genome were being expressed in non-permissive cell lines alongside the S glycoprotein from SARS-CoV-2 using the latest transcriptomic and proteomic approaches. We recently developed a pipeline to analyse and characterise transcriptomic data derived from direct RNA sequencing (dRNAseq) on Oxford Nanopore devices, and we have used this pipeline to analyse the transcriptome of adenovirus [24] and SARS-CoV-2 [25] during productive viral replication. Potentially of significance for adenovirus-based vaccine delivery systems, we showed that in a wild type human adenovirus replication cycle in human cells, there is a wide range of non-standard transcripts produced as a result of numerous atypical splicing and polyadenylation events [24]. Because of this, we also wanted to understand the transcriptomic repertoire of the S glycoprotein transcript in the ChAdOx1 vector backbone. This is critical because during SARS-CoV-2 replication the S glycoprotein transcript, like all SARS-CoV-2 transcripts, is made in the cytoplasm in the absence of host splicing and polyadenylation machinery [26]. Indeed, in the ChAdOx1 nCoV-19 vector, expression of the S glycoprotein is enhanced by inclusion of an intron between the transcription start site and the S glycoprotein ORF. Inclusion of introns in gene expression cassettes is thought to promote engagement of the nascent mRNA with the host splicing, polyadenylation and mRNA export machinery and has been previously shown to significantly enhance transgene protein expression in adenoviral vectors [27]. We wanted to understand if this addition of an intron before the S glycoprotein ORF could enable unwanted splicing events. We chose to examine the transcriptomic and proteomic repertoire of ChAdOx1 nCoV-19 in two lung-derived human cell lines. MRC5 cells are genetically normal, slow growing and with a finite lifespan (although not primary), and A549 cells are immortal and cancerous. Thus, the cell types are both human and lung derived but quite distinct in their growth properties, and we reasoned this could result in distinct and informative transcriptomic expression profiles from the same vector.

We determined that transcription of the S glycoprotein gene is, as expected, by far the dominant transcript generated in both non-permissive cell lines tested. Interestingly, we noted distinct transcriptomic repertoires from the ChAdOx1 vector were possible in the two non-permissive cell lines and that at a very low-level aberrant splicing and polyadenylation of the SARS-CoV-2 S glycoprotein transcript did occur in all settings. Proteomically, we were able to detect SARS-CoV-2 S glycoprotein and several ChAdOx1 proteins as well as using phosphoproteomics to confirm that the SARS-CoV-2 S protein is phosphorylated as we have previously reported [25]. Our findings are likely to be reflected in all adenovirus E1/E3-deleted replication-defective type vectors irrespective of the parent adenovirus and represent a comprehensive and informative analysis of the transcriptomic repertoire of this important class of virus vaccine vector.

Methods

Virus and cells

Human MRC-5 cells (a genetically normal male human lung fibroblast-like line), A549 cells (a human male lung epithelial-like continuous line derived from carcinomatous tissue) and HEK293 cells (human embryonic kidney epithelial cell line immortalised by human adenovirus E1 region) were obtained from the European Collection of Authenticated Cell Cultures (ECACC). The cells were cultured in DMEM supplemented with 10% foetal bovine serum, 100 U/ml penicillin and 100 μg/ml streptomycin. The ChAdOx1 vector and the ChAdOx1 nCoV-19 vaccine vector have been previously reported [1, 5, 17], but briefly, the vector is derived from the Y25 chimpanzee adenovirus with E1 and E3 regions deleted and E4 region ORF 4, ORF 6 and ORF6/7 replaced with their equivalent HuAd5 counterparts to improve vector yield. The S glycoprotein ORF inserted into the vector (Genbank accession number YP_009724390.1) has been human codon optimised and contains the human tissue plasminogen activator (tPA) signal peptide at the N-terminus. Transcription of the S glycoprotein is driven by a CMV immediate promoter containing intron A [28] and the tetracycline-regulated operon sequence, and the polyadenylation signal is the bovine growth hormone (BGH) polyadenylation sequence. After reaching confluence, the cells were infected with ChAdOx1 nCoV-19 at a multiplicity of 10 infectious units per cell to ensure infections were synchronous. Duplicate flasks of infected cells were harvested at 24, 48 and 72 h post-infection (hpi) for MRC5-5 and A549 cells. One flask was used to extract total RNA whilst the duplicate was used to extract total protein content. For the HEK293 cells, just one flask was infected for total RNA extraction.

RNA extraction and sequencing

Total RNA was extracted from the infected cells using TRIzol™ reagent (#15596026, Ambion) at 1 ml of reagent per 107 cells but with two additional washes of extracted RNA with 70% ethanol prior to storage at − 80 °C under 70% ethanol. Once resuspended in H2O, extracted RNA was immediately enriched for polyadenylated RNA and then immediately sequenced as we have described previously for dRNAseq of human adenovirus-infected cells [24]. As before, the SQK-RNA002 kits and MIN106D R9 version flow cells (Oxford Nanopore Technologies) were used following the manufacturer’s protocols exactly. Between 0.8 and 1.5 million QC-passed reads were typically obtained per flow cell over 48 h.

Data analysis and characterisation of viral transcripts

As before, our previously described ORF-centric data analysis pipeline was used to characterise the RNA derived from the ChAdOx1 nCoV-19 genome [24]. Briefly, the transcripts were mapped to the viral genome with minimap2 and the mapping data was used to try to identify commonly used transcription start and termination locations alongside the splice acceptor/donor sites. Once this is complete, the software then assigns each transcript to a “transcript group” depending on its pattern of transcription start, transcription termination and splice acceptor/donor sites and counts how many transcripts belong to each transcript group. This data is then used by a second in-house script to generate pseudo transcripts based on the genomic sequences alongside a table of features predicted to be present on the vector genome (e.g. ORFS and predicted transcription start sites, Additional file 1). For each transcript group, this script then determines which features (if any) are contained within this transcript group. In addition, it notes which one (if any) of the known ORFs noted in the features table are 5′ proximal for each transcript group. In this manner, transcripts with different structures that still code for the same given protein can be counted together. This script also produced GFF files describing the structure of each transcript group and also produces a GFF file describing the dominant transcript type coding for each ORF detailed in the table of features provided. Finally, the pipeline also produces a list of proteins that are 5′ proximal for any given transcript group but are not identical to the list of known proteins predicted from the features table, referred to as the “proteins not known” list.

Characterisation of human transcripts

For each dataset, minimap2 [29] was used to map the transcripts to a list of human transcripts from Ensembl covering all human transcripts but with the fasta header modified so that each transcript included ENSG, ENST and ENSP information (Additional file 2). In this way, dRNAseq transcripts mapped to any given human transcript could be unambiguously assigned to its gene group and to the protein coded by that human transcript. An in-house script counted the number of times dRNAseq transcripts mapped to each curated Ensemble transcript and compared transcript mapping abundance between the mock and each of the time points after adjusting for the total number of mapped transcripts at each time point. This was collated as a simple table of log2 fold changes between mapping abundance of the mock sample transcripts time point and the mapping abundance of the transcripts at each time point. Further statistical evaluation was not applied since there is currently no consensus on how to use dRNAseq data to robustly infer gene expression changes or even the minimum depth of reads required [30].

Total and phosphoproteome analysis

Protein lysates were prepared from MRC5 or A549 cells only following mock infection, or infection with ChAdOx1 nCoV-19 for 24 h, 48 h or 72 h as previously described [25].

Aliquots of 100 μg of each sample were digested with trypsin (2.5 μg trypsin per 100 μg protein; 37 °C, overnight) and labelled with Tandem Mass Tag (TMT) ten plex reagents according to the manufacturer’s protocol (Thermo Fisher Scientific, Loughborough, LE11 5RG, UK), and the labelled samples pooled according to cell line.

For the total proteome analysis, an aliquot of 50 μg of the pooled sample was desalted using a SepPak cartridge according to the manufacturer’s instructions (Waters, Milford, MA, USA). Eluate from the SepPak cartridge was evaporated to dryness and resuspended in buffer A (20 mM ammonium hydroxide, pH 10) prior to fractionation by high pH reversed-phase chromatography using an Ultimate 3000 liquid chromatography system (Thermo Fisher Scientific). In brief, the sample was loaded onto an XBridge BEH C18 Column (130 Å, 3.5 μm, 2.1 mm × 150 mm, Waters, UK) in buffer A and peptides eluted with an increasing gradient of buffer B (20 mM ammonium hydroxide in acetonitrile, pH 10) from 0 to 95% over 60 min. The resulting fractions (15 in total) were evaporated to dryness and resuspended in 1% formic acid prior to analysis by nano-LC MSMS using an Orbitrap Fusion Lumos mass spectrometer (Thermo Scientific).

For the phosphoproteome analysis, the remainder of the TMT-labelled pooled sample was also desalted using a SepPak cartridge (Waters, Milford, MA, USA). Eluate from the SepPak cartridge was evaporated to dryness and subjected to TiO2-based phosphopeptide enrichment according to the manufacturer’s instructions (Pierce). The flow-through and washes from the TiO2-based enrichment were then subjected to FeNTA-based phosphopeptide enrichment according to the manufacturer’s instructions (Pierce). The phospho-enriched samples were again evaporated to dryness and then resuspended in 1% formic acid prior to analysis by nano-LC MSMS using an Orbitrap Fusion Lumos mass spectrometer (Thermo Scientific).

High pH RP fractions (total proteome analysis) or the phospho-enriched fractions (phosphoproteome analysis) were further fractionated using an Ultimate 3000 nano-LC system in line with an Orbitrap Fusion Lumos mass spectrometer (Thermo Scientific). In brief, peptides in 1% (vol/vol) formic acid were injected onto an Acclaim PepMap C18 nano-trap column (Thermo Scientific). After washing with 0.5% (vol/vol) acetonitrile 0.1% (vol/vol) formic acid, peptides were resolved on a 250 mm × 75 μm Acclaim PepMap C18 reverse-phase analytical column (Thermo Scientific) over a 150-min organic gradient, using 7 gradient segments (1–6% solvent B over 1 min, 6–15% B over 58 min, 15–32%B over 58 min, 32–40%B over 5 min, 40–90%B over 1 min, held at 90%B for 6 min and then reduced to 1%B over 1 min) with a flow rate of 300 nl min−1. Solvent A was 0.1% formic acid and solvent B was aqueous 80% acetonitrile in 0.1% formic acid. Peptides were ionised by nano-electrospray ionisation at 2.0 kV using a stainless-steel emitter with an internal diameter of 30 μm (Thermo Scientific) and a capillary temperature of 300 °C.

All spectra were acquired using an Orbitrap Fusion Lumos mass spectrometer controlled by Xcalibur 3.0 software (Thermo Scientific) and operated in data-dependent acquisition mode using an SPS-MS3 workflow. FTMS1 spectra were collected at a resolution of 120,000, with an automatic gain control (AGC) target of 400,000 and a max injection time of 100 ms. Precursors were filtered with an intensity threshold of 5000, according to charge state (to include charge states 2–7) and with monoisotopic peak determination set to peptide. Previously interrogated precursors were excluded using a dynamic window (60 s ± 10 ppm). The MS2 precursors were isolated with a quadrupole isolation window of 0.7 m/z. ITMS2 spectra were collected with an AGC target of 10,000, max injection time of 70 ms and CID collision energy of 35%.

For FTMS3 analysis, the Orbitrap was operated at 30,000 resolution with an AGC target of 50,000 and a max injection time of 105 ms. Precursors were fragmented by high-energy collision dissociation (HCD) at a normalised collision energy of 60% to ensure maximal TMT reporter ion yield. Synchronous precursor selection (SPS) was enabled to include up to 5 MS2 fragment ions in the FTMS3 scan.

Proteomics data analysis

The raw data files were processed using Proteome Discoverer software v2.1 (Thermo Scientific) and searched against a bespoke human database with the fasta headers amended to contain ENSG, ENST and ENSP notation of target human protein concatenated with a custom ChAdOx1 nCoV-19 protein database consisting of a list of “known proteins” predicted to exist by homology with known adenovirus proteins alongside the “proteins not known” list predicted by the transcriptome analysis pipeline described above (Additional file 3). Searches against this database and against an in-house “common contaminants” database were performed using the SEQUEST HT algorithm. Peptide precursor mass tolerance was set at 10 ppm, and MS/MS tolerance was set at 0.6 Da. Search criteria included oxidation of methionine (+ 15.995 Da), acetylation of the protein N-terminus (+ 42.011 Da) and methionine loss plus acetylation of the protein N-terminus (− 89.03 Da) as variable modifications and carbamidomethylation of cysteine (+ 57.021 Da) and the addition of the TMT mass tag (+ 229.163) to peptide N-termini and lysine as fixed modifications. For the phosphoproteome analysis, phosphorylation of serine, threonine and tyrosine (+ 79.966 Da) was also included as a variable modification. Searches were performed with full tryptic digestion and a maximum of 2 missed cleavages were allowed. The reverse database search option was enabled and all data was filtered to satisfy a false discovery rate (FDR) of 1%. Quantitative changes in phosphopeptide abundance were normalised with respect to total protein abundance.

Data integration

Once the quantitative data for proteome, phosphoproteome and transcriptome were collected, an in-house script was used to integrate the data producing a seamless dataset where for any given gene, data is provided on the number of transcripts mapping to that gene at any time point alongside log2 fold changes in protein or phosphoprotein abundance (if available). For the purposes of noting changes in phosphopeptide abundance, if multiple phosphorylation sites were identified for any given protein, in the integrated dataset, the largest fold change (positive or negative) is the one reported.

Results

Overview of sequencing data outputs

Table 1 illustrates the number of sequence reads obtained for each timepoint and cell line combination alongside the number of reads that mapped to the human transcriptome, the ChAdOx1 nCoV-19 genome or the region of HuAd5 present in HEK293 cells. For MRC-5 cells, there was a peak of reads mapping to the ChAdOx1 genome at 48 hpi but overall the proportion of reads does not vary markedly. For A549 cells, there was a slight decline over time but notably there were approximately fivefold fewer reads as a proportion of human mapped transcripts mapping to the ChAdOx1 genome compared to MRC-5 cells. As expected, the proportion of reads that mapped to the ChAdOx1 genome in HEK293 cells, which are permissive for viral genome replication, was significantly higher. We also detected over a thousand reads mapping to the E1 and pIX region in HEK293 cells; this region covers just over 4000 bp from HuAd5 that is integrated into the HEK293 genome [31]. We utilised our previously published ORF-centric analysis pipeline to further characterise the transcripts allowing a qualitative and quantitative understanding of gene expression as we have previously done for wild type isolates of both human adenovirus type 5 and SARS-CoV-2 [24, 25]. Briefly, this analysis pipeline uses the genome sequence to correct the error-prone transcripts generated by the nanopore device. The dRNAseq data is biased towards the 3′ polyadenylated end, so we further constrained our analysis to transcripts whose mapped start site was within 50 bp of viral promoters as described for our analysis of wild type HuAd5 and SARS-CoV-2. After correction, each transcript was scanned to determine the 5′-most ORF which is compared to the locations and sequences of known ORFs and other features supplied by the user. In this way, the pipeline can identify and count the presence of known features on transcripts and ultimately provide quantitative data on the abundance of transcripts that code for each ORF. Moreover, this analysis illustrates the structure of the dominant transcript that codes for each ORF as well as characterising the structure of minor transcripts that could also code for any given ORF (Additional file 4: Table s1, Additional file 5: Table S2 and Additional file 6: Table S3).

Table 1 Read counts for individual dRNAseq datasets

Comparison of ChAdOx1 nCoV-19 adenoviral gene expression levels in non-permissive and permissive cells

Strikingly, in MRC-5 cells, the ChAdOx1 vector backbone genes (i.e., excluding the S glycoprotein gene) were hardly transcribed at any timepoint, but in A549 cells, we saw a wide range of transcripts that code for almost every protein predicted to be expressed by the ChAdOx1 backbone (Fig. 1b, c and Table 2). The pattern of dominant transcripts seen in A549 cells was broadly typical of mammalian adenoviruses [32] and was similar to the expression pattern of the virus grown in permissive HEK293 cells (compare Fig. 1c, d). Looking quantitatively, Table 2 shows that in A549 cells a large number of different ChAdOx1 vector backbone genes were expressed to one degree or another with the adenovirus protein DBP being dominant. For MRC-5 cells, transcripts that code for DBP and E4 ORFs 2 and 3 could be detected, but only at 48 hpi. However in A549 cells, in addition to DBP/pTP coding E2 transcripts, transcripts that code for all the E4 ORFs except E4 ORF6 could be detected, which could explain the detectable levels of major late transcript proteins in A549 cells. Notably, ChAdOx1 vector backbone gene expression declined over the time course analysed.

Table 2 List of ORFs searched for and the transcript frequency for each ORF

Comparing the expression of ChAdOx1 gene expression in non-permissive vs permissive HEK293 cells illustrated how deletion of the E1 region affects ChAdOx1 gene expression in A549 cells where some vector backbone gene expression has occurred (Table 2). One key difference was the level of expression of the fibre (L5) transcript which was the dominant transcript in HEK293 cells (9.65% of transcripts mapped to ChAdOx1 nCoV19) after the S glycoprotein transcript but is expressed only at very low levels in A549 cells (0.46% of transcripts mapped to ChAdOx1 nCoV19 at its peak at 48 hpi, Table 2). Furthermore, transcripts that code for all the E4 region ORFs were readily detected, indicating that replacing the Y25 ORFs 4, 6 and 6/7 with their HuAd5 equivalents did not ablate the expression of the remaining E4 ORF transcripts (Fig. 1 and Table 2).

Analysis of the transcriptome of ChAdOx1 nCoV-19 in HEK293 cells also allowed the observation of transcription from the region of HuAd5 that is integrated into HEK293 cells (Fig. 1e). Notably, there is some expression of HuAd5 pIX from this region (Table 2) which was approximately 100-fold lower than the ChAdOx1-derived pIX protein. In addition, in the non-permissive cell lines, there were transcripts that code for protein IX, which appear to be derived from RNA pol II and have ignored the polyA signal after the S glycoprotein and instead used the pIX poly A signal and spliced out the coding region for S glycoprotein (Fig. 1b, c). In permissive HEK293 cells, the dominant transcript coding for pIX appeared to revert to the expected transcript for pIX with the transcription start site just upstream of the initiating AUG (Fig. 1d). Notably, the type of transcript seen in the non-permissive cells was still present in HEK293 cells, but it was no longer the dominant transcript coding for pIX (Additional file 6: Table S3).

Expression of SARS-CoV-2 S glycoprotein transcripts

In both permissive and non-permissive cell lines, the most abundant ChAdOx1 nCoV-19 virus vector transcript, as expected, was the transcript for expression of the SARS-CoV-2 S glycoprotein. For MRC5 cells, the S transcript accounts for over 90% of transcripts from the ChAdOx1 nCoV19 genome; this is reduced to between 75 and 90% in A549 cells. However, in HEK293 cells (where there is productive viral replication), only 17% of transcripts code for the S glycoprotein relative to transcripts that code for adenovirus ORFs. This reflects the large numbers of adenovirus transcripts generated during productive replication. Notably, the S glycoprotein transcript is still the dominant transcript even in HEK293 cells. However, a number of transcripts were identified in all cell lines at low levels which appear to arise from aberrant splicing events. We also found evidence that occasionally some transcripts were spliced to other elements (notably pIX); these would still code for SARS-CoV-2 S glycoprotein but would be polycistronic (Fig. 2). We also saw limited evidence of transcripts that extended beyond the pIX polyA signal and deep into the ChAdOx1 vector backbone in HEK293 cells when there was active viral replication. However, it is important to place these rare transcripts in context—the vast majority of transcripts starting at the transcription start site for the S glycoprotein message would generate mRNA that codes for S glycoprotein (Table 3).

Fig. 2
figure 2

Transcript heterogeneity. This figure illustrates the structure of various transcripts observed at very low levels (less than 10 transcripts per dataset) all of which initiate at the transcription start site that drives expression of the S glycoprotein. For orientation, the structure of the dominant S glycoprotein transcript is shown in red in the examples selected from the MRC5 dataset. Three sets of example transcripts are shown from MRC5, A549 and HEK293 cells. The transcripts shown arise from splicing events that utilise only the canonical GU-AG splice donor/acceptor pair present in over 90% of eukaryotic splicing events. As can be seen, some of the transcripts arise from unintended splice site usage as well as from missing the intended BGH polyadenylation signal and utlilising other polyadenylation signals further down the vector genome

Table 3 Counts of transcripts arising from the expected transcription start site for S glycoprotein transcripts and their assignment

Additional polyadenylation site for ChAdOx1 preVI transcripts

Typically, in human adenoviruses, the late transcripts code for most of the structural viral proteins. The late transcripts typically originate at the major late promoter and include three exons (known as the tripartite leader) and optional i-leader exon and are grouped into five classes called L1 to L5 based on which of the five major polyadenylation sites they use [32, 33]. The transcript for preVI is usually regarded as part of the L3 group of transcripts which share a common polyadenylation site but utilise different splice acceptor/donor pairs to place one of three different ORFs (preVI, hexon or 23K) proximal to the 5′ cap for translation. In the case of ChAdOx1 replicating in HEK293 cells, although there were transcripts coding for preVI that would fit this pattern (e.g. utilising the same polyadenylation site as transcripts for hexon and 23K), the dominant transcript for preVI instead utilises an additional polyadenylation site upstream of the start codon for the hexon protein (Fig. 3a). This was also the case in A549 cells (Fig. 1b) where 12 transcripts utilised this additional polyadenylation site compared to just 4 using the classical L3 polyadenylation site (Additional file 5: Table S2). Examining the sequence of the viral genome in the region of this additional polyadenylation site reveals a GU-rich region preceded by a classical polyadenylation signal [34] (Fig. 3b)—this signal is not present in the equivalent region of HuAd5 (Fig. 3c) where no such preVI-specific polyadenylation site has been observed [24]. However, this additional polyadenylation site is present in human adenovirus type 4 (HuAd4, Fig. 3c) which is also a member of the group E mammalian adenoviruses like chimpanzee adenovirus Y25 from which ChaAdOx1 was derived.

Fig. 3
figure 3

Novel polyadenylation site usage for preVI. a The transcript structure of transcripts coding for three proteins classically considered to be L3 proteins viz preVI, hexon and 23K proteinase. The three vertical boxes on the left represent the tripartite exons, labelled TPL1, TPL2 and TPL3, and normally included in all major late transcripts generated by mammalian adenoviruses. The L3 transcripts would classically share the same L3 polyadenylation site (indicated on the diagram). In this diagram, the dominant transcripts seen in ChAdOx1 nCoV-19-infected 293 cells for each ORF are coloured red. This figure includes, in black, the transcript structure for observed transcripts that would both code for preVI and fit the canonical transcript structure for an L3 transcript, utilising the L3 polyadenylation site. In the case of 293 cells, however, there are only 192 copies of this transcript compared to 403 copies of the novel preVI transcript indicated in red at the bottom of a. b The sequences at the proposed polyadenylation site on ChAdOx1 showing the location of the polyadenylation signal and the GU-rich region that is often present downstream of a polyadenylation signal. c The equivalent regions in the genome of human adenoviruses 5 (a species C adenovirus) and human adenovirus type 4 which is, like ChAdOx1, a species E adenovirus. Taken together, b, c and d illustrate the similarities between the two species E adenovirus genomes and their distinction from a species C adenovirus genome

Proteomic detection of SARS-CoV-2 S glycoprotein and ChAdOx1 proteins

As expected, we were able to detect a range of S glycoprotein peptides in MRC5 and A549 cells as well as one phosphorylation site (Ser 1292) derived from the C-terminal portion of the glycoprotein which is internal to the viral particle (Table 4 and Additional files 7, 8, 9, 10, 11, 12, 13, 14: Tables S4, S5, S6, S7, S8, S9, S10, S11). For the vector backbone expressed proteins, we were only able to detect E4ORF3 in MRC5-infected cells whereas in A549 cells in addition to E4ORF3 we also detected DBP and hexon (Table 4). Interestingly in both MRC5 and A549 cells, we detected a peptide (SYLTPGDSSSGWTAGAAAY, aa248-265) apparently derived from a single transcript designated mRNA#962. However, this transcript has only one copy present in the transcriptome of the A549-infected cells and a similar transcript was not detected in the MRC-5 cells. Alternatively, we believe it is more likely that the detected peptide resulted from proteolytic activity at aa265 on full-length S glycoprotein either from intracellular processes or during sample harvest and preparation. In support of this, an analysis of the proteomics data to include semi tryptic peptides did identify the additional peptide (i.e. on the C-terminal side of an expected cleavage at aa265) which would be consistent with proteolytic activity at this site. The list of possible proteins predicted by translation of the 5′-most ORF of all the observed transcripts is very large and varied, and they were all used in the initial searches for peptides. However, this peptide was the only one that could not be predicted by the standard list of expected ChAdOx1 proteins. The proteomic analysis supported the transcriptomic analysis suggesting that the ChAdOx1 nCoV-19 vaccine does not make additional unexpected proteins. Notably, in both A549 and MRC5 cells, we observed some of the highest fold increases over time for the S glycoprotein as expected (Additional file 7: Table S4 and Additional file 9: Table S6).

Table 4 Peptides identified as derived from ChAdOx1 nCoV-19 in each cell line used

Analysis of the cellular total transcriptome, proteome and phosphoproteome

Whilst our primary focus was on the transcriptome, proteome and phosphoproteome of the ChAdOx1 nCoV-19 vaccine, we also collected data on the host cell counterparts which are summarised in Additional file 15: Table S12 and Additional file 16: Table S13. The two cell lines appear to have distinct patterns of response to the vaccine virus. We collated lists of proteins that were increased or decreased in abundance at least twofold at each time point and utilised STRING pathway analysis to determine if cellular pathways were apparently overrepresented (Additional file 15: Table S12 and Additional file 16: Table S13). Whilst the host cell responses to the vaccine are varied, our STRING analysis suggested that in MRC-5 cells there was an apparent enrichment for upregulated proteins involved in the unfolded protein response (e.g. HSPA6 and SLC5A3) and ER stress (e.g. HSPA5, DNAJB9 and HERPUD1), which was not the case for A549 cells. In A549 cells, there was some evidence that ChAdOx1 nCoV-19 infection affected proteins involved in ribosome biogenesis and host mRNA splicing (lists of proteins analysed are detailed in Additional file 15: Table S12 and Additional file 16: Table S13).

Discussion

Recombinant adenovirus vaccines have been developed over many years and have shown great promise as safe and effective vaccine platforms for mass vaccination programmes. Despite their widespread development, this is the first study to directly and comprehensively survey the transcriptomic repertoire of a replication-defective adenovirus in a non-permissive host cell. A significant advantage of using simian-based adenoviruses like ChAdOx1 (a species E non-human adenovirus) to study the transcriptome of E1-deleted adenoviruses is the absence of replication-competent adenoviruses (RCAs) [35]. RCAs can arise from homology between the HuAd5 (a species C human adenovirus) sequences present in the HEK293 cells used to produce the vaccine vector which extend beyond the E1 region into the pIX gene (as shown in Fig. 1e) and sequences present in the recombinant vector itself (e.g. pIX gene in the vector), which are frequently based on HuAd5. Whilst the emergence of an RCA could confound such an analysis in other adenoviral backgrounds, we have never detected RCA in any preparation of ChAdOx1-based vaccine vectors. This is presumably due to insufficient homology between the species E chimpanzee adenoviruses and the sequences in HEK293 cells [36]. Given the significant advantages in using simian-based vectors as a human vaccine platform [37], we wanted to deepen our understanding of these virus vectors on a molecular level. We were especially keen to determine that no additional unanticipated transcripts or proteins were being made. Any one of such proteins could be antigenic with unintended consequences (e.g. generating auto-immune responses for example).

The two host cell lines chosen for this study have distinct properties despite both being derived from male human lungs. The A549 cell line is an immortal cancerous cell line with an average of 66 chromosomes and a deletion in the CDKN2A locus leaving the cell line defective in the p53/ARFp14/MDM2 pathway. By contrast, MRC5 cells are genetically normal with a finite cell passage capability. Thus, the distinct transcriptional profiles seen from the ChAdOx1 nCoV-19 vaccine could result from the A549 cells being genetically defective and immortal or it may be connected to unrelated differences in the intracellular environment. Previous work with human adenovirus-based vectors has suggested that a cellular E1A-like activity is in some manner connected to pre-existing levels of heat shock proteins, with promyelocytic leukaemia protein (PML) also being implicated [38, 39]. Notably, in the mock transcriptomic data, there are higher levels of sequenced HSPA1A and HSPA1B transcripts in A549 cells compared to MRC-5 cells prior to infection with ChAdOx1 nCoV-19 with the reverse being true for PML transcripts (Additional file 15: Table S12 and Additional file 16: Table S13). We noted that in MRC5 cells, levels of transcripts for HSPA1A and HSPA1B show increases over time and for HSPA1B there is a corresponding protein increase (Additional file 15: Table S12 and Additional file 16: Table S13), whereas in A549 cells, there is a decrease in sequenced transcripts for HSPA1A and HSPA1B with no significant change in protein levels for HSPA1A. This supports previous research which implies heat shock (and thus increased expression of heat shock proteins) does not overcome the lack of E1A expression [39] as in MRC5 cells we do not see increasing amounts of ChAdOx1 vector backbone transcription. In addition, NFκB has been implicated in enabling HuAd5-based vector backbone expression in non-permissive cell lines [40]. However, no marked differences were observed between the two cell lines in the transcriptomic abundance data for NFKB1 in our datasets (Additional file 15: Table S12 and Additional file 16: Table S13). A deeper understanding of the intracellular environments that allow the ChAdOx1 vector to overcome the lack of E1A, even in a limited fashion, could lead to better ways to prevent this which in turn may lead to improved transgene expression. The transcriptome and proteome of the two cell lines appeared to respond differently to the vaccine vector. However, whether these differences are driven by or a result of the distinct adenovirus vector expression patterns observed on the vector genome is yet to be determined.

That we observe quite distinct vaccine backbone expression in cell lines from the same tissue site is notable. However, care is needed not to overinterpret this study which was designed to investigate the potential repertoire of transcription from a respiratory virus-based vaccine vector in two different human respiratory cell lines. This vaccine vector is currently administered intramuscularly (other routes of administration including respiratory ones may also prove effective longer term), and so in further analysis, it might be useful to examine primary cells and/or cell lines derived from muscles or even biopsy material. In addition, there is clearly the potential for other cell types to become infected after intramuscular administration that may also be relevant to the adaptive immune response.

We have previously analysed the proteome of purified adenovirus particles, both wild type HuAd5- and HuAd5-based recombinant vaccine vectors, including a sample of adenovirus manufactured to clinical grade. We were able to detect non-structural proteins DBP, 100K and E4 14.7K protein in purified virus particles in addition to the expected viral structural proteins [41]. If any cells in a vaccinee did express the full range of adenovirus proteins over time similar to the A549 transcriptomic profile, then T-cell responses to the virus vector could both derive from the incoming viral proteins and from subsequent low-level expression of any of the remaining virus vector backbone genes. That some cells allow low-level expression of ChAdOx1 vector backbone transcripts would be consistent with data from human adenovirus vector studies where additional deletions in E2 and E4 or helper-dependent adenovirus vectors (where essentially all the vector backbone genes are removed) were shown to afford longer transgene expression in vivo [13, 15, 42,43,44]. The continuing low-level expression of vector backbone genes is therefore likely to be the main driver of immune-mediated clearance of cells infected with E1/E3-deleted adenoviruses. The repercussions in a vaccinee of broad low-level expression as observed in A549 cells is difficult to evaluate at this stage, but in principle, it suggests that repeated re-administration of a particular vaccine vector beyond the usual prime-boost vaccination regimen should be examined in further detail. In addition, low-level expression of vector backbone genes in A549 cells may be related to our finding that A549 cells express lower levels of S glycoprotein mRNA as a proportion of total mRNA than MRC5 cells (Table 1).

Our proteomics analysis is consistent with the transcriptomics data in that the S glycoprotein is readily detected by MS/MS and that we detect a slightly wider array of ChAdOx1 proteins in A549 cells compared to MRC-5 cells. We also detect evidence of phosphorylation of the S glycoprotein as we have noted previously in our work on SARS-CoV-2 [25]; the significance of this phosphorylation is unclear at this time. We do not detect any evidence of additional unexpected proteins being expressed by the ChAdOx1 nCoV19 vaccine despite a broad search of proteins coded for by all transcripts that map to the ChAdOx1 nCoV-19 genome.

We previously analysed the transcriptomic repertoire of wild type HuAd5 in MRC-5 cells showing that during a productive infection there is significant low-level heterogeneity in the usage of both splice sites and polyadenylation signals [24]. Here, we show that as with the human virus there is similar low-level expression of transcripts with aberrant splice site and polyadenylation signal usage in cell lines infected with the recombinant ChAdOx1 nCoV-19 virus in replication permissive and non-permissive cell lines. Critically, the overwhelming majority of S glycoprotein transcripts have the expected structure with only a small minority of transcripts originating from the S glycoprotein transcription start site being unable to express the S glycoprotein. The most common issue seems to be failure to use the polyadenylation signal immediately after the S glycoprotein ORF. This leads to, for example, usage of the pIX polyadenylation site and the opportunity for splicing events that may disrupt S glycoprotein expression. Other S transcripts with aberrant splice patterns or polyadenylation site usage were detected, but these were minor compared to the dominant transcript and we did not observe peptides that could arise from such aberrant transcripts. In future vector designs, this approach should highlight such issues potentially allowing changes to vector design to minimise issues identified.

Notably, in our analysis of wild type HuAd5 transcripts, we proposed that a key aspect of HuAd5 evolution would be the exploration of different splice acceptor/donor sites and initiating codon combinations. Here, analysis of expression from the ChAdOx1 backbone identified alternative usage of polyadenylation signals. Typically, the polyadenylation signal for the preVI protein transcript would be shared with the transcripts for hexon and 23K—collectively known as the L3 group since they share this polyadenylation signal. Protein preVI uses two polyadenylation signals, one that is the L3 location used by the hexon and 23K transcripts and a second just after its own stop codon. That this second one is not observed in human adenovirus type 5 is likely due to the lack of a strong polyadenylation signal as shown in Fig. 3. However, this arrangement seems to be present both in HuAd4 and chimpanzee adenovirus Y25 (which the ChAdOx1 vector is based on) which are both species E adenoviruses. Indeed, HuAd4 is closely enough related to chimpanzee adenoviruses to suggest that the human virus may have jumped into chimpanzees at some point in the past [45]. What advantage, if any, the two polyadenylation sites offer the ChAdOx1 virus is not clear, but the use of the polyadenylation site after the preVI ORF would by necessity reduce the influence of splice site selection in determining whether preVI, hexon or 23K would be expressed from transcripts utilising the canonical L3 polyadenylation site. This could, in principle, bias gene expression in favour of preVI relative to hexon or 23K. Whether this additional polyadenylation site is being selected for or against remains an open question.

Examination of the HEK293 cell data revealed, as expected, expression of HuAd5 transcripts corresponding to E1 region genes dominated by E1b19K, E1a12S and E1a13S as we saw for wild type HuAd5 gene expression at 16hpi infection [24]. As noted in the results, in HEK293 cells, we observed a small number of HuAd5 pIX transcripts, approximately 100-fold fewer than ChadOx1 pIX transcripts. However, with 240 copies of protein IX per virus capsid [46,47,48], it is possible that the ChAdOx1 virus particles grown in HEK293 cells (or derivatives) contain one or two copies of HuAd5-pIX per viral particle assuming the two proteins are equivalently interchangeable on the virus particle.

Conclusions

Our analyses provide valuable insight into the transcriptomic and proteomic repertoire of an important vaccine, providing confirmation that the vaccine vector’s transcriptome is essentially as intended in these cell lines. However, care should be taken when extrapolating the results of the transcriptomic analysis using single cell types to the complex multicellular differentiated cell-rich environment these types of vaccines encounter when they are administered (by whatever route). Recently, an in-depth multi-omic analysis of the herpes virus genome revealed that an oncolytic herpes virus licenced in 2015 under the name Imlygic is in fact deleted for an additional third gene rather than the two intended because of a previously unknown ORF present in the deleted region [49]. Whilst there is no suggestion that this has been a problem, it highlights the importance of utilising state-of-the-art and unbiased approaches to survey genetically modified viruses intended for clinical use. Finally, we argue that this kind of analysis is relatively straightforward and should be routinely incorporated into the early stages of future viral vector evaluation pipelines to allow a robust understanding of the transcriptomic potential of engineered viral vectors.

Availability of data and materials

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository [50] with the dataset identifier PXD023286, https://www.ebi.ac.uk/pride/archive/projects/PXD023286 [51].

All the fastq files are available from the Short Read Archive with the identifier PRJNA688412, https://www.ncbi.nlm.nih.gov/sra/PRJNA688412 [52].

References

  1. Folegatti PM, Ewer KJ, Aley PK, Angus B, Becker S, Belij-Rammerstorfer S, Bellamy D, Bibi S, Bittaye M, Clutterbuck EA, et al. Safety and immunogenicity of the ChAdOx1 nCoV-19 vaccine against SARS-CoV-2: a preliminary report of a phase 1/2, single-blind, randomised controlled trial. Lancet. 2020;396:467–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Logunov DY, Dolzhikova IV, Zubkova OV, Tukhvatullin AI, Shcheblyakov DV, Dzharullaeva AS, Grousova DM, Erokhova AS, Kovyrshina AV, Botikov AG, et al. Safety and immunogenicity of an rAd26 and rAd5 vector-based heterologous prime-boost COVID-19 vaccine in two formulations: two open, non-randomised phase 1/2 studies from Russia. Lancet. 2020;396(10255):887–97.

  3. Mercado NB, Zahn R, Wegmann F, Loos C, Chandrashekar A, Yu J, Liu J, Peter L, McMahan K, Tostanoski LH, et al. Single-shot Ad26 vaccine protects against SARS-CoV-2 in rhesus macaques. Nature. 2020.

  4. Zhu FC, Guan XH, Li YH, Huang JY, Jiang T, Hou LH, Li JX, Yang BF, Wang L, Wang WJ, et al. Immunogenicity and safety of a recombinant adenovirus type-5-vectored COVID-19 vaccine in healthy adults aged 18 years or older: a randomised, double-blind, placebo-controlled, phase 2 trial. Lancet. 2020;396:479–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Hassan AO, Kafai NM, Dmitriev IP, Fox JM, Smith BK, Harvey IB, Chen RE, Winkler ES, Wessel AW, Case JB, et al. A single-dose intranasal ChAd vaccine protects upper and lower respiratory tracts against SARS-CoV-2. Cell. 2020.

  6. Davison AJ, Benko M, Harrach B. Genetic content and evolution of adenoviruses. J Gen Virol. 2003;84:2895–908.

    Article  CAS  PubMed  Google Scholar 

  7. Wold WS, Toth K. Adenovirus vectors for gene therapy, vaccination and cancer gene therapy. Curr Gene Ther. 2013;13:421–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Graham FL, Smiley J, Russell WC, Nairn R. Characteristics of a human cell line transformed by DNA from human adenovirus type 5. J Gen Virol. 1977;36:59–74.

    Article  CAS  PubMed  Google Scholar 

  9. Saha B, Parks RJ. Human adenovirus type 5 vectors deleted of early region 1 (E1) undergo limited expression of early replicative E2 proteins and DNA replication in non-permissive cells. PLoS One. 2017;12:e0181012.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Gorziglia MI, Kadan MJ, Yei S, Lim J, Lee GM, Luthra R, Trapnell BC. Elimination of both E1 and E2 from adenovirus vectors further improves prospects for in vivo human gene therapy. J Virol. 1996;70:4173–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Rittner K, Schultz H, Pavirani A, Mehtali M. Conditional repression of the E2 transcription unit in E1-E3-deleted adenovirus vectors is correlated with a strong reduction in viral DNA replication and late gene expression in vitro. J Virol. 1997;71:3307–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Shimizu K, Sakurai F, Machitani M, Katayama K, Mizuguchi H. Quantitative analysis of the leaky expression of adenovirus genes in cells transduced with a replication-incompetent adenovirus vector. Mol Pharm. 2011;8:1430–5.

    Article  CAS  PubMed  Google Scholar 

  13. Yang Y, Nunes FA, Berencsi K, Furth EE, Gonczol E, Wilson JM. Cellular immunity to viral antigens limits E1-deleted adenoviruses for gene therapy. Proc Natl Acad Sci U S A. 1994;91:4407–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Nakai M, Komiya K, Murata M, Kimura T, Kanaoka M, Kanegae Y, Saito I. Expression of pIX gene induced by transgene promoter: possible cause of host immune response in first-generation adenoviral vectors. Hum Gene Ther. 2007;18:925–36.

    Article  CAS  PubMed  Google Scholar 

  15. Morsy MA, Gu M, Motzel S, Zhao J, Lin J, Su Q, Allen H, Franlin L, Parks RJ, Graham FL, et al. An adenoviral vector deleted for all viral coding sequences results in enhanced safety and extended expression of a leptin transgene. Proc Natl Acad Sci U S A. 1998;95:7866–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Bangari DS, Mittal SK. Development of nonhuman adenoviruses as vaccine vectors. Vaccine. 2006;24:849–62.

    Article  CAS  PubMed  Google Scholar 

  17. Dicks MD, Spencer AJ, Edwards NJ, Wadell G, Bojang K, Gilbert SC, Hill AV, Cottingham MG. A novel chimpanzee adenovirus vector with low human seroprevalence: improved systems for vector derivation and comparative immunogenicity. PLoS One. 2012;7:e40385.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. van Doremalen N, Haddock E, Feldmann F, Meade-White K, Bushmaker T, Fischer RJ, Okumura A, Hanley PW, Saturday G, Edwards NJ, et al. A single dose of ChAdOx1 MERS provides protective immunity in rhesus macaques. Sci Adv. 2020;6:eaba8399.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Stedman A, Wright D, Wichgers Schreur PJ, Clark MHA, Hill AVS, Gilbert SC, Francis MJ, van Keulen L, Kortekaas J, Charleston B, Warimwe GM. Safety and efficacy of ChAdOx1 RVF vaccine against Rift Valley fever in pregnant sheep and goats. NPJ Vaccines. 2019;4:44.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. López-Camacho C, Abbink P, Larocca RA, Dejnirattisai W, Boyd M, Badamchi-Zadeh A, Wallace ZR, Doig J, Velazquez RS, Neto RDL, et al. Rational Zika vaccine design via the modulation of antigen membrane anchors in chimpanzee adenoviral vectors. Nat Commun. 2018;9:2441.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Coughlan L, Sridhar S, Payne R, Edmans M, Milicic A, Venkatraman N, Lugonja B, Clifton L, Qi C, Folegatti PM, et al. Heterologous two-dose vaccination with simian adenovirus and poxvirus vectors elicits long-lasting cellular immunity to influenza virus A in healthy adults. EBioMedicine. 2018;29:146–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Antrobus RD, Coughlan L, Berthoud TK, Dicks MD, Hill AV, Lambe T, Gilbert SC. Clinical assessment of a novel recombinant simian adenovirus ChAdOx1 as a vectored vaccine expressing conserved influenza A antigens. Mol Ther. 2014;22:668–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Voysey M, Clemens SAC, Madhi SA, Weckx LY, Folegatti PM, Aley PK, Angus B, Baillie VL, Barnabas SL, Bhorat QE, et al. Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: an interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK. Lancet. 2020.

  24. Donovan-Banfield I, Turnell AS, Hiscox JA, Leppard KN, Matthews DA. Deep splicing plasticity of the human adenovirus type 5 transcriptome drives virus evolution. Commun Biol. 2020;3:124.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Davidson AD, Williamson MK, Lewis S, Shoemark D, Carroll MW, Heesom KJ, Zambon M, Ellis J, Lewis PA, Hiscox JA, Matthews DA. Characterisation of the transcriptome and proteome of SARS-CoV-2 reveals a cell passage induced in-frame deletion of the furin-like cleavage site from the spike glycoprotein. Genome Med. 2020;12:68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Wang Y, Grunewald M, Perlman S. Coronaviruses: an updated overview of their replication and pathogenesis. Methods Mol Biol. 2020;2203:1–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Matthews DA, Cummings D, Evelegh C, Graham FL, Prevec L. Development and use of a 293 cell line expressing lac repressor for the rescue of recombinant adenoviruses expressing high levels of rabies virus glycoprotein. J Gen Virol. 1999;80(Pt 2):345–53.

    Article  CAS  PubMed  Google Scholar 

  28. Chapman BS, Thayer RM, Vincent KA, Haigwood NL. Effect of intron A from human cytomegalovirus (Towne) immediate-early gene on heterologous expression in mammalian cells. Nucleic Acids Res. 1991;19:3979–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34:3094–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Soneson C, Yao Y, Bratus-Neuenschwander A, Patrignani A, Robinson MD, Hussain S. A comprehensive examination of Nanopore native RNA sequencing for characterization of complex transcriptomes. Nat Commun. 2019;10:3359.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Louis N, Evelegh C, Graham FL. Cloning and sequencing of the cellular-viral junctions from the human adenovirus type 5 transformed 293 cell line. Virology. 1997;233:423–9.

    Article  CAS  PubMed  Google Scholar 

  32. Akusjarvi G. Temporal regulation of adenovirus major late alternative RNA splicing. Front Biosci. 2008;13:5006–15.

    Article  CAS  PubMed  Google Scholar 

  33. Chow LT, Broker TR. The spliced structures of adenovirus 2 fiber message and the other late mRNAs. Cell. 1978;15:497–510.

    Article  CAS  PubMed  Google Scholar 

  34. Beaudoing E, Freier S, Wyatt JR, Claverie JM, Gautheret D. Patterns of variant polyadenylation signal usage in human genes. Genome Res. 2000;10:1001–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Lochmuller H, Jani A, Huard J, Prescott S, Simoneau M, Massie B, Karpati G, Acsadi G. Emergence of early region 1-containing replication-competent adenovirus in stocks of replication-defective adenovirus recombinants (delta E1 + delta E3) during multiple passages in 293 cells. Hum Gene Ther. 1994;5:1485–91.

    Article  CAS  PubMed  Google Scholar 

  36. Tatsis N, Tesema L, Robinson ER, Giles-Davis W, McCoy K, Gao GP, Wilson JM, Ertl HC. Chimpanzee-origin adenovirus vectors as vaccine carriers. Gene Ther. 2006;13:421–9.

    Article  CAS  PubMed  Google Scholar 

  37. Morris SJ, Sebastian S, Spencer AJ, Gilbert SC. Simian adenoviruses as vaccine vectors. Future Virol. 2016;11:649–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Atwan Z, Wright J, Woodman A, Leppard KN. Promyelocytic leukemia protein isoform II inhibits infection by human adenovirus type 5 through effects on HSP70 and the interferon response. J Gen Virol. 2016;97:1955–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Imperiale MJ, Kao HT, Feldman LT, Nevins JR, Strickland S. Common control of the heat shock gene and early adenovirus genes: evidence for a cellular E1A-like activity. Mol Cell Biol. 1984;4:867–74.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Machitani M, Sakurai F, Wakabayashi K, Nakatani K, Shimizu K, Tachibana M, Mizuguchi H. NF-κB promotes leaky expression of adenovirus genes in a replication-incompetent adenovirus vector. Sci Rep. 2016;6:19922.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Alqahtani A, Heesom K, Bramson JL, Curiel D, Ugai H, Matthews DA. Analysis of purified wild type and mutant adenovirus particles by SILAC based quantitative proteomics. J Gen Virol. 2014;95:2504–11.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Rhee EG, Blattman JN, Kasturi SP, Kelley RP, Kaufman DR, Lynch DM, La Porte A, Simmons NL, Clark SL, Pulendran B, et al. Multiple innate immune pathways contribute to the immunogenicity of recombinant adenovirus vaccine vectors. J Virol. 2011;85:315–23.

    Article  CAS  PubMed  Google Scholar 

  43. Coughlan L. Factors which contribute to the immunogenicity of non-replicating adenoviral vectored vaccines. Front Immunol. 2020;11:909. https://www.frontiersin.org/article/10.3389/fimmu.2020.00909.

  44. Lusky M, Christ M, Rittner K, Dieterle A, Dreyer D, Mourot B, Schultz H, Stoeckel F, Pavirani A, Mehtali M. In vitro and in vivo biology of recombinant adenovirus vectors with E1, E1/E2A, or E1/E4 deleted. J Virol. 1998;72:2022–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Purkayastha A, Ditty SE, Su J, McGraw J, Hadfield TL, Tibbetts C, Seto D. Genomic and bioinformatics analysis of HAdV-4, a human adenovirus causing acute respiratory disease: implications for gene therapy and vaccine vector development. J Virol. 2005;79:2559–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. van Oostrum J, Burnett RM. Molecular composition of the adenovirus type 2 virion. J Virol. 1985;56:439–48.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Reddy VS, Nemerow GR. Structures and organization of adenovirus cement proteins provide insights into the role of capsid maturation in virus entry and infection. Proc Natl Acad Sci U S A. 2014;111:11715–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kundhavai Natchiar S, Venkataraman S, Mullen TM, Nemerow GR, Reddy VS. Revised crystal structure of human adenovirus reveals the limits on protein IX quasi-equivalence and on analyzing large macromolecular complexes. J Mol Biol. 2018;430:4132–41.

    Article  CAS  PubMed  Google Scholar 

  49. Whisnant AW, Jürges CS, Hennig T, Wyler E, Prusty B, Rutkowski AJ, L’Hernault A, Djakovic L, Göbel M, Döring K, et al. Integrative functional genomics decodes herpes simplex virus 1. Nat Commun. 2020;11:2038.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, Inuganti A, Griss J, Mayer G, Eisenacher M, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47:D442–50.

    Article  CAS  PubMed  Google Scholar 

  51. Almuqrin A, Davidson A D., Williamson MK, Lewis P, Heesom K, Morris S, Gilbert SC, Matthews DA: ChAdOx infection of human cell lines A549 MRC-5. Project PXD023286. PRIDE database. https://www.ebi.ac.uk/pride/archive/projects/PXD023286. 2021.

  52. Almuqrin A, Davidson A D., Williamson MK, Lewis P, Heesom K, Morris S, Gilbert SC, Matthews DA: ChAdOx1 nCoV19 infection of human cell lines. BioProject PRJNA688412. NCBI Sequence Read Archive https://www.ncbi.nlm.nih.gov/sra/PRJNA688412. 2021.

Download references

Acknowledgements

We would like to thank Prof A. Finn for valuable discussions prior to undertaking the work.

Funding

D.A.M. and A.D.D. were supported by the BBSRC (grant BB/M02542X/1) and by the US Food and Drug Administration grant number HHSF223201510104C “Ebola Virus Disease: correlates of protection, determinants of outcome and clinical management” amended to incorporate urgent COVID-19 studies. SJG and SM were funded through VaxHub from the Engineering and Physical Sciences Research Council (EP/R013756/1).

Author information

Authors and Affiliations

Authors

Contributions

DAM, AA, SM and SJG led the conception and design of the work and led the manuscript drafting. DAM, MKW and ADD performed infections and harvesting of ChAdOx1 nCoV-19-infected cells. AA extracted and sequenced cellular mRNA. KH led proteomics and phosphoproteomics data acquisition and analysis with PAL providing additional analysis and proteomics dataset integration. DAM integrated proteomics and transcriptomics datasets. All authors contributed to manuscript drafts and figure preparation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to David A. Matthews.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

SG is the co-founder of Vaccitech (co-inventors of this vaccine) and named as an inventor on a patent covering the use of ChAdOx1-vectored vaccines and a patent application covering this SARS-CoV-2 vaccine. The remaining authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Features file. A text file delimiting the locations of known features on the ChAdOx1 nCoV-19 genome.

Additional file 2.

Known human transcripts. A fasta file of human transcripts from the EBI data repository and used for mapping of sequence data to the human transcriptome.

Additional file 3.

Proteins search list. The protein search list used in the proteomics analysis.

Additional file 4: Table S1.

MRC5 transcript groups. An excel spreadsheet describing the structure of all ChAdOx1 nCoV19 derived transcript groups identified and how many transcripts belong to each group. In addition, it also describes what features are present on each transcript group. This list is the combined data from all time points in MRC5 cells.

Additional file 5: Table S2.

A549 transcript groups. An excel spreadsheet describing the structure of all ChAdOx1 nCoV19 derived transcript groups identified and how many transcripts belong to each group. In addition, it also describes what features are present on each transcript group. This list is the combined data from all time points in A549 cells.

Additional file 6: Table S3.

HEK293 transcript groups. An excel spreadsheet describing the structure of all ChAdOx1 nCoV19 derived transcript groups identified and how many transcripts belong to each group. In addition, it also describes what features are present on each transcript group. This list is the combined data from all time points in 293 cells.

Additional file 7: Table S4.

A549 proteins identified. List of all proteins identified in A549 cells or A549 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 8: Table S5.

A549 phosphoproteins identified. List of all phosphoproteins identified in A549 cells or A549 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 9: Table S6.

MRC5 proteins identified. List of all proteins identified in MRC5 cells or MRC5 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 10: Table S7.

MRC5 phosphoproteins identified. List of all phosphoproteins identified in MRC5 cells or MRC5 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 11: Table S8.

A549 Peptide Spectral Matches. List of all peptide spectral matches (PSMs) identified in A549 cells or A549 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 12: Table S9.

A549 phospho Peptide Spectral Matches. List of all phospho peptide spectral matches (PSMs) identified in A549 cells or A549 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 13: Table S10.

MRC5 Peptide Spectral Matches. List of all peptide spectral matches (PSMs) identified in MRC5 cells or MRC5 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 14: Table S11.

MRC5 phospho Peptide Spectral Matches. List of all phospho peptide spectral matches (PSMs) identified in MRC5 cells or MRC5 cells infected for various times with ChAdOx1 nCoV-19.

Additional file 15: Table S12.

Integrated transcriptomics and proteomics table, A549 cells. Integrated list of transcriptomic, proteomic and phosphoproteomic changes over time in A549 cells.

Additional file 16: Table S13.

Integrated transcriptomics and proteomics table, MRC5 cells. Integrated list of transcriptomic, proteomic and phosphoproteomic changes over time in MRC5 cells.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Almuqrin, A., Davidson, A.D., Williamson, M.K. et al. SARS-CoV-2 vaccine ChAdOx1 nCoV-19 infection of human cell lines reveals low levels of viral backbone gene transcription alongside very high levels of SARS-CoV-2 S glycoprotein gene transcription. Genome Med 13, 43 (2021). https://doi.org/10.1186/s13073-021-00859-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13073-021-00859-1