Cloning and expression pattern of a heat shock cognate protein 70 gene in ticks (Haemaphysalis flava)

Comparative analysis of the anticoagulant activities and immunogenicity of HSC70 and HSC70TKD of Haemaphysalis flava

BACKGROUND

Haemaphysalis flava is a hematophagous ectoparasite that acquires the nutrition needed for development and reproduction by sucking blood and digesting the blood meal. During blood-sucking and blood-meal digestion, the prevention of blood coagulation is important for this tick. Previous studies have shown that heat shock cognate 70 (HSC70) protein has certain anticoagulant activities, but its immunogenicity remains unclear. Also, whether the mutation of individual bases of the TKD-like peptide of HSC70 through the overlap extension method can change its anticoagulant activities and immunogenicity remains to be investigated.

METHODS

The gene encoding the HSC70 protein was cloned from a complementary DNA library synthesized from H. flava. The coding gene of the TKD-like peptide of HSC70 was mutated into a TKD peptide coding gene (HSC70^TKD) using the overlap extension method. Escherichia coli prokaryotic expression plasmids were constructed to obtain the recombinant proteins of HSC70 (rHSC70) and HSC70^TKD (rHSC70^TKD). The purified rHSC70 and rHSC70^TKD were evaluated at different concentrations for anticoagulant activities using four in vitro clotting assays. Emulsifying recombinant proteins with complete and incomplete Freund's adjuvants were subcutaneously immunized in Sprague Dawley rats. The serum antibody titers and serum concentrations of interferon-gamma (IFN-γ) and interleukin-4 (IL-4) were detected using an indirect enzyme-linked immunosorbent assay to assess the immunogenicity of rHSC70 and rHSC70^TKD.

RESULTS

The open reading frame of HSC70 was successfully amplified and found to have a length of 1958 bp. The gene encoding the TKD-like peptide of HSC70 was artificially mutated, with the 1373-position adenine (A) of the original sequence mutated into guanine (G), the 1385-position cytosine (C) mutated into G and the 1386-position G mutated into C. rHSC70 and rHSC70^TKD that fused with His-tag were obtained using the expression plasmids pET-28a-HSC70 and pET-28a-HSC70^TKD, respectively. rHSC70 and rHSC70^TKD prolonged the thrombin time (TT) and reduced the fibrinogen (FIB) content in the plasma, but did not affect the prothrombin time (PT) or activated partial thromboplastin time (APTT) when compared to the negative control. Interestingly, the ability of rHSC70^TKD to prolong the TT and reduce the FIB content in the plasma was better than that of rHSC70. The specific antibody titers of both rHSC70 and rHSC70^TKD in rat serum reached 1:124,000 14 days after the third immunization. The serum concentration of IFN-γ in the rHSC70^TKD group was higher than that in the rHSC70 group. The rHSC70 group has the highest serum concentration of IL-4, and the serum concentration of IL-4 in the rHSC70^TKD group was higher than that in the negative group.

CONCLUSIONS

rHSC70 and rHSC70^TKD exhibited anticoagulant activities by prolonging the TT and reducing the FIB content in vitro. rHSC70^TKD had better anticoagulant activities than rHSC70. Both rHSC70 and rHSC70^TKD had good immunogenicity and induced humoral and cellular immunity.

Protein profiling of hemolymph in Haemaphysalis flava ticks

Background

Tick hemolymph bathes internal organs, acts as an exchange medium for nutrients and cellular metabolites, and offers protection against pathogens. Hemolymph is abundant in proteins. However, there has been limited integrated protein analysis in tick hemolymph thus far. Moreover, there are difficulties in differentiating tick-derived proteins from the host source. The aim of this study was to profile the tick/host protein components in the hemolymph of Haemaphysalis flava.

Methods

Hemolymph from adult engorged H. flava females was collected by leg amputation from the Erinaceus europaeus host. Hemolymph proteins were extracted by a filter-aided sample preparation protocol, digested by trypsin, and assayed by liquid chromatography–tandem mass spectrometry (LC–MS/MS). MS raw data were searched against the UniProt Erinaceidae database and H. flava protein database for host- and tick-derived protein identification. Protein abundance was further quantified by intensity-based absolute quantification (iBAQ).

Results

