Selection of reference genes for quantitative RT-PCR studies in Rhipicephalus (Boophilus) microplus and Rhipicephalus appendiculatus ticks and determination of the expression profile of Bm86

Towards modelling tick-virus interactions using the weakly pathogenic Sindbis virus: Evidence that ticks are competent vectors

Most tick-borne viruses (TBVs) are highly pathogenic and require high biosecurity, which severely limits their study. We found that Sindbis virus (SINV), predominantly transmitted by mosquitoes, can replicate in ticks and be subsequently transmitted, with the potential to serve as a model for studying tick-virus interactions. We found that both larval and nymphal stages of Rhipicephalus haemaphysaloides can be infected with SINV-wild-type (WT) when feeding on infected mice. SINV replicated in two species of ticks (R. haemaphysaloides and Hyalomma asiaticum) after infecting them by microinjection. Injection of ticks with SINV expressing enhanced Green Fluorescent Protein (eGFP) revealed that SINV-eGFP specifically aggregated in the tick midguts for replication. During blood-feeding, SINV-eGFP migrated from the midguts to the salivary glands and was transmitted to a new host. SINV infection caused changes in expression levels of tick genes related to immune responses, substance transport and metabolism, cell growth and death. SINV mainly induced autophagy during the early stage of infection; with increasing time of infection, the level of autophagy decreased, while the level of apoptosis increased. During the early stages of infection, the transcript levels of immune-related genes were significantly upregulated, and then decreased. In addition, SINV induced changes in the transcription levels of some functional genes that play important roles in the interactions between ticks and tick-borne pathogens. These results confirm that the SINV-based transmission model between ticks, viruses, and mammals can be widely used to unravel the interactions between ticks and viruses.

From the fat body to the hemolymph: Profiling tick immune and storage proteins through transcriptomics and proteomics

Ticks are blood-feeding arachnids that are known to transmit various pathogenic microorganisms to their hosts. During blood feeding, ticks activate their metabolism and immune system to efficiently utilise nutrients from the host's blood and complete the feeding process. In contrast to insects, in which the fat body is known to be a central organ that controls essential metabolic processes and immune defense mechanisms, the function of the fat body in tick physiology is still relatively unexplored. To fill this gap, we sought to uncover the repertoire of genes expressed in the fat body associated with trachea (FB/Tr) by analyzing the transcriptome of individual, partially fed (previtellogenic) Ixodes ricinus females. The resulting catalog of individual mRNA sequences reveals a broad repertoire of transcripts encoding proteins involved in nutrient storage and distribution, as well as components of the tick immune system. To gain a detailed insight into the secretory products of FB/Tr specifically involved in inter-tissue transport and humoral immunity, the transcriptomic data were complemented with the proteome of soluble proteins in the hemolymph of partially fed female ticks. Among these proteins, the hemolipoglyco-carrier proteins were predominant. When comparing immune peptides and proteins from the fat body with those produced by hemocytes, we found that the fat body serves as a unique producer of certain immune components. Finally, time-resolved transcriptional regulation of selected immune transcripts from the FB/Tr was examined in response to experimental challenges with model microbes and analyzed by RT-qPCR. Overall, our data show that the fat body of ticks, similar to insects, is an important metabolic tissue that also plays a remarkable role in immune defense against invading microbes. These findings improve our understanding of tick biology and its impact on the transmission of tick-borne pathogens.

Biochemical characterization of a novel sphingomyelinase-like protein from the Rhipicephalus microplus tick

Sphingomyelinase D is a toxin present in venomous spiders and bacteria and is associated with infection symptoms in patients affected by spider bites. It was observed that in Ixodes scapularis ticks, sphingomyelinase-like protein secreted in saliva can modulate the host immune response, affecting the transmission of flavivirus to the host via exosomes. In this work, a sphingomyelinase D-like protein (RmSMase) from R. microplus, a tick responsible for economic losses and a vector of pathogens for cattle, was investigated. The amino acid sequence revealed the lack of important residues for enzymatic activity, but the recombinant protein showed sphingomyelinase D activity. RmSMase shows Ca2+ and Mg2+ dependence in acidic pH, differing from IsSMase, which has Mg2+ dependence in neutral pH. Due to the difference between RmSMase and other SMases described, the data suggest that RmSMase belongs to SMase D class IIc. RmSMase mRNA transcription levels are upregulated during tick feeding, and the recombinant protein was recognized by host antibodies elicited after heavy tick infestation, indicating that RmSMase is present in tick saliva and may play a role in the tick feeding process.

Stable internal reference genes for quantitative RT-PCR analyses in Rhipicephalus microplus during embryogenesis

Studies on the transcriptional control of gene expression are crucial to understand changes in organism's physiological or cellular conditions. To obtain reliable data on mRNA amounts and the estimation of gene expression levels, it is crucial to normalize the target gene with one or more internal reference gene(s). However, the use of constitutive genes as reference genes is controversial, as their expression patterns are sometimes more complex than previously thought. In various arthropod vectors, including ticks, several constitutive genes have been identified by studying gene expression in different tissues and life stages. The cattle tick Rhipicephalus microplus is a major vector for several pathogens and is widely distributed in tropical and subtropical regions globally. Tick developmental physiology is an essential aspect of research, particularly embryogenesis, where many important developmental events occur, thus the identification of stable reference genes is essential for the interpretation of reliable gene expression data. This study aimed to identify and select R. microplus housekeeping genes and evaluate their stability during embryogenesis. Reference genes used as internal control in molecular assays were selected based on previous studies. These genes were screened by quantitative PCR (qPCR) and tested for gene expression stability during embryogenesis. Results demonstrated that the relative stability of reference genes varied at different time points during the embryogenesis. The GeNorm tool showed that elongation factor 1α (Elf1a) and ribosomal protein L4 (Rpl4) were the most stable genes, while H3 histone family 3A (Hist3A) and ribosomal protein S18 (RpS18) were the least stable. The NormFinder tool showed that Rpl4 was the most stable gene, while the ranking of Elf1a was intermediate in all tested conditions. The BestKeeper tool showed that Rpl4 and cyclophilin A (CycA) were the more and less stable genes, respectively. These data collectively demonstrate that Rpl4, Elf1a, and GAPDH are suitable internal controls for normalizing qPCR during R. microplus embryogenesis. These genes were consistently identified as the most stable in various analysis methods employed in this study. Thus, findings presented in this study offer valuable information for the study of gene expression during embryogenesis in R. microplus.

Mining gene expression data for rational identification of novel drug targets and vaccine candidates against the cattle tick, Rhipicephalus microplus

Control of complex parasites via vaccination remains challenging, with the current combination of vaccines and small drugs remaining the choice for an integrated control strategy. Studies conducted to date, are providing evidence that multicomponent vaccines will be needed for the development of protective vaccines against endo- and ectoparasites, though multicomponent vaccines require an in-depth understanding of parasite biology which remains insufficient for ticks. With the rapid development and spread of acaricide resistance in ticks, new targets for acaricide development also remains to be identified, along with novel targets that can be exploited for the design of lead compounds. In this study, we analysed the differential gene expression of Rhipicephalus microplus ticks that were fed on cattle vaccinated with a multi-component vaccine (Bm86 and 3 putative Bm86-binding proteins). The data was scrutinised for the identification of vaccine targets, small drug targets and novel pathways that can be evaluated in future studies. Limitations associated with targeting novel proteins for vaccine and/or drug design is also discussed and placed into the context of challenges arising when targeting large protein families and intracellular localised proteins. Lastly, this study provide insight into how Bm86-based vaccines may reduce successful uptake and digestion of the bloodmeal and overall tick fecundity.

More than Three Decades of Bm86: What We Know and Where to Go

Tick and tick-borne disease control have been a serious research focus for many decades. In a global climate of increasing acaricide resistance, host immunity against tick infestation has become a much-needed complementary strategy to common chemical control. From the earliest acquired resistance studies in small animal models to proof of concept in large production animals, it was the isolation, characterization, and final recombinant protein production of the midgut antigen Bm86 from the Australian cattle tick strain of Rhipicephalus (Boophilus) microplus (later reinstated as R. (B.) australis) that established tick subunit vaccines as a viable alternative in tick and tick-borne disease control. In the past 37 years, this antigen has spawned numerous tick subunit vaccines (either Bm86-based or novel), and though we are still describing its molecular structure and function, this antigen remains the gold standard for all tick vaccines. In this paper, advances in tick vaccine development over the past three decades are discussed alongside the development of biotechnology, where existing gaps and future directives in the field are highlighted.

Molecular identification and functional analysis of Niemann-Pick type C2 proteins，carriers for semiochemicals and other hydrophobic compounds in the brown dog tick, Rhipicephalus linnaei

Ticks are important vectors of many pathogens with tremendous impact on human and animal health. Studies of semiochemical interactions and mechanisms underlying chemoreception can provide important tools in tick management. Niemann-Pick type C2 (NPC2) proteins have been proposed as one type of chemoreceptor in arthropods. Here, we cloned two NPC2 genes in the brown dog tick, Rhipicephalus linnaei, the tropical lineage previously named R. sanguineus sensu lato and characterized them functionally. R.linNPC2a and R.linNPC2b genes were found to be expressed at each developmental stage with the highest level in adult males. By using quantitative real-time PCR we revealed expression in multiple tissues, including midgut, ovary, salivary glands and legs. Ligand binding analysis revealed that R.linNPC2b bound a wide spectrum of compounds, with β-ionone, α-amylcinnamaldehyde, 2-nitrophenol and benzaldehyde displaying the strongest binding affinity (Ki < 10 μM), whereas R.linNPC2a showed a more narrow ligand binding range, with intermediate binding affinity to α-amylcinnamaldehyde and 2-nitrophenol (Ki < 20 μM). Molecular docking indicated that the amino acid residue Phe89, Leu77 and Val131 of R.linNPC2a and Phe70, Leu132 and Phe73 of R.linNPC2b could bind multiple ligands. These residues might thus play a key role in the identification of the volatiles. Our results contribute to the understanding of olfactory mechanisms of R. linnaei and can offer new pathways towards new management strategies.

