Journal of Virology & Antiviral ResearchISSN: 2324-8955

All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

Research Article, J Virol Antivir Res Vol: 7 Issue: 2

First Identification of Orungo Virus (Orbivirus) from a Febrile Patient in South-Eastern Senegal, 2015

Martin Faye1*, Arame Ba1, Yakhouba Kane1, Andre Pouedewou Bedekelabou1, Mamadou Diop2, Moussa Moise Diagne1, Gamou Fall1, Oumar Faye1, Amadou Alpha Sall1 and Ousmane Faye1

1Virology Pole, Institut Pasteur of Dakar, 220 Dakar, Senegal

2Bioinformatics and Statistics group G4, Institut Pasteur of Dakar, 220 Dakar, Senegal

*Corresponding Author : Dr Martin Faye, PhD
36 avenue Pasteur, Institut Pasteur of Dakar, Dakar, Senegal
Tel: +221.33 839 92 23
Fax:
+221.33 839 92 00
E-mail:
martin.faye@pasteur.sn

Received: October 01, 2018 Accepted: November 16, 2018 Published: November 23, 2018

Citation: Faye M, Ba A, Kane Y, Bedekelabou AP, Diop M, et al. (2018) First Identification of Orungo Virus (Orbivirus) from a Febrile Patient in South- Eastern Senegal, 2015. J Virol Antivir Res 7:2. doi: 10.4172/2324-8955.1000187

Abstract

Orbiviruses are of medical interest for humans and livestock worldwide. We isolated a novel Orungo virus strain from in a febrile patient living in Southeastern Senegal. The full-length genome was determined using MiSeq next-generation sequencing and the phylogenetic relationship to other orbivirus groups was analyzed. Inoculation experiments in mammalian Vero cells and mice were conducted and a specific real time reverse transcription-polymerase chain reaction assay was developed.

This is the first isolation of Orungo virus in a human in Senegal in addition to the prototype strain. This new characterized Orungo virus genome presents ten segments with the typical genomic organization of orbiviruses. However, important amino acid gaps were identified on it VP1 and VP3 segments. Phylogenetic analysis confirms that this Senegalese Orungo isolate is closely related to Culicoides-borne orbiviruses. In addition, this isolate is able to infect mammalian Vero cells and is lethal for mice when inoculated intra-cerebrally. The developed reverse transcription-polymerase chain reaction assay can be used for rapid and specific detection of Orungo virus.

This study would be helpful not only for future genomic studies to understand the biological and pathogenic properties of Orungo virus, but also for better understanding of its evolutionary history of orbiviruses and assessment of their potential public health impact.

Keywords: First identification; Orungo virus; orbivirus; Human; Senegal

Introduction

Orungo virus (ORUV) belongs to the genus orbivirus in the family Reoviridae. Orbiviruses present a worldwide distribution, a large genetic diversity and are responsible for diseases of veterinary and human importance transmitted by a wide range of vectors including Culicoides and mosquito species [1]. Orbiviruses have a nonenveloped, segmented, double-stranded RNA (dsRNA) genome coding 11 distinct proteins (VP1 to VP7 and NS1 to NS4) organized in 10 segments, the NS4 overlapping with the VP6 protein. Each segment is flanked by short 5’ and 3’ non-coding regions (NCRs) [1,2]. Although vaccines for Bluetongue virus (BTV), African horse sickness virus (AHSV) and Epizootic hemorrhagic disease virus (EHDV) exist in animals, no antiviral drugs are currently available for humans [3,4]. ORUV is a mosquito-transmitted orbivirus first isolated in 1959 from Anopheles funestus mosquitoes in Uganda and later from 8 humans and Aedes dentalus mosquitoes in Nigeria [5,6]. ORUV presents a large diversity with four distinct serotypes (ORUV-1 to ORUV-4) that are widely distributed in tropical Africa [7]. It is responsible of a febrile disease with symptoms including headache, myalgia, nausea and vomiting and induces antibodies production between day 7 and day 14 post infection [5,7,8]. In addition, ORUV kills less than 10-days old mice and hamsters following intra-cranial inoculation and cause cytopathic effect (CPE) in Vero and BHK-21. However, it didn’t infect adult mice, hamsters and chickens by subcutaneous inoculation [8].

