Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Simultaneous circulation of two West Nile virus lineage 2 clades and Bagaza virus in the Zambezi region, Namibia

  • Heiko D. Guggemos,

    Roles Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Institute of Virology, Charité - Universitätsmedizin Berlin, corporate member of Free University Berlin, Humboldt-University Berlin, and Berlin Institute of Health, Germany, German Center for Infection Research (DZIF), Berlin, Germany

    ⨯
  • Matthias Fendt,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliations Institute of Virology, Charité - Universitätsmedizin Berlin, corporate member of Free University Berlin, Humboldt-University Berlin, and Berlin Institute of Health, Germany, German Center for Infection Research (DZIF), Berlin, Germany

    ⨯
  • Christian Hieke,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliation Institute of Virology, Charité - Universitätsmedizin Berlin, corporate member of Free University Berlin, Humboldt-University Berlin, and Berlin Institute of Health, Germany

    ⨯
  • Verena Heyde,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliation Institute of Virology, Charité - Universitätsmedizin Berlin, corporate member of Free University Berlin, Humboldt-University Berlin, and Berlin Institute of Health, Germany

    ⨯
  • John K. E. Mfune,

    Roles Resources, Writing – review & editing

    Affiliation Department of Biological Sciences, University of Namibia, Windhoek, Namibia

    ⨯
  • Christian Borgemeister,

    Roles Conceptualization, Funding acquisition, Project administration, Writing – review & editing

    Affiliation Center for Development Research (ZEF), University of Bonn, Bonn, Germany

    ⨯
  • Sandra Junglen

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing

    Sandra.junglen@charite.de

    Affiliations Institute of Virology, Charité - Universitätsmedizin Berlin, corporate member of Free University Berlin, Humboldt-University Berlin, and Berlin Institute of Health, Germany, German Center for Infection Research (DZIF), Berlin, Germany

    ⨯

Abstract

Flaviviruses include a great diversity of mosquito-borne arboviruses with epidemic potential and high global disease burden. Several flaviviruses are circulating in southern Africa affecting humans and livestock, among them West Nile virus (WNV) and Wesselsbron virus. Despite their high relevance, no arbovirus surveillance study has been conducted for more than 35 years in Namibia. In this study we assessed the diversity of flaviviruses circulating in mosquitoes in the densely populated, semi-tropical Zambezi region of north-eastern Namibia. In total, 10,206 mosquitoes were sampled in Bwabwata and Mudumu national parks and Mashi and Wuparo conservancies and screened for flavivirus infections. A high infection rate with insect-specific flaviviruses was found with 241 strains of two previously known and seven putative novel insect-specific flaviviruses. In addition, we identified ten strains of WNV in the main vector Cx. univittatus sampled in the Mashi conservancy. Surprisingly, the strains fell into two different clades of lineage 2, 2b and 2d. Further, three strains of Bagaza Virus (BAGV) were found in Cx. univittatus mosquitoes originating from Mudumu national park. Assessment of BAGV growth in different cell lines showed high replication rates in mosquito and duck cells and about 100,000fold lower replication in human, primate and rodent cells. We demonstrate a wide genetic diversity of flaviviruses is circulating in mosquitoes in the Zambezi region. Importantly, WNV and BAGV can cause outbreaks including severe disease and mortality in humans and birds, respectively. Future studies should focus on WNV and BAGV geographic distribution, as well as on their potential health impacts in and the associated social and economic implications for southern Africa.

Author summary

Mosquitoes serve as vectors for the transmission of infectious diseases. Some of the most important mosquito-borne arboviruses belong to the genus Flavivirus, which can induce severe disease in humans and livestock. Surveillance of vector populations provide information on circulating arboviruses and may help to identify local outbreaks. Here we sampled mosquitoes over three wet seasons in the densely populated, semi-tropical Zambezi region of north-eastern Namibia and tested them for infections with flaviviruses. We observed simultaneous circulation of two different West Nile virus clades in the main vector species Cx. univittatus. Humans infected with West Nile virus can develop flu-like symptoms or in rare cases meningoencephalitis. Further, we detected Bagaza virus in Cx. univittatus from another locality and season. Bagaza virus infects birds leading to high mortality rates and may also infect humans. Our data suggest that both viruses are endemic in the Zambezi region and may affect human health and well-being in Namibia.

Citation: Guggemos HD, Fendt M, Hieke C, Heyde V, Mfune JKE, Borgemeister C, et al. (2021) Simultaneous circulation of two West Nile virus lineage 2 clades and Bagaza virus in the Zambezi region, Namibia. PLoS Negl Trop Dis 15(4): e0009311. https://doi.org/10.1371/journal.pntd.0009311

Editor: Doug E. Brackney, Connecticut Agricultural Experiment Station, UNITED STATES

Received: December 21, 2020; Accepted: March 16, 2021; Published: April 2, 2021

Copyright: © 2021 Guggemos et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All sequence files are available from the NCBI GenBank database (accession number(s) MW383507, MW383508, MW383509, MW436414 - MW436420, MW672101).

Funding: This work was supported by the German Research Foundation (DFG, Collaborative Research Center 228: Future Rural Africa [TRR 228/1], grant to SJ and CB). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The genus Flavivirus within the family Flaviviridae contains some of the most dangerous and life-threatening mosquito-borne viruses like Dengue virus, Yellow fever virus, Japanese encephalitis virus, Zika virus and West Nile virus (WNV) [1]. Land use and socio-economic changes, as well as global trade and travel fueled the spread of these and other flaviviruses to new geographic regions where they caused epidemics and often became endemic [2]. For instance, Dengue virus is now endemic in >100 countries and affects about 400 million people each year [3].

Several flaviviruses are endemic in southern Africa, with WNV and Wesselsbron virus (WSLV) among the most prevalent ones [4]. WNV is a mosquito-borne flavivirus that was first isolated in 1937 from a patient with febrile symptoms in the West Nile district of Uganda [5]. The virus is maintained in nature in a mosquito-bird transmission cycle where it can be transmitted to humans and horses by a bite of an infected mosquito. WNV is geographically the most widely distributed encephalitic flavivirus and has caused outbreaks in the Americas, Africa, Europe and Asia [6–10]. Approximately 80% of human infections with WNV are asymptomatic. Symptomatic patients usually show flu-like symptoms, such as fever, fatigue, headache, nausea, muscle pain and rash [11,12]. Less than 1% of WNV patients develop severe symptoms, such as meningoencephalitis that can lead to death [13]. WNV can cause sporadic outbreaks, usually after heavy rainfall, resulting in a sharp vector population build-up [14]. In 1974, the Karoo and Northern Cape provinces in South Africa experienced a large WNV outbreak that affected ten thousand of humans though no cases of encephalitis and no deaths were recorded [15–17]. In South Africa WNV infections have been observed in a range of wild and domesticated animals and seroprevalence rates of 24%, 18% and 30% against WNV were detected in humans, cattle, and sheep, respectively using an in-house ELISA indicating that flaviviruses are highly prevalent across the country [18,19].

For WSLV, first isolated in the South African town of Wesselsbron in 1955, similar seroprevalence rates in humans and cattle were observed in South Africa [18,20]. WSLV causes arthralgia, myalgia and fever in humans [21,22]. Symptoms in cattle and sheep are more severe and similar to Rift Valley fever with abortion rates of up to 20% [23–25].

Bagaza virus (BAGV) is a mosquito-borne flavivirus pathogenic to birds and has only recently been detected in dead monal pheasants near Pretoria in South Africa [26]. Bagaza virus shares high sequence similarity with Israel turkey meningoencephalitis virus (ITV), which was first isolated in 1958 from domesticated turkeys in Israel [27] and was detected in South African turkeys in 1978 [28]. Based on the high sequence identity, both viruses have been suggested to belong to the same species [29]. Typical symptoms in susceptible bird species include weight loss, weakness, and apathy with mortality rates as high as 30% in partridges and pheasants [30,31]. BAGV was first isolated in 1966 from Culex mosquitoes collected in Bagaza, Central African Republic [32]. Subsequently, it has been detected in mosquitoes collected in Senegal and Mauritania [33–35] and India [36]. BAGV neutralizing antibodies were later detected in encephalitic patients from India suggesting that the virus can also infect humans [36]. BAGV was first isolated from vertebrates in 2010 from dead red-legged partridges and pheasants in southern Spain [31].

