Next Article in Journal
Emerging Human Babesiosis with “Ground Zero” in North America
Next Article in Special Issue
Development of a Well-Characterized Rhesus Macaque Model of Ebola Virus Disease for Support of Product Development
Previous Article in Journal
Characterization of the Habitat- and Season-Independent Increase in Fungal Biomass Induced by the Invasive Giant Goldenrod and Its Impact on the Fungivorous Nematode Community
Previous Article in Special Issue
Taï Forest Virus Does Not Cause Lethal Disease in Ferrets
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 2
10.3390/microorganisms9020438
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogenetic Characterization of Crimean-Congo Hemorrhagic Fever Virus Detected in African Blue Ticks Feeding on Cattle in a Ugandan Abattoir

by
Eddie M. Wampande
1,
Peter Waiswa
1,
David J. Allen
2,
Roger Hewson
2,3,
Simon D. W. Frost
4,5 and
Samuel C. B. Stubbs
4,6,*
1
College of Veterinary Medicine, Animal Resources and Biosecurity, Makerere University, Kampala, Uganda
2
Department of Infection Biology, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK
3
Virology and Pathogenesis Group, Public Health England, Porton Down, Salisbury SP4 0JG, UK
4
Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK
5
Microsoft Research, Redmond, Washington, DC 98052, USA
6
Department of Veterinary Medicine, University of Cambridge, Cambridge CB3 0ES, UK
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(2), 438; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9020438
Submission received: 27 January 2021 / Revised: 10 February 2021 / Accepted: 15 February 2021 / Published: 20 February 2021
(This article belongs to the Special Issue Hemorrhagic Fever Viruses: Pathogenesis and Countermeasures)
Browse Figures
Review Reports Versions Notes

Abstract

:
Crimean-Congo hemorrhagic fever virus (CCHFV) is the most geographically widespread of the tick-borne viruses. However, African strains of CCHFV are poorly represented in sequence databases. In addition, almost all sequence data collected to date have been obtained from cases of human disease, while information regarding the circulation of the virus in tick and animal reservoirs is severely lacking. Here, we characterize the complete coding region of a novel CCHFV strain, detected in African blue ticks (Rhipicephalus (Boophilus) decoloratus) feeding on cattle in an abattoir in Kampala, Uganda. These cattle originated from a farm in Mbarara, a major cattle-trading hub for much of Uganda. Phylogenetic analysis indicates that the newly sequenced strain belongs to the African genotype II clade, which predominantly contains the sequences of strains isolated from West Africa in the 1950s, and South Africa in the 1980s. Whilst the viral S (nucleoprotein) and L (RNA polymerase) genome segments shared >90% nucleotide similarity with previously reported genotype II strains, the glycoprotein-coding M segment shared only 80% nucleotide similarity with the next most closely related strains, which were derived from ticks in Western India and Northern China. This genome segment also displayed a large number of non-synonymous mutations previously unreported in the genotype II strains. Characterization of this novel strain adds to our limited understanding of the natural diversity of CCHFV circulating in both ticks and in Africa. Such data can be used to inform the design of vaccines and diagnostics, as well as studies exploring the epidemiology and evolution of the virus for the establishment of future CCHFV control strategies.

1. Introduction

Crimean-Congo hemorrhagic fever virus (CCHFV) is a zoonotic pathogen primarily transmitted between vertebrate hosts through the bite of infected ticks, as well as through direct contact with infected livestock, ticks, and bodily fluids [1]. CCHFV is the most geographically widespread of the tick-borne viruses; endemic throughout Africa, Asia, Eastern Europe, and the Middle East, it circulates in an enzootic cycle between animals and ticks, sporadically emerging to cause severe outbreaks of human disease [2]. Outbreak fatality rates of up to 30% have been reported, although CCHF symptoms can range in severity, from mild febrile illness to hemorrhagic fever and multi-organ failure [2].
CCHFV is classified within the genus Orthonairovirus (family Nairoviridae) and carries a tripartite RNA genome of negative polarity. The three genome segments are labelled small (S), medium (M), and large (L), according to their respective lengths [3]. The ~1.7 kb S segment encodes the viral nucleoprotein (NP) and, in the reverse orientation, an overlapping non-structural (NSs) protein. The M segment is approximately 5.3 kb and encodes a single polyprotein precursor, which is post-translationally cleaved into two envelope glycoproteins (Gn and Gc), several secreted glycoproteins (GP38, GP85, and GP160), and the non-structural M protein, NSm. The pre-Gn precursor also contains a hyper-variable region known as the mucin-like domain (MLD). The ~12.1 kb L segment encodes the viral RNA-dependent RNA polymerase.
Strains of CCHFV display a high degree of genetic variability for an arthropod-borne virus, with nucleotide diversity as high as 20% for the S segment and 31% for the M segment [4]. Genotypes are assigned to each genomic segment independently based on phylogenetic analysis, with each genotype largely defined by a distinct geographical region: three are African (genotypes I, II, and III), two Asian (genotypes IVa and IVb), and two European (genotypes V and VI). However, limited sampling and uneven representation of CCHFV sequences has hindered phylogeographic and evolutionary analyses [2,5]. This is particularly true for African strains, which are poorly represented in sequence databases in comparison to strains from Europe, Asia, and the Middle East. Furthermore, the majority of CCHFV strains reported to date have been derived from cases of severe human disease, while insufficient sampling from tick vectors and reservoir hosts is almost certainly limiting our knowledge of the natural diversity of CCHFV and the potential for spillover of unanticipated viral variants.
Recent detection of an African genotype III strain following an outbreak in Spain [6], and the presence of African genotype II M segments in viruses obtained from pools of Hyalomma spp. ticks in India and China [7,8] has highlighted the need for a more complete genomic dataset in order to better understand the epidemiology of the virus, including viral transmission routes and the extent of re-assortment among genotypes. The potential for such re-assortment is particularly important, as many molecular tests for CCHFV are only capable of detecting a restricted number of closely related genotypes [9]. It also remains to be seen whether recently reported vaccine candidates, including one based on the S and M segments of the Hoti strain (genotype V) [10] and another based on the M segment of the IbAr102000 strain (genotype III) [11], can elicit a protective immune response against viruses belonging to other more distantly related genotypes.