In Africa, without BTV circulation notified in Algeria and Morocco in 2008 and 2009, respectively [5], no major orbivirus outbreak has been yet reported. Indeed, although high prevalence of ORUV neutralizing antibodies was reported in humans in CAR (88%) in 1979 and Southeastern Gabon (58%) in 1982, in Gambian population (up to 23%) in 1978-1979 and in food animals of Northern Nigeria (up to 30%) in 1984 [9-11], reported clinical human cases still rare. Herein, we report the first identification and isolation of ORUV from a patient in Southeastern Senegal in 2015 during arboviruses surveillance and the results of the comparative genome and phylogenetic analysis and mice experiments.

Materials and Methods

Sample collection and sequencing

In February 2015, during arboviruses surveillance in Kédougou region, Southeastern Senegal, a blood sample was obtained from a 31- years-old man with a 3-days history of fever (38 °C) and headache (SH264542). The patient lived in the department of Kédougou at the onset of disease and had not traveled abroad Figure 1. The sample was analyzed at WHO collaborating centre (WHO-CC) for arboviruses and hemorrhagic fevers at the Institut pasteur of Dakar (IPD) using molecular and serological diagnostic methods for arboviruses and a neutralizing antibodies test for yellow fever virus was positive. The sample was later processed for high-throughput RNA sequencing at IPD. RNA was extracted from the serum using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany). Briefly, viral RNA was fragmented and the Nextera XT library Prep kit (Illumina, San Diego, CA, USA) was used for strands synthesis, adapter ligation and libraries amplification according to the manufacturer’s recommandations (Illumina, San Diego, CA, USA). The amplified DNA library was analyzed using the 2100 Bioanalyzer instrument (Agilent Technologies, Santa Clara, CA, USA). The purified sequencing library was quantified using the KAPA library quantification kit (Kapa Biosystems, Wilmington, DE, USA) and subsequently sequenced on the Illumina MiSeq platform with 151 bp paired-end reads using the Miseq reagents kit v3 according to the manufacturer’s protocol (Illumina, San Diego, CA, USA). Because the protocol followed made use of randomly primed cDNA synthesis, it provided a comprehensive and quantitative view of all RNA present in the sample, thus enabling the characterization of a potential present virus. Illumina sequence reads were quality trimmed by Prinseq-lite, and de novo genome assembly was performed with IDBA (Iterative de Bruijn graph De Novo Assembler).

Figure 1: Sample collection site in Southeastern Senegal. The department of Kédougou is highlighted in brown and bordered by the departments of Salémata and Saraya colored in pink and purple, respectively.

Comparative genome analysis

Overlapping nucleotide sequences obtained from the 10 dsRNA genome segments of this human ORUV were assembled using the Unipro UGENE software [12]. Multiple alignments of nucleotide sequences were carried out by using Muscle algorithm [13] within Unipro UGENE software. Based on these alignments, we investigated the genetic properties of this senegalese ORUV isolate, such as segments length and amino acid distances between available ORUV coding regions from Uganda in 1959 (isolate UGMP 359, Genbank acc. no. JQ610675-84) and sequences of the new characterized ORUV isolate.We also assessed amino acid differences in coding regions.

Phylogenetic tree inference

A Bayesian phylogenetic analysis for estimation of data quality and selection of the best-fit nucleotide substitution model were performed using Model Finder test [14] on the IQ-TREE 1.6.3 web server [15]. The maximum likelihood (ML) tree was then constructed with sequences of VP1 (RNA-dependent RNA polymerase) and VP3 (T2, Major subcore Protein) proteins representing the two most conserved genes between orbiviruses. Tree was inferred using Fast Tree v2.1.7 [16] with the best fit nucleotide substitution model to our sequence data. Nodes were labeled with local support values, which were computed with the Shimodaira-Hasegawa test (SH-like) for 5000 replications. Topology was visualized by FigTree v.1.4.2.

Viral isolation in suckling mice and vertebrate cells