Identification and expression of short neuropeptide F and its receptors in the tick Ixodes ricinus

In Europe, the tick Ixodes ricinus is the most important vector of numerous pathogens that are transmitted during blood feeding on their vertebrate hosts. To elucidate mechanisms controlling blood intake and associated transmission of pathogens we identified and described expression of short neuropeptide F (sNPF) and its receptors which are known to regulate feeding in insects. Using in situ hybridization (ISH) and immunohistochemistry (IHC) we stained numerous neurons producing sNPF in the central nervous system (CNS; synganglion), while a few peripheral neurons were detected anteriorly to the synganglion, and on the surface of the hindgut and leg muscles. Apparent sNPF expression was also found in enteroendocrine cells individually scattered in anterior lobes of the midgut. In silico analyses and BLAST search for sNPF receptors revealed two putative G protein-coupled receptors (sNPFR1 and sNPFR2) in the I. ricinus genome. Aequorin-based functional assay in CHO cells showed that both receptors were specific and sensitive to sNPF in nanomolar concentrations. Increased expression levels of these receptors in the gut during blood intake suggest that sNPF signaling may be involved in regulation of feeding and digestion processes of I. ricinus.

Cymbopogon citratus and Citral Overcome Doxorubicin Resistance in Cancer Cells via Modulating the Drug's Metabolism, Toxicity, and Multidrug Transporters

Multidrug resistance (MDR) is the major complex mechanism that causes the failure of chemotherapy, especially with drugs of natural origin such as doxorubicin (DOX). Intracellular drug accumulation and detoxification are also involved in cancer resistance by reducing the susceptibility of cancer cells to death. This research aims to identify the volatile composition of Cymbopogon citratus (lemon grass; LG) essential oil and compare the ability of LG and its major compound, citral, to modulate MDR in resistant cell lines. The composition of LG essential oil was identified using gas chromatography mass spectrometry (GC-MS). In addition, a comparison of the modulatory effects of LG and citral, performed on breast (MCF-7/ADR), hepatic (HepG-2/ADR), and ovarian (SKOV-3/ADR) MDR cell lines, were compared to their parent sensitive cells using the MTT assay, ABC transporter function assays, and RT-PCR. Oxygenated monoterpenes (53.69%), sesquiterpene hydrocarbons (19.19%), and oxygenated sesquiterpenes (13.79%) made up the yield of LG essential oil. α-citral (18.50%), β-citral (10.15%), geranyl acetate (9.65%), ylangene (5.70), δ-elemene (5.38%), and eugenol (4.77) represent the major constituents of LG oil. LG and citral (20 μg/mL) synergistically increased DOX cytotoxicity and lowered DOX dosage by >3-fold and >1.5-fold, respectively. These combinations showed synergism in the isobologram and CI < 1. DOX accumulation or reversal experiment confirmed that LG and citral modulated the efflux pump function. Both substances significantly increased DOX accumulation in resistant cells compared to untreated cells and verapamil (the positive control). RT-PCR confirmed that LG and citral targeted metabolic molecules in resistant cells and significantly downregulated PXR, CYP3A4, GST, MDR1, MRP1, and PCRP genes. Our results suggest a novel dietary and therapeutic strategy combining LG and citral with DOX to overcome multidrug resistance in cancer cells. However, these results should be confirmed by additional animal experiments before being used in human clinical trials.

Identification and Characterization of Rhipicephalus microplus ATAQ Homolog from Haemaphysalis longicornis Ticks and Its Immunogenic Potential as an Anti-Tick Vaccine Candidate Molecule

Although vaccines are one of the environmentally friendly means to prevent the spread of ticks, there is currently no commercial vaccine effective against Haemaphysalis longicornis ticks. In this study, we identified, characterized, localized, and evaluated the expression patterns, and tested the immunogenic potential of a homologue of Rhipicephalus microplus ATAQ in H. longicornis (HlATAQ). HlATAQ was identified as a 654 amino acid-long protein present throughout the midgut and in Malpighian tubule cells and containing six full and one partial EGF-like domains. HlATAQ was genetically distant (homology < 50%) from previously reported ATAQ proteins and was expressed throughout tick life stages. Its expression steadily increased (p < 0.001) during feeding, reached a peak, and then decreased slightly with engorgement. Silencing of HlATAQ did not result in a phenotype that was significantly different from the control ticks. However, H. longicornis female ticks fed on a rabbit immunized with recombinant HlATAQ showed significantly longer blood-feeding periods, higher body weight at engorgement, higher egg mass, and longer pre-oviposition and egg hatching periods than control ticks. These findings indicate that the ATAQ protein plays a role in the blood-feeding-related physiological processes in the midgut and Malpighian tubules and antibodies directed against it may affect these tissues and disrupt tick engorgement and oviposition.

Function and evolution of the aquaporin IsAQP1 in the Lyme disease vector Ixodes scapularis

Ticks are important vectors of pathogenic viruses, bacteria, and protozoans to humans, wildlife, and domestic animals. Due to their life cycles, ticks face significant challenges related to water homeostasis. When blood-feeding, they must excrete water and ions, but when off-host (for stretches lasting several months), they must conserve water to avoid desiccation. Aquaporins (AQPs), a family of membrane-bound water channels, are key players in osmoregulation in many animals but remain poorly characterized in ticks. Here, we bioinformatically identified AQP-like genes from the deer tick Ixodes scapularis and used phylogenetic approaches to map the evolution of the aquaporin gene family in arthropods. Most arachnid AQP-like sequences (including those of I. scapularis) formed a monophyletic group clustered within aquaglycerolporins (GLPs) from bacteria to vertebrates. This gene family is absent from insects, revealing divergent evolutionary paths for AQPs in different hematophagous arthropods. Next, we sequenced the full-length cDNA of I. scapularis aquaporin 1 (IsAQP1) and expressed it heterologously in Xenopus oocytes to functionally characterize its permeability to water and solutes. Additionally, we examined IsAQP1 expression across different life stages and adult female organs. We found IsAQP1 is an efficient water channel with high expression in salivary glands prior to feeding, suggesting it plays a role in osmoregulation before or during blood feeding. Its functional properties are unique: unlike most GLPs, IsAQP1 has low glycerol permeability, and unlike most AQPs, it is insensitive to mercury. Together, our results suggest IsAQP1 plays an important role in tick water balance physiology and that it may hold promise as a target of novel vector control efforts.

Tick Diversity and Distribution of Hard (Ixodidae) Cattle Ticks in South Africa

Ticks are amongst the important ectoparasites where livestock are concerned, as they adversely affect the animals through bloodsucking. In tropical and subtropical countries, they transmit pathogens such as babesiosis, theileriosis, ehrlichiosis, and anaplasmosis in cattle, causing a reduction in production rate and significant concomitant economic losses. Ticks affect 80% of the cattle population across the world, with an estimated economic loss of USD 20–30 billion per year. In South Africa, economic losses in the livestock industry caused by ticks and tick-borne diseases are estimated to exceed USD 33 million per year (ZAR 500 million). There are seven major genera of ixodid ticks in Southern Africa (i.e., Amblyomma, Dermacentor, Haemaphysalis, Hyalomma, Ixodes and Rhipicephalus). The environment in which a tick lives is made up of all the various biological and abiotic factors that are either necessary or unnecessary for its life. The areas where various ticks have been found have been documented in many publications. Using these data, maps of possible species’ habitats can be made. Historical records on tick distribution may be incorrect due to identification mistakes or a change in the tick’s name. All the sources used to generate the maps for this review were unpublished and came from a wide range of sources. To identify tick species and the accompanying microbial ecosystems, researchers are increasingly adopting tick identification methods including 16S and 18S rDNA gene sequencing. Indeed, little is known about the genetic alterations that give important traits, including the predilection for tick hosts, transmission, and acaricide resistance. Opportunities for exploring these changes in tick populations and subpopulations are provided by advancements in omics technologies. The literature on the variety of ixodid ticks, their direct and indirect effects, and control methods in South Africa is compiled in this review.

Differential Expression of Immune Genes in the Rhipicephalus microplus Gut in Response to Theileria equi Infection

Rhipicephalus microplus is the only tick species known to serve as a biological vector of Theileria equi for horses and other equids in Brazil. The protozoan T. equi is one of the causal agents of equine piroplasmosis, a major threat in horse breeding systems. Vector competence is closely linked to the pathogens' ability to evade tick defense mechanisms. However, knowledge of tick immune response against infections by hemoparasites of the Theileria genus is scarce. In the present study, the expression of genes involved in immune signaling pathways of R. microplus adults' guts when challenged with a high or low parasitic load of T. equi was evaluated. This research demonstrates divergences in the immune gene expression pattern linked to T. equi infection in R. microplus since the Toll, IMD, and JNK signaling pathways were transcriptionally repressed in the guts of adult ticks infected with T. equi. Moreover, the results showed that different infectious doses of T. equi induce differential gene expression of key components of immune signaling cascades in R. microplus gut, suggesting a link between the intensity of infection and the activation of tick immunity response. The present study adds knowledge to elucidate the gut immune signaling response of R. microplus to T. equi infection. In addition, the generated data can serve as a basis for further investigations to develop strategies for controlling and preventing equine piroplasmosis.

