Vaccine
The optimized DNA sequence encoding SARS-CoV-2 IgE-spike was created using Inovio’s proprietary in silico Gene Optimization Algorithm to enhance expression and immunogenicity [12]. The optimized DNA sequence was synthesized, digested with BamHI and XhoI, and cloned into the expression vector pGX0001 under the control of the human cytomegalovirus immediate-early promoter and a bovine growth hormone polyadenylation signal.
Animals
Eighteen rhesus macaques of Indian origin (Macaca mulatta) from a UK breeding colony were used in this study. Study groups comprised three males and three females of each species and all were adults aged between 2.5 and 3.5 years of age and weighing > 4 Kg at time of challenge. Prior to the start of the experiment, socially compatible animals were randomly assigned to challenge groups, to minimize bias. Animals were housed in compatible social groups, in cages in accordance with the UK Home Office Code of Practice for the Housing and Care of Animals Bred, Supplied or Used for Scientific Procedures (2014) and National Committee for Refinement, Reduction and Replacement (NC3Rs) Guidelines on Primate Accommodation, Care and Use, August 2006. Housing prior and for the duration of challenge is described in [14]. All experimental work was conducted under the authority of a UK Home Office approved project license that had been subject to local ethical review at PHE Porton Down by the Animal Welfare and Ethical Review Body (AWERB) and approved as required by the Home Office Animals (Scientific Procedures) Act 1986. Animals were sedated by intramuscular (IM) injection with ketamine hydrochloride (Ketaset, 100 mg mL-1, Fort Dodge Animal Health Ltd, Southampton, UK; 10mg kg-1) for procedures requiring removal from their housing. None of the animals had been used previously for experimental procedures.
Animals received 1 mg of SARS-CoV-2 DNA vaccine, INO-4800, by intradermal injection at day 28 only (1 dose group) or 0 and 28 (2 dose group) followed by an EP treatment using the CELLECTRA 2000® Adaptive Constant Current Electroporation Device with a 3P array (Inovio Pharmaceuticals).
Serum and heparinised whole blood were collected whilst animals were sedated at bi-weekly intervals during the vaccination phase. Nasal and throat swabs were also collected on the day of challenge on D56. After challenge, nasal swabs, throat swabs and serum were collected at 1, 3, 5 dpc and at cull (6, 7 or 8 dpc – staggered due to the high level of labour involved in procedures), with heparinised whole blood collected at 3 dpc and at cull. Nasal and throat swabs were obtained as described [14].
Clinical observations
Animals were monitored multiple times per day for behavioral and clinical changes. Behavior was evaluated for contra-indicators including depression, withdrawal from the group, aggression, changes in feeding patterns, breathing pattern, respiration rate and cough. Animals were observed and scored as follows for activity and health throughout the study. Key: Activity Level: A0 = Active & Alert; A1 = Only active when stimulated by operator; A2=Inactive even when stimulated/Immobile; H=Healthy; S=Sneeze, C=Cough, Nd=Nasal Discharge, Od=Ocular Discharge, Rn=Respiratory Noises, Lb=Laboured breathing, L=Lethargy, Di=Diarrhoea, Ax = Loss of Appetite, Dx = Dehydration, RD=Respiratory Distress. Animal body weight, temperature and haemoglobin levels were measured and recorded throughout the study.
Viruses and Cells
SARS-CoV-2 Victoria/01/2020[25] was generously provided by The Doherty Institute, Melbourne, Australia at P1 after primary growth in Vero/hSLAM cells and subsequently passaged twice at PHE Porton Down in Vero/hSLAM cells [ECACC 04091501]. Infection of cells was with ~0.0005 MOI of virus and harvested at day 4 by dissociation of the remaining attached cells by gentle rocking with sterile 5 mm borosilicate beads followed by clarification by centrifugation at 1,000 x g for 10 mins. Whole genome sequencing was performed, on the P3 challenge stock, using both Nanopore and Illumina as described previously [26]. Virus titre of the challenge stocks was determined by plaque assay on Vero/E6 cells [ECACC 85020206]. Cell lines were obtained from the European Collection of Authenticated Cell Cultures (ECACC) PHE, Porton Down, UK. Cell cultures were maintained at 37oC in Minimum essential medium (MEM) (Life Technologies, California, USA) supplemented with 10% fetal bovine serum (FBS) (Sigma, Dorset, UK) and 25 mM HEPES (Life Technologies, California, USA). In addition, Vero/hSLAM cultures were supplemented with 0.4 mg mL-1 of geneticin (Invitrogen) to maintain the expression plasmid. Challenge substance dilutions were conducted in phosphate buffer saline (PBS). Inoculum (5 x 106 PFU) was delivered by intratracheal route (2 ml) and intranasal instillation (1.0 ml total, 0.5 ml per nostril).
Clinical signs and in-life imaging by computerised tomography
CT scans were performed two weeks before and five days after challenge with SARS-CoV2. CT imaging was performed on sedated animals using a 16 slice Lightspeed CT scanner (General Electric Healthcare, Milwaukee, WI, USA) in both the prone and supine position and scans evaluated by a medical radiologist expert in respiratory diseases (as described previously [14]). To provide the power to discriminate differences between individual NHP’s with low disease volume (i.e. <25% lung involvement), a refined score system was designed in which scores were attributed for possession of abnormal features characteristic of COVID in human patients (COVID pattern score) and for the distribution of features through the lung (Zone score). The COVID pattern score was calculated as sum of scores assigned for the number of nodules identified, and the possession and extent of GGO and consolidation according to the following system: Nodule(s): Score 1 for 1, 2 for 2 or 3, 3 for 4 or more; GGO: each affected area was attributed with a score according to the following: Score 1 if area measured < 1 cm, 2 if 1 to 2 cm, 3 if 2 -3 cm, 4 if > 3 cm and scores for each area of GGO were summed to provide a total GGO score; Consolidation: each affected area was attributed with a score according to the following: 1 if area measured < 1 cm, 2 if 1 to 2 cm, 3 if 2 -3 cm, 4 if > 3 cm. Scores for each area of consolidation are summed to provide a total consolidation score. To account for estimated additional disease impact on the host of consolidation compared to GGO, the score system was weighted by doubling the score assigned for consolidation. To determine the zone score, the lung was divided into 12 zones and each side of the lung divided (from top to bottom) into three zones: the upper zone (above the carina), the middle zone (from the carina to the inferior pulmonary vein), and the lower zone (below the inferior pulmonary vein). Each zone was further divided into two areas: the anterior area (the area before the vertical line of the midpoint of the diaphragm in the sagittal position) and the posterior area (the area after the vertical line of the mid-point of the diaphragm in the sagittal position). This results in 12 zones in total where a score of one is attributed to each zone containing structural changes. The COVID pattern score and the zone are summed to provide the Total CT score.
Post-mortem examination and histopathology
Animals were euthanized at 3 different time-points, in groups of six (one male and one female from each group) at 6, 7 and 8 dpc. The bronchial alveolar lavage fluid (BAL) was collected at necropsy from the right lung. The left lung was dissected prior to BAL collection and used for subsequent histopathology and virology procedures. At necropsy nasal and throat swabs, heparinised whole blood and serum were taken alongside tissue samples for histopathology. Samples from the left cranial and left caudal lung lobe together with spleen, kidney, liver, mediastinal and axillary lymph nodes, small intestine (duodenum), large intestine (colon), trachea, larynx inoculation site and draining lymph node, were fixed by immersion in 10% neutral-buffered formalin and processed routinely into paraffin wax. Four µm sections were cut and stained with hematoxylin and eosin (H&E) and examined microscopically. A lung histopathology scoring system [14] was used to evaluate lesions affecting the airways and the parenchyma. Three tissue sections from each left lung lobe were used to evaluate the lung histopathology. In addition, samples were stained using the RNAscope technique to identify the SARS-CoV-2 virus RNA in lung tissue sections. Briefly, tissues were pre-treated with hydrogen peroxide for 10 mins (RT), target retrieval for 15 mins (98-102°C) and protease plus for 30 mins (40°C) (Advanced Cell Diagnostics). A V-nCoV2019-S probe (SARS-CoV-2 Spike gene specific) was incubated on the tissues for two hours at 40°C. In addition, samples were stained using the RNAscope technique to identify the SARS-CoV-2 virus RNA. Amplification of the signal was carried out following the RNAscope protocol using the RNAscope 2.5 HD Detection kit – Red (Advanced Cell Diagnostics, Biotechne). All H&E and ISH stained slides were digitally scanned using a Panoramic 3D-Histech scanner and viewed using CaseViewer v2.4 software. The presence of viral RNA by ISH was evaluated using the whole lung tissue section slides. Digital image analysis was performed in RNAscope labelled slides to ascertain the percentage of stained cells within the lesions, by using the Nikon-NIS-Ar software package.
Viral load quantification by RT-qPCR
RNA was isolated from nasal swabs and throat swabs. Samples were inactivated in AVL (Qiagen) and ethanol. Downstream extraction was then performed using the BioSprint™96 One-For-All vet kit (Indical) and Kingfisher Flex platform as per manufacturer’s instructions. Tissues were homogenized in Buffer RLT+ betamercaptoethanol (Qiagen). Tissue homogenate was then centrifuged through a QIAshredder homogenizer (Qiagen) and supplemented with ethanol as per manufacturer’s instructions. Downstream extraction from tissue samples was then performed using the BioSprint™96 One-For-All vet kit (Indical) and Kingfisher Flex platform as per manufacturer’s instructions.
Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) targeting a region of the SARS-CoV-2 nucleocapsid (N) gene was used to determine viral loads and was performed using TaqPath™ 1-Step RT-qPCR Master Mix, CG (Applied Biosystems™), 2019-nCoV CDC RUO Kit (Integrated DNA Technologies) and QuantStudio™ 7 Flex Real-Time PCR System. Sequences of the N1 primers and probe were: 2019-nCoV_N1-forward, 5’ GACCCCAAAATCAGCGAAAT 3’; 2019-nCoV_N1-reverse, 5’ TCTGGTTACTGCCAGTTGAATCTG 3’; 2019-nCoV_N1-probe, 5’ FAM-ACCCCGCATTACGTTTGGTGGACC-BHQ1 3’. The cycling conditions were: 25°C for 2 minutes, 50°C for 15 minutes, 95°C for 2 minutes, followed by 45 cycles of 95°C for 3 seconds, 55°C for 30 seconds. The quantification standard was in vitro transcribed RNA of the SARS-CoV-2 N ORF (accession number NC_045512.2) with quantification between 1 and 6 log copies µL-1. Positive swab and fluid samples detected below the limit of quantification (LoQ) of 4.11 log copies mL-1, were assigned the value of 5 copies µL-1, this equates to 3.81 log copies mL-1, whilst undetected samples were assigned the value of < 2.3 copies µL-1, equivalent to the assay’s lower limit of detection (LoD) which equates to 3.47 log copies mL-1. Positive tissue samples detected below the limit of quantification (LoQ) of 4.76 log copies mL-1 were assigned the value of 5 copies µL-1, this equates to 4.46 log copies g-1, whilst undetected samples were assigned the value of < 2.3 copies µL-1, equivalent to the assay’s lower limit of detection (LoD) which equates to 4.76 log copies g-1.
Subgenomic RT-qPCR was performed on the QuantStudio™ 7 Flex Real-Time PCR System using TaqMan™ Fast Virus 1-Step Master Mix (Thermo Fisher Scientific) and oligonucleotides as specified by Wolfel et al.[27], with forward primer, probe and reverse primer at a final concentration of 250 nM, 125 nM and 500 nM respectively. Sequences of the sgE primers and probe were: 2019-nCoV_sgE-forward, 5’ CGATCTCTTGTAGATCTGTTCTC 3’; 2019-nCoV_sgE-reverse, 5’ ATATTGCAGCAGTACGCACACA 3’; 2019-nCoV_sgE-probe, 5’ FAM- ACACTAGCCATCCTTACTGCGCTTCG-BHQ1 3’. Cycling conditions were 50°C for 10 minutes, 95°C for 2 minutes, followed by 45 cycles of 95°C for 10 seconds and 60°C for 30 seconds. RT-qPCR amplicons were quantified against an in vitro transcribed RNA standard of the full length SARS-CoV-2 E ORF (accession number NC_045512.2) preceded by the UTR leader sequence and putative E gene transcription regulatory sequence described by Wolfel et al [27]. Positive samples detected below the lower limit of quantification (LLOQ) were assigned the value of 5 copies µL-1, whilst undetected samples were assigned the value of ≤0.9 copies µL-1, equivalent to the assays lower limit of detection (LLOD). For nasal swab, throat swab and BAL samples extracted samples this equates to an LLOQ of 4.11 log copies mL-1 and LLOD of 3.06 log copies mL-1. For tissue samples this equates to an LLOQ of 4.76 log copies g-1 and LLOD of 3.71 log copies g-1.
Plaque reduction neutralisation test
Neutralising virus titres were measured in heat-inactivated (56°C for 30 minutes) serum samples. SARS-CoV-2 was diluted to a concentration of 1.4 x 103 pfu mL-1 (70 pfu in 50 µl) and mixed 50:50 in 1% FCS/MEM with doubling serum dilutions from 1:10 to 1:320 in a 96-well V-bottomed plate. The plate was incubated at 37°C in a humidified box for one hour to allow the antibody in the serum samples to neutralise the virus. The neutralised virus was transferred into the wells of a washed 96-well plate containing virus-susceptible VERO/E6 cells, allowed to adsorb at 37°C for a further hour, and overlaid with 1% CMC in 10% FCS/MEM. After a 24-hour incubation at 37°C in a humified box, the plates were formaldehyde-fixed before infected cell foci were detected by immunostaining with anti-SARS-CoV-2 RBD spike protein, a polyclonal rabbit horseradish peroxidase conjugate and TrueBlue Substrate. The percentage reduction of foci in serum compared to virus only control was calculated using SoftMax Pro v7.0.3GxP and neutralization titres (ND50) are reported as the serum dilution that neutralized 50% of the virus foci.
Antigen Binding ELISA
Recombinant SARS-CoV-2 Spike- and RBD-specific IgG responses were determined by ELISA. A full-length trimeric and stabilised version of the SARS-CoV-2 Spike protein was supplied by Lake Pharma (#46328). Recombinant SARS-CoV-2 Receptor-Binding-Domain (319-541) Myc-His was developed and kindly provided by MassBiologics.
High-binding 96-well plates (Nunc Maxisorp, 442404) were coated with 50 µl per well of 2 µg mL-1 Spike trimer (S1+S2) or RBD in 1x PBS (Gibco) and incubated overnight at 4°C. The ELISA plates were washed and blocked with 5% Foetal Bovine Serum (FBS, Sigma, F9665) in 1x PBS/0.1% Tween 20 for 1 hour at room temperature. Serum collected from animals after vaccination had a starting dilution of 1:50 followed by 8 two-fold serial dilutions. Post-challenge samples were inactivated in 0.5% triton and had a starting dilution of 1:100 followed by 8 three-fold serial dilutions. Serial dilutions were performed in 10% FBS in 1x PBS/0.1% Tween 20. After washing the plates, 50 µl per well of each serum dilution was added to the antigen-coated plate in duplicate and incubated for 2 hours at room temperature. Following washing, anti-monkey IgG conjugated to HRP (Invitrogen, PA1-84631) was diluted (1: 10,000) in 10% FBS in 1X PBS containing 0.1% Tween 20 and 100 µl per well was added to the plate. Plates were then incubated for 1 hour at room temperature. After washing, 1 mg mL-1 O-Phenylenediamine dihydrochloride solution (Sigma P9187) was prepared and 100 µl per well were added. The development was stopped with 50 µl per well 1M Hydrochloric acid (Fisher Chemical, J/4320/15) and the absorbance at 490 nm was read on a Molecular Devices versamax plate reader using Softmax (version 7.0). Titers were determined using the endpoint titer determination method. For each sample, an endpoint titer is defined as the reciprocal of the highest sample dilution that gives a reading (OD) above the cut-off. The cut-off was determined for each experimental group as the mean OD + 3SD of naïve samples.
Peripheral Blood Mononuclear Cell isolation and resuscitation
PBMCs were isolated from whole blood anticoagulated with heparin (132 Units per 8 720 ml blood) (BD Biosciences, Oxford, UK) using standard methods. PBMCs isolated from tissues were stored at −180 °C. For resuscitation PBMCs were thawed, washed in R10 medium (consisting of RPMI 1640 supplemented with 2 mM L-glutamine, 50 U/ml penicillin- 50 μg mL-1 streptomycin, and 10% heat-inactivated FBS) with 1 U mL-1 of DNase (Sigma), and resuspended in R10 medium and incubated at 37°C 5% CO2 overnight.
ELISpot
An IFNγ ELISpot assay was used to estimate the frequency and IFNγ production capacity of SARS-CoV-2-specific T cells in PBMCs using a human/simian IFNγ kit (MabTech, Nacka. Sweden), as described previously[28]. The cells were assayed at 2 x 105 cells per well. Cells were stimulated overnight with SARS-CoV-2 peptide pools spanning the ECD spike protein. Five peptide pools were 748 used, comprising of 15mer peptides, overlapping by 9 amino acids. Phorbol 12-myristate (Sigma) (100 ng mL-1) and ionomycin (CN Biosciences, 753 Nottingham, UK) (1 mg mL-1) were used as a positive control. Results were calculated and reported as spot forming units (SFU) per million cells. All SARS-CoV-2 peptides were assayed in duplicate and media only wells subtracted to give the antigen-specific SFU. ELISPOT plates were analysed using a CTL scanner and software (CTL, Germany) and further analysis carried out using GraphPad Prism (GraphPad Software, USA).
Statistics
All statistical analyses were performed using GraphPad Prism 7 or 8 software (La Jolla, CA). These data were considered significant if p <0.05. The type of statistical analysis performed is detailed in the figure legend. No samples or animals were excluded from the analysis.