Our findings provide a rationale for the continued testing of investigational therapies, confirm the effectiveness of the diagnostic assays used in the region, and establish a paradigm for the use of reverse genetics to inform response activities in an outbreak.

The phylogenetic analysis of whole-genome sequences from each Ebola virus outbreak suggests there are at least two Ebola virus strains in DR Congo, which have independently crossed into the human population. The Ituri Ebola strain initially grew slower than the Makona strain, yet reached similar mean yields of 3 × 10 7 50% tissue culture infectious dose by 72 h infection in Huh-7 cells. Ituri Ebola virus was similar to Makona in its susceptibility to inhibition by remdesivir and to neutralisation by monoclonal antibodies from ZMapp and other monoclonal antibodies. Remdesivir inhibited Ituri Ebola virus at a 50% effective concentration (EC 50 ) of 12nM (with a selectivity index of 303) and Makona Ebola virus at 13nM (with a selectivity index of 279). The Zmapp monoclonal antibodies 2G4 and 4G7 neutralised Ituri Ebola virus with a mean EC 50 of 0·24 μg/mL and 0·48 μg/mL, and Makona Ebola virus with a mean EC 50 of 0·45 μg/mL and 0·2 μg/mL. The Xpert Ebola and US Centers for Disease Control and Prevention real-time RT-qPCR diagnostic assays detected Ituri and Makona Ebola virus sequences with similar sensitivities and efficiencies, despite primer site binding mismatches in the Ituri Ebola virus.

Using Ebola virus sequences provided by organisations in DR Congo and a reverse genetics system, we generated an authentic Ebola virus from the ongoing outbreak in Ituri and North Kivu provinces. To relate this virus to other Ebola viruses in DR Congo, we did a phylogenetic analysis of representative complete Ebola virus genome sequences from previous outbreaks. We evaluated experimental therapies being tested in clinical trials in DR Congo, including remdesivir and ZMapp monoclonal antibodies, for their ability to inhibit the growth of infectious Ituri Ebola virus in cell culture. We also tested diagnostic assays for detection of the Ituri Ebola virus sequence.

The ongoing Ebola virus outbreak in the Ituri and North Kivu Provinces of the Democratic Republic of the Congo, which began in July, 2018, is the second largest ever recorded. Despite civil unrest, outbreak control measures and the administration of experimental therapies and a vaccine have been initiated. The aim of this study was to test the efficacy of candidate therapies and diagnostic tests with the outbreak strain Ituri Ebola virus. Lacking a virus isolate from this outbreak, a recombinant Ituri Ebola virus was compared with a similarly engineered Makona virus from the 2013–16 outbreak.

Because of the containment requirements, surrogate assays and in silico analyses are often used to predict biological activities for EBOV. Although these can be informative, we sought to provide a more definitive evaluation of candidate therapies and diagnostic assays by testing them using infectious EBOV-Ituri. Lacking an isolate from this outbreak, the aim of this study was to rescue recombinant, infectious EBOV-Ituri from plasmid DNA, and to use this virus for our evaluations. The International Committee on Taxonomy of Viruses has established guidelinesfor Ebola nomenclature and declared EBOV to indicate the Ebola Zaire species only. Ituri is a variant of the Ebola Zaire species; we have referred to other Ebola virus species by name (eg, Bundibugyo).

Despite security concerns in the area, scientists with the Institut National de Recherche Biomédicale (INRB; Kinshasa, Democratic Republic of the Congo) and the United States Army Medical Research Institute of Infectious Diseases (USAMRIID; Fort Detrick, MD, USA) established a sequencing facility in DR Congo in May, 2018. They succeeded in generating high-quality, accurate sequences from the Ituri-North Kivu outbreak, which they shared with the scientific community. Their work allowed the early identification of Zaire EBOV as the causative agent of the outbreak. In partnership with local and international collaborators, analyses of these outbreak sequences elucidated early transmission chains and identified the index case to be from the Mangina health district.

Our data support the continued testing of remdesivir and ZMapp in DR Congo. Our work illustrates the utility of reverse-genetic systems to generate outbreak strains, allowing the testing of investigational therapies and diagnostic assays against emerging virus sequences to provide a rational basis for response activities that may have otherwise been absent.

