Rift Valley Fever Reemergence after 7 Years of Quiescence, South Africa, May 2018.
Communicable disease surveillance and outbreak investigation activities of the National Institute for Communicable Diseases (Johannesburg, South Africa) are approved by the Human Research Ethics Committee of the University of the Witwatersrand, Johannesburg, South Africa (M160667). In mid-May 2018, an outbreak of RVF in sheep on a farm in Free State Province, South Africa, was reported, followed by 4 probable cases in humans detected by RVFV serology (6). The affected farm is located in Jacobsdal District, a farming community close to the border of Northern Cape Province, where sheep are the main livestock species. In addition to the 6 patients sampled on May 21, 2018, described previously (6), another 4 were sampled on June 4, 2018. These patients experienced headache, muscle pain, fever, body ache, rigors, and nausea, as reported previously (6). A recent history of influenza-like illness was reported for only 2 of these 4 patients. All 4 lived and worked on the farm and were involved in high-risk activities, such as slaughtering, autopsying, disposal and burial of carcasses, or handling of raw meat. We obtained follow-up samples from all 10 patients for paired serologic testing (Appendix Table 1, https://wwwnc.cdc.gov/EID/ article/25/2/18-1289-App1.pdf).
We performed the serologic assays, hemagglutination inhibition assay, virus neutralization test, and IgM ELISA with all serial serum samples collected from all 10 patients (7) and real-time reverse transcription PCR (RT-PCR) (8) on the serum fractions of clotted blood collected from the first 6 patients with suspected cases described previously (6). We extracted nucleic acid from EDTA whole blood samples collected from the initial 4 patients with probable cases using the MagMax Pathogen RNA/DNA Kit (Applied Biosystems, https://www.thermofisher.com) and then tested by RT-PCR. We determined the partial genome sequences of viruses from RT-PCR-positive whole blood samples using sequence-independent single-primer amplification combined with sequencing in triplicate on an Illumina MiSeq (https://www.illumina.com) and raw data processing, as described previously (9). After quality and host filtering and using a requirement of [greater than or equal to] 3x coverage per base, we mapped raw reads to reference sequences representing the RVFV large (L), medium (M), and small (S) segments (GenBank accession nos. KX611605-7). We concatenated sequence fragments of segments, prepared alignments in MEGA6 (https://www.megasoftware.net), performed phylogenetic analyses using RAxML version 8.2.10 (http:// evomics.org/learning/phylogenetics/raxml), and visualized trees with Figtree version 1.4.3 (http://tree.bio.ed.ac.uk/ software/figtree).
Of 10 patients sampled, 8 seroconverted after 2 or 3 serial bleeds, as evidenced by a 4-fold increase in the hemagglutination inhibition assay or virus neutralization test titers (Appendix Table 1), and had RVFV-specific IgM, confirming their recent RVFV infection status. We detected RVFV RNA in EDTA whole blood samples of 3 of 4 patients sampled 7 days after estimated symptom onset (Appendix Table 1). Sequence-independent single-primer amplification sequencing yielded sequence fragments of the M and L segments in 1 (SA344-18) of 3 samples (GenBank accession nos. MH753234-41). We concatenated these partial sequence fragments (Appendix Table 2) and attained 86% (3,341/3,885 nt) of the M segment and 77.7% (4,975/6,404 nt) of the L segment. Only a single fragment was obtained of the S segment (776 nt, 45.9%) spanning nucleotides 18-793 (GenBank accession no. MH753233). Partial sequence fragments of the L segment from another sample (SA343-18) were also obtained; this sequence had a 165-bp overlap with sample SA344-18 (at nucleotides 3,276-3,440) and a single-nucleotide mismatch (A3305T, 99.4% identical) but was not included in phylogenetic analyses because of its close identity to SA344-18 and small fragment size. We prepared alignments with similarly concatenated L and M or partial S sequences from GenBank (Appendix Table 3) and a separate alignment with a 490-nt portion of the M segment of SA344-18 and sequences available from GenBank from a previous study, including sequences from the 2010 RVF outbreak in South Africa (5).
In phylogenetic analyses of the L and M segments, isolate SA344-18 grouped with lineage E (Figure 1; Appendix Figure 1), according to the lineage nomenclature of Bird et al. (4). In analyses of all segments, this RVFV isolate had the closest relationship to Beijing-01, an isolate from a 2016 case exported to China from Angola (10), and to a 2009 sheep isolate from Kakamas, Northern Cape Province, South Africa (11) (Figures 1, 2; Appendix Figures 1, 2). The E lineage also included 25010-24, an isolate from a camel sampled during an outbreak in Mauritania in 2010 (Appendix Figure 1) (11). Nodes of the partial S segment tree were poorly supported, probably because the fragments were small (Appendix Figure 2). The overall pairwise nucleotide differences were 3.1% for L, 3.6% for M, and 3.1% for S (nonstructural protein gene), similar to values reported previously (4,5).
According to the partial 490-nt M segment tree and lineage nomenclature of Grobbelaar et al. (5), SA344-18 is lineage K (Figure 2). All 2010 RVFV sequences from humans in South Africa, except SA184-10, group in lineage H, distant from SA344-18. Isolate SA184-10 (lineage K) was obtained from a patient who had a needle-stick injury while vaccinating sheep with Smithburn vaccine (5). The 490-nt partial M sequence of SA184-10 is identical to that of the 2009 Kakamas sheep isolate (Figure 2). The only RVFV sequence from a human in Northern Cape Province (SA404-09, GenBank accession no. HM587096) grouped in lineage H (data not shown) (5).
The major RVF epidemic in South Africa during 20102011 included mostly lineage H RVFV isolates, and the smaller outbreaks during the 2 preceding years included lineage C (5). The virus detected in Free State Province in 2018 groups in the distant lineage K, despite extensive transmission of lineage H during 2010-2011 in the same province. The 2018 isolate, SA344-18, detected in the absence of a major epidemic, is closely related to viruses from isolated cases, such as Kakamas-2009 and Beijing-01. These isolates group in the lineage that includes Smithburn vaccine viruses and its parent strain, Entebbe, isolated in Uganda in 1944 (4). Camel isolate 25010-24 also falls in this lineage, despite most other RVFV isolates from Mauritania grouping with distant lineages (11).
Intensified surveillance during interepidemic periods would enhance our knowledge of the genetic evolution of RVFV. Isolated outbreaks, similar to the one reported here, probably have occurred more frequently than recorded. Establishment of sentinel herds of susceptible animals and increased vector surveillance could yield isolates more frequently and lead to a better characterization of the strains maintained at low levels between epidemics. Our results indicate an ongoing activity and evolution of RVFV during interepizootic periods and highlight the importance of a cryptic transmission cycle that enables establishment of endemicity, which at times inevitably emerges in the form of an explosive outbreak.
The work presented here was partly funded by the Poliomyelitis Research Foundation (grant no. PRF12/63) and partly through the routine arbovirus reference diagnostic work of the Centre for Emerging Zoonotic and Parasitic Diseases, National Institute for Communicable Diseases.
Dr. Jansen van Vuren is a medical scientist at the National Institute for Communicable Diseases in Johannesburg, South Africa. His research interests are diagnostics and viral hemorrhagic and arboviral diseases.
(1.) Nanyingi MO, Munyua P, Kiama SG, Muchemi GM, Thumbi SM, Bitek AO, et al. A systematic review of Rift Valley fever epidemiology 1931-2014. Infect Ecol Epidemiol. 2015;5:28024. http://dx.doi.org/10.3402/iee.v5.28024
(2.) Archer BN, Thomas J, Weyer J, Cengimbo A, Landoh DE, Jacobs C, et al. Epidemiologic investigations into outbreaks of Rift Valley fever in humans, South Africa, 2008-2011. Emerg Infect Dis. 2013;19:1918-25. http://dx.doi.org/10.3201/eid1912.121527
(3.) Pienaar NJ, Thompson PN. Temporal and spatial history of Rift Valley fever in South Africa: 1950 to 2011. Onderstepoort J Vet Res. 2013;80:384. http://dx.doi.org/10.4102/ojvr.v80i1.384
(4.) Bird BH, Khristova ML, Rollin PE, Ksiazek TG, Nichol ST. Complete genome analysis of 33 ecologically and biologically diverse Rift Valley fever virus strains reveals widespread virus movement and low genetic diversity due to recent common ancestry. J Virol. 2007;81:2805-16. http://dx.doi.org/10.1128/ JVI.02095-06
(5.) Grobbelaar AA, Weyer J, Leman PA, Kemp A, Paweska JT, Swanepoel R. Molecular epidemiology of Rift Valley fever virus. Emerg Infect Dis. 2011;17:2270-6. http://dx.doi.org/10.3201/ eid1712.111035
(6.) Jansen van Vuren P, Kgaladi J, Patharoo V, Ohaebosim P, Msimang V, Nyokong B, et al. Human cases of Rift Valley fever in South Africa, 2018. Vector Borne Zoonotic Dis. 2018;vbz.2018.2357. http://dx.doi.org/10.1089/vbz.2018.2357
(7.) Paweska JT, Burt FJ, Swanepoel R. Validation of IgG-sandwich and IgM-capture ELISA for the detection of antibody to Rift Valley fever virus in humans. J Virol Methods. 2005;124:173-81. http://dx.doi.org/10.1016/j.jviromet.2004.11.020
(8.) Drosten C, Gottig S, Schilling S, Asper M, Panning M, Schmitz H, et al. Rapid detection and quantification of RNA of Ebola and Marburg viruses, Lassa virus, Crimean-Congo hemorrhagic fever virus, Rift Valley fever virus, dengue virus, and yellow fever virus by real-time reverse transcription-PCR. J Clin Microbiol. 2002;40:2323-30. http://dx.doi.org/10.1128/ JCM.40.7.2323-2330.2002
(9.) Ladner JT, Wiley MR, Mate S, Dudas G, Prieto K, Lovett S, et al. Evolution and spread of Ebola virus in Liberia, 2014-2015. Cell Host Microbe. 2015;18:659-69. http://dx.doi.org/10.1016/ j.chom.2015.11.008
(10.) Liu J, Sun Y, Shi W, Tan S, Pan Y, Cui S, et al. The first imported case of Rift Valley fever in China reveals a genetic reassortment of different viral lineages. Emerg Microbes Infect. 2017;6:e4. http://dx.doi.org/10.1038/emi.2016.136
(11.) El Mamy AB, Lo MM, Thiongane Y, Diop M, Isselmou K, Doumbia B, et al. Comprehensive phylogenetic reconstructions of Rift Valley fever virus: the 2010 northern Mauritania outbreak in the Camelus dromedarius species. Vector Borne Zoonotic Dis. 2014;14:856-61. http://dx.doi.org/10.1089/ vbz.2014.1605
Address for correspondence: Janusz T. Paweska, Center for Emerging Zoonotic and Parasitic Diseases, National Institute for Communicable Diseases of the National Health Laboratory Service, 1 Modderfontein Rd, Sandringham 2131, Johannesburg, South Africa; email: firstname.lastname@example.org
Author affiliation: National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
Caption: Figure 1. Maximum-likelihood tree showing the phylogeny of Rift Valley fever virus isolate SA344-18, collected in South Africa in May 2018, on the basis of the concatenated large segment. Lineage names according to the nomenclature of Bird et al. (4) are indicated. Maximum-likelihood analysis was performed in RAxML version 8.2.10 (http://evomics.org/learning/phylogenetics/raxml); 10,000 bootstrap replicates were performed. Bootstrap values are shown at nodes. Scale bar indicates nucleotide changes per base pair. BFA, Burkina Faso; CAR, Central African Republic; EGY, Egypt; GIN, Guinea; KEN, Kenya; MAD, Madagascar; MAU, Mauritania; MAY, Mayotte; PRC, China; RSA, South Africa; SAU, Saudi Arabia; SUD, Sudan; TAN, Tanzania; UGA, Uganda; ZIM, Zimbabwe.
Caption: Figure 2. Maximum-likelihood tree showing the phylogeny of Rift Valley fever virus isolate SA344-18, collected in South Africa in May 2018, on the basis of a 490-nt fragment of the medium segment. Lineage names according to the nomenclature of Grobbelaar et al. (5) are indicated. Maximum-likelihood analysis was performed in RAxML version 8.2.10 (http://evomics. org/learning/phylogenetics/ raxml); 100,000 bootstrap replicates were performed. Bootstrap values are shown at the nodes. Scale bar indicates nucleotide changes per base pair. ANG, Angola; CAR, Central African Republic; EGY, Egypt; GIN, Guinea; KEN, Kenya; MAD, Madagascar; MAU, Mauritania; NAM, Namibia; PRC, China; RSA, South Africa; SAU, Saudi Arabia; SEN, Senegal; SOM, Somalia; UGA, Uganda; ZAM, Zambia; ZIM, Zimbabwe.
|Printer friendly Cite/link Email Feedback|
|Author:||van Vuren, Petrus Jansen; Kgaladi, Joe; Msimang, Veerle; Paweska, Janusz T.|
|Publication:||Emerging Infectious Diseases|
|Date:||Feb 1, 2019|
|Previous Article:||Bat Influenza A(HL18NL11) Virus in Fruit Bats, Brazil.|
|Next Article:||Tick-Borne Encephalitis Virus Antibodies in Roe Deer, the Netherlands.|