Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
QTL Mapping of Tiller Number in Korean Japonica Rice Varieties
Previous Article in Journal
Mitogenomic Characterization of Cameroonian Endemic Coptodon camerunensis (Cichliformes: Cichlidae) and Matrilineal Phylogeny of Old-World Cichlids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morphological and Molecular Identification of Hard Ticks in Hainan Island, China

1
Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 570100, China
2
Laboratory of Tropical Veterinary Medicine and Vector Biology, School of Life Sciences, Hainan University, Haikou 570228, China
3
International School of Public Health and One Health, Hainan Medical University, Haikou 571199, China
4
Provincial Engineering Technology Research Center for Biological Vector Control, Guangzhou 510080, China
5
Key Laboratory of Tropical Disease Control, Ministry of Education, Sun Yat-Sen University, Guangzhou 510080, China
6
Hainan Provincial Bureau of Disease Prevention and Control, Haikou 570100, China
*
Authors to whom correspondence should be addressed.
Genes 2023, 14(8), 1592; https://doi.org/10.3390/genes14081592
Submission received: 6 July 2023 / Revised: 29 July 2023 / Accepted: 4 August 2023 / Published: 6 August 2023

Abstract

:
Ticks are small, blood-sucking arachnids, known vectors of various diseases, and found throughout the world. They are distributed basically in almost all regions of China. At present, there is not much information regarding tick species on Hainan Island. They were subjected to morphological identification and imaging on an individual basis. Molecular phylogenetic analyses, based on cox1 and 16S rRNA genes, were utilized to identify the species and determine their approximate phylogenetic origin and genetic diversity. The genomic DNA of tick species was extracted, and cytochrome oxidase subunit 1 (cox1) and 16S ribosomal RNA (rRNA) genes were amplified and sequenced. The identification of five tick species, namely Rhipicephalus microplus, Rhipicephalus sanguineus, Rhipicephalus haemaphysaloides, Haemaphysalis cornigera and Haemaphysalis mageshimaensis, was carried out by morphological analysis. When employing the cox1 and 16S rRNA phylogenetic tree, all isolates of R. microplus from Hainan Island were classified as clade A and B, respectively. R. sanguineus was recognized as a member of the tropical lineage by phylogenetic analysis on the cox1 and 16S rRNA genes. Three phylogenetic groups of R. haemaphysaloides were recognized and found to be related closely to strains from China. H. cornigera and H. mageshimaensis formed one phylogenetic group, presumably from tick strains prevalent in Japan and China. The haplotype network analysis indicated that R. microplus is classed into 26 and 6 haplotypes, which correspond to cox1 and 16S rRNA gene assemblages, respectively. In addition, four cox1 haplotypes were detected in R. sanguineus. This is the first evidence that suggests genetic diversity, host range and geographical distribution of hard ticks in Hainan Island, China.

1. Introduction

Ticks are blood-sucking ectoparasitic arthropods that play an essential role in transmitting various pathogens to humans and animals worldwide [1,2], they are the second most significant transmitters of infectious diseases after mosquitoes [3,4]. Most of the pathogenic organisms that inhabit ticks include viruses, bacteria (especially rickettsiae and spirochetes), protozoa and helminths, with their diversity increasing over the past three decades [5,6]. Tick-borne illnesses have been reported in nearly all Provinces/Autonomous Regions/Municipalities (P/A/M) in China in recent years, and the incidence rate is on the rise [7,8]. Roughly 800 tick species and 18 genera have been recognized globally, whereas, 9 genera have been identified in China, with 111 and 14 species of hard (Ixodidae) and soft ticks (Argasidae), respectively [9,10]. Morphological identification refers to the process of visually identifying tick species based on their morphological characteristics, such as body shape, size, color and pattern [11,12]. This involves the use of specialized keys or taxonomic guides that provide detailed descriptions and illustrations of different tick species, as well as their geographic distribution, host preferences and pathogen transmission potential [13]. However, morphological identification is insufficient in separating between related species complexes, particularly when the specimens are at immature stages, physically damaged or engorged [14]. Therefore, molecular techniques such as sequencing of the cytochrome oxidase subunit 1 (cox1) gene and 16S ribosomal RNA (16S rRNA) gene can be used to characterize tick species [15,16]. The cox1 gene is used commonly for species identification and phylogenetic studies [17], while the 16S rRNA gene is a useful marker for determining the evolutionary relationships between different tick species [18]. These techniques can be used in combination with morphological identification to provide a more accurate and comprehensive tick species characterization [19,20].
Hainan Island is located in the southernmost part of China and it has a hot climate. Climatic conditions have shaped rich animal and plant resources on the island and provided a suitable environment for the reproduction of ticks. Previous research on tick species distribution on domestic and small wild animals in Hainan has been investigated [21,22]. However, these studies were conducted over a long period of time, in which climate change and changes in the range of ticks occurred, and no update on data for ticks has been carried out in Hainan Island in the past ten years. At present, there is little information regarding tick species on Hainan Island. Therefore, this study conducted a survey on the distribution of ticks in Hainan Island from July to December 2022 in order to report the morphological features and genetic diversity of field-collected ticks within hard ticks collected from dogs, goats and cattle in the area. Differentiations of hard ticks were characterized using morphological features and nucleotide analysis of the cox1 mitochondrial gene and 16S rRNA gene.

2. Materials and Methods

2.1. Study Sites and Tick Collection

Tick samples were collected in Hainan Island, China, from July to December 2022. Hainan Island is located in southern China, and it has a hot climate, sufficient sunshine and no snow throughout the year. Its climate is mostly tropical with between 1500 and 2000 millimeters of average annual precipitation and average temperature ranging from 23 to 26 °C, and it sits on a geographical coordinate of 20°1′ N latitude and 110°20′ E longitude. On the island, a total of 858 adult and nymph ticks were collected from cattle, dogs and goats in 24 sites and 12 districts of cowshed farms, dog houses and farms and goat farms, respectively (Figure 1 and Supplementary Table S1). Sampling points, geographical coordinates, tick samples, hosts, habitat and time were recorded. Ticks were collected directly from animal hosts by using steel forceps, before placing them in sterile tubes for transportation to the laboratory at low temperature. Then, they were rinsed 3 times with 70% ethanol, rinsed once with phosphate-buffered saline (PBS), air dried and kept at −80 °C until further processing.

2.2. Morphological Identification and Photography