Ebola virus vaccines, therapeutics, and diagnostic assays have been developed on the basis of strains from past outbreaks. Using our reverse genetics system and Ebola virus sequences provided by organisations in DR Congo, we generated an authentic Ebola virus from the ongoing outbreak in Ituri and North Kivu provinces. We offer evidence that the Ebola-Ituri virus is phylogenetically distinct from past outbreak viruses, is inhibited by experimental antiviral compounds and monoclonal antibodies being tested in DR Congo, and can be detected by diagnostic assays used in the field.

We searched PubMed using the terms “ebola” or “filovirus” combined with “DRC” or “outbreak” for articles published up to April 14, 2019, with no restrictions on language or publication date. Reports were mostly commentaries and focused on the socioeconomic aspects of the outbreak in the Democratic Republic of the Congo, particularly the difficulties with community engagement, public trust for participating in contact tracing, vaccination, and reporting to Ebola treatment centres. We found no research articles reporting the sensitivity of the Ituri-North Kivu Ebola virus to experimental therapies and its detection by diagnostic assays used in the current outbreak. After we had submitted our manuscript, two articles were released in The Lancet Infectious Diseases describing the Ebola virus sequences obtained from both the Équateur Province and Ituri-North Kivu 2018 outbreaks in DR Congo. Mbala-Kingebeni and colleagues sequenced Ebola virus locally and used the results to track transmission and test deployed diagnostics and medical countermeasures in silico.

In DR Congo, an extensive immunisation campaign with the experimental vaccine rVSV-ZEBOV is underway, but effective treatments for Ebola virus disease are urgently needed. In North Kivu, a randomised clinical trial (registered with ClinicalTrials.gov , number NCT03719586 ) has been initiated to compare mortality among patients receiving four different investigational EBOV therapies: monoclonal antibody cocktails ZMapp (developed by Mapp Biopharmaceuticals, San Diego, CA, USA) and REGN-EB3 (developed by Regeneron Pharmaceuticals, Eastview, NY, USA), single monoclonal antibody mAb114 (developed by the National Institutes of Health, Bethesda, MD, USA), and the small molecule remdesivir (also known as GS-5734, developed by Gilead Sciences, Foster City, CA, USA). The monoclonal antibodies bind to the EBOV glycoprotein (GP) to neutralise the virus, while remdesivir is a nucleotide analogue prodrug that inhibits the viral polymerase.

Viruses of the genus Ebolavirus cause sporadic outbreaks of severe haemorrhagic fever, with case fatality of up to 90%.On July 30, 2018, the Ministry of Health of the Democratic Republic of the Congo declared their tenth outbreak of Ebola virus disease, in the Mabalako health district in North Kivu. This was just weeks after declaring the end of the ninth outbreak in the western town of Bikoro, Équateur Province, 780 miles away. Responding to any filovirus outbreak is challenging and requires a combination of clinical, laboratory, and epidemiological approaches to halt virus transmission. Community engagement and education are key to enabling contact tracing and safe and dignified burials. The mineral-rich area of eastern DR Congo is experiencing a humanitarian crisis that has lasted more than a decade, hampering response activities. Despite best efforts of the Ministry of Health, Médecins Sans Frontières, WHO, and many other organisations, the ongoing outbreak in the Ituri and North Kivu provinces is now the second largest and longest Ebola virus (EBOV) outbreak ever recorded, with 2181 cases and 1459 deaths (as of June 17, 2019).

The sponsor of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Statistical significance was determined using the Holm-Sidak method with α=0·05 for the means of each virus titre each timepoint after infection. To calculate the IC 50 , we did a non-linear regression analysis with a variable slope and four parameters. The means of the IC 50 s were calculated for at least two independent experiments and the standard deviations reported.

The highest posterior density interval with 95% of the probability of the posterior distribution is reported for the Bayesian analysis of genetic sequence data by the BEAST program. The number of independent samples from the posterior distribution was estimated in the effective sample size in Tracer. The maximum clade credibility tree was calculated from the posterior sample with the maximum product of posterior probabilities of each clade in TreeAnnotator.

