The Complete Mitochondrial Genome and Expression Profile of Mitochondrial Protein-Coding Genes in the Bisexual and Parthenogenetic Haemaphysalis longicornis

Comparative Analysis of Bisexual and Parthenogenetic Populations in Haemaphysalis Longicornis

Haemaphysalis longicornis, a three-host tick with a wide host range, is widely distributed in different countries and regions. It stands out among ticks due to its unique feature of having both parthenogenetic and bisexual populations. Despite their morphological resemblance, the characteristics of the parthenogenetic population have been overlooked. In this comprehensive study, we systematically compared the similarities and differences between these two populations. Our investigation revealed that the parthenogenetic H. longicornis, widely distributed in China, was found in ten provinces, surpassing the previously reported distribution. Notably, individuals from the parthenogenetic population exhibited a prolonged blood-feeding duration during the larval and nymph stages compared to their bisexual counterparts. Additionally, the life cycle of the parthenogenetic population was observed to be longer. A flow cytometry analysis indicated a DNA content ratio of approximately 2:3 between the bisexual and parthenogenetic populations. A phylogenetic analysis using whole mitochondrial genome sequences resulted in the separation of the phylogenetic tree into two distinct branches. A molecular analysis unveiled a consistent single T-base deletion at nucleotide 8497 in the parthenogenetic population compared to the bisexual population. Both populations displayed high viral infection capability and significant resistance to ivermectin. Intriguingly, despite these differences, the parthenogenetic population exhibited a similar life cycle to the bisexual population, retaining the ability to transmit pathogens such as Severe fever with thrombocytopenia syndrome virus (SFTSV) and Heartland Virus (HRTV). These findings contribute to a deeper understanding of the distinct characteristics and similarities between different populations of H. longicornis, laying the foundation for future research in this field.

Characterization of the complete mitochondrial genomes of five hard ticks and phylogenetic implications

Ticks are blood-sucking ectoparasites with significant medical and veterinary importance, capable of transmitting bacteria, protozoa, fungi, and viruses that cause a variety of human and animal diseases worldwide. In the present study, we sequenced the complete mitochondrial (mt) genomes of five hard tick species and analyzed features of their gene contents and genome organizations. The complete mt genomes of Haemaphysalis verticalis, H. flava, H. longicornis, Rhipicephalus sanguineus and Hyalomma asiaticum were 14855 bp, 14689 bp, 14693 bp, 14715 bp and 14722 bp in size, respectively. Their gene contents and arrangements are the same as those of most species of metastriate Ixodida, but distinct from species of genus Ixodes. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes with two different computational algorithms (Bayesian inference and maximum likelihood) revealed the monophylies of the genera Rhipicephalus, Ixodes and Amblyomma, however, rejected the monophyly of the genus Haemaphysalis. To our knowledge, this is the first report of the complete mt genome of H. verticalis. These datasets provide useful mtDNA markers for further studies of the identification and classification of hard ticks.

Complete Mitogenomes of Ticks Ixodes acutitarsus and Ixodes ovatus Parasitizing Giant Panda: Deep Insights into the Comparative Mitogenomic and Phylogenetic Relationship of Ixodidae Species

Ticks rank second in the world as vectors of disease. Tick infestation is one of the factors threatening the health and survival of giant pandas. Here, we describe the mitogenomes of Ixodes acutitarsus and Ixodes ovatus parasitizing giant pandas, and perform comparative and phylogenetic genomic analyses on the newly sequenced and other available mitogenomes of hard ticks. All six newly determined mitogenomes contain a typical gene component and share an ancient Arthropoda gene arrangement pattern. Our study suggests that I. ovatus is a species complex with high genetic divergence, indicating that different clades of I. ovatus represent distinct species. Comparative mitogenomic analyses show that the average A + T content of Ixodidae mitogenomes is 78.08%, their GC-skews are strongly negative, while AT-skews fluctuate around 0. A large number of microsatellites are detected in Ixodidae mitogenomes, and the main microsatellite motifs are mononucleotide A and trinucleotide AAT. We summarize five gene arrangement types, and identify the trnY-COX1-trnS1-COX2-trnK-ATP8-ATP6-COX3-trnG fragment is the most conserved region, whereas the region near the control region is the rearrangement hotspot in Ixodidae mitogenomes. The phylogenetic trees based on 15 genes provide a very convincing relationship (Ixodes + (Robertsicus + ((Bothriocroton + Haemaphysalis) + (Amblyomma + (Dermacentor + (Rhipicentor + (Hyalomma + Rhipicephalus))))))) with very strong supports. Remarkably, Archaeocroton sphenodonti is embedded in the Haemaphysalis clade with strong supports, resulting in paraphyly of the Haemaphysalis genus, so in-depth morphological and molecular studies are essential to determine the taxonomic status of A. sphenodonti and its closely related species. Our results provide new insights into the molecular phylogeny and evolution of hard ticks, as well as basic data for population genetics assessment and efficient surveillance and control for the giant panda-infesting ticks.

The Ovarian Development Genes of Bisexual and Parthenogenetic Haemaphysalis longicornis Evaluated by Transcriptomics and Proteomics

