martes, 19 de noviembre de 2019

Bagaza Virus in Himalayan Monal Pheasants, South Africa, 2016–2017 - Volume 25, Number 12—December 2019 - Emerging Infectious Diseases journal - CDC

Bagaza Virus in Himalayan Monal Pheasants, South Africa, 2016–2017 - Volume 25, Number 12—December 2019 - Emerging Infectious Diseases journal - CDC

Issue Cover for Volume 25, Number 12—December 2019



Volume 25, Number 12—December 2019
Dispatch

Bagaza Virus in Himalayan Monal Pheasants, South Africa, 2016–2017

Jumari Steyn1, Elizabeth M. Botha1, Carina Lourens, Jacobus A.W. Coetzer, and Marietjie VenterComments to Author 
Author affiliations: University of Pretoria, Pretoria, South Africa

Abstract

Bagaza virus (BAGV) has not been reported in birds in South Africa since 1978. We used phylogenetic analysis and electron microscopy to identify BAGV as the likely etiology in neurologic disease and death in Himalayan monal pheasants in Pretoria, South Africa. Our results suggest circulation of BAGV in South Africa.
The flavivirus genus of family Flaviviridae consists of 53 virus species, including arboviruses of medical and veterinary relevance, such as West Nile virus and Bagaza virus (BAGV). BAGV was isolated in 1966 from Culex mosquitoes in the Bagaza district of Central African Republic (1). In 1978, BAGV was isolated from turkeys with clinical signs similar to Israel turkey meningoencephalitis virus (ITV) in South Africa (2). BAGV infection causes neurologic disease in avian species, especially turkeys and other members of the Phasianidae family; 1 report suggests that BAGV and ITV are the same viral species (3).
BAGV also has been detected in various mosquito species in western Africa (4,5), India (6), and the Arabian Peninsula (7) and in wild partridges in Spain (8). No evidence of the virus has been reported in other parts of Africa. Zoonotic transmission was reported in India after patients with acute encephalitis demonstrated 15% positivity for BAGV neutralizing antibodies (6). We report detection of BAGV in fatalities in Himalayan monal pheasants in South Africa during 2016–2017.

The Study