Twenty microliter of the filtrated serum was inoculated by intracerebral route into 10 Swiss suckling mice (1–2 days old) at animal laboratory of the IPD authorized by the national ethical committee in Senegal to use virus isolation in suckling mice for surveillance, routine diagnostic and research of viruses. Infected mice were followed for 21 days post infection (dpi). In addition, 150 μl of filtrated serum were inoculated onto confluent Vero cells (African Green Monkey Kidney) monolayer (25 cm) in Leibovitz 15 (L-15) growth medium supplemented with 5% fetal bovine serum (FBS) (both GibcoBRL, Grand Island, NY, USA), and 10 mM each of Penicillin and Streptomycin (Sigma, Gmbh, Germany). Inoculated cells were incubated at 37 °C for a maximum of 7 days and cells were collected until observation of a CPE. In order to assess ORUV infection, cells were tested by indirect immunofluorescence assay (IFA) using specific hyper-immune mouse ascitic fluid for different orbiviruses including Orungo, Palyam, Acado, Kindia, Gomoka, Andesibe, Pata, Ndelle, Ngoupe and Sanar, as previously [17].

RT-PCR for orbivirus detection

Viral RNA was extracted from initial serum, culture supernatant from infected Vero cells and suckling mice brain using the QIAamp viral RNA mini kit (Qiagen, Heiden, Germany) according to manufacturer’s instructions. Extracted RNA was tested by real-time reverse transcription-polymerase chain reaction (RT-PCR) as previously described [18] using Bluetongue virus and Corriparta virus as positive controls. Subsequently, 5 μL of each PCR products were analyzed by gel electrophoresis on 1 % agarose gels stained with ethidium bromide to check the size of amplified fragments by comparison to a DNA molecular weight marker (HyperLadder™ 1kb, Bioline, USA). The DNA bands from the PCR amplification were purified (QIAquick Gel Extraction Kit, Qiagen, Heiden, Germany) and sequenced (Beckmann Coulter, UK).

Development of a real time RT-PCR for Orungo virus detection

In order to design an assay capable of detecting ORUV, both primers and TaqMan probe were designed on VP1 segment of the ORUV sequence available from Uganda (GenBank Ac.No.JQ610675) using Primer3web® software (version 4.0.0, Whitehead Institute for Biomedical Research). To avoid non-specific cross-reactions with others orbiviruses, primers and probe were validated by BLAST analysis on NCBI. Oligonucleotides were synthetized by TIBMol-Biol (Berlin, Germany) (Table 1).

Name Type Length Sense Sequence 5'-3' Gene Position GC (%) Tm (°C)b Product size (nt)
ORUV_F3022 primer 21 sense CAAGTGAGAGCATTGCCACTG VP1pol 3022-3042 50.00 59.04 238
ORUV_P3052 probe 21 sense 6FAM-TCGCCTTTCAATCTCTCCGCA--BBQ VP1pol 3052-3072 52.38 67.48
ORUV_R3259 primer 21 Antisense CACTCGCGTAGAGAGATCAAC VP1pol 3259-3279 55.00 59.14

Table 1: Oligonucleotide sequences of primers and probe designed for Orungo virus detection.

Briefly, real-time RT-PCR was performed using the Quantitect Probe RT-PCR kit (Qiagen, Hilden, Germany) in a 20 μl reaction mixture containing 2x QuantiTect probe RT-PCR master mix, 0.25 μl of QuantiTect RT mix, 400 nM PCR primers, 200 nM TaqMan probe and 5 μl RNA in RNase-free water. Positive control containing master mix with ORUV RNA and negative control consisting of master mix with sterile RNase-free water were included in each run. The reaction was carried out on a 7500 Fast Real Time system cycler or on a 7500 Real Time system cycler (Applied Biosystems, Foster City, US) using the following temperature profile: 1 cycle of reverse transcription at 50 °C for 10 min, 1 cycle of denaturation at 95 °C for 15 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. In addition, a total of 11 other major representatives of orbiviruses including Bluetongue, Corriparta, Palyam, Acado, Kindia, Gomoka, Andesibe, Pata, Ndelle, Ngoupe and Sanar and obtained from the archive of the WHO-CC at IPD were also used for specificity testing. The initial serum, culture supernatant from infected Vero cells and suckling mice brain were tested in duplicates by real time RT-PCR.

Plaque assay and mice experiments