To test GeneXpert, virus RNA was extracted from TriPure isolation reagent (Sigma-Aldrich, St Louis, MO, USA) with Direct-zol RNA purification kits (Zymo Research, Irvine, CA, USA) and MagMax Viral RNA isolation kits (ThermoFisher). The extracted RNA was spiked into matrix and the Xpert Ebola assay was done according to manufacturer's instructions on the GeneXpert XVI system (Cepheid, Sunnyvale, CA, USA). The CDC EBOV NP and VP40 real-time RT-qPCR assays were done with SuperScript III reverse transcriptase and HiFi Platinum Taq one-step kit (ThermoFisher) in a Bio-Rad Touch CFX 96 thermal cycler (Bio-Rad, Hercules, CA, USA) for the following cycling conditions: reverse transcription at 50°C for 15 min, inactivation of reverse transcription at 95°C for 2 min, 40 cycles of 95°C for 30 s to denature, and 55°C for 60 s for annealing and probe detection.

For compound testing, fluorescence readings from quadruplicate wells at each compound concentration were taken. Analyses of signals from antibody-inhibition experiments were done similarly, but with triplicate wells for each antibody concentration. Background signals (no virus) were deducted from fluorescence readings, and data normalised to vehicle-only control. GraphPad Prism (version 8.1.2; GraphPad Software, La Jolla, CA, USA) was used to generate concentration–response plots. A four-parameter equation was used to fit semi-log plots of the data and derive the 50% effective concentration (EC 50 ) or 50% cytotoxic concentration (CC 50 ). The selectivity index was defined as the CC 50 divided by the EC 50 . The data shown are from at least two independent experiments.

Antibody neutralisation activity was measured in Huh-7 cells. Antibodies were diluted in DMEM with an initial dilution of 1:100 and a serial four-fold dilution for eight points. A constant amount of virus, calculated for a multiplicity of infection of 0·4, was added to all antibody dilutions and allowed to incubate for 2 h at 37°C, then added to Huh-7 cells seeded the day before at 10 000 cells per well in 96-well black-walled, CellCarrier Ultra plates. After 24 h, the medium was removed, the cells were washed with phosphate-buffered saline, and fixed in formalin. NucBlue was added to stain for nuclei and ZsG fluorescence measured per cell using an Operetta high-content imaging system with a 10× objective. Cell segmentation and ZsG intensity were quantified using Harmony (v4.1).

For fluorescence microscopy, Huh-7 cells were seeded, compounds were diluted and added, and cells were infected. 3 days after infection, cells were fixed with 10% formalin and nuclei were stained with NucBlue (ThermoFisher, Waltham, MA, USA). Cells were visualised using an Operetta high-content imaging system with a 10× objective and were quantified using Harmony (v4.1; PerkinElmer, Waltham, MA, USA).

Assays to measure the anti-EBOV activities of compounds were done as described previously.Briefly, Huh-7 cells were seeded in 96-well black-walled, CellCarrier Ultra plates (PerkinElmer) at 10 000 cells per well in FluoroBrite DMEM supplemented with 1× nonessential amino acids, 1× glutamax, and 5% FCS. The next day, unless noted otherwise, two-fold serial dilutions of compounds were generated in dimethylsulfoxide (DMSO), and then diluted in FluoroBrite media with 10% FCS before they were added to cells, to yield a final DMSO concentration of 0·5%. 2 h after the addition of compounds, the cells were infected with EBOV reporter viruses at a multiplicity of infection of 0·2. 3 days after infection, ZsG fluorescence was measured with a BioTek Synergy H1 plate reader (BioTek Instruments, Winooski, VT, USA) with excitation at 493 nm and emission at 505 nm, 6 mm read-height, and gain of 80.

Huh-7 cells were seeded at 5 × 10 5 cells per well in a six-well plate. The following day, cells were infected at a multiplicity of infection of 0·2. Samples were taken at various time-points (0, 24, 48, 72, 96, and 144 h post-infection), and the EBOV titre in each sample was assessed by a 50% tissue culture infectious dose assay in Huh-7 cells.