Usutu virus (USUV) and Banzi virus (BANV) are also endemic in the region and were both first isolated close to the Usutu river in South Africa [37–39]. Finally, there are historic reports of DENV epidemics in southern Africa in the late 19th and early 20th century, with most recent reports from northern Mozambique dating back to 1984 [40]. Since then, only few reports on DENV in the region have been published, primarily originating from imported cases, and hence DENV is most likely not circulating in this part of the continent [41].

Little is known on the current incidence of flaviviruses in Namibia. WNV, WSLV, BANV, Spondwenni and Yellow fever have been detected in the local population in serosurveys across the north-eastern Zambezi region in the 1950s and 1960s with antibody prevalence rates of 12%, 23%, 9%, 4 and 2% in neutralization assays [42]. In the 1980s, rates of 29% WNV and 5% WSLV seroprevalence were recorded by hemagglutination inhibition assays in the same region [43], and WNV was detected in a febrile patient in Namibia’s Ovamboland in 1989 [44]. Since independence in 1990, only one mosquito-borne arbovirus surveillance study was conducted in Namibia, focusing on the hot and semi-arid central region of the country, with a seroprevalence rate not exceeding 8% against flaviviruses reported in human sera by ELISA testing [45]. Yet since the 1980s, no comparable data has been gathered for the rather densely populated, semi-tropical regions in the north of the country [45].

In this study we sought to assess the diversity of flaviviruses circulating in mosquitoes in the Zambezi region in north-eastern Namibia. To obtain an overview on the genetic diversity of endemic viruses, we sought to sample a broad diversity of mosquito species from various land-use types.

Material and methods

Ethics statement

A research permit to conduct this research was received from the National Commission on Research, Science and Technology (permit number RPIV00442018).

Mosquito sampling and identification

Mosquitoes were sampled in north-eastern Namibia in the Sachinga Livestock Development Centre in May and June 2018, in Mudumu National Park and Wuparo Conservancy in November and December 2018, in Bwabwata National Park in February and March 2019, and in Mashi Conservancy in March and April 2019 and in February and March 2020 (Fig 1). Four sites each in the National Parks were selected to represent sites with medium and high pressure of elephants. In the conservancies, we selected four sites each in rangeland, agricultural field and in the conservancy core area to cover all existing land-use types of the conservancies. In total, two BG Sentinel traps (Biogents AG, Regensburg, Germany), two CDC Light traps and one CDC Gravid trap (both John W. Hock Company, Gainesville, United States of America) were operated from dusk till dawn for five consecutive days at each site (total trapping effort 1,480 trap nights). BG Sentinel traps were baited with BG-Lure (Biogents AG, Regensburg, Germany), 1-Octen-3-ol and worn socks on different days. The combination of trap types and attractants was used in order to sample a high diversity of mosquito species. Mosquitoes were stored in liquid nitrogen in the field. Mosquito species were identified morphologically in the laboratory using standard keys [46–48].

thumbnail
Download:
Fig 1. Schematic map of mosquito-sampling sites in Namibia.

Samples were retrieved from Bwabwata and Mudumu National Parks, Wuparo and Mashi Conservancies and Sachinga Livestock Development Centre (Sachinga LDC). The map base layer is based on https://www.naturalearthdata.com, and the map’s details are based on [94].

https://doi.org/10.1371/journal.pntd.0009311.g001

Viral RNA extraction and PCR-screening

