Identification of the Bcl-2 and Bax Homologs From the Tick Rhipicephalus Haemaphysaloides and Their Function in the Degeneration of Tick Salivary Glands

doi:10.21203/rs.3.rs-154051/v1

Download PDF

Research

Identification of the Bcl-2 and Bax Homologs From the Tick Rhipicephalus Haemaphysaloides and Their Function in the Degeneration of Tick Salivary Glands

https://doi.org/10.21203/rs.3.rs-154051/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 04 Aug, 2021

Read the published version in Parasites & Vectors →

You are reading this latest preprint version

Background

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

Methods

The full-length cDNA of the RhBcl-2 and RhBax genes were obtained by transcriptome analysis and RhBcl-2 and RhBax were expressed in E. coli. Two molecules containing conserved BH (Bcl-2 family homology) domains were identified and named RhBcl-2 and RhBax. After protein purification and mouse immunizations, specific polyclonal antibodies (PcAb) were created in response to the recombinant protein. Reverse transcription quantitative PCR (RT-qPCR) and western blot were used to detect the existence of RhBcl-2 and RhBax in ticks. TUNEL assays were used to determine the apoptosis level in the salivary glands at different feeding times after gene silencing. The interactions between RhBcl-2 and RhBax were identified by co-transfection and GST pull-down assays. The RT-qPCR assay demonstrated that the transcription of RhBax genes increased significantly during the engorged periods of the tick developmental stages (engorged larval, nymph, and adult female ticks).

Results

Transcriptional levels of RhBcl-2 and RhBax in the salivary glands increased more significantly than other tissues post engorgement. RhBcl-2 treatment significantly restrained tick feeding. In contrast, RhBax interference had no effect on tick feeding. TUNEL staining showed that apoptosis levels were significantly reduced after interfering with RhBcl-2 expression. Co-transfection and GST pull-down assays showed that RhBcl-2 and RhBax could combine with each other, but failed to combine without the BH3 domain.

Conclusions

This study identified the roles of RhBcl-2 and RhBax in tick salivary gland degeneration and clarifies their interactions.

Parasitology

Salivary gland degeneration

Apoptosis

Bcl-2

Bax

Ticks are obligatory blood-feeding arthropods and vectors of many pathogens [1]. Ticks usually stay attached to their hosts and feed for several days or even weeks [2]. The salivary glands of ticks mediate diverse functions that ensure tick survival [3]. Tick-borne pathogens (TBP) are transmitted to the host by the saliva used during tick feeding [4].

The salivary glands of female ticks degenerate rapidly within four days after engorgement [5]. Salivary gland degeneration begins at the early stage of rapid feeding and may be caused by apoptosis and autophagy related to the increase of ecdysone in the hemolymph [6, 7]. During degeneration, the granular acini undergo DNA fragmentation and caspase enzyme activity increases [8, 9]. According to the comprehensive analysis of the global protein changes, some proteins associated with apoptosis and autophagy changed and some proteins linked to the degradation of DNA and proteins were consistently upregulated [10]. However, the mechanism of these effects is unknown.

Apoptosis is an important part of programmed cell death (PCD). Apoptosis is programmed to transmit death induction signals and carry out biochemical processes of cell decomposition by a group of specific interacting proteins [11]. For embryonic development, tissue homeostasis, and organism stability in multicellular organisms, apoptosis is indispensable and is a normal physiological process [12]. During apoptosis, a decrease in cell volume, cytoplasmic condensation, mitochondrial membrane permeabilization, DNA fragmentation, and chromatin condensation followed by nuclear fragmentation and cytoplasmic membrane blebbing take place, which leads to complete cell division. In the late stage of apoptosis, apoptotic bodies are produced due to cell fragmentation [13].

Apoptosis is caused by sequential activation of cysteine proteases of the caspase family in two distinct, but converging pathways [14]. The extrinsic pathway, also known as the death receptor pathway, activates Caspase-8 when the “death receptor” Fas binds to its ligand FasL, and then recruits FADD proteins and caspase-8 enzymes in the cytoplasm to form a death complex DISC [15]. The intrinsic pathway (also termed the mitochondrial pathway) can be triggered by the pro-apoptotic proteins of the Bcl-2 family. These cause mitochondrial membrane permeabilization to release mitochondrial cytochrome c into the cytoplasm, further activating caspases to induce apoptosis [16]. These initiator caspases can cleave and activate the effector caspases (Caspase-3, Caspase-6, and Caspase-7) that mediate cellular demolition by cleaving multiple critical cellular proteins [17].

The Bcl-2 (B-cell lymphoma-2) protein family is involved in apoptosis via the intrinsic pathway. Bcl-2 was first discovered via the chromosome translocation that hallmarks human follicular lymphoma, and it was subsequently found to promote cell survival as an inhibitory apoptosis protein of the Bcl-2 protein family [18]. Its homologs Bcl-X, Bcl-w, Mcl-1, and A1 were also found to inhibit apoptosis, but Bax, Bak, and Bok promoted apoptosis. BH (Bcl-2 homology) domain is a vital functional domain for the interaction of the Bcl-2 family proteins. The BH3 domain was a key for Bax, Bak, and Bok to promote apoptosis [19]. Bim, Bid, Bad, Puma, and Noxa only had the BH3 domain, so they were termed only-BH3 proteins.

