Abstract
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.
Introduction
Ticks transmit a wide range of zoonosis, such as tick-borne encephalitis (Grabowski et al., 2017), Lyme disease (Ehrmann et al., 2018), tick-borne spotted fever (Han et al., 2018), anaplasmosis (Parola and Raoult, 2001; Eisen and Eisen, 2018), and babesiosis (Kotal et al., 2015). These pathogens are transmitted via tick saliva to their vertebrate host when ticks feed on host. Transovarial transmission from adult female tick to its offspring may also occur in some tick species and tick-borne diseases (Danielova et al., 2002; Harris et al., 2017; Umemiya-Shirafuji et al., 2017; Jongejan et al., 2018; Moore et al., 2018). Tick salivary glands secrete immunosuppressant proteins including P36 protein that regulate the host immune system during blood feeding (Bergman et al., 2000; Alarcon-Chaidez et al., 2003; Konnai et al., 2009; Anatriello et al., 2010; Wang et al., 2017). Although P36 homolog proteins from Ixodid ticks are not highly conserved, all ticks share the same key site and structural features with conserved region “IDKGMLSPF”. Based on 3D structures of P36 from Dermacentor andersoni, a predicted conserved antigenic region was identified, located within the exposed loop with a potential of binding immunomodulating ligands, including glycerol and lactose (Oyugi et al., 2018). In our previous study, we found that RH36 from Rhipicephalus haemaphysaloides, P36 homolog, suppressed the T-lymphocyte mitogen-driven proliferation of splenocytes in vitro and in vivo, as well as the expression of several cytokines, including IL-2, IL-12, and TNF-α (Wang et al., 2017). P36 homologs expressed mainly in the salivary glands of partially fed ticks (Alarcon-Chaidez et al., 2003; Konnai et al., 2009) and RH36 protein production correlated with blood-feeding success and oviposition (Wang et al., 2017). Little information is known about the expression and interaction of RH36 molecules in the salivary glands and reproductive organs, including the ovary, while the mechanism by which the RH36 gene regulates ovary development is also unknown.
Ticks with a huge reproductive capacity produce hundreds of eggs, while the key of egg maturation is vitellogenesis, a process in which massive amounts of vitellin (Vn) precursors (vitellogenin, Vg) are expressed in fat bodies and midguts, and then secreted into hemolymph, and subsequently accumulated by developing oocytes via vitellogenin receptor-mediated clathrin-dependent endocytosis (Rosell and Coons, 1992; Esmeralda Parra-Peralbo, 2011; Khalil et al., 2011). After blood feeding, tick vitellogenesis is activated by a cooperative action of nutritional amino acid/target of rapamycin/S6 kinase (AA/TOR/S6K) (Umemiya-Shirafuji et al., 2012) and steroid hormone (20-hydroxyecdysone, 20E) pathways, closely correlated with mating (Thompson et al., 2005). In addition, heme and iron metabolism during tick blood feeding are essential for tick reproduction and egg fertility (Hajdusek et al., 2009; Galay et al., 2013; Perner et al., 2016). Genome and proteomics analyses of the parasitic processes seem to be unique in ticks, including prolonged feeding, new methods of hemoglobin digestion, heme detoxification, vitellogenesis and transmission of pathogens (Gulia-Nuss et al., 2016; Barrero et al., 2017). We hypothesized that the RH36 molecules participated in the integration between reproduction and nutrition and investigated the mechanism of RH36 in oviposition by proteomics analyses.
The significance of our study is revelation of the molecular mechanism of immunosuppressor RH36 modulating ovary development and oviposition in R. haemaphysaloides ticks. The main hypothesis addressed in this study was that RH36 may regulate tick vitellogenesis, while heat shock protein 70 (HSP70) was found to be down-regulated in the RH36 gene-silenced ticks by proteomic analysis, which regulated tick blood feeding and the expression of vitellogenin. This study illustrates the critical function of RH36 in the process of ovary development and tick reproduction, which are required for blood meal completion.
Materials and Methods
Ethics Statement
The protocols (shvri-ra-2017070578) were approved by the Institutional Animal Care and Use Committee of the Shanghai Veterinary Research Institute, and authorized by the Animal Ethical Committee of Shanghai Veterinary Research Institute.
Ticks and Sample Preparation
Ticks were obtained from a colony of R. haemaphysaloides maintained on the ears of New Zealand white rabbits at Shanghai Veterinary Research Institute, CAAS. Off-host ticks were kept in an incubator at 25°C and 95% humidity. To detect dynamic expression of RH36 and other proteins, midguts, salivary glands, ovaries, hemolymph, and fat bodies were dissected from unfed female ticks, 5-day-fed female ticks, 7-day-fed female ticks, engorging ticks and ticks on day 3 or 10 after engorgement. Furthermore, unfed female adult ticks (N = 80 females per group, three replicates) were microinjected with RH36 and luciferase double-stranded RNA (dsRNA), respectively, as described previously (Wang et al., 2017). Two individual five-day fed female ticks from each group were collected to characterize RH36 mRNA levels after RH36 RNA interference (RNAi). Blood-feeding ticks from RH36 or luciferase gene silencing groups (14 ticks per group) were removed from rabbits at 5 (indicated by RF5 and LF5, respectively) or 7 days (indicated by RF7 and LF7). Engorged ticks (14 ticks per group) were also collected from RH36 or luciferase gene silencing groups at day 0 (denoted by RE0 and LE0), day 3 (denoted by RE3 and LE3), day 7 (denoted by RE7 and LE7) and day 10 (denoted by RE10 and LE10) after engorgement. The ovaries, hemolymph, fat bodies and salivary glands were dissected from 7 ticks per group for RNA or protein extraction and washed in PBS to remove hemolymph -related cells.
Total RNA was extracted from ovaries using TRIZOL (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. Total RNA quality was evaluated using Agilent 2100 Bioanalyzer RNA Nano Chip (Agilent Technologies, Santa Clara, CA, USA). Tick ovaries with SDT buffer (containing 4% SDS, 100 mM DTT and 100 mM Tris-HCl, pH 7.6) and quartz sand (another 1/4 inch ceramic bead MP 6540-424 for tissue samples) were homogenized by MP homogenizer (24 × 2, 6.0M/S, 60s, twice) and subsequently sonicated and boiled for 15 min. After it was centrifuged at 10,000 rpm for 20 min, the supernatant was filtered with 0.22 μm filters. Total protein was also extracted from tick salivary glands, fat bodies and hemolymph according to the above instruction, respectively. The filtrate was quantified with the BCA Protein Assay Kit (Bio-Rad, USA). These RNA and protein samples were stored at −80°C until use.
Antibody Against RH36 Preparation and Expression Level of RH36 Protein by Western Blot
The B cell linear epitope of RH36 protein was predicted using the Immune Epitope Database (IEDB) of http://tools.iedb.org/bcell/result/. The peptide (LTDDYWKQGEHPSGEYPIS) was selected randomly and synthesized with peptide purify of at least 95% and coupled with KLH by GL Biochem (Shanghai, China), and emulsified with an equal volume of Freund's complete adjuvant (Sigma-Aldrich, St. Louis, MO, USA) used to intraperitoneally inject Balb/c mice. Two booster injections of RH36 peptide were given with the same procedure but using Freund's incomplete adjuvant. One week after the final immunization, sera were collected following centrifugation. Antibodies against RH36 peptides were used to detect expression levels by Western blot. According to standard techniques, total proteins were separated on a 12% SDS-PAGE gel and transferred to PVDF membranes (0.45 μm, Millipore, Billerica, MA, USA), and the membranes were immunoblotted with RH36 antibody (1:200). The monoclonal antibody of β-actin against mouse (Proteintech, Rosemont, USA, Cat. No. 60008-1-Ig) was used as the control. IgG goat anti-mouse or goat anti-rabbit antibodies conjugated with HRP were used for secondary antibodies (1:8,000, Bethyl Laboratories, Inc., Montgomery, Texas, USA). The membranes were visualized using automatic chemiluminescence image analysis system (Tanon-5200, Tanon Science & Technology Co., Ltd, Shanghai, China). Densitometry used to analyze indicated protein bands was determined by Image J (NIH), normalized with the control (Schneider et al., 2012).
Antibody Against Vitellin Preparation and Expression Level of Vitellin Protein by Western Blot
Crude eggs were homogenized and subsequently centrifuged at 10,000 rpm for 20 min, the supernatant was applied onto a Sepharose CL-4B gel filtration column and diethyl-amino-ethyl (DEAE)-cellulose column in turn (Yang et al., 2015). Polyclonal antibody and monoclonal antibody against vitellin was generated from Balb/c mice (Wang et al., 2017). Splenocytes from immunized mice were fused to myeloma cells, SP2/0, using polyethylene glycol. Hybridomas were selected in HAT medium and the hybridoma culture supernatants were screened by enzyme-linked immunosorbent assay (ELISA) in 96-microtiter plates and Western blot (Yang et al., 2015).
Expression of vitellin was detected in the ovaries, hemolymph, fat bodies and salivary glands by Western blot according to the above description. The monoclonal antibodies of β-actin (Proteintech, Cat. No. 60008-1-Ig) and GAPDH (Proteintech, Cat. No. 60004-1-Ig) against mouse were used as the control. The bands were visualized using automatic chemiluminescence image analysis system (Tanon-5200, Tanon Science & Technology Co., Ltd., Shanghai, China). Heat shock protein 70 rabbit antibody (HSP70, Cell Signaling Technology, Boston, USA, Cat. No. 3177) was also used to detect the protein expression level in the different organs during different blood-feeding ticks by Western blot. The anti-GAPDH antibody against rabbit (Cell Signaling Technology, Cat. No. 2118) was used as the control.
Immunofluorescence Assays (IFA)
For immunofluorescence, ovaries were obtained from adult female ticks at 10 days after engorgement as described above, fixed with 4% paraformaldehyde in 0.2 M sodium cacodylate buffer, and dehydrated in a graded series of ethanol and embedded in paraffin. Sections (4 μm) were prepared and mounted on glass slides. Paraffin was removed from the sections with xylene and the sections were hydrated by successive 5-min washes with a graded series of 100, 80, 75, and 65% ethanol. The slides were treated with EDTA antigen repair buffer for 30 min at 37°C, washed with PBST and incubated with 2% bovine serum albumin (BSA; Sigma-Aldrich) in PBST for 30 min at room temperature. The slides were then incubated for 14 h at 4°C with primary antibodies against vitellin diluted 1:5–1:200 in 2% BSA mixture and developed for 30 min with goat-anti-mouse IgG conjugated with FITC (Sigma-Aldrich) (diluted 1:500 in 2% BSA mixture) after 3 washes in PBST buffer. The slides were incubated for 20 min with DAPI (1:200, 2% BSA dilution) after 3 washes in PBST buffer and mounted in Fluoromount Aqueous Mounting Medium (Sigma-Aldrich). Samples were observed and photographed with 10X objective, 1.0X zoom and the excitation of 488 nm (FITC) and 405 nm (DAPI) under a Zeiss inverted fluorescence microscope (Zeiss, Germany).
Proteomic Data Collection and Bioinformatics Analysis
Quantitative proteomic analysis by tandem mass tag (TMT) technology was performed by Shanghai Applied Protein Technology (APT, Shanghai, China). Each 100 μg peptide mixture digested with trypsin (Promega, Madison, WI, USA) was labeled using TMT reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer's instructions. The 6-plex TMT Labeling were fractionated into 15 strong cation exchange (SCX). LC-MS/MS analysis (1H) was performed on a Q Exactive mass spectrometer (Thermo Fisher Scientific) that was coupled to Easy nLC (Proxeon Biosystems, now Thermo Fisher Scientific) for 60 min. The mass spectrometer was operated in positive ion mode. MS data was acquired using a data-dependent top 10 method dynamically choosing the most abundant precursor ions from the survey scan (300–1,800 m/z) for HCD fragmentation. Automatic gain control (AGC) target was set to 3e6, and maximum inject time to 10 ms. Dynamic exclusion duration was 40.0 s. Survey scans were acquired at a resolution of 70,000 at m/z 200 and resolution for HCD spectra was set to 17,500 at m/z 200 and isolation width was 2 m/z. Normalized collision energy was 30 eV and the underfill ratio, which specifies the minimum percentage of the target value likely to be reached at maximum fill time, was defined as 0.1%. The instrument was run with peptide recognition mode enabled. And results were compared to a Uniprot database containing all sequences from Ixodidae (January 1, 2017) using MASCOT engine (Matrix Science, London, UK; version 2.2) embedded into Proteome Discoverer 1.4. The MASCOT identified parameters were set in Supplementary Table 1. Peptide identification was validated using false discovery rate (FDR ≤ 0.01). The standard of screening differentially expressed proteins was Fold change >1.2 and a P-value < 0.05 was considered statistically significant. Sequences of differentially expressed proteins were also searched using blastp tool against the online Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://geneontology.org/) and were subsequently mapped to pathways in KEGG. GO and KEGG pathway enrichment analysis were applied based on the Fisher's exact test. All assays were performed in three replicates.
RNA Sequencing and Analysis
After total RNA was treated with DNase I, mRNA was enriched by oligo (dT) with beads and fragmented into fragments of 100–400 bp. First-strand cDNA was synthesized from the mRNA fragments. And second-strand cDNA was then synthesized using DNA polymerase I and RNaseH. Short fragments were purified, end-repaired, connected with signal nucleotide A (adenine) addition and ligated with library-specific adaptors. The cDNA library was constructed and tested in qualification by Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System, and then sequenced using Illumina HiSeq 4000.
Tick ovary transcriptome de novo assembly was performed by Beijing Genomics Institute (BGI, Shenzhen, China). After sequencing, low-quality, adaptor-polluted and high content of unknown base (N) reads were filtered to get clean reads. Clean reads were assembled by Trinity assembles software to obtain high-quality contigs (Yu et al., 2015). T gicl was used to perform clustering and eliminate redundant data in the assembled transcripts to obtain unique genes. Subsequently, expression analysis, heterozygous SNP detection, and functional annotation were performed. Differentially expression genes (DEGs) were further analyzed with gene ontology functional enrichment and KEGG pathway functional enrichment. The false discovery rate (FDR ≤ 0.001) was used to determine the threshold P value in multiple tests. Results were the mean of three replicates.
HSP70 RNA Interference (RNAi) for Gene Knockdown in Ticks
The dsRNA was synthesized using the T7 RiboMAX™ RNAi System (Promega, Madison, WI, USA) according to the manufacturer's recommendations. Oligonucleotide primers containing T7 promoter sequences (in italics at 5′-end) were listed in Supplementary Table 2. A dsRNA targeting luciferase used as a negative control was subjected to the same PCR amplification protocol using luciferase-specific primers (Yu et al., 2013). These dsRNA were maintained at −80°C until use. Unfed adult ticks (N = 50 females per group, two independent groups) were microinjected with 0.5 μL HSP70 dsRNA (about 1 μg RNA) at the base of the fourth right leg of the ventral surface of ticks. Control ticks were microinjected with unrelated luciferase dsRNA. After microinjection, ticks were maintained in an incubator at 25°C and 95% humidity for 24 h and then allowed to feed on rabbits' ears. Four female ticks per group were collected at 5 days of blood feeding for RNA and protein extraction to characterize gene knockdown by real-time quantitative PCR and Western blot with respect to luciferase control. Remaining ticks were allowed to feed until engorgement, which tick mortality, engorgement rate and engorgement weight was determined in individual engorged female ticks.
Determination of Tick mRNA Levels by Real-time Quantitative PCR
The primers and amplicon length of target genes used for real-time quantitative PCR are listed in Supplementary Table 3. The first-stranded cDNA was synthesized from 1 μg RNAs with the PrimeScript RT reagent kit (Perfect Real Time) (Takara, Shiga, Japan) in the following program: 37°C for 15 min, 85°C for 7 s, and finally 4°C. Quantitative real-time PCR was performed using SYBR Premix Ex Taq (Takara, Japan) with QuantStudio 5 (ABI, Waltham, USA) in 20 μL reaction mixtures including 0.4 μL ROX Dye II. The reactions were incubated for 30 s at 95°C, followed by 40 cycles at 95°C for 5 s, 60°C for 34 s, and analysis of a melting curve. Tick elongation factor 1α (ELF1α) was used as an internal control because of its stability in R. microplus and R. appendiculatus ticks (Nijhof et al., 2009). All reactions were performed in triplicate. The 2−ΔCt method was used to calculate relative expression levels.
Statistical Analysis
The body weights of engorged ticks were expressed as the mean ± standard error of mean (SEM) and analyzed using Student t-test. Differences in the attachment rate, engorgement rate, and death rate of ticks in different groups were tested using the Chi-square test. Significance of quantitative real-time PCR results was expressed as means ± SEM and determined with the one-way analysis of variance (ANOVA), and statistical analysis was conducted using software SPSS version 20.0. P < 0.05 and 0.01 were considered to be statistically significant (*) and extremely significant (**), respectively.
Results
The Immunomodulatory Protein RH36 Regulated the Complex Changes of Vitellin
RH36 protein was expressed exclusively in the salivary glands, not in the midguts, ovaries or fat bodies of 5-day-fed ticks and engorged ticks (Figure 1A). Additionally, RH36 protein was expressed at higher levels in 7-day-fed and saturated ticks and subsequently decreased after engorgement, which was highly consistent with the time of tick ovarian development (Figure 1B). To confirm our hypothesis that RH36 modulating tick ovary development via vitellogenesis, four vitellogenin genes (named RHVg1, RHVg2, RHVg3, and RHVg4, respectively) were cloned from ovarian transcriptome from the post-engorgement R. haemaphysaloides female ticks (Supplementary Figure 1). These sequences are highly homologous in Haemaphysalis longicornis and Dermacentor variabilis (Supplementary Figure 2). RHVg1, RHVg2, RHVg3, RHVg4 sequences were available in the GenBank databases under the accession numbers: MK584634, MK584635, MK584636, MK584637, respectively. In addition, RHVg1 (Supplementary Figure 3A), RHVg2 (Supplementary Figure 3B) and RHVg3 (Supplementary Figure 3C) genes were widely distributed in the different development stages of ticks and RHVg4 (Supplementary Figure 3D) was detected only in engorged adults not larvae and nymphs by real-time quantitative PCR and semi-quantitative PCR (Supplementary Figure 3E). What's more, RHVg2 and RHVg4 was expressed in the fat body and ovary of engorged adult ticks, respectively (Supplementary Figure 3F).
Figure 1
Vitellin was purified from crude egg extraction by gel filtration (Sepharose CL-4B) and ion exchange chromatography on DEAE-cellulose and identified by possible polypeptide bands with molecular weights of 43.0, 47.3, 64.4, 72.6, 79, 85.9, 94, 102.3, 108.7, 136.5, and 161.5 kDa, respectively (Supplementary Figure 4A), in which two polypeptides (136.5, 72.7 kDa) were identified by LC-MS/MS (Supplementary Tables 4, 5). Monoclonal antibodies were generated and had specific immunological reaction with two polypeptides (136.5, 72.7 kDa) of vitellin (Supplementary Figure 4B). Vitellin proteins were highly expressed in fat body, hemolymph and ovary of adult ticks after engorgement (Figure 2).
Figure 2
After RH36 gene silencing, there was no difference for four vitellogenin genes in the fat body of adult ticks (Figures 3A,C,D). There was one exception: RHVg2 gene was up-regulated significantly at 7 days after engorgement (p < 0.01, Figure 3B). In the ovaries, although there was no difference in the transcript level of RHVg1 (Figure 4A), RHVg2 (Figure 4B) and RHVg3 (Figure 4C) genes from blood feeding to post-engorgement, RHVg4 gene was down-regulated significantly at 7 days or 10 days post engorgement (p < 0.01, Figure 4D). In addition, vitellin (the band of 136 kDa) protein decreased in the ovary of adult ticks at 3 or 10 days after engorgement, but not at 5 days of blood feeding (Figure 4E). Compared to mature oocytes from luciferase gene-silenced adult ticks, immunofluorescence showed that vitellin protein was less expressed in the cytoplasm of immature oocytes in response to RH36 gene-silencing (Figure 4F).
Figure 3
Figure 4
The Importance of HSP70 in the Ovary Development in Response to RH36 Depletion
To identify which proteins regulate ovary development, ovaries from adult female ticks at 5 days blood feeding and 10 days post-engorgement were selected to conduct TMT proteomics analysis. Because of significant difference in total proteins found in ovaries between blood feeding and post-engorgement (Supplementary Figure 5), the TMT Mass-Tagging mass spectrometric (2, 6-plex) results were analyzed using MASCOT engine version 2.2. Ovary proteomics analysis resulted in the identification of 5,326 and 1,685 unique proteins in ticks at 5 days feeding and 10 days post-engorgement respectively. In the 5-day-fed ovaries, there were 44 proteins statistically significant (p < 0.05 and fold change > 1.2). Among these proteins, 25 and 19 proteins were up- and down-regulated, respectively (Figure 5A and Supplementary Table 6). In the day 10 post-engorgement ovaries, RH36 gene knockdown resulted in 208 proteins with statistically significant production differentiation (p < 0.05 and fold change > 1.2). Among these proteins, 74 were up-regulated and 134 were down-regulated (Figure 5A and Supplementary Table 7). According to GO and KEGG analyses, differentially expressed proteins mainly correlated with cellular process, developmental process and metabolic process with function of lipid transporter activity, catalytic activity and binding in the ovaries at 5 days feeding (Figure 5B) and 10 days post-engorgement (Figure 5C), in which vitellogenin (Supplementary Figure 6A) and HSP70 (Supplementary Figure 6B) involved in protein export and protein processing in endoplasmic reticulum, respectively (Figures 5D,E).
Figure 5
RH36 Gene Regulated the Expression of HSP70 Protein
Since HSP70 protein was very important in tick blood-sucking and ovary maturation, we obtained HSP70 gene sequence (GenBank accession numbers: MK584638) from the ovary transcriptome (Table 1) and found that the open reading frame (ORF) of HSP70 was 1,980 bp long, encoding a deduced protein with 660 amino acids in length, a predicted molecular mass of 72.5 kDa and a pI of 5.32. It was also found that HSP70 protein included three signature sequences and the extreme C-terminal endoplasmic reticulum targeting sequence (KDEL, Supplementary Figure 7A). As shown in Supplementary Figure 7B, HSP70 from R. haemaphysaloides was highly conserved with over 83% identity in arthropods and mammals, in which it was highly homologous in Haemaphysalis flava and Ixodes scapularis with the identity of 96.21 and 91.06%, respectively. HSP70 protein was widely distributed in the hemolymph, ovary, salivary gland and fat body of post-engorgement ticks, with highest levels in the ovary at 7 days blood feeding and in the circulating hemolymph at 3 days post-engorgement (Figure 6).
Table 1
| Unigenes | Subject_id | Gene Name | Species | Identity | Align length | Miss_match | E_value | Score |
|---|---|---|---|---|---|---|---|---|
| Unigene10768_All | ACA84007.1 | Heat shock 70 kDa protein 5 | Haemaphysalis longicornis | 100 | 21 | 0 | 1.01E-11 | 45.4394 |
| CL2077.Contig1_All | ACA84007.1 | Heat shock 70 kDa protein 5 | Haemaphysalis longicornis | 97.51 | 401 | 10 | 0 | 761.14 |
| Unigene11678_All | XP_013782730.1 | Heat shock 70 kDa protein 4-like | Limulus polyphemus | 53.91 | 831 | 383 | 0 | 869.766 |
| Unigene13971_All | P19120.2 | Heat shock 70 kDa protein 8 | Bos taurus | 93.68 | 174 | 11 | 3.24E-86 | 321.627 |
| Unigene16158_All | XP_012264348.1 | Heat shock 70 kDa protein cognate 5 | Athalia rosae | 79.97 | 694 | 139 | 0 | 1061.21 |
| Unigene18598_All | XP_001012263.1 | Heat shock 70 kDa protein | Tetrahymena thermophila | 76 | 75 | 18 | 6.00E-24 | 114.39 |
| Unigene19775_All | XP_007453496.1 | Heat shock 70 kDa protein | Lipotes vexillifer | 61.64 | 219 | 84 | 3.08E-68 | 263.462 |
| Unigene23136_All | XP_013788178.1 | Heat shock 70 kDa protein 12A-like | Limulus polyphemus | 72.4 | 616 | 170 | 0 | 920.613 |
| Unigene48580_All | AIS39468.1 | Heat shock 70 kDa protein | Haemaphysalis flava | 100 | 20 | 0 | 2.31E-10 | 45.4394 |
| Unigene48580_All | DAA34432.1 | Heat shock 70 kDa protein | Amblyomma variegatum | 95 | 20 | 1 | 5.35E-10 | 41.5874 |
| Unigene56143_All | XP_010125825.1 | Heat shock 70 kDa protein | Chlamydotis macqueenii | 88.61 | 79 | 9 | 3.87E-34 | 148.288 |
| CL434.Contig1_All | XP_014239635.1 | Major heat shock 70 kDa protein Ba-like | Cimex lectularius | 78.74 | 635 | 135 | 0 | 1012.29 |
| CL2742.Contig1_All | XP_013778668.1 | Heat shock 70 kDa protein cognate 4 | Limulus polyphemus | 89.68 | 649 | 67 | 0 | 1162.52 |
List of heat shock 70 kDa protein in the tick ovaries transcriptome.
Figure 6
After RH36 gene knocking down, there was no difference for HSP70 gene expression in the fat body, salivary gland, midgut, ovary and hemolymph of 5-day-fed and 7-day-fed adult ticks (Figures 7A,B, respectively). On the day of engorgement, HSP70 gene was down-regulated significantly in the ovary, midgut and salivary gland of adult ticks (p < 0.05, Figure 7C), coinciding with RH36 down-regulation expression. At 3 days post-engorgement, HSP70 gene also decreased markedly in the fat body, ovary, midgut, salivary gland, and hemolymph of adult ticks (p < 0.05, Figure 7D). In addition, HSP70 gene expression reduced observably in the ovary of adult ticks at 10 days post-engorgement (p < 0.01, Figure 7E). From tick blood-sucking to engorgement, HSP70 protein level decreased in the ovaries of adult ticks from RH36 RNAi groups by Western blot (Figure 7F). These results suggested that RH36 promoted expression of HSP70 protein in the post-engorged adult ticks.
Figure 7
Impacts of HSP70 Gene Knockdown on Tick Blood Feeding
HSP70 mRNA decreased by 73.86% (Figure 8A) and HSP70 protein was significantly down-regulated after microinjecting HSP70 dsRNA (p < 0.05, Figure 8B). After HSP70 gene silencing, mean of 48 h attachment rate and average engorgement rate was 52.22 and 2.56%, respectively. Meanwhile, these parameters showed 88.89 and 80.39% in the luciferase RNAi groups, respectively, which revealed that down-regulation of HSP70 gene decreased tick blood feeding and engorgement rate (Table 2). Although ticks attached to rabbit ears after HSP70 gene silencing (Figure 8C), the engorged body weight of ticks decreased and death rate increased compared to the luciferase RNAi groups (p < 0.05, Table 2 and Figure 8D), which shows that HSP70 gene facilitates abundance of tick blood feeding and protects ticks from death. After HSP70 gene knockdown, two vitellogenins (RHVg1 and RHVg2, Figure 9A) and vitellin (Figure 9B) were observably down-regulated (p < 0.05), indicating that HSP70 protein regulated expression of vitellogenin proteins via promoting tick blood feeding.
Figure 8
Table 2
| Groups | 48 h attachment rate (%)a | Engorgement rate (%)a | Engorged body weight (mg)b | Death rate (%)a |
|---|---|---|---|---|
| HSP70 | 52.22 (47/90) | 2.56 (2/78) | 208.95 ± 201.53 (2 ticks) | 97.44 (76/78) |
| L | 88.89 (80/90) | 80.39 (41/51) | 324.45 ± 85.24 | 19.61 (10/51) |
Effect of HSP70 gene knockdown on tick blood feeding and death rate.
Significantly different using the chi-square test.
Significantly different using Student's t-test.
Figure 9
RH36 Gene Played an Important Role in the Integration of Nutrition Uptake and Reproduction
As shown in Figure 10, RH36 gene expressed-uniquely in salivary gland controlled tick blood feeding and modulated tick ovary development and oviposition by regulating the expression of vitellin in ovary, while HSP70 protein was down-regulated significantly in the ovaries of post-engorged ticks after RH36 gene silencing by real-time quantitative PCR and Western blot (p < 0.05). HSP70 gene positively modulated tick blood feeding and ovary development via nutrition uptake. Therefore, we speculated that RH36 gene silencing hindered tick blood feeding and oocytes maturity by inhibiting vitellin uptake mediated by HSP70 protein. In addition, HSP70 gene silencing considerably increased tick death rate during blood feeding.
Figure 10
Discussion
RH36 not only is an immunomodulatory protein that suppressed the proliferation of splenocytes and expression of several cytokines such as IL-2, IL-12, and TNF-α, but also modulates blood feeding and oviposition in hard ticks (Wang et al., 2017). Moreover, RH36 was induced by blood feeding and expressed specifically in salivary glands, but not in other organs such as ovary and hemolymph (Figure 1). How did RH36 regulate tick ovary development indirectly? We speculated that RH36 was involved in the integration of nutrition and reproduction in ticks, in which vitellogenin is critical to tick ovary maturation since it serves as the nutritional source for tick survivability and reproduction (Xavier et al., 2018). Vitellogenin and vitellin are structurally, biochemically and immunologically similar in the majority of insects, including ticks, which results in McAb against vitellin recognized vitellogenins (Kaufman, 2004; Yang et al., 2015). Four vitellogenin genes (Supplementary Figure 1) and McAb against vitellin (Figure 4) were obtained to test the dynamics of vitellogenesis. Vitellogenin genes were widely distributed in fat body, midgut, hemolymph and ovaries of fed and engorged adult ticks (Figure 2 and Supplementary Figure 3). It is reported that vitellogenin is synthesized in the fat body and midgut of female ticks after mating and transported through the hemolymph, captured by surface receptors called vitellogenin receptors (VgRs), and endocytosed into developing oocytes within the ovaries (Rosell and Coons, 1992; Khalil et al., 2011; Mitchell et al., 2019). RH36 gene did not regulate transcripts of RHVg1, RHVg3, and RHVg4 genes in the fat body of adult ticks (Figures 3A,C,D), but RHVg2 gene in fat body was up-regulated at 7 days post-engorgement (Figure 3B), which may be stimulated by down-regulation of vitellin protein in the ovaries of 7-day-engorged ticks (Figure 4D). RHVg4 gene and vitellin protein decreased significantly in the post-engorged ovaries after RH36 gene silencing (p < 0.05, Figures 4A–E), and vitellin protein was less expressed in the cytoplasm of immature oocytes compared to mature oocytes at 10 days after engorgement (Figure 4F), which suggested RH36 gene silencing attenuated vitellin expression in ovaries at post-engorgement. Due to the absence of available genome for Rhipicephalus haemaphysaloides ticks, some ovarian-specific proteins may be lost by searching Ixodidae database resulting in less proteins correlating with reproductive process in proteomic analysis (Figure 5C). Differentially-expressed proteins relating to endocytosis played an important role in vitellogenesis (Figures 5C,E). It is suggested that RH36 may modulate vitellin uptake via endocytosis. Clathrin heavy chain involved in endocytic pathway regulates expression of vitellogenin in R. haemaphysaloides (Kuang et al., 2020). Vitellin is a hemoglycolipoprotein consisting of many bands (James and Oliver, 1997; Yang et al., 2015). Similarly, there were 11 polypeptides of vitellin identified by proteomic analysis in the R. haemaphysaloides ticks and vitellin antibody in this study recognized two bands of 136 and 72 kDa (Supplementary Figure 4) that is different from the peptides of vitellogenin in ovary proteomics analysis (Supplementary Tables 4, 5), resulting in inconsistent phenomenon.
The ovarian transcriptome from post-engorged R. haemaphysaloides females was analyzed and revealed that the transcripts related to lipid metabolism and other metabolisms were differentially expressed. These transcripts may play an important role in ovary development after tick engorgement. The metabolism processes could be involved in the increase of fat body cells observing in the vitellogenesis/oogenesis processes (Denardi et al., 2009). In addition, the abundance of metabolic transcripts could be related to changes in the reproductive tissues during the oogenesis/vitellogenesis processes, such as ovary size growth and oocyte maturation (Moreira et al., 2017), which is consistent with the result that RH36 gene silencing affects tick ovary development and oocytes maturity (Wang et al., 2017).
In the ovary proteomic analysis of ticks fed for 5 days, HSP70 was a typical stress protein in response to micro-injection of RH36 dsRNA (Villar et al., 2010). In the 10 days post-engorged ovary, HSP70 was inhibited significantly after RH36 gene silencing (p < 0.05, Figures 5A, 7) and related closely to protein processing in endoplasmic reticulum (Figure 5E). Recent work has shown that heat shock proteins are secreted to the cell exterior partially in response to stress, we hypothesize that RH36 may modulate ovary development indirectly via mediator HSP70 protein. On the one hand, HSP-70 like protein in salivary gland of ticks participates in variable fibrinogenolysis and accelerates engorged quantity upon blood feeding (Vora et al., 2017), immunomodulatory protein RH36 may be dependent on HSP70 protein to modulate tick blood feeding and nutrient supply (Wang et al., 2017). On the other hand, HSP70, as a endoplasmic reticulum chaperone, exhibits both pro-inflammatory and anti-inflammatory properties, depending on the context in which they encounter responding immune cells (Pockley and Henderson, 2018). HSP70 may have a synergistic effect with RH36 protein in salivary gland of adult ticks. However, the mechanism of HSP70 gene regulating tick oocytes development was still to be further identified. Furthermore, it is still unknown that RH36 gene regulates tick oviposition by accelerating cell death via HSP70.
After HSP70 gene silencing, the mean of 48 h attachment rate, average engorgement rate and expression of vitellin significantly decreased (p < 0.05, Table 2 and Figure 9). This indicates that the down-regulation of HSP70 gene decreases tick blood feeding and nutrition uptake. In ticks, blood feeding triggered salivary gland degeneration (L'Amoreaux et al., 2003) and ovary maturation (Friesen and Kaufman, 2009), suggesting vitellin protein reduced by HSP70 gene silencing inhibits tick ovary development. Starvation activates cell apoptosis and increases tick death. Consistent with this, HSP70 gene silencing dramatically increased tick death rate (Table 2 and Figure 8). We speculated that nutrition deficiency accelerated HSP70-mediated cell damage in ticks (Figure 10). The expression of HSP with antioxidant activity extends lifespan via regulating the downstream protein autophagy-related gene 7 (Morrow et al., 2004; Liao et al., 2008; Chen et al., 2012; Sarup et al., 2014; Vos et al., 2016), while HSP depletion reduces the number of viviparous offspring and simultaneously increases the number of premature nymphs, suggesting an unexpected role in aphid embryogenesis and eclosion (Will et al., 2017). Therefore, it is suggested that RH36 gene regulates vitellogenin uptake via HSP70 protein and extends lifespan of ovary to improve oviposition rate.
Statements
Data availability statement
The Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GIJA00000000.
Ethics statement
The protocols (shvri-ra-2017070578) were approved by the Institutional Animal Care and Use Committee of the Shanghai Veterinary Research Institute, and authorized by the Animal Ethical Committee of Shanghai Veterinary Research Institute.
Author contributions
FW planned the experiments and wrote the paper. FW, YW, GW, and CK performed the experiments. FW, HZ, and JZ analyzed the data. JC and YZ contributed to the reagents or ticks.
Funding
This study was supported by the National Natural Science Foundation of China (NSFC; grant No. 31572512).
Acknowledgments
We would like to thank Beijing Genomics Institute (BGI, Shenzhen, China) for helping with RNA sequencing and Shanghai Applied Protein Technology (APT, shanghai, China) for helping with quantitative proteomic analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00093/full#supplementary-material
References
1
Alarcon-ChaidezF. J.Muller-DobliesU. U.WikelS. (2003). Characterization of a recombinant immunomodulatory protein from the salivary glands of Dermacentor andersoni. Parasite Immunol.25, 69–77. 10.1046/j.1365-3024.2003.00609.x
2
AnatrielloE.RibeiroJ. M.De Miranda-SantosI. K.BrandaoL. G.AndersonJ. M.ValenzuelaJ. G.et al. (2010). An insight into the sialotranscriptome of the brown dog tick, Rhipicephalus sanguineus. BMC Genomics11:450. 10.1186/1471-2164-11-450
3
BarreroR. A.GuerreroF. D.BlackM.MccookeJ.ChapmanB.SchilkeyF.et al. (2017). Gene-enriched draft genome of the cattle tick Rhipicephalus microplus: assembly by the hybrid Pacific Biosciences/Illumina approach enabled analysis of the highly repetitive genome. Int. J. Parasitol. 47, 569–583. 10.1016/j.ijpara.2017.03.007
4
BergmanD. K.PalmerM. J.CaimanoM. J.RadolfJ. D.WikelS. K. (2000). Isolation and molecular cloning of a secreted immunosuppressant protein from Dermacentor andersoni salivary gland. J. Parasitol.86, 516–525. 10.1645/0022-3395(2000)086[0516:IAMCOA]2.0.CO;2
5
ChenS. F.KangM. L.ChenY. C.TangH. W.HuangC. W.LiW. H.et al. (2012). Autophagy-related gene 7 is downstream of heat shock protein 27 in the regulation of eye morphology, polyglutamine toxicity, and lifespan in Drosophila. J. Biomed. Sci.19:52. 10.1186/1423-0127-19-52
6
DanielovaV.HolubovaJ.PejcochM.DanielM. (2002). Potential significance of transovarial transmission in the circulation of tick-borne encephalitis virus. Folia Parasitol.49, 323–325. 10.14411/fp.2002.060
7
DenardiS. E.BecharaG. H.Camargo-MathiasM. I. (2009). Fat body cells of Amblyomma cajennense partially engorged females (Acari: Ixodidae) and their role on vitellogenesis process. Exp. Parasitol.121, 213–218. 10.1016/j.exppara.2008.10.019
8
EhrmannS.RuytsS. C.Scherer-LorenzenM.BauhusJ.BrunetJ.CousinsS.et al. (2018). Habitat properties are key drivers of Borrelia burgdorferi (s.l.) prevalence in Ixodes ricinus populations of deciduous forest fragments. Parasit. Vectors. 11:23. 10.1186/s13071-017-2590-x
9
EisenR. J.EisenL. (2018). The blacklegged tick, Ixodes scapularis: an increasing public health concern. Trends Parasitol.34, 295–309. 10.1016/j.pt.2017.12.006
10
Esmeralda Parra-PeralboJ. C. (2011). Drosophila lipophorin receptors mediate the uptake of neutral lipids in oocytes and imaginal disc cells by an endocytosis-independent mechanism. PLoS Genet.7:e1001297. 10.1371/journal.pgen.1001297
11
FriesenK. J.KaufmanW. R. (2009). 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.55, 936–942. 10.1016/j.jinsphys.2009.06.007
12
GalayR. L.AungK. M.Umemiya-ShirafujiR.MaedaH.MatsuoT.KawaguchiH.et al. (2013). Multiple ferritins are vital to successful blood feeding and reproduction of the hard tick Haemaphysalis longicornis. J. Exp. Biol.216, 1905–1915. 10.1242/jeb.081240
13
GrabowskiJ. M.TsetsarkinK. A.LongD.ScottD. P.RosenkeR.SchwanT. G.et al. (2017). Flavivirus infection of Ixodes scapularis (black-legged tick) ex vivo organotypic cultures and applications for disease control. MBio8, e01255–e01217. 10.1128/mBio.01255-17
14
Gulia-NussM.NussA. B.MeyerJ. M.SonenshineD. E.RoeR. M.WaterhouseR. M.et al. (2016). Genomic insights into the Ixodes scapularis tick vector of Lyme disease. Nat. Commun.7, 1569–1576. 10.1038/ncomms10507
15
HajdusekO.SojkaD.KopacekP.BuresovaV.FrantaZ.SaumanI.et al. (2009). Knockdown of proteins involved in iron metabolism limits tick reproduction and development. Proc. Natl. Acad. Sci. U.S.A. 106, 1033–1038. 10.1073/pnas.0807961106
16
HanR.YangJ.NiuQ.LiuZ.ChenZ.KanW.et al. (2018). Molecular prevalence of spotted fever group rickettsiae in ticks from Qinghai province, northwestern China. Infect. Genet. Evol. 57, 1–7. 10.1016/j.meegid.2017.10.025
17
HarrisE. K.VerhoeveV. I.BanajeeK. H.MacalusoJ. A.AzadA. F.MacalusoK. R. (2017). Comparative vertical transmission of Rickettsia by Dermacentor variabilis and Amblyomma maculatum. Ticks Tick Borne Dis. 8, 598–604. 10.1016/j.ttbdis.2017.04.003
18
JamesA. M.OliverJ. H.Jr. (1997). Purification and partial characterization of vitellin from the black-legged tick, Ixodes scapularis. Insect Biochem. Mol. Biol. 27, 639–649. 10.1016/S0965-1748(97)00038-6
19
JongejanF.SuB. L.YangH. J.BergerL.BeversJ.LiuP. C.et al. (2018). Molecular evidence for the transovarial passage of Babesia gibsoni in Haemaphysalis hystricis (Acari: Ixodidae) ticks from Taiwan: a novel vector for canine babesiosis. Parasit. Vectors. 11:134. 10.1186/s13071-018-2722-y
20
KaufmanW. R. (2004). Assuring paternity in a promiscuous world: are there lessons for ticks among the insects?Parasitology129, 145–160. 10.1017/S0031182004004846
21
KhalilS. M.DonohueK. V.ThompsonD. M.JeffersL. A.AnanthapadmanabanU.SonenshineD. E.et al. (2011). Full-length sequence, regulation and developmental studies of a second vitellogenin gene from the American dog tick, Dermacentor variabilis. J. Insect Physiol. 57, 400–408. 10.1016/j.jinsphys.2010.12.008
22
KonnaiS.NakajimaC.ImamuraS.YamadaS.NishikadoH.KodamaM.et al. (2009). Suppression of cell proliferation and cytokine expression by HL-p36, a tick salivary gland-derived protein of Haemaphysalis longicornis. Immunology126, 209–219. 10.1111/j.1365-2567.2008.02890.x
23
KotalJ.LanghansovaH.LieskovskaJ.AndersenJ. F.FrancischettiI. M.ChavakisT.et al. (2015). Modulation of host immunity by tick saliva. J. Proteomics128, 58–68. 10.1016/j.jprot.2015.07.005
24
KuangC.WangF.ZhouY.CaoJ.ZhangH.GongH.et al. (2020). Molecular characterization of clathrin heavy chain (Chc) in Rhipicephalus haemaphysaloides and its effect on vitellogenin (Vg) expression via the clathrin-mediated endocytic pathway. Exp. Appl. Acarol.80, 71–89. 10.1007/s10493-019-00438-5
25
L'AmoreauxW. J.JunaidL.TrevidiS. (2003). Morphological evidence that salivary gland degeneration in the American dog tick, Dermacentor variabilis (Say), involves programmed cell death. Tissue Cell35, 95–99. 10.1016/S0040-8166(02)00109-X
26
LiaoP. C.LinH. Y.YuhC. H.YuL. K.WangH. D. (2008). The effect of neuronal expression of heat shock proteins 26 and 27 on lifespan, neurodegeneration, and apoptosis in Drosophila. Biochem. Biophys. Res. Commun.376, 637–641. 10.1016/j.bbrc.2008.08.161
27
MitchellR. DIIISonenshineD. E.Perez De LeonA. A. (2019). Vitellogenin receptor as a target for tick control: a mini-review. Front. Physiol.10:618. 10.3389/fphys.2019.00618
28
MooreT. C.PulscherL. A.CaddellL.Von FrickenM. E.AndersonB. D.GonchigooB.et al. (2018). Evidence for transovarial transmission of tick-borne rickettsiae circulating in Northern Mongolia. PLoS Negl. Trop. Dis. 12:e0006696. 10.1371/journal.pntd.0006696
29
MoreiraH. N. S.BarcelosR. M.VidigalP. M. P.KleinR. C.MontandonC. E.MacielT. E. F.et al. (2017). A deep insight into the whole transcriptome of midguts, ovaries and salivary glands of the Amblyomma sculptum tick. Parasitol. Int.66, 64–73. 10.1016/j.parint.2016.10.011
30
MorrowG.BattistiniS.ZhangP.TanguayR. M. (2004). Decreased lifespan in the absence of expression of the mitochondrial small heat shock protein Hsp22 in Drosophila. J. Biol. Chem. 279, 43382–43385. 10.1074/jbc.C400357200
31
NijhofA. M.BalkJ. A.PostigoM.JongejanF. (2009). Selection of reference genes for quantitative REAL-TIME PCR studies in Rhipicephalus (Boophilus) microplus and Rhipicephalus appendiculatus ticks and determination of the expression profile of Bm86. BMC Mol. Biol.10:112. 10.1186/1471-2199-10-112
32
OyugiM. O.KinyuaJ. K.MagiriE. N.KigoniM. W.CostaE. P.GithakaN. W. (2018). In silico characterization and structural modeling of Dermacentor andersoni p36 immunosuppressive protein. Adv. Bioinformatics. 2018:7963401. 10.1155/2018/7963401
33
ParolaP.RaoultD. (2001). Ticks and tickborne bacterial diseases in humans: an emerging infectious threat. Clin. Infect. Dis. 32, 897–928. 10.1086/319347
34
PernerJ.SobotkaR.SimaR.KonvickovaJ.SojkaD.OliveiraP. L.et al. (2016). Acquisition of exogenous haem is essential for tick reproduction. Elife. 5:e12318. 10.7554/eLife.12318
35
PockleyA. G.HendersonB. (2018). Extracellular cell stress (heat shock) proteins-immune responses and disease: an overview. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 373:20160522. 10.1098/rstb.2016.0522
36
RosellR.CoonsL. B. (1992). The role of the fat body, midgut and ovary in vitellogenin production and vitellogenesis in the female tick, Dermacentor variabilis. Int. J. Parasitol. 22, 341–349. 10.1016/S0020-7519(05)80012-8
37
SarupP.SorensenP.LoeschckeV. (2014). The long-term effects of a life-prolonging heat treatment on the Drosophila melanogaster transcriptome suggest that heat shock proteins extend lifespan. Exp. Gerontol. 50, 34–39. 10.1016/j.exger.2013.11.017
38
SchneiderC. A.RasbandW. S.EliceiriK. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods. 9, 671–675. 10.1038/nmeth.2089
39
ThompsonD. M.KhalilS. M.JeffersL. A.AnanthapadmanabanU.SonenshineD. E.MitchellR. D.et al. (2005). In vivo role of 20-hydroxyecdysone in the regulation of the vitellogenin mRNA and egg development in the American dog tick, Dermacentor variabilis (Say). J. Insect Physiol.51, 1105–1116. 10.1016/j.jinsphys.2005.05.011
40
Umemiya-ShirafujiR.BoldbaatarD.LiaoM.BatturB.RahmanM. M.KubokiT.et al. (2012). Target of rapamycin (TOR) controls vitellogenesis via activation of the S6 kinase in the fat body of the tick, Haemaphysalis longicornis. Int. J. Parasitol. 42, 991–998. 10.1016/j.ijpara.2012.08.002
41
Umemiya-ShirafujiR.HattaT.OkuboK.SatoM.MaedaH.KumeA.et al. (2017). Transovarial persistence of Babesia ovata DNA in a hard tick, Haemaphysalis longicornis, in a semi-artificial mouse skin membrane feeding system. Acta Parasitol. 62, 836–841. 10.1515/ap-2017-0100
42
VillarM.AyllonN.BusbyA. T.GalindoR. C.BlouinE. F.KocanK. M.et al. (2010). Expression of heat shock and other stress response proteins in ticks and cultured tick cells in response to Anaplasma spp. infection and heat shock. Int. J. Proteomics2010:657261. 10.1155/2010/657261
43
VoraA.TaankV.DuttaS. M.AndersonJ. F.FishD.SonenshineD. E.et al. (2017). Ticks elicit variable fibrinogenolytic activities upon feeding on hosts with different immune backgrounds. Sci. Rep. 7:44593. 10.1038/srep44593
44
VosM. J.CarraS.KanonB.BosveldF.KlaukeK.SibonO. C.et al. (2016). Specific protein homeostatic functions of small heat-shock proteins increase lifespan. Aging Cell15, 217–226. 10.1111/acel.12422
45
WangF.LuX.GuoF.GongH.ZhangH.ZhouY.et al. (2017). The immunomodulatory protein RH36 is relating to blood-feeding success and oviposition in hard ticks. Vet. Parasitol. 240, 49–59. 10.1016/j.vetpar.2017.03.017
46
WillT.SchmidtbergH.SkaljacM.VilcinskasA. (2017). Heat shock protein 83 plays pleiotropic roles in embryogenesis, longevity, and fecundity of the pea aphid Acyrthosiphon pisum. Dev. Genes Evol. 227, 1–9. 10.1007/s00427-016-0564-1
47
XavierM. A.TirloniL.PintoA. F. M.DiedrichJ. K.YatesJ. R.MulengaA.et al. (2018). A proteomic insight into vitellogenesis during tick ovary maturation. Sci. Rep.8:4698. 10.1038/s41598-018-23090-2
48
YangX.YuZ.HeY.XuX.GaoZ.WangH.et al. (2015). Purification of vitellin and dynamics of vitellogenesis in the parthenogenetic tick Haemaphysalis longicornis (Acari: Ixodidae). Exp. Appl. Acarol.65, 377–388. 10.1007/s10493-014-9866-z
49
YuX.GongH.ZhouY.ZhangH.CaoJ.ZhouJ. (2015). Differential sialotranscriptomes of unfed and fed Rhipicephalus haemaphysaloides, with particular regard to differentially expressed genes of cysteine proteases. Parasit. Vectors8:597. 10.1186/s13071-015-1213-7
50
YuY.CaoJ.ZhouY.ZhangH.ZhouJ. (2013). Isolation and characterization of two novel serpins from the tick Rhipicephalus haemaphysaloides. Ticks Tick Borne Dis.4, 297–303.
Summary
Keywords
Rhipicephalus haemaphysaloides, RH36, vitellogenesis, RNA interference, HSP70 protein
Citation
Wang F, Wang Y, Wang G, Zhang H, Kuang C, Zhou Y, Cao J and Zhou J (2020) Ovary Proteome Analysis Reveals RH36 Regulates Reproduction via Vitellin Uptake Mediated by HSP70 Protein in Hard Ticks. Front. Cell. Infect. Microbiol. 10:93. doi: 10.3389/fcimb.2020.00093
Received
02 November 2019
Accepted
21 February 2020
Published
10 March 2020
Volume
10 - 2020
Edited by
Lucas Tirloni, National Institute of Allergy and Infectious Diseases (NIAID), United States
Reviewed by
Camila Fernanda Waltero Orjuela, Federal University of Rio de Janeiro, Brazil; Marina Amaral Xavier, Federal University of Rio Grande Do Sul, Brazil
Updates
Copyright
© 2020 Wang, Wang, Wang, Zhang, Kuang, Zhou, Cao and Zhou.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinlin Zhou jinlinzhou@shvri.ac.cn
This article was submitted to Parasite and Host, a section of the journal Frontiers in Cellular and Infection Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.