2. Materials and Methods

Seventy ticks were collected from cattle, goats, and a pig in an abattoir in Kampala, Uganda between September and December 2019 (Table S1) as part of a pilot study into the virome of ticks in Uganda. Ticks were immediately frozen in liquid nitrogen following collection and stored at −80 °C until ready for processing. Individual ticks were morphologically identified using taxonomic keys [12] before being crushed using a sterile micropestle and resuspended in 250 µL of phosphate-buffered saline (PBS). Tissue debris was pelleted by low-speed centrifugation (400× g) and total nucleic acid (TNA) was extracted from 200 µL of the supernatant using the High Pure Viral Nucleic Acid extraction kit (Roche, Basel, Switzerland) following the manufacturer’s instructions. The extracted nucleic acid was then lyophilized before shipping at ambient temperature to the University of Cambridge, UK for further analysis.
Upon arrival in the UK, TNA was reconstituted in 20 µL of RNase-free water (RFW) supplemented with 2 U of RNaseOUT (Thermo Fisher Scientific, Waltham, MA, USA). Two microliters of reconstituted TNA was subjected to reverse transcription in 20 µL reactions using Superscript III reverse transcriptase (Thermo Fisher Scientific) and primed with 50 ng random hexamers (Thermo Fisher Scientific), as per the enzyme manufacturer’s instructions. Two microliters from each reverse transcription reaction was used as direct input to screen for the presence of nairovirus RNA by (RT-)PCR. PCR was performed in 20 µL reactions, using Phusion DNA polymerase (New England Biolabs, Ipswich, MA, USA) and 0.5 μM pan-nairovirus primers targeting the S segment [13], set up according to the enzyme manufacturer’s instructions. PCR cycling conditions were as follows: incubation at 98 °C for 2 min, followed by 45 cycles of 98 °C for 5 s, 60 °C for 20 s, and 72 °C for 20 s. A final extension was performed at 72 °C for a further 5 min, followed by incubation at 4 °C. The resulting PCR products were visualized by gel electrophoresis, and amplicons of the expected size (~400 bp) were individually purified and confirmed to be CCHFV by Sanger dideoxy sequencing.
Five RT-PCR-positive ticks from Mbarara were selected at random for further characterization by Illumina metagenomic sequencing. Five microliters of TNA from each tick were pooled for a total volume of 25 μL before being treated with 2 U of TURBO DNase (Thermo Fisher Scientific), and incubating for 20 min at 37 °C. DNase-treated RNA was purified using RNA Clean & Concentrator-5 spin columns (Zymo Research, Irvine, CA, USA), eluted in 15 µL of RFW and quantified using the Qubit High Sensitivity RNA assay (Thermo Fisher Scientific ). Illumina sequencing libraries were prepared from 200 ng of purified RNA using the Zymo-Seq RiboFree Total RNA Library kit (Zymo Research), which incorporates a ribosomal RNA depletion step. The resulting library was sequenced on the Illumina NextSeq 500 platform, using the 300-cycle mid-output kit (v2.5, Illumina San Diego, CA, USA). Resulting paired-end reads were imported into CLC Genomics Workbench v7.5.1 (Qiagen, Hilden, Germany) and sequences sharing 95% sequence identity over 35% of their length were assembled de novo. Contiguous sequences corresponding to the CCHFV S, M, and L segments were identified using BLASTn against the redundant nucleotide sequence database with default settings [14].
The CCHFV genome segments (S, M, and L) were independently characterized through phylogenetic analysis. MAFFT v7.427 [15] was used to align the open reading frame of each genome segment against a representative set of CCHFV sequences from GenBank, including members of each genotype (African, Asian, European, etc.) and all available sequences belonging to the genotype to which the Mbarara strain was found to belong (genotype II). For the S and L segments, complete open reading frame (ORF) sequences were analyzed. For the M segment, an 857 bp region at the 5′ end corresponding to the MLD was excluded from phylogenetic analysis due to the extreme variability of the region, which produced sub-optimal alignments. Inclusion of this region did not significantly alter phylogenetic clustering, but bootstrap values were adversely affected. Following trimming of the M segment alignment, the final datasets were 1449, 4323, and 11,838 nt long and included 50, 56, and 46 sequences for S, M, and L segments, respectively.
Maximum likelihood phylogenies were calculated using IQTree v1.6.11 [16], employing 1000 iterations of the ultrafast bootstrap estimation [17]. Optimal substitution models for each alignment were selected using the ModelFinder algorithm [18] according to Bayesian information criterion scores (GTR+F+I+G4 for S and L segments and GTR+F+R3 for the M segment). Nucleotide identity matrices were generated for each segment using CLUSTAL Omega [19] and intra-genotype nucleotide sequence similarity plots were generated using EMBOSS v6.6.0 [20] Plotcon function across a window size l = 4. Midpoint rooted, maximum likelihood phylogenies and sequence similarity plots were visualized in R v3.6.3 [21] using the ggplot2 [22] and ggtree [23] packages. Genotype II M segment sequences were also analyzed for evidence of recombination using the RDP5 tool [24] and for evidence of episodic diversifying selection using the adaptive branch-site random effects likelihood (aBRSEL) method [25]. Nucleotide alignments of genotype II sequences were also translated into predicted amino acid alignments in Seaview v4.7 [26] and polymorphic and informative sites were identified using the ape v5.0 package [27] in R v3.6.3.