The Bcl-2 family proteins regulate apoptosis mainly by regulating mitochondrial membrane permeability, working by themselves or in combination with other proteins. For example, Bax forms oligomers on the mitochondrial membrane that perturb membrane integrity and promote apoptosis, while Bcl-2 and Bcl-xL inhibit their effects by binding to Bax to form heterodimers [20, 21]. BH3-only protein can play a direct or indirect pro-apoptotic function. It can directly bind and activate Bax and Bak proteins to promote apoptosis, and it can indirectly promote the expression of Bax and Bak proteins to cause apoptosis by inhibiting the action of the anti-proteins [22, 23].

Bcl-2 family proteins exist in both mammalian and non-mammalian species. The nematode Caenorhabditis elegans has two members of the Bcl-2 family, EGL-1 and CED-9, that regulate cell death [24]. Two Bcl-2 family proteins which have high homology with Bok found in Drosophila, Debcl and Buffy [6]. Using RNA interference, Ayllon [25] injected Bcl-2 dsRNA into Ixodes scapularis and showed feeding disruption. This indicated that Bcl-2 plays an important role in the intrinsic apoptosis of ticks. However, the mechanism by which Bcl-2 family proteins regulate apoptosis in ticks is unknown. In this study, we identified Bcl-2 and Bax from the tick Rhipicephalus haemaphysaloides and investigated their functions in the degeneration of tick salivary glands at different feeding times.

2.1 Tick Feeding and tissue collection

Adult Rhipicephalus haemaphysaloides Supino 1897 (Ixodida: Ixodidae) were collected from Wuhan in Hubei Province, China. The ticks were fed on the ears of New Zealand white rabbits (SLAC, Shanghai Institutes for Biological Science, CAS), and maintained in artificial climate incubators of the Shanghai Veterinary Research Institute [26]. Rabbits and mice were maintained following the approved guidelines of the Institutional Animal Care and Use Committee of the Shanghai Veterinary Research Institute (IACUC approval number: SHVRI-SZ-20190726-02; SHVRI-SZ-20191024-01). Tick tissues were rapidly dissected, washed with phosphate-buffered saline (PBS, PH 7.4, with 0.14 M NaCl and 0.0027 M KCl, 0.01 M phosphate buffer; Gibco, Life Technologies), and placed in PBS or TRIzol (Invitrogen, Carlsbad, CA, USA) reagent at − 80°C.

2.2 RNA extraction and cDNA synthesis

RNA was derived from ticks dissected at different development stages (eggs, unfed and fed larvae, nymphs, and adults) and different tissues (salivary glands, midguts, and ovaries of female ticks), and preserved in TRIzol reagent. The HiScript® III RT SuperMix for qPCR (+ gDNA wiper) kit (Vazyme Biotech, Nanjing, China) was used, according to manufacturer protocol, to synthesis cDNA.

2.3 Cloning, sequence analysis, and expression of RhBcl-2 and RhBax