Based on morphological characteristics, the ticks were recognized under a light stereomicroscope Leica S9E (Leica Microsystems, Heerbrugg, Switzerland), and their species were identified using a previously reported taxonomic key [13,23]. The genus, species, stage and gender of the ticks were sorted and counted. Out of the 858 ticks examined, 199 (2–30 specimens per site) were sampled and photographed, with a total of 15,000 pictures taken with a stereomicroscope Leica MZ10F connected to a computer, and the processing and measurement of the images were carried out using Leica Application Suite X version 5.0 software (Leica Microsystems, Heerbrugg, Switzerland). All of the selected ticks were observed, and the following characteristics were taken: dorsal whole-body view; ventral whole-body view; coxae; scutum; dorsal capitulum view; ventral capitulum view; genital opening; anal aperture; spiracular plates. Approximately 5–15 captures were taken in each position.

2.3. DNA Extraction and PCR Amplification of the Tick cox1 and 16S rRNA Genes

After morphological study, 199 selected specimens from different areas, hosts and tick species were subjected to molecular study. The QIAamp Blood and Tissue Kit (Qiagen, Hilden, Germany) was utilized in accordance with the manufacturer’s protocol to perform DNA extraction. Individual ticks were rinsed with PBS and air dried for 5 min on sterile paper. Only semi and fully engorged ticks were sliced separately into small pieces along the length of the body by a sterile scalpel blade. They were then placed in a 2 mL microcentrifuge tube filled with glass beads (900 mg, 0.1 and 3 mm in diameter), together with 180 microliters of ATL buffer and 20 µL of proteinase K, before homogenizing while shaking for 3 minutes in a tissue grinder at 60 HZ for extracting DNA. The extracted DNA was frozen at −20 °C until used.
Individual specimens for each tick species were subjected to polymerase chain reaction (PCR) to amplify regions of the cox1 and 16S rRNA genes. PCR conditions and primers have been summarized in Supplementary Table S2. The PCR reaction mix (final volume 25 µL) included 12.5 µL of 2X TaKaRa Ex Premier™ DNA Polymerase Dye plus (Takara, Japan), 1.5 µL of 5 µM forward and reverse primers, 1–4 µL of DNA template (approximately 100 ng of genomic DNA) and 5.5–8.5 µL of ddH2O. A negative no-template water control was included in all of the PCR runs. Each amplified product was loaded onto 1.5% agarose gel (Biowest, Shanghai, China) using five microliter volumes per sample and the resulting bands were visualized on a gel documentation system (JUNYI, Beijing, China). After purification using the QIAquick® PCR Purification Kit (Qiagen, Hilden, Germany), by following the manufacturer’s instructions, the PCR amplified products were sent subsequently to the Sangon Company (Guangzhou, China) for Sanger bidirectional sequencing.

2.4. Phylogenetic Analyses

Sequences of cox1 and 16S rRNA genes were edited manually and assembled, and bidirectional consensus sequences were originated using BioEdit software version 7.2.6.1 [24]. The sequences were compared with those available in GenBank using the Basic Local Alignment Search Tool (BLAST) available at the National Center for Biotechnology Information (NCBI) website “https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 4 March 2023)”. The multiple sequences were performed using the default parameters of the ClustalW multiple alignment tool in MEGA 11. The phylogenetic analysis was performed using the maximum likelihood (ML) or neighbor-joining (NJ) methods, based on the general time reversible (GTR) model [25] or Kimura’s two-parameter model (K2P) [26] in MEGA 11 [27]. Bootstrap values were estimated for 1000 replicates. The phylogenetic trees were edited using adobe illustration (AI) software.

2.5. DNA Polymorphism Analysis

DNA sequences from each tick species and gene fragment (cox1 and 16S rRNA) were defined as sequence sets that estimated the genetic differentiation and gene flow by using DnaSP software version 6.12.03 [28]. The number of sequences (N), number of haplotypes (Hn), number of segregating sites (n), haplotype diversity (Hd), average number of nucleotide differences (k), nucleotide diversity (π) and FST values between pairs of populations (pairwise FST) were generated. The haplotype networks were created to evaluate the ancestral relationship between detected haplotypes using the Median Joining Network [29] in PopART software version 1.7 [30].

3. Results

3.1. Morphological Features of Tick Species

A total of 858 ticks were collected from 94 domestic animals, including 74 cattle, 11 dogs and 9 goats (Supplementary Table S1). Most of the ticks collected from hosts were in the adult stage (n = 720), from which 550 were female (76.39%) and 170 male (23.61%). There were 138 female nymphs (16.08%) from the total number of samples. The percentage of unfed, semi-engorged and fully engorged ticks was found to be 25.41% (n = 218), 22.84% (n = 196) and 51.75% (n = 444), respectively. The tick species with the highest prevalence was Rhipicephalus microplus 70.62% (n = 606), followed by R. sanguineus sensu lato (s.l.) tropical lineage (R. linnaei) 28.79% (n = 247), R. haemaphysaloides 0.35% (n = 3), H. cornigera 0.12% (n = 1) and H. mageshimaensis 0.12% (n = 1). R. sanguineus (R. linnaei) ticks were collected from dogs, R. haemaphysaloides, H. cornigera and H. mageshimaensis from cattle and R. microplus from cattle and goats. Morphologically, five species of ticks belonged to two genera, namely Rhipicephalus and Haemaphysalis. The capitulum (mouthparts) is the anterior portion of the body used to distinguish between these genera. It is angulated on the sides with triangular porose areas in the genus Rhipicephalus. The palpi are short and broad with no transverse ridges (Figure 2A–F). The capitulum of the genus Haemaphysalis is not angulated on the sides, but has large porose areas that are longitudinal and distant. The palpi is short and its second segment has an acute outward basal prolongation (Figure 2G–J). The distinguishing morphological features of R. microplus, R. sanguineus (R. linnaei), R. haemaphysaloides, H. cornigera and H. mageshimaensis are presented in Supplementary Table S3. The tick species in this study could be distinguished based on color, body, capitulum, coxae, genital opening, scutum (female), anal plates (male), accessory adanal plates (male) and spiracular plate shapes (Supplementary Table S3 and Figures S1–S3).

3.2. Molecular Identification and Classification of Ticks by Nucleotide BLAST