Individual mosquito specimens were homogenized in 500 μl phosphate-buffered saline (PBS) using ceramic beads and a Tissue Lyser (QIAGEN, Hilden, Germany) and subsequently arranged into pools of 10 according to species and sampling location using 100 μl of the individual mosquito homogenate. The rest of the mosquito homogenate was retained for further analysis based on individual specimens. RNA was extracted from 200 μl of pooled mosquito homogenates using the MagNA Pure 96 DNA and Viral NA Small Volume Kit (Roche Diagnostics, Mannheim, Germany). The RNA was eluted in 100 μl of provided elution buffer. Subsequently cDNA synthesis was performed using SuperScript IV reverse transcriptase (Thermo Fisher Scientific GmbH, Dreieich, Germany) and random hexamer primers (Integrated DNA Technologies Germany GmbH, Munich, Germany). Mosquito pools were screened for flaviviruses using generic PCR assays targeting the RNA-dependent RNA polymerase gene according to published protocols [49,50]. Nucleotide sequences were obtained by Sanger sequencing (Microsynth Seqlab GmbH, Göttingen, Germany), analysed in Geneious 9.1.8 (Biomatters ltd., Auckland, New Zealand) and compared to GenBank using the Basic Local Alignment Search Tool (BLAST) (https://blast.ncbi.nlm.nih.gov/Blast.cgi). In addition, all mosquito pools were tested for WNV [51], Dengue virus [52], Zika virus [53], Yellow fever virus [54] and BAGV (forward 5’-TCCAGGGAAGACCAGAGAGG-3’, reverse 5’-AGGCTTCAGCAATCCTTCCC-3’, probe 5’-FAM-TGCGTTGAACACCTTTACCA-ZEN-3’) using virus specific real-time-PCRs. RNA was extracted from single mosquito specimen of virus-positive pools and tested individually for infection using virus-specific RT-qPCRs as described above.

Genome sequencing and analysis

A fragment of the viral E protein (1,336 nucleotides) of all detected WNV strains was amplified by nested PCR using the primers WNV-E F1 5’-GAGGGAGTGTCTGGAGCTAC-3’, WNV-E R1 5’-CTGTCACGGGCATTGATTCC-3’; WNV-E F2 5’-TGGGTTGATCTGGTACTGGA-3’, WNV-E R2 5’-GAAGTCCCTGTGTGATCCAG-3’. Full WNV genomes of two selected strains (M6848NA-2020 and M6646-NA-2020) were sequenced by amplifying overlapping PCR amplicons of 800–1,000 nucleotides in length using generic nested primers based on an alignment of WNV lineage 2 sequences (primer sequences will be published elsewhere). For full genome sequencing of the BAGV strain MP314-NA-2018 nested primer pairs were designed based on an alignment of all currently available BAGV genomes aiming to amplify overlapping fragments of approximately 1,500 nucleotides (S1 Table). Amplicons were sequenced by Sanger sequencing and genomes were assembled using closely related strains as reference genomes in Geneious 9.1.8 (Biomatters ltd., Auckland, New Zealand).

Phylogenetic analyses

Nucleotide and amino acid sequences of RdRp and E genes, as well as entire open reading frames (ORFs) were aligned with related viral sequences in Geneious using the CLUSTAL W [55] and MAFFT algorithms [56]. Phylogenetic trees were inferred by the maximum-likelihood (ML) method with the best suitable substitution matrix as identified by MEGA-X [57] using PhyML [58]. The phylogeny of all detected flaviviruses was inferred using the GTR model with a fixed Gamma shape parameter of 0.36 and an estimated fraction of invariable sites. The phylogenetic relationship based on WNV E genes was inferred with the GTR model using a fraction of 0.46 invariable sites and a fixed Gamma shape parameter of 1.2. The phylogenetic relationship based on WNV full genomes was inferred with the JC69 model using an estimated fraction of invariable sites and an estimated Gamma shape parameter. The BAGV phylogeny based on full BAGV genomes was inferred using the TN93 model, a fixed Gamma shape parameter of 0.27, and an estimated fraction of invariable sites and the BAGV tree based on RdRP sequence fragments was inferred with the TN93 model, a fixed Gamma shape parameter of 0.34, and an estimated fraction of invariable sites. Confidence testing was performed based on 1,000 bootstrap iterations.

Genetic barcoding for mosquito species identification and blood-meal analysis

From selected mosquito specimens a fragment of the cytochrome c oxidase I (COI) gene was amplified based on generic primers for invertebrates [59]. To identify the vertebrate feeding source of selected mosquitoes, a fragment of the COI gene was amplified using generic vertebrate specific primers [60]. Obtained PCR-products were sequenced by Sanger sequencing (Microsynth Seqlab GmbH, Göttingen, Germany) and compared to GenBank using BLAST (https://blast.ncbi.nlm.nih.gov/).

Virus isolation in cell culture

For virus isolation attempts, mosquito (C6/36, Aedes albopictus; CXT, Culex tarsalis) and mammalian (VeroE6, African green monkey) cell lines were either inoculated with an aliquot of the virus-positive mosquito-pool or an aliquot of the individual virus-positive mosquito homogenate as described previously [61]. Cultures were checked daily for occurrence of a cytopathic effect (CPE). An aliquot of the supernatant (100 μl) was used for infection of fresh cells eight days post infection (passage 1). This procedure was repeated four times. The cell culture supernatants of each passage were in addition checked for viral replication by real-time PCR as described below.

Virus stocks were prepared from infectious cell culture supernatant of MP312-NA-2018 passage 4, MP314-NA-2018 passage 4, and MP370-NA-2018 passage 5. The number of infectious particles was determined by Tissue culture Infectious Dose 50 (TCID50) end-point dilution assay [62,63].

Virus growth kinetics

Different vertebrate (VeroE6, African green monkey; BHK-21, hamster; KN-R, cattle; DF-1, chicken [64]; LMH, chicken [65]; AGE1.CR, duck [66]; HEK293, human) and insect (C6/36, Aedes albopictus; CXT, Culex tarsalis) cell lines were infected in duplicates at a multiplicity of infection (MOI) of 0.1 and 0.01, respectively. Aliquots of cell culture supernatants were sampled every 24 hours for five days and virus copy numbers were determined by virus specific RT-qPCR using a plasmid-based standard dilution series. Briefly, the 114 nucleotide BAGV qPCR amplicon was amplified by conventional PCR using infectious supernatant of MP314-NA-2018 and subsequently cloned into the TOPO TA vector using the TOPO TA cloning kit according to the manufacturer’s protocol (Thermo Fisher Scientific GmbH, Dreieich, Germany). An overnight culture of a positive clone was prepared in LB medium supplemented with ampicillin (100 μg/mL) and plasmid DNA was extracted with the QIAprep Spin Miniprep kit (QIAGEN, Hilden, Germany). Plasmid concentration was measured using a NanoDrop spectrophotometer (Thermo Fisher Scientific GmbH, Dreieich, Germany) and a tenfold serial dilution series from 1x101 – 1x106 genome copies/ml was used for quantification.

Results

Mosquito collection and flavivirus screening

A total of 10,206 mosquitoes were collected in national parks and conservancies in the Zambezi region, Namibia, from May 2018 to March 2020 (Fig 1). In total, 42 different species were collected with Culex univittatus (3032, 30%), Coquillettidia metallica (1851, 18%) and Mansonia uniformis (1652, 16%) as the most abundant species. About three quarters of the sample set represented female (n = 7840, 76.8%) and about one quarter male mosquitoes (n = 2181, 21.4%) and for 185 (1.8%) mosquitoes the sex could not be identified due to poor sample condition. The average catch size was 4.3 mosquitoes per catch bag. Most of the mosquitoes were collected in the Wuparo conservancy and low catch rates were associated with low precipitation rates (Fig 2). After individual homogenization, all mosquitoes were sorted into 1,021 pools, comprising 10 individuals each.

thumbnail
Download:
Fig 2. Overview of collected mosquitoes in the Zambezi region, Namibia.

Number of sampled mosquitoes per locality and collection time is given by bar charts according to trap type and amount of rainfall at the time of collection [69]. Number of collected Cx. univittatus mosquitoes per locality is shown by a dashed line.

https://doi.org/10.1371/journal.pntd.0009311.g002

Screening of all pools with a generic flavivirus RT-PCR yielded 246 viral RdRP sequences falling in eleven clusters. Two sequences from two pools of female Cx. univittatus mosquitoes showed 96.9% and 97.9% identity to the South African WNV strains HM147822 and EF429197, respectively, and three sequences from Culex sp. mosquitoes showed 99.3–100% identity to BAGV indicating the detection of WNV and BAGV in mosquitoes from Namibia. The term strain is in the following referred to as a variant of a given virus species. All other 241 sequences showed identities to insect-specific flaviviruses (ISF). We detected 10 strains of Culex flavivirus in Culex spp. mosquitoes and 60 strains of Mosquito flavivirus AMH010516 in Mansonia spp. and Cx. univittatus mosquitoes (Table 1). The other 171 sequences assembled in seven groups which showed pairwise identities of 70–77% to known ISF suggesting that they most likely belong to seven distinct species (Table 1). In addition, all mosquito pools were tested negative for DENV, ZIKV and YFV in specific qPCR assays.

thumbnail
Download:
Table 1. Detected insect-specific flaviviruses in the Zambezi region, Namibia.

https://doi.org/10.1371/journal.pntd.0009311.t001

Phylogenetic analyses based on one representative strain of every identified virus species and selected flaviviruses confirmed the identification of WNV, BAGV and nine ISFs (Fig 3). WNV M6646-NA-2020 and BAGV MP314-NA-2018 grouped together with Dengue virus 2 (NC001474) in a monophyletic clade comprising mosquito-borne flaviviruses in sister relationship to the clade of ISFs. The Namibian ISFs did not group according to their geographic origin, but rather formed host species associated clades.

thumbnail
Download:
Fig 3. Phylogenetic relationship of flaviviruses detected in the Zambezi region.

The phylogenetic tree is based on RdRP gene sequences of 153 nucleotides. Branch support is shown at each branch. Virus sequences detected in this work are shown in red. Gen Bank accession numbers are marked in grey.

https://doi.org/10.1371/journal.pntd.0009311.g003

Characterization of WNV

Following the identification of WNV in two mosquito pools, all pools were screened for WNV infections by specific RT-qPCR and individual mosquitoes of the eight WNV-positive pools were tested by specific RT-qPCR. WNV was identified in ten individual female Cx. univittatus mosquitoes with viral genome copy numbers of 1.3 x 105–6.8 x 107 per mosquito (Table 2). All attempts to isolate the WNV strains from the pooled or individual mosquito homogenates in cell culture failed despite the detection of high viral loads in the mosquitoes. Mosquito species were confirmed by genetic barcoding. We further examined the WNV-positive mosquito extracts (n = 10) for remnants of vertebrate blood meal sources. We identified human mitochondrial DNA in the mosquito that was infected with WNV strain M6848-NA-2020. No vertebrate DNA was found in the other mosquitoes.

thumbnail
Download:
Table 2. WNV strains identified in the Mashi conservancy in this study.

https://doi.org/10.1371/journal.pntd.0009311.t002

E-gene sequence analysis and phylogenetic inference revealed a grouping of the ten identified WNV strains with lineage 2 sequences which widely occur across the African continent [67] (Fig 4A). Surprisingly, the detected WNV strains did not group together but fell within two different clades of WNV lineage 2. Nine of the ten WNV strains that originated from the agricultural field clustered with clade 2d, whereas strain M6646-NA-2020 that originated from rangeland grouped with clade 2b. Clade 2b has so far only been detected in South Africa in 1958 and in Cyprus in 1968 [67]. A complete genome of one strain from each clade was sequenced. Phylogenetic analyses based on complete coding sequences confirmed grouping with clades 2b and 2d and showed that M6848-NA-2020 and M6646-NA-2020 each were most closely related to viruses detected in South Africa (Fig 4B). Strain M6646-NA-2020 contained a deletion of 18 nucleotides in the 3’UTR region which was not present in any of the closely related WNV strains (Fig 5).

thumbnail
Download:
Fig 4. Phylogenetic relationship of detected WNV strains.

The phylogenetic tree is based on a 1,032 nucleotide alignment of E-protein sequences of all detected WNV strains and selected global WNV strains (A) or on an alignment of the entire open reading frame (ORF) of WNV M6848-NA-2020 and M6646-NA-2020 and selected global WNV strains (B). Branch support is given for each branch. Virus sequences identified in this study are shown in red. Reference sequences are marked with their accession number in grey and the two-letter code of their countries of origin. AU, Australia; CY, Cyprus; CZ, Czechia; ES, Spain; GR, Greece; IN, India; MG, Madagascar; RU, Russia; UG, Uganda; US, United States of America; ZA, South Africa. The phylogenies are rooted to a WNV lineage 4 strain from Russia (AY277251).

https://doi.org/10.1371/journal.pntd.0009311.g004

thumbnail
Download:
Fig 5. Schematic depiction of WNV 3’ UTR sequences.

Strains of different lineages were aligned, and insertions or deletions are shown. Numbers indicate genomic positions. Sequences identified in this study are shown in red.

https://doi.org/10.1371/journal.pntd.0009311.g005

Characterization of BAGV

DNA barcoding showed that the three BAGV-positive mosquito pools consisted of Cx. univittatus mosquitoes. Virus isolation attempts were successful for all three strains MP312-NA-2018, MP314-NA-2018, and MP370-NA-2018 in mosquito cells (C6/36). A CPE was observed in C6/36 cells, 24-, 10-, and 14-days post infection (dpi) for strains MP312-NA-2018, MP314-NA-2018, and MP370-NA-2018, respectively. Growth analyses in C6/36 cells revealed similar growth of the three strains and extremely high genome copy numbers of ca. 1x1012 RNA copies/ml five dpi (Fig 6). For further assessment of virus growth and putative host tropism, an array of different cell lines was infected with MP314-NA-2018. The mosquito cell line CXT was also permissive albeit genome copy numbers were approximately two log factors lower compared to C6/36 cells. Since BAGV is a bird pathogen, we infected three different avian cell lines. BAGV was only observed to replicate in AGE1.CR cells derived from muscovy ducks but not in LMH and DF-1 cells derived from chicken. Replication of BAGV in AGE1.CR was almost as productive as in C6/36 cells. BAGV growth was the least-well in mammalian cells with maximum genome copy numbers ca. 106-fold less than in C6/36 cells. However, similar viral RNA copy numbers were reached in non-human primate (Vero), human (HEK293T) and rodent (BHK-21) cells. No growth was observed in cattle-derived cells (KNR).

thumbnail
Download:
Fig 6. BAGV growth in vertebrate and insect cell lines.

The mosquito cell line C6/36 (derived from Aedes Stegomyia albopictus) was infected with the three BAGV strains MP312-NA-2018, MP314-NA-2018, and MP370-NA-2018 in duplicates at an MOI of 0.01. Mammalian cell lines BHK-21 (hamster), HEK293T (human), VeroE6 (monkey) and KNR (cattle) were infected in duplicates with BAGV strain MP314-NA-2018 at an MOI of 0.1. Avian cell lines AGE1.CR (duck), LMH (chicken) and DF1 (chicken) and the insect cell line CXT (derived from Culex tarsalis) were infected in duplicates with BAGV strain MP314-NA-2018 at an MOI of 0.01. A sample of cell culture supernatant was taken every 24 hours for five consecutive days and viral copy numbers were determined by RT-qPCR.

https://doi.org/10.1371/journal.pntd.0009311.g006

The genome of BAGV strain MP314-NA-2018 was sequenced from infectious cell culture supernatant. The full sequence was uploaded to GenBank under accession number MW672101. Full genome based phylogenetic tree inference showed that BAGV MP314-NA-2018 was most closely related to BAGV strain LC318701 from Zambia (Fig 7A). In further phylogenetic analyses, which included the partial NS5 sequence of a BAGV strain (MN329584) that was isolated from monal pheasants in South Africa in 2016, BAGV MP314-NA-2018 grouped with this sequence (Fig 7B). As MP314-NA-2018 is clearly distinct from LC318701 from Zambia (97.5% divergence) and from MN329584 from South Africa (99.2% divergence based on 1041 nucleotides of viral RdRP sequence), it seems to represent a strain endemic to Namibia.

thumbnail
Download:
Fig 7. Phylogenetic inference of detected BAGV strains.

Phylogenetic relationship of MP314-NA-2018, Israel turkey meningoencephalomyelitis virus (ITV) and selected global BAGV strains based on the entire RdRp gene (A). The tree is based on an alignment of 1,078 nt to include the partial NS5 sequence of a BAGV strain isolated from monal pheasants in South Africa in 2019 (B). Branch support is marked at each branch. Virus sequences detected in this work are shown in red. Reference sequences from GenBank contain the virus name, the accession number in grey and the two-letter code of their countries of origin. CF, Central African Republic; CI, Côte d’Ivoire; ES, Spain; IN, India; IL, Israel; SN, Senegal; ZA, South Africa; ZM, Zambia.

https://doi.org/10.1371/journal.pntd.0009311.g007

Discussion

Information on arbovirus circulation in vector populations in the tropics remain sparse, although surveillance in these regions is of particular importance since arboviruses are primarily endemic in tropical and subtropical regions, often causing local outbreaks there. Here, we screened mosquitoes from the Zambezi region of north-eastern Namibia for infection with flaviviruses and report the detection of the vertebrate-pathogenic arboviruses WNV and BAGV, as well as the identification of 241 strains of two previously known and seven putative novel insect-specific flaviviruses. This is the first detection of WNV and BAGV in insect vectors in Namibia and also the first detection of WNV in Namibia in more than 30 years.

We detected WNV at the end of the yearly rainy season in March 2020 but not at the end of 2018 nor during or at the end of the rainy season in 2019. Heavy precipitation is a key driver for an upsurge in mosquito populations and WNV infections, especially in dryer landscapes where mosquitoes only find breeding places after heavy rainfall and in the wet season [68]. Such heavier rainfalls may have contributed to the detection of WNV in 2020. With 253.9 mm precipitation in the Zambezi region was notably higher between January and March 2020 compared to the 98.3 mm and 76.1 mm recorded for the same time period in 2018 and 2019, respectively [69].

WNV is transmitted by a large variety of mosquito species [70]. Mosquitoes of the Culex pipiens L complex are the primary WNV vectors in North America and Europe. In contrast, Cx. univittatus has been identified as the most important WNV vector in South Africa [71]. We exclusively detected WNV in female Cx. univittatus mosquitoes, confirming the importance of this species as WNV vector in southern Africa. Of note, WNV was only found in the Mashi conservancy in 2020 but not a single WNV infection was detected in Mudumu National Park and Wuparo Conservancy in 2018 where the absolute numbers of collected Cx. univittatus mosquitoes peaked (Fig 2). Such a discrete temporal and spatial distribution pattern of WNV has been observed in long term surveys in the USA [72].

We observed a WNV prevalence of 0.1% (10/10,206 mosquitoes). Other studies conducted in South Africa between 1956 to 1980 found WNV infection rates in mosquitoes of less than 0.026% [71]. However, infection rates in Cx. univittatus can reach up to 39.0/1000 or 9.6/1000 during WNV epizootics in the same region [71]. Notably, all the ten detected WNV strains fell into two different lineage 2 clades, with one strain grouping with clade 2b and nine strains grouping with 2d. The 2b strain (M6646-NA-2020) originated from rangeland and the 2d strains from agricultural fields in the Mashi conservancy that were about 9 km apart from each other, further emphasizing that WNV is endemic and diversified in Namibia. Identity values of 95.9–99.9% in pairwise alignment of partial E-protein fragments of the nine 2d sequences support a genetically diversified set of strains. This becomes particularly evident when comparing the observed genetic divergence rate of maximal 4% to that found in other studies. For example, nucleotide divergences in mosquitoes were found to range between 1–2% in the Montréal area of Canada from 2004–2016 [73], and below 1% in the US, 0.41%–0.72% in Texas [74], 0.26%–0.5% in Chicago [75], and 0.3% in Connecticut [76]. Higher divergence rates of 0.16–3% were found in mosquitoes from Israel [77]. Overall, the detection of WNV strains from two clades and the observed genetic nucleotide divergence of about 4% in 10 mosquitoes from one area at a single time point demonstrate the circulation of a genetically diverse WNV population in the Zambezi region of Namibia.

We could further show that the mosquito specimen that was infected with M6848-NA-2020 took a blood-meal on a human host indicating human exposure to WNV. We could not identify any vertebrate DNA remaining from blood feeding in the other nine WNV-positive mosquitoes. Although Cx. univittatus is primarily ornithophilic, it has been shown to also feed on humans and was found in human dwellings [71]. Maintenance of WNV in its enzootic cycle in southern Africa involves ornithophilic mosquitoes and a wide range of southern African avian species [71]. Taken together, the high WNV detection rate and genetic diversity, as well as human host sources indicate the public health relevance of WNV in Namibia. However, there is a general lack of knowledge on the ecology and epidemiology of WNV in Africa including the study region of this paper. Other knowledge gaps include unknown distribution and incidence patterns of WNV in Africa, as well as on species involved in its natural amplification cycles. Future studies should focus on endemic transmission cycles and geographic distribution of WNV, as well as human infection rates with WNV.

To our knowledge, so far only two other clade 2b strains have been detected. One was found in a bird in South Africa in 1958 and the other in a bird in Cyprus in 1968 [67]. We discovered an 18-nucleotide deletion in the 3’ UTR of strain M6646-NA-2020 that was not present in the other strains of clade 2b. The 3’UTR region of WNV is known to harbor a variety of sequence elements that play a role in maintenance of correct viral RNA secondary structures, interaction with cellular proteins and protecting the viral genome from degradation by RNases [78,79]. These highly conserved, functional regions are located in the distal part of flavivirus UTRs, close to the genome 3’-end. However, the downstream part of the flaviviral 3’ UTR that closely follows the 3’-end of the CDS, is well known for its high variability and typically contains insertions and deletions [80].

Unfortunately, we were not successful in isolating any of the WNV strains in cell culture although genome copies of up to 6.8 x 107 per mosquito were found. This is most probably due to thawing of the mosquito samples during their transport to the laboratory in Germany as the dryshipper was damaged on the flight. In contrast, all detected BAGV strains were successfully isolated in cell culture from mosquitoes. These samples were transported previously, and no damage was observed.

Three BAGV strains were uncovered in Cx. univittatus mosquitoes from Mudumu National Park in the Zambezi region. BAGV has previously been isolated from several species of the Culex genus [32,34,35] including Cx. univittatus [33]. BAGV is maintained in nature by transmission between wild birds of the Phasianidae family and Culex spp. [81]. Susceptible species also include farmed game birds like turkeys and wild game birds like pheasants and partridges [26–28,31]. Several members of the Phasianidae family, like Peliperdix coqui, Scleroptila gutturalis, Dendroperdix sephaena, Pternistis adspersus, and Pternistis swainsonii are endemic in north-eastern Namibia [82]. However, little is known on BAGV maintenance in the Zambezi region. Screening of potential natural bird reservoir species for antibodies against BAGV could provide further insight into endemic transmission cycles, geographic distribution, and virus prevalence. In this study, growth analyses revealed great differences between cells derived from hosts involved in the natural virus transmission cycle, mosquito and bird cell lines, and mammalian cell lines derived from primates, rodents, cattle, and humans. Viral genome copy numbers were on average 100,000 times higher in mosquito and bird cells than in mammalian cells, except for cattle cells which were not susceptible. The mammalian cell lines BHK-21 and VeroE6 have been reported to support BAGV replication [30]. We found that the human cell line HEK293T enabled BAGV replication to a similar extent than BHK-21 and VeroE6 cells. Possibly BAGV is capable of infecting humans as neutralizing antibodies were found in encephalitic patients from India [36]. Such susceptibility of human cells for BAGV infection would indicate that humans may not be refractory for BAGV infection. We further found that muscovy duck cells (Cairina moschata, Anatidae) support BAGV replication to similar levels as mosquito cells, whereas chicken-derived cells were not susceptible. Although muscovy ducks are endemic to Central America, these findings may imply that BAGV may infect a wider range of bird species than partridges and pheasants. In Spain common wood pigeons (Columba palumbus, Columbidae) can also be infected by BAGV [83]. Thus, testing a wide variety of different bird species of the Zambezi region will help to identify the susceptible species and provide insight into virus maintenance in nature.

The high detection rate of insect-specific flaviviruses (95%, 241 of 254 pools positive) in the sample set confirmed findings from other studies, providing further evidence that ISFs are widely spread in mosquito populations worldwide, now including southern Africa [84–87]. Co-infection studies have shown that infection with insect-specific viruses can alter replication and transmission of arboviruses in mosquitoes [88]. Notably, an inhibition of WNV replication mediated by a variety of ISFs in mosquitoes has been reported [89–92]. In contrast, Culex flavivirus, that was also found in mosquitoes of this study (Fig 3), seems to be positively associated with WNV infections in mosquitoes [93]. We have not observed any co-infections with multiple flaviviruses in our samples. However, in vitro co-infection experiments with arboviruses and ISFs detected in this study could identify if these viruses could have an effect on vector competence or transmission rates.

In conclusion, we demonstrated that a wide genetic diversity of flaviviruses is circulating in mosquitoes in the Zambezi region of north-eastern Namibia. Importantly, we identified two clusters of WNV belonging to two different clades and BAGV. WNV and BAGV can cause outbreaks, including severe disease and mortality in humans and birds, respectively. Further studies are needed to assess their health impact, as well as their social and economic impact in southern Africa.

Supporting information

S1 Table. Primers used for BAGV genome amplification.

https://doi.org/10.1371/journal.pntd.0009311.s001

(DOCX)

Acknowledgments

We thank the Namibian Ministry of Environment, Forestry and Tourism (MET) and the National Commission on Research Science and Technology (NCRST) for granting permission to conduct this research. We are grateful to the staff of Mudumu and Bwabwata National Parks and Wuparo and Mashi Conservancies for their support, especially Mrs. Ellen Simataa, Mr. Matambo Singwangwa, Mrs. Manga Sangwali and Mr. Nicky Pudeho. We thank Marco Marklewitz for help with organizing field work and Anne Kopp for assistance in the laboratory.

References

  1. 1. Pierson TC, Diamond MS. The continued threat of emerging flaviviruses. Nat Microbiol. 2020;5(6):796–812. Epub 2020/05/06. pmid:32367055; PubMed Central PMCID: PMC7696730.
  2. 2. Campos GS, Bandeira AC, Sardi SI. Zika Virus Outbreak, Bahia, Brazil. Emerg Infect Dis. 2015;21(10):1885–6. Epub 2015/09/25. pmid:26401719; PubMed Central PMCID: PMC4593454.
  3. 3. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–7. Epub 2013/04/09. pmid:23563266; PubMed Central PMCID: PMC3651993.
  4. 4. Mensah E, Zowalaty MEE. Arboviruses in South Africa, known and unknown. 2018;13(11):787–802.
  5. 5. Smithburn KC, Hughes TP, Burke AW, Paul JH. A Neurotropic Virus Isolated from the Blood of a Native of Uganda1. 1940;s1-20(4):471–92. https://doi.org/10.4269/ajtmh.1940.s1-20.471.
  6. 6. Le Guenno B, Bougermouh A, Azzam T, Bouakaz R. West Nile: a deadly virus? Lancet. 1996;348(9037):1315. Epub 1996/11/09. pmid:8909403.
  7. 7. Platonov AE, Shipulin GA, Shipulina OY, Tyutyunnik EN, Frolochkina TI, Lanciotti RS, et al. Outbreak of West Nile virus infection, Volgograd Region, Russia, 1999. Emerg Infect Dis. 2001;7(1):128–32. Epub 2001/03/27. pmid:11266303; PubMed Central PMCID: PMC2631674.
  8. 8. Murgue B, Murri S, Zientara S, Durand B, Durand JP, Zeller H. West Nile outbreak in horses in southern France, 2000: the return after 35 years. Emerg Infect Dis. 2001;7(4):692–6. Epub 2001/10/05. pmid:11585534; PubMed Central PMCID: PMC2631744.
  9. 9. Weinberger M, Pitlik SD, Gandacu D, Lang R, Nassar F, Ben David D, et al. West Nile fever outbreak, Israel, 2000: epidemiologic aspects. Emerg Infect Dis. 2001;7(4):686–91. Epub 2001/10/05. pmid:11585533; PubMed Central PMCID: PMC2631774.
  10. 10. May FJ, Davis CT, Tesh RB, Barrett AD. Phylogeography of West Nile virus: from the cradle of evolution in Africa to Eurasia, Australia, and the Americas. J Virol. 2011;85(6):2964–74. Epub 2010/12/17. pmid:21159871; PubMed Central PMCID: PMC3067944.
  11. 11. Watson JT, Pertel PE, Jones RC, Siston AM, Paul WS, Austin CC, et al. Clinical characteristics and functional outcomes of West Nile Fever. Ann Intern Med. 2004;141(5):360–5. Epub 2004/09/09. pmid:15353427.
  12. 12. Chancey C, Grinev A, Volkova E, Rios M. The global ecology and epidemiology of West Nile virus. Biomed Res Int. 2015;2015:376230. Epub 2015/04/14. pmid:25866777; PubMed Central PMCID: PMC4383390.
  13. 13. Hayes EB, Sejvar JJ, Zaki SR, Lanciotti RS, Bode AV, Campbell GL. Virology, pathology, and clinical manifestations of West Nile virus disease. Emerg Infect Dis. 2005;11(8):1174–9. Epub 2005/08/17. pmid:16102303; PubMed Central PMCID: PMC3320472.
  14. 14. Soverow JE, Wellenius GA, Fisman DN, Mittleman MA. Infectious disease in a warming world: how weather influenced West Nile virus in the United States (2001–2005). Environmental health perspectives. 2009;117(7):1049–52. Epub 2009/08/06. pmid:19654911; PubMed Central PMCID: PMC2717128.
  15. 15. McIntosh B, Jupp P, Dos Santos I, Meenehan GJSAJoS. Epidemics of West Nile and Sindbis viruses in South Africa with Culex (Culex) univittatus Theobald as vector. 1976;72(10):295–300.
  16. 16. Venter M, Human S, Zaayman D, Gerdes GH, Williams J, Steyl J, et al. Lineage 2 west nile virus as cause of fatal neurologic disease in horses, South Africa. Emerg Infect Dis. 2009;15(6):877–84. Epub 2009/06/16. pmid:19523285; PubMed Central PMCID: PMC2727306.
  17. 17. Venter M, Swanepoel R. West Nile virus lineage 2 as a cause of zoonotic neurological disease in humans and horses in southern Africa. Vector Borne Zoonotic Dis. 2010;10(7):659–64. Epub 2010/09/22. pmid:20854018.
  18. 18. Mathengtheng L, Burt FJ. Use of envelope domain III protein for detection and differentiation of flaviviruses in the Free State Province, South Africa. Vector Borne Zoonotic Dis. 2014;14(4):261–71. Epub 2014/04/03. pmid:24689894.
  19. 19. Steyn J, Botha E, Stivaktas VI, Buss P, Beechler BR, Myburgh JG, et al. West Nile Virus in Wildlife and Nonequine Domestic Animals, South Africa, 2010–2018. Emerg Infect Dis. 2019;25(12):2290–4. Epub 2019/11/20. pmid:31742510; PubMed Central PMCID: PMC6874268.
  20. 20. Weyer J, Thomas J, Leman PA, Grobbelaar AA, Kemp A, Paweska JT. Human cases of Wesselsbron disease, South Africa 2010–2011. Vector Borne Zoonotic Dis. 2013;13(5):330–6. Epub 2013/03/12. pmid:23473219.
  21. 21. Oymans J, van Keulen L, Wichgers Schreur PJ, Kortekaas J. Early Pathogenesis of Wesselsbron Disease in Pregnant Ewes. Pathogens. 2020;9(5). Epub 2020/05/18. pmid:32414152; PubMed Central PMCID: PMC7281025.
  22. 22. Diagne MM, Faye M, Faye O, Sow A, Balique F, Sembene M, et al. Emergence of Wesselsbron virus among black rat and humans in Eastern Senegal in 2013. One Health. 2017;3:23–8. Epub 2017/06/16. pmid:28616499; PubMed Central PMCID: PMC5454166.
  23. 23. Blackburn NK, Swanepoel R. An investigation of flavivirus infections of cattle in Zimbabwe Rhodesia with particular reference to Wesselsbron virus. J Hyg (Lond). 1980;85(1):1–33. Epub 1980/08/01. pmid:6300226; PubMed Central PMCID: PMC2134009.
  24. 24. Coetzer JA, Barnard B. Hydrops amnii in sheep associated with hydranencephaly and arthrogryposis with wesselsbron disease and rift valley fever viruses as aetiological agents. 1977.
  25. 25. Coetzer J, Theodoridis AJOJVR. Clinical and pathological studies in adult sheep and goats experimentally infected with Wesselsbron disease virus. 1982;49:19–22.
  26. 26. Steyn J, Botha EM, Lourens C, Coetzer JAW, Venter M. Bagaza Virus in Himalayan Monal Pheasants, South Africa, 2016–2017. Emerg Infect Dis. 2019;25(12):2299–302. Epub 2019/11/20. pmid:31742517; PubMed Central PMCID: PMC6874265.
  27. 27. Komarov A, Kalmar EJVR. A hitherto undescribed disease-turkey meningo-encephalitis. 1960;72:257–61.
  28. 28. Barnard BJ, Buys SB, Du Preez JH, Greyling SP, Venter HJ. Turkey meningo-encephalitis in South Africa. Onderstepoort J Vet Res. 1980;47(2):89–94. Epub 1980/06/01. pmid:6251413.
  29. 29. Fernandez-Pinero J, Davidson I, Elizalde M, Perk S, Khinich Y, Jimenez-Clavero MA. Bagaza virus and Israel turkey meningoencephalomyelitis virus are a single virus species. J Gen Virol. 2014;95(Pt 4):883–7. Epub 2014/01/25. pmid:24457974.
  30. 30. Llorente F, Perez-Ramirez E, Fernandez-Pinero J, Elizalde M, Figuerola J, Soriguer RC, et al. Bagaza virus is pathogenic and transmitted by direct contact in experimentally infected partridges, but is not infectious in house sparrows and adult mice. Vet Res. 2015;46:93. Epub 2015/09/05. pmid:26338714; PubMed Central PMCID: PMC4559182.
  31. 31. Aguero M, Fernandez-Pinero J, Buitrago D, Sanchez A, Elizalde M, San Miguel E, et al. Bagaza virus in partridges and pheasants, Spain, 2010. Emerg Infect Dis. 2011;17(8):1498–501. Epub 2011/08/02. pmid:21801633; PubMed Central PMCID: PMC3381565.
  32. 32. Digoutte JP. Bagaza (BAG). 1978;27(2_Part_2):376–7. https://doi.org/10.4269/ajtmh.1978.27.376.
  33. 33. Gordon SW, Tammariello RF, Linthicum KJ, Dohm DJ, Digoutte JP, Calvo-Wilson MA. Arbovirus Isolations from Mosquitoes Collected During 1988 in the Senegal River Basin. 1992;47(6):742–8. ://doi.org/10.4269/ajtmh.1992.47.742.
  34. 34. Traore-Lamizana M, Zeller HG, Mondo M, Hervy JP, Adam F, Digoutte JP. Isolations of West Nile and Bagaza viruses from mosquitoes (Diptera: Culicidae) in central Senegal (Ferlo). J Med Entomol. 1994;31(6):934–8. Epub 1994/11/01. pmid:7815413.
  35. 35. Diallo M, Nabeth P, Ba K, Sall AA, Ba Y, Mondo M, et al. Mosquito vectors of the 1998–1999 outbreak of Rift Valley Fever and other arboviruses (Bagaza, Sanar, Wesselsbron and West Nile) in Mauritania and Senegal. Med Vet Entomol. 2005;19(2):119–26. Epub 2005/06/17. pmid:15958020.
  36. 36. Bondre VP, Sapkal GN, Yergolkar PN, Fulmali PV, Sankararaman V, Ayachit VM, et al. Genetic characterization of Bagaza virus (BAGV) isolated in India and evidence of anti-BAGV antibodies in sera collected from encephalitis patients. J Gen Virol. 2009;90(Pt 11):2644–9. Epub 2009/07/03. pmid:19570951.
  37. 37. Smithburn KC, Paterson HE, Heymann CS, Winter PA. An agent related to Uganda S virus from man and mosquitoes in South Africa. S Afr Med J. 1959;33:959–62. Epub 1959/11/14. pmid:13832180.
  38. 38. Williams MC, Woodall JP. An Epidemic of an Illness Resembling Dengue in the Morogoro District of Tanganyika. East Afr Med J. 1964;41:271–5. Epub 1964/06/01. pmid:14182630.
  39. 39. Braack L, Gouveia de Almeida AP, Cornel AJ, Swanepoel R, de Jager C. Mosquito-borne arboviruses of African origin: review of key viruses and vectors. Parasit Vectors. 2018;11(1):29. Epub 2018/01/11. pmid:29316963; PubMed Central PMCID: PMC5759361.
  40. 40. McIntosh BM. Mosquito-borne virus diseases of man in southern Africa. S Afr Med J. 1986;Suppl:69–72. Epub 1986/10/11. pmid:3022406.
  41. 41. Were F. The dengue situation in Africa. Paediatr Int Child Health. 2012;32 Suppl 1:18–21. Epub 2012/06/08. pmid:22668445; PubMed Central PMCID: PMC3381440.
  42. 42. Kokernot RH, Szlamp EL, Levitt J, McIntosh BM. Survey for antibodies against arthropod-borne viruses in the sera of indigenous residents of the Caprivi Strip and Bechuanaland Protectorate. Trans R Soc Trop Med Hyg. 1965;59(5):553–62. Epub 1965/09/01. pmid:5893148.
  43. 43. Joubert JJ, van der Merwe CA, Lourens JH, Lecatsas G, Siegruhn C. Serological markers of hepatitis B virus and certain other viruses in the population of eastern Caprivi, Namibia. Trans R Soc Trop Med Hyg. 1991;85(1):101–3. Epub 1991/01/01. pmid:2068734.
  44. 44. Burt FJ, Grobbelaar AA, Leman PA, Anthony FS, Gibson GV, Swanepoel R. Phylogenetic relationships of southern African West Nile virus isolates. Emerg Infect Dis. 2002;8(8):820–6. Epub 2002/07/27. pmid:12141968; PubMed Central PMCID: PMC2732512.
  45. 45. Noden BH, Musuuo M, Aku-Akai L, van der Colf B, Chipare I, Wilkinson R. Risk assessment of flavivirus transmission in Namibia. Acta Trop. 2014;137:123–9. Epub 2014/05/29. pmid:24865792.
  46. 46. Edwards FWJMotERI-CA Pupae. Mosquitoes of the Ethiopian Region. III.-Culicine adults and pupae. 1941.
  47. 47. Gillies MT, De Meillon BJTAoAsotS The Anophelinae of Africa south of the Sahara (Ethiopian zoogeographical region). 1968.
  48. 48. Jupp PG. Mosquitoes of Southern Africa: culicinae and toxorhynchitinae: Ekogilde Publishers; 1996.
    • 49. Crochu S, Cook S, Attoui H, Charrel RN, De Chesse R, Belhouchet M, et al. Sequences of flavivirus-related RNA viruses persist in DNA form integrated in the genome of Aedes spp. mosquitoes. J Gen Virol. 2004;85(Pt 7):1971–80. Epub 2004/06/26. pmid:15218182.
    • 50. Moureau G, Temmam S, Gonzalez JP, Charrel RN, Grard G, de Lamballerie X. A real-time RT-PCR method for the universal detection and identification of flaviviruses. Vector Borne Zoonotic Dis. 2007;7(4):467–77. Epub 2007/11/21. pmid:18020965.
    • 51. Eiden M, Vina-Rodriguez A, Hoffmann B, Ziegler U, Groschup MH. Two new real-time quantitative reverse transcription polymerase chain reaction assays with unique target sites for the specific and sensitive detection of lineages 1 and 2 West Nile virus strains. J Vet Diagn Invest. 2010;22(5):748–53. Epub 2010/09/03. pmid:20807934.
    • 52. 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(7):2323–30. Epub 2002/06/29. pmid:12089242; PubMed Central PMCID: PMC120575.
    • 53. Corman VM, Rasche A, Baronti C, Aldabbagh S, Cadar D, Reusken CB, et al. Assay optimization for molecular detection of Zika virus. Bull World Health Organ. 2016;94(12):880–92. Epub 2016/12/21. pmid:27994281; PubMed Central PMCID: PMC5153932.
    • 54. Domingo C, Patel P, Yillah J, Weidmann M, Mendez JA, Nakoune ER, et al. Advanced yellow fever virus genome detection in point-of-care facilities and reference laboratories. J Clin Microbiol. 2012;50(12):4054–60. Epub 2012/10/12. pmid:23052311; PubMed Central PMCID: PMC3503008.
    • 55. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22(22):4673–80. Epub 1994/11/11. pmid:7984417; PubMed Central PMCID: PMC308517.
    • 56. Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30(14):3059–66. Epub 2002/07/24. pmid:12136088; PubMed Central PMCID: PMC135756.
    • 57. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35(6):1547–9. Epub 2018/05/04. pmid:29722887; PubMed Central PMCID: PMC5967553.
    • 58. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel OJSb. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. 2010;59(3):307–21.
    • 59. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9. Epub 1994/10/01. pmid:7881515.
    • 60. Alcaide M, Rico C, Ruiz S, Soriguer R, Munoz J, Figuerola J. Disentangling vector-borne transmission networks: a universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals. PLoS One. 2009;4(9):e7092. Epub 2009/09/22. pmid:19768113; PubMed Central PMCID: PMC2740869.
    • 61. Junglen S, Kopp A, Kurth A, Pauli G, Ellerbrok H, Leendertz FH. A new flavivirus and a new vector: characterization of a novel flavivirus isolated from uranotaenia mosquitoes from a tropical rain forest. J Virol. 2009;83(9):4462–8. Epub 2009/02/20. pmid:19224998; PubMed Central PMCID: PMC2668441.
    • 62. Reed LJ, Muench HJAjoe. A simple method of estimating fifty per cent endpoints. 1938;27(3):493–7.
    • 63. Marklewitz M, Zirkel F, Rwego IB, Heidemann H, Trippner P, Kurth A, et al. Discovery of a unique novel clade of mosquito-associated bunyaviruses. J Virol. 2013;87(23):12850–65. Epub 2013/09/27. pmid:24067954; PubMed Central PMCID: PMC3838119.
    • 64. Himly M, Foster DN, Bottoli I, Iacovoni JS, Vogt PK. The DF-1 chicken fibroblast cell line: transformation induced by diverse oncogenes and cell death resulting from infection by avian leukosis viruses. Virology. 1998;248(2):295–304. Epub 1998/08/29. pmid:9721238.
    • 65. Kawaguchi T, Nomura K, Hirayama Y, Kitagawa T. Establishment and characterization of a chicken hepatocellular carcinoma cell line, LMH. Cancer Res. 1987;47(16):4460–4. Epub 1987/08/15. pmid:3607775.
    • 66. Jordan I, Vos A, Beilfuss S, Neubert A, Breul S, Sandig V. An avian cell line designed for production of highly attenuated viruses. Vaccine. 2009;27(5):748–56. Epub 2008/12/17. pmid:19071186.
    • 67. McMullen AR, Albayrak H, May FJ, Davis CT, Beasley DWC, Barrett ADT. Molecular evolution of lineage 2 West Nile virus. J Gen Virol. 2013;94(Pt 2):318–25. Epub 2012/11/09. pmid:23136360; PubMed Central PMCID: PMC3709619.
    • 68. Paz S. Climate change impacts on West Nile virus transmission in a global context. Philos Trans R Soc Lond B Biol Sci. 2015;370(1665). Epub 2015/02/18. pmid:25688020; PubMed Central PMCID: PMC4342965.
    • 69. World Weather Online [cited 2020 09/12]. Available from: http://www.worldweatheronline.com/.
      • 70. Ciota AT. West Nile virus and its vectors. Curr Opin Insect Sci. 2017;22:28–36. Epub 2017/08/15. pmid:28805636.
      • 71. Jupp PG. The ecology of West Nile virus in South Africa and the occurrence of outbreaks in humans. Ann N Y Acad Sci. 2001;951:143–52. Epub 2002/01/19. pmid:11797772.
      • 72. Dunphy BM, Kovach KB, Gehrke EJ, Field EN, Rowley WA, Bartholomay LC, et al. Long-term surveillance defines spatial and temporal patterns implicating Culex tarsalis as the primary vector of West Nile virus. Scientific reports. 2019;9(1):6637. Epub 2019/05/01. pmid:31036953; PubMed Central PMCID: PMC6488619.
      • 73. Therrien C, Fournier É, Ludwig A, Ménard J, Charest H, Martineau C. Phylogenetic analysis of West Nile virus in Quebec, Canada, 2004–2016: Co-circulation of distinct variants harbouring conserved amino acid motifs in North America. Virology. 2019;537:65–73. Epub 2019/08/30. pmid:31465892.
      • 74. Mann BR, McMullen AR, Swetnam DM, Salvato V, Reyna M, Guzman H, et al. Continued evolution of West Nile virus, Houston, Texas, USA, 2002–2012. Emerg Infect Dis. 2013;19(9):1418–27. Epub 2013/08/24. pmid:23965756; PubMed Central PMCID: PMC3810927.
      • 75. Amore G, Bertolotti L, Hamer GL, Kitron UD, Walker ED, Ruiz MO, et al. Multi-year evolutionary dynamics of West Nile virus in suburban Chicago, USA, 2005–2007. Philos Trans R Soc Lond B Biol Sci. 2010;365(1548):1871–8. Epub 2010/05/19. pmid:20478882; PubMed Central PMCID: PMC2880109.
      • 76. Armstrong PM, Vossbrinck CR, Andreadis TG, Anderson JF, Pesko KN, Newman RM, et al. Molecular evolution of West Nile virus in a northern temperate region: Connecticut, USA 1999–2008. Virology. 2011;417(1):203–10. Epub 2011/07/05. pmid:21723580; PubMed Central PMCID: PMC3152594.
      • 77. Lustig Y, Hindiyeh M, Orshan L, Weiss L, Koren R, Katz-Likvornik S, et al. Mosquito Surveillance for 15 Years Reveals High Genetic Diversity Among West Nile Viruses in Israel. The Journal of infectious diseases. 2016;213(7):1107–14. Epub 2015/11/26. pmid:26597260.
      • 78. Hussain M, Torres S, Schnettler E, Funk A, Grundhoff A, Pijlman GP, et al. West Nile virus encodes a microRNA-like small RNA in the 3’ untranslated region which up-regulates GATA4 mRNA and facilitates virus replication in mosquito cells. Nucleic Acids Res. 2012;40(5):2210–23. Epub 2011/11/15. pmid:22080551; PubMed Central PMCID: PMC3300009.
      • 79. Ng WC, Soto-Acosta R, Bradrick SS, Garcia-Blanco MA, Ooi EE. The 5’ and 3’ Untranslated Regions of the Flaviviral Genome. Viruses. 2017;9(6). Epub 2017/06/08. pmid:28587300; PubMed Central PMCID: PMC5490814.
      • 80. Gritsun DJ, Jones IM, Gould EA, Gritsun TS. Molecular archaeology of Flaviviridae untranslated regions: duplicated RNA structures in the replication enhancer of flaviviruses and pestiviruses emerged via convergent evolution. PLoS One. 2014;9(3):e92056. Epub 2014/03/22. pmid:24647143; PubMed Central PMCID: PMC3960163.
      • 81. Elizalde M, Cano-Gomez C, Llorente F, Perez-Ramirez E, Casades-Marti L, Aguilera-Sepulveda P, et al. A Duplex Quantitative Real-Time Reverse Transcription-PCR for Simultaneous Detection and Differentiation of Flaviviruses of the Japanese Encephalitis and Ntaya Serocomplexes in Birds. Front Vet Sci. 2020;7:203. Epub 2020/05/07. pmid:32373639; PubMed Central PMCID: PMC7186316.
      • 82. Brown C, Thomson N, Boorman MJBO. Checklist and analysis of the birds of Namibia as at 31 January 2016. 2017;8:20: 1–153.
      • 83. Gamino V, Gutierrez-Guzman AV, Fernandez-de-Mera IG, Ortiz JA, Duran-Martin M, de la Fuente J, et al. Natural Bagaza virus infection in game birds in southern Spain. Vet Res. 2012;43:65. Epub 2012/09/13. pmid:22966904; PubMed Central PMCID: PMC3483237.
      • 84. Kopp A, Gillespie TR, Hobelsberger D, Estrada A, Harper JM, Miller RA, et al. Provenance and geographic spread of St. Louis encephalitis virus. mBio. 2013;4(3):e00322–13. Epub 2013/06/14. pmid:23760463; PubMed Central PMCID: PMC3685209.
      • 85. Cholleti H, Hayer J, Abilio AP, Mulandane FC, Verner-Carlsson J, Falk KI, et al. Discovery of Novel Viruses in Mosquitoes from the Zambezi Valley of Mozambique. PLoS One. 2016;11(9):e0162751. Epub 2016/09/30. pmid:27682810; PubMed Central PMCID: PMC5040392.
      • 86. Colmant AMG, Hobson-Peters J, Bielefeldt-Ohmann H, van den Hurk AF, Hall-Mendelin S, Chow WK, et al. A New Clade of Insect-Specific Flaviviruses from Australian Anopheles Mosquitoes Displays Species-Specific Host Restriction. mSphere. 2017;2(4). Epub 2017/07/18. pmid:28713857; PubMed Central PMCID: PMC5506557.
      • 87. Guzman H, Contreras-Gutierrez MA, Travassos da Rosa APA, Nunes MRT, Cardoso JF, Popov VL, et al. Characterization of Three New Insect-Specific Flaviviruses: Their Relationship to the Mosquito-Borne Flavivirus Pathogens. Am J Trop Med Hyg. 2018;98(2):410–9. Epub 2017/10/11. PubMed Central PMCID: PMC5929187. pmid:29016330
      • 88. Öhlund P, Lundén H, Blomström AL. Insect-specific virus evolution and potential effects on vector competence. Virus genes. 2019;55(2):127–37. Epub 2019/01/12. pmid:30632016; PubMed Central PMCID: PMC6458977.
      • 89. Kent RJ, Crabtree MB, Miller BR. Transmission of West Nile virus by Culex quinquefasciatus say infected with Culex Flavivirus Izabal. PLoS neglected tropical diseases. 2010;4(5):e671. Epub 2010/05/11. pmid:20454569; PubMed Central PMCID: PMC2864301.
      • 90. Hobson-Peters J, Yam AW, Lu JW, Setoh YX, May FJ, Kurucz N, et al. A new insect-specific flavivirus from northern Australia suppresses replication of West Nile virus and Murray Valley encephalitis virus in co-infected mosquito cells. PLoS One. 2013;8(2):e56534. Epub 2013/03/06. pmid:23460804; PubMed Central PMCID: PMC3584062.
      • 91. Kenney JL, Solberg OD, Langevin SA, Brault AC. Characterization of a novel insect-specific flavivirus from Brazil: potential for inhibition of infection of arthropod cells with medically important flaviviruses. J Gen Virol. 2014;95(Pt 12):2796–808. Epub 2014/08/26. pmid:25146007; PubMed Central PMCID: PMC4582674.
      • 92. Colmant AMG, Hall-Mendelin S, Ritchie SA, Bielefeldt-Ohmann H, Harrison JJ, Newton ND, et al. The recently identified flavivirus Bamaga virus is transmitted horizontally by Culex mosquitoes and interferes with West Nile virus replication in vitro and transmission in vivo. PLoS neglected tropical diseases. 2018;12(10):e0006886. Epub 2018/10/26. pmid:30356234; PubMed Central PMCID: PMC6200184.
      • 93. Newman CM, Cerutti F, Anderson TK, Hamer GL, Walker ED, Kitron UD, et al. Culex flavivirus and West Nile virus mosquito coinfection and positive ecological association in Chicago, United States. Vector borne and zoonotic diseases (Larchmont, NY). 2011;11(8):1099–105. Epub 2011/01/22. pmid:21254845.
      • 94. Kalvelage L, Revilla Diez J, Bollig M. How much remains? Local value capture from tourism in Zambezi, Namibia. Tourism Geographies. 2020:1–22.