In April 2016, two Himalayan monal pheasants (Lophophorus impejanus) and 1 tragopan pheasant (Tragopan melanocephalus) suddenly died on a property northeast of Pretoria, Gauteng Province, South Africa. In June 2017, this property had another 4 monal pheasants that displayed signs of lethargy and ataxia and died within a day. Around the same time, a residence in the northern suburbs of Pretoria had 5 monal and 2 tragopan pheasants that exhibited neurologic signs and died. That residence had another incidence in 2018 when a monal pheasant exhibited neurologic disease. Also in 2018, a monal pheasant was found dead in North West Province, South Africa. The cause of these deaths was unknown. All the birds were adults that were locally bred from parents imported from Belgium >2 years before.
Brain tissue from the 16 birds was sent to the Department of Veterinary Tropical Diseases (DVTD), University of Pretoria (Pretoria, South Africa), for virus isolation and to the Centre for Viral Zoonoses (CVZ), University of Pretoria, for zoonotic arbovirus investigations. At the CVZ, we extracted RNA from the brain tissues by using the RNeasy Mini Kit (QIAGEN, https://www.qiagen.comExternal Link) according to manufacturer’s instructions under Biosafety Level 3 conditions. We used nested real-time reverse transcription Pan-Flavi assay targeting the nonstructural coding gene 5 (NS5) (9) to identify the etiologic agent (10,11). To obtain a larger NS5 gene segment, we performed additional PCR using SuperScript III/Platinum Taq Mix (Invitrogen, https://www.thermofisher.comExternal Link) and the MAMD (9) and FLAVI-2 (10) primers with the following cycling conditions: 50°C for 30 min; 94°C for 15 min; 35 cycles of 94°C for 45 s, 50°C for 45 s, 72°C for 1 min; and 72°C for 10 min. We successfully obtained a larger NS5 gene segment for phylogenetic analyses, but only for 4 positive birds.
We assembled and edited sequence data by using CLC Main WorkBench (https://www.qiagenbioinformatics.comExternal Link) and performed multiple sequence alignments using the online version of MAFFT (http://mafft.cbrc.jp/alignment/server/index.htmlExternal Link) with default parameters. We used MEGA 6.06 (https://www.megasoftware.netExternal Link) to view, edit, and truncate the datasets. We downloaded reference sequences for the flavivirus genus from GenBank (12). We conducted maximum likelihood analysis in RAxML (13), invoking the autoMRE bootstopping function applying a general time-reversible plus gamma model with default 4 rate categories on both datasets. We performed an analysis on the longer dataset by using BEAST version 1.8 (http://beast.community) and a relaxed log-normal clock, general time-reversible plus gamma model, and default priors to generate a maximum clade credibility tree (MCC). We ran a Markov chain Monte Carlo analysis for 106 generations, saving every 1,000th tree. We estimated effective sample size by using Tracer version 1.6 (http://tree.bio.ed.ac.uk/software/tracerExternal Link) with an effective sample size value >200. We used TreeAnnotator version 1.8 (http://beast.community) to generate the MCC tree and discarded 15% as burn-in. We displayed bootstrap and posterior probabilities on RAxML topology.
We performed virus isolation on all PCR-positive samples. We inoculated brain tissue supernatant onto a confluent monolayer of baby hamster kidney fibroblast cells (BHK-21 line) in 25 cm2 tissue culture flasks and incubated at 37°C for 1 h. Then we added Dulbecco’s Minimum Essential Medium (ThermoFisher Scientific, https://www.thermofisher.comExternal Link) containing 2% fetal bovine serum and 0.1 mg/mL gentamycin. We harvested cells and supernatant when 80% of the cell monolayer showed cytopathic effect and sent these to the Electron Microscopy Unit of The University of Pretoria and to the CVZ for molecular identification.
Thumbnail of Maximum-likelihood phylogram of BAGV isolated in samples from Himalayan monal pheasants (black dots), South Africa, 2016–2017. Phylogram represents partial (1,079 nt) nonstructural coding gene 5 (NS5; taxa = 30). Bootstrap support with values >60 indicated on branches with posterior probabilities >0.95 from a maximum clade credibility tree. BAGV strains from this study are available in GenBank under the following accession nos.: ZRU349/17/6, no. MN329586; ZRU350/17/1, no. MN32
Figure 1. Maximum-likelihood phylogram of BAGV isolated in samples from Himalayan monal pheasants (black dots), South Africa, 2016–2017. Phylogram represents partial (1,079 nt) nonstructural coding gene 5 (NS5; taxa = 30). Bootstrap support with values...
The Pan-Flavi assay targeting the NS5 gene resulted in amplicons of the expected size in 8/13 (61.5%, 95% CI 35.1%–88.0%) Himalayan monal pheasants but not in the 3 tragopan pheasants tested. Neurologic signs were reported before death in 7/8 (87.5%, 95% CI 0.5%–66.5%) positive birds, but 1/8 (12.5%, 95% CI 0–2.1%) was found dead and its clinical signs are unknown. Inqaba Biotec (https://www.inqababiotec.co.zaExternal Link) performed Sanger sequencing; we confirmed all positive samples as BAGV using phylogenetic analysis of the flavivirus genus NS5 PCR regions at CVZ. Phylogenetic analyses were based on partial nucleotide sequences of NS5 from genomic positions 9091–9280 (166 nt) and 9030–10109 (1,079 nt) and were used to compare the identified strains with other flaviviruses. Analyses confirmed the molecular results from all 8 Himalayan monals as BAGV in the Ntaya virus group with a bootstrap value of 92, sister to BAGV strains from Spain (bootstrap value 68) (Appendix Figure). The 4 strains for which we amplified a larger region (1,079 nt) formed 2 well-supported sister groups, both with a bootstrap value of 100 and phylogenetic probability of 1, with nucleotide similarities of 97.7%–99.7%, and highest nucleotide identity (96.7%–97.7%) to strain Zambia_Zmq13mz26 (GenBank accession no. LC318701.1) isolated from a mosquito (Figure 1).
Thumbnail of Electron microscopy of Bagaza virus isolated in samples from Himalayan monal pheasants, South Africa, 2016–2017. A) Circles indicate occasional particles with size range and approximate morphology of Flaviviridae observed in samples ZRU350_17_1 and ZRU350_17_2. Scale bars indicate 200 nm. B) A few isolated fringed isometric particles of 40–65 nm (top row) and free-lying smooth-surfaced particles of 25–40 nm (bottom row) of suspected Flaviviridae observed in sample ZRU349_17_6. C) A
Figure 2. Electron microscopy of Bagaza virus isolated in samples from Himalayan monal pheasants, South Africa, 2016–2017. A) Circles indicate occasional particles with size range and approximate morphology of Flaviviridae observed in samples...
Electron microscopy on 3 BAGV cultures (sample nos. ZRU350_17_1, ZRU350_17_2, and ZRU349_17_6) (Appendix Table) from 2017 confirmed the presence of Flaviviridae particles (Figure 2). We observed fringed isometric and free-lying smooth-surfaced particles typical of Flaviviridae (Figure 2, panels B and C).