Ten-fold serial dilutions (10-1 to 10-11) of the serum sample were prepared in maintenance medium L-15 supplemented with 5% FBS and virus titer in initial serum was determined in plaque forming units by milliliter (PFU/mL) by inoculation onto Vero cells as previously described [19]. Furthermore, 20 μl Vero-grown ORUV containing 103 PFU/mL of the virus was inoculated into 4 Swiss adult mice (up to 4-weeks old) by group through intracerebral, intraperitoneal and subcutaneous route. Control groups were inoculated with 20 μl of 1x phosphate buffer saline (PBS). The animals were observed daily for illness and followed for 21 dpi and the experiments were undertaken two times. Finally, sera, brain, liver and spleen from animals euthanized at 21 dpi were tested in duplicates for ORUV infection by real time RT-PCR.

Results

Comparative genome analysis

Surprisingly, de novo assembly of the sample reads revealed an Orungo virus genome from this human serum (ORUV; SH264542) and assembled segments sequences were deposited in GenBank (Accession numbers: MH477874-83). The likelihood of sample contamination is small, since the sequencing run hadn’t included any material containing ORUV nor none of the additional RNA samples which were treated at the same time provided sequence reads that corresponded to ORUV. Analysis of complete coding sequences of the new characterized ORUV was performed at amino acid level including previously available sequences from Uganda in 1959 (Genbank accession number JQ610675-84) into multiple sequences alignment. The protein length was determined with respect to coding regions size and amino acid distances were ranged from 0 to 6%, with the highest distances observed at proteins VP1, VP3 and VP6 (Table 2).

Segment Comparative analysis with available ORUV sequences Genbank Accession number Nucleotide length Protein names: Amino acid length Amino acid distance Genbank Accession number
Segment 1 JQ610675 3756 VP1 : 1252/1302 6% MH477874
Segment 2 JQ610676 2880 VP2 : 960/960 2% MH477875
Segment 3 JQ610677 2649 VP3 :883/906 4% MH477876
Segment 4 JQ610678 1911 VP4 : 637/637 1% MH477877
Segment 5 JQ610679 1638 NS1 : 546/546 0% MH477878
Segment 6 JQ610680 1557 VP5 : 519/519 2% MH477879
Segment 7 JQ610681 1050 VP7 : 350/350 0% MH477880
Segment 8 JQ610682 1029 NS2 : 343/343 1% MH477881
Segment 9 JQ610683 858 VP6 : 286/286 6% MH477882
Segment 10 JQ610684 681 NS3 : 227/227 0% MH477883

Table 2: Comparative analysis of coding regions lengths of the new human ORUV isolate.

In addition, the highest numbers of non-conservative amino acid substitutions were located on proteins VP1, VP2, VP3 and VP6 with 18, 11, 9 and 9 mutations, respectively. Interestingly, sequences of the new characterized ORUV isolate showed gaps of 50 and 23 amino acids at VP1 and VP3, respectively when compared with sequences from Uganda (Table 3).

Segment Protein names Amino acid position Ugandan ORUV isolate Senegalese ORUV isolate Presence of gaps
Segment 1 VP1 12 R I gap of 50 aa between aa positions
263-312
184 S P
188 S P
197 F S
205 M I
315 F V
339 M L
557 S N
633 L F
648 N Y
650 K N
677 T A
732 N S
738 I F
822 I F
827 D Y
921 P L
1011 T A
Segment 2 VP2 183 S N  
211 L F
246 H D
289 Q R
478 T A
571 P S
633 N S
712 T A
788 N S
795 L S
898 H D
905 S P
Segment 3 VP3 8 G R gap of 23 aa between aa positions
180-202
31 G S
55 A P
61 L W
239 A T
366 A T
378 G C
386 T I
425 T I
Segment 4 VP4 476 H Y  
Segment 5 NS1 - - -  
Segment 6 VP5 360 Y N  
363 P T
490 F L
495 L M
508 P Q
514 L R
Segment 7 VP7 - - -  
Segment 8 NS2 164 D G
235 R C
Segment 9 VP6/NS4 15 H D  
38 P S
43 D G
78 G R
87 T I
130 L P
185 L F
192 P S
261 T I
Segment 10 NS3 - - -  