RhBcl-2 and RhBax primers were designed based on a comparison of the salivary gland transcriptomes of starved and engorged R. haemaphysaloides [27, 28]. A BLAST analysis of the translation products deduced from the open reading frames (ORFs) was performed. SignalP 4.1 (http://www.cbs.dtu.dk/services/SignalP/) [29] and ExPASy (http://web.expasy.org/compute pi/) [30] were used for the Signal peptide analysis and isoelectric point (PI) prediction. We aligned RhBcl-2 and RhBax with the Bcl-2 and Bax protein sequences of other species using Genetyx ver. 6 (GENETYX, Tokyo, Japan). For phylogenetic analysis, the alignment of the sequences was performed using the MUSCLE algorithm [31] and inferred using the Maximum likelihood method with the default settings in MEGA X software [32]. Bootstrap support values were estimated using 500 bootstrap replicates [28].

2.4 RT-qPCR analyses

The expression levels of RhBcl-2 and RhBax were examined in ticks from different development stages (eggs, unfed and fed larvae, nymphs, and adults) and different tissues (salivary glands, midguts, and ovaries of female ticks). After RhBcl-2 and RhBax primers (Table S1) designed using Primer Premier 5, the cDNAs of the above stages and tissues were analyzed by RT-qPCRs which were conducted using ChamQ Universal SYBR qPCR Master Mix (Vazyme) green and gene-specific (Table S1) primers with a QuantStudio 5 PCR System (Applied Biosystems, Austin, TX, USA). The RT-qPCR process was 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s, followed by analysis of the melting curve. All samples were analyzed in triplicate. The data used elongation factor-1(ELF1A, Genbank accession no. AB836665) as an internal control [31], and this was used to analyze the relative gene expression in each sample by the 2^−△Ct method [32, 33].

2.5 Expression of recombinant RhBcl-2 or RhBax and antibody production

Specific RhBcl-2 and RhBax primers (Table S2) were designed in a pET-30a or PGEX-4T-1 vector. The RhBcl-2 and RhBax amplified-PCR products were purified and digested with BamHI, EcoRI, or XhoI (New England Biolabs, USA) and ligated into pET-28a or PGEX-4T-1(Invitrogen) using In-Fusion HD Cloning Kits (Takara Clontech, Mountain View, CA, USA) [28]. These recombinants were transformed and expressed in Escherichia coli (E. coli) BL21 (DE3) strain (TIANGEN, Beijing, China). The strains were grown until the OD₆₀₀ reached 0.8 at 37°C. After isopropyl β-D thiogalacto pyranoside (IPTG) was added to 1 mM, the protein expression was induced at 25°C for 20 h. The recombinant, containing His-tagged protein, was purified by affinity chromatography using Ni-NTA His resin (Thermo Fisher Scientific, Waltham, MA, USA) by gravity. The recombinant containing the GST-tagged protein was purified by affinity chromatography using resin and gravity.

2.6 GST-Pull Down of RhBcl-2 and RhBax

To confirm whether RhBcl-2 interacts with RhBax, 1 mM IPTG induced expression for 8 h of pET-30a-RhBax, PGEX-4T-RhBcl-2, and PGEX-4T-1 recombinant E. coli BL21. The supernatant was separated using GST agarose (Merck, Darmstadt, Germany) according to manufacturer instructions. Anti-GST antibody (Proteintech, 215 Rosemont, IL, USA) and goat anti-mouse secondary antibody (Proteintech) were used to confirm that GST/GST-RhBcl-2 was separated successfully, and anti-His antibody (Proteintech) was used to detect the target protein.

2.7 Cells and Transient Co-transfection assays

HEK 293 cells were maintained in Dulbecco’s Modified Eagle Medium (DMEM, Gibco), supplemented with 8% heat-inactivated fetal bovine serum (Biological Industries, Kibbutz Beit Haemek, Israel) and 1% penicillin (Gibco) at 37°C [28].

The full-length ORF of RhBcl-2 was inserted into the p3×Flag-CMV-14 vector (MiaoLing Plasmid Sharing Platform, Wuhan, China) with Flag tag at the N-terminal with gene specific primers (Table S2). RhBax was inserted into the pcmv-HA vector using same method. The fragment of RhBax deleted BH3 domain (Rhbax- ΔBH3) and BH1 domain (Rhbax-ΔBH1) was amplified by splicing by overlap extension (SOE) of recombinant PCR with specific primers (Table S3), and they were inserted into the pcmv-HA vector [34]. Transfection, using Lipofectamine™ 3000 Transfection Reagent (Invitrogen), was performed according to manufacturer protocol with a DNA to Lipofectamine ratio of 1: 2 w/v. The HEK 293 cells were transformed with 2.5 µg/well of plasmid or co-transfected with the equivalent amount of two different plasmids in 6-well plates [28].

2.8 Antibody generation and western blotting

Recombinant RhBcl-2 and RhBax were expressed and purified according to previously described protocols [35]. Polyclonal antibodies (pcAb) against recombinant RhBcl-2 and RhBax were developed in mice. We tested the specific recognition of the recombinant protein by the antibodies using western blots. We predicted the epitopes of RhBcl-2 and RhBax online (http://www.iedb.org/), and synthesized Keyhole limpet hemocyanin (KLH) [36] coupled polypeptides based on the epitope amino acid sequences (RhBcl-2: GLQWNTCPPLPRPSK, RhBax: STPTHEETREE). Polypeptides were dissolved in PBS and Freund’s adjuvant (complete and incomplete; Invitrogen). Equal volumes of these solutions were then emulsified together and injected into six to eight weeks old BALC/c mice provided by the Shanghai Laboratory Animals Center (Shanghai Institutes for Biological Science, Chinese Academy of Sciences, Shanghai, China). The mixtures were injected three times at two weeks intervals. The sera were collected on the seventh day after the third immunization.

Total proteins of different development stage ticks were extracted using Tris-buffered saline (TBS; 10 mM Tris–HCl, pH 7.5; 150 mM NaCl with 1 mM phenylmethanesulfonyl fluoride). A Bradford Protein Assay Kit (Beyotime Tiangen) was used to determine the protein concentration. Proteins were separated by 10% sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE, 12%, Genescript, Nanjing, China) and transferred onto polyvinylidene difluoride (PDVF) membranes. After blocking, the proteins were incubated at 4°C overnight with primary antibodies. The sera anti RhBcl-2 and RhBax were used to detect RhBcl-2 and RhBax in protein extracts, and anti-tubulin (1:1000, Proteintech, Rosemont, IL, USA) primary antibody was used as a constitutive control to normalize the signal from the target protein. The membranes were rinsed five times in tris-buffered saline with Tween-20 (TBST) and incubated with a secondary antibody (goat anti-mouse, 1:5000, Cell Signaling Technology, USA) for 1 h at room temperature. The membranes were rinsed five times by TBST. Images were captured by ChemiDoc Touch (Bio-rad, Hercules, CA, USA) [28].

2.9 RNAi of RhBcl-2 or RhBax

The RNAi experiments were designed against RhBcl-2 and RhBax genes. For the design of RNAi primers, RhBcl-2 and RhBax sequences were screened by Primer Premier 5. Specific primers (Table S4) containing the T7 polymerase promoter sequence were used for PCR amplification. The amplicons were then purified to obtain templates for double-stranded RNA synthesis using the T7 RiboMAX Express RNAi system (Promega, Madison, WI, USA). The unrelated gene Luciferase dsRNA was synthesized using the same methods described previously and used as the negative control. Unfed female ticks (N = 20 females per group) were microinjected with approximately 1 µg of dsRNA. Control ticks were injected with unrelated Luciferase dsRNA. After dsRNA injection, female ticks were held in a humidity chamber for 1 day after which they were allowed to feed on rabbit ears. The biological parameters analyzed were attachment rate at 24 h, number of engorged ticks, death rate, and weight. RNA was used to characterize gene knockdown by real-time RT-PCR.

We interfered with the salivary gland using dsRNA in vitro. The salivary gland of female ticks fed on rabbit for 5 days was dissected and placed into complete L15 medium with 1% penicillin-streptomycin [37]. Each well received 5 µg of dsRNA and incubated at 27°C with no CO₂ for 48 h. The control group received unrelated Luciferase dsRNA. Each treatment group included the salivary glands from 10 ticks. After incubation, a pair of the salivary glands in each group was fixed in 4% paraformaldehyde for TUNEL staining or fixed in 2.5% glutaraldehyde for at least 24 h at 4°C for TEM. The remaining salivary glands were used for western blots and RNA.

2.10 TUNEL staining

Salivary glands were fixed in 4% formalin and paraffin-embedded. Sections of salivary glands were mounted on microscope slides. Tissue sections were then deparaffinized, washed in 100% ethanol, and rehydrated. Samples were washed with PBS. After antigen retrieval with 0.1% Triton X-100, the tissues were incubated for 1 h with 1:9 TdT mixed with fluorescent-labeled dUTP at 37°C, following the instructions of the Roche In Situ Cell Death Detection Kit, POD (Roche, Mannheim, Germany). After washing the Sect. 2–3 times with PBS, the sections were stained with 1 µg/mL 4’, 6’-diamidino-2-phenylindole (DAPI, Invitrogen) in distilled H₂O for 20 min [28]. After washing, the sections were mounted using Lab Vision™ PermaFluor™ (Invitrogen) medium under glass coverslips, then viewed and photographed on a Pannoramic DESK Digital Slide Scanner (3D HISTECH, Budapest, Hungary) [28].

2.11 Data analysis

All statistical analyses were performed using GraphPad Prism 6.0 software (Graph Pad Software Inc., San Diego, CA, USA). Mean ± standard error (SEM) values were calculated for three separate experiments, and two-tailed Student’s t tests were used to identify significant differences between groups. P < 0.05 was considered to be statistically significant.

3.1 Sequencing analysis of RhBcl-2 and RhBax

We used nucleic acid sequences obtained from RNAsEq. Specific cloning primers were designed according to the predicted sequences of R. haemaphysaloides Bcl-2 and Bax (Table S2). The ORF regions of Bcl-2 and Bax genes were cloned from the cDNA of fully engorged R. haemaphysaloides salivary glands and named RhBcl-2 and RhBax. RhBcl-2 ORF has 576 bp and encodes a protein of 192 amino acid residues with a deduced molecular weight (MW) and theoretical isoelectric point (PI) of 21.4 kDa and 4.87, respectively. RhBax ORF has 648 bp and encodes a protein with 216 amino acid residues. It has a deduced MW of 24.6 kDa and a PI of 5.28. RhBcl-2 had the conserved BH1, BH2, BH3, and BH4 domains, but RhBax had the conserved BH1 and BH3 domains.

3.2 Transcription and translation of RhBcl-2 and RhBax profiles in different stages, tissues, and feeding periods

The cDNA of eggs, larvae (unfed and engorged), and nymphs (unfed and engorged) were subjected to RT-qPCR to evaluate the expression profiles of RhBcl-2 and RhBax genes during the developmental stages. The cDNA of unfed adults (females and males), fed adults (females and males), and engorged female ticks were used to determine the sex-specific profile of RhBcl-2 and RhBax genes. The transcription level of RhBax increased after adult female feeding.

After microdissection, the cDNA of salivary glands, ovary, and midgut of adult female ticks at different feeding times were analyzed for RhBcl-2 and RhBax gene expression profiles. RhBcl-2 and RhBax genes were expressed in different tissues during tick feeding (including the early-feeding period, the fast-feeding period, and the end of the feeding period). However, RhBcl-2 and RhBax had the same expression trends in different organs.

Western blot analysis used anti-sera of RhBcl-2 and RhBax. RhBax was detected during the fast-feeding time (fed 5 d to 7 d) to post engorgement and RhBcl-2 was detected at all stages.

3.3 Expression of RhBcl-2 and RhBax

The coding sequences of RhBcl-2 and RhBax were cloned into prokaryotic expression vectors (pET-28a and pGEX-4T-1) to produce recombinant RhBcl-2 and RhBax. All of the recombinant proteins were expressed as inclusion bodies in E. coli.

After solubilization and purification, His-RhBax (molecular weight: ~38 kDa) and GST-RhBcl-2 (molecular weight: ~50 kDa) were obtained, respectively. The MWs of RhBcl-2 were likely the corresponding mammalian Bcl-2 molecules. Specific polyclonal antibodies (PcAb) were created in response to the recombinant protein. Western blot analysis revealed that the two sera were able to identify the full recombinant proteins, and the fragment sizes were consistent with those observed by SDS-PAGE.

3.4 RhBcl-2 binds with RhBax by BH3 domain

To verify that RhBcl-2 interacts with RhBax, an in-vitro GST-Pull down assay was used to confirm the interaction between His-RhBax and GST-RhBcl-2. For this experiment, two recombinant proteins with either a GST-tag or His-tag, respectively (His-RhBax and GST-RhBcl-2), were produced from E. coli BL21. Immunoblotting showed that the GST and GST-RhBcl-2 were successfully pulled down with GST beads and that the GST-RhBcl-2 interacted directly with His-RhBax. GST as a negative control confirmed that the interaction of GST-RhBcl-2 and His-RhBax was specific.

Co-transfection of HEK 293 cells with plasmids containing RhBcl-2 and RhBax genes resulted in co-expression of the proteins. The results also showed that RhBcl-2 and RhBax could interact with each other, but failed to combine without the BH3 domain (Fig. 8).

3.5 RNAi assay of RhBcl-2 and RhBax

To identify the function of RhBcl-2 and RhBax, the genes of RhBcl-2 and RhBax were in vivo targets of RNA interference. The RT-qPCR analysis of ticks showed decreased levels of RhBcl-2 and RhBax in RhBcl-2-dsRNA and RhBax-dsRNA injected groups compared to luciferase-dsRNA injected controls (RhRhBcl-2: t₍₃₎ = 3.967, P = 0.0166; RhBax: t₍₃₎ = 4.153, P = 0.0142). Compared to the control groups, ticks injected with RhBax-dsRNA sere similar to the control groups, but ticks injected with RhBcl-2-dsRNA were significantly different. However, there was no significant difference in the attachment rate between the experiment and control groups.

TUNEL assays showed that salivary gland interference by RhBcl-2 dsRNA increased apoptosis levels, but RhBax was similar to the control.

Discussion

R. haemaphysaloides is a three-host hard tick widely distributed in China. It is a carrier of several human pathogens [38, 39, 40]. The salivary glands of female ticks undergo degeneration after engorgement[3]. Apoptosis plays a role in the process while caspase genes and apoptosis-related genes are upregulated [41]. We identified and characterized two Bcl-2 family proteins (Bcl-2 and Bax) in R. haemaphysaloides and found that they are involved in salivary gland degeneration.

The Bcl-2 family proteins are also found in insects such as Drosophila, Bombyx mori, and Anopheline mosquitoes, and are involved in growth and metamorphosis [24, 42, 43, 44]. Ayllon [25] found that Bcl-2 may play a key role in the intrinsic apoptosis of Ixodes scapularis, but other Bcl-2 family proteins and the apoptosis mechanism are unknown. We found that RhBcl-2 and RhBax are involved in salivary glands apoptosis degeneration in R. haemaphysaloides. Sequence analysis demonstrated that RhBcl-2 and RhBax have conserved sequences characteristic of BH domains. The BH3 domain was a key for Bax, Bak, and Bok to promote apoptosis and combine anti-apoptosis proteins with pro-apoptosis proteins [45, 46]. A GST pull-down assay verified that RhBcl-2 and RhBax combined with each other. Co-transfection assays also showed that RhBcl-2 and RhBax could interact with each other but did not combine without the BH3 domain.

Salivary gland degeneration is triggered after the completion of female tick feeding [47]. The RT-qPCR results showed that RhBcl-2 and RhBax had a high level of transcription in salivary glands during tick engorgement as determined by western blot. These data suggest that the roles of RhBcl-2 and RhBax are related to apoptosis, and degeneration, of the salivary glands.

RNAi was used to study the physiological roles of RhBcl-2 and RhBax. TUNEL staining was used to evaluate the rate of DNA fragmentation in degenerated tick salivary glands. Compared to the control group, the TUNEL positive staining rate increased significantly by interfering with RhBcl-2. These results showed that RhBcl-2 functions in preventing apoptosis during degeneration of the salivary glands. Bax is a pro-apoptosis protein in the Bcl-2 protein family [48]. However, interfering with RhBax had no significant effect compared to the control group. This suggests that other apoptosis pathways help regulate the apoptosis of the salivary glands.

RhBcl-2 silencing produced the most significant phenotype alteration in ticks during blood feeding compared to RhBax and the control groups. After injection with dsRhBcl-2, the female ticks were unable to feed and engorge until they died. These data illustrated the importance of RhBcl-2 during blood feeding. Bcl-2 inhibited Bax to protect cell survival, but had other affects [49, 50]. Consequently, RhBcl-2 functioned to protect cell survival, promote normal organ functioning, and ensure that ticks grew in the stage of blood feeding. The function of the RhBax protein remains unclear.

In conclusion, two Bcl-2 family molecules were identified in ticks and the roles of RhBcl-2 and RhBax in tick physiology were studied. These findings allowed for advanced studies of tick apoptosis. RhBcl-2 and RhBax, like Bcl-2 and Bax in mammals, combine with each other to regulate apoptosis in the tick salivary glands. The knockdown of RhBcl-2 in vivo inhibited tick growth and blood feeding and demonstrated that RhBcl-2 has the potential to be a tick vaccine protein.

TBPs, Tick-borne pathogens; PCD, Programmed cell death; PcAb, polyclonal antibodies; RT-qPCR, reverse transcription quantitative PCR; PBS, Phosphate-buffered saline; TBS; 10 mM Tris–HCl, pH 7.5; 150 mM NaCl with 1 mM phenylmethanesulfonyl fluoride; SDS-PAGE, Sodium dodecyl sulfate polyacrylamide gel electrophoresis; IPTG, Isopropyl β-D-1-Thiogalactopyranoside; DAPI, 4', 6'-diamidino-2-phenylindole; TUNEL, Terminal deoxynucleotidyl transferase dUTP Nick-End Labeling; BH, Bcl-2 homology

Ethics approval and consent to participate

The treatment of rabbits and mice was approved by the Institutional Animal Care and Use Committee of the Shanghai Veterinary Research Institute (IACUC approval number: SHVRI-SZ-20190726-02; SHVRI-SZ-20191024-01), and authorized by the Animal Ethical Committee of Shanghai Veterinary Research Institute.

Consent for publication

Not applicable.

Availability of data and material

The datasets supporting the conclusions of this article are uploaded as Additional files.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was supported by a grant from “Central Public-interest Scientific Institution Basal Research Fund (No. Y2021GH01-3)” and Strategic International Collaborative Research Project (JPJ008837) Promoted by the Ministry of Agriculture, Forestry and Fisheries, Japan.

Authors' contributions

SH conceived the study, performed the statistical analyses and participated in the dynamic detection of the target genes. SH drafted the manuscript. YW and ZX participated in the design of the study and the interpretation of the data. YZ, HZ and JC participated in the experiments on animals. JZ participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgments

We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.

Jongejan F, Uilenberg G. The global importance of ticks. Parasitology. 2004;129:S3-S14;
Sonenshine DE. Biology of ticks. New York: Oxford University Press; 1991.
Simo L, Kazimirova M, Richardson J, Bonnet SI. The Essential Role of Tick Salivary Glands and Saliva in Tick Feeding and Pathogen Transmission. Front Cell Infect Microbiol. 2017;7:281;
Kazimirova M, Stibraniova I. Tick salivary compounds: their role in modulation of host defences and pathogen transmission. Front Cell Infect Microbiol. 2013;3:43;
Francischetti IM, Sa-Nunes A, Mans BJ, Santos IM, Ribeiro JM. The role of saliva in tick feeding. Front Biosci (Landmark Ed). 2009;14:2051-88;
Abdelwahid E, Rolland S, Teng X, Conradt B, Hardwick JM, White K. Mitochondrial involvement in cell death of non-mammalian eukaryotes. Biochim Biophys Acta. 2011;1813 4:597-607;
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 1:33-42;
L'Amoreaux WJ, Junaid L, Trevidi S. Morphological evidence that salivary gland degeneration in the American dog tick, Dermacentor variabilis (Say), involves programmed cell death. Tissue Cell. 2003;35 2:95-9;
Freitas DR, Rosa RM, Moura DJ, Seitz AL, Colodel EM, Driemeier D, et al. Cell death during preoviposition period in Boophilus microplus tick. Vet Parasitol. 2007;144 3-4:321-7;
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. 2018;9:1943;
Maghsoudi N, Zakeri Z, Lockshin RA. Programmed cell death and apoptosis--where it came from and where it is going: from Elie Metchnikoff to the control of caspases. Exp Oncol. 2012;34 3:146-52;
Kerr JF, Wyllie AH, Currie AR. Apoptosis: a basic biological phenomenon with wide-ranging implications in tissue kinetics. Br J Cancer. 1972;26 4:239-57;
Serrano-Heras G, Diaz-Maroto I, Castro-Robles B, Carrion B, Perona-Moratalla AB, Gracia J, et al. Isolation and Quantification of Blood Apoptotic Bodies, a Non-invasive Tool to Evaluate Apoptosis in Patients with Ischemic Stroke and Neurodegenerative Diseases. Biol Proced Online. 2020;22:17;
Shi Y. Mechanical aspects of apoptosome assembly. Curr Opin Cell Biol. 2006;18 6:677-84;
Yin F, Zhou H, Fang Y, Li C, He Y, Yu L, et al. Astragaloside IV alleviates ischemia reperfusion-induced apoptosis by inhibiting the activation of key factors in death receptor pathway and mitochondrial pathway. J Ethnopharmacol. 2020;248:112319;
Santos LC, Vogel R, Chipuk JE, Birtwistle MR, Stolovitzky G, Meyer P. Mitochondrial origins of fractional control in regulated cell death. Nat Commun. 2019;10 1:1313;
Adams JM, Cory S. Bcl-2-regulated apoptosis: mechanism and therapeutic potential. Curr Opin Immunol. 2007;19 5:488-96;
Tsujimoto Y, Cossman J, Jaffe E, Croce CM. Involvement of the bcl-2 gene in human follicular lymphoma. Science. 1985;228 4706:1440-3;
Hartman ML, Czyz M. BCL-w: apoptotic and non-apoptotic role in health and disease. Cell Death Dis. 2020;11 4:260;
Gross A, Jockel J, Wei MC, Korsmeyer SJ. Enforced dimerization of BAX results in its translocation, mitochondrial dysfunction and apoptosis. Embo J. 1998;17 14:3878-85;
Minn AJ, Swain RE, Ma A, Thompson CB. Recent progress on the regulation of apoptosis by Bcl-2 family members. Adv Immunol. 1998;70:245-79;
Strasser A. The role of BH3-only proteins in the immune system. Nat Rev Immunol. 2005;5 3:189-200;
Deng J, Carlson N, Takeyama K, Dal Cin P, Shipp M, Letai A. BH3 profiling identifies three distinct classes of apoptotic blocks to predict response to ABT-737 and conventional chemotherapeutic agents. Cancer Cell. 2007;12 2:171-85;
Igaki T, Miura M. Role of Bcl-2 family members in invertebrates. Biochim Biophys Acta. 2004;1644 2-3:73-81;
Ayllon N, Villar M, Galindo RC, Kocan KM, Sima R, Lopez JA, et al. Systems biology of tissue-specific response to Anaplasma phagocytophilum reveals differentiated apoptosis in the tick vector Ixodes scapularis. PLoS Genet. 2015;11 3:e1005120;
Zhou J, Gong H, Zhou Y, Xuan X, Fujisaki K. Identification of a glycine-rich protein from the tick Rhipicephalus haemaphysaloides and evaluation of its vaccine potential against tick feeding. Parasitol Res. 2006;100 1:77-84;
Yu X, Gong H, Zhou Y, Zhang H, Cao J, Zhou J. Differential sialotranscriptomes of unfed and fed Rhipicephalus haemaphysaloides, with particular regard to differentially expressed genes of cysteine proteases. Parasit Vectors. 2015;8:597;
Wang Y, Hu S, Tuerdi M, Yu X, Zhang H, Zhou Y, et al. Initiator and executioner caspases in salivary gland apoptosis of Rhipicephalus haemaphysaloides. Parasit Vectors. 2020;13 1:288;
Cruz LM, Trefflich S, Weiss VA, Castro MAA. Protein Function Prediction. Methods Mol Biol. 2017;1654:55-75;
Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40 Web Server issue:W597-603;
Nijhof AM, Balk JA, Postigo M, Jongejan F. Selection of reference genes for quantitative RT-PCR studies in Rhipicephalus (Boophilus) microplus and Rhipicephalus appendiculatus ticks and determination of the expression profile of Bm86. BMC Mol Biol. 2009;10:112;
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25 4:402-8;
Ginzinger DG. Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Exp Hematol. 2002;30 6:503-12;
Kim HS, Smithies O. Recombinant fragment assay for gene targetting based on the polymerase chain reaction. Nucleic Acids Res. 1988;16 18:8887-903;
Wang Y, Li Z, Zhou Y, Cao J, Zhang H, Gong H, et al. Specific histamine binding activity of a new lipocalin from Hyalomma asiaticum (Ixodidae) and therapeutic effects on allergic asthma in mice. Parasit Vectors. 2016;9 1:506;
Greenfield EA, DeCaprio J, Brahmandam M. Making Weak Antigens Strong: Cross-Linking Peptides to KLH with Maleimide. Cold Spring Harb Protoc. 2018;2018 10:pdb prot100016;
Grabowski JM, Tsetsarkin KA, Long D, Scott DP, Rosenke R, Schwan TG, et al. Flavivirus Infection of Ixodes scapularis (Black-Legged Tick) Ex Vivo Organotypic Cultures and Applications for Disease Control. MBio. 2017;8 4;
Chen Z, Yang X, Bu F, Yang X, Yang X, Liu J. Ticks (acari: ixodoidea: argasidae, ixodidae) of China. Exp Appl Acarol. 2010;51 4:393-404;
Yu Z, Wang H, Wang T, Sun W, Yang X, Liu J. Tick-borne pathogens and the vector potential of ticks in China. Parasit Vectors. 2015;8:24;
Zhang J, Liu Q, Wang D, Li W, Beugnet F, Zhou J. Epidemiological survey of ticks and tick-borne pathogens in pet dogs in south-eastern China. Parasite. 2017;24:35;
Yu X, Zhou Y, Cao J, Zhang H, Gong H, Zhou J. Caspase-1 participates in apoptosis of salivary glands in Rhipicephalus haemaphysaloides. Parasit Vectors. 2017;10 1:225;
Brachmann CB, Jassim OW, Wachsmuth BD, Cagan RL. The Drosophila bcl-2 family member dBorg-1 functions in the apoptotic response to UV-irradiation. Curr Biol. 2000;10 9:547-50;
Colussi PA, Quinn LM, Huang DC, Coombe M, Read SH, Richardson H, et al. Debcl, a proapoptotic Bcl-2 homologue, is a component of the Drosophila melanogaster cell death machinery. J Cell Biol. 2000;148 4:703-14;
Igaki T, Kanuka H, Inohara N, Sawamoto K, Nunez G, Okano H, et al. Drob-1, a Drosophila member of the Bcl-2/CED-9 family that promotes cell death. Proc Natl Acad Sci U S A. 2000;97 2:662-7;
Chou JJ, Li H, Salvesen GS, Yuan J, Wagner G. Solution structure of BID, an intracellular amplifier of apoptotic signaling. Cell. 1999;96 5:615-24;
McDonnell JM, Fushman D, Milliman CL, Korsmeyer SJ, Cowburn D. Solution structure of the proapoptotic molecule BID: a structural basis for apoptotic agonists and antagonists. Cell. 1999;96 5:625-34;
Scopinho Furquim KC, Bechara GH, Camargo Mathias MI. Death by apoptosis in salivary glands of females of the tick Rhipicephalus sanguineus (Latreille, 1806) (Acari: Ixodidae). Exp Parasitol. 2008;119 1:152-63;
Oltvai ZN, Milliman CL, Korsmeyer SJ. Bcl-2 heterodimerizes in vivo with a conserved homolog, Bax, that accelerates programmed cell death. Cell. 1993;74 4:609-19;
Brooks C, Dong Z. Regulation of mitochondrial morphological dynamics during apoptosis by Bcl-2 family proteins: a key in Bak? Cell Cycle. 2007;6 24:3043-7;
Bettaieb A, Dubrez-Daloz L, Launay S, Plenchette S, Rebe C, Cathelin S, et al. Bcl-2 proteins: targets and tools for chemosensitisation of tumor cells. Curr Med Chem Anticancer Agents. 2003;3 4:307-18;

Table 1 Effect of knocking-down RhBcl-2 and RhBax on tick feeding behavior.

Test groups	Attachment rate at 24h (%)	Engorgement rate (%)	Death rate (%)	egg laying rate (%)	Hatchability rate (%)
dsLuciferase	76.67%	73.33%	20.00%	65.00%	65.00%
dsRhBcl-2	85.00%	0.00%	95.00%	0.00%	0.00%
dsRhBax	70.00%	61.67%	33.33%	48.33%	48.33%

Download PDF

Journal Publication

published 04 Aug, 2021

Read the published version in Parasites & Vectors →

Editorial decision: Major Revision
17 Feb, 2021
Review #1 received at journal
15 Feb, 2021
Review #2 received at journal
15 Feb, 2021
Reviewer #2 agreed at journal
02 Feb, 2021
Reviewer #1 agreed at journal
01 Feb, 2021
Reviewers invited by journal
30 Jan, 2021
First submitted to journal
19 Jan, 2021
Editor invited by journal
19 Jan, 2021
Editor assigned by journal
19 Jan, 2021
Submission checks completed at journal
19 Jan, 2021

You are reading this latest preprint version

Identification of the Bcl-2 and Bax Homologs From the Tick Rhipicephalus Haemaphysaloides and Their Function in the Degeneration of Tick Salivary Glands

Status:

Journal Publication

Version 1

Abstract

Figures

1 Introduction

2 Material And Methods

2.1 Tick Feeding and tissue collection

2.2 RNA extraction and cDNA synthesis

2.3 Cloning, sequence analysis, and expression of RhBcl-2 and RhBax

2.4 RT-qPCR analyses

2.5 Expression of recombinant RhBcl-2 or RhBax and antibody production

2.6 GST-Pull Down of RhBcl-2 and RhBax

2.7 Cells and Transient Co-transfection assays

2.8 Antibody generation and western blotting

2.9 RNAi of RhBcl-2 or RhBax

2.10 TUNEL staining

3 Results

3.1 Sequencing analysis of RhBcl-2 and RhBax

3.2 Transcription and translation of RhBcl-2 and RhBax profiles in different stages, tissues, and feeding periods

3.3 Expression of RhBcl-2 and RhBax

3.4 RhBcl-2 binds with RhBax by BH3 domain

3.5 RNAi assay of RhBcl-2 and RhBax

Discussion

Abbreviations

Declarations

Ethics approval and consent to participate

Acknowledgments

References

Tables

Supplementary Files

Status:

Journal Publication

Version 1