The tick Haemaphysalis longicornis has two reproductive groups: a bisexual group (HLBP) and a parthenogenetic group (HLPP). The comparative molecular regulation of ovarian development in these two groups is unexplored. We conducted transcriptome sequencing and quantitative proteomics on the ovaries of HLBP and HLPP, in different feeding stages, to evaluate the molecular function of genes associated with ovarian development. The ovarian tissues of HLBP and HLPP were divided into three feeding stages (early-fed, partially-fed and engorged). A total of 87,233 genes and 2,833 proteins were annotated in the ovary of H. longicornis in the different feeding stages. The differentially expressed genes (DEGs) of functional pathway analysis indicated that Lysosome, MAPK Signaling Pathway, Phagosome, Regulation of Actin Cytoskeleton, Endocytosis, Apoptosis, Insulin Signaling Pathway, Oxidative Phosphorylation, and Sphingolipid Metabolism were most abundant in the ovary of H. longicornis in the different feeding stages. Comparing the DEGs between HLBP and HLPP revealed that the ABC Transporter, PI3K-Akt Signaling Pathway and cAMP Signaling Pathway were the most enriched and suggested that the functions of signal transduction mechanisms may have changed during ovarian development. The functions of the annotated proteome of ovarian tissues were strongly correlated with the transcriptome annotation results, and these were further validated using quantitative polymerase chain reaction (qPCR). In the HLBP, the expression of cathepsin L, secreted proteins and glycosidase proteins was significantly up-regulated during feeding stages. In the HLPP, the lysozyme, yolk proteins, heat shock protein, glutathione S transferase, myosin and ATP synthase proteins were up-regulated during feeding stages. The significant differences of the gene expression between HLBP and HLPP indicated that variations in the genetic background and molecular function might exist in the two groups. These results provide a foundation for understanding the molecular mechanism and exploring the functions of genes in the ovarian development of different reproductive groups of H. longicornis.

Complete mitochondrial genome of Asian longhorned tick, Haemaphysalis longicornis, Neumann, 1901 (Acari: Ixodida: Ixodidae) identified in the United States

Haemaphysalis longicornis (Ixodida: Ixodidae), the Asian longhorned tick, which is native to temperate East Asia, has been recently detected in the northeastern region of the United States, drawing concerns about its potential impact on the US animal and public health sectors. Knowledge about the genetic features of H. longicornis found in the US is limited. Therefore, we sequenced the complete mitochondrial genome (mt-genome) from two H. longicornis ticks recently collected in the State of New York, USA, in 2020. These ticks were morphologically identified and tested for tick-borne pathogens at the Connecticut Veterinary Medical Diagnostic Laboratory (Storrs, CT). The mt-genome was 14,694 bp in length and encoded 37 genes, including 13 protein-coding genes, 22 transfer RNAs, and two ribosomal RNAs. Phylogenetic analysis showed that the mt-genome clustered with those of other H. longicornis identified in China. The mt-genome sequence was 99.7% identical to a H. longicornis mt-genome (GenBank: MK439888) collected in China. The cox1 gene haplotype in these ticks belonged to the H1 type, which is the dominant haplotype present in central NJ and Staten Island, NY. The complete mt-genome data are needed to provide insights into genetic changes and phylogenetic studies of H. longicornis ticks.

First glimpse into the origin and spread of the Asian longhorned tick, Haemaphysalis longicornis, in the United States

Established populations of Asian longhorned ticks (ALT), Haemaphysalis longicornis, were first identified in the United States (US) in 2017 by sequencing the mitochondrial cytochrome c oxidase subunit I (cox1) 'barcoding' locus followed by morphological confirmation. Subsequent investigations detected ALT infestations in 12, mostly eastern, US states. To gain information on the origin and spread of US ALT, we (1) sequenced cox1 from ALT populations across 9 US states and (2) obtained cox1 sequences from potential source populations [China, Japan and Republic of Korea (ROK) as well as Australia, New Zealand and the Kingdom of Tonga (KOT)] both by sequencing and by downloading publicly available sequences in NCBI GenBank. Additionally, we conducted epidemiological investigations of properties near its initial detection locale in Hunterdon County, NJ, as well as a broader risk analysis for importation of ectoparasites into the area. In eastern Asian populations (China/Japan/ROK), we detected 35 cox1 haplotypes that neatly clustered into two clades with known bisexual versus parthenogenetic phenotypes. In Australia/New Zealand/KOT, we detected 10 cox1 haplotypes all falling within the parthenogenetic cluster. In the United States, we detected three differentially distributed cox1 haplotypes from the parthenogenetic cluster, supporting phenotypic evidence that US ALT are parthenogenetic. While none of the source populations examined had all three US cox1 haplotypes, a phylogeographic network analysis supports a northeast Asian source for the US populations. Within the United States, epidemiological investigations indicate ALT can be moved long distances by human transport of animals, such as horses and dogs, with smaller scale movements on wildlife. These results have relevant implications for efforts aimed at minimizing the spread of ALT in the United States and preventing additional exotic tick introductions.

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.