Proteins extracted from hemolymph unevenly varied in size with intense bands between 100 and 130 kDa. In total, 312 proteins were identified in the present study. Therein 40 proteins were identified to be host-derived proteins, of which 18 were high-confidence proteins. Top 10 abundant host-derived proteins included hemoglobin subunit-α and subunit-β, albumin, serotransferrin-like, ubiquitin-like, haptoglobin, α-1-antitrypsin-like protein, histone H2B, apolipoprotein A-I, and C3-β. In contrast, 169 were high-confidence tick-derived proteins. These proteins were classified into six categories based on reported functions in ticks, i.e., enzymes, enzyme inhibitors, transporters, immune-related proteins, muscle proteins, and heat shock proteins. The abundance of Vg, microplusin and α-2-macroglobulin was the highest among tick-derived proteins as indicated by iBAQ.

Conclusions

Numerous tick- and host-derived proteins were identified in hemolymph. The protein profile of H. flava hemolymph revealed a sophisticated protein system in the physiological processes of anticoagulation, digestion of blood meal, and innate immunity. More investigations are needed to characterize tick-derived proteins in hemolymph.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05287-7.

Cloning, expression, and function of ferritins in the tick Haemaphysalis flava

The full-length cDNA of two ferritins of Haemaphysalis flava were cloned after which recombinant Hf-FER1 and Hf-FER2 were expressed and their function was analyzed. In addition, RNA interference (RNAi) based on the injection of Hf-fer1 or Hf-fer2 dsRNA into fully engorged female ticks was performed. The cDNA encoding Hf-FER1 is 834 bp in length. It contains an iron-responsive element in the 5' untranslated region and encodes 174 amino acid residues. The full-length cDNA of Hf-FER2 contains 696 bp and encodes 199 amino acids, including a putative signal peptide sequence. Hf-FER1 and Hf-FER2 both have the ferroxidase iron center and the ferrihydrite nucleation center. The evolutionary relationship of Hf-FER1 and Hf-FER2 was established, and the predicted quaternary structures were assembled as typical spherical shells composed of 24 subunits which was demonstrated by nature PAGE. Real-time PCR showed that Hf-fer1 and Hf-fer2 were expressed in all developmental stages, with the highest expression in fully engorged females. The expression of Hf-fer1 and Hf-fer2 were relatively high in unfed larvae. Hf-fer1 was expressed in all tissues and was especially abundant in the salivary glands of fully engorged females. In contrast, the highest levels of Hf-fer2 were found in the midgut of fully engorged females, and no expression was found in the salivary glands of this life stage. Both recombinant Hf-FER1 and Hf-FER2 had iron-binding capabilities. Silencing of both Hf-fer1 and Hf-fer2 affected fecundity. Compared to the control, the percentage of ticks that laid eggs in the Hf-fer1 and Hf-fer2 RNAi groups was 73.3% and 66.7%, respectively. The silenced ticks that laid eggs had lower egg weight to body weight ratios, and the eggs had abnormal morphologies. The hatchability of eggs with normal morphology in the Hf-fer1 and Hf-fer2 silenced groups was 47.8% and 22.8%, respectively, which was significantly different from the control group (P < 0.005). These findings indicate that Hf-FER1 and Hf-FER2 play important roles in the iron storage of H. flava.

The complex evolution of the metazoan HSP70 gene family

The metazoan 70-kDa heat shock protein (HSP70) family contains several members localized in different subcellular compartments. The cytosolic members have been classified into inducible HSP70s and constitutive heat shock cognates (HSC70s), but their distinction and evolutionary relationship remain unclear because of occasional reports of “constitutive HSP70s” and the lack of cross-phylum comparisons. Here we provide novel insights into the evolution of these important molecular chaperones. Phylogenetic analyses of 125 full-length HSP70s from a broad range of phyla revealed an ancient duplication that gave rise to two lineages from which all metazoan cytosolic HSP70s descend. One lineage (A) contains a relatively small number of genes from many invertebrate phyla, none of which have been shown to be constitutively expressed (i.e., either inducible or unknown). The other lineage (B) included both inducible and constitutive genes from diverse phyla. Species-specific duplications are present in both lineages, and Lineage B contains well-supported phylum-specific clades for Platyhelminthes, Rotifera, Nematoda, Porifera/Cnidaria, and Chordata. Some genes in Lineage B have likely independently acquired inducibility, which may explain the sporadic distribution of “HSP70” or “HSC70” in previous phylogenetic analyses. Consistent with the diversification history within each group, inducible members show lower purifying selection pressure compared to constitutive members. These results illustrate the evolutionary history of the HSP70 family, encouraging us to propose a new nomenclature: “HSP70 + subcellular localization + linage + copy number in the organism + inducible or constitutive, if known.” e.g., HSP70cA1i for cytosolic Lineage A, copy 1, inducible.