A total of 199 specimens: R. microplus (n = 152), R. sanguineus (R. linnaei) (n = 42), R. haemaphysaloides (n = 3), H. cornigera (n = 1) and H. mageshimaensis (n = 1) were used. All nucleotide sequences obtained in this study were submitted to GenBank with Accession Numbers: OQ704485‒OQ704683 and OQ725381‒OQ725579 for cox1 and 16S rRNA, respectively. Partial cox1 and the 16S rRNA DNA sequence of all ticks were 96–100% of sequence matching identities in GenBank using BLAST (Supplementary Table S4). The closest sequence to those of R. microplus, based on cox1 genes of 99–100% identity, was from Columbia (KT906178.1), Haikou, China (MK685985.1), Kenya (KX228549.1) and Benin (MT249801.1), and to those of R. sanguineus (R. linnaei) was from Angola (MF425995.1). The closest sequence to those of R. haemaphysaloides with 96% identity, H. cornigera with 99% identity and H. mageshimaensis with 99% identity, based on cox1 genes, was from Yingtan, China (OP050242.1), Ganzhou, China (OM368283.1) and Haikou, China (NC062163.1), respectively. The closest sequence to those of R. microplus, based on 16S rRNA genes of 99–100% identity, was from Thailand (KC170742.1) and Mozambique (EU918187.1). The closest sequence to those of R. sanguineus (R. linnaei) with 100% identity, R. haemaphysaloides with 98% identity, H. cornigera with 99% identity and H. mageshimaensis with 98% identity, based on 16S rRNA genes, was from Thailand (KC170744.1), Taiwan (AY972533.1), Ganzhou, China (OM368283.1) and Haikou, China (NC062163.1), respectively.

3.2.1. Genetic Distances and Phylogenetic Analyses for R. microplus

Genetic analysis of R. microplus complex comprises 5 taxa, including R. annulatus, R. australis and R. microplus Clade A, B and C, based on the cox1 gene [15]. Clade A comprised R. microplus from Colombia, Kenya, South Africa, Brazil, Thailand, etc., Clade B contained R. microplus from China, and Clade C consisted of R. microplus from India, Pakistan, Bangladesh and Myanmar. Clade A and C exhibited a sibling association with R. australis and R. annulatus, respectively, indicating a close relationship. In this study, all of the R. microplus ticks and five isolates from Colombia, China, Kenya, Brazil and Thailand were clustered within Clade A (Figure 3). There was variation in the mean K2P distances between and within groups from 0.050 to 0.104 and 0.000 to 0.016, respectively (Supplementary Table S5).
The 16S rRNA genes revealed 4 taxa of R. microplus complex comprising R. australis, R. annulatus and R. microplus Clade A and B. The R. microplus isolates from India, Pakistan and China formed Clade A, while those from this study were Clade B, which clustered together with Thailand, China, Taiwan, Malaysia, Africa and South America (Figure 4). The mean K2P distances between and within groups varied from 0.012 to 0.048 and 0.000 to 0.012, respectively (Supplementary Table S6).

3.2.2. Genetic Distances and Phylogenetic Analyses for R. sanguineus

The cox1 and 16S rRNA phylogenetic trees of R. sanguineus revealed three genetic lineages, including tropical, temperate and southeast European lineages. Both genes of all the R. sanguineus ticks in this study were clustered into one clade of tropical lineage, and closely related to those from Thailand, China, India, Cuba, Brazil, Iraq, Angola, Pakistan and Sub-Saharan Africa (Figure 5). R. sanguineus ticks in the southeast European lineage showed a sister relationship with the tropical lineage and were closer genetically to R. turanicus, while those in the temperate lineage formed a distinct clade away from the tropical lineage. There was variation in the mean K2P distances between and within groups of cox1 and 16S rRNA, ranging from 0.090 to 0.121 and 0.030 to 0.071 and 0.000 to 0.055 and 0.000 to 0.027, respectively (Supplementary Tables S7 and S8).

3.2.3. Genetic Distances and Phylogenetic Analyses for R. haemaphysaloides

The phylogenetic tree of R. haemaphysaloides, based on cox1 and 16S rRNA genes, consisted of three groups. Group 1 comprised R. haemaphysaloides ticks from China (Hunan, Sichuan and Yunnan) and Thailand, which were related closely to R. sanguineus and R. turanicus. Group 2 consisted of ticks from Pakistan, India and Sri Lanka, and Group 3 included ticks from this study and China (Taiwan, Yangxin, Yingtan, Ganzhou and China–Myanmar border) (Figure 6). There was variation observed in the mean K2P distances between and within groups of cox1 and 16S rRNA, with a range of 0.086 to 0.150 and 0.055 to 0.073, and 0.000 to 0.038 and 0.000 to 0.049, respectively (Supplementary Tables S9 and S10).

3.2.4. Genetic Distances and Phylogenetic Analyses for Haemaphysalis spp.

H. cornigera and H. mageshimaensis ticks, and four sequences of cox1 and 16S rRNA genes from this study formed a single clade and were analyzed for phylogenetic trees (Figure 7). According to phylogenetic analysis, based on cox1 genes, H. cornigera and H. mageshimaensis ticks in this study were clustered together in Ganzhou, China and Haikou, China, respectively (Figure 7A). The mean K2P distances between and within species varied from 0.000 to 0.003 and 0.000 to 0.003, respectively (Supplementary Tables S11 and S12). Similarly, 16S rRNA genes of the H. cornigera ticks in this study were clustered alongside the isolates from Japan and Ganzhou, China, while H. mageshimaensis ticks in this study were clustered together with those collected from Japan and China (Taiwan and Haikou) (Figure 7B). The mean K2P distances for the cox1 gene was 0.005, and they varied from 0.000 to 0.003 for the 16S rRNA gene (Supplementary Tables S11 and S12).

3.3. Genetic Diversity and Haplotype Analyses