3. Results

A total of 70 ticks were collected from cattle (n = 48), goats (n = 21), and a pig (n = 1) (Table S1). The ticks were morphologically identified as Rhipicephalus appendiculatus (n = 16), Rhipicephalus (Boophilus) decoloratus (n = 40), and Amblyomma variegatum (n = 14). Nairovirus RNA was detected by RT-PCR in 8 of the 70 ticks (11.4%). These 8 ticks were all identified as the African blue tick R. decoloratus, and all were picked from cattle. Seven of the CCHFV-positive ticks were sampled from cattle originating from a farm in Mbarara, Western Uganda, and one from a farm in Nakasongola in Central Uganda, as determined based on cattle movement permits. Sanger dideoxy sequencing of the RT-PCR products confirmed the presence of CCHFV and revealed all 8 of the ~400 bp fragments of the S segment to be identical.
Metagenomic sequencing was performed on a pool of five randomly selected (RT-)PCR-positive R. decoloratus ticks picked from Mbarara cattle, producing 27,745,476 150 bp paired end reads. A contiguous sequence corresponding to the S segment was 1624 nt in length, with an open reading frame (ORF) of 1449 nt (482 AA), the M segment was 5341 nt (ORF length: 5102 nt; 1700 AA), and the L segment was 11,838 nt (ORF length: 11,838 nt; 3945 AA). The average coverage depth was 6189×, 3882×, and 1833× for S, M, and L, respectively. The full contiguous sequences have been deposited in NCBI GenBank under accession numbers MW452933-MW452935.
Phylogenetic analyses revealed that all three segments of the Mbarara CCHFV strain clustered within the African genotype II (Africa 2) clade, together with strains from the Democratic Republic of Congo, Uganda, and South Africa (Figure 1, Figure 2 and Figure 3). The S and L genome segments were highly similar to previously reported genotype II strains, sharing between 91.79–94.34% nucleotide identity (98.76–99.38% amino acid identity) and 91.57–94.12% nucleotide identity (97.49–98.56% amino acid identity), respectively (Figure S1).
In comparison to the S and L segments, the M segment of the Mbarara strain displayed reduced similarity to genotype II sequences, sharing between 77.54% and 79.62% nucleotide identity (80.45–83.40% amino acid identity). Those sequences with the greatest similarity to the Mbarara strain belonged to several strains from Western India (Gujarat in 2016 and Rajasthan in 2019) and Northern China (Xinjiang in 1979, 2004 and 2016) (Figure S1). As expected, sequence conservation of the M segment was lowest (~40% similarity) in the 5′ region that encompasses the signal peptide and mucin-like domains (Figure 4). However, the variability in genotype II sequences was not limited to the MLD alone; the regions coding for the secreted glycoprotein GP38 and the surface glycoproteins Gn and Gc also displayed a notable degree of intra-genotype nucleotide variability (80–85% similarity). Finally, analysis of the genotype II M segment sequences for recombination and diversifying episodic selection revealed no evidence for either (p > 0.05).
Alignment of genotype II predicted amino acid sequences revealed 164 residues unique to the Mbarara strain (singletons) (Table 1). These novel residues were located throughout all three genomic segments, at a frequency ranging from average to above average for the genotype. A majority of singletons (148/164) were observed in the M segment, in which Mbarara strain displayed the greatest number of unique residues across every gene, including the MLD, the non-structural NSm, and the glycoproteins GP38, Gn, and Gc.