The Makona ZsG-VP40 plasmid based on EBOV strain H.sapiens-wt/LIB/2014/Makona-201403007 (GenBank accession KP178538) was described previously.The ZsGreen open reading frame and porcine teschovirus 2A self-cleaving peptide were cloned at the amino-terminus of VP40. Human codon-optimised genes for nucleo-protein (NP), VP35, VP30, and L (EBOV replication proteins) were cloned into pcDNA3.1. The full-length Ituri virus was based on sequence 18FHV089 (GenBank accession MK007329). This sequence in GenBank was incomplete, missing nine nucleotides at the 3ʹ end and seven at the 5ʹ end. To replace these, we used the sequence from EBOV strain H.sapiens-wt/LIB/2014/Makona-201403007, which we had previously confirmed using 3ʹ and 5ʹ rapid amplification of cDNA ends, and found the genomic ends to be conserved in Zaire EBOV. The complete Ituri virus genome was cloned behind a T7 promoter followed by a hepatitis δ virus ribozyme and T7 terminator. Briefly, infectious virus was rescued following transfection of Huh-7 cells seeded in 12-well plates with 1 μg pEBOV, 0·5 μg pC-L (or pC-L inactive), 0·5 μg pC-NP, 0·05 μg pC-VP35, 0·05 μg pC-VP30, and 1 μg of codon-optimised pC-T7. Supernatants were harvested 4 days after transfection, clarified by low speed centrifugation, and passaged twice on Huh-7 cells in 12-well plates. Recombinant viruses were sequenced by next-generation sequencing on the MiniSeq platform (Illumina, San Diego, CA, USA) and found to contain no changes after passage.

Complete genome sequences of Zaire EBOV were obtained from GenBank ( appendix p 2 ). The EBOV genomic sequences from the Ituri-North Kivu outbreak were generously made available before publication through GenBank, by INRB-USAMRIID scientists. Full-length EBOV genomic sequences were aligned using MUSCLE.Maximum likelihood was used to infer phylogeny among the two species of Ebola virus found in DR Congo, Zaire and Bundibugyo, using GARLI (v2.01) with the general time-reversible nucleotide substitution model with discrete γ distribution justified by jMODELTEST with Akaike information criterion.For analysis of the Zaire EBOV, a Bayesian coalescent phylogenetic analysis was used to infer a molecular clock to estimate the dates of origin. Evolutionary rates were estimated using Bayesian Markov chain Monte Carlo inference in BEAST (v.1.10.4).A Hasegawa-Kishino-Yano substitution model with γ distribution was chosen along with an uncorrelated relaxed clock with log-normal distribution and a constant population size model. A final analysis was run for 500 million chain lengths and analysed by Tracer (v.1.7.1) to achieve an effective sample size of over 200. TreeAnnotator (v.1.10.4) estimated the maximum clade credibility tree with a 10% burn-in, and FigTree (v1.4.4) was used to edit the final phylogenetic tree.

Human hepatocellular carcinoma Huh-7 cells (Apath, New York, NY, USA) were cultured in Dulbecco's modified Eagle medium (DMEM) supplemented with 1× nonessential amino acids and 10% fetal calf serum (FCS) at 37°C in 5% CO. The EBOV-Makona-ZsGreen1 reporter virus expressing ZsG was generated as described previously.Monoclonal antibodies 13F6 (0201–022), c13C6 (0201–023), c6D8 (0201–021), 4F3 (0201–020), and KZ52 (0260–001) were from IBT Bioservices (Rockville, MD, USA). The ZMapp prototype monoclonal antibodies 4G7 and 2G4 were donated by Gary Kobinger (Infectious Disease Research Center, Laval University, Quebec City, QC, Canada). Unlabelled IgG controls were obtained from Southern Biotech (Birmingham, AL, USA). Antibodies derived from convalescent patients were donated by Carl Davis and Rafi Ahmed (Emory University, Atlanta, GA, USA).

This study was designed to generate infectious EBOV-Ituri with which to evaluate candidate antiviral therapies and existing diagnostic tests. We generated a recombinant virus expressing a reporter protein, Zoanthus green fluorescent protein (ZsG), such that virus growth was quantifiable by measuring fluorescence in infected cells. Growth and inhibition of EBOV-Ituri was compared with that of similarly engineered EBOV-Makona. In compliance with US Centers for Disease Control and Prevention (CDC) policy, this study was reviewed and approved by the CDC Institutional Biosafety Committee before its initiation, and it is reported in accordance with the US Government Policy for Oversight of Life Sciences Dual Use Research of Concern.