Cloning of four HSPA multigene family members in Haemaphysalis flava ticks

The heat shock protein 70 (HSPA) family and their genes have been studied in ticks and are considered as possible antigen candidates for the development of anti-tick vaccines. However, knowledge about their members, structure and function in ticks is incomplete. Based on our transcriptomic data, the full length of four HSPA genes in Haemaphysalis flava (Acari: Ixodidae) was cloned via rapid amplification of cDNA ends. The open reading frame of HSPA2A, HSPA2B, HSPA5 and HSPA9 was 1920, 1911, 1983 and 2088 bp in length, respectively. Three family signatures and one localization motif were in the encoding proteins. HSPA2A and HSPA2B were predicted to be located at cytoplasm/nucleus, whereas HSPA5 and HSPA9 were at endoplasmic reticulum and mitochondria, respectively. In silico simulation demonstrated that those proteins had distinct numbers of α-helixes, extended strands and coils, and different antigenic epitopes. Expression of HSPA5 and HSPA9 in the salivary gland was significantly higher in partially-engorged female adult ticks than the fully-engorged (P < 0.01) as shown by a quantitative polymerase chain reaction. Our data indicated that H. flava ticks had at least four HSPA genes encoding proteins with different cellular locations, structures and expression profiles, suggesting their diverse roles in tick biology.

A survey of proteins in midgut contents of the tick, Haemaphysalis flava, by proteome and transcriptome analysis

Tick blood meals are stored and digested in their midguts. Blood digestion is complex, and many proteins are involved. Study of the tick-derived proteins in the midgut content may aid in the discovery of active molecules that would be useful for anti-tick vaccines. We analyzed the midgut content proteomes of partially engorged female Haemaphysalis flava, fully engorged female H. flava, and hedgehog serum using liquid chromatography tandem-mass spectrometry and label-free quantitation. In this study, high-confidence protein profiling of tick midgut content was determined. Based on the search against our in-house transcriptome database, the 28 high-confidence proteins were identified. Of these, 17 were identified as tick-derived, and the rest were of unspecified origin (proteins that could not be differentiated as host-derived or tick-derived proteins). The function of these midgut content proteins identified here may involve nutrient transportation, anti-coagulation, erythrocyte lysis, detoxification, lipid metabolism, and immunization. The presence of hemoglobin suggested that the red blood cells were lysed in the gut lumen. The midgut contents contain a large amount of fibrinogen and it has the ability to clot immediately. The midgut contained mostly host-derived proteins, and these host proteins provide rich nutrients for tick development and reproduction. However, some intracellular proteins were also identified, suggesting the possibility of shedding of the midgut epithelium and ingestion of saliva during feeding. This finding advances our understanding of the digestive mechanism and will be useful in the screening of vaccine antigens.

De novo transcriptome sequencing and comparative analysis of Haemaphysalis flava Neumann, 1897 at larvae and nymph stages

Haemaphysalis flava Neumann, 1897 (H. flava) is of public health significance due to its capability of transmitting several pathogens such as Rickettsia, Ehrlichia, Bartonella and Francisella tularensis. However, lack of complete genome, transcriptome and proteome information limits the understanding of the biology of H. flava. Here, the total RNA of H. flava was collected separately at larvae and nymph stages and analyzed with high-throughput RNA sequencing technology. The obtained data were assembled and annotated based on the near origin species in the Nr database. The functions of the unigenes were annotated and classified by seven databases, including Nr, Nt, Pfam, KOG, Swiss-Prot, GO and KEGG. A total of 61,850,967 and 79,579,368 clean reads were obtained with a data bulk of 9.28 G and 11.94 G in larvae and nymph stages, respectively. The number of unigenes was 440,896, with 48.6% of them being matched to the Nr database and 51.4% remaining unknown. Additionally, 1,776,404 SNPs were identified in the unigenes. Differential analysis revealed 80 differentially expressed genes (DEGs), including 56 up-regulated genes and 24 down-regulated genes in the nymph versus larvae. qPCR confirmed 4 of the 56 up-regulated genes and 4 of the down-regulated genes. KEGG analysis of the DEGs showed that aldehyde dehydrogenase and sorbitol dehydrogenase, two up-regulated unigenes in nymph versus larvae, were both matched to the top three enriched pathways: "chloroalkane and chloroalkene degradation", "fatty acid degradation" and "glycolysis and gluconeogenesis". This is the first report on the whole transcriptome of H. flava at larvae and nymph stages. This study contributes to the understanding of H. flava at the gene expression level in different developmental stages and provides a theoretical basis for the development of vaccines against H. flava.

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.