Conclusions

We detected BAGV in the offspring of monal pheasants imported from Belgium to South Africa. We sequenced BAGV strains and found they monophyletically clustered with strains from Spain rather than strains from West Africa. However, nucleotide similarities in the large gene segment were highest when compared with a strain from Zambia that was isolated from a Cx. quinquefasciatus mosquito (GenBank accession no. LC318701.1; 14), an endemic species in South Africa that could be a BAGV vector. We noted 2 distinct monophyletic clusters of BAGV, a cluster composed of strains from West Africa and older strains and a cluster containing the newly sequenced birds with BAGV from Spain and more recent strains that could indicate several circulating strains or genotypes.
We used virus isolation and electron microscopy results to confirm the etiology of the agent as a flavivirus. The causative link between the clinical symptoms of the monal pheasants and evidence of BAGV infection should be regarded with caution because we did not exclude other possible infectious and noninfectious etiologies. However, detection of BAGV in the brain suggests crossing of the blood–brain barrier and exclusion of other flaviviruses, arboviruses, and orthobunyaviruses suggests BAGV as a probable cause. Future work will focus on next-generation sequencing to obtain full genomes because initial attempts were unsuccessful. More data are needed to determine the endemicity of BAGV and the reservoir host and vectors of BAGV in South Africa and to define the seroprevelance of these infections in birds and possibly in humans.
Ms. Steyn is a virologist and PhD candidate at the Centre for Viral Zoonoses at the University of Pretoria, Pretoria, South Africa. Her primary research focuses on investigating arboviruses with zoonotic potential at human–animal interface areas. Mrs. Botha was a master’s student and research assistant in the Centre for Viral Zoonoses at the University of Pretoria, and currently is employed at a private pathology laboratory. Her primary interest is flaviviruses.

Acknowledgments

We would like to attribute this paper to the memory of Chris Kingsley, who submitted some of the cases to the project, and recognize him for his work in bird conservation. Dr. Kingsley sadly passed away before the Bagaza virus findings could be published. We would also like to acknowledge Louwtjie Snyman for his help with the phylogenetic analyses.
This study was cleared by section 20 (12/11/1/1) approval through the Department of Agriculture, Forestry and Fisheries (clearance no. V057-15) and by the animal ethics committee (clearance no. H12-16) of the University of Pretoria (UP) and the PhD research committee. The work was funded through UP Zoonotic Arbo- and Respiratory Virus Program income-generated funds. J. Steyn received doctoral scholarships from the National Research Foundation (grant no. 95175), the Meat Industry Trust (grant no. IT8114/98) and the Poliomyelitis Research Foundation (grant no. 15/112), as well as the US Centers for Disease Control and Prevention cooperative agreement with the University of Pretoria (no. 5 NU2GGH001874-02-00).

