Salivary gland proteome analysis of developing adult female Haemaphysalis longicornis ticks: molecular motor and TCA cycle-related proteins play an important role throughout development
Parasites & Vectors volume 12, Article number: 613 (2019)
Ticks are notorious blood-feeding arthropods that can spread a variety of deadly diseases. The salivary gland is an important organ for ticks to feed on blood, and this organ begins to develop rapidly when ixodid ticks suck blood. When these ticks reach a critical weight, the salivary glands stop developing and begin to degenerate. The expression levels of a large number of proteins during the development and degeneration of salivary glands change, which regulate the biological functions of the salivary glands. Furthermore, to the best of our knowledge, there are only a few reports on the role of molecular motor and TCA cycle-related proteins in the salivary glands of ticks.
We used iTRAQ quantitative proteomics to study the dynamic changes in salivary gland proteins in female Haemaphysalis longicornis at four feeding stages: unfed, partially fed, semi-engorged and engorged. Using bioinformatics methods to analyze the dynamic changes of a large number of proteins, we found that molecular motor and TCA cycle-related proteins play an important role in the physiological changes of the salivary glands. The results of RNAi experiments showed that when dynein, kinesin, isocitrate dehydrogenase and citrate synthase were knocked down independently, the weight of the engorged female ticks decreased by 63.5%, 54.9%, 42.6% and 48.6%, respectively, and oviposition amounts decreased by 83.1%, 76.0%, 50.8%, and 55.9%, respectively, and the size of type III acini of females salivary glands decreased by 35.6%, 33.3%, 28.9%, and 20.0%, respectively.
The results showed that the expression of different types of proteins change in different characteristics in salivary glands during the unfed to engorged process of female ticks. Corresponding expression changes of these proteins at different developmental stages of female ticks are very important to ensure the orderly development of the organ. By analyzing these changes, some proteins, such as molecular motor and TCA cycle-related proteins, were screened and RNAi carried out. When these mRNAs were knocked down, the female ticks cannot develop normally. The research results provide a new protein target for the control of ticks and tick-borne diseases.
Ticks are harmful external parasitic arthropods that transmit a variety of pathogenic microorganisms to their hosts while sucking blood. These pathogens can lead to the development of many lethal diseases in hosts, such as Lyme disease, human anaplasmosis, babesiosis, and tick-borne encephalitis . To adapt to a blood-feeding lifestyle, ticks have evolved special body structures and blood-feeding organs with unique functions. Among them, the salivary glands are important for successful blood-feeding, because they secrete saliva that includes large amounts of powerful pharmacologically active proteins, which can inhibit host blood coagulation [2, 3], immune response  and vascular repair . Salivary glands are also the main organs involved in tick-borne diseases, because saliva contains a large number of pathogens [6, 7].
The salivary glands of different female ixodidae tick species exhibit similar developmental characteristics. The salivary glands of female ticks are composed of different types of acini (I, II, III) with ducts connected to each other . Salivary glands have the ability to regulate water balance during the unfed stage. When a tick finds a host and pre-attachment occurs, some genes are upregulated in the salivary glands [9,10,11]. The development and functional adjustment of the salivary glands are more related to blood-feeding. When the tick starts to suck blood, the salivary glands begin to develop rapidly, and expressions of a large number of genes are upregulated rapidly. However, the salivary glands begin to degenerate when the female tick reaches a critical weight [12,13,14,15]. Salivary gland degeneration is caused by the ecdysteroid hormone , which is also regulated by many proteins .
The coordination of various proteins plays a decisive role in the development and degeneration of salivary glands. Analyzing global changes in salivary gland proteins is critical to develop more effective molecular control methods. Quantitative proteomics was used to rapidly monitor global proteins in the development of salivary glands. Some studies on the proteome of tick salivary glands of Hyalomma and Amblyomma genus have been reported [18,19,20], but there are no studies on global proteins in the feeding cycle of female adult of tick salivary glands. In the present study, we focused on the expression changes of global proteins in salivary glands during the feeding cycle of female adult Haemaphysalis longicornis, providing deeper insights into the molecular mechanisms of tick salivary gland development and degeneration. Haemaphysalis longicornis is widely distributed in East Asia and Oceania  and is very harmful to animal husbandry and human health. This species is a major vector of multiple pathogens, such as Babesia, Rickettsia, Theileria, and severe fever with thrombocytopenia syndrome virus (SFTSV) [22, 23].
The widely used isobaric tags for relative and absolute quantification (iTRAQ) quantitative proteomics method has many advantages and can be used for large-scale and accurate quantitative studies of multicomponent samples. In the present study, we determined the protein expression kinetics of the salivary glands of adult H. longicornis females during the four feeding stages: unfed, partially fed, semi-engorged and engorged. In total, we identified 5831 salivary gland proteins in the four feeding stages. The expression of many proteins related to energy production and material transport was significantly increased in the salivary glands of female ticks after they started to suck blood. The increased expression of these proteins provides sufficient energy and materials for the rapid development of the salivary glands. These proteins were then individually knocked down with RNAi.
The results showed that after knockdown of the proteins dynein, kinesin, isocitrate dehydrogenase, and citrate synthase, female ticks were unable to suck blood normally. In particular, 52.5% and 57.5% of ticks died after silencing of dynein and kinesin. There are few reports on these four proteins regulating the function of tick salivary glands. Dynein and kinesin, which are molecular motor proteins, play an important role in mitosis, meiosis and transport of cellular cargo [24,25,26]. In addition, isocitrate dehydrogenase and citrate synthase are indispensable proteins of the tricarboxylic acid (TCA) cycle . These proteins provide sufficient energy during blood-feeding to enable rapid expansion of the tick body.
Tick breeding and protein extraction
Haemaphysalis longicornis were captured from Xiaowutai Mountain National Nature Reserve of China. Ticks were allowed to feed on the ears of New Zealand white rabbits by using earmuffs made of white cloth. Rabbits were kept in constant temperature culture room (25 ± 1 °C). When not sucking blood, the ticks were cultured in artificial climate incubators at 25 ± 1 °C and 75% relative humidity. Female adult ticks in four feeding stages, unfed (approximately 1.8 ± 0.1 mg, 0.32 ± 0.02 cm), partially fed (approximately 2.5 days post-feeding, 11.9 ± 0.4 mg, 0.48 ± 0.02 cm), mated semi-engorged (approximately 4 days post-feeding, 40.1 ± 0.8 mg, 0.67 ± 0.02 cm) and engorged (approximately 5.5 days post-feeding, 249.2 ± 6.8 mg, 0.99 ± 0.02 cm), were collected for this study. Haemaphysalis longicornis is a rather small tick, therefore, in order to obtain enough proteins for this experiment, approximately 500 adult female ticks were needed per group for experiments (approximately 6000 ticks in total). Approximately 50 New Zealand white rabbits (Oryctolagus cuniculus) were used to feed the ticks and the salivary glands were dissected from the ticks at each feeding stage. After being washed with PBS (0.01 M), the salivary glands were deposited immediately in PBS buffer containing cOmplete™ protease inhibitor cocktail (Roche, Mannheim, Germany) and quickly frozen at − 80 °C. When enough salivary glands were collected, protein extraction was performed as previously described . The salivary glands were fully ground in a glass homogenizer then transferred into a 50 ml centrifuge tube and centrifuged for 20 min (4 °C, 12,000×g). Then, supernatant was collected and an equal volume of Tris-saturated phenol (pH 7.8) was added to the supernatant. The sample was mixed thoroughly and centrifuged for 20 min (4 °C, 12,000×g). Then, an equal volume of 50 mM Tris-HCl (pH 8.0) was added and the solution mixed thoroughly, followed by centrifugation for 20 min (4 °C, 12,000×g). Then 0.1 M ammonium acetate in methanol was added to the solution after the supernatant was removed and stored at − 20 °C overnight. Then the solution was centrifuged for 20 min (4 °C, 12,000×g), and the supernatant removed. The protein pellets were washed again with methanol and were lyophilized and stored at −80 °C.
Protein digestion and iTRAQ labeling
Protein digestion was performed as previously described . Salivary gland protein samples (200 μg) from each blood-feeding stage were reduced with 10 mM dithiothreitol and then alkylated with 20 mM iodoacetamide. Then, the samples were processed using a filter-aided sample preparation (FASP) protocol . Enzymatic digestion was performed with sequencing-grade modified trypsin (1:2 w/w, Promega, Madison, USA) at 37 °C (strictly controlled) for 12 h. After digestion, the peptides were eluted (centrifugation at 12,000×g for 20 min) from the ultrafiltration membrane with iTRAQ dissolution buffer (AB SCIEX, Redwood, USA). The concentrations of the trypsin-digested peptides were detected with a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, USA) and an LC-MS system (Thermo Fisher Scientific) so that the components could be adjusted to equal concentrations. Enzyme efficiency was also monitored by LC-MS. Figure 1 shows the experimental workflow using iTRAQ. Each sample was labeled with iTRAQ (AB SCIEX) 4-plex reagents (114–117) according to the manufacturer’s instructions, respectively. After deionized water was added to stop the reaction, the four labeled samples were mixed together for further analysis.
High-pH reversed-phase (RP) fractionation
The iTRAQ-labeled peptide mix samples were separated by high-pH (pH = 10) C18 (Agela; 5 μm particle size, 100 Å pore size, 0.46 × 25 cm) reversed-phase high-performance liquid chromatography (C18 RP-HPLC). The peptides were eluted at a gradient flow of 1 ml/min, and the concentration of solvent B (solvent A: H2O containing 5 mM ammonium formate, pH 10.0; solvent B: acetonitrile (ACN) containing 5 mM ammonium formate, pH 10.0) in the elution solvent was raised to 60% (v/v) over 70 min. A tube of elution components was collected every minute and randomly combined. As a result, 30 eluted fractions were ultimately obtained for further LC-MS analysis.
Peptide mixtures in each eluted fraction were analyzed with a Q Exactive HF (Thermo Fisher Scientific) mass spectrometer coupled online to a nanoACQUITY UPLC M-Class system (waters, Milford, USA) as previously described . Each sample was desalted with a C18 RP trap column (5 μm particle size, 100 Å pore size, 180 μm ID × 20 mm length; Waters, Milford, USA) and then separated with a C18 RP analytical column (1.8 μm particle size, 100 μm ID × 150 mm length; Waters) at a flow rate of 300 nl/min using a linear gradient (0–40%) of solvent B (solvent A: 99.9% H2O + 0.1% formic acid; solvent B: 99.9% ACN + 0.1% formic acid) over 75 min. LC-MS data were acquired in the data-dependent acquisition mode. All mass spectrometry parameters were set as previously described . Three biological replicates were performed for iTRAQ analysis.
Protein identification and iTRAQ quantification
Thirty files (raw data) identified by mass spectrometry (MS) were combined into one MS/MS dataset and then analyzed using the SEQUEST algorithm embedded in Proteome Discoverer (version 2.1) software (Thermo Fisher Scientific). The H. longicornis proteins database derived from transcriptome sequencing (GenBank: GHLT00000000). Ticks were raised in laboratory conditions since they were collected from Xiaowutai. The ticks used in this study have been bred in the laboratory for more than 3 generations and fed on the ears of New Zealand white rabbits. In order to eliminate any contamination, host O. cuniculus and human sequences were used as a contaminated database for proteomic searching. Search parameters were applied as previously described : (i) a precursor mass tolerance of 10 ppm was used; (ii) a fragment mass tolerance of 0.5 Da was used; (iii) trypsin was set as the enzyme, and 2 missed cleavages were allowed; (iv) methionine oxidation was set as the variable modification; (v) the iTRAQ 4-plex reagents (N terminus and lysine residues, 144.102 Da) were defined, and carbamidomethylation of cysteine residues was set as the fixed modification; (vi) high-energy collision-induced dissociation (HCD) was chosen as the activation type; (vii) the “total peptide amount” tab was checked to normalize the protein levels by the protein ratio median; and (viii) a decoy database search was simultaneously performed to estimate the false discovery rate (FDR). The q-value was used as the criterion to evaluate FDR, and the target FDR was set to 0.01. Only quantified proteins with at least two unique peptides and high confidence (FDR < 0.01) were considered for further analysis. Protein quantification was based on the normalized spectrum abundance factor (NSAF). Proteins with expression changes greater than 1.5-fold were considered to be upregulated or downregulated.
To further reveal the functions of the proteins, a number of bioinformatics analyses were carried out. The GProX platform was used to cluster the salivary gland proteins with similar expression patterns in the 4 tick feeding stages . The number of clusters was set to 4, and a fixed regulation threshold (upper limit of 0.58 and lower limit of − 0.58, corresponding to the original ratios of approximately 1.5 and 0.67) was used [31, 32]. The minimal membership for the plot was set as 0.5. Other parameters were set to default values. PANTHER classification (http://www.pantherdb.org/) was used to carry out a Gene Ontology (GO) functional annotation. The pathways of the differentially expressed proteins (P-values < 0.05) were carried out using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.kegg.jp/kegg/) and KOBAS 3 software (http://kobas.cbi.pku.edu.cn/).
T7 RiboMAXTM Express RNAi System (Promega) was used to synthesize double-stranded RNA (dsRNA) according to the manufacturer’s instructions. Specific H. longicornis nucleotide sequence targeting positions 13456–13965, 104–515, 235–763, and 649–1291 of the dynein heavy chain (GenBank: GHLT01015731), the kinesin heavy chain (GenBank: GHLT01011136), isocitrate dehydrogenase (GenBank: GHLT01012103), and citrate synthase (GenBank: GHLT01017326), respectively, were selected for dsRNA synthesis. After the target mRNA sequences were cloned and sequenced, the correct target cDNA sequences for RNAi were used for dsRNA synthesis. The primers used to synthesize dsRNA included the T7 promoter sequence (underlined) were as follows: 5′-TAA TAC GAC TCA CTA TAG GCA GTG GGT GAC GGA TTT CTC-3′ and 5′-TAA TAC GAC TCA CTA TAG GGG AAC ACA GCA GAG CGA CAC-3′ (dynein heavy chain), 5′-TAA TAC GAC TCA CTA TAG GGC AAG TTC ATC GTC AAG TTC-3′ and 5′-TAA TAC GAC TCA CTA TAG GAC TCT GTT CTT GTC CTC GTG-3′ (kinesin heavy chain), 5′-TAA TAC GAC TCA CTA TAG GAT CCC GCA GAA GGC TAT T-3′ and 5′-TAA TAC GAC TCA CTA TAG GCG TAC AAG TTG GGC ATC AC-3′ (isocitrate dehydrogenase), 5′-TAA TAC GAC TCA CTA TAG GGC CAC CAT CTA CCG CAA CC-3′ and 5′-TAA TAC GAC TCA CTA TAG GGG GAC ACT CCA AAC AGC ACC-3′ (citrate synthase). Green fluorescent protein (GFP, GenBank: KX247384.1) dsRNA was synthesized as the control . Twenty unfed female H. longicornis were injected with 0.5–1 µl (4 µg/µl) dsRNA using a Hamilton syringe (33-gauge needle) in the lower right quadrant of the body . After injection, the ticks were cultured in a humid incubator for 24 h to recover, then were allowed to feed on the ears of New Zealand white rabbits using cloth earmuffs. The degree of mRNA knockdown was confirmed by analysis of mRNA expression levels using RT-qPCR. A digital microscope (DVM6; Leica, Oskar, Germany) was used to observe morphological changes of the salivary gland acini. Ten RNAi ticks at each feeding stage were used to observe the shape and size of salivary glands. Other physiological changes of female ticks after RNAi were measured: (i) tick mortality rate; (ii) time required for engorgement (starting with biting the host) and the weight of the engorged tick; (iii) number of eggs; and (vi) egg hatching rate.
RT-qPCR quantitative analysis
In order to observe the expression of transcription levels of these four proteins during the four feeding stages: unfed, partially fed, semi-engorged, and engorged. We performed RT-qPCR quantitative analysis of mRNA from different periods. In addition, RT-qPCR analysis was performed at 48 h after the unfed ticks were injected with dsRNA. Total RNA of salivary glands of 10 ticks were extracted using TRIzol Reagent (Invitrogen, Carlsbad, USA). After treatment with DNase, the cDNA was synthesized using the PrimeScript™ II 1st Strand cDNA Synthesis Kit according to the manufacturer’s protocol (Takara, Tokyo, Japan). The primers used for RT-qPCR were as follows: 5′-CTC CGG GTT GAA CAG GC-3′ and 5′-CAG TGG GTG ACG GAT TTC T-3’ (dynein heavy chain), 5′-GTG CTC ATG GGC TAC AA-3′ and 5′-GGG TAT GAT GCC TTG ACT-3′ (kinesin heavy chain), 5′-GCA TCC AGG AAC ATC GC-3′ and 5′-GCC GTC GGA CAT TCT CAT-3′ (isocitrate dehydrogenase), 5′-CAC CAT CTA CCG CAA CCT G-3′ and 5′-GGT CGT CGT ATC CGA GCA T-3′ (citrate synthase). Actin was used as an internal reference gene. The 2−ΔΔCq equation  was used to calculate the gene expression silencing rate.
The coefficients of variation (CV) were used to calculate the degree of data dispersion between repeated experiments (CV = SD/Mean × 100%), proteins abundance ratios with CV < 20% were selected for further analysis. R2 was used to analyze correlations between repeated experiments. If > 0.9, it means that the data we obtained is reproducible and accurate. The student’s t-test (Statistica 6.0) was used to analyze the quantitative results of RT-qPCR. If P-values < 0.05, it was considered that there was a significant difference in mRNA expression levels between different groups.
Protein identification and quantitation
Figure 1 shows the experimental workflow. After trypsin digestion, mass spectrometry identification and data analysis, the proteins were identified, and quantitative results were obtained (Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3, Additional file 4: Table S4). The first, second and third biological replicates identified 5059 (experiment 1), 5526 (experiment 2) and 5584 (experiment 3) proteins, respectively. Only quantified proteins with at least 2 unique peptides and high confidence (FDR < 1%) were considered for further analysis. The Venn diagram shows the overlap of proteins from three biological replicates, a total of 3667 high-confidence unique proteins were in all three experiments (Fig. 2a); 2507 high-confidence proteins were shared among all three experiments (Fig. 2a), accounting for 68.4% of the total proteins. The coefficients of variation (CV) for the partially fed:unfed protein abundance ratio (115:114), the mated semi-engorged:unfed ratio (116:114), and the engorged:unfed ratio (117:114) were calculated. Proteins whose abundance ratios with CV < 20% were selected for further analysis (Additional file 4: Table S4). log2 correlation coefficients for protein abundance ratios from the different experimental replicates were calculated. The results showed that the relative changes in protein abundance were correlated (R2 > 0.9) between different experiment replicates (Additional file 5: Figure S1).
Cluster analysis of differentially expressed proteins
After eliminating the null values in the quantitative results, we selected proteins with CV < 20% from among the 2507 identified proteins shared among all three experiments for further analysis. The construction of clusters was based on all the quantified proteins from the 115/114, 116/114 and 117/114 comparisons (Fig. 2b). We divided these proteins into 5 clusters according to the trends in the protein expression changes in the 4 different feeding stages (Fig. 2c, Additional file 6: Table S5). Proteins with values higher than 1.5 or lower than 0.67 were considered to be upregulated or downregulated, which corresponded to log 2 values of 0.58 and − 0.58, respectively.
Proteins with similar expression patterns were grouped together. Overall, 435 proteins exhibited no significant changes in expression (Cluster 0). The expression patterns of the upregulated and downregulated proteins were generally divided into five clusters by GProX clustering analysis: Cluster 1 contained 364 proteins; Cluster 2 contained 174 proteins; Cluster 3 contained 176 proteins; Cluster 4 contained 141 proteins; and Cluster 5 contained 229 proteins. There were many different types of proteins in each cluster; for example, the overall trend of protein expression in Cluster 1 from the unfed stage to engorgement was upregulation. Then, the protein expression remained largely unchanged: dynein heavy chain (115/114 = 1.8; 116/114 = 1.8; 117/114 = 1.9); kinesin heavy chain (115/114 = 1.6; 116/114 = 1.6; 117/114 = 1.4); isocitrate dehydrogenase (115/114 = 1.9; 116/114 = 1.8; 117/114 = 1.9); citrate synthase (115/114 = 1.7; 116/114 = 1.5; 117/114 = 1.8); angiotensin-converting enzyme (115/114 = 2.3; 116/114 = 2.5; 117/114 = 3.3). The overall trend of protein expression in Cluster 2 from the unfed stage to engorgement was continuous upregulation. The proteins in Cluster 2 were mainly enzymes, such as enolase (115/114 = 2.1; 116/114 = 2.0; 117/114 = 3.2). In addition, Cluster 2 also contained certain structural proteins, such as alpha tubulin (115/114 = 3.0; 116/114 = 3.2; 117/114 = 6.0). The overall trend of protein expression in Cluster 4 from the unfed stage to engorgement was continuous downregulation. The proteins in Cluster 4 included some phosphatases, such as deoxyuridine 5′-triphosphate nucleotidohydrolase (115/114 = 0.9; 116/114 = 0.8; 117/114 = 0.5) and phosphatidylinositol-4-phosphate 5-kinase (115/114 = 0.8; 116/114 = 0.8; 117/114 = 0.6). In addition, Cluster 4 also contained certain binding proteins, such as polyadenylate-binding protein-interacting protein (115/114 = 0.9; 116/114 = 0.7; 117/114 = 0.6).
GO function annotation
As shown in Fig. 3, the GO functions were annotated for all the differentially expressed proteins in the following comparisons: partially fed:unfed (115:114); mated semi-engorged:partially fed (116:115); and engorged:mated semi-engorged (117:116). The proteins were all grouped into three major functional groups: biological processes; cellular components; and molecular functions. Further enrichment analysis was conducted in these three categories (Fig. 3), the results showed that for the proteins upregulated and downregulated in partially fed ticks compared with unfed ticks, a total of 11 terms (means different kinds of biological functions) and 9 terms were enriched in the biological process category, respectively (Fig. 3a), 7 terms and 6 terms were enriched in the cellular component category, respectively (Fig. 3b), and 8 terms and 7 terms were enriched in the molecular function category, respectively (Fig. 3c). Upregulated and downregulated proteins in mated semi-engorged ticks compared with partially fed ticks show that a total of 6 terms and 6 terms were enriched in the biological process category, respectively (Fig. 3d), 3 terms and 3 terms were enriched in the cellular component category, respectively (Fig. 3e), and 2 terms and 3 terms were enriched in the molecular function category, respectively (Fig. 3f). Upregulated and downregulated proteins in engorged ticks compared with mated semi-engorged ticks show that a total of 11 terms and 10 terms were enriched in the biological process category, respectively (Fig. 3g), 7 terms and 4 terms were enriched in the cellular component category, respectively (Fig. 3h), and 7 terms and 5 terms were enriched in the molecular function category, respectively (Fig. 3i). GO annotation was also performed for the differentially expressed proteins in Clusters 1–5 (Additional file 7: Figure S2).
KEGG pathway analysis
The differentially expressed proteins in the 5 different Clusters were subjected to KEGG pathway enrichment analysis. In Clusters 1–4, and 5, 66, 30, 22, 15, and 35 terms were mapped to KEGG pathways, respectively (Additional file 8: Table S6, Additional file 9: Table S7, Additional file 10: Table S8, Additional file 11: Table S9, Additional file 12: Table S10). Terms with P-values < 0.05 were used to draw column diagrams. After screening, the number of terms was reduced to 34, 10, 3, 5, and 15 for Clusters 1–5, respectively (Fig. 4). The “RNA transport” term was distributed in all four Clusters except Cluster 2; the term “endocytosis” was distributed in Clusters 1, 2, and 4; the term “metabolic pathways” was distributed in Clusters 1, 2, and 5; and the terms “ribosome” and “protein processing in endoplasmic reticulum” were distributed in Clusters 1, 4 and 5.
Phenotypes associated with RNAi
After bioinformatics analysis, we found that the protein expression levels of dynein, kinesin, isocitrate dehydrogenase, and citrate synthase were upregulated. In addition, we analyzed the mRNA of these four proteins at the transcriptional level during the four feeding stages of salivary gland development by RT-qPCR (Fig. 5, Additional file 13: Table S11). It was found that the expression of transcription level and protein level was slightly different, but the expression change trend was the same. These proteins may directly affect the physiological function of salivary glands, so we further analysis their function using RNAi. We found that ticks could not suck blood normally after RNAi of dynein, kinesin, isocitrate dehydrogenase, and citrate synthase, and the body length of females represented 45.7 ± 6.1%, 51.2 ± 7.3%, 55.5 ± 8.0% and 53.5 ± 8.6% (mean ± SD) of the body length of females in the control group after blood-feeding for 4 days, respectively. There were significant differences in body length of females after RNAi of dynein, kinesin, isocitrate dehydrogenase and citrate synthase mRNA (Studentʼs t-test: t(2) = 10.34, P = 0.0092; t(2) = 7.45, P = 0.0175; t(2) = 5.16, P = 0.0354; and t(2) = 6.23, P = 0.0247, respectively) (Fig. 6). Female ticks had a higher mortality rate after RNAi of dynein (47.5%) and kinesin (42.5%). In addition, the development of the salivary glands was greatly affected after RNAi injection. Figure 7 shows the morphological changes of tick salivary glands following injection with different dsRNA. The development of the salivary gland acini in ticks injected with RNA targeting the four genes clearly did not progress to the same extent as that in ticks injected with control GFP dsRNA. The size of acini in the RNAi groups was significantly smaller than that in the control group. After dynein, kinesin, isocitrate dehydrogenase, and citrate synthase mRNA were interfered, the size of type III acini of females were about 60.7 ± 2.7%, 66.7 ± 3.5%, 71.1 ± 3.6% and 80.0 ± 2.1% (mean ± SD) of the size observed in the control group after blood-feeding for 5 days, respectively. There were statistically significant differences in the size of type III acini after RNAi of dynein, kinesin, isocitrate dehydrogenase and citrate synthase mRNA (Studentʼs t-test: t(2) = 34.00, P = 0.0008; t(2) = 26.11, P = 0.0014; t(2) = 20.19, P = 0.0024; and t(2) = 8.39, P = 0.0139, respectively). In addition, in the RNAi groups, the weights of the engorged female ticks and the egg hatching rates decreased compared with those of the control group (Table 1). Table 2 shows the qPCR results confirming that RNAi knocked down the mRNA of these four genes. The gene expression silencing rate was compared between female ticks RNAi with test genes dsRNA and GFP control dsRNA. The Cq values of the qPCR are provided in Additional file 14: Table S12 and Additional file 15: Table S13. After the injection of dsRNA from the target mRNA, the silencing rates of dynein, kinesin, isocitrate dehydrogenase, and citrate synthase in the salivary glands of female ticks reached 78.3%, 83.3%, 73.2% and 83.1%, respectively. Meanwhile, the silencing rates in the whole ticks were 81.2%, 85.6%, 69.1% and 86.7%, respectively.
In the present study, we developed a quantitative iTRAQ proteomics strategy to investigate differential protein expression in salivary glands and thus to further investigate both the development and the functional adjustment of salivary glands in the different feeding stages of female H. longicornis. Protein regulation during changes in salivary gland function in female ticks was systematically analyzed. The results of the present study showed that 615 salivary gland proteins were upregulated in partially fed females compared with unfed females (Additional file 4: Table S4).
Tick saliva contains a variety of functional proteins and peptides, such as vasodilators, anticoagulants, inhibitors of platelet aggregation, and immunomodulators that are essential for successful blood-feeding . In addition, it was found that proteins in the saliva of ticks also included heme/iron metabolism-related proteins, oxidation/detoxification proteins, a variety of proteinase inhibitors, and so on . To complete a blood meal at the feeding site ticks secrete proteins produced in the salivary glands that inhibit the host immune response [38, 39]. In the present study, the expression of a protein related to immunosuppressant protein (Da-p36), which can suppress the concanavalin A-induced in vitro proliferation of normal murine T-lymphocytes , was higher in the salivary glands of semi-engorged and engorged ticks than in those of unfed ticks. Although the change in this protein was not very stable among the three experiments, the data indicated that the expression of this protein was upregulated in the salivary glands of semi-engorged ticks compared with those of unfed ticks. Furthermore, the angiotensin-converting enzyme-like protein found in ticks may have functions similar to those of angiotensin-converting enzyme in mammals , which can control blood pressure by regulating fluid volume . This control would allow the tick to continuously feed on host blood.
We further found that two sequences annotated as Kunitz-type serine protease inhibitors were upregulated when ticks sucked blood . However, the homology of these two sequences with those of Kunitz-type serine protease inhibitors was relatively low. Therefore, we cannot reach any definite conclusions about changes in the expression of proteins with anticoagulant functions in tick salivary glands . We believe that the main reason why some small peptides such as antimicrobial peptides [45, 46] were not detected is that we used the FASP method. Although the FASP method can produce pre-labeled samples that are sufficiently clean for iTRAQ labeling and is thus currently an ideal pre-treatment method for iTRAQ, peptides smaller than 10 kDa (sometimes slightly larger than 10 kDa) are washed away by centrifugation before trypsin digestion. In addition, the sequences of functional peptides are very short, which ultimately causes the number of digested peptides that can be identified by MS to be less than that of peptides from digested proteins. Therefore, we did not obtain quantitative information about functional peptides.
The salivary glands of unfed female ticks are undeveloped and protein expression is low . When ticks were exposed to the host but not allowed to feed, some proteins are expressed when female ticks are stimulated to start feeding [9,10,11]. Once a tick starts sucking blood, salivary gland begins to develop. In the present study, we found that during constant tick blood-feeding the expression of large amounts of proteins increase, with the expression of thrombospondin upregulated in partially fed female salivary glands compared with unfed female salivary glands. We found that the thrombospondin sequence provided in GenBank (XP_002434432.1) showed the greatest homology with thrombospondin 3. Thrombospondin 1 has been confirmed to have antiangiogenic activity in mammals . Although thrombospondin 3 in mammals has not been found to have the same function as thrombospondin 1 , it is possible that thrombospondin 3 in ticks could have the same function as thrombospondin 1. The expression of the angiotensin-converting enzyme-like and thrombospondin 3 proteins did not differ significantly among the partially fed, semi-engorged, and engorged stage, which seems to indicate that these proteins are maintained at constant levels in the salivary glands during these periods.
Coat protein complexes I and II (COPI and COPII) can transport proteins between the rough endoplasmic reticulum and the Golgi apparatus , which is an essential process in protein synthesis. These two proteins were both upregulated in partially fed compared with unfed female ticks, indicating an acceleration of protein synthesis and transport in the partially fed ticks (Additional file 4: Table S4). In contrast, there was a significant decrease in the expression of COPI in engorged compared with semi-engorged ticks that was associated with a decrease in protein expression during salivary gland degeneration (Additional file 4: Table S4). In partially fed ticks, salivary gland cell division produces large numbers of new cells. Tubulin, a cytoskeleton component , was upregulated at this stage, which reflects the massive proliferation and growth of salivary gland cells.
When the salivary glands develop, cell mitosis and meiosis accelerate, and processes that transport cellular cargo become faster. Therefore, developing salivary gland cells require more microtubule motor proteins to transport proteins and organelles . The activity of dynein and kinesin is powered by ATP hydrolysis [52, 53]; these proteins can move along microtubules and transport cargo through interactions with tubulin . In the present study, we found that the expression of these two motor proteins was indeed upregulated in partially fed females compared with unfed females. However, at other blood-feeding stages, there were no differences in the expression of these two motor proteins, indicating that the expression and transport of total proteins in the other three stages may be relatively stable. Through RNAi, we proved that the reduction of these two proteins had a marked effect on the physiological activity of tick salivary glands.
Because a large number of cells are rapidly proliferating, the energy supply during salivary gland development is particularly important. Our results showed that a series of enzymes associated with ATP production, including various enzymes involved in glycolysis, the TCA cycle, the pentose phosphate pathway, and oxidative phosphorylation in respiratory chains, were clearly upregulated. Additional file 4: Table S4 shows the upregulated enzymes involved in glycolysis , such as 3-phosphoglycerate kinase, phosphoglycerate mutase, glucose-6-phosphate isomerase, and enolase. Significantly upregulated enzymes involved in the TCA cycle included pyruvate dehydrogenase, citrate synthase, isocitrate dehydrogenase, succinyl-CoA synthetase, succinate dehydrogenase, fumarase, and malate dehydrogenase. Through RNAi, we proved that a lack of citrate synthase and isocitrate dehydrogenase had a marked effect on the physiological activity of tick salivary glands. Furthermore, the enzyme 6-phosphogluconate dehydrogenase, which is involved in the pentose phosphate pathway, was also upregulated. The upregulation of all of these enzymes in ticks should lead to increased NADPH production for the synthesis of lipids, fatty acids, and nucleotides . In addition, the expression of some proteins involved in the respiratory chain protein complexes in the inner mitochondrial membrane, such as NADH-ubiquinone oxidoreductase, succinate dehydrogenase, cytochrome b-c1 complex subunit 2, and F0F1-type ATP synthase, increased significantly. The results indicate that when ticks enter the blood-feeding stage, the TCA cycle and other energy-metabolism pathways in their salivary gland cells become more active so that the ticks can immediately convert carbohydrates, fats, and proteins from the blood of the host into a large amount of ATP via the respiratory chain. Ultimately, most organs in ticks were able to develop rapidly after a short duration of blood-feeding.
The mortality of female ticks with dynein and kinesin (47.5% and 42.5%) RNAi was significantly higher than that of the control group (12.5%). In addition, with the significant decrease of engorged weight of female ticks after RNAi, the egg-laying rate, egg quality and hatching rate of female ticks also decreased significantly. The physiological changes of female ticks after RNAi of isocitrate dehydrogenase and citrate synthase were not as obvious as those after RNAi of dynein and kinesin. It is thought that a decrease in these two enzymes involved in the TCA cycle reduces the ATP provided for tick development, but the ticks can obtain ATP in other ways . The motor proteins dynein and kinesin play a decisive role in the transport of substances [58, 59] in the cell of tick, so their reduced expression after RNAi can affect the physiological activity of the tick even more.
A small number of salivary gland proteins appeared to be differentially expressed between semi-engorged and partially fed ticks. According to statistical analysis, 53 proteins, including 12 hypothetical proteins, were upregulated and 34 proteins, including 8 hypothetical proteins, were downregulated in semi-engorged ticks compared with partially fed ticks. The expression levels and change fold of proteins between these two stages are all small, which indicate the function of salivary glands in these two periods is relatively balanced and that the expression levels of the proteins are relatively stable. The female tick enters a rapid feeding stage after mating , and the body weight can be increased 100-fold when they are engorged . Previous studies have shown that the expression level of a large number of proteins in ovaries of female ticks changes after mating . In the present study, we found that the expression levels of a large number of proteins in salivary glands also changed after mating (Additional file 4: Table S4), which suggests that energy production should be increased in this stage. However, there were no significant increases in the expression of enzymes related to glycolysis, the TCA cycle, or the pentose phosphate pathway in the salivary glands at this stage. There was little change in the expression of proteins associated with energy metabolism. We speculate that during this period of time, to achieve rapid blood-feeding, the female tick primarily uses energy for development of the salivary glands and other organs, such as the midgut, and that the energy demand of the salivary glands does not increase as much.
The protein expression levels of dynein, kinesin, isocitrate dehydrogenase, and citrate synthase in the salivary glands of female ticks are upregulated. Quantitative results of RT-qPCR showed that the transcription level of these proteins was consistent with the protein level during the development of salivary glands, but the difference was only slight in fold change, which is a frequent occurrence . Although iTRAQ is a more ideal method for protein quantification, and the data generated by us have good reproducibility, it cannot guarantee the accurate quantification of each protein, which is caused by various reasons. Therefore, it may be necessary to use other more advanced quantitative methods to make up for the shortcomings of iTRAQ in the future.
In this study, the dynamic changes of salivary gland proteins expression in H. longicornis from unfed to engorgement were systematically analyzed. This will provide a deeper understanding of the synergistic patterns of proteins in the whole process of salivary gland development to degeneration. Meanwhile, through RNAi on dynein, kinesin, isocitrate dehydrogenase, and citrate synthase, the study investigated the functions and effects of the four proteins on ticks and their salivary glands. The research results provide new protein targets for the control of ticks and tick-borne diseases. The results also provide new protein targets for controlling ticks which can have profound implications for the prevention of tick-borne diseases.
Availability of data and materials
Data supporting the conclusions of this article are included within the article. The proteomic dataset has been submitted to ProteomeXchange via the PRIDE database (accession number: PXD013678).
isobaric tags for relative and absolute quantification
filter-aided sample preparation
high-energy collision-induced dissociation
false discovery rate
normalized spectrum abundance factor
Kyoto Encyclopedia of Genes and Genomes
green fluorescent protein
coefficients of variation
- COPI and COPII:
Coat protein complexes I and II
Estrada-Pena A, Jongejan F. Ticks feeding on humans: A review of records on human-biting Ixodoidea with special reference to pathogen transmission. Exp Appl Acarol. 1999;23:685–715.
Iwanaga S, Okada M, Isawa H, Morita A, Yuda M, Chinzei Y. Identification and characterization of novel salivary thrombin inhibitors from the ixodidae tick, Haemaphysalis longicornis. Eur J Biochem. 2003;270:1926–34.
Tang J, Fang Y, Han Y, Bai X, Yan X, Zhang Y, et al. YY-39, a tick anti-thrombosis peptide containing RGD domain. Peptides. 2015;68:99–104.
Kotál J, Stergiou N, Buša M, Chlastáková A, Beránková Z, Řezáčová P, et al. The structure and function of Iristatin, a novel immunosuppressive tick salivary cystatin. Cell Mol Life Sci. 2019;76:2003–13.
Islam MK, Tsuji N, Miyoshi T, Alim MA, Huang X, Hatta T, et al. The Kunitz-like modulatory protein haemangin is vital for hard tick blood-feeding success. PLoS Pathog. 2009;5:e1000497.
Artigas-Jerónimo S, Estrada-Peña A, Cabezas-Cruz A, Alberdi P, Villar M, de la Fuente J. Modeling modulation of the tick regulome in response to Anaplasma phagocytophilum for the identification of new control targets. Front Physiol. 2019;10:462.
Lejal E, Moutailler S, Šimo L, Vayssier-Taussat M, Pollet T. Tick-borne pathogen detection in midgut and salivary glands of adult Ixodes ricinus. Parasites Vectors. 2019;12:152.
Nunes PH, Bechara GH, Camargo-Mathias MI. Morphological changes in the salivary glands of Amblyomma cajennense females (Acari: Ixodidae) in different feeding stages on rabbits at first infestation. Exp Appl Acarol. 2008;45:199–209.
Mulenga A, Blandon M, Khumthong R. The molecular basis of the Amblyomma americanum tick attachment phase. Exp Appl Acarol. 2007;41:267–87.
Lew-Tabor AE, Moolhuijzen PM, Vance ME, Kurscheid S, Valle MR, Jarrett S, et al. Suppressive subtractive hybridization analysis of Rhipicephalus (Boophilus) microplus larval and adult transcript expression during attachment and feeding. Vet Parasitol. 2010;167:304–20.
Tirloni L, Kim TK, Pinto AF, Yates JR, da Silva Vaz I, Jr Mulenga A. Tick-host range adaptation: changes in protein profiles in unfed adult Ixodes scapularis and Amblyomma americanum saliva stimulated to feed on different hosts. Front Cell Infect Microbiol. 2017;7:517.
Weiss LB, Reuben Kaufman W. The relationship between ‘critical weight’ and 20-hydroxiecdysone in the female ixodid tick, Amblyomma hebraeum. J Insect Physiol. 2001;47:1261–7.
Nunes ET, Bechara GH, Saito KC, Denardi SE, Oliveira PR, Mathias MI. Morphological, histological, and ultrastructural characterization of degenerating salivary glands in females of the cattle-tick Rhipicephalus (Boophilus) microplus (CANESTRINI, 1887) (Acari: Ixodidae). Micron. 2005;36:437–47.
Friesen KJ, Kaufman WR. Salivary gland degeneration and vitellogenesis in the ixodid tick Amblyomma hebraeum: surpassing a critical weight is the prerequisite and detachment from the host is the trigger. J Insect Physiol. 2009;55:936–42.
Nodari EF, Roma GC, Furquim KC, De Oliveira PR, Bechara GH, Camargo-Mathias MI. Degenerative process and cell death in salivary glands of Rhipicephalus sanguineus (Latreille, 1806) (Acari: Ixodidae) semi-engorged female exposed to the acaricide permethrin. Microsc Res Tech. 2012;75:1012–8.
Mao H, Kaufman WR. Profile of the ecdysteroid hormone and its receptor in the salivary gland of the adult female tick, Amblyomma hebraeum. Insect Biochem Mol Biol. 1999;29:33–42.
Scopinho FKC, Bechara GH, Camargo MMI. Death by apoptosis in salivary glands of females of the tick Rhipicephalus sanguineus (Latreille, 1806) (Acari: Ixodidae). Exp Parasitol. 2008;119:152–63.
Bensaoud C, Aounallah H, Sciani JM, Faria F, Chudzinski-Tavassi AM, Bouattour A, et al. Proteomic informed by transcriptomic for salivary glands components of the camel tick Hyalomma dromedarii. BMC Genom. 2019;20:675.
Di Venere M, Fumagalli M, Cafiso A, De Marco L, Epis S, Plantard O, et al. Ixodes ricinus and its endosymbiont Midichloria mitochondrii: A comparative proteomic analysis of salivary glands and ovaries. PLoS ONE. 2015;23(10):e0138842.
Esteves E, Maruyama SR, Kawahara R, Fujita A, Martins LA, Righi AA, et al. Analysis of the salivary gland transcriptome of unfed and partially fed Amblyomma sculptum ticks and descriptive proteome of the saliva. Front Cell Infect Microbiol. 2017;7:476.
Zheng H, Yu Z, Zhou L, Yang X, Liu J. Seasonal abundance and activity of the hard tick Haemaphysalis longicornis (Acari: Ixodidae) in North China. Exp Appl Acarol. 2012;56:133–41.
Chen X, Xu S, Yu Z, Guo L, Yang S, Liu L, et al. Multiple lines of evidence on the genetic relatedness of the parthenogenetic and bisexual Haemaphysalis longicornis (Acari: Ixodidae). Infect Genet Evol. 2014;21:308–14.
Luo LM, Zhao L, Wen HL, Zhang ZT, Liu JW, Fang LZ, et al. Haemaphysalis longicornis ticks as reservoir and vector of severe fever with thrombocytopenia syndrome virus in China. Emerg Infect Dis. 2015;21:1770–6.
Tsuda A, Mulenga A, Sugimoto C, Nakajima M, Ohashi K, Onuma M. cDNA cloning, characterization and vaccine effect analysis of Haemaphysalis longicornis tick saliva proteins. Vaccine. 2001;19:4287–96.
Vale RD. The molecular motor toolbox for intracellular transport. Cell. 2003;112:467–80.
Kiyomitsu T, Cheeseman IM. Chromosome and spindle pole-derived signals generate an intrinsic code for spindle position and orientation. Nat Cell Biol. 2012;14:311–7.
Kay J, Weitzman PD. Krebs’ citric acid cycle: half a century and still turning. London: Biochemical Society; 1987. p. 25.
Wang H, Zhang X, Wang X, Zhang B, Wang M, Yang X, et al. Comprehensive analysis of the global protein changes that occur during salivary gland degeneration in female ixodid ticks Haemaphysalis longicornis. Front Physiol. 2019;22:1943.
Wiśniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6:359–62.
Rigbolt KT, Vanselow JT, Blagoev B. GProX, a user-friendly platform for bioinformatics analysis and visualization of quantitative proteomics data. Mol Cell Proteom. 2011;10(O110):007450.
Cao J, Wang T, Wang Q, Zheng X, Huang L. Functional insights into protein acetylation in the hyperthermophilic archaeon Sulfolobus islandicus. Mol Cell Proteom. 2019;18:1572–87.
Carne NA, Bell S, Brown AP, Määttä A, Flagler MJ, Benham AM. Reductive stress selectively disrupts collagen homeostasis and modifies growth factor-independent signaling through the MAPK/Akt pathway in human dermal fibroblasts. Mol Cell Proteom. 2019;18:1123–37.
Gong H, Umemiya R, Zhou J, Liao M, Zhang H, Jia H, et al. Blocking the secretion of saliva by silencing the HlYkt6, gene in the tick Haemaphysalis longicornis. Insect Biochem Mol Biol. 2009;39:372–81.
Kocan KM, Zivkovic Z, Blouin EF, Naranjo V, Almazán C, Mitra R, et al. Silencing of genes involved in Anaplasma marginale-tick interactions affects the pathogen developmental cycle in Dermacentor variabilis. BMC Dev Biol. 2009;9:42.
Verna E. More attention should be paid on the interpretation of gene expression data. World J Gastroenterol. 2012;18:3181–2.
Maritz-Olivier C, Stutzer C, Jongejan F, Neitz AW, Gaspar AR. Tick anti-hemostatics: targets for future vaccines and therapeutics. Trends Parasitol. 2007;23:397–407.
Tirloni L, Islam MS, Kim TK, Diedrich JK, Yates JR, Pinto AF, et al. Saliva from nymph and adult females of Haemaphysalis longicornis: a proteomic study. Parasites Vectors. 2015;8:338.
Wikel SK, Ramachandra RN, Bergman DK. Tick-induced modulation of the host immune response. Int J Parasitol. 1994;24:59–66.
Lu S, Soares TS, Vaz Junior IS, Lovato DV, Tanaka AS. Rmcystatin, a cysteine protease inhibitor from Rhipicephalus microplus hemocytes involved in immune response. Biochimie. 2014;106:17–23.
Bergman DK, Palmer MJ, Caimano MJ, Radolf JD, Wikel SK. Isolation and molecular cloning of a secreted immunosuppressant protein from Dermacentor andersoni salivary gland. J Parasitol. 2000;86:516–25.
Jarmey JM, Riding GA, Pearson RD, McKenna RV, Willadsen P. Carboxydipeptidase from Boophilus microplus: a “concealed” antigen with similarity to angiotensin-converting enzyme. Insect Biochem Mol Biol. 1995;25:969–74.
Ehlers MR, Riordan JF. Angiotensin-converting enzyme: new concepts concerning its biological role. Biochemistry. 1989;28:5311–8.
Cao J, Shi L, Zhou Y, Gao X, Zhang H, Gong H, et al. Characterization of a new Kunitz-type serine protease inhibitor from the hard tick Rhipicephalus hemaphysaloides. Arch Insect Biochem Physiol. 2013;84:104–13.
Ibrahim MA, Masoud HMM. Thrombin inhibitor from the salivary gland of the camel tick Hyalomma dromedarii. Exp Appl Acarol. 2018;74:85–97.
Sun T, Pan W, Song Y, Zhang J, Wang J, Dai J. Functional characterization of two defensins, HlDFS1 and HlDFS2, from the hard tick Haemaphysalis longicornis. Parasites Vectors. 2017;10:455.
Yada Y, Talactac MR, Kusakisako K, Hernandez EP, Galay RL, Andoh M, et al. Hemolymph defensin from the hard tick Haemaphysalis longicornis attacks Gram-positive bacteria. J Invertebr Pathol. 2018;156:14–8.
Slovák M, Hajnická V, Labuda M, Fuchsberger N. Comparison of the protein profiles of salivary gland extracts derived from three species of unfed and partially fed ixodid ticks analysed by SDS-PAGE. Folia Parasitol. 2000;47:67–71.
Haviv F, Bradley MF, Kalvin DM, Schneider AJ, Davidson DJ, Majest SM, et al. Thrombospondin-1 mimetic peptide inhibitors of angiogenesis and tumor growth: design, synthesis, and optimization of pharmacokinetics and biological activities. J Med Chem. 2005;48:2838–46.
Qabar AN, Bullock J, Matej L, Polverini P. Expression and characterization of novel thrombospondin 1 type I repeat fusion proteins. Biochem J. 2000;346:147–53.
Letourneur F, Gaynor EC, Hennecke S, Démollière C, Duden R, Emr SD, et al. Coatomer is essential for retrieval of dilysine-tagged proteins to the endoplasmic reticulum. Cell. 1994;79:1199–307.
Hachouf-Gheras S, Besson MT, Bosquet G. Identification and developmental expression of a Bombyx mori alpha-tubulin gene. Gene. 1998;208:89–94.
Gibbons IR. Dynein ATPases as microtubule motors. J Biol Chem. 1988;263:15837–40.
Schnitzer MJ, Block SM. Kinesin hydrolyses one ATP per 8-nm step. Nature. 1997;388:386–90.
Rice S, Lin AW, Safer D, Hart CL, Naber N, Carragher BO, et al. A structural change in the kinesin motor protein that drives motility. Nature. 1999;402:778–84.
O’Neill LA, Hardie DG. Metabolism of inflammation limited by AMPK and pseudo-starvation. Nature. 2013;493:346–55.
Ma L, Cheng Q. Inhibiting 6-phosphogluconate dehydrogenase reverses doxorubicin resistance in anaplastic thyroid cancer via inhibiting NADPH-dependent metabolic reprogramming. Biochem Biophys Res Commun. 2018;498:912–7.
Moraes J, Galina A, Alvarenga PH, Rezende GL, Masuda A, da Silva Vaz I Jr, et al. Glucose metabolism during embryogenesis of the hard tick Boophilus microplus. Comp Biochem Physiol A Mol Integr Physiol. 2007;146:528–33.
Cyr JL, Brady ST. Molecular motors in axonal transport. Cellular and molecular biology of kinesin. Mol Neurobiol. 1992;6:137–55.
Schnapp BJ, Reese TS. Dynein is the motor for retrograde axonal transport of organelles. Proc Natl Acad Sci USA. 1989;86:1548–52.
Weiss BL, Kaufman WR. Two feeding-induced proteins from the male gonad trigger engorgement of the female tick Amblyomma hebraeum. Proc Natl Acad Sci USA. 2004;101:5874–9.
Reuben Kaufman W. Gluttony and sex in female ixodid ticks: how do they compare to other blood-sucking arthropods? J Insect Physiol. 2007;53:264–73.
Xavier MA, Tirloni L, Pinto AFM, Diedrich JK, Yates JR 3rd, Mulenga A, et al. A proteomic insight into vitellogenesis during tick ovary maturation. Sci Rep. 2018;8:4698.
Lu P, Vogel C, Wang R, Yao X, Marcotte EM. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2007;25:117–24.
The authors would like to thank Ecology Post-doctoral Research Station of Hebei Normal University.
This work was supported by the Natural Science Foundation for Excellent Youth Scholars of Hebei Province of China (No. C2017205135), Natural Science Fund for Distinguished Young Scholars of Hebei Normal University (No. L2017J04), Hebei Normal University Doctor Initial Science Foud (No. L2017B13). Youth Foundation of Hebei Province Department of Education Fund (No. QN2017092).
Ethics approval and consent to participate
This study was approved by the Animal Ethics Committee of Hebei Normal University (protocol number: 165031) as complying with the animal protection law of the People’s Republic of China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Statistics for protein quantification information of 115:114.
. Statistics for protein quantification information of 116:114.
Statistics for protein quantification information of 117:114.
Statistics for protein quantification information of identified proteins shared among all three biological replicates.
Overview of the reproducibility of iTRAQ quantitative proteomics. a-c Correlation coefficients for the partially fed:unfed abundance ratios among 3 replicates; d-f Correlation coefficients for the semi-engorged:unfed abundance ratios among 3 replicates; g-i Correlation coefficients for the engorged:unfed abundance ratios among 3 replicates.
Raw data obtained by cluster analysis using the GProX platform.
GO functional annotations for the differentially expressed proteins in Clusters 1–5. BP: biological process; CC: cellular component; MF: molecular function.
Raw data for the KEGG pathway enrichment analysis of differentially expressed proteins in Cluster 1.
Raw data for the KEGG pathway enrichment analysis of differentially expressed proteins in Cluster 2.
Raw data for the KEGG pathway enrichment analysis of differentially expressed proteins in Cluster 3.
Raw data for the KEGG pathway enrichment analysis of differentially expressed proteins in Cluster 4.
Raw data for the KEGG pathway enrichment analysis of differentially expressed proteins in Cluster 5
Cq values of the RT-qPCR for 4 mRNA in 4 different stages.
Cq values of the RT-qPCR in the salivary glands of female ticks.
Cq values of the RT-qPCR in the whole female ticks.
About this article
Cite this article
Ren, S., Zhang, B., Xue, X. et al. Salivary gland proteome analysis of developing adult female Haemaphysalis longicornis ticks: molecular motor and TCA cycle-related proteins play an important role throughout development. Parasites Vectors 12, 613 (2019). https://doi.org/10.1186/s13071-019-3864-2