Low Genetic Polymorphism in the Immunogenic Sequences of Rhipicephalus microplus Clade C

Rhipicephalus microplus tick highly affects the veterinary sector throughout the world. Different tick control methods have been adopted, and the identification of tick-derived highly immunogenic sequences for the development of an anti-tick vaccine has emerged as a successful alternate. This study aimed to characterize immunogenic sequences from R. microplus ticks prevalent in Pakistan. Ticks collected in the field were morphologically identified and subjected to DNA and RNA extraction. Ticks were molecularly identified based on the partial mitochondrial cytochrome C oxidase subunit (cox) sequence and screened for piroplasms (Theileria/Babesia spp.), Rickettsia spp., and Anaplasma spp. PCR-based pathogens-free R. microplus-derived cDNA was used for the amplification of full-length cysteine protease inhibitor (cystatin 2b), cathepsin L-like cysteine proteinase (cathepsin-L), glutathione S-transferase (GST), ferritin 1, 60S acidic ribosomal protein (P0), aquaporin 2, ATAQ, and R. microplus 05 antigen (Rm05Uy) coding sequences. The cox sequence revealed 100% identity with the nucleotide sequences of Pakistan's formerly reported R. microplus, and full-length immunogenic sequences revealed maximum identities to the most similar sequences reported from India, China, Cuba, USA, Brazil, Egypt, Mexico, Israel, and Uruguay. Low nonsynonymous polymorphisms were observed in ATAQ (1.5%), cathepsin-L (0.6%), and aquaporin 2 (0.4%) sequences compared to the homologous sequences from Mexico, India, and the USA, respectively. Based on the cox sequence, R. microplus was phylogenetically assembled in clade C, which includes R. microplus from Pakistan, Myanmar, Malaysia, Thailand, Bangladesh, and India. In the phylogenetic trees, the cystatin 2b, cathepsin-L, ferritin 1, and aquaporin 2 sequences were clustered with the most similar available sequences of R. microplus, P0 with R. microplus, R. sanguineus and R. haemaphysaloides, and GST, ATAQ, and Rm05Uy with R. microplus and R. annulatus. This is the first report on the molecular characterization of clade C R. microplus-derived immunogenic sequences.

Periviscerokinin (Cap2b; CAPA) receptor silencing in females of Rhipicephalus microplus reduces survival, weight and reproductive output

BACKGROUND

The cattle fever tick, Rhipicephalus (Boophilus) microplus, is a vector of pathogens causative of babesiosis and anaplasmosis, both highly lethal bovine diseases that affect cattle worldwide. In Ecdysozoa, neuropeptides and their G-protein-coupled receptors play a critical integrative role in the regulation of all physiological processes. However, the physiological activity of many neuropeptides is still unknown in ticks. Periviscerokinins (CAP2b/PVKs) are neuropeptides associated with myotropic and diuretic activities in insects. These peptides have been identified only in a few tick species, such as Ixodes ricinus, Ixodes scapularis and R. microplus, and their cognate receptor only characterized for the last two.

METHODS

Expression of the periviscerokinin receptor (Rhimi-CAP2bR) was investigated throughout the developmental stages of R. microplus and silenced by RNA interference (RNAi) in the females. In a first experiment, three double-stranded (ds) RNAs, named ds680-805, ds956-1109 and ds1102-1200, respectively, were tested in vivo. All three caused phenotypic effects, but only the last one was chosen for subsequent experiments. Resulting RNAi phenotypic variables were compared to those of negative controls, both non-injected and dsRNA beta-lactamase-injected ticks, and to positive controls injected with beta-actin dsRNA. Rhimi-CAP2bR silencing was verified by quantitative reverse-transcriptase PCR in whole females and dissected tissues.

RESULTS

Rhimi-CAP2bR transcript expression was detected throughout all developmental stages. Rhimi-CAP2bR silencing was associated with increased female mortality, decreased weight of surviving females and of egg masses, a delayed egg incubation period and decreased egg hatching (P < 0.05).

CONCLUSIONS

CAP2b/PVKs appear to be associated with the regulation of female feeding, reproduction and survival. Since the Rhimi-CAP2bR loss of function was detrimental to females, the discovery of antagonistic molecules of the CAP2b/PVK signaling system should cause similar effects. Our results point to this signaling system as a promising target for tick control.

Pyrokinin receptor silencing in females of the southern cattle tick Rhipicephalus (Boophilus) microplus is associated with a reproductive fitness cost

BACKGROUND

Rhipicephalus microplus is the vector of deadly cattle pathogens, especially Babesia spp., for which a recombinant vaccine is not available. Therefore, disease control depends on tick vector control. However, R. microplus populations worldwide have developed resistance to available acaricides, prompting the search for novel acaricide targets. G protein-coupled receptors (GPCRs) are involved in the regulation of many physiological processes and have been suggested as druggable targets for the control of arthropod vectors. Arthropod-specific signaling systems of small neuropeptides are being investigated for this purpose. The pyrokinin receptor (PKR) is a GPCR previously characterized in ticks. Myotropic activity of pyrokinins in feeding-related tissues of Rhipicephalus sanguineus and Ixodes scapularis was recently reported.

METHODS

The R. microplus pyrokinin receptor (Rhimi-PKR) was silenced through RNA interference (RNAi) in female ticks. To optimize RNAi, a dual-luciferase assay was applied to determine the silencing efficiency of two Rhimi-PKR double-stranded RNAs (dsRNA) prior to injecting dsRNA in ticks to be placed on cattle. Phenotypic variables of female ticks obtained at the endpoint of the RNAi experiment were compared to those of control female ticks (non-injected and beta-lactamase dsRNA-injected). Rhimi-PKR silencing was verified by quantitative reverse-transcriptase PCR in whole females and dissected tissues.

RESULTS

The Rhimi-PKR transcript was expressed in all developmental stages. Rhimi-PKR silencing was confirmed in whole ticks 4 days after injection, and in the tick carcass, ovary and synganglion 6 days after injection. Rhimi-PKR silencing was associated with an increased mortality and decreased weight of both surviving females and egg masses (P < 0.05). Delays in repletion, pre-oviposition and incubation periods were observed (P < 0.05).

CONCLUSIONS

Rhimi-PKR silencing negatively affected female reproductive fitness. The PKR appears to be directly or indirectly associated with the regulation of female feeding and/or reproductive output in R. microplus. Antagonists of the pyrokinin signaling system could be explored for tick control.

Myotropic Activities of Tick Pyrokinin Neuropeptides and Analog in Feeding Tissues of Hard Ticks (Ixodidae)

Neuropeptides regulate many important physiological processes in animals. The G protein-coupled receptors of corresponding small neuropeptide ligands are considered promising targets for controlling arthropod pests. Pyrokinins (PKs) are pleiotropic neuropeptides that, in some insect species, stimulate muscle contraction and modulate pheromone biosynthesis, embryonic diapause, and feeding behavior. However, their function remains unknown in ticks. In this study, we reported the myotropic activity of tick endogenous PKs and a PK agonist analog, PK-PEG₈ (MS[PEG₈]-YFTPRLa), on feeding tissues of two tick species representing the family Ixodidae lineages, namely, Prostriata (Ixodes scapularis) and Metastriata (Rhipicephalus sanguineus). First, we predicted the sequences of two periviscerokinins (PVK), one with a derived ending RNa and five PKs encoded by the CAPA peptide precursor from R. sanguineus and found the encoded PKs were identical to those of R. microplus identified previously. The pharynx-esophagus of both tick species responded with increased contractions to 10 μM of the endogenous PK as well as to PK-PEG₈ but not to the scrambled PK peptide, as expected. A dose-dependent myotropic activity of the PK-PEG₈ was found for both tick species, validating the analog activity previously found in the pyrokinin recombinant receptor assay. In agreement with the tissue activity elicited, we quantified the relative transcript abundance of R. sanguineus PK receptor in unfed female ticks and found it was the highest in the feeding tissues extracted from the capitulum and lowest in the reproductive tissue. This is the first report of the activity of pyrokinins in ticks. These findings strongly indicate the potential role of PKs in regulating tick blood feeding and therefore, making the tick PK receptor a potential target for interference.

The role of complement in the tick cellular immune defense against the entomopathogenic fungus Metarhizium robertsii

Entomopathogenic fungi (EPF) have been widely explored for their potential in the biological control of insect pests and as an environmentally friendly alternative to acaricides for limiting tick infestation in the field. The arthropod cuticle is the main barrier against fungal infection, however, an understanding of internal defense mechanisms after EPF intrusion into the invertebrate hemocoel is still rather limited. Using an infection model of the European Lyme borreliosis vector Ixodes ricinus with the EPF Metarhizium robertsii, we demonstrated that ticks are capable of protecting themselves to a certain extent against mild fungal infections. However, tick mortality dramatically increases when the capability of tick hemocytes to phagocytose fungal conidia is impaired. Using RNAi-mediated silencing of tick thioester-containing proteins (TEPs), followed by in vitro and/or in vivo phagocytic assays, we found that C3-like complement components and α2-macroglobulin pan-protease inhibitors secreted to the hemolymph play pivotal roles in M. robertsii phagocytosis.