The CDC Ebola assay approved by the US Food and Drugs Administration (FDA) for emergency use is a real-time RT-qPCR assay that detects two targets, the NP and VP40 EBOV genes. Aligning the primer and probe sequences against the EBOV-Ituri sequence, we found a mismatch in the middle of the NP2 reverse primer and in the VP40 reverse primer close to the 3ʹ end, at position 654 in the VP40 open reading frame ( appendix p 7 ). We evaluated the CDC Ebola assays alongside a VP40 assay with a redesigned primer (654G) completely complementary to the EBOV-Ituri sequence; both assays recognised EBOV-Ituri and EBOV-Makona sequences with similar sensitivities and efficiencies over the seven-log dilution series tested ( figure 3 appendix p 7 ).

(A) Results from the CDC NP2 real-time RT-qPCR assay, using dilution series of Ituri or Makona Ebola virus RNA. (B) Results from the CDC VP40 real-time RT-qPCR assay and the modified VP40 assay that accounts for the mismatch in the NP2 reverse primer (654G), using dilution series of Ituri or Makona Ebola virus RNA. For each assay, a representative of at least three independent experiments is shown. CDC=US Centers for Disease Control and Prevention.

Aligning the sequences available for the Ituri-North Kivu outbreak, we found a mismatch in one of the Xpert Ebola primer-binding sites in the NP gene (a T in the forward primer that binds 880-902 of the NP gene for EBO-Makona is substituted with C in the EBOV-Ituri sequence; appendix p 7 ). We isolated virus RNA from the EBOV-Ituri and EBOV-Makona viruses. When viral RNA was amplified using the Xpert Ebola assay, we found that both viruses had similar cycle threshold values ( table 3 ), indicating that this mismatch did not significantly affect detection of the EBOV-Ituri sequence.

EBOV-Ituri and EBOV-Makona were also similarly susceptible to inhibition by EBOV patient monoclonal antibodies 5.6.1A2, 9.6.3D6, and 2.1.1D5, whereas 2.1.1D7 did not neutralise either virus ( table 2 appendix pp 5–6 ).Other therapeutic antibodies being tested in the 2018 Democratic Republic of the Congo outbreak are monoclonal antibody mAb114 and antibody cocktail REGN-EB3 (comprised of R3470, R3471, and R3479 monoclonal antibodies).The cocrystal structure of EBOV-Mayinga GP and mAb114 implicate contact residues in GP important for mAb114 binding; these contact residues are conserved in the EBOV-Ituri sequence ( appendix p 4 ).Convalescent sera collected at Yambuku in 1976, as well as rabbit and non-human primate control sera, similarly inhibited the EBOV-Ituri and EBOV-Makona viruses ( table 2 appendix pp 5–6 ).

Within Zaire EBOV species, there is conservation of the GP amino acid sequence with only 31 differences between EBOV Ituri, Bikoro, Makona, and Mayinga, despite their isolations spanning 42 years ( appendix p 4 ). However, there is greater diversity between EBOV species, which might limit the activity of EBOV monoclonal antibodies. We found that antibodies KZ52, 4G7, and 2G4 could inhibit both EBOV-Ituri and EBOV-Makona in a concentration-dependent manner. The ZMapp monoclonal antibody cocktail is composed of 2G4, 4G7, and 13C6; earlier versions contained mAbs 6D8 and 13F6.KZ52 is considered the prototypical monoclonal antibody isolated from a patient in the 1995 Kikwit outbreak. The monoclonal antibodies 13C6, 6D8, and 13F6 do not inhibit either EBOV-Ituri or EBOV-Makona; these glycan cap and mucin binding antibodies have been shown to protect non-human primates from Ebola virus disease, but do not neutralise EBOV in cell culture ( table 2 appendix pp 4–6 ).

EBOV-Ituri and EBOV-Makona were equally sensitive to inhibition by remdesivir, with mean ECvalues of 13 nM (SD 2) for EBOV-Ituri and 12 nM (2) for EBOV-Makona ( table 1 appendix p 3 ). Mean CCvalues were 3·6 μM (SD 2·9) with a selectivity index of around 300 for both viruses, consistent with previous findings ( table 1 ).We also tested a panel of compounds previously reported to inhibit EBOV for their ability to inhibit EBOV-Ituri. Inhibition of infection by most compounds, including ribavirin, the selective oestrogen receptor modulator inhibitors toremifene and tamoxifen, and the calcium channel blocker verapamil was similar for EBOV-Ituri and EBOV-Makona ( table 1 ). The anti-arrhythmic drug amiodarone hydrochloride, was previously shown to block EBOV entry in vitro and was given to patients in west Africa,but its clinical benefits remain unclear. We found that inhibition of infection by amiodarone hydrochloride was similar for EBOV-Ituri and EBOV-Makona ( table 1 ).

