Mitochondrial phylogeography and population structure of the cattle tick Rhipicephalus appendiculatus in the African Great Lakes region

Phylogeography and demographic dynamics of Rhipicephalus microplus from North India

Rhipicephalus microplus, a hematophagous vector prevalent in the tropics and subtropics, is responsible for huge economic losses throughout the globe. However, the taxonomy of the tick species, especially prevalent in north India and south China has been challenged in the recent past. The present study attempted to assess the cryptic status of R. microplus ticks of north India based on two mitochondrial markers; the 16S rRNA and cox1 gene sequences. The phylogenetic tree corresponding to both markers demonstrated the presence of three distinct genetic assemblages/ clades of R. microplus. The present study isolates (n = five and seven for the cox1 and 16S rRNA gene sequences, respectively) from north India along with other isolates from India assorted in the R. microplus clade C sensu. Based on the median joining network analysis corresponding to the 16S rRNA gene sequence, 18 haplotypes were recorded, exhibiting a stellate shape, which was indicative of rapid population expansion. For the cox1 gene, the haplotypes corresponding to clades A, B and C were distantly placed with two exceptions. While performing the population structure analysis, low nucleotide (0.04745 ± 0.00416 and 0.01021 ± 0.00146) and high haplotype diversities (0.913 ± 0.032 and 0.794 ± 0.058) were recorded for the different clades of R. microplus based on the cox1 and 16S rRNA mitochondrial markers, respectively. Eventually, high genetic differentiation and low gene flow were recorded among the different clades. A negative value for the neutrality indices (Tajima's D = -1.44125, Fu's Fs = -4.879, Fu and Li's D = -2.78031 and Fu and Li's F = -2.75229) corresponding to the 16S rRNA gene for the overall dataset evinced an expansion of population size. Based on the detailed studies, it was inferred that the R. microplus tick species circulating in north India belonged to clade C sensu, similar to that of the species prevalent in the other parts of the country as well as in the Indian subcontinent.

Dermacentor reticulatus - a tick on its way from glacial refugia to a panmictic Eurasian population

The ornate dog tick (Dermacentor reticulatus) shows a recently expanding geographic distribution. Knowledge on its intraspecific variability, population structure, rate of genetic diversity and divergence, including its evolution and geographic distribution, is crucial to understand its dispersal capacity. All such information would help to evaluate the potential risk of future spread of associated pathogens of medical and veterinary concern. A set of 865 D. reticulatus ticks was collected from 65 localities across 21 countries, from Portugal in the west to Kazakhstan and southern Russia in the east. Cluster analyses of 16 microsatellite loci were combined with nuclear (ITS2, 18S) and mitochondrial (12S, 16S, COI) sequence data to uncover the ticks' population structures and geographical patterns. Approximate Bayesian computation was applied to model evolutionary relationships among the found clusters. Low variability and a weak phylogenetic signal showing an east-west cline were detected both for mitochondrial and nuclear sequence markers. Microsatellite analyses revealed three genetic clusters, where the eastern and western cluster gradient was supplemented by a third, northern cluster. Alternative scenarios could explain such a tripartite population structure by independent formation of clusters in separate refugia, limited gene flow connected with isolation by distance causing a "bipolar pattern", and the northern cluster deriving from admixture between the eastern and western populations. The best supported demographic scenario of this tick species indicates that the northern cluster derived from admixture between the eastern and western populations 441 (median) to 224 (mode) generations ago, suggesting a possible link with the end of the Little Ice Age in Europe.

Comparative mitogenomics elucidates the population genetic structure of Amblyomma testudinarium in Japan and a closely related Amblyomma species in Myanmar

Ticks are the second most important vector capable of transmitting diseases affecting the health of both humans and animals. Amblyomma testudinarium Koch 1844 (Acari: Ixodidae), is a hard tick species having a wide geographic distribution in Asia. In this study, we analyzed the composition of A. testudinarium whole mitogenomes from various geographical regions in Japan and investigated the population structure, demographic patterns, and phylogeographic relationship with other ixodid species. In addition, we characterized a potentially novel tick species closely related to A. testudinarium from Myanmar. Phylogeographic inference and evolutionary dynamics based on the 15 mitochondrial coding genes supported that A. testudinarium population in Japan is resolved into a star‐like haplogroup and suggested a distinct population structure of A. testudinarium from Amami island in Kyushu region. Correlation analysis using Mantel test statistics showed that no significant correlation was observed between the genetic and geographic distances calculated between the A. testudinarium population from different localities in Japan. Finally, demographic analyses, including mismatch analysis and Tajima’s D test, suggested a possibility of recent population expansion occurred within Japanese haplogroup after a bottleneck event. Although A. testudinarium has been considered widespread and common in East and Southeast Asia, the current study suggested that potentially several cryptic Amblyomma spp. closely related to A. testudinarium are present in Asia.

Molecular survey of cattle ticks in Burundi: First report on the presence of the invasive Rhipicephalus microplus tick