REDOX IMBALANCE INDUCES REMODELING OF GLUCOSE METABOLISM IN RHIPICEPHALUS MICROPLUS EMBRYONIC CELL LINE

Carbohydrate metabolism not only functions in supplying cellular energy but also has an important role in maintaining physiological homeostasis and in preventing oxidative damage caused by reactive oxygen species. Previously, we showed that arthropod embryonic cell lines have high tolerance to H₂O₂ exposure. Here, we describe that Rhipicephalus microplus tick embryonic cell line (BME26) employs an adaptive glucose metabolism mechanism that confers tolerance to hydrogen peroxide at concentrations too high for other organisms. This adaptive mechanism sustained by glucose metabolism remodeling promotes cell survival and redox balance in BME26 cell line after millimolar H₂O₂ exposure. The present work shows that this tick cell line could tolerate high H₂O₂ concentrations by initiating a carbohydrate-related adaptive response. We demonstrate that gluconeogenesis was induced as a compensation strategy that involved, among other molecules, the metabolic enzymes NADP-ICDH, G6PDH, and PEPCK. We also found that this phenomenon was coupled to glycogen accumulation and glucose uptake, supporting the pentose phosphate pathway to sustain NADPH production and leading to cell survival and proliferation. Our findings suggest that the described response is not atypical, being also observed in cancer cells, which highlights the importance of this model to all proliferative cells. We propose that these results will be useful in generating basic biological information to support the development of new strategies for disease treatment and parasite control.

Characterization of a single-domain von Willebrand factor type C protein (HaSVC) from the salivary gland of the tick Hyalomma asiaticum

As a widely distributed arthropod and vector for various pathogens, Hyalomma asiaticum presents great risk and potential losses in animal husbandry. Effective measures, including the use of vaccines, are necessary for controlling ticks and tick-borne diseases. A concise understanding of the tick-host interaction associated molecules and pathways is required for vaccine development. In the present study, a protein containing a single-domain von Willebrand factor type C (HaSVC) was isolated from H. asiaticum and was subjected to functional identification. As a result, the full-length sequence of the HaSVC (506 bp) gene was obtained, which putatively encodes 100 amino acids with a predicted molecular mass of 11 kDa, excluding the 23-amino acid signal peptide. HaSVC contains 8 cysteines to form 4 disulfide bonds. The native HaSVC protein was detected in multiple tick organs. HaSVC neither attenuated the anti-coagulation process nor directly affected the blood feeding of adult ticks. However, the purified recombinant protein HaSVC (rHaSVC/GST) significantly increased the proliferation of mice spleen cells. This might suggest a regulatory function for HaSVC on inflammation, thus providing new information that may explain the "crosstalk" between ticks and hosts.

Haem-responsive gene transporter enables mobilization of host haem in ticks

Ticks, notorious blood-feeders and disease-vectors, have lost a part of their genetic complement encoding haem biosynthetic enzymes and are, therefore, dependent on the acquisition and distribution of host haem. Solute carrier protein SLC48A1, aka haem-responsive gene 1 protein (HRG1), has been implicated in haem transport, regulating the availability of intracellular haem. HRG1 transporter has been identified in both free-living and parasitic organisms ranging from unicellular kinetoplastids, nematodes, up to vertebrates. However, an HRG1 homologue in the arthropod lineage has not yet been identified. We have identified a single HRG1 homologue in the midgut transcriptome of the tick Ixodes ricinus, denoted as IrHRG, and have elucidated its role as a haem transporter. Data from haem biosynthesis-deficient yeast growth assays, systemic RNA interference and the evaluation of gallium protoporphyrin IX-mediated toxicity through tick membrane feeding clearly show that IrHRG is the bona fide tetrapyrrole transporter. We argue that during evolution, ticks profited from retaining a functional hrg1 gene in the genome because its protein product facilitates host haem escort from intracellularly digested haemoglobin, rendering haem bioavailable for a haem-dependent network of enzymes.

Identification of the Bcl-2 and Bax homologs from Rhipicephalus haemaphysaloides and their function in the degeneration of tick salivary glands

Background

The salivary glands of female ticks degenerate rapidly by apoptosis and autophagy after feeding. Bcl-2 family proteins play an important role in the apoptosis pathways, but the functions of these proteins in ticks are unclear. We studied Bcl-2 and Bax homologs from Rhipicephalus haemaphysaloides and determined their functions in the degeneration of the salivary glands.

Methods

Two molecules containing conserved BH (Bcl-2 family homology) domains were identified and named RhBcl-2 and RhBax. After protein purification and mouse immunization, specific polyclonal antibodies (PcAb) were created in response to the recombinant proteins. Reverse transcription quantitative PCR (RT-qPCR) and western blot were used to detect the presence of RhBcl-2 and RhBax in ticks. TUNEL assays were used to determine the level of apoptosis in the salivary glands of female ticks at different feeding times after gene silencing. Co-transfection and GST pull-down assays were used to identify interactions between RhBcl-2 and RhBax.

Results

The RT-qPCR assay revealed that RhBax gene transcription increased significantly during feeding at all tick developmental stages (engorged larvae, nymphs, and adult females). Transcriptional levels of RhBcl-2 and RhBax increased more significantly in the female salivary glands than in other tissues post engorgement. RhBcl-2 silencing significantly inhibited tick feeding. In contrast, RhBax interference had no effect on tick feeding. TUNEL staining showed that apoptosis levels were significantly reduced after interference with RhBcl-2 expression. Co-transfection and GST pull-down assays showed that RhBcl-2 and RhBax could interact but not combine in the absence of the BH3 domain.

Conclusions

This study identified the roles of RhBcl-2 and RhBax in tick salivary gland degeneration and finds that the BH3 domain is a key factor in their interactions.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04879-z.

Engorgement of Rhipicephalus haemaphysaloides ticks blocked by silencing a protein inhibitor of apoptosis

Inhibitors of apoptosis (IAPs) are regulators of cell death and may play a role in the salivary glands of ticks during blood-feeding. We cloned the open reading frame (ORF) sequence of the IAP gene in Rhipicephalus haemaphysaloides (RhIAP). The RhIAP ORF of 1887 bp encodes a predicted protein of 607 amino acids, which contains three baculovirus IAP repeat domains and a RING finger motif. A real-time PCR assay showed that RhIAP mRNA was expressed in all the tick developmental stages (eggs, larvae, nymphs, and adults) and in all tissues examined (midgut, ovary, salivary glands, fat body, and hemolymph). Western blot showed that the protein level of RhIAP in salivary glands increased during tick blood-feeding and decreased towards the end of tick engorgement. RhIAP gene silencing in vitro experiments with salivary glands demonstrated that RhIAP could be effectively knocked down within 48 h after dsRNA treatment, and as a consequence, salivary glands displayed apoptotic morphology. RhIAP gene silencing also inhibited tick blood-feeding and decreased the engorgement rate. These data suggest that RhIAP might be a suitable RNAi target for tick control.

Selection and Validation of Reference Genes for Gene Expression Analysis in Tuta absoluta Meyrick (Lepidoptera: Gelechiidae)

Simple Summary

Reference genes are critical for standardizing expression data of RT-qPCR across samples of organisms under different experimental conditions. However, most commonly used reference genes may not be stably expressed leading to a risk of misinterpretation of results. In our study, nine reference genes were evaluated in Tuta absoluta (a destructive pest of tomato) at different developmental stages, tissues, 20-hydroxyecdysone (20E) and insecticide treatments. Finally, the expression profile of indicator gene EcR after 20E treatment was evaluated to verify the accuracy of the results. This study is essential for improving accuracy and reliability to normalize gene expression data in T. absoluta and provides a useful strategy for other insects.

Abstract

The tomato leaf miner, Tuta absoluta is a destructive pest of tomato. The leaf-mining activities of its larvae can cause significant yield losses. Real-time quantitative polymerase chain reaction (RT-qPCR) is commonly used to measure gene expression, and the selection of stable reference genes for calibration and standardization is critical for accurate use of RT-qPCR. We studied the stable expression of nine common housekeeping genes in T. absoluta. These were examined at different developmental stages, in larval tissues, as well as those induced by exposure to 20E and insecticides. Four dedicated algorithms (geNorm, BestKeeper, NormFinder, and ΔCt method) and online tool (RefFinder) were used to analyze and rank the tested reference genes. Based on the standardized gene expression data of target gene ecdysone receptor (EcR), the applicability of specific reference genes was verified. The results clarify that the optimal internal reference genes vary greatly under different experimental conditions. GAPDH and RPS11 were the best reference genes for developmental stages; RPL28 and RPL10 for different tissues; EF1α and RPL28 for 20E treatment; EF1α and RPL7A for insecticide treatments. The most suitable reference genes in all experimental conditions are EF1α and RPL28.

Identification and Characterization of Immunodominant Proteins from Tick Tissue Extracts Inducing a Protective Immune Response against Ixodes ricinus in Cattle