4. Discussion

Here we describe the detection and characterization of a genotype II strain of CCHFV, present in African blue ticks feeding on cattle from a farm in Mbarara, Uganda. To our knowledge, CCHFV has only rarely been described in this species of tick (R. decoloratus): twice in Senegal (1972 and 1975) and twice in Uganda (1978 and 2015) [1,28]. It was not possible to determine whether the virus was derived from within the tick or the host bloodmeal as the specimen was homogenized during the nucleic acid extraction process. However, the tick appears to be the more likely source, as a limited number of experimental studies have found CCHFV viremia to be transient and short-lived in cattle [29]. As a single-host tick species, R. decoloratus is unlikely to play a major role in zoonotic transmission of the virus. However, it is possible that these ticks may act as a reservoir for the virus, maintaining infection of the animal host and facilitating the infection of vector species (e.g., Hyalomma spp.) through co-feeding [30]. Additionally, several two- and three-host Rhipicephalus species have also been shown to support the presence of CCHFV [30]. Therefore, given that Rhipicephalus ticks are common throughout Africa, further investigation into their role in CCHFV epizootiology and epidemiology is clearly warranted.
Previous reports of high seroprevalence in livestock and animal handlers throughout Africa (e.g., Bukbuk et al. 2016 [31]), suggest that CCHFV-infected ticks and/or animals can often be found in abattoirs, highlighting the potential for focusing surveillance and control measures around these facilities. Notably, the majority of CCHFV-infected ticks sampled from the abattoir were picked from cattle originating from Mbarara. This district, situated in the cattle corridor of Western Uganda, is a major source of beef and dairy products for the country. Indeed, cattle herds in Nakasongola have been largely sourced from Mbarara, following their depletion during rebel occupation. This may account for the surprisingly high degree of similarity between the Mbarara CCHFV strain and the strain detected in a tick picked from a cow from Nakasongola, and suggests that further sampling of tick CCHFV strains from other regions of Uganda is required in order to better examine transmission dynamics within the country. Possible alternative explanations for the similarity between the CCHFV strains detected in Mbarara and Nakasongola ticks include host-switching of the tick upon arrival at the abattoir, or that the CCHFV-infected tick population originated from the abattoir itself.
Two of the Mbarara strain genomic segments—the nucleocapsid (S) and RNA polymerase (L)—showed a high degree of nucleotide conservation, sharing >90% similarity with other African genotype II sequences, even though these were derived from strains isolated several decades ago. A similar observation was made for viruses from Pakistan and Nigeria [31], suggesting a tendency for these genomic segments to remain genetically stable over time. In contrast, sequence diversity was greatest in the surface-glycoprotein-coding (M) segment, which shared only 78% nucleotide similarity with strains previously reported from Africa. Instead, the M segment sequences bearing the greatest (~80%) nucleotide similarity to the Mbarara strain belonged to two strains sequenced directly from Hyalomma spp. ticks in Western India (Gujarat in 2016 and Rajasthan in 2019), as well as several strains from Northern China (Xinjiang). The Chinese strains, one isolated from a rodent (Euchoreutes naso) in 1979, and two from Hyalomma asiaticum ticks in 2004 and 2016, were first isolated and then passaged in suckling mice prior to sequencing. This should be noted, as the virus may have acquired mutations during adaptation to its mammalian host and therefore may not be representative of the virus where it was initially isolated. Interestingly, the M segments of both the Indian and Chinese strains clustered within the African genotype II clade, however their corresponding S and L segments belonged to the Asian genotype IV, indicating that the M segment had likely been obtained through re-assortment with an African virus. Such a process is made possible through intercontinental movement of CCHFV-infected ticks on livestock and migratory birds [32]. It is unclear whether the relatively low homology of the Mbarara strain M segment to previously reported genotype II strains is due to its novel lineage or simply the result of several decades of evolution. Testing this theory is made difficult by the paucity of closely related sequences. However, the stark dissimilarity to previously reported genotype II M segment sequences, in contrast to the S and L segments, implies recombination or re-assortment with an unidentified viral lineage. Unfortunately, many methods designed to identify instances of recombination and re-assortment are reliant upon the inclusion of the parent lineages in the analysis, which is not possible with the current small sample of genotype II strains.
Amino acid residues unique to the genotype II clade were detected throughout the Mbarara strain genome, although much of the observed diversity, at both nucleotide and amino acid levels, was in the genomic M segment. Unsurprisingly, a high proportion of this variability could be attributed to the region encoding the MLD, a highly glycosylated protein of unproven function but hypothesized to be involved in shielding viral epitopes against the host immune response, based on the role of a similar domain in filoviruses [33]. This hypothesis is supported by a recent epitope-mapping study, which reported a high degree of immune reactivity by CCHFV survivors against the MLD [34]. However, it should be noted that intra-genotype diversity was not exclusively limited to this region of known hyper-variability. A high degree of nucleotide and amino acid variation was also observed in the genes encoding the viral glycoproteins, particularly the secreted glycoprotein GP38. The function of GP38 is similarly unknown, however it is believed to act as a chaperone, assisting in PreGn folding [3]. Interestingly, GP38 has recently shown promise as an antigen following the use of anti-GP38 monoclonal antibodies to protect mice against lethal challenge by CCHFV [35]. Knowledge of the spectrum of GP38 variants in animal and vector reservoirs may therefore prove useful for the design of future vaccines.