Table 3: Summary of non-conservative amino-acid changes in the various putative coding regions encoded by the Senegalese ORUV genome relative to those of the isolate from Uganda in 1959.

Phylogenetic analyses

The General Time-Reversible with a Gamma distribution of 4 categories rate (GTR+F+G4) was the best nucleotide substitution model for our sequences data presenting scores of 103,978.11 and 95,928.15, 103,590.62 and 95,566.50, 103,592.20 and 95,567.87 for Bayesian information criterion score (BIC), Akaike information criterion score (AIC) and corrected AIC score (AICc) criteria using our VP1 and VP3 sequence datasets, respectively. The maximum likelihood (ML) trees was inferred using FastTree v2.1.7 [16] and phylogenetic analysis revealed that this new characterized ORUV (ORUV_SN_SH264542) belongs to Culicoides-borne orbiviruses group, clustering closely with the ORUV isolate sampled in Uganda in 1959 (Figure 2).

Figure 2: Maximum Likelihood (ML) trees inferred on complete sequences from VP1 (A) and VP3 (B) proteins of orbivirus groups. Phylogenetic trees were constructed using FastTree and nodes were labeled with local support values computed using the Shimodaira- Hasegawa (SH) test for 5000 bootstrap replications. The new characterized Senegalese Orungo virus isolate belonging to the Culicoides-borne orbiviruses group was highlighted in red on these trees.

Viral isolation and identification

We observed mortality at 4 dpi in suckling mice inoculated with 20 μl of the initial serum. In addition, the inoculated Vero cells developed a CPE within 6 dpi, indicating virus growth. Immunofluorescence assay showed that this serum cross-reacted only with ORUV specific hyper-immune ascitic fluid and not with other tested orbiviruses important in African context (Figure 3).

Figure 3: Assessment of Orungo virus infection by Immunofluorescence assay (IFA). Specificity of the Orungo virus infection was assessed using specific immune ascetic fluids for ten different orbiviruses on cells inoculated with the serum and positive and negative controls (from 1 to 10).

RT-PCR amplification

The designed primer pair and probe exhibited no cross-reactivity with other species of orbivirus genus. The initial serum, culture supernatant from infected Vero cells and suckling mice brain were tested positive by conventional and real time RT-PCR and confirmed the ORUV positive serum (ORUV_SN_SH264542). However, the virus was not detected in any sera and tissues from adult mice euthanized at 21 dpi (Table 4).

Samples RT-PCR results (Ct values)
SH264542 initial serum 28.15 ± 0.20
SH264542-vero-cells-supernatant 22.62 ± 0.56
SH264542-Suckling-mice-brain tissues 25.59 ± 0.72
Orungo virus 23.50 ± 0.43
Bluetongue virus Negative
Corriparta virus Negative
Palyam virus Negative
Acado virus Negative
Kindia virus Negative
Gomoka virus Negative
Andesibe virus Negative
Pata virus Negative
Ndelle virus Negative
Ngoupe virus Negative
Sanar virus Negative
Control group Spleen tissues Negative
Control group Liver tissues Negative
Control group Brain tissues Negative
Control group Sera Negative
IC group Spleen tissues Negative
IC group Liver tissues Negative
IC group Brain tissues Negative
IC group Sera Negative
IP group Spleen tissues Negative
IP group Liver tissues Negative
IP group Brain tissues Negative
IP group Sera Negative
IM group Spleen tissues Negative
IM group Liver tissues Negative
IM group Brain tissues Negative
IM group Sera Negative

Table 4: Specificity assessment of Orungo virus assay and tropism evaluation Viral characterization.

The plaque assay determined a titer of 5.75 x 104 PFU/ml in the initial serum (ORUV_SN_SH264542) and no CPE was observed in the negative control culture (Figure 4). The ORUV inoculum was tested negative for mycoplasma sequences using PCR (data not shown). However, no illness was observed at 21 dpi after intracerebral, subcutaneous and intraperitoneal inoculation of adult mice.

Figure 4: Orungo virus plaque assay: A1: negative control (cell culture medium), A2: initial human serum. A3-D6: 10-fold serial dilutions in duplex wells down to 10-11 of ORUV human serum on Vero cell culture monolayers. The initial human serum as well as the Vero cells was previously confirmed by PCR to be negative for mycoplasma sequences (data not shown).