A total of 398 cox1 and 16S rRNA gene sequences from five tick species were analyzed (Supplementary Table S4). Twenty-six haplotypes of the cox1 sequences and 6 of the 16S rRNA sequences were identified from the 17 populations of R. microplus ticks in Hainan Island, with a haplotype diversity (Hd) of 0.840 and 0.109, respectively. Fifty-two sequences of each gene from the six populations of R. sanguineus (R. linnaei) ticks were analyzed, and four haplotypes of the cox1 sequences were detected with a low haplotype diversity (Hd = 0.216). In contrast, the 16S rRNA sequences showed no variation site and they generated only 1 haplotype, while two, one and one haplotype were identified from 1 population of R. haemaphysaloides, H. cornigera and H. mageshimaensis ticks, respectively (Supplementary Table S13). From this study, only R. microplus ticks, based on cox1 genes, generated a higher haplotype diversity with the Hd ranging from 0.200 to 1.000 (mean = 0.840) and nucleotide diversity (π) ranging from 0.00171 to 0.00389 (mean = 0.00334). The pairwise FST values among populations ranged from –0.5128 to 0.822 (Supplementary Table S14). In 136 comparisons, 86 showed a significantly higher genetic differentiation (p = 0.05).
The median-joining haplotype network, based on cox1 and/or 16S rRNA sequences of R. microplus and R. sanguineus (R. linnaei) in each population is shown in Supplementary Figure S4A–C and Supplementary Table S15. In the frequency haplotype analysis of R. microplus ticks, based on cox1 genes, H1 was the most frequent haplotype, represented by 40 from 142 sequences (28.2%) and a shared haplotype with 13 geographic populations. H9 was positioned centrally on the network map, while the rest of the shared haplotypes were clustered around their respective centers (Supplementary Figure S4A). H1 was the most frequent haplotype detected in the 16S rRNA sequences, represented by 134 from all sequences (94.4%), while the remaining shared haplotypes formed small cluster centers (Supplementary Figure S4B). H1 was also the most frequent haplotype of cox1 sequences in R. sanguineus (R. linnaei) ticks, represented by 46 from 52 sequences (88.5%), while the remaining shared haplotypes formed small cluster centers (Supplementary Figure S4C). The high percentage of haplotype frequency in the 16S rRNA sequences of R. microplus, and cox1 sequences of R. sanguineus (R. linnaei), suggested that these populations were relatively stable and capable of adapting to a variety of environments. CLHK was the locality of R. microplus ticks, with the highest number of haplotypes (Hn = 8 for cox1 and Hn = 3 for 16S rRNA sequences), while the locality of BWJ, XTDZ and QZ had an equal number of haplotypes (Hn = 2) for R. sanguineus (R. linnaei) cox1 sequences.

4. Discussion

The R. microplus species had the highest prevalence in this study and is found in many parts of the world, including South and Central America, Africa, Asia and Australia [31,32]. It is a major livestock pest, particularly affecting cattle. The R. microplus complex, based on cox1 and 16S rRNA genes, includes 5 and 4 distinct taxa, respectively, including R. annulatus, R. australis and R. microplus (clade A–C for the cox1 gene and clade A–B for the 16S rRNA gene) [15,31,33]. However, the cox1 gene has a greater intraspecies resolution within the R. microplus complex when compared to that of the internal transcribed spacer 2 (ITS2), 16S and 12S genes [15,33]. A previous study revealed that, apart from R. australis, all the other taxa were detected in China [20]. However, in this study the R. microplus ticks in Hainan Island were clustered within clade A and B, based on cox1 and 16S rRNA genes, respectively, together with isolates from Asia, Africa and South America [32]. R. sanguineus complex, based on molecular markers, is classified into four potential lineages, including temperate (R. sanguineus s.s.), southeastern Europe, tropical (R. sanguineus s.l.) and Afrotropical (R. afranicus) [34,35,36,37]. Šlapeta et al. (2021) [37] confirmed that Rhipicephalus linnaei, which is synonymous with R. sanguineus, was regarded as the oldest name and suggested that its adoption should refer to R. sanguineus sensu lato as the “tropical lineage”. In addition, Šlapeta et al. (2022) [38] elucidated further regarding the recognition of the tropical lineage (R. sanguineus s.l.) and officially designated it as Rhipicephalus linnaei. In this study, the phylogenetic analysis of cox1 and 16S rRNA genes identified R. sanguineus in Hainan Island as being clustered into one clade of tropical lineage in the temperate zone of a country (Thailand, China, India, Cuba, Brazil, Iraq, Angola, Pakistan and Sub-Saharan Africa), similar to previous reports [15,39,40]. A study by Zemtsova et al., 2016 [41] revealed that R. sanguineus of tropical lineage occurs in geographical areas with an annual mean temperature of greater than 20 °C, whereas the temperate lineage presents in geographical areas with an annual mean temperature of less than 20 °C. This supports the presence of R. sanguineus in this study site, where the average temperature ranges from 23 to 26 °C. The genetic diversity of R. haemaphysaloides, based on the cox1 gene, had not been investigated previously. Three genetic groups were recognized, based on this gene, and R. haemaphysaloides from this study was found to belong to one of the three. The result for the 16S rRNA gene was similar to that in a previous study by Li et al. (2018) [20], who recognized three genetic groups, from which R. haemaphysaloides was found in group one that belonged to the isolated tick from China. This study found two tick samples in the genus Haemaphysalis in Changliu town, Haikou district, namely H. cornigera and H. mageshimaensis, which were collected from cattle. A review of the geographical distribution of ticks in China found that H. cornigera was distributed in Fujian, Taiwan, Guangdong, Guangxi and Hainan, and the natural host of this species was buffaloes and cattle [10]. Moreover, a report by Doi et al. (2020) [42] revealed that H. cornigera is the dominant tick in wildlife (sika deer) on Niijima Island, Japan. The natural hosts for H. mageshimaensis are livestock, wild mammals and birds, which have been reported in China [10,43] and Japan [44]. Although there is little information on nucleotide sequence in BLAST searches and limited determination of evolutionary relationships of these species, the phylogenetic analysis in this study revealed that both H. cornigera and H. mageshimaensis were clustered and related together with isolates from China and Japan, respectively, with the Chinese isolates based on cox1 and 16S rRNA genes.
Moreover, this study also investigated the genetic diversity of different tick species in local areas. The cox1 gene marker presents greater resolution than the 16S rRNA gene marker in revealing the genetic variability of R. microplus and R. sanguineus (R. linnaei) in Hainan Island. From the 16S rRNA gene sequences, only six haplotypes were identified from R. microplus ticks in this study, while no variation site in R. sanguineus (R. linnaei) was observed, and only one haplotype generated. On the other hand, cox1 gene sequences revealed 26 and 4 haplotypes of R. microplus and R. sanguineus (R. linnaei), respectively. However, only R. microplus ticks generated a higher haplotype diversity, based on the cox1 gene. Similar to previous studies, several COI (cox1) haplotypes were revealed with high genetic differences in R. microplus [15], and the haplotype frequency in R. microplus was higher than in R. sanguineus [15,16]. Low et al. (2018) [45] revealed four COI haplotypes with slightly higher variation in R. sanguineus in Malaysia, but 16S rRNA sequences were highly conserved, which was similar to R. sanguineus (R. linnaei) in Hainan Island in this study. In the median-joining haplotype network, H1 was the most frequent haplotype for R. microplus and R. sanguineus (R. linnaei), based on 16S rRNA and cox1 sequences, respectively, with only a few shared haplotypes of all geographic populations in Hainan Island. This suggested that haplotype H1 is an ancestor of R. sanguineus (R. linnaei) and R. microplus and other haplotypes were shared over time. In contrast, the cox1 haplotype of R. microplus from seventeen sampled populations had high genetic diversity and revealed highly abundant shared haplotypes. It is indicated that ticks may have the ability to adapt to different environments and there was frequent gene communication between individuals or populations [46]. However, the remaining three tick species from this study were not investigated, due to the smaller number of tick specimens and species being represented by only one population.