5. Conclusions

A severely neglected pathogen, the lack of genomic data from African CCHFV makes exploration of the evolutionary history of this novel strain difficult. In addition, a lack of viral data from animal and tick reservoirs has highlighted the need to better characterize the natural variation of CCHFV strains in circulation outside of severe human cases of disease. Indeed, to our knowledge only one other CCHFV strain has been reported from African ticks (IbAr102000). This genotype III virus was isolated in 1966 from a Hyalomma tick taken from a camel in Nigeria. Notably, the virus was passaged multiple times over several decades prior to sequencing, likely introducing a number of mutations throughout this period. It is therefore uncertain how representative this sequence is of the original virus, particularly a tick virus adapting to mammalian cell culture. In conclusion, the generation of additional CCHFV genome sequence directly from ticks, such as that presented here, is vital in obtaining a more complete understanding of CCHFV ecology, transmission, and evolution—all of which are necessary for the effective design and deployment of future vaccines and molecular diagnostic tests.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2076-2607/9/2/438/s1. Figure S1: Pairwise nucleotide identity comparisons for CCHFV genome segments S, M, and L; Table S1: List of collected tick specimens.

Author Contributions

Conceptualization, E.M.W. and S.C.B.S.; Data curation S.C.B.S.; Formal analysis, S.C.B.S.; Funding acquisition, E.M.W., D.J.A., R.H., S.D.W.F., and S.C.B.S.; Investigation, E.M.W., P.W., and S.C.B.S.; Methodology, E.M.W., P.W., and S.C.B.S.; Project administration, S.C.B.S.; Resources, D.J.A. and S.D.W.F.; Visualization, S.C.B.S.; Writing—original draft, S.C.B.S.; Writing—review and editing, E.M.W., D.J.A., R.H., and S.D.W.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the ALBORADA Trust, grant number RG G101892, awarded to S.C.B.S. and E.M.W.. E.M.W. is supported by the Government of the Republic of Uganda through Makerere University Research and Innovations Fund. S.C.B.S. is supported in part by an unrestricted gift from Microsoft Research. D.J.A. and R.H. are supported by the Medical Research Council, grant number MR/T029196/1.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Science and Technology Review Board of the Uganda National Council (code NS 673 on 31st July 2019).

Informed Consent Statement

Not applicable.

Data Availability Statement

All consensus sequences generated and analyzed for this study were uploaded to NCBI GenBank under accession numbers MW452933–MW452935. Raw sequence data have been deposited in ENA under study PRJEB42607. Nucleotide and amino acid multiple sequence alignments and maximum likelihood phylogenies for each segment are available at: DOI: 10.5281/zenodo.4455021.

Acknowledgments