Discussion and Conclusion

Our study is the first in Senegal identifying human infection caused by Orungo virus, an orbivirus responsible for diseases of veterinary and human importance. We found evidence of locally acquired infection with ORUV in Senegal through antigenic typing and sequencing, which raises the issue of the risk of circulation of the ORUV and human infection in Southeastern Senegal. These data also highlight the potential usefulness of high-throughput sequencing, particularly untargeted approaches, for pathogen surveillance [20]. Our phylogenetic analysis indicated a strong VP1 and VP3 sequence homology between this Senegalese ORUV isolate and the isolate from Uganda, despite evidence of amino acid differences observed in these major proteins such as gaps. Cumulated with evident gene exchanges (reassortment) occurring in orbiviruses replication cycle [2,21-23], these sequences modifications may lead to differences in viral fitness and cause more or less severity in illness presentation. Thus, more studies need to be carried out on a large number of samples from a broader spectrum of diseases other than Yellow fever suspicion, with complete clinical records, to study the pathogenic determinants of this virus. In addition, other forces driving orbiviruses evolution such as intragenic recombination reported for Bluetongue virus [24], could be also assessed in future studies focusing on ORUV diversity.

Although the Senegalese ORUV isolate showed lethality to suckling mice when injected intra-cerebrally, results of pathogenicity experiments in adult mice indicated that ORUV illness can’t occur from the intra-cerebral, intra-peritoneal and intra-muscular routes of injection as in a previous study [5]. In addition, viral absence in tested tissues could be explained by ORUV tropism different to brain, liver and spleen or viral clearance by the host’s immune system. However, it could be interesting to perform further experiments with higher viral doses in less than 10-days old mice inoculated by subcutaneous route or in other animal models such as sheep.

Since both the vectors and suitable hosts (e.g., Culicoides) for orbiviruses are present in Senegal [25-27], more studies of this viral genus in this locality, including serosurveillance, are needed. Increasing levels of population movement inside Senegal and between Southeastern Senegal and other counties which features a rich diversity of arthropod species implicated in the maintenance and transmission of arboviruses of medical interest [28,29] may provide opportunities for pathogens to expand their geographical distribution. Cumulatively, this isolation from febrile human indicates that ORUV may be a health risk for humans and other mammals, and future strategies of vaccine or antiviral development against orbivirus genus in Africa should consider these possibilities.

Despite a previous suggestion that the dsRNA viruses have evolved by co-evolution with their respective hosts [30], recent studied provided evidence for co-evolution of the orbiviruses with their arthropod vectors rather than their vertebrate host [2,31]. Thus, it could be interesting to assess also vector competency of potential ORUV vector groups circulating in Southeastern Senegal [3,32].

Since ORUV could emerge or re-emerge in geographically distinct locations, outside its original geographical boundaries as other orbiviruses or with new clinical onsets, this new molecular diagnostic assay have been developed to assist in broad identification of ORUV and diagnosis in livestock and/or wildlife, and to enhance a rapid response to this emerging pathogen, if necessary. Our complete genome could provide a resource for further studies on biological and pathogenic properties and its possible association with human diseases, even during outbreaks or circulations of arbovirus species such as Yellow fever virus.

Acknowledgment

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors and was only supported by the Institut Pasteur of Dakar. We acknowledge colleagues of virology pole of Institute Pasteur de Dakar, Senegal for sharing serum and supportive information necessary for establishment and accomplishment of this study. “The authors declare no conflict of interest."

Author Contributions

Martin Faye, Amadou Alpha Sall and Ousmane Faye conceived and designed the experiments. Martin Faye, Arame Ba, Andre Pouwedeou Bedekelabou, and Yakhouba Kane performed the experiments. Martin Faye, Mamadou Diop and Moussa Moise Diagne analyzed the data. Gamou Fall, Oumar Faye, Amadou Alpha Sall and Ousmane Faye contributed reagents/materials. All authors wrote the paper. All authors read and approved the final manuscript.

References

international publisher, scitechnol, subscription journals, subscription, international, publisher, science

Track Your Manuscript

Awards Nomination

Media Partners