5. Conclusions

This study employed morphological and molecular approaches to characterize and identify the genetic relationships of field-collected ticks from eight hundred and fifty-eight hard ticks collected from cattle, dogs and goats in 24 locations around Hainan Island. Based on morphological features, five tick species were identified, namely Rhipicephalus microplus, R. sanguineus (R. linnaei), R. haemaphysaloides, Haemaphysalis cornigera and H. mageshimaensis and confirmed subsequently by molecular analysis of cox1 and 16S rRNA genes. It was concluded that cox1 and 16S rRNA genes were useful markers for verifying species identification of hard ticks, according to analyses of the sequence data obtained. The analyses established a dependable DNA reference database that could be utilized in forensic entomology not only in China but also in other countries where these species are present.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14081592/s1, Figure S1: Morphological features of the adult female; Figure S2: Morphological features of the adult male; Figure S3: Morphological features of the female nymph; Figure S4: Median-joining haplotype network based on cox1 gene and 16S rRNA gene sequences of R. microplus and cox1 gene sequences of R. sanguineus (R. linnaei); Table S1: Tick specimens collected between July and December 2022 in Hainan Island; Table S2: Details of primers and thermal sequences used in PCR assays for tick species; Table S3: Comparison of distinguishing morphological features of R. microplus, R. sanguineus (R. linnaei), R. haemaphysaloides, H. cornigera and H. mageshimaensis in this study; Table S4: Detailed descriptions of the ticks analyzed in the current study; Table S5: Genetic distances between R. microplus and other tick species of cox1 sequences were determined using the Kimura 2-Parameter model; Table S6: Genetic distances between R. microplus and other tick species of 16S rRNA sequences were determined using the Kimura 2-Parameter model; Table S7: Genetic distances between R. sanguineus and other tick species of cox1 sequences were determined using the Kimura 2-Parameter model; Table S8: Genetic distances between R. sanguineus and other tick species of 16S rRNA sequences were determined using the Kimura 2-Parameter model; Table S9: Genetic distances between R. haemaphysaloides and other tick species of cox1 sequences were determined using the Kimura 2-Parameter model; Table S10: Genetic distances between R. haemaphysaloides and other tick species of 16S rRNA sequences were determined using the Kimura 2-Parameter model; Table S11: Genetic distances between Haemaphysalis spp. in present study and other tick species of cox1 sequences were determined using the Kimura 2-Parameter model; Table S12: Genetic distances between Haemaphysalis spp. in present study and other tick species of 16S rRNA sequences were determined using the Kimura 2-Parameter model; Table S13: Number of haplotypes (Hn), haplotype diversity (Hd) and nucleotide diversity (π), based on cox1 and 16S rRNA gene sequences of tick populations in Hainan Island; Table S14: Pairwise FST values among populations of R. microplus, based on data of the cox1 gene; Table S15: Distribution of cox1 and 16S rRNA haplotype sequences obtained in this study across sample sites.

Author Contributions

Conceptualization, J.I., X.L., Q.H., Z.-Y.L. and T.C.; methodology, J.I., X.L., Q.H., Z.-Y.L. and T.C.; software, J.I.; formal analysis, J.I., X.L., Q.H., Z.-Y.L. and T.C; data curation, J.I., X.L. and Q.H.; writing—original draft preparation, J.I.; writing—review and editing, Q.H., Z.-Y.L. and T.C.; project administration, T.C.; funding acquisition, J.I., Q.H., Z.-Y.L. and T.C. All authors have read and agreed to the published version of the manuscript.

Funding

The research project was supported by the National Natural Science Foundation of China (grant number 81060096 and 81560202 to T.C.), Hainan Natural Science Foundation Innovative Research Team Projects (grant number 2016CXTD010 to T.C.), High-level Talents of the Hainan Natural Science Foundation (grant number 820RC759 to T.C.), Hainan Key Scientific Research Projects (grant number 14A110060 to T.C.), Key Educational Reform Projects of Hainan Province (grant number Hnjg2019ZD-15 to T.C.), Key R&D Plan Projects of Hainan Province (grant number ZDYF2021SHFZ091 to T.C.), Postdoctoral Fund of Hainan General Hospital (grant number 2022BSH01 to J.I.), Hainan Province Science and Technology Special Fund (grant number ZDKJ2021035 to Z.L. and Q.H.) and National Natural Science Foundation of China (grant number U22A20363 to Q.H.). The Construction Project of Hainan Province Clinical Medical Center and Science, Innovation Platform for Academicians of Hainan Province provided support for this work. The study design, data collection and analysis, data interpretation and manuscript writing were not influenced by the funding bodies.

Institutional Review Board Statement

The animal study was reviewed and approved by the Medical Ethics Committee of Hainan General Hospital [No.: Med-Eth-Re [2023] 42]. Written informed consent was obtained from the owners for the participation of their animals in this study.

Informed Consent Statement

Not applicable.

Data Availability Statement

The genetic information obtained in this study is currently being submitted to the National Center for Biotechnology Information (NCBI).

Acknowledgments