We thank. Robert Aruho, wildlife Veterinary officer, and Dorothy Kirumira, Monitoring and Research LMNP, Uganda Wildlife Authority; the contribution of Rose Ademun Okurut, Commissioner, Animal Health Ministry of Agriculture, Animal Industry and Fisheries (MAAIF); the area Veterinarian (Grace Assiimwe) of Kiruhura District (part of original greater Mbarara); the cattle traders and the staff of the Uganda Meat Packers’ Abattoir; and especially Joseph Amagoro, for the collection of the ticks used in this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Hoogstraal, H. Review Article 1: The Epidemiology of Tick-Borne Crimean-Congo Hemorrhagic Fever in Asia, Europe, and Africa23. J. Med Èntomol. 1979, 15, 307–417. [Google Scholar] [CrossRef] [PubMed]
  2. Bente, D.A.; Forrester, N.L.; Watts, D.M.; McAuley, A.J.; Whitehouse, C.A.; Bray, M. Crimean-Congo hemorrhagic fever: History, epidemiology, pathogenesis, clinical syndrome and genetic diversity. Antivir. Res. 2013, 100, 159–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Zivcec, M.; Scholte, F.E.M.; Spiropoulou, C.F.; Spengler, J.R.; Bergeron, É. Molecular Insights into Crimean-Congo Hemorrhagic Fever Virus. Viruses 2016, 8, 106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Deyde, V.M.; Khristova, M.L.; Rollin, P.E.; Ksiazek, T.G.; Nichol, S.T. Crimean-Congo Hemorrhagic Fever Virus Genomics and Global Diversity. J. Virol. 2006, 80, 8834–8842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lukashev, A.N.; Klimentov, A.S.; Smirnova, S.E.; Dzagurova, T.K.; Drexler, J.F.; Gmyl, A.P. Phylogeography of Crimean Congo Hem-orrhagic Fever Virus. PLoS ONE 2016, 11, e0166744. [Google Scholar] [CrossRef] [PubMed]
  6. Ramírez de Arellano, E.; Hernández, L.; Goyanes, M.J.; Arsuaga, M.; Fernández, C.A.; Negredo, A.; Paz Sánchez-Seco, M. Phylo-genetic characterization of Crimean-Congo hemorrhagic fever virus, Spain. Emerg. Infect. Dis. 2017, 23, 2078–2080. [Google Scholar] [CrossRef] [Green Version]
  7. Sahay, R.R.; Dhandore, S.; Yadav, P.D.; Chauhan, A.; Bhatt, L.; Garg, V.; Gupta, N.; Nyayanit, D.A.; Shete, A.M.; Singh, R.; et al. Detection of African genotype in Hyalomma tick pools during Crimean Congo hemorrhagic fever outbreak, Rajasthan, India. Virus Res. 2020, 286, 198046. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Shen, S.; Fang, Y.; Liu, J.; Su, Z.; Liang, J.; Zhang, Z.; Wu, Q.; Wang, C.; Abudurexiti, A.; et al. Isolation, Characterization, and Phylogenetic Analysis of Two New Crimean-Congo Hemorrhagic Fever Virus Strains from the Northern Region of Xinjiang Province, China. Virol. Sin. 2018, 33, 74–86. [Google Scholar] [CrossRef]
  9. Gruber, C.E.M.; Bartolini, B.; Castilletti, C.; Mirazimi, A.; Hewson, R.; Christova, I.; Avšič, T.; Grunow, R.; Papa, A.; Sánchez-Seco, M.P.; et al. Geographical Variability Affects CCHFV Detection by RT–PCR: A Tool for In-Silico Evaluation of Molecular Assays. Viruses 2019, 11, 953. [Google Scholar] [CrossRef] [Green Version]
  10. Hawman, D.W.; Ahlén, G.; Appelberg, K.S.; Meade-White, K.; Hanley, P.W.; Scott, D.; Monteil, V.; Devignot, S.; Okumura, A.; Weber, F.; et al. A DNA-based vaccine protects against Crimean-Congo haemorrhagic fever virus disease in a Cynomolgus macaque model. Nat. Microbiol. 2021, 6, 187–195. [Google Scholar] [CrossRef] [PubMed]
  11. Buttigieg, K.R.; Dowall, S.D.; Findlay-Wilson, S.; Miloszewska, A.; Rayner, E.; Hewson, R.; Carroll, M.W. A Novel Vaccine against Crimean-Congo Haemorrhagic Fever Protects 100% of Animals against Lethal Challenge in a Mouse Model. PLoS ONE 2014, 9, e91516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Walker, A.R.; Bouattour, A.; Camicas, J.-L.; Estrada-Peña, A.; Horak, I.G.; Latif, A.A.; Pegram, R.G.; Preston, P.M. Ticks of Domestic Animals in Africa: A Guide to Identification of Species, 1st ed.; Bioscience Reports: Edinburgh, UK, 2003. [Google Scholar]
  13. Lambert, A.J.; Lanciotti, R.S. Consensus Amplification and Novel Multiplex Sequencing Method for S Segment Species Identification of 47 Viruses of the Orthobunyavirus, Phlebovirus, and Nairovirus Genera of the Family Bunyaviridae. J. Clin. Microbiol. 2009, 47, 2398–2404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  15. 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, 3059–3066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  17. Minh, B.Q.; Nguyen, M.A.T.; Von Haeseler, A. Ultrafast Approximation for Phylogenetic Bootstrap. Mol. Biol. Evol. 2013, 30, 1188–1195. [Google Scholar] [CrossRef] [PubMed]
  18. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phy-logenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  19. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [Green Version]
  20. Rice, P.; Longden, I.; Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16, 276–277. [Google Scholar] [CrossRef]
  21. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Compu-ting: Vienna, Austria, 2020. [Google Scholar]
  22. Wickham, H. Ggplot2: Elegrant Graphics for Data Analysis; Springer: New York, NY, USA, 2009. [Google Scholar]
  23. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T. ggtree: An r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  24. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murrell, B. RDP5: A computer program for analysing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. 2020. [Google Scholar] [CrossRef]
  25. Smith, M.D.; Wertheim, J.O.; Weaver, S.; Murrell, B.; Scheffler, K.; Pond, S.L.K. Less Is More: An Adaptive Branch-Site Random Effects Model for Efficient Detection of Episodic Diversifying Selection. Mol. Biol. Evol. 2015, 32, 1342–1353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Gouy, M.; Guindon, S.; Gascuel, O. SeaView Version 4: A Multiplatform Graphical User Interface for Sequence Alignment and Phylogenetic Tree Building. Mol. Biol. Evol. 2009, 27, 221–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Paradis, E.; Schliep, K. ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformation 2018, 35, 526–528. [Google Scholar] [CrossRef]
  28. Balinandi, S.; Patel, K.; Ojwang, J.; Kyondo, J.; Mulei, S.; Tumusiime, A.; Lubwama, B.; Nyakarahuka, L.; Klena, J.D.; Lutwama, J.; et al. Investigation of an isolated case of human Crimean–Congo hemorrhagic fever in Central Uganda, 2015. Int. J. Infect. Dis. 2018, 68, 88–93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Causey, O.R.; Kemp, G.E.; Madbouly, M.H.; David-West, T.S. Congo Virus from Domestic Livestock, African Hedgehog, and Arthropods in Nigeria. Am. J. Trop. Med. Hyg. 1970, 19, 846–850. [Google Scholar] [CrossRef]
  30. Gargili, A.; Estrada-Peña, A.; Spengler, J.R.; Lukashev, A.; Nuttall, P.A.; Bente, D.A. The role of ticks in the maintenance and trans-mission of Crimean-Congo hemorrhagic fever virus: A review of published field and laboratory studies. Antiviral Res. 2017, 144, 93–119. [Google Scholar] [CrossRef]
  31. Bukbuk, D.N.; Dowall, S.D.; Lewandowski, K.; Bosworth, A.; Baba, S.S.; Varghese, A.; Watson, R.J.; Bell, A.; Atkinson, B.; Hewson, R. Se-rological and Virological Evidence of Crimean-Congo Haemorrhagic Fever Virus Circulation in the Human Population of Borno State, Northeastern Nigeria. PLoS Negl. Trop. Dis. 2016, 10, e0005126. [Google Scholar] [CrossRef] [Green Version]
  32. Lindeborg, M.; Barboutis, C.; Ehrenborg, C.; Fransson, T.; Jaenson, T.G.; Lindgren, P.-E.; Lundkvist, Å.; Nyström, F.; Salaneck, E.; Waldenström, J.; et al. Migratory Birds, Ticks, and Crimean-Congo Hemorrhagic Fever Virus. Emerg. Infect. Dis. 2012, 18, 2095–2097. [Google Scholar] [CrossRef]
  33. Fusco, M.L.; Hashiguchi, T.; Cassan, R.; Biggins, J.E.; Murin, C.D.; Warfield, K.L.; Li, S.; Holtsberg, F.W.; Shulenin, S.; Vu, H.; et al. Protective mAbs and Cross-Reactive mAbs Raised by Immunization with Engineered Marburg Virus GPs. PLoS Pathog. 2015, 11, e1005016. [Google Scholar]
  34. Fritzen, A.; Risinger, C.; Korukluoglu, G.; Christova, I.; Hitzeroth, A.C.; Viljoen, N.; Burt, F.J.; Mirazimi, A.; Blixt, O. Epitope-mapping of the glycoprotein from Crimean-Congo hemorrhagic fever virus using a microarray approach. PLoS Neglected Trop. Dis. 2018, 12, e0006598. [Google Scholar] [CrossRef] [PubMed]
  35. Golden, J.W.; Shoemaker, C.J.; Lindquist, M.E.; Zeng, X.; Daye, S.P.; Williams, J.A.; Liu, J.; Coffin, K.M.; Olschner, S.; Flusin, O.; et al. GP38-targeting monoclonal antibodies protect adult mice against lethal Crimean-Congo hemorrhagic fever virus infection. Sci. Adv. 2019, 5, eaaw9535. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Microorganisms 09 00438 g001 550