A recent research study on prevalence of tick-borne pathogens in Burundi reported high prevalence and endemicity of Theileria parva, Anaplasma marginale and Babesia bigemina infections in cattle. Detailed information about tick species infesting animals, their distribution and genetic diversity in Burundi is outdated and limited. This study therefore assessed the prevalence and genetic diversity of tick species infesting cattle across agroecological zones (AEZs) in Burundi. A cross-sectional study on the occurrence of tick species was conducted in 24 districts of Burundi between October and December 2017. Differential identification and characterization of ticks collected was conducted using tick morphological keys and molecular tools (cox1 and 12S rRNA gene). Chi-square test was used to test for association between agroecological zones and the prevalence of tick species. Phylogenetic relationships were inferred using bayesian and maximum likelihood algorithms. A total of 483 ticks were collected from the five AEZs sampled. Six tick species comprising of Rhipicephalus appendiculatus, R. sanguineus, R. evertsi evertsi, R. microplus, R. decoloratus and Amblyomma variegatum were observed. Rhipicephalus appendiculatus were the most prevalent ticks (~45%). A total of 138 specimens (28%) were found to be Rhipicephalus microplus, suggesting an emerging threat for cattle farmers. Twelve R. appendiculatus cox1 haplotypes were obtained from 106 specimens that were sequenced. Two cox1 haplotypes of R. microplus which clustered into previously reported Clade A were observed. Rhipicephalus sanguineus and R. evertsi evertsi ticks, the vectors of numerous zoonotic pathogens, were collected from cattle, which constitute a high risk for public health. These findings reveal an overlapping distribution of tick vectors in Burundi. The design of ticks and tick-borne diseases control strategies should consider the distribution of different vectors across the AEZs particularly the presence of the highly invasive R. microplus tick in Burundi and the potential risk of introducing the pathogenic Babesia bovis.

A review of recent research on Theileria parva: Implications for the infection and treatment vaccination method for control of East Coast fever

The infection and treatment (ITM) live vaccination method for control of Theileria parva infection in cattle is increasingly being adopted, particularly in Maasai pastoralist systems. Several studies indicate positive impacts on human livelihoods. Importantly, the first detailed protocol for live vaccine production at scale has recently been published. However, quality control and delivery issues constrain vaccination sustainability and deployment. There is evidence that the distribution of T. parva is spreading from endemic areas in East Africa, North into Southern Sudan and West into Cameroon, probably as a result of anthropogenic movement of cattle. It has also recently been demonstrated that in Kenya, T. parva derived from cape buffalo can 'breakthrough' the immunity induced by ITM. However, in Tanzania, breakthrough has not been reported in areas where cattle co-graze with buffalo. It has been confirmed that buffalo in northern Uganda national parks are not infected with T. parva and R. appendiculatus appears to be absent, raising issues regarding vector distribution. Recently, there have been multiple field population genetic studies using variable number tandem repeat (VNTR) sequences and sequencing of antigen genes encoding targets of CD8+ T-cell responses. The VNTR markers generally reveal high levels of diversity. The antigen gene sequences present within the trivalent Muguga cocktail are relatively conserved among cattle transmissible T. parva populations. By contrast, greater genetic diversity is present in antigen genes from T. parva of buffalo origin. There is also evidence from several studies for transmission of components of stocks present within the Muguga cocktail, into field ticks and cattle following induction of a carrier state by immunization. In the short term, this may increase live vaccine effectiveness, through a more homogeneous challenge, but the long-term consequences are unknown.

The Rhipicephalus appendiculatus tick vector of Theileria parva is absent from cape buffalo (Syncerus caffer) populations and associated ecosystems in northern Uganda

Rhipicephalus appendiculatus is the major tick vector of Theileria parva, an apicomplexan protozoan parasite that causes the most economically important and lethal disease of cattle in East and central Africa. The African cape buffalo (Syncerus caffer) is the major wildlife host of T. parva from southern Uganda and Kenya to southern Africa. We show herein that R. appendiculatus appears to be absent from the two largest national parks in northern Uganda. Syncerus caffer is common in both of these national parks, specifically Murchison falls (MFNP) and Kidepo Valley (KVNP). We re-confirmed the previously reported absence of T. parva in buffalo sampled in the two northern parks based on RLB data using a nested PCR based on the T. parva p104 gene. By contrast, T. parva-infected R. appendiculatus ticks and parasite-infected buffalo were present in Lake Mburo (LMNP) in South central Uganda. This suggests that the distribution of R. appendiculatus, which is predicted to include the higher rainfall regions of northern Uganda, may be limited by additional, as yet unknown factors.

Genetic and antigenic variation of the bovine tick-borne pathogen Theileria parva in the Great Lakes region of Central Africa

Background