Ixodes ricinus is the main vector of tick-borne diseases in Europe. An immunization trial of calves with soluble extracts of I. ricinus salivary glands (SGE) or midgut (ME) previously showed a strong response against subsequent tick challenge, resulting in diminished tick feeding success. Immune sera from these trials were used for the co-immunoprecipitation of tick tissue extracts, followed by LC-MS/MS analyses. This resulted in the identification of 46 immunodominant proteins that were differentially recognized by the serum of immunized calves. Some of these proteins had previously also drawn attention as potential anti-tick vaccine candidates using other approaches. Selected proteins were studied in more detail by measuring their relative expression in tick tissues and RNA interference (RNAi) studies. The strongest RNAi phenotypes were observed for MG6 (A0A147BXB7), a protein containing eight fibronectin type III domains predominantly expressed in tick midgut and ovaries of feeding females, and SG2 (A0A0K8RKT7), a glutathione-S-transferase that was found to be upregulated in all investigated tissues upon feeding. The results demonstrated that co-immunoprecipitation of tick proteins with host immune sera followed by protein identification using LC-MS/MS is a valid approach to identify antigen–antibody interactions, and could be integrated into anti-tick vaccine discovery pipelines.

Mialostatin, a Novel Midgut Cystatin from Ixodes ricinus Ticks: Crystal Structure and Regulation of Host Blood Digestion

The hard tick Ixodes ricinus is a vector of Lyme disease and tick-borne encephalitis. Host blood protein digestion, essential for tick development and reproduction, occurs in tick midgut digestive cells driven by cathepsin proteases. Little is known about the regulation of the digestive proteolytic machinery of I. ricinus. Here we characterize a novel cystatin-type protease inhibitor, mialostatin, from the I. ricinus midgut. Blood feeding rapidly induced mialostatin expression in the gut, which continued after tick detachment. Recombinant mialostatin inhibited a number of I. ricinus digestive cysteine cathepsins, with the greatest potency observed against cathepsin L isoforms, with which it co-localized in midgut digestive cells. The crystal structure of mialostatin was determined at 1.55 Å to explain its unique inhibitory specificity. Finally, mialostatin effectively blocked in vitro proteolysis of blood proteins by midgut cysteine cathepsins. Mialostatin is likely to be involved in the regulation of gut-associated proteolytic pathways, making midgut cystatins promising targets for tick control strategies.

Functional characterization of the insulin signaling pathway in the hard tick Ixodes ricinus

Ticks are blood-feeding arachnids transmitting a variety of pathogens to humans and animals. A unique trait in tick physiology is their ability to engorge and digest large amounts of host blood, ensuring their high reproductive potential. Activation of the blood digestive machinery in the tick gut, as well as processes controlling maturation of ovaries, are triggered upon blood meal uptake by still largely unknown mechanisms. Sensing of the nutritional status in metazoan organisms is facilitated by the evolutionarily conserved Insulin Signaling Pathway (ISP) and the interlinked Target of Rapamycin (TOR) pathway. Recently, we have identified three components of these pathways in the hard tick Ixodes ricinus midgut transcriptome, namely a putative insulin receptor (InR), and the downstream intracellular serine/threonine kinases AKT and TOR. In this study, we primarily focus on the molecular and functional characterization of the I. ricinus insulin receptor (IrInR), the first InR characterized in Chelicerates. A phylogenetic analysis across the major Arthropod lineages demonstrated that ticks possess only one gene encoding an InR-related molecule. Tissue expression profiling by quantitative PCR in semi-engorged I. ricinus females revealed that the IrInR, as well as AKT (IrAKT) and TOR (IrTOR) are expressed in various organs, with the highest expression being detected in ovaries. We have further evaluated the impact of RNAi-mediated knock-down (KD) of IrInR, IrAKT, and IrTOR on tick blood-feeding and reproductive capacity. Weights of engorged IrInR KD females and laid egg clutches were reduced compared to the control group, and these quantitative parameters clearly correlated with the efficiency of RNAi-KD achieved in individual ticks. The most striking phenotype was observed for IrAKT KD that impaired tick feeding and completely aborted egg production. A recombinant extracellular fragment of the IrInR α-subunit was used to produce antibodies in experimental rabbits to assess its potential as a protective antigen against tick feeding and reproduction. Our data clearly indicate the functionality of the ISP in ticks and demonstrate the need for further investigation of specific roles played by the endogenous insulin-like peptides in tick physiological processes.

Selection of reference genes for quantitative PCR analysis in poultry red mite (Dermanyssus gallinae)

Poultry red mites (PRMs, Dermanyssus gallinae) are harmful ectoparasites that affect farmed chickens and cause serious economic losses in the poultry industry worldwide. Acaricides are used for PRM control; however, some PRMs have developed acaricide-resistant properties, which have indicated the need for different approaches for PRM control. Therefore, it is necessary to elucidate the biological status of PRMs to develop alternative PRM control strategies. Quantitative polymerase chain reaction (qPCR) allows analysis of the biological status at the transcript level. However, reference genes are preferable for accurate comparison of expression level changes given the large variation in the quality of the PRM samples collected in each farm. This study aimed to identify candidate reference genes with stable expression levels in the different blood feeding states and life stages of PRMs. First, we selected candidates based on the following criteria: sufficient expression intensity and no significant expression difference between fed and starved states. We selected and characterized seven candidate reference genes. Among them, we evaluated the gene expression stability between the starved and fed states using RefFinder; moreover, we compared their expression levels in each life-stage and identified two reference genes, Elongation factor 1-alpha (ELF1A)-like and apolipophorins-like. Finally, we evaluated the utility of the candidates as reference genes, and the use of ELF1A-like and apolipophorins-like successfully normalized ATP synthase subunit g -like gene expression. Thus, ELF1A-like and apolipophorins-like could be suitable reference genes in PRMs.

Characterization of the Rhipicephalus (Boophilus) microplus Sialotranscriptome Profile in Response to Theileria equi Infection

This study intends to characterize the sialotranscriptome profile of Rhipicephalus (Boophilus) microplus in response to Theileria equi and identify genes of interest with differential genomic expression, indicating relevant targets in the tick–protozoan interactions. The experimental design consisted of RNA sequencing from uninfected and T. equi-infected R. microplus salivary glands (SGs) to obtain transcriptomic profiles for characterization and comparison. A total of 288,952 transcripts were obtained from both tick profiles, 3456 transcripts (p < 0.05) differentially expressed in response to T. equi infection. The uninfected SGs’ registered 231,179 transcripts, of which 155,359 were annotated. The most transcribed sequences were female-specific histamine binding protein and lipocalins. Regarding the T. equi-infected SGs, from the 238,964 assembled transcripts, 163,564 were annotated. The most transcribed sequences were histone demethylase JARID1 and Y-box-binding protein. Five transcripts (cystatin, arginase, nuclear factor κB kinase inhibitor subunit β (IκB), IκB delta, lysosomal-trafficking regulator, and reeler protein) presented the gene ontology (GO) category “response to protozoan” and were exclusively displayed in the T. equi-infected profile. The transcriptome of T. equi was also analyzed, registering 4728 hits. The study’s genetic and molecular information would be of great value for future studies and biotechnological applications envisaging disease control.

ATG5 is instrumental in the transition from autophagy to apoptosis during the degeneration of tick salivary glands

Female tick salivary glands undergo rapid degeneration several days post engorgement. This degeneration may be caused by the increased concentration of ecdysone in the hemolymph during the fast feeding period and both autophagy and apoptosis occur. In this work, we first proved autophagy-related gene (ATG) and caspase gene expression peaks during degeneration of the tick salivary glands. We explored the regulatory role of Rhipicephalus haemaphysaloides autophagy-related 5 (RhATG5) in the degeneration of tick salivary glands. During the fast feeding phase, RhATG5 was cleaved and both calcium concentration and the transcription of Rhcalpains increased in the salivary glands. Recombinant RhATG5 was cleaved by μ-calpain only in the presence of calcium; the mutant RhATG5^191-199Δ was not cleaved. Treatment with 20-hydroxyecdysone (20E) led to programmed cell death in the salivary glands of unfed ticks in vitro, RhATG8-phosphatidylethanolamine (PE) was upregulated in ticks treated with low concentration of 20E. Conversely, RhATG8-PE decreased and Rhcaspase-7 increased in ticks treated with a high concentration of 20E and transformed autophagy to apoptosis. High concentrations of 20E led to the cleavage of RhATG5. Calcium concentration and expression of Rhcalpains were also upregulated in the tick salivary glands. RNA interference (RNAi) of RhATG5 in vitro inhibited both autophagy and apoptosis of the tick salivary glands. RNAi of RhATG5 in vivo significantly inhibited the normal feeding process. These results demonstrated that high concentrations of 20E led to the cleavage of RhATG5 by increasing the concentration of calcium and stimulated the transition from autophagy to apoptosis.

Ticks are well-known pathogen vectors which transmitted virus, bacterial and protozoan. They are considered to be second only to mosquitoes as global vectors of human diseases. Most tick-borne pathogens (TBPs) are transmitted to hosts through tick bites assisted by saliva. Control of ticks has been achieved primarily by the application of acaricides, a method that has drawbacks such as environmental contamination and selection of pesticide-resistant ticks. Understanding the tick physiological characteristics is the key step for this objective; however, there are knowledge gap remained in tick physiology. Tick salivary glands rapidly degenerate and disappear within 4 days post engorgement. In this research, we are focused on tick salivary glands rapidly degeneration within 4 days post engorgement, and made several highlights findings: The first work demonstrated that 20E promotes both autophagy and apoptosis during tick salivary gland degeneration; RhATG5 is the first reported ATG5 homologue in ticks; RhATG5 play an important role in both autophagy and apoptosis during the degeneration of tick salivary glands.

Transcriptome analysis of the innate immune system of Hyalomma asiaticum