Figure 1. Phylogenetic analysis of the Crimean-Congo hemorrhagic fever virus (CCHFV) S genome segment. Maximum likelihood phylogenies were estimated from complete open reading frame (ORF) nucleotide sequences where available. GenBank accession numbers, strain, country of origin, and year of sampling are listed in taxon names. Numbers at the tree nodes represent bootstrap support values estimated using the ultrafast bootstrap approximation. The Mbarara strain is highlighted in red with a triangular node tip.
Figure 1. Phylogenetic analysis of the Crimean-Congo hemorrhagic fever virus (CCHFV) S genome segment. Maximum likelihood phylogenies were estimated from complete open reading frame (ORF) nucleotide sequences where available. GenBank accession numbers, strain, country of origin, and year of sampling are listed in taxon names. Numbers at the tree nodes represent bootstrap support values estimated using the ultrafast bootstrap approximation. The Mbarara strain is highlighted in red with a triangular node tip.
Microorganisms 09 00438 g001
Microorganisms 09 00438 g002 550
Figure 2. Phylogenetic analysis of the CCHFV M genome segment. Maximum likelihood phylogenies were calculated from nucleotide sequences. The hyper-variable (857 bp) mucin-like domain was trimmed from the 5′ end prior to analysis to facilitate accurate alignment, particularly between genotypes. GenBank accession numbers, strain, country of origin, and year of sampling are listed in taxon names. Numbers at the tree nodes represent bootstrap support values estimated using the ultrafast bootstrap approximation. The Mbarara strain is highlighted in red with a triangular node tip.
Figure 2. Phylogenetic analysis of the CCHFV M genome segment. Maximum likelihood phylogenies were calculated from nucleotide sequences. The hyper-variable (857 bp) mucin-like domain was trimmed from the 5′ end prior to analysis to facilitate accurate alignment, particularly between genotypes. GenBank accession numbers, strain, country of origin, and year of sampling are listed in taxon names. Numbers at the tree nodes represent bootstrap support values estimated using the ultrafast bootstrap approximation. The Mbarara strain is highlighted in red with a triangular node tip.
Microorganisms 09 00438 g002
Microorganisms 09 00438 g003 550
Figure 3. Phylogenetic analysis of the CCHFV L genome segment. Maximum likelihood phylogenies were estimated from complete ORF nucleotide sequences where available. GenBank accession numbers, strain, country of origin, and year of sampling are listed in taxon names. Numbers at the nodes represent bootstrap support values estimated using the ultrafast bootstrap approximation. The Mbarara strain is highlighted in red with a triangular node tip.
Figure 3. Phylogenetic analysis of the CCHFV L genome segment. Maximum likelihood phylogenies were estimated from complete ORF nucleotide sequences where available. GenBank accession numbers, strain, country of origin, and year of sampling are listed in taxon names. Numbers at the nodes represent bootstrap support values estimated using the ultrafast bootstrap approximation. The Mbarara strain is highlighted in red with a triangular node tip.
Microorganisms 09 00438 g003
Microorganisms 09 00438 g004 550
Figure 4. Nucleotide conservation between the glycoprotein precursor (GPC) ORFs of M segments of genotype II strains. Nucleotide sequences of genotype II M ORFs (n = 12) were aligned, and the average nucleotide similarity was calculated using a sliding window size of 4. Variability over the context of all mature proteins that are processed from GPC is highlighted.
Figure 4. Nucleotide conservation between the glycoprotein precursor (GPC) ORFs of M segments of genotype II strains. Nucleotide sequences of genotype II M ORFs (n = 12) were aligned, and the average nucleotide similarity was calculated using a sliding window size of 4. Variability over the context of all mature proteins that are processed from GPC is highlighted.
Microorganisms 09 00438 g004
Table
Table 1. Predicted amino acid differences in the coding region of the Mbarara strain of CCHFV in comparison to previously sequenced African genotype II strains.
Table 1. Predicted amino acid differences in the coding region of the Mbarara strain of CCHFV in comparison to previously sequenced African genotype II strains.
Genome SegmentGeneNumber of Sequences (Full Length)Gene Length (AA)Polymorphic Sites * (Informative Sites †)Median Singleton ‡ Count (Range) §Mbarara Strain Singletons% Sites Unique to Mbarara Strain
SNP and Ns8 (5)48211 (3)1 (0–2)10.2
MSignal peptide10 (10)2016 (13)0 (0–1)525.0
MLD10 (10)242206 (167)4 (1–29)7631.4
GP3810 (10)27180 (55)7 (4–14)228.1
Gn11 (10)32258 (41)1 (0–4)123.7
NSm11 (11)19741 (27)1.5 (0–3)105.1
Gc12 (10)643103 (53)2 (0–12)233.6
LRdRp7 (5)3945153 (32)7 (0–60)320.8
S, M, and LAll-6144668 (391)-1642.7
* Sites with at least two alleles. † Sites with at least two alleles, where each of those alleles is found in at least two sequences. ‡ Sites displaying amino acid residues unique to the genotype II clade. § Genotype II strains, excluding Mbarara strain.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wampande, E.M.; Waiswa, P.; Allen, D.J.; Hewson, R.; Frost, S.D.W.; Stubbs, S.C.B. Phylogenetic Characterization of Crimean-Congo Hemorrhagic Fever Virus Detected in African Blue Ticks Feeding on Cattle in a Ugandan Abattoir. Microorganisms 2021, 9, 438. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9020438

AMA Style

Wampande EM, Waiswa P, Allen DJ, Hewson R, Frost SDW, Stubbs SCB. Phylogenetic Characterization of Crimean-Congo Hemorrhagic Fever Virus Detected in African Blue Ticks Feeding on Cattle in a Ugandan Abattoir. Microorganisms. 2021; 9(2):438. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9020438

Chicago/Turabian Style

Wampande, Eddie M., Peter Waiswa, David J. Allen, Roger Hewson, Simon D. W. Frost, and Samuel C. B. Stubbs. 2021. "Phylogenetic Characterization of Crimean-Congo Hemorrhagic Fever Virus Detected in African Blue Ticks Feeding on Cattle in a Ugandan Abattoir" Microorganisms 9, no. 2: 438. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9020438

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop