Next Article in Journal
Hemagglutinin Gene Variation Rate of H9N2 Avian Influenza Virus by Vaccine Intervention in China
Previous Article in Journal
Robust AAV Genotyping Based on Genetic Distances in Rep Gene That Are Maintained by Ubiquitous Recombination
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 5
10.3390/v14051041
viruses-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

Jingmen Tick Virus in Ticks from Kenya

by
Edwin O. Ogola
1,2,
Anne Kopp
3,4,
Armanda D. S. Bastos
2,
Inga Slothouwer
3,4,
Marco Marklewitz
3,4,
Dorcus Omoga
1,
Gilbert Rotich
1,
Caroline Getugi
1,
Rosemary Sang
1,5,
Baldwyn Torto
1,2,
Sandra Junglen
3,4,*,† and
David P. Tchouassi
1,*,†
1
International Centre of Insect Physiology and Ecology, Nairobi P.O. Box 30772-00100, Kenya
2
Department of Zoology and Entomology, University of Pretoria, Private Bag 20, Pretoria 0028, South Africa
3
Berlin Institute of Health, Institute of Virology, Charité—Universitätsmedizin Berlin, Corporate Member of Free University Berlin, Humboldt-University Berlin, Chariteplatz 1, 10117 Berlin, Germany
4
German Centre for Infection Research (DZIF), Associated Partner Site Charité, 10117 Berlin, Germany
5
Kenya Medical Research Institute (KEMRI), Off Raila Odinga Way, Nairobi P.O. Box 54840-00200, Kenya
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this article.
Viruses 2022, 14(5), 1041; https://doi.org/10.3390/v14051041
Submission received: 15 April 2022 / Revised: 9 May 2022 / Accepted: 10 May 2022 / Published: 13 May 2022
(This article belongs to the Section Invertebrate Viruses)
Browse Figures
Review Reports Versions Notes

Abstract