Ticks are considered to be the second most important vectors of human infectious diseases. The innate immune system is the key factor that affects its vector competence. Hyalomma asiaticum is the primary vector of Crimean-Congo hemorrhagic fever virus (CCHFV). However, the immune system of H. asiaticum remains virtually unknown. Here, a high throughput full-length mRNA sequencing method was adopted to define the immunotranscriptome of H. asiaticum infected with the fungal pathogen Beauveria bassiana and gram-negative bacterium Enterobacter cloacae. The analysis yielded 22,300 isoforms with an average length of 3233 bps. In total, 68 potential immunity-related genes were identified based on similarity to the homologs known to be involved in immunity. These included most members of the Toll and JAK/STAT signaling pathways, but not the IMD signaling pathway. Moreover, two copies of Dicer-2 and five copies of Argonaute-2 were detected. These genes are postulated to be involved in the RNA interference (RNAi) pathway, which is an important defense against RNA viruses. Overall, this study provides the foundation for understanding the immune response of H. asiaticum to CCHFV.

Sequence identification and expression profile of seven Dermacentor marginatus glutathione S-transferase genes

Dermacentor marginatus is a widespread tick species and a vector of many pathogens in Eurasia. Due to the medical importance of D. marginatus, control measures are needed for this tick species. Currently tick control approaches rely mostly on acaricide application, whereas wrong and irrational acaricide use may result in drug resistance and residue problems. Vaccination as an alternative approach for tick control has been proven to be effective towards some tick species. However, immunization against D. marginatus has not yet reached satisfactory protection. The effort of in silico based analysis could predict antigenicity and identify candidates for anti-tick vaccine development. We carried out an in silico analysis of D. marginatus glutathione S-transferases (DmGSTs) in order to identify blood-feeding induced GSTs as antigens that can be used in anti-tick vaccine development. Phylogenetic analysis, linear B-cell epitope prediction, homology modeling, and conformational B-cell epitope mapping on the GST models were performed to identify highly antigenic DmGSTs. Relative gene expressions of the seven GSTs were profiled through real-time quantitative PCR (RT-qPCR) to outline GSTs up-regulated during blood feeding. The phylogenetic analysis indicated that the seven GSTs belonged to four classes of GST, including one in epsilon-class, one in zeta-class, one in omega-class, and four in mu-class. Linear B-cell epitope prediction revealed mu-class GSTs share similar conserved antigenic regions. The conformational B-cell epitope mapped on the homology model of the GSTs displayed that GSTs of mu-class showed stronger antigenicity than that of other classes. RT-qPCR revealed DmGSTM1 and DmGSTM2 were positively related to blood feeding. In sum, the data suggest that DmGSTM1 and DmGSTM2 could be tested for potential anti-tick vaccine trials.

Molecular evidence confirms occurrence of Rhipicephalus microplus Clade A in Kenya and sub-Saharan Africa

Background

The tick vector Rhipicephalus microplus which transmits Babesia spp. and rickettsial pathogens has not been reported in Kenya since 1998. More recently, the pathogenic Babesia bovis has been detected in cattle blood DNA. The status of R. microplus in Kenya remains unknown. This study employed morphological and molecular tools to characterize R. microplus originating from Kenya and assess the genetic relationships between Kenyan and other African R. microplus genotypes.

Methods

Ticks were collected in south-eastern Kenya (Kwale County) from cattle and characterized to investigate the existence of R. microplus. Genetic and phylogenetic relationships between the Kenyan and other annotated R. microplus reference sequences was investigated by analysis of the cytochrome c oxidase subunit 1 (cox1) gene. To further characterize Kenyan ticks, we generated low coverage whole genome sequences of two R. microplus, one R. decoloratus and R. appendiculatus. A B. bovis specific TaqMan probe qPCR assay was used to detect B. bovis in gDNA from R. microplus ticks.

Results

Occurrence of R. microplus was confirmed in Kwale County, Kenya. The Kenyan R. microplus cox1 sequences showed very high pairwise identities (> 99%) and clustered very closely with reference African R. microplus sequences. We found a low genetic variation and lack of geographical sub-structuring among the African cox1 sequences of R. microplus. Four complete mitochondrial (mt) genomes for two R. microplus, one R. decoloratus and one R. appendiculatus were assembled from next generation sequence data. The mitochondrial genome sequences of the two Kenyan R. microplus ticks clustered closely with reference genome sequences from Brazil, USA, Cambodia and India forming R. microplus Clade A. No B. bovis was detected in the Kwale R. microplus DNA.

Conclusions

These findings confirm the presence of R. microplus in Kenya and suggest that R. microplus Clade A is prevalent in cattle in sub-Saharan Africa. These and other recent findings of widespread occurrence of R. microplus in Africa provide a strong justification for urgent surveillance to determine and monitor the spread of R. microplus and vector competence of Boophilus ticks for B. bovis in Africa, with the ultimate goal of strategic control.

A cysteine protease of Babesia microti and its interaction with tick cystatins

Babesiosis is a tick-borne protozoonosis caused by Babesia, which can cause fever, hemolytic anemia, hemoglobinuria, and even death. Babesia microti is a parasite found in rodents and can be pathogenic to humans. In this study, the full-length cDNA of a B. microti cysteine protease (BmCYP) was expressed and the recombinant rBmCYP protein analyzed and characterized. BmCYP is encoded by an ORF of 1.3 kb, with a predicted molecular weight of 50 kDa and a theoretical pI of 8.5. The amino acid sequence of BmCYP exhibits an identity of 32.9 to 35.2% with cysteine proteases of Babesia ovis, Babesia bovis, and Theileria, respectively. The results of the proteinase assays show that rBmCYP has cysteine protease enzymatic activity. In addition, we demonstrate that tick cystatins rRhcyst-1 and rRhcyst-2 were able to effectively inhibit the activity of rBmCYP; the inhibition rates were 57.2% and 30.9%, respectively. Tick cystatins Rhcyst-1 and Rhcyst-2 were differentially expressed in ticks that fed on Babesia-infected mice relative to non-infected control ticks. Our results suggest that BmCYP is a functional enzyme with cysteine protease enzymatic activity and may be involved in tick-B. microti interactions.

Catalogue of stage-specific transcripts in Ixodes ricinus and their potential functions during the tick life-cycle

Background

The castor bean tick Ixodes ricinus is an important vector of several clinically important diseases, whose prevalence increases with accelerating global climate changes. Characterization of a tick life-cycle is thus of great importance. However, researchers mainly focus on specific organs of fed life stages, while early development of this tick species is largely neglected.

Methods

In an attempt to better understand the life-cycle of this widespread arthropod parasite, we sequenced the transcriptomes of four life stages (egg, larva, nymph and adult female), including unfed and partially blood-fed individuals. To enable a more reliable identification of transcripts and their comparison in all five transcriptome libraries, we validated an improved-fit set of five I. ricinus-specific reference genes for internal standard normalization of our transcriptomes. Then, we mapped biological functions to transcripts identified in different life stages (clusters) to elucidate life stage-specific processes. Finally, we drew conclusions from the functional enrichment of these clusters specifically assigned to each transcriptome, also in the context of recently published transcriptomic studies in ticks.

Results

We found that reproduction-related transcripts are present in both fed nymphs and fed females, underlining the poorly documented importance of ovaries as moulting regulators in ticks. Additionally, we identified transposase transcripts in tick eggs suggesting elevated transposition during embryogenesis, co-activated with factors driving developmental regulation of gene expression. Our findings also highlight the importance of the regulation of energetic metabolism in tick eggs during embryonic development and glutamate metabolism in nymphs.

Conclusions

Our study presents novel insights into stage-specific transcriptomes of I. ricinus and extends the current knowledge of this medically important pathogen, especially in the early phases of its development.

Initiator and executioner caspases in salivary gland apoptosis of Rhipicephalus haemaphysaloides

BACKGROUND

Apoptosis is fundamental in maintaining cell balance in multicellular organisms, and caspases play a crucial role in apoptosis pathways. It is reported that apoptosis plays an important role in tick salivary gland degeneration. Several different caspases have been found in ticks, but the interactions between them are currently unknown. Here, we report three new caspases, isolated from the salivary glands of the tick Rhipicephalus haemaphysaloides.

METHODS

The full-length cDNA of the RhCaspases 7, 8 and 9 genes were obtained by transcriptome, and RhCaspases 7, 8 and 9 were expressed in E. coli; after protein purification and immunization in mice, specific polyclonal antibodies (PcAb) were created in response to the recombinant protein. Reverse-transcription quantitative PCR (RT-qPCR) and western blot were used to detect the existence of RhCaspases 7, 8 and 9 in ticks. TUNEL assays were used to determine the apoptosis level in salivary glands at different feeding times after gene silencing. The interaction between RhCaspases 7, 8 and 9 were identified by co-transfection assays.

RESULTS

The transcription of apoptosis-related genes in R. haemaphysaloides salivary glands increased significantly after tick engorgement. Three caspase-like molecules containing conserved caspase domains were identified and named RhCaspases 7, 8 and 9. RhCaspase8 and RhCaspase9 contain a long pro-domain at their N-terminals. An RT-qPCR assay demonstrated that the transcription of these three caspase genes increased significantly during the engorged periods of the tick developmental stages (engorged larval, nymph, and adult female ticks). Transcriptional levels of RhCaspases 7, 8 and 9 in salivary glands increased more significantly than other tissues post-engorgement. RhCaspase9-RNAi treatment significantly inhibited tick feeding. In contrast, knockdown of RhCaspase7 and RhCaspase8 had no influence on tick feeding. Compared to the control group, apoptosis levels were significantly reduced after interfering with RhCaspase 7, 8 and 9 expressions. Co-transfection assays showed RhCaspase7 was cleaved by RhCaspases 8 and 9, demonstrating that RhCaspases 8 and 9 are initiator caspases and RhCaspase7 is an executioner caspase.