Data are n or mean (SD), where n is the selectivity index. The mean was calculated from at least two independent determinations. The selectivity index was defined as the CC 50 divided by the EC 50 . EBOV=Ebola virus. CC 50 =50% cytotoxic concentration. EC 50 =50% effective concentration.

We generated plasmid DNA corresponding to an EBOV-Ituri sequence and included the ZsGreen reporter gene ( figure 2A appendix p 3 ). Under biosafety level 4 containment, we rescued a recombinant, infectious EBOV-Ituri reporter virus. We compared the replication of EBOV-Ituri and EBOV-Makona (similarly engineered to express ZsG) in Huh-7 cells ( figure 2 ). The yields of EBOV-Ituri were marginally but significantly lower than those of EBOV-Makona at early time-points, and increased seven-fold between 24 h and 48 h. At 24 h post-infection the mean yield was 4·64 × 1050% tissue culture infectious dose (TCID; SD 0) for EBOV-Makona and 3·21 × 10TCID(7·425 × 10) for EBOV-Ituri (p=0·96), and at 48 h post-infection the mean yield was 3·75 × 10TCID(0) for EBOV-Makona and 2·56 × 10TCID(8·56 × 10) for EBOV-Ituri (p=0·00068). At 72 h and later, virus yields reached similar mean yields of 3 × 10TCID. In addition, ZsG was readily detected when expressed by EBOV-Ituri and could be used as a reporter for screening of antiviral inhibitors ( appendix p 3 ).

(A) Schematic of the recombinant EBOV-Ituri genome in viral complementary sense with reporter ZsG expressed with the gene for virus membrane protein VP40 and separated by the gene for the self-cleaving 2A peptide. (B) Time course of viral growth. Huh-7 cells were infected at a multiplicity of infection of 0·2. Samples of media were harvested at varying time points and virus titres were measured. Data represent the mean of three biological replicates and the error bars represent the SD. EBOV=Ebola virus. TCID 50 =50% tissue culture infectious dose.

We did a phylogenetic analysis of whole-genome sequences from each EBOV outbreak ( figure 1 ); EBOV established four clades. Clade 1 includes viruses from the first identified EBOV outbreak in 1976–77 in the remote villages of Yambuku and Bonduni, north of the Congo River ( figure 1 ).Clade 4 includes the Makona variant viruses from the 2013–16 west African epidemic in Guinea, Liberia, and Sierra Leone. A simultaneous EBOV outbreak also occurred in 2014 in Boende (Tshuapa Province, Democratic Republic of the Congo), but these viruses are phylogenetically distinct, forming part of clade 3. EBOV from the 2018 outbreak in Bikoro (Democratic Republic of the Congo) was also in clade 3. The sequences from the current Ituri-North Kivu outbreak establish a distinct clade 2, along with viruses from a small outbreak in the northern health zone of Likati in DR Congo from 2017.

Bayesian coalescent phylogenetic analysis with Markov chain Monte Carlo sampling was used to infer a molecular clock to estimate the dates of origin. (A) Maximum-likelihood analysis of full-length genomes showing Ebola virus sequences cluster separately from Bundibugyo virus. Scale bar is substitutions per site. (B) Maximum clade credibility tree of the Bayesian analysis of the Ebola-Ituri sequence with those from other Ebola virus outbreaks. The most recent common ancestor, with years since 2018, and posterior probability values are shown at the nodes, and the 95% highest posterior density interval is shown in parentheses. The evolutionary rate was estimated to be 4·6 × 10 −4 nucleotide substitutions per site per year (2·8 × 10 −4 – 6·3 × 10 −4 highest posterior density), typical for an RNA virus with a 19 kb genome. The Ebola-Ituri genome used in this study is marked with a star. (C) Locations of Ebola virus outbreaks within the Democratic Republic of the Congo.

Discussion

4 Mbala-Kingebeni P

Aziza A

Di Paola N

et al. Medical countermeasures during the 2018 Ebola virus disease outbreak in the North Kivu and Ituri Provinces of the Democratic Republic of the Congo: a rapid genomic assessment. 25 Carroll SA

Towner JS

Sealy TK