Theileria parva causes East Coast fever (ECF), one of the most economically important tick-borne diseases of cattle in sub-Saharan Africa. A live immunisation approach using the infection and treatment method (ITM) provides a strong long-term strain-restricted immunity. However, it typically induces a tick-transmissible carrier state in cattle and may lead to spread of antigenically distinct parasites. Thus, understanding the genetic composition of T. parva is needed prior to the use of the ITM vaccine in new areas. This study examined the sequence diversity and the evolutionary and biogeographical dynamics of T. parva within the African Great Lakes region to better understand the epidemiology of ECF and to assure vaccine safety. Genetic analyses were performed using sequences of two antigen-coding genes, Tp1 and Tp2, generated among 119 T. parva samples collected from cattle in four agro-ecological zones of DRC and Burundi.

Results

The results provided evidence of nucleotide and amino acid polymorphisms in both antigens, resulting in 11 and 10 distinct nucleotide alleles, that predicted 6 and 9 protein variants in Tp1 and Tp2, respectively. Theileria parva samples showed high variation within populations and a moderate biogeographical sub-structuring due to the widespread major genotypes. The diversity was greater in samples from lowlands and midlands areas compared to those from highlands and other African countries. The evolutionary dynamics modelling revealed a signal of selective evolution which was not preferentially detected within the epitope-coding regions, suggesting that the observed polymorphism could be more related to gene flow rather than recent host immune-based selection. Most alleles isolated in the Great Lakes region were closely related to the components of the trivalent Muguga vaccine.

Conclusions

Our findings suggest that the extensive sequence diversity of T. parva and its biogeographical distribution mainly depend on host migration and agro-ecological conditions driving tick population dynamics. Such patterns are likely to contribute to the epidemic and unstable endemic situations of ECF in the region. However, the fact that ubiquitous alleles are genetically similar to the components of the Muguga vaccine together with the limited geographical clustering may justify testing the existing trivalent vaccine for cross-immunity in the region.

Mitochondrial and nuclear multilocus phylogeny of Rhipicephalus ticks from Kenya

The morphological diversity of African ticks of the genus Rhipicephalus and subgenus Boophilus have been studied in detail. However, their taxonomy remains poorly resolved with limited molecular studies performed to improve inter-species discrimination. Herein, ribosomal cytochrome c oxidase I (COI), 12S ribosomal DNA (12S rDNA) and nuclear ribosomal DNA internal transcriber spacer 2 (ITS2) were analyzed in Rhipicephalus tick populations in Kenya. While the morphological and molecular criteria separated R. e. evertsi, R. pulchellus and R. appendiculatus from other members of the genus, except the morphologically similar sibling species R. zambeziensis, this was not the case for other tick populations. COI sequences of Rhipicephalus ticks from Ruma National Park (RNP) in Southwestern Kenya, that were morphologically similar to R. praetextatus/R. simus, a formed distinct clade and barcode gap group. 12S rDNA haplotypes of this population were 99% identical to a GenBank accession of R. muhsamae which is thought to be endemic in West and Central Africa. However, the ITS2 locus indicated that the RNP samples were genetically closest to ticks identified morphologically as R. praetextatus. The COI and 12S rDNA haplotype sequences of R. praetextatus clustered closely with R. simus reference sequences though the two species occurred in distinct barcode gap groups. Our results suggest that the R. simus/R. praetextatus/R. muhsamae comprise a closely related tick species complex found across sub-Saharan Africa and includes the yet to be described RNP population. More studies on the biology, ecology and genomics of all life stages of tick species in the complex may clarify their taxonomic status. A continent-wide study that combines morphology, DNA marker sequencing and emerging methods, such as mass spectrometry and whole-genome resequencing may reveal the diversity and distribution of taxa within the genus Rhipicephalus in sub-Saharan Africa.

Mitochondrial Gene Heterogeneity and Population Genetics of Haemaphysalis longicornis (Acari: Ixodidae) in China

INTRODUCTION

Haemaphysalis longicornis is an important ectoparasite of domestic and wild animals that can transmit many pathogens including viruses, fungi, bacteria and protozoa.

MATERIALS AND METHODS

In this study, we examined genetic variation and population genetics in three mitochondrial (mt) genes [cox1 (cytochrome c subunit 1), rrnL (large subunit ribosomal RNA) and nad5 (NADH dehydrogenase 5)] among four H. longicornis populations from China.

RESULTS

The sizes of the partial sequences of cox1, rrnL and nad5 were 776 bp, 409 bp, 510 bp, respectively. Among the obtained sequences, we identified 22 haplotypes for cox1, 2 haplotypes for rrnL and 17 haplotypes for nad5. Low gene flow and significant genetic differentiation (66.2%) were detected among H. longicornis populations. There was no rapid expansion event in the demographic history of four H. longicornis populations in China. In addition, phylogenetic analyses confirmed that all the Haemaphysalis isolates were H. longicornis which were segregated into two major clades.

CONCLUSION

The mt DNA genes provide a potential novel genetic marker for molecular epidemiology of H. longicornis and assist in the control of tick and tick-borne diseases in humans and animals.