CONCLUSIONS

To the best of our knowledge, this is the first study to identify initiator and executioner caspases in ticks, confirm the interaction among them, and associate caspase activation with tick salivary gland degeneration.

Poor Unstable Midgut Microbiome of Hard Ticks Contrasts With Abundant and Stable Monospecific Microbiome in Ovaries

Culture-independent metagenomic methodologies have enabled detection and identification of microorganisms in various biological systems and often revealed complex and unknown microbiomes. In many organisms, the microbiome outnumbers the host cells and greatly affects the host biology and fitness. Ticks are hematophagous ectoparasites with a wide host range. They vector a number of human and animal pathogens and also directly cause major economic losses in livestock. Although several reports on a tick midgut microbiota show a diverse bacterial community, in most cases the size of the bacterial population has not been determined. In this study, the microbiome was quantified in the midgut and ovaries of the ticks Ixodes ricinus and Rhipicephalus microplus before, during, and after blood feeding. Although the size of bacterial community in the midgut fluctuated with blood feeding, it was overall extremely low in comparison to that of other hematophagous arthropods. In addition, the tick ovarian microbiome of both tick species exceeded the midgut 16S rDNA copy numbers by several orders of magnitude. This indicates that the ratio of a tick midgut/ovary microbiome represents an exception to the general biology of other metazoans. In addition to the very low abundance, the tick midgut diversity in I. ricinus was variable and that is in contrast to that found in the tick ovary. The ovary of I. ricinus had a very low bacterial diversity and a very high and stable bacterial abundance with the dominant endosymbiont, Midichloria sp. The elucidation of this aspect of tick biology highlights a unique tissue-specific microbial-invertebrate host interaction.

Identification and validation of internal reference genes for real-time quantitative polymerase chain reaction-based studies in Hyalomma anatolicum ticks

Crimean-Congo hemorrhagic fever (CCHF) is an emerging tick-borne viral disease caused by the orthonairovirus CCHF virus (CCHFV). Ticks of the genus Hyalomma are the viral reservoir and they represent the main vector transmitting the virus to their hosts during blood feeding. However, how CCHFV replicates in its natural arthropod host cells and the nature of virus/host interactions are still largely unknown. With the aim of developing tools for use in this field, we identified and validated expression of four candidate endogenous control tick genes in a Hyalomma anatolicum-derived cell line. These genes will be useful for normalization of viral/cellular transcripts in infection/expression studies or as internal controls in molecular epidemiology surveys of pathogens transmitted by Hyalomma ticks.

Ovary Proteome Analysis Reveals RH36 Regulates Reproduction via Vitellin Uptake Mediated by HSP70 Protein in Hard Ticks

Ticks are blood-sucking vector arthropods, which play an important role in transmitting pathogens between humans and animals. RH36 is an immunomodulatory protein expressed in the salivary glands, but not other organs, of partially fed Rhipicephalus haemaphysaloides ticks, and it reaches its peak on the day of tick engorgement. RH36 gene silencing inhibited tick blood feeding and induced a significant decrease in tick oviposition, indicating that another function of immunosuppressor RH36 was regulating tick reproduction. Why did RH36 protein expressed uniquely in the salivary gland regulate tick reproduction? RH36 regulated positively the expression of vitellogenin in ovary, which indicated RH36 protein played an important role in the integration of nutrition and reproduction. According to proteomic analysis, heat shock protein 70 (HSP70) was significantly down-regulated in the immature ovary of post-engorged ticks. In addition, gene silencing of HSP70 not only inhibited tick blood-sucking and the expression of vitellogenin, but also increased tick death rate. These results suggested RH36 affected tick vitellogenin uptake and then regulated ovary cell maturation by modulating the expression of HSP70 protein, and finally controlled tick oviposition.

Expression analysis of glutathione S-transferases and ferritins during the embryogenesis of the tick Haemaphysalis longicornis

In the tick life cycle, embryogenesis is the only stage of development wherein no blood meal is required. Nevertheless, even in the absence of a blood meal, which is the source of nutrients as well as the ferrous iron and heme that could cause oxidative stress in ticks, malondialdehyde (MDA) has been reported to increase during this period. Additionally, the knockdown of some oxidative stress-related molecules such as ferritin has resulted in abnormal eggs and embryonic death. Here, we investigate the gene and protein expression profiles of the identified glutathione S-transferases (GSTs) and ferritins (Fers) of the tick H. longicornis during embryogenesis through quantitative reverse transcription polymerase chain reaction (RT-qPCR) and Western blotting, respectively. We also confirm the lipid peroxidation and ferrous iron concentration level using a thiobarbituric acid reactive substances (TBARS) assay. Finally, we attempt to correlate these findings with the events occurring by establishing a staging process in H. longicornis embryos. Lipid peroxidation increased during the course of embryogenesis, as does the amount of GST proteins. On the other hand, the GST genes have high expression at the 1st day post-oviposition, during the early stage of embryogenesis and at day 10 during the period wherein the germ band is observable. Fer gene expression also starts to increase at day 10 and peaks at day 15. In the ferritin proteins, only the secretory ferritin (Fer2) is detected and constitutively expressed during embryogenesis. Events occurring during embryogenesis, such as energy production and iron metabolism for cellular proliferation and differentiation cause oxidative stress in the embryo. To counteract oxidative stress, it is possible that the embryo may utilize oxidative stress-related molecules such as GSTs and Fer2, which could be either maternally or embryo-derived.

Quantitative Proteomics Identifies Metabolic Pathways Affected by Babesia Infection and Blood Feeding in the Sialoproteome of the Vector Rhipicephalus bursa

The negative impact of ticks and tick-borne diseases on animals and human health is driving research to discover novel targets affecting both vectors and pathogens. The salivary glands are involved in feeding and pathogen transmission, thus are considered as a compelling target to focus research. In this study, proteomics approach was used to characterize Rhipicephalusbursa sialoproteome in response to Babesiaovis infection and blood feeding. Two potential tick protective antigens were identified and its influence in tick biological parameters and pathogen infection was evaluated. Results demonstrate that the R. bursa sialoproteome is highly affected by feeding but infection is well tolerated by tick cells. The combination of both stimuli shifts the previous scenario and a more evident pathogen manipulation can be suggested. Knockdown of ub2n led to a significative increase of infection in tick salivary glands but a brusque decrease in the progeny, revealing its importance in the cellular response to pathogen infection, which is worth pursuing in future studies. Additionally, an impact in the recovery rate of adults (62%), the egg production efficiency (45.75%), and the hatching rate (88.57 %) was detected. Building knowledge on vector and/or pathogen interplay bridges the identification of protective antigens and the development of novel control strategies.

A novel type 1 cystatin involved in the regulation of Rhipicephalus microplus midgut cysteine proteases

Rhipicephalus microplus is a cattle ectoparasite found in tropical and subtropical regions around the world with great impact on livestock production. R. microplus can also harbor pathogens, such as Babesia sp. and Anaplasma sp. which further compromise cattle production. Blood meal acquisition and digestion are key steps for tick development. In ticks, digestion takes place inside midgut cells and is mediated by aspartic and cysteine peptidases and, therefore, regulated by their inhibitors. Cystatins are a family of cysteine peptidases inhibitors found in several organisms and have been associated in ticks with blood acquisition, blood digestion, modulation of host immune response and tick immunity. In this work, we characterized a novel R. microplus type 1 cystatin, named Rmcystatin-1b. The inhibitor transcripts were found to be highly expressed in the midgut of partially and fully engorged females and they appear to be modulated at different days post-detachment. Purified recombinant Rmcystatin-1b displayed inhibitory activity towards typical cysteine peptidases with high affinity. Moreover, rRmcystatin-1b was able to inhibit native R. microplus cysteine peptidases and RNAi-mediated knockdown of the cystatin transcripts resulted in increased proteolytic activity. Moreover, rRmcystatin-1b was able to interfere with B. bovis growth in vitro. Taken together our data strongly suggest that Rmcystatin-1b is a regulator of blood digestion in R. microplus midgut.

Expression dynamics of autophagy-related genes in the cattle tick Rhipicephalus microplus during embryonic development and under increasing larval starvation

Rhipicephalus microplus is a hematophagous ectoparasite that significantly affects parasitized cattle. As a one-host tick its entire life cycle consists of free-living and parasitic forms. Its extraordinary ability to survive during prolonged off-host periods has been related to the process of cytoplasmic degradation called autophagy. In order to deepen our understanding of this process during R. microplus non-parasitic stages, we determined the expression dynamics of a set of five autophagy-related genes (ATG genes) during embryonic development and over an increasing larval starvation period of 50 days. We found two apparent successive waves of ATG genes transcriptional activation, which paralleled key embryonic changes such as cellularization and organogenesis, as well as nutrient utilization. Moreover, during increasing larval starvation, ATG genes were up-regulated cyclically every 10-15 days. Taken together, our results suggest that autophagy is playing a major role in embryo development and energy metabolism during starvation in R. microplus.

TOR as a Regulatory Target in Rhipicephalus microplus Embryogenesis