Cloning and Expression Profile of Glyceraldehyde-3-Phosphate Dehydrogenase in Haemaphysalis flava (Acari: Ixodidae)

Haemaphysalis flava (Acari: Ixodidae) harbors pathogenic microorganisms and transfers these to hosts during blood feeding. Proteomic analysis in the midgut contents of H. flava detected glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and contig 1683 was retrieved as a GAPDH gene fragment by searching our previous transcriptomic library. In the study, the 5' and 3' ends of contig 1683 were cloned by rapid amplification of cDNA ends (RACE) and a full length, 1340 bp cDNA of Hf-GAPDH was obtained. The open-reading frame had 999 bp and coded for 333 amino acids. Hf-GAPDH was predicted to have an N-terminal NAD binding domain and a C-terminal glyceraldehyde dehydrogenase catalytic domain. The molecular structure of Hf-GAPDH was analyzed and the evolutionary relationship also established. The GAPDH protein sequence was conserved among ticks. The expression pattern of Hf-GAPDH, analyzed by real-time PCR, significantly differed among life phases, feeding stages, and tissues. As the ticks grew, the expression level of Hf-GAPDH was up-regulated. The expression levels of Hf-GAPDH in salivary glands and midguts from half-engorged ticks were lower than the same tissues from engorged ticks. This study will provide reference data for the follow-up verification of the GAPDH-related function and the feasibility as a potential anti-tick vaccine.

HSC70 from Haemaphysalis flava (Acari: Ixodidae) exerts anticoagulation activity in vitro

Ticks and tick-borne diseases are major global health threats. During blood feeding, ticks insert their hypostomes into hosts and inject an array of anticoagulant molecules to maintain fluidity of the blood-meal. These anticoagulant molecules may provide insights into understanding the feeding biology of ticks and to develop vaccines against infestations. In Haemaphysalis flava, the heat shock cognate 70 (HSC70), a member of the heat shock protein (HSP) family, is differentially expressed in salivary glands at different levels of engorgement during blood feeding. However, its function in ticks is largely not known. The present study was designed to explore the possible effects of HSC70 on the plasma. The open reading frame (ORF) of HSC70 was expressed in a prokaryotic system, and recombinant HSC70 (rHSC70) was purified and characterized. The anticoagulation activity of rHSC70 was estimated by measuring prothrombin time (PT), activated partial thromboplastin time (APTT), thrombin time (TT) and fibrinogen (FIB) with/without its inhibitor, VER155008. The results demonstrated that rHSC70 from H. flava extended TT (P < 0.001) and FIB clotting times (>300 s), but showed little effect on PT and APTT. Adding an inhibitor reversed anticlotting effects of rHSC70 on TT and FIB. These data indicate that rHSC70 is an anticoagulant agent, and the anticlotting activity likely attributes to the inhibition of thrombin and the transformation of fibrinogen into fibrin.

Proteomic profiling of the midgut contents of Haemaphysalis flava

Scant information is available regarding the proteins involved in blood meal processing in ticks. Here, we aimed to highlight the midgut proteins involved in preventing blood meal coagulation, and in facilitating intracellular digestion in the tick Haemaphysalis flava. Proteins were extracted from the midgut contents of fully engorged and partially engorged ticks. We used liquid chromatography tandem-mass spectrometry (LC-MS/MS) analysis to identify 131 unique peptides, and 102 proteins. Of these, 15 proteins, each with at least two unique peptides, were recognized with high confidence. We also retrieved 18 unigenes from our previous published transcriptomic libraries of the midguts and salivary glands of H. flava, and inferred the primary structures of nine proteins and fragments of five proteins. There were 23 and 21 unique proteins in the midgut contents of fully engorged and partially engorged ticks, respectively. We detected 58 shared proteins in the midgut contents of both fully engorged and partially engorged ticks. Of these, seven were significantly differentially expressed between fully engorged and partially engorged ticks: actin, calmodulin, elongation factor-1α, hsp90, multifunctional chaperone, tubulin α, and tubulin β. Our results demonstrated that the proteome of the midgut contents, combined with the transcriptome of the midgut, was a viable method for the reinforcement of protein identification. This method will facilitate further study of blood meal processing by ticks, as well as the identification of clues for tick infestation control. The existence of numerous proteins detected in the midgut contents also highlight the complexity of blood digestion in ticks; this area is in need of further investigation.