:
Jingmen tick virus (JMTV) is an arbovirus with a multisegmented genome related to those of unsegmented flaviviruses. The virus first described in Rhipicephalus microplus ticks collected in Jingmen city (Hubei Province, China) in 2010 is associated with febrile illness in humans. Since then, the geographic range has expanded to include Trinidad and Tobago, Brazil, and Uganda. However, the ecology of JMTV remains poorly described in Africa. We screened adult ticks (n = 4550, 718 pools) for JMTV infection by reverse transcription polymerase chain reaction (RT-PCR). Ticks were collected from cattle (n = 859, 18.88%), goats (n = 2070, 45.49%), sheep (n = 1574, 34.59%), and free-ranging tortoises (Leopard tortoise, Stigmochelys pardalis) (n = 47, 1.03%) in two Kenyan pastoralist-dominated areas (Baringo and Kajiado counties) with a history of undiagnosed febrile human illness. Surprisingly, ticks collected from goats (0.3%, 95% confidence interval (CI) 0.1–0.5), sheep (1.8%, 95% CI 1.2–2.5), and tortoise (74.5%, 95% CI 60.9–85.4, were found infected with JMTV, but ticks collected from cattle were all negative. JMTV ribonucleic acid (RNA) was also detected in blood from tortoises (66.7%, 95% CI 16.1–97.7). Intragenetic distance of JMTV sequences originating from tortoise-associated ticks was greater than that of sheep-associated ticks. Phylogenetic analyses of seven complete-coding genome sequences generated from tortoise-associated ticks formed a monophyletic clade within JMTV strains from other countries. In summary, our findings confirm the circulation of JMTV in ticks in Kenya. Further epidemiological surveys are needed to assess the potential public health impact of JMTV in Kenya.

1. Introduction

Jingmen tick virus (JMTV) is an emerging human pathogen first described in Rhipicephalus microplus ticks in Jingmen city (Hubei Province, China) in 2010 [1,2]. The public health importance of the virus has been demonstrated following its detection in patients presenting with mild to severe disease, and in ticks collected from the patients [2]. Jingmen tick virus was implicated in fatal human cases associated with Crimean–Congo hemorrhagic fever virus (CCHFV) in Kosovo in 2019 [3]. Variants of the virus have been detected in non-human primates in Uganda [4]. In Finland, Russia, and China, a JMTV-like virus known as Alongshan virus (ALSV) was recovered in patients with febrile illness and in ticks [5,6,7]. JMTV has also been detected in mosquitoes and in vertebrates including cattle and rodents in China [1,8,9]. However, the ecology of JMTV remains poorly described in Africa.
JMTV and related viruses have so far not been classified by the International Committee on Taxonomy of Viruses (ICTV) and are referred to as jingmenviruses [10,11]. Jingmenviruses form two phylogenetic groups. The viruses of one group are transmitted by blood-feeding arthropods and are vertebrate pathogenic (arbo-jingmenviruses) such as JMTV, Mogiana tick virus (MGTV), and Yanggou tick virus (YGTV) [2,10,11,12]. The second group contains viruses that only infect arthropods (arthropod-specific jingmenviruses) exemplified by the Guaico Culex virus (GCXV), Wuhan aphid virus (WHAV), and Shuangao insect virus (SHIV) groups [8,10,12]. Jingmenviruses have a genome comprising four segments of positive-sense single-stranded RNA and are distantly related to flaviviruses [1,10,11,12]. JMTV has a multisegmented genome comprising three monocistronic segments (segments 1, 2, and 3) and two separate open reading frames (ORF) in one segment (segment 4) [1,10,12]. The non-structural proteins of segment 1 (NSP1) and segment 3 (NSP2) seem to be homologous to those of the genus flavivirus (NS5, RNA-dependent RNA polymerase (RdRP) and NS3, helicase/protease) [1,10,12]. The close phylogenetic relationship of these two proteins provides an unusual evolutionary link between the non-segmented flaviviruses and the segmented jingmenviruses.
To date, there is no information regarding the presence of JMTV infections in vertebrates or ticks in Kenya. Towards predictive surveillance for arboviral threats, this study sought to assess the presence and distribution of JMTV in adult ticks in Kenya. Ticks collected from livestock hosts such as cattle, goats, and sheep and free-ranging tortoises from two pastoralist-predominant areas, Baringo and Kajiado counties were assessed.

2. Material and Methods

2.1. Ethical Approval

All experimental procedures were approved by the Kenya Medical Research Institute Scientific and Ethics Review Unit (KEMRI-SERU) (SERU protocol number 3312). Local authorities and village elders were consulted before the study was initiated and during each sampling visit. Additionally, informed oral consent was obtained from owners and/or caretakers before sampling domestic animals. Tick sampling from livestock was performed in accordance with ethical guidelines.

2.2. Study Sites, Tick Sampling and Morphological Identification

Adult ticks were collected from seven sites in two counties, Baringo and Kajiado, within the Great Rift Valley Region of Kenya. These counties (Figure 1) were selected because they have previous history of arboviral infections [13,14,15,16]. Four sampling sites were included in Baringo County, namely Ntepes, Sandai, Logumgum, and Kapkuikui, and three sites in Kajiado County, namely Oloisinyai, Oldorko, and Soweto (Figure 1). Baringo and Kajiado counties are home to Ruko and Olkirimatian conservancies, respectively, and have a semi-arid ecology. Both counties are inhabited by nomadic pastoralist communities that keep livestock (sheep, cattle, and goats). Both areas have a history of undiagnosed febrile human illness with poorly described aetiologies.
Adult ticks were collected from cattle (Bos taurus), goats (Capra hircus), and sheep (Ovis aries), as well as from tortoises (Leopard tortoise, Stigmochelys pardalis) in Baringo only, between August 2019 and July 2020 after the rains when arboviral activities are expected to be high [18]. Individual livestock specimens (cattle, goats, and sheep) were restrained to allow manual picking of attached adult ticks from the animal’s skin. The ticks were initially taken to a temporary field laboratory and cryopreserved in liquid nitrogen for transportation and subsequent storage in −80 °C freezers in the laboratory at icipe in Nairobi. Following surface sterilization with 70% ethanol and washing with deionized water to remove foreign particles from animals, the collected ticks were identified morphologically to species using established morphologic keys [19,20,21,22]. The ticks were pooled according to species, host, sex, and sampling sites in groups of up to eight individuals and stored in −80 °C freezers until further processing. Ticks that could not be identified to the species level were not pooled and analyzed individually.

2.3. Tortoise Sampling

Tortoise sampling was carried out in September 2019 at Kapkuikui, Baringo county after tick collection from livestock on a single occasion on August 2019. Free-ranging tortoises were physically restrained to allow blood collection from either the dorsal coccygeal vein or brachial sinus by a team including a registered animal health technician in accordance with the animal ethics protocol approved by KEMRI-SERU. About 0.5 mL of whole blood was collected aseptically using 2 mL BD Vacutainer with Ethylenediaminetetraacetic acid (EDTA) (Becton Dickinson, Sunnyvale, CA, USA) or approximately 50 µL whole blood spotted on FTA Classic Card Whatman (Merck KGaA; Darmstadt, Germany). The samples on 2 mL BD Vacutainer were cryopreserved in liquid nitrogen for transportation while the samples on FTA cards were stored at room temperature during transportation. All blood samples were stored at −80 °C during transfer to the laboratory at icipe in Nairobi, until further processing.

2.4. Tick Homogenization and Molecular Identification

Before genomic DNA extraction, ticks were frozen in liquid nitrogen and mechanically homogenized in 1.5 mL microcentrifuge tubes with zirconia beads (2.0 mm and 0.1 mm diameter) for 45 s in Mini-Beadbeater-16 (Biospec, Bartlesville, OK, USA). One mL DPBS (Dulbeccos phosphate-buffered saline) was added to each tick homogenate prior to centrifugation at 2500 revolutions per minute (rpm) at 4 °C for 10 min in a bench top centrifuge (Eppendorf 5430 R). Following the manufacturer’s protocol, genomic DNA was extracted from the pellet of Amblyomma ticks (that could not be identified at the species level by morphology) using the ISOLATE II Genomic DNA extraction kit (Bioline, London, UK). The supernatant was preserved at −80 °C for RNA extraction and virus isolation. Three PCR assays targeting the 16S ribosomal ribonucleic acid (rRNA), internal transcribed spacer 2 (ITS2), and cytochrome oxidase 1 (CO1) genes were employed for amplification (Supplementary Table S1) using the thermal cycling conditions described earlier [23,24,25]. The PCR products were examined on an ethidium bromide-stained 1.5% agarose gel and amplicons of the correct size purified for Sanger sequencing (Macrogen, Amsterdam, The Netherlands) using ExoSAP-IT (USB Corporation, Cleveland, OH, USA).

2.5. RNA Extraction and PCR Screening

Total RNA was extracted from 140 µL of tick homogenates (homogenization described above) and from tortoise whole blood using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The tick homogenates included both samples homogenized in pools and those homogenized individually. The remaining homogenates were preserved at −80 °C for virus isolation. Reverse transcription was carried out with 5 µL of the isolated RNA using the High Capacity cDNA Reverse Transcription (RT) kit (LifeTechnologies, Carlsbad, CA, USA) and 600 μM non-ribosomal random primers [26]. Samples were screened for JMTV infection by pan-JMTV RT-PCR using primers targeting the NS5 gene (Supplementary Table S1). PCRs were set up in a final volume of 25 µL containing 2.5 units of Mytaq DNA polymerase (Bioline, London, UK). The touch down thermal cycling conditions involved initial denaturation for 3 min at 95 °C, followed by 10 cycles of 20 s at 95 °C, 30 s at 64–56 °C and 30 s at 72 °C, and 35 cycles of 20 s at 95 °C, 30 s at 56 °C, and 30 s at 72 °C. The final extension lasted 5 min at 72 °C. The post-PCR analysis and Sanger sequencing were carried out as described above.

2.6. Virus Isolation and Quantification of Viral Genome Copies

Virus isolation attempts from JMTV-positive tick homogenates was performed in Vero E6 (Ceropithecus aethiops) and C6/36 (Aedes albopictus) cells as previously described following filtration and antibiotic treatment approaches [27]. Briefly, before inoculation, cell culture medium was removed from the cells seeded in a 48-well plate (Nunc, Roskilde, Denmark) and 150 µL of cell line-specific medium without additives added. In the filtered approach, 100 µL of each sample (supernatant of tick homogenate) was passed through 0.45 µm sterile membrane filters (Merck Millipore, Billerica, MA, USA) using a syringe and 50 µL inoculated onto each cell line. Following inoculation, cells were allowed to adsorb for an hour and maintained in 300 µL media. Gibco Dulbecco’s modified Eagle’s medium (DMEM), used for maintaining Vero E6 cells, was supplemented with 5% fetal calf serum (FCS; Gibco, Thermo Fisher Scientific, Darmstadt, Germany) and 1% l-glutamine for Vero E6 cells. Gibco Leibovitz’s L-15 medium (L-15), used for maintaining C6/36 cells, was supplemented with 5% FCS. The antibiotic approach involved inoculating 50 µL of unfiltered tick homogenate and the addition of 100 U/mL penicillin together with 100 ug/mL amphotericin B. Vero E6 cells were incubated in a 5% CO2 atmosphere at 37 °C while C6/36 cells were incubated at 28 °C without CO2. The cells were monitored regularly for cytopathic effects (CPE) and four blind passages on fresh cells were carried out after every seven days.
During each passage, 75 µL of cell culture supernatant was collected for viral RNA isolation in the MagNa Pure 96 extraction system (Roche Diagnostics, Rotkreuz, Switzerland) followed by cDNA synthesis using Superscript IV reverse transcriptase and random hexamer primers (Invitrogen, Karlsruhe, Germany) as described by the manufacturer. A quantitative TaqMan Real-time PCR (qRT-PCR) using JMTV-specific primers and probe (Supplementary Table S1) was established to measure the amount of viral genome copies in the collected cell culture supernatants.

2.7. Library Preparation and Next-Generation Sequencing

Viral RNA was extracted from JMTV-positive tick homogenates using QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) as described above without adding carrier RNA. Viral RNA was further extracted from infectious cell culture supernatant using the NucleoSpin® RNA Virus Kit (Macherey-Nagel, Düren, Germany) following the manufacturer’s protocol. The extracted RNA was used to sequence JMTV full genomes on an Illumina MiSeq Next-generation sequencing (NGS) platform as described earlier with modification [28] in that DNA libraries were prepared using the KAPA HyperPlus kit (Roche Diagnostics, Rotkreuz, Switzerland). After de-multiplexing and extracting the raw data in fastq format, the resulting paired reads were quality trimmed and filtered to remove Illumina adapters using BBDuk (filter = 27, trimk = 30; http://jgi.doe.gov/data-and-tools/bb-tools/, accessed on 25 January 2022). The read coverage was then normalized, employing the kmer-based normalization toolBBNorm with a target option of 40 and minimum depth of 6. After trimming and normalization, the duplicate reads were binned using Dedupe. Viral sequence identification involved mapping paired reads against JMTV segments S1, S2, S3, and S4 as well as by de novo assembly using Spades v3.11.1 (http://cab.spbu.ru/software/spades/, accessed on 25 January 2022). Genome analysis was performed using Geneious Prime, and nucleotides (nt) and deduced amino acids (aa) sequences were queried against the GenBank database using blastn and blastx searches [29,30].

2.8. Phylogenetic Analysis and Genome Characterization

The partial sequences (ITS2 and 16S rRNA) for ticks and fragments of NS5 genes for JMTV were analyzed in Geneious Prime and queried in GenBank [30]. Using the L-INS-i algorithm implemented in MAFFT, nucleotide sequences generated in this study were aligned with related sequences and the alignment was edited manually to adjust regions aligned ambiguously [29,31]. Maximum likelihood (ML) phylogenetic trees were estimated using PhyML v. 2.2.4 with the best fit model determined by Modeltest implemented in MEGA-X version 10.2.5 with nodal support being assessed through 1000 bootstrap replications [32].
Full JMTV genomes were analyzed using Geneious prime [30]. In summary, InterProScan implemented in Geneious prime was used to predict the location of transmembrane domains and signal peptides [29,33]. The location of potential N-glycosylation sites in the four segments were determined in the NetNGlyc v1.0 server (http://www.cbs.dtu.dk/services/NetNGlyc/, accessed on 25 January 2022). JMTV phylogenetic trees were inferred from MAFFT alignment using PhyML v. 2.2.4 with General-time-reversible (GTR) substitution models employing 1000 bootstrap replicates. Trees were midpoint rooted.

2.9. Statistical Analysis

The relative abundance of tick species collected from different hosts was estimated binomially and evaluated by Chi square tests at 95% confidence intervals using R version 4.1.2 [34]. To estimate JMTV prevalence in individual tick species, a frequentist model was used to implement maximum likelihood (ML) analysis in an online platform, EpiTools epidemiological calculator (http://epitools.ausvet.com.au/, accessed on 4 November 2021) [24,35,36,37].

2.10. Sequence Accession Numbers

The ITS2 and 16S rRNA gene fragments sequences were deposited in GenBank under the accession numbers ON212401–ON212405 and ON220154–ON220159, respectively. The JMTV NS5 gene sequences were deposited in GenBank under the accession numbers ON158817–ON158867. JMTV coding complete genome sequences were deposited in GenBank under the accession numbers ON186499–ON186526 (Supplementary Table S2).

3. Results

3.1. Tick Collection from Livestock and Tortoises

A total of 4550 ticks were collected from different hosts including cattle, goats, sheep, and tortoises from the study sites (Figure 1). The ticks comprised thirteen species in three genera including Rhipicephalus (Rh.) appendiculatus, Rh. evertsi evertsi, Rh. pulchellus, Hyalomma (Hy.) marginatum, Hy. rufipes, Hy. truncatum, Hy. impeltatum, Hy. albiparmatum, Amblyomma (Am.) gemma, Am. nuttalli, Am. sparsum, Am. Variegatum and Am. lepidum (Supplementary Table S3). Thirty-five Amblyomma ticks collected from tortoises could not be identified at the species level by morphology nor by sequencing of ITS2 and 16S rRNA gene fragments. These ticks are denoted as Am. sp., henceforth. The generated sequences showed 99%–100% pairwise nucleotide identities within the 16S rRNA gene fragment, and 98%–100% pairwise nucleotide identities within the ITS2 fragment. Amplification of the CO1 gene was not successful. Phylogenetic analysis of ITS2 fragment sequences (1008 bp nucleotides in length) revealed that the Am. sp. ticks are most closely related to Am. marmoreum and Am. loculosum, showing pairwise nucleotide identities of 93%–94% to both species (Supplementary Figure S1). Analysis of the 16S rRNA gene fragment (402 bp nucleotides in length) showed a similar relationship of the collected Am. sp. ticks to Am. marmoreum (93%–94% nt similarity). These data suggest that the collected Am. sp. ticks are distinct from formally recognized Amblyomma tick species for which homologous 16S rRNA and ITS2 data are available (Supplementary Figure S1) and most likely represent a so-far unrecognized tick species [38].

3.2. JMTV Infection in Ticks and Tortoises

Pan-JMTV PCR and/or qRT-PCR screening of 718 tick pools (≤ 8 ticks/pool) yielded sixty-seven positive samples. Assuming only one positive tick per pool, this translates to an overall estimated individual-level JMTV prevalence of 1.5% (67/718, 95% CI 1.2–1.9). The highest individual-level JMTV prevalence was observed for Am. sparsum (85.7%, 6/7) and Am. sp. (77.1%, 27/35), both collected from tortoises sampled in Kapkuikui (Baringo county), while Rh. appendiculatus (0.5%,1/30) collected from sheep in Oloisinyai (Kajiado county) had the lowest prevalence (Table 1). Other positive tested tick species included Rh. evertsi evertsi and Hy. truncatum (Table 1). Jingmen tick virus infection was also observed in ticks collected from small ruminants (sheep and goats) while no JMTV was detected in ticks collected from cattle at all sampling sites (Table 1).
Out of the sixty-seven JMTV-positive samples, fifty showed 93%–94% nt identity to JMTV detected in a non-human primate from Uganda (KX377513.1) in the NS5 gene (489 bp; positions 1733–2222) as shown in Figure 2. The other seventeen samples were identified by qRT-PCR and no sequence information from these samples is available (Table 2). Due to the high prevalence of JMTV in ticks collected from tortoises in Kapkuikui, whole blood was sampled from three tortoises from the same sampling locality a month after the initial sampling and screened for JMTV infection by Pan-JMTV PCR and qRT-PCR. Two of the three whole blood samples collected from tortoises tested positive using qRT-PCR (concentration = 1.58–18.2copies/µL). From one sample, a sequence fragment could be amplified using the Pan-JMTV PCR (nucleotide sequence accession number ON158817, Table 2). The sequence showed 97%–100% nt identity to JMTV detected in ticks collected from tortoise (Supplementary Figure S2) and grouped within the tick-associated JMTV-derived sequences (Figure 2).
Phylogenetic analyses based on the NS5 gene sequences showed that the sequences of this study form a monophyletic clade closely related to JMTV identified from non-human primates in Uganda. The diversity of tortoise-associated JMTV was much greater than that of sheep-associated JMTV (Figure 2). JMTV detected in Rh. appendiculatus collected from sheep clustered together in a distinct monophyletic clade embedded in tortoise-tick-associated viruses, suggesting a tick and/or host specific clustering pattern (Figure 2).

3.3. Virus Isolation and Genome Organization

Virus isolation was attempted using representative JMTV-positive samples (n = 42). JMTV genome copies were detected in cell culture supernatants of seven samples following two cell culture passages. No cytopathic effects (CPE) were present. Cell culture supernatants comprised Rh. appendiculatus (KT125) collected from sheep, from Am. sparsum (MT293), and from Am. sp. (samples MT299, MT304, MT305, MT308, and MT314) both collected from tortoise hosts. However, no virus replication was measured after the third passage in both Vero E6 and C6/36 cells. Multiple attempts to further passage the virus after two cell culture passages were not successful.
Coding complete genomes were obtained by NGS from two Am. sp. ticks (MT304 and MT299) inoculated into C6/36 and Vero E6 cells respectively, as well as from two Am. sparsum (MT290 and MT293) and three Am. sp. (MT297, MT308, and MT 328) homogenates. Each genome consisted of four genome segments and had a size of approximately 11 kb excluding the non-coding untranslated regions (UTRs). The genomes showed a typical segmented JMTV genome organization comprising three monocistronic segments (segments 1, 2, and 3) and two separate ORFs in segment 4 (Figure 3A). A 924-amino acid polypeptide on segment 1 was predicted to code for non-structural protein 1 (NSP1). The protein is homologous to the flavivirus NS5 protein and contains an RNA-dependent RNA polymerase (RdRp) motif (Figure 3A). The other motif present in segment 1 included S-adenosyl-L-methionine-dependent methyltransferases (Figure 3A). Further, a transmembrane domain (position 16–35) and two N-glycosylation sites at position 335 and 519 were observed (Figure 3A). Like other JMTVs, motifs A, B, and C on NSP1 were present in all sequences found in this study (Figure 3B). Analysis of RdRp genes did not reveal any unique nonsynonymous nucleotide substitutions compared to other JMTV sequences. Segment 2 was predicted to encode the structural glycoprotein VP1 and contained a potential N-glycosylation site at position 171 and 225 (Figure 3A). Segment 3 was predicted to code for the NSP2 protein, which shares functional homology with the flavivirus NS2b–NS3 complex. Further, it contained P-loop structures with nucleoside triphosphate hydrolases and RNA helicase motifs (Figure 3A). Signal peptides were also identified in NSP2. Similarly, an ATP-binding site (PGAGKTR) and DEAD-box helicase domain on NSP2 were conserved in all sequences from this study (Figure 3C). Analysis of protease/helicase genes did not reveal any unique nonsynonymous nucleotide substitutions to other JMTVs. The two separate ORFs in segment 4 were predicted to encode two polypeptide viral glycoproteins (VP) 2 and 3 (Figure 3A).
Phylogenetic analyses based on all segments showed that the Kenyan sequences formed a monophyletic clade closely related to JMTV identified in a non-human primate from Uganda and in ticks and rodents from China and ticks from Laos (Figure 4).
Notably, all viruses sequenced in this study were closely related to each other at RdRp and Protease/helicase genes with 99%–100% deduced aa identity (Supplementary Figure S3).

4. Discussion

JMTV has recently been found in many countries, but little is known about its presence in Africa. The only report of JMTV in Africa details its detection in a non-human primate from Uganda. Here, we describe the detection of JMTV in several tick species collected from sheep, goats, and tortoises. JMTV genome copies were also detected in whole blood collected from tortoises. This represents the first record of JMTV in Kenya and in a reptilian host.
JMTV was found in tick species that were not associated with JMTV before. Positive tick species included Rh. appendiculatus, Rh. evertsi evertsi, and Hy. truncatum associated with goats and sheep, as well as Am. sparsum, Am. nuttalli, and Am. sp. associated with tortoises. JMTV has been detected in other studies in Rh. microplus in France, China, Brazil, and Trinidad and Tobago [11,39,40]. Other tick species that have been shown to harbor JMTV include Am. variegatum originating from France, Hy. marginatum, and Haemaphysalis inermis originating from Turkey, showing that JMTV infects species of the genera Rhipicephalus, Amblyomma, Haemaphysalis, and Hyalomma [11,41]. The highest JMTV prevalence was detected in Am. sparsum (85.7%, 95% CI 50.6–99.1) and Am. sp. (77.1%, 95% CI 61.6–88.8) collected from tortoises in Kapkuikui. Interestingly, we did not find any tick collected from cattle positive for JMTV, although cattle from China have been shown to be infected with JMTV [1,8]. Instead, ticks infesting sheep and goats were found positive with a percentage of 1.8% (95% CI 1.2–2.5) and 0.3% (95% CI 0.1–0.5), respectively. Together with earlier findings, these results indicate that JMTV is a potential emerging tick-borne pathogen with a wide geographic distribution and infecting a wide diversity of tick species.
Integral to arbovirus disease transmission dynamics is the presence of a pathogen, a competent vector, and a susceptible animal host. In the frame of susceptible animal hosts inhabiting areas in Kenya with history of arboviral disease, tortoises are among widely distributed tick hosts. Tortoises are known to be long-lived reptiles and have low species richness, with only one or two species living in sympatry [42,43]. Common viral infections in tortoise include ranaviruses and herpesviruses [44,45,46]. However, the role of reptiles in arbovirus transmission dynamics is barely investigated despite them living near humans and domestic animals and being parasitized by ticks such as Am. gemma, which is also known to feed on humans and livestock and Am. falsomarmoreum [23,47]. In some cases, arbovirus transmission is facilitated by the presence of a reservoir/maintenance host where the virus circulates latently as in the case of Crimean–Congo hemorrhagic fever virus (CCHFV) in birds and small mammals [48,49]. Screening of tortoises’ blood for JMTV infection revealed the presence of JMTV RNA (nucleotide sequence accession number: ON158817; Table 2) indicative of JMTV presence in both the reptile and ticks from the same geographic location, revising our understanding of JMTV host range. Further, phylogenetic analyses of the NS5 gene revealed that the diversity of tortoise-associated JMTV was much greater than that of sheep-associated JMTV, suggesting that a much wider JMTV diversity may exist in tortoises (Figure 2). More sequence information would be required for a comprehensive phylogenetic analysis to investigate if (ancestral) JMTV variants may have spread to mammals including sheep. Although the feeding pattern of Am. sp. is not known, other tortoise-associated tick species reporting JMTV infection including Am. sparsum and Am. nuttalli have been shown to feed on livestock and humans, respectively [47,50]. Overall, the findings suggest that tortoises could become infected with JMTV. However, as the present study involved only a small number of tortoises, a focused JMTV survey including antibody screening could provide more insights into JMTV maintenance in tortoise and transmission between tortoise and ticks.
JMTV isolation attempts suggested only replication for two cell culture passages in Vero E6 and C6/36 cells. While the detection of JMTV RNA in cell culture supernatant may be an indication of virus replication in these cells, it is also possible that the JMTV RNA detections in cell culture supernatants were carryover from tick homogenates as speculated by Kobayashi et al. [51]. Earlier virus isolation attempts have shown mixed outcomes with inconsistent replication in C6/36 as well as in DH82 (canine macrophage) cells, and persistent growth in BME/CTVM23 (Rh. microplus embryo-derived cell line) [1,2,51]. Other jingmenviruses such as ALSV has been shown to replicate in Vero and IRE/CTVM19 (Ixodes ricinus) cell lines [52,53]. In the present study, JMTV was detected from seven out of forty-two cell culture supernatants only up to the second blind passages, highlighting the need to establish a more efficient JMTV isolation approach.
Seven complete-coding genomes of JMTV were sequenced and characterized. A comprehensive analysis of coding regions showed a genome structure and functional organization comparable to those of other reported JMTVs. The RNA helicase in segment 3 was comparable to flavivirus NS3 helicase responsible for viral RNA capping and synthesis [54]. The segment plays an important role in polyprotein processing and genome replication with NS3 involved in ATP-dependent RNA or dsRNA unwinding activities [1]. Further, similarities in genome segments organization, amino acid motifs, and signal peptides were noted. In addition, the monophyletic grouping of the current samples suggests that like other JMTVs, the segmented genome may be a product of a single evolutionary event [10]. However, a lot remains to be investigated on the evolutionary history of JMTV especially following its detection in tortoises, an animal known to be long living and having low species richness.

5. Conclusions

Our findings showed circulation of JMTV in several tick species in Kenya. Virus-positive ticks were collected from sheep, goats, and tortoises, with the latter showing the highest JMTV infection rates. Interestingly, the diversity of tortoise-associated JMTV was much greater than that of sheep-associated JMTV. Further epidemiological surveys including antibody screening of livestock and humans are necessary to appraise the potential health risk posed by JMTV in the study areas and beyond.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v14051041/s1, Figure S1. Phylogenetic relationship of Amblyomma ticks collected from tortoises (that could not be identified at the species level using morphological keys), and other Amblyomma congeners, Figure S2. JMTV distance matrix, Figure S3. Jingmenviruses distance matrixes, Table S1: Primers, probes and PCR conditions used in the study [38,55,56], Table S2: Sequences generated in this study submitted to GenBank, Table S3: Relative abundance of ticks sampled from different hosts.

Author Contributions

Conceptualization, D.P.T., S.J., A.D.S.B., B.T., R.S. and E.O.O.; methodology, E.O.O., D.P.T. and M.M.; software, E.O.O., S.J. and A.K.; validation, E.O.O., D.P.T., A.K., M.M. and S.J.; formal analysis, E.O.O., D.O., I.S., G.R. and C.G.; investigation, E.O.O., D.O., I.S., G.R. and C.G.; resources, S.J., D.P.T., B.T., R.S. and A.D.S.B.; data curation, E.O.O., S.J., D.P.T. and A.D.S.B.; writing—original draft preparation, review and editing, E.O.O., D.P.T., S.J., A.D.S.B., B.T., R.S., A.K., G.R., D.O., C.G., M.M. and I.S.; supervision, D.P.T., S.J., B.T., R.S. and A.D.S.B.; funding acquisition, S.J., D.P.T., B.T. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

The work was funded by the Deutsche Forschungsgemeinschaft (JU 2857/9-1 and JU 2857/9-2 to S.J.) and the German Center for Infection Research (DZIF), Germany (TTU 01.801). Edwin O. Ogola was supported by a German Academic Exchange Service (DAAD) through the icipe ARPPIS-DAAD scholarship and a UP postgraduate bursary. David P. Tchouassi is supported by a Wellcome Trust International Intermediate Fellowship (222005/Z/20/Z) and from the Norad-funded project Combatting Arthropod Pests for better Health, Food and Climate Resilience (CAP-Africa; project number RAF-3058 KEN-18/0005). We gratefully acknowledge the financial support for this research by the following organizations and agencies: Swedish International Development Cooperation Agency (Sida), Swiss Agency for Development and Cooperation (SDC), Federal Democratic Republic of Ethiopia, and the Government of the Republic of Kenya. The views expressed herein do not necessarily reflect the official opinion of the donors. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Sequences generated were deposited to GenBank under accession numbers ON158817–ON158867, ON186499–ON186526, ON220154–ON220159 and ON212401–ON212405. Other data presented in the study are available in the article and as supplements.

Acknowledgments

We are grateful for the technical support of Tatenda Chiuya, Josephine Osalla, Kevin Kidambasi and Dickens Ondifu of icipe’s ML-EID laboratory, and Christian Hieke and Verena Hyde of Institute of Virology, Charité Universitätsmediz in Berlin.

Conflicts of Interest

The authors declare that there is no conflict of interests.

Disclosure Statement

No potential conflict of interest was reported by the author(s).

References

  1. Qin, X.C.; Shi, M.; Tian, J.H.; Lin, X.D.; Gao, D.Y.; He, J.R.; Wang, J.B.; Li, C.X.; Kang, Y.J.; Yu, B.; et al. A Tick-Borne Segmented RNA Virus Contains Genome Segments Derived from Unsegmented Viral Ancestors. Proc. Natl. Acad. Sci. USA 2014, 111, 6744–6749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Jia, N.; Liu, H.B.; Ni, X.B.; Bell-Sakyi, L.; Zheng, Y.C.; Song, J.L.; Li, J.; Jiang, B.G.; Wang, Q.; Sun, Y.; et al. Emergence of Human Infection with Jingmen Tick Virus in China: A Retrospective Study. eBioMedicine 2019, 43, 317–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Emmerich, P.; Jakupi, X.; von Possel, R.; Berisha, L.; Halili, B.; Günther, S.; Cadar, D.; Ahmeti, S.; Schmidt-Chanasit, J. Viral Metagenomics, Genetic and Evolutionary Characteristics of Crimean-Congo Hemorrhagic Fever Orthonairovirus in Humans, Kosovo. Infect. Genet. Evol. 2018, 65, 6–11. [Google Scholar] [CrossRef] [PubMed]
  4. Ladner, J.T.; Wiley, M.R.; Beitzel, B.; Auguste, A.J.; Dupuis, A.; Lindquist, M.E.; Sibley, S.D.; Kota, K.P.; Fetterer, D.; Eastwood, G.; et al. A Multicomponent Animal Virus Isolated from Mosquitoes. Physiol. Behav. 2017, 176, 139–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Kuivanen, S.; Levanov, L.; Kareinen, L.; Sironen, T.; Jääskeläinen, A.J.; Plyusnin, I.; Zakham, F.; Emmerich, P.; Schmidt-Chanasit, J.; Hepojoki, J.; et al. Detection of Novel Tick-Borne Pathogen, Alongshan Virus, in Ixodes Ricinus Ticks, South-Eastern Finland, 2019. Eurosurveillance 2019, 24, 1900394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kholodilov, I.S.; Litov, A.G.; Klimentov, A.S.; Belova, O.A.; Polienko, A.E.; Nikitin, N.A.; Shchetinin, A.M.; Ivannikova, A.Y.; Bell-Sakyi, L.; Yakovlev, A.S.; et al. Isolation and Characterisation of Alongshan Virus in Russia. Viruses 2020, 12, 362. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, Z.D.; Wang, B.; Wei, F.; Han, S.Z.; Zhang, L.; Yang, Z.T.; Yan, Y.; Lv, X.L.; Li, L.; Wang, S.C.; et al. A New Segmented Virus Associated with Human Febrile Illness in China. N. Engl. J. Med. 2019, 380, 2116–2125. [Google Scholar] [CrossRef]
  8. Guo, J.J.; Lin, X.D.; Chen, Y.M.; Hao, Z.Y.; Wang, Z.X.; Yu, Z.M.; Lu, M.; Li, K.; Qin, X.C.; Wang, W.; et al. Diversity and Circulation of Jingmen Tick Virus in Ticks and Mammals. Virus Evol. 2020, 6, veaa051. [Google Scholar] [CrossRef]
  9. Yu, Z.M.; Chen, J.T.; Qin, J.; Guo, J.J.; Li, K.; Xu, Q.Y.; Wang, W.; Lu, M.; Qin, X.C.; Zhang, Y.Z. Identification and Characterization of Jingmen Tick Virus in Rodents from Xinjiang, China. Infect. Genet. Evol. 2020, 84, 104411. [Google Scholar] [CrossRef]
  10. Shi, M.; Lin, X.-D.; Vasilakis, N.; Tian, J.-H.; Li, C.-X.; Chen, L.-J.; Eastwood, G.; Diao, X.-N.; Chen, M.-H.; Chen, X.; et al. Divergent Viruses Discovered in Arthropods and Vertebrates Revise the Evolutionary History of the Flaviviridae and Related Viruses. J. Virol. 2016, 90, 659–669. [Google Scholar] [CrossRef] [Green Version]
  11. Temmam, S.; Bigot, T.; Chrétien, D.; Gondard, M.; Pérot, P.; Pommelet, V. Insights into the Host Range, Genetic Diversity, and Geographical Distribution of Jingmenviruses. Ecol. Evol. Sci. 2019, 4, e00645-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Garry, C.E.; Garry, R.F. Proteomics Computational Analyses Suggest That the Envelope Glycoproteins of Segmented Jingmen Flavi-like Viruses Are Class II Viral Fusion Proteins (β-Penetrenes) with Mucin-like Domains. Viruses 2020, 12, 260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tigoi, C.; Lwande, O.; Orindi, B.; Irura, Z.; Ongus, J.; Sang, R. Seroepidemiology of Selected Arboviruses in Febrile Patients Visiting Selected Health Facilities in the Lake/River Basin Areas of Lake Baringo, Lake Naivasha, and Tana River, Kenya. Vector-Borne Zoonotic Dis. 2015, 15, 124–132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ajamma, Y.U.; Onchuru, T.O.; Ouso, D.O.; Omondi, D.; Masiga, D.K.; Villinger, J. Vertical Transmission of Naturally Occurring Bunyamwera and Insect-Specific Flavivirus Infections in Mosquitoes from Islands and Mainland Shores of Lakes Victoria and Baringo in Kenya. PLoS Negl. Trop. Dis. 2018, 12, e0006949. [Google Scholar] [CrossRef] [PubMed]
  15. Tchouassi, D.P.; Marklewitz, M.; Chepkorir, E.; Zirkel, F.; Agha, S.B.; Tigoi, C.C.; Koskei, E.; Drosten, C.; Borgemeister, C.; Torto, B.; et al. Sandfly-Associated Phlebovirus with Evidence of Neutralizing Antibodies in Humans, Kenya. Emerg. Infect. Dis. 2019, 25, 681–690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Marklewitz, M.; Tchouassi, D.P.; Hieke, C.; Heyde, V.; Torto, B.; Sang, R.; Junglen, S. Insights into the Evolutionary Origin of Mediterranean Sandfly Fever Viruses. mSphere 2020, 5, e00598-20. [Google Scholar] [CrossRef]
  17. QGIS Development Team. Welcome to the QGIS Project! Available online: https://qgis.org/en/site/ (accessed on 11 June 2019).
  18. Harrison, B.A.; Whitt, P.B.; Roberts, L.F.; Lehman, J.A.; Lindsey, N.P.; Nasci, R.S.; Hansen, G.R. Rapid Assessment of Mosquitoes and Arbovirus Activity after Floods in Southeastern Kansas, 2007. J. Am. Mosq. Control Assoc. 2009, 25, 265–271. [Google Scholar] [CrossRef]
  19. John, G.; Matthysse, J.G.M. The Ixodid Ticks of Uganda Together with Species Pertinent to Uganda Because of Their Present Known Distribution; Entomological Society of America: Annapolis, MA, USA, 1987. [Google Scholar]
  20. Okello-Onen, J.; Hassan, S.M.; Essuman, S. Taxonomy of African Ticks: An Identification Manual, 1st ed.; International Centre of Insect Physiology and Ecology, African Postgraduate Programme in Insect Science: Nairobi, Kenya, 1999. [Google Scholar]
  21. Walker, J.; Keirans, J.E.; Horak, I.G. Rhipicephalud (Acari, Ixodidae)—A Guide to Brown Ticks of the World; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  22. 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 Municipal Abattoir for Their Technical Support. Africa: A Guide to Identification of Tick Species; Bioscience Reports: Edinburgh, UK, 2003. [Google Scholar]
  23. Wang’ang’a Oundo, J.; Villinger, J.; Jeneby, M.; Ong’amo, G.; Otiende, M.Y.; Makhulu, E.E.; Musa, A.A.; Ouso, D.O.; Wambua, L. Pathogens, Endosymbionts, and Blood-Meal Sources of Host-Seeking Ticks in the Fast-Changing Maasai Mara Wildlife Ecosystem. PLoS ONE 2020, 15, e0228366. [Google Scholar] [CrossRef]
  24. Chiuya, T.; Masiga, D.K.; Falzon, L.C.; Bastos, A.D.S.; Fèvre, E.M.; Villinger, J. Tick-Borne Pathogens, Including Crime-an-Congo Haemorrhagic Fever Virus, at Livestock Markets and Slaughterhouses in Western Kenya. Transbound. Emerg. Dis. 2020, 68, 2429–2445. [Google Scholar] [CrossRef]
  25. Getange, D.; Bargul, J.L.; Kanduma, E.; Collins, M.; Bodha, B.; Denge, D.; Chiuya, T.; Githaka, N.; Younan, M.; Fèvre, E.M.; et al. Ticks and Tick-Borne Pathogens Associated with Dromedary Camels (Camelus dromedarius) in Northern Kenya. Microorganisms 2021, 9, 1414. [Google Scholar] [CrossRef]
  26. Endoh, D.; Mizutani, T.; Kirisawa, R.; Maki, Y.; Saito, H.; Kon, Y.; Morikawa, S.; Hayashi, M. Species-Independent Detection of RNA Virus by Representational Difference Analysis Using Non-Ribosomal Hexanucleotides for Reverse Transcription. Nucleic Acids Res. 2005, 33, e65. [Google Scholar] [CrossRef] [PubMed]
  27. Junglen, S.; Kopp, A.; Kurth, A.; Pauli, G.; Ellerbrok, H.; Leendertz, F.H. 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, 4462–4468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Marklewitz, M.; Dutari, L.C.; Paraskevopoulou, S.; Page, R.A.; Loaiza, J.R.; Junglen, S. Diverse Novel Phleboviruses in Sandflies from the Panama Canal Area, Central Panama. J. Gen. Virol. 2019, 100, 938–949. [Google Scholar] [CrossRef] [PubMed]
  29. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An Integrated and Extendable Desktop Software Platform for the Organization and Analysis of Sequence Data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  30. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 2014, 215, 403–410. [Google Scholar] [CrossRef]
  31. Katoh, K.; Misawa, K.; Kuma, K.I.; 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] [Green Version]
  32. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [Green Version]
  33. Jones, P.; Binns, D.; Chang, H.Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-Scale Protein Function Classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [Green Version]
  34. R: The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 4 November 2021).
  35. Christopher, J.; Williams, C.M.M. A Critique of Methods of Sampling and Reporting Pathogens in Populations of Fish. J. Aquat. Anim. Health 2001, 13, 300–309. [Google Scholar]
  36. Cowling, D.W.; Gardner, I.A.; Johnson, W.O. Comparison of Methods for Estimation of Individual-Level Prevalence Based on Pooled Samples. Prev. Vet. Med. 1999, 39, 211–225. [Google Scholar] [CrossRef]
  37. Sergeant, E.S.G. Epitools Epidemiological Calculators. Available online: https://epitools.ausvet.com.au/ (accessed on 4 November 2021).
  38. Lv, J.; Wu, S.; Zhang, Y.; Chen, Y.; Feng, C.; Yuan, X.; Jia, G.; Deng, J.; Wang, C.; Wang, Q.; et al. Assessment of Four DNA Fragments (COI, 16S RDNA, ITS2, 12S RDNA) for Species Identification of the Ixodida (Acari: Ixodida). Parasites Vectors 2014, 7, 93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Sameroff, S.; Tokarz, R.; Charles, R.A.; Jain, K.; Oleynik, A.; Che, X.; Georges, K.; Carrington, C.V.; Lipkin, W.I.; Oura, C. Viral Diversity of Tick Species Parasitizing Cattle and Dogs in Trinidad and Tobago. Sci. Rep. 2019, 9, 10421. [Google Scholar] [CrossRef] [Green Version]
  40. De Souza, W.M.; Fumagalli, M.J.; Torres Carrasco, A.D.O.; Romeiro, M.F.; Modha, S.; Seki, M.C.; Gheller, J.M.; Daffre, S.; Nunes, M.R.T.; Murcia, P.R.; et al. Viral Diversity of Rhipicephalus Microplus Parasitizing Cattle in Southern Brazil. Sci. Rep. 2018, 8, 16315. [Google Scholar] [CrossRef] [PubMed]
  41. Dinçer, E.; Hacioǧlu, S.; Kar, S.; Emanet, N.; Brinkmann, A.; Nitsche, A.; Özkul, A.; Linton, Y.M.; Ergünay, K. Survey and Characterization of Jingmen Tick Virus Variants. Viruses 2019, 11, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Hailey, A.; Coulson, I.M. Habitat Association of Tortoises Geochelone Pardalis and Kinixys Spekii in Sengwa Wildlife Research Area, Zimbabwe. J. Herpetol. 1995, 5, 305–309. [Google Scholar]
  43. Demaya, G.S.; Benansio, J.S.; Petrozzi, F.; Di Vittorio, M.; Dendi, D.; John, E.F.; Luiselli, L. Habitat Use and Spatial Niche Overlap of Sympatric Savannah Tortoises at Multiple Spatial Scales in South Sudan. J. Arid Environ. 2020, 183, 104287. [Google Scholar] [CrossRef]
  44. Marschang, R.E. Viruses Infecting Reptiles. Viruses 2011, 3, 2087–2126. [Google Scholar] [CrossRef] [Green Version]
  45. Allender, M.C.; Abd-Eldaim, M.; Schumacher, J.; Mcruer, D.; Christian, L.S.; Kennedy, M. PCR Prevalence of Ranavirus in Free-Ranging Eastern Box Turtles (Terrapene carolina carolina) at Rehabilitation Centers in Three Southeastern US States. J. Wildl. Dis. 2011, 47, 759–764. [Google Scholar] [CrossRef]
  46. Marschang, R.E.; Kolesnik, E.; Mittenzwei, F. Detection of Testudinid Herpesvirus Type 4 in a Leopard Tortoise (Stigmochelys pardalis). Tierärztl Prax. 2016, 44, 283–286. [Google Scholar] [CrossRef]
  47. Omondi, D.; Masiga, D.K.; Fielding, B.C.; Kariuki, E.; Ajamma, Y.U.; Mwamuye, M.M.; Ouso, D.O.; Villinger, J. Molecular Detection of Tick-Borne Pathogen Diversities in Ticks from Livestock and Reptiles along the Shores and Adjacent Islands of Lake Victoria and Lake Baringo, Kenya. Front. Vet. Sci. 2017, 4, 73. [Google Scholar] [CrossRef] [Green Version]
  48. Ergönül, Ö. Crimean-Congo Haemorrhagic Fever. Lancet Infect. Dis. 2006, 6, 203–214. [Google Scholar] [CrossRef]
  49. Spengler, J.R.; Estrada-Pena, A.; Garrison, A.R.; Schmaljohn, C.; Spiropoulou, C.F.; Bergeron, E.; Bente, D. A Chronological Review of Experimental Infection Studies of the Role of Wild Animals and Livestock in the Maintenance and Transmission of Crimean-Congo Hemorrhagic Fever Virus. Antiviral Res. 2016, 135, 31–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Pietzsch, M.; Quest, R.; Hillyard, P.D.; Medlock, J.M.; Leach, S. Importation of Exotic Ticks into the United Kingdom via the International Trade in Reptiles. Exp. Appl. Acarol. 2006, 38, 59–65. [Google Scholar] [CrossRef] [PubMed]
  51. Kobayashi, D.; Kuwata, R.; Kimura, T.; Shimoda, H.; Fujita, R.; Faizah, A.N.; Kai, I.; Matsumura, R.; Kuroda, Y.; Watanabe, S.; et al. Detection of Jingmenviruses in Japan with Evidence of Vertical Transmission in Ticks. Viruses 2021, 13, 2547. [Google Scholar] [CrossRef]
  52. Kholodilov, I.S.; Belova, O.A.; Morozkin, E.S.; Litov, A.G.; Ivannikova, A.Y.; Makenov, M.T.; Shchetinin, A.M.; Aibulatov, S.V.; Bazarova, G.K.; Bell-Sakyi, L.; et al. Geographical and Tick-Dependent Distribution of Flavi-Like Alongshan and Yanggou Tick Viruses in Russia. Viruses 2021, 13, 458. [Google Scholar] [CrossRef]
  53. Wang, Z.D.; Wang, W.; Wang, N.N.; Qiu, K.; Zhang, X.; Tana, G.; Liu, Q.; Zhu, X.Q. Prevalence of the Emerging Novel Alongshan Virus Infection in Sheep and Cattle in Inner Mongolia, Northeastern China. Parasites Vectors 2019, 12, 450. [Google Scholar] [CrossRef] [Green Version]
  54. Ferron, F.; Rancurel, C.; Longhi, S.; Cambillau, C.; Henrissat, B.; Canard, B. VaZyMolO: A Tool to Define and Classify Modularity in Viral Proteins. J. Gen. Virol. 2005, 86, 743–749. [Google Scholar] [CrossRef]
  55. Brahma, R.K.; Dixit, V.; Sangwan, A.K.; Doley, R. Identification and Characterization of Rhipicephalus (Boophilus) microplus and Haemaphysalis bispinosa Ticks (Acari: Ixodidae) of Northeast India by ITS2 and 16S RDNA Sequences and Morphological Analysis. Exp. Appl. Acarol. 2014, 62, 253–265. [Google Scholar] [CrossRef]
  56. Chitimia, L.; Lin, R.Q.; Cosoroaba, I.; Braila, P.; Song, H.Q.; Zhu, X.Q. Molecular Characterization of Hard and Soft Ticks from Romania by Sequences of the Internal Transcribed Spacers of Ribosomal DNA. Parasitol. Res. 2009, 105, 907–911. [Google Scholar] [CrossRef]
Viruses 14 01041 g001 550
Figure 1. Map showing tick collection sites in Kenya. The red dots indicate sampling points. (A), location of Baringo county; (B), location of Kajiado county. The maps were created in the open source GIS software QGIS 2.12 using GPS co-ordinates and shape files derived from Natural Earth (http://www.naturalearthdata.com/, a free GIS data source, accessed on 20 October 2021) and Africa Open data (https://africaopendata.org/dataset/kenya-counties-shapefile/, license Creative Commons, accessed on 20 October 2021) [17].
Figure 1. Map showing tick collection sites in Kenya. The red dots indicate sampling points. (A), location of Baringo county; (B), location of Kajiado county. The maps were created in the open source GIS software QGIS 2.12 using GPS co-ordinates and shape files derived from Natural Earth (http://www.naturalearthdata.com/, a free GIS data source, accessed on 20 October 2021) and Africa Open data (https://africaopendata.org/dataset/kenya-counties-shapefile/, license Creative Commons, accessed on 20 October 2021) [17].
Viruses 14 01041 g001
Viruses 14 01041 g002 550
Figure 2. Phylogenetic relationship of detected JMTV sequences. The ML phylogenetic analysis is based on a 489-nucleotide fragment of the NS5 gene. Sequences detected in the study and a JMTV reference sequence from Uganda (KX377513.1) were aligned using MAFFT and tree inferred using PhyML v. 2.2.4 with GTR substitution models employing 1000 bootstrap replicates. Trees are midpoint rooted. Only bootstrap values exceeding 50% are shown. Species information of infected ticks and the hosts from which ticks have been collected are indicated.
Figure 2. Phylogenetic relationship of detected JMTV sequences. The ML phylogenetic analysis is based on a 489-nucleotide fragment of the NS5 gene. Sequences detected in the study and a JMTV reference sequence from Uganda (KX377513.1) were aligned using MAFFT and tree inferred using PhyML v. 2.2.4 with GTR substitution models employing 1000 bootstrap replicates. Trees are midpoint rooted. Only bootstrap values exceeding 50% are shown. Species information of infected ticks and the hosts from which ticks have been collected are indicated.
Viruses 14 01041 g002
Viruses 14 01041 g003 550
Figure 3. Genome structure of JMTV. (A) genome organization; grey regions represent non cytoplasmic domains while light grey regions show cytoplasmic domains; (B) conserved JMTV RNA-dependent RNA polymerase; motifs A, B, and C represent highly conserved regions, light blue boxes highlight regions of functional significance; (C) protease/helicase motifs of JMTV and viruses identified in the present study.
Figure 3. Genome structure of JMTV. (A) genome organization; grey regions represent non cytoplasmic domains while light grey regions show cytoplasmic domains; (B) conserved JMTV RNA-dependent RNA polymerase; motifs A, B, and C represent highly conserved regions, light blue boxes highlight regions of functional significance; (C) protease/helicase motifs of JMTV and viruses identified in the present study.
Viruses 14 01041 g003
Viruses 14 01041 g004 550
Figure 4. Phylogentic analyses of JMTV complete coding sequences. The maximum likelihood phylogenetic analyses were based on nucleotide sequences of segements 1, 2, 3, and 4. Host information and geographic origin of the sequences are indicated. Phylogenetic trees were inferred from MAFFT alignment using PhyML v. 2.2.4 with General-time-reversible (GTR) substitution models employing 1000 bootstrap replicates. Trees are midpoint rooted. JMTV sequences identified in this study are highlighed in red. Bootstrap support values of more than 65 are shown.
Figure 4. Phylogentic analyses of JMTV complete coding sequences. The maximum likelihood phylogenetic analyses were based on nucleotide sequences of segements 1, 2, 3, and 4. Host information and geographic origin of the sequences are indicated. Phylogenetic trees were inferred from MAFFT alignment using PhyML v. 2.2.4 with General-time-reversible (GTR) substitution models employing 1000 bootstrap replicates. Trees are midpoint rooted. JMTV sequences identified in this study are highlighed in red. Bootstrap support values of more than 65 are shown.
Viruses 14 01041 g004
Table
Table 1. Estimated individual-level JMTV prevalence in ticks collected from different hosts.
Table 1. Estimated individual-level JMTV prevalence in ticks collected from different hosts.
CountySampling SiteTick SpeciesGoats (n)Sheep (n)Tortoises (n) Cattle (n)
KajiadoOloisinyaiRh. appendiculatus0 (0/4)0.5 (1/30)0 (0/0)0 (0/5)
BaringoNtepesRh. appendiculatus0 (0/84)4.4 (16/55)0 (0/0)0 (0/0)
SandaiRh. evertsi evertsi0 (0/0)42.3 (1/2)0 (0/0)0 (0/2)
Rh. appendiculatus0.5 (1/28)0 (0/7)0 (0/0)0 (0/20)
Hy. truncatum0 (0/2)29.3 (1/2)0 (0/0)0 (0/3)
LogumgumHy. truncatum15.5 (1/3)0 (0/2)0 (0/0)0 (0/0)
Rh. appendiculatus1.5 (4/35)3.3 (7/30)0 (0/0)0 (0/0)
KapkuikuiAm. sp.0 (0/0)0 (0/0)77.1 (27/35)0 (0/0)
Am. sparsum0 (0/0)0 (0/0)85.7 (6/7)0 (0/0)
Am. nuttalli0 (0/0)0 (0/0)40.0 (2/5)0 (0/0)
Total (n) 0.3 (6/156)1.8 (26/128)74.5 (35/47)0 (0/30)
n: proportion of positive pools.
Table
Table 2. Characteristics of JMTV-positive ticks and tortoises in Baringo and Kajiado counties, Kenya.
Table 2. Characteristics of JMTV-positive ticks and tortoises in Baringo and Kajiado counties, Kenya.
CodeCountySampling SiteTime of Sample CollectionSpeciesPool SizeHostSequence Length (nt)GenBank Accession No
MT146BaringoNtepes10 August 2019Rh. appendiculatus8♀Sheep594ON158858
MT147 Ntepes10 August 2019Rh. appendiculatus8♀Sheep595ON158857
MT148 Ntepes10 August 2019Rh. appendiculatus8♀Sheep602ON158854
MT149 Ntepes10 August 2019Rh. appendiculatus8♀Sheep594ON158864
MT150 Ntepes10 August 2019Rh. appendiculatus8♀Sheep591ON158852
MT151 Ntepes10 August 2019Rh. appendiculatus8♀Sheep589ON158862
MT152 Ntepes10 August 2019Rh. appendiculatus8♀Sheep593ON158853
MT153 Ntepes10 August 2019Rh. appendiculatus8♀Sheep586ON158856
MT154 Ntepes10 August 2019Rh. appendiculatus8♀Sheep604ON158860
MT155 Ntepes10 August 2019Rh. appendiculatus8♀Sheep588ON158859
MT156 Ntepes10 August 2019Rh. appendiculatus8♀Sheep591ON158851
MT157 Ntepes10 August 2019Rh. appendiculatus8♀Sheep591ON158855
MT158 Ntepes10 August 2019Rh. appendiculatus8♀Sheep590ON158865
MT159 Ntepes10 August 2019Rh. appendiculatus8♀Sheep595ON158861
MT160 Ntepes10 August 2019Rh. appendiculatus8♀Sheep596ON158863
MT161 Ntepes10 August 2019Rh. appendiculatus8♀Sheep592ON158850
MT290 Kapkuikui10 August 2019Am. sparsum1♂Tortoise577ON158844
MT291 Kapkuikui10 August 2019Am. sparsum1♂Tortoise578ON158847
MT292 Kapkuikui10 August 2019Am. sparsum1♂Tortoise580ON158831
MT293 Kapkuikui10 August 2019Am. sparsum1♂Tortoise561ON158846
MT294 Kapkuikui10 August 2019Am. sparsum1♂Tortoise590ON158837
MT295 Kapkuikui10 August 2019Am. sparsum1♂Tortoise539ON158845
MT297 Kapkuikui10 August 2019Amblyomma sp.1♂Tortoise559ON158835
MT298 Kapkuikui10 August 2019Amblyomma sp.1♂Tortoise538ON158841
MT299 Kapkuikui10 August 2019Amblyomma sp.1♂Tortoise532ON158822
MT300 Kapkuikui10 August 2019Amblyomma sp.1♂Tortoise531ON158842
MT302 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise580ON158834
MT303 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise559ON158840
MT304 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise559ON158839
MT305 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise489ON158843
MT306 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise524ON158836
MT307 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise516ON158849
MT308 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise520ON158838
MT309 Kapkuikui12 August 2019Am. nuttalli1♀Tortoise562ON158866
MT314 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise577ON158829
MT317 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise578ON158833
MT319 Kapkuikui12 August 2019Amblyomma sp.1♂Tortoise562ON158827
MT320 Kapkuikui13 August 2019Am. nuttalli1♂Tortoise603ON158832
MT323 Kapkuikui13 August 2019Amblyomma sp.1♂Tortoise588ON158821
MT324 Kapkuikui13August 2019Amblyomma sp.1♂Tortoise584ON158830
MT327 Kapkuikui14 August 2019Amblyomma sp.1♂Tortoise531ON158848
MT328 Kapkuikui14 August 2019Amblyomma sp.1♂Tortoise520ON158820
MT329 Kapkuikui14 August 2019Amblyomma sp.1♂Tortoise531ON158867
MT330 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise514ON158826
MT331 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise519ON158823
MT332 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise565ON158824
MT333 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise568ON158825
MT334 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise575ON158819
MT335 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise575ON158818
MT336 Kapkuikui15 August 2019Amblyomma sp.1♂Tortoise571ON158828
MT4 Logumgum15 August 2019Rh. appendiculatus8♂Sheep‡‡
MT8 Logumgum15 August 2019Rh. appendiculatus8♂Sheep‡‡
MT19 Logumgum15 August 2019Rh. appendiculatus8♀Sheep‡‡
MT23 Logumgum15 August 2019Rh. appendiculatus8♀Sheep‡‡
MT26 Logumgum15 August 2019Rh. appendiculatus8♀Sheep‡‡
MT29 Logumgum15 August 2019Rh. appendiculatus8♀Sheep‡‡
MT31 Logumgum15 August 2019Rh. appendiculatus8♀Sheep‡‡
MT42 Logumgum15 August 2019Rh. appendiculatus8♀Goat‡‡
MT54 Logumgum15 August 2019Rh. appendiculatus8♂Goat‡‡
MT55 Logumgum15 August 2019Rh. appendiculatus8♂Goat‡‡
MT61 Logumgum15 August 2019Hy. truncatum2♂Goat‡‡
MT62 Logumgum15 August 2019Rh. appendiculatus8♀Goat‡‡
MT101 Sandai15 August 2019Rh. appendiculatus8♂Goat‡‡
MT136 Sandai15 August 2019Hy. truncatum2♂Sheep‡‡
MT144 Sandai15 August 2019Rh. evertsi evertsi2♂Sheep‡‡
MT313 Kapkuikui12 Auguast 2019Amblyomma sp.1♂Tortoise‡‡
T3 * Kapkuikui29 September 2019 1♀ 596ON158817
T2 * Kapkuikui29 September 2019 1♀ ‡‡
KT125KajiadoOloisinyai19 July 2020Rh. appendiculatus8♂Sheep‡‡
♀: female; ♂: male; nt: nucleotides; *: sequence was derived from vertebrate; ‡: positive by qPCR only, no sequence information available.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ogola, E.O.; Kopp, A.; Bastos, A.D.S.; Slothouwer, I.; Marklewitz, M.; Omoga, D.; Rotich, G.; Getugi, C.; Sang, R.; Torto, B.; et al. Jingmen Tick Virus in Ticks from Kenya. Viruses 2022, 14, 1041. https://doi.org/10.3390/v14051041

AMA Style

Ogola EO, Kopp A, Bastos ADS, Slothouwer I, Marklewitz M, Omoga D, Rotich G, Getugi C, Sang R, Torto B, et al. Jingmen Tick Virus in Ticks from Kenya. Viruses. 2022; 14(5):1041. https://doi.org/10.3390/v14051041

Chicago/Turabian Style

Ogola, Edwin O., Anne Kopp, Armanda D. S. Bastos, Inga Slothouwer, Marco Marklewitz, Dorcus Omoga, Gilbert Rotich, Caroline Getugi, Rosemary Sang, Baldwyn Torto, and et al. 2022. "Jingmen Tick Virus in Ticks from Kenya" Viruses 14, no. 5: 1041. https://doi.org/10.3390/v14051041

APA Style

Ogola, E. O., Kopp, A., Bastos, A. D. S., Slothouwer, I., Marklewitz, M., Omoga, D., Rotich, G., Getugi, C., Sang, R., Torto, B., Junglen, S., & Tchouassi, D. P. (2022). Jingmen Tick Virus in Ticks from Kenya. Viruses, 14(5), 1041. https://doi.org/10.3390/v14051041

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