Embryogenesis is a metabolically intensive process carried out under tightly controlled conditions. The insulin signaling pathway regulates glucose homeostasis and is essential for reproduction in metazoan model species. Three key targets are part of this signaling pathway: protein kinase B (PKB, or AKT), glycogen synthase kinase 3 (GSK-3), and target of rapamycin (TOR). While the role of AKT and GSK-3 has been investigated during tick embryonic development, the role of TOR remains unknown. In this study, TOR and two other downstream effectors, namely S6 kinase (S6K) and eukaryotic translation initiation factor 4E-binding protein 1 (4E-BP1), were investigated in in vitro studies using the tick embryonic cell line BME26. First, we show that exogenous insulin can stimulate TOR transcription. Second, TOR chemical inhibition led to a decrease in BME26 cell viability, loss of membrane integrity, and downregulation of S6K and 4E-BP1 transcription. Conversely, treating BME26 cells with chemical inhibitors of AKT or GSK-3 did not affect S6K and 4E-BP1 transcription, showing that TOR is specifically required to activate its downstream targets. To address the role of TOR in tick reproduction, in vivo studies were performed. Analysis of relative transcription during different stages of tick embryonic development showed different levels of transcription for TOR, and a maternal deposition of S6K and 4E-BP1 transcripts. Injection of TOR double-stranded RNA (dsRNA) into partially fed females led to a slight delay in oviposition, an atypical egg external morphology, decreased vitellin content in eggs, and decreased larval hatching. Taken together, our data show that the TOR signaling pathway is important for tick reproduction, that TOR acts as a regulatory target in Rhipicephalus microplus embryogenesis and represents a promising target for the development of compounds for tick control.

A combination of antibodies against Bm86 and Subolesin inhibits engorgement of Rhipicephalus australis (formerly Rhipicephalus microplus) larvae in vitro

Background

Rhipicephalus microplus is a hard tick species that has a high impact on cattle health and production in tropical and subtropical regions. Recently, ribosomal DNA and morphological analysis resulted in the reinstatement of R. australis as a separate species from R. microplus. Both feed on cattle and can transmit bovine pathogens such as Anaplasma and Babesia species. The current treatment with acaricides is becoming increasingly less effective due to the emergence of resistant tick strains. A promising alternative can be found in the form of anti-tick vaccines. The available commercial vaccines can be used to control tick infestation, but the lack of a knockdown effect (> 90% reduction in tick numbers as seen with effective acaricides) hampers its widespread use, hence higher efficacious vaccines are needed. Instead of searching for new protective antigens, we investigated the efficacy of vaccines that contain more than one (partially) protective antigen. For screening vaccine formulations, a previously developed in vitro feeding assay was used in which R. australis larvae are fed sera that were raised against the candidate vaccine antigens. In the present study, the efficacy of the Bm86 midgut antigen and the cytosolic Subolesin (SUB) antigen were evaluated in vitro.

Results

Antiserum against recombinant Bm86 (rBm86) partially inhibited larval engorgement, whereas antiserum against recombinant SUB (rSUB) did not have any effect on feeding of larvae. Importantly, when larvae were fed a combination of antiserum against rBm86 and rSUB, a synergistic effect on significantly reducing larval infestations was found. Immunohistochemical analysis revealed that the rBm86 antiserum reacted with gut epithelium of R. australis larvae, whereas the antiserum against rSUB stained salivary glands and rectal sac epithelium.

Conclusions

Combining anti-Bm86 and anti-subolesin antibodies synergistically reduced R. australis larval feeding in vitro. Rhipicephalus australis is a one host tick, meaning that the larvae develop to nymphs and subsequently adults on the same host. Hence, this protective effect could be even more pronounced when larvae are used for infestation of vaccinated cattle, as the antibodies could then affect all three developmental stages. This will be tested in future in vivo experiments.

Genetic diversity of Bm86 sequences in Rhipicephalus (Boophilus) microplus ticks from Mexico: analysis of haplotype distribution patterns

Background

Ticks are a problem for cattle production mainly in tropical and subtropical regions, because they generate great economic losses. Acaricides and vaccines have been used to try to keep tick populations under control. This has been proven difficult given the resistance to acaricides and vaccines observed in ticks. Resistance to protein rBm86-based vaccines has been associated with the genetic diversity of Bm86 among the ectoparasite’s populations. So far, neither genetic diversity, nor spatial distribution of circulating Bm86 haplotypes, have been studied within the Mexican territory. Here, we explored the genetic diversity of 125 Bm86 cDNA gene sequences from R. microplus from 10 endemic areas of Mexico by analyzing haplotype distribution patterns to help in understanding the population genetic structure of Mexican ticks.

Results

Our results showed an average nucleotide identity among the Mexican isolates of 98.3%, ranging from 91.1 to 100%. Divergence between the Mexican and Yeerongpilly (the Bm86 reference vaccine antigen) sequences ranged from 3.1 to 7.4%. Based on the geographic distribution of Bm86 haplotypes in Mexico, our results suggest gene flow occurrence within different regions of the Mexican territory, and even the USA.

Conclusions

The polymorphism of Bm86 found in the populations included in this study, could account for the poor efficacy of the current Bm86 antigen based commercial vaccine in many regions of Mexico. Our data may contribute towards designing new, highly-specific, Bm86 antigen vaccine candidates against R. microplus circulating in Mexico.

Folate pathway modulation in Rhipicephalus ticks in response to infection

Folate pathways components were demonstrated to be present in RNA-sequencing data obtained from uninfected and pathogen-infected Rhipicephalus ticks. Here, PCR and qPCR allowed the identification of folate-related genes in Rhipicephalus spp. ticks and in the tick cell line IDE8. Genes coding for GTP cyclohydrolase I (gch-I), thymidylate synthase (ts) and 6-pyrovoyltetrahydropterin (ptps) were identified. Differential gene expression was evaluated by qPCR between uninfected and infected samples of four biological systems, showing significant upregulation and largest fold-change for the gch-I gene in the majority of the biological systems, supporting the selection for functional analysis by RNAi silencing. Efficient knockdown of the gch-I gene in uninfected and Ehrlichia canis-infected IDE8 cells showed no detectable impact on the capacity of the bacteria to invade or replicate in the tick cells. Overall, this work demonstrated an increase in the expression of some folate-related genes, though not always statistically significantly, in the presence of infection, suggesting gene expression modulation of these pathways, either as a tick response to an invader or manipulation of the tick cell machinery by the pathogens to their advantage. This discovery points to folate pathways as interesting targets for further studies.

Transcriptome and Proteome Response of Rhipicephalus annulatus Tick Vector to Babesia bigemina Infection

A system biology approach was used to gain insight into tick biology and interactions between vector and pathogen. Rhipicephalus annulatus is one of the main vectors of Babesia bigemina which has a massive impact on animal health. It is vital to obtain more information about this relationship, to better understand tick and pathogen biology, pathogen transmission dynamics, and new potential control approaches. In ticks, salivary glands (SGs) play a key role during pathogen infection and transmission. RNA sequencing obtained from uninfected and B. bigemina infected SGs obtained from fed female ticks resulted in 6823 and 6475 unigenes, respectively. From these, 360 unigenes were found to be differentially expressed (p < 0.05). Reversed phase liquid chromatography-mass spectrometry identified a total of 3679 tick proteins. Among them 406 were differently represented in response to Babesia infection. The omics data obtained suggested that Babesia infection lead to a reduction in the levels of mRNA and proteins (n = 237 transcripts, n = 212 proteins) when compared to uninfected controls. Integrated transcriptomics and proteomics datasets suggested a key role for stress response and apoptosis pathways in response to infection. Thus, six genes coding for GP80, death-associated protein kinase (DAPK-1), bax inhibitor-1 related (BI-1), heat shock protein (HSP), heat shock transcription factor (PHSTF), and queuine trna-ribosyltransferase (QtRibosyl) were selected and RNA interference (RNAi) performed. Gene silencing was obtained for all genes except phstf. Knockdown of gp80, dapk-1, and bi-1 led to a significant increase in Babesia infection levels while hsp and QtRibosyl knockdown resulted in a non-significant decrease of infection levels when compared to the respective controls. Gene knockdown did not affect tick survival, but engorged female weight and egg production were affected in the gp80, dapk-1, and QtRibosyl-silenced groups in comparison to controls. These results advanced our understanding of tick-Babesia molecular interactions, and suggested new tick antigens as putative targets for vaccination to control tick infestations and pathogen infection/transmission.

Molecular characterization, tissue-specific expression, and RNA knockdown of the putative heat shock cognate 70 protein from Rhipicephalus haemaphysaloides

Heat shock cognate 70-kDa protein (RH-Hsc70) was identified from a cDNA library synthesized from the sialotranscriptomes of unfed and fed Rhipicephalus haemaphysaloides. The RH-Hsc70 open reading frame is 1950 bp long and encodes a protein that is 649 amino acids in length, with a predicted molecular weight of 71.1 kDa and a theoretical pI of 5.43. RH-Hsc70 exhibits 98% amino acid identity with Hsc70 in Haemaphysalis flava and 83% identity with Hsc70 in arthropods and mammals. RH-Hsc70 was mainly expressed in nymphs and adult ticks, not in larvae. Real-time quantitative PCR analysis indicated that RH-Hsc70 mRNA expression was induced by blood feeding in adult ticks. In addition, RH-Hsc70 gene expression was higher in the ovaries of fed adult ticks than that in the midguts, salivary glands, and fat bodies of unfed or fed adult ticks. RH-Hsc70 gene knockdown inhibited tick blood feeding, significantly decreased tick engorgement rate, and increased tick death rate. These data illustrate the importance of RH-Hsc70 in tick blood feeding and aging, which makes it a promising candidate for the development of anti-tick vaccines.