et al. Molecular evolution of viruses of the family filoviridae based on 97 whole-genome sequences. The current EBOV outbreak in Ituri and North Kivu Provinces in DR Congo that began in July, 2018, is now the second largest ever recorded. Response to this outbreak continues to be extremely challenging, with a precarious security situation hampering case finding, contact tracing, decontamination, and other standard response activities. INRB-USAMRIID scientists have succeeded, however, in generating EBOV sequences from samples collected from patients in the region. We did a phylogenetic analysis with the Ituri-North Kivu and other EBOV sequences from past outbreaks and, consistent with the recent report by Mbala-Kingebeni and colleagues,found the Ituri viruses to be an independent introduction into the human population, unrelated to the previous Bikoro outbreak which ended just weeks before this tenth outbreak began. DR Congo is endemic for at least two species of Ebola virus: Zaire and Bundibugyo. Our analysis of the 2018 outbreaks suggests that there are at least two EBOV populations circulating separately in the virus' natural reservoir. The viruses from Bikoro in the west and Ituri in the east belong to different clades (Bikoro EBOV belongs to clade 3 and Ituri EBOV to clade 2) and they only share a distant common ancestor 39 years ago. Interestingly, for all clades, the Bayesian analysis suggests the variants of EBOV to be circulating 1–2 years before each outbreak. This is inferred from the most recent common ancestor at each node of an outbreak ( figure 1 ). Each outbreak is thought to be a separate introduction into the human population from an unknown natural reservoir. We do not know if the host is dispersed and virus transmitted among a widespread population or if virus persists in groups of mobile hosts, which must contend with geographical features, such as rivers or mountains. There might have been multiple transmission events from the host to other intermediate species before Ebola came into contact with humans. Within the EBOV species, there is 96% nucleotide identity among the viruses sequenced from 1976 to 2018. This is remarkable conservation for a negative-sense RNA virus, since RNA viruses typically have a rapid evolutionary rate. The exceptionally low diversity of all EBOV indicates genetic bottlenecks must be continually restraining the virus genome sequence.

26 Diehl WE

Lin AE

Grubaugh ND

et al. Ebola virus glycoprotein with increased infectivity dominated the 2013–2016 epidemic. , 27 Urbanowicz RA

McClure CP

Sakuntabhai A

et al. Human adaptation of Ebola virus during the west African outbreak. , 28 Dietzel E

Schudt G

Krähling V

Matrosovich M

Becker S Functional characterization of adaptive mutations during the west African Ebola virus outbreak. 29 Marzi A

Chadinah S

Haddock E

et al. Recently identified mutations in the Ebola virus-Makona genome do not alter pathogenicity in animal models. Although sequence analysis can guide the epidemiology of an EBOV outbreak, it can be difficult to predict functional aspects of a virus variant, such as sensitivity to antivirals and antibody treatments and correlates of protection in vaccine studies. There have been no EBOV isolates available to the scientific community from the past four outbreaks in DR Congo. We rescued a recombinant, infectious EBOV based on a recent circulating virus sequence in the Ituri-North Kivu region; this is only the third variant of EBOV recovered from a molecular clone. The EBOV-Ituri virus replicated slightly more slowly than EBOV-Makona, the virus responsible for the 2013–16 west Africa outbreak. Although the difference in growth was statistically significant, this marginal difference might not be biologically relevant and we have not confirmed if this difference is seen in other cell lines. Additionally, we confirmed the sequence of virus stocks used in this study but cannot exclude the presence of minor variants below the level of detection available with the Illumina method that might affect the growth of the viral population. The GP mutation Ala82Val proposed to be associated with increased infectivity and mortality in the west Africa epidemicis not present in the Ituri sequence. Although the Ala82Val change did not alter pathogenicity in animal models,its absence might explain the somewhat reduced replication of EBOV-Ituri relative to EBOV-Makona in cell culture.

Because of the stringent laboratory containment requirements for working with infectious EBOV, surrogate assays are often used to study small molecules and antibodies to infer anti-EBOV activities. Such assays include using minigenomes for studying viral replication and pseudotype viruses for studying cell entry, and often recapitulate what is observed with the authentic virus. However, findings from experiments using surrogate assays should be confirmed with authentic virus whenever possible. For example, the predominant GP form made by authentic virus is the secreted GP; the full-length form of GP is synthesised at lower frequency and is only expressed when the viral polymerase performs the non-templated addition of an extra adenosine. The secreted GP may bind antibodies, but is not present in pseudotype virus assays. Other assays make use of GP without a mucin-like domain or with a cleaved glycan cap, but might not accurately represent findings with an intact GP. Several investigational therapies are being tested in the current Ituri-North Kivu outbreak despite a lack of evidence for effectiveness of these therapies against the circulating virus. Here, we show that the investigational drug remdesivir and monoclonal antibodies present in the antibody cocktail ZMApp, as well as other small molecules and antibodies with therapeutic potential (humanised monoclonal antibodies, patient monoclonal antibodies, and convalescent patient sera), inhibited EBOV-Ituri similarly to the more commonly used EBOV-Makona isolate. Importantly, the humanised monoclonal antibodies present in the ZMapp cocktail blocked EBOV-Ituri and EBOV-Makona infection with a similar potency.

The Xpert Ebola assay, run on GeneXpert instruments from Cepheid , has received approval for emergency use both from WHO and FDA, and is being used extensively for diagnosis in the current outbreak. The assay uses self-contained cartridges for RNA extraction and real-time RT-qPCR on two different EBOV target genes, NP and GP. As with all assays that detect nucleic acid, sequence changes might impair assay performance, resulting in reduced sensitivity or even false-negative test results. It is therefore crucial to confirm assay performance on newly emerging outbreak sequences.

4 Mbala-Kingebeni P

Aziza A

Di Paola N

et al. Medical countermeasures during the 2018 Ebola virus disease outbreak in the North Kivu and Ituri Provinces of the Democratic Republic of the Congo: a rapid genomic assessment. The Xpert Ebola assay contains a mismatch in one of its NP primer-binding sites. Despite this, Mbala-Kingebeni and colleagueswere able to use a panel of EBOV-Ituri samples to show that the mean difference in cycle threshold values between the NP and GP assays was similar to the mean difference previously published for EBOV-Makona samples from the 2013–16 outbreak in Sierra Leone. Consistent with these observations, we showed that the Xpert Ebola assay detects the EBOV-Ituri sequence with similar sensitivity to EBOV-Makona, over a range of dilutions. Additionally, we showed that the CDC Ebola real-time RT-qPCR assay, which has a different mismatch, also detects both viruses with similar sensitivity. The satisfactory performance of these assays, despite the changes in their respective target sequences, is encouraging for their use in detecting future virus strains, but assay performance in any future outbreak should be confirmed using sequences specific to that outbreak.

One limitation of our study is that to facilitate the quantification of infection, we compared EBOV-Ituri and EBOV-Makona viruses encoding the fluorescent reporter protein ZsG. The inclusion of the ZsGreen gene means that these viruses are not ideal for testing in animal models to compare pathogenicity, for example, where it would be preferable to use viruses without an exogenous gene. Work to generate such a full-length EBOV-Ituri virus is underway. Additionally, although it seems unlikely, we cannot exclude the possibility that the inclusion of the ZsGreen gene might affect the growth of EBOV-Makona or EBOV-Ituri differently, which could affect the comparison of their growth kinetics.

Our results could help to inform outbreak response activities. Our data support the continued testing of remdesivir and ZMapp in DR Congo. We suggest that other investigational therapies such as monoclonal antibodies mAb114 and REGN-EB3 also be tested against the EBOV-Ituri strain as soon as possible, to provide a rational basis for their continued administration. We also show the validity of using authentic virus to test the performance of current diagnostic assays, and empirically prove the effect a primer mismatch might have on the efficiency of the assay.

Our work further shows the utility of viral genome sequencing during an outbreak. Such sequences could be used to supplement traditional epidemiological measures, identify chains of transmission, and help contact tracing. They could also be used to detect any known viral variants resistant to therapy and potential signatures of host adaptation. As technology progresses and rapid sequencing of virus strains in low-resource settings in the field becomes more and more feasible, if isolates are unavailable, we recommend a policy whereby reverse genetics is used to generate outbreak strains as a standard practice. As we have shown in this study, such recombinant viruses can be used to confirm the efficacy of investigational therapies and the sensitivity of diagnostic assays for use in the field. This testing might provide a rationale for the prioritisation of candidate therapeutics and could identify issues with diagnostic assays that are rectifiable by altering primer sequences.