The authors would like to thank all of the participants in this study as well as the staff at the Hainan General Hospital and Hainan University for their support. We would also like to thank the Agriculture and Rural Bureau of Danzhou City, Ledong Li Autonomous County, Lingao County, Lingshui Li Autonomous County, Qionghai City, Sanya City, Wanning City, Wuzhishan City and Xiuying District, Hainan, China, as well as all veterinarians for facilitating and assisting in the collection of tick samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wikel, S.K. Tick modulation of host immunity: An important factor in pathogen transmission. Int. J. Parasitol. 1999, 29, 851–859. [Google Scholar] [CrossRef] [PubMed]
  2. Parola, P.; Raoult, D. Ticks and tickborne bacterial diseases in humans: An emerging infectious threat. Clin. Infect. Dis. 2001, 32, 897–928. [Google Scholar] [CrossRef] [Green Version]
  3. de la Fuente, J.; Estrada-Pena, A.; Venzal, J.M.; Kocan, K.M.; Sonenshine, D.E. Overview: Ticks as vectors of pathogens that cause disease in humans and animals. Front. Biosci. 2008, 13, 6938–6946. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Dantas-Torres, F.; Chomel, B.B.; Otranto, D. Ticks and tick-borne diseases: A One Health perspective. Trends Parasitol. 2012, 28, 437–446. [Google Scholar] [CrossRef]
  5. Kilpatrick, A.M.; Randolph, S.E. Drivers, dynamics, and control of emerging vector-borne zoonotic diseases. Lancet 2012, 380, 1946–1955. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Fang, L.Q.; Liu, K.; Li, X.L.; Liang, S.; Yang, Y.; Yao, H.W.; Sun, R.X.; Sun, Y.; Chen, W.J.; Zuo, S.Q.; et al. Emerging tick-borne infections in mainland China: An increasing public health threat. Lancet Infect. Dis. 2015, 15, 1467–1479. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, G.; Huang, J.; Gong, X. Aishajiang·Abula: Hard tick and tick-borne diseases. Xinjiang Xu Mu Ye 2010, 9, 48–50. [Google Scholar]
  8. Zhao, J.; Wang, H.; Wang, Y. Regional distribution profiles of tick-borne pathogens in China. Chin. J. Vector Biol. Control 2012, 5, 445–447. [Google Scholar]
  9. Bian, C.; Liu, C.Y.; Zhang, C.Y.; Wang, K.; Li, X.D. Analysis of multiple infections of tick-borne diseases such as forest encephalitis, Lyme disease and spotted fever. World Latest Med. Inf. 2017, 17, 137. [Google Scholar]
  10. Zhang, Y.K.; Zhang, X.Y.; Liu, J.Z. Ticks (Acari: Ixodoidea) in China: Geographical distribution, host diversity, and specificity. Arch. Insect Biochem. Physiol. 2019, 102, e21544. [Google Scholar] [CrossRef] [Green Version]
  11. Abdel-Shafy, S.; El Namaky, A.H.; Zeina, H.A.A.; Ghazy, A.A. Description of immature stages of the brown dog tick Rhipicephalus Sanguineus (Acari: Ixodidae) using scanning electron microscopy. Int. J. Vet. Sci. Res. 2015, 1, 12–30. [Google Scholar]
  12. Amrutha, B.M.; Kumar, K.G.A.; Kurbet, P.S.; Varghese, A.; Deepa, C.K.; Pradeep, R.K.; Nimisha, M.; Asaf, M.; Juliet, S.; Ravindran, R.; et al. Morphological and molecular characterization of Rhipicephalus microplus and Rhipicephalus annulatus from selected states of southern India. Ticks Tick. Borne Dis. 2023, 14, 102086. [Google Scholar] [CrossRef] [PubMed]
  13. Walker, A.; Bouattour, A.; Camicas, J.L.; Estrada-Peña, A.; Horak, I.; Latif, A.; Pegram, R.G.; Preston, P.M. Ticks of Domestic Animals in Africa: A Guide to Identification of Species; Bioscience Reports: Edinburgh, UK, 2003; pp. 1–227. [Google Scholar]
  14. Paguem, A.; Manchang, K.; Kamtsap, P.; Renz, A.; Schaper, S.; Dobler, G.; Bakkes, D.K.; Chitimia-Dobler, L. Ticks and rickettsiae associated with wild animals sold in bush meat markets in Cameroon. Pathogens 2023, 12, 348. [Google Scholar] [CrossRef] [PubMed]
  15. Low, V.L.; Tay, S.T.; Kho, K.L.; Koh, F.X.; Tan, T.K.; Lim, Y.A.; Ong, B.L.; Panchadcharam, C.; Norma-Rashid, Y.; Sofian-Azirun, M. Molecular characterisation of the tick Rhipicephalus microplus in Malaysia: New insights into the cryptic diversity and distinct genetic assemblages throughout the world. Parasites Vectors 2015, 8, 341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Tantrawatpan, C.; Vaisusuk, K.; Chatan, W.; Pilap, W.; Suksavate, W.; Andrews, R.H.; Petney, T.N.; Saijuntha, W. Genetic diversity and phylogenetic analyses of ixodid ticks infesting cattle in northeast Thailand: The discovery of Rhipicephalus microplus clade C and the rarely detected R. haemaphysaloides. Exp. Appl. Acarol. 2022, 86, 535–548. [Google Scholar] [CrossRef]
  17. Sándor, A.D.; Milchev, B.; Takács, N.; Kontschán, J.; Szekeres, S.; Hornok, S. Five ixodid tick species including two morphotypes of Rhipicephalus turanicus on nestlings of Eurasian eagle owl (Bubo bubo) from south-eastern Bulgaria. Parasites Vectors 2021, 14, 334. [Google Scholar] [CrossRef]
  18. Kamani, J. Molecular identification and genetic analysis of Rhipicephalus sanguineus sensu lato of dogs in Nigeria, West Africa. Exp. Appl. Acarol. 2021, 85, 277–289. [Google Scholar] [CrossRef]
  19. Hornok, S.; Kontschán, J.; Estrada-Peña, A.; de Mera, I.G.; Tomanović, S.; de la Fuente, J. Contributions to the morphology and phylogeny of the newly discovered bat tick species, Ixodes ariadnae in comparison with I. vespertilionis and I. simplex. Parasites Vectors 2015, 8, 47. [Google Scholar] [CrossRef] [Green Version]
  20. Li, L.H.; Zhang, Y.; Wang, J.Z.; Li, X.S.; Yin, S.Q.; Zhu, D.; Xue, J.B.; Li, S.G. High genetic diversity in hard ticks from a China-Myanmar border county. Parasites Vectors 2018, 11, 469. [Google Scholar] [CrossRef]
  21. Sun, D.W.; Wang, S.Q.; Zhao, W. Seasonal succession of ticks in small wild animals in Chengmai county of Hainan province. Chin. J. Public Health 2012, 28, 1620. [Google Scholar]
  22. Lin, L.H.; Bai, Z.J.; Cai, J.; Fang, M.Y. Investigation and analysis of epiparasitic ticks on domestic animals in Qiongzhong area, Hainan province. Chin. J. Vector Biol. Control 2008, 6, 544–545. [Google Scholar]
  23. Communicable Disease Center. Pictorial Keys to Arthropods, Reptiles, Birds, and Mammals of Public Health Significance; CDC: Doha, Qatar, 1966; pp. 1–196. [Google Scholar]
  24. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  25. Tavaré, S. Some probabilistic and statistical problems in the analysis of DNA sequences. Lect. Math. Life Sci. (Am. Math. Soc.) 1986, 17, 57–86. [Google Scholar]
  26. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  27. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  28. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  29. Bandelt, H.J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  30. Leigh, J.W.; Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  31. Roy, B.C.; Estrada-Peña, A.; Krücken, J.; Rehman, A.; Nijhof, A.M. Morphological and phylogenetic analyses of Rhipicephalus microplus ticks from Bangladesh, Pakistan and Myanmar. Ticks Tick Borne Dis. 2018, 9, 1069–1079. [Google Scholar] [CrossRef]
  32. Kanduma, E.G.; Emery, D.; Githaka, N.W.; Nguu, E.K.; Bishop, R.P.; Šlapeta, J. Molecular evidence confirms occurrence of Rhipicephalus microplus Clade A in Kenya and sub-Saharan Africa. Parasites Vectors 2020, 13, 432. [Google Scholar] [CrossRef]
  33. Burger, T.D.; Shao, R.; Barker, S.C. Phylogenetic analysis of mitochondrial genome sequences indicates that the cattle tick, Rhipicephalus (Boophilus) microplus, contains a cryptic species. Mol. Phylogenetics Evol. 2014, 76, 241–253. [Google Scholar] [CrossRef]
  34. Hornok, S.; Sándor, A.D.; Tomanović, S.; Beck, R.; D’Amico, G.; Kontschán, J.; Takács, N.; Görföl, T.; Bendjeddou, M.L.; Földvári, G.; et al. East and west separation of Rhipicephalus sanguineus mitochondrial lineages in the Mediterranean Basin. Parasites Vectors 2017, 10, 39. [Google Scholar] [CrossRef] [Green Version]
  35. Nava, S.; Beati, L.; Venzal, J.M.; Labruna, M.B.; Szabó, M.P.J.; Petney, T.; Saracho-Bottero, M.N.; Tarragona, E.L.; Dantas-Torres, F.; Silva, M.M.S.; et al. Rhipicephalus sanguineus (Latreille, 1806): Neotype designation, morphological re-description of all parasitic stages and molecular characterization. Ticks Tick Borne Dis. 2018, 9, 1573–1585. [Google Scholar] [CrossRef] [PubMed]
  36. Bakkes, D.K.; Chitimia-Dobler, L.; Matloa, D.; Oosthuysen, M.; Mumcuoglu, K.Y.; Mans, B.J.; Matthee, C.A. Integrative taxonomy and species delimitation of Rhipicephalus turanicus (Acari: Ixodida: Ixodidae). Int. J. Parasitol. 2020, 50, 577–594. [Google Scholar] [CrossRef]
  37. Šlapeta, J.; Chandra, S.; Halliday, B. The “tropical lineage” of the brown dog tick Rhipicephalus sanguineus sensu lato identified as Rhipicephalus linnaei (Audouin, 1826). Int. J. Parasitol. 2021, 51, 431–436. [Google Scholar] [CrossRef]
  38. Šlapeta, J.; Halliday, B.; Chandra, S.; Alanazi, A.D.; Abdel-Shafy, S. Rhipicephalus linnaei (Audouin, 1826) recognised as the “tropical lineage” of the brown dog tick Rhipicephalus sanguineus sensu lato: Neotype designation, redescription, and establishment of morphological and molecular reference. Ticks Tick Borne Dis. 2022, 13, 102024. [Google Scholar] [CrossRef]
  39. Dantas-Torres, F.; Latrofa, M.S.; Annoscia, G.; Giannelli, A.; Parisi, A.; Otranto, D. Morphological and genetic diversity of Rhipicephalus sanguineus sensu lato from the New and Old Worlds. Parasites Vectors 2013, 6, 213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Caetano, R.L.; Vizzoni, V.F.; Bitencourth, K.; Carriço, C.; Sato, T.P.; Pinto, Z.T.; De Oliveira, S.V.; Amorim, M.; Voloch, C.M.; Gazeta, G.S. Ultrastructural morphology and molecular analyses of tropical and temperate “Species” of Rhipicephalus sanguineus sensu lato (Acari: Ixodidae) in Brazil. J. Med. Entomol. 2017, 54, 1201–1212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Zemtsova, G.E.; Apanaskevich, D.A.; Reeves, W.K.; Hahn, M.; Snellgrove, A.; Levin, M.L. Phylogeography of Rhipicephalus sanguineus sensu lato and its relationships with climatic factors. Exp. Appl. Acarol. 2016, 69, 191–203. [Google Scholar] [CrossRef] [Green Version]
  42. Doi, K.; Nishida, K.; Kato, T.; Hayama, S.I. Effects of introduced sika deer (Cervus nippon) and population control activity on the distribution of Haemaphysalis ticks in an island environment. Int. J. Parasitol. Parasites Wildl. 2020, 11, 302–307. [Google Scholar] [CrossRef]
  43. Tian, J.; Hou, X.; Ge, M.; Xu, H.; Yu, B.; Liu, J.; Shao, R.; Holmes, E.C.; Lei, C.; Shi, M. The diversity and evolutionary relationships of ticks and tick-borne bacteria collected in China. Parasites Vectors 2022, 15, 352. [Google Scholar] [CrossRef] [PubMed]
  44. Tsunoda, T.; Tatsuzawa, S. Questing height of nymphs of the bush tick, Haemaphysalis longicornis, and its closely related species, H. mageshimaensis: Correlation with body size of the host. Parasitology 2004, 128, 503–509. [Google Scholar] [CrossRef] [PubMed]
  45. Low, V.L.; Prakash, B.K. First genetic characterization of the brown dog tick Rhipicephalus sanguineus sensu lato in Peninsular Malaysia. Exp. Appl. Acarol. 2018, 75, 299–307. [Google Scholar] [CrossRef]
  46. Gui, Z.; Wu, L.; Cai, H.; Mu, L.; Yu, J.F.; Fu, S.Y.; Si, X.Y. Genetic diversity analysis of Dermacentor nuttalli within Inner Mongolia, China. Parasites Vectors 2021, 14, 131. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Collection sites of Rhipicephalus and Haemaphysalis ticks in Hainan Island, with tick numbers and species recorded at each study area. BWJ: Beiwuji, Meilan District; CLHK: Changliu Town, Xiuying District; CM: Chengmai County; GLYHK: Guilinyang, Meilan District; LD: Ledong Li Autonomous County; LS: Lingshui Li Autonomous County; MAHK: Meian Town, Xiuying District; MLCM: Meilang Village, Chengmai County; MYCM: Meiyang Village, County; PTHK: Potousan Village, Meilan District; QBHK: Qibi Village, Meilan District; QH: Qionghai City; QZ: Qiongzhong Li and Miao Autonomous County; SCDZ: Shancun Town, Danzhou City; SY: Sanya City; WN: Wanning City; WZS: Wuzhishan City; XHDZ: Xihuacha, Danzhou City; XTDZ: Xitian Village, Danzhou City; XYHK: Xiangtang Village, Xiuying District; ZHDZ: Zhonghe Town, Danzhou City; ZZLS: Zaozai Town, Lingshui Li Autonomous County.