References

  1. Digoutte  JPBagaza (BAG) strain: Dak Ar B 209. Am J Trop Med Hyg1978;27:3767DOIExternal Link
  2. Barnard  BJBuys  SBDu Preez  JHGreyling  SPVenter  HJTurkey meningo-encephalitis in South Africa. Onderstepoort J Vet Res1980;47:8994.PubMedExternal Link
  3. Fernández-Pinero  JDavidson  IElizalde  MPerk  SKhinich  YJiménez-Clavero  MABagaza virus and Israel turkey meningoencephalomyelitis virus are a single virus species. J Gen Virol2014;95:8837DOIExternal LinkPubMedExternal Link
  4. Diallo  MNabeth  PBa  KSall  AABa  YMondo  Met 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 Entomol2005;19:11926DOIExternal LinkPubMedExternal Link
  5. Traore-Lamizana  MZeller  HGMondo  MHervy  JPAdam  FDigoutte  JPIsolations of West Nile and Bagaza viruses from mosquitoes (Diptera: Culicidae) in central Senegal (Ferlo). J Med Entomol1994;31:9348DOIExternal LinkPubMedExternal Link
  6. Bondre  VPSapkal  GNYergolkar  PNFulmali  PVSankararaman  VAyachit  VMet al. Genetic characterization of Bagaza virus (BAGV) isolated in India and evidence of anti-BAGV antibodies in sera collected from encephalitis patients. J Gen Virol2009;90:26449DOIExternal LinkPubMedExternal Link
  7. Camp  JVKaruvantevida  NChouhna  HSafi  EShah  JNNowotny  NMosquito biodiversity and mosquito-borne viruses in the United Arab Emirates. Parasit Vectors2019;12:153DOIExternal LinkPubMedExternal Link
  8. Agüero  MFernández-Pinero  JBuitrago  DSánchez  AElizalde  MSan Miguel  Eet al. Bagaza virus in partridges and pheasants, Spain, 2010. Emerg Infect Dis2011;17:1498501DOIExternal LinkPubMedExternal Link
  9. Zaayman  DHuman  SVenter  MA highly sensitive method for the detection and genotyping of West Nile virus by real-time PCR. J Virol Methods2009;157:15560DOIExternal LinkPubMedExternal Link
  10. Scaramozzino  NCrance  JMJouan  ADeBriel  DAStoll  FGarin  DComparison of flavivirus universal primer pairs and development of a rapid, highly sensitive heminested reverse transcription-PCR assay for detection of flaviviruses targeted to a conserved region of the NS5 gene sequences. J Clin Microbiol2001;39:19227DOIExternal LinkPubMedExternal Link
  11. Ayers  MAdachi  DJohnson  GAndonova  MDrebot  MTellier  RA single tube RT-PCR assay for the detection of mosquito-borne flaviviruses. J Virol Methods2006;135:2359DOIExternal LinkPubMedExternal Link
  12. Benson  DAKarsch-Mizrachi  ILipman  DJOstell  JSayers  EW. GenBank. Nucleic Acids Res. 2010;38(suppl_1):D46–51.
  13. Stamatakis  ARAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics2014;30:13123DOIExternal LinkPubMedExternal Link
  14. Orba  YHang’ombe  BMMweene  ASSasaki  MEshita  YSawa  H. Bagaza virus genomic RNA, nearly complete genome, strain: Zmq13mz26. GenBank;2017.
Figures
Cite This Article

DOI: 10.3201/eid2512.190756
Original Publication Date: 11/6/2019
1These authors contributed equally to this article.

No hay comentarios:

Publicar un comentario