Figure 1. Collection sites of Rhipicephalus and Haemaphysalis ticks in Hainan Island, with tick numbers and species recorded at each study area. BWJ: Beiwuji, Meilan District; CLHK: Changliu Town, Xiuying District; CM: Chengmai County; GLYHK: Guilinyang, Meilan District; LD: Ledong Li Autonomous County; LS: Lingshui Li Autonomous County; MAHK: Meian Town, Xiuying District; MLCM: Meilang Village, Chengmai County; MYCM: Meiyang Village, County; PTHK: Potousan Village, Meilan District; QBHK: Qibi Village, Meilan District; QH: Qionghai City; QZ: Qiongzhong Li and Miao Autonomous County; SCDZ: Shancun Town, Danzhou City; SY: Sanya City; WN: Wanning City; WZS: Wuzhishan City; XHDZ: Xihuacha, Danzhou City; XTDZ: Xitian Village, Danzhou City; XYHK: Xiangtang Village, Xiuying District; ZHDZ: Zhonghe Town, Danzhou City; ZZLS: Zaozai Town, Lingshui Li Autonomous County.
Genes 14 01592 g001
Figure 2. Morphological features of the capitulum of the adult female. R. microplus (A,D), R. sanguineus (R. linnaei) (B,E), R. haemaphysaloides (C,F), H. cornigera (G,I) and H. mageshimaensis (H,J) ticks collected in Hainan Island. (AC,G,H): Dorsal view; (DF,I,J): Ventral view. Bc, basis capitulum; Cap, capitulum; Ch, chelicerae; Hyp, hypostome; P, palps; Pa, porose area.
Figure 2. Morphological features of the capitulum of the adult female. R. microplus (A,D), R. sanguineus (R. linnaei) (B,E), R. haemaphysaloides (C,F), H. cornigera (G,I) and H. mageshimaensis (H,J) ticks collected in Hainan Island. (AC,G,H): Dorsal view; (DF,I,J): Ventral view. Bc, basis capitulum; Cap, capitulum; Ch, chelicerae; Hyp, hypostome; P, palps; Pa, porose area.
Genes 14 01592 g002
Figure 3. Neighbor-joining (NJ) phylogenetic tree of R. microplus based on the cox1 gene sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Figure 3. Neighbor-joining (NJ) phylogenetic tree of R. microplus based on the cox1 gene sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Genes 14 01592 g003
Figure 4. Neighbor-joining (NJ) phylogenetic tree of R. microplus based on the 16S rRNA gene sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Figure 4. Neighbor-joining (NJ) phylogenetic tree of R. microplus based on the 16S rRNA gene sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Genes 14 01592 g004
Figure 5. Neighbor-joining (NJ) phylogenetic tree of R. sanguineus based on the cox1 gene (A) and 16S rRNA gene (B) sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Figure 5. Neighbor-joining (NJ) phylogenetic tree of R. sanguineus based on the cox1 gene (A) and 16S rRNA gene (B) sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Genes 14 01592 g005
Figure 6. Maximum Likelihood (ML) phylogenetic tree of R. haemaphysaloides based on the cox1 gene (A) and 16S rRNA gene (B) sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Figure 6. Maximum Likelihood (ML) phylogenetic tree of R. haemaphysaloides based on the cox1 gene (A) and 16S rRNA gene (B) sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Genes 14 01592 g006
Figure 7. Maximum Likelihood (ML) phylogenetic tree of H. cornigera and H. mageshimaensis based on the cox1 gene (A) and 16S rRNA gene (B) sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Figure 7. Maximum Likelihood (ML) phylogenetic tree of H. cornigera and H. mageshimaensis based on the cox1 gene (A) and 16S rRNA gene (B) sequences. Evolutionary analyses were conducted in MEGA XI. Bootstrap values (1000 replications) are shown on the branches. Sequences generated from this study are in red label.
Genes 14 01592 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Intirach, J.; Lv, X.; Han, Q.; Lv, Z.-Y.; Chen, T. Morphological and Molecular Identification of Hard Ticks in Hainan Island, China. Genes 2023, 14, 1592. https://doi.org/10.3390/genes14081592

AMA Style

Intirach J, Lv X, Han Q, Lv Z-Y, Chen T. Morphological and Molecular Identification of Hard Ticks in Hainan Island, China. Genes. 2023; 14(8):1592. https://doi.org/10.3390/genes14081592

Chicago/Turabian Style

Intirach, Jitrawadee, Xin Lv, Qian Han, Zhi-Yue Lv, and Tao Chen. 2023. "Morphological and Molecular Identification of Hard Ticks in Hainan Island, China" Genes 14, no. 8: 1592. https://doi.org/10.3390/genes14081592

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