Differential sialotranscriptomes of unfed and fed Rhipicephalus haemaphysaloides, with particular regard to differentially expressed genes of cysteine proteases
© Yu et al. 2015
Received: 27 August 2015
Accepted: 12 November 2015
Published: 17 November 2015
Rhipicephalus haemaphysaloides, a hard tick, is a common ectoparasite and can be found in many countries. It is recognized as the primary vector of bovine babesiosis in the south of China. During blood feeding, the tick’s salivary glands secret numerous essential multifunctional proteins. In this study, a R. haemaphysaloides salivary gland transcriptome was described following the production and analysis of the transcripts from the two cDNA libraries of unfed and fed female ticks. The study focused on the differentially expressed genes and cysteine proteases, which play essential roles in the tick life cycle, that were detected most commonly in the up-regulation libraries.
The sialotranscriptome was assembled and analyzed though bioinformatic tools and the cysteine protease which is differentially expressed form sialotranscriptome were confirmed by Real-time PCR in salivary glands and different developments of ticks.
On the basis of sequence similarities with other species in various databases, we analyzed the unfed and fed sialotranscriptome of R. haemaphysaloides to identify the differentially expressed proteins secreted from the salivary glands during blood feeding and to investigate their biological functions. There were 25,113 transcripts (35 % of the total assembled transcripts) that showed significant similarity to known proteins with high BLAST from other species annotated. In total, 88 % and 89 % of the sequencing reads could be mapped back to assembled sequences in the unfed and fed library, respectively. Comparison of the abundance of transcripts from similar contigs of the two salivary gland cDNA libraries allowed the identification of differentially expressed genes. In total, there were 1179 up-regulated genes and 574 down-regulated genes found by comparing the two libraries. Twenty-five predicted cysteine proteases were screened from the transcript databases, whereas only six protein molecules were confirmed by gene cloning and molecular expression in E.coli which all belonged to the cysteine protease family. Bioinformatic evolutionary analysis showed the relationship of cysteine proteases in ticks with those of other species, suggesting the origin and conservation of these genes. Analysis of sequences from different tick species indicated the further relationships among the proteases, suggesting the closely related function of these genes. Thus, we confirmed their changes in unfed, fed and engorged ticks and salivary glands. The dynamic changes revealed their important roles in the tick life cycle.
Our survey provided an insight into the R. haemaphysaloides sialotranscriptome. The dynamic changes of cysteine proteases in ticks will assist further study of these proteases, which may contribute to the development of anti-tick vaccines or drugs, as well as improving understanding of the roles of cysteine proteases in the tick life cycle.
KeywordsSialotranscriptome Rhipicephalus haemaphysaloides Cysteine proteases Dynamic change
As a kind of blood-feeding arthropod, ticks can transmit viruses, bacteria and protozoa in their meals . Up to now, there are over 800 species described worldwide which are classified in Ixodidae (hard ticks) and Argasidae (soft ticks) primarily . Rhipicephalus haemaphysaloides is a three-host tick belonging to the Ixodidae and is widely distributed in China, India, and other South Asian countries . This tick is a major vector of bovine babesiosis in China  and can also transmit the Kyasanur Forest disease virus .
Within the blood-feeding, ticks possess salivary glands that secrete bioactive substances, which can exhibit a range of pharmacological properties to thwart the host defense mechanisms in response to attachment [6, 7]. The components of the saliva are of major importance for the tick’s survival, helping it feed and evade host defenses, hemostatic factors and the inflammatory response . Proteases are one of the most important components of tick saliva and essential for the life cycle of the ectoparasite.
Cysteine proteases are ancient conserved proteases that are involved in different physiological processes . Most of these proteases belong to the papain-like superfamily and are associated with the development of hematophagous arthropod ectoparasites . Ticks express cysteine peptidases with important roles in physiological events that are crucial to the ectoparasitic lifestyle, including the digestion of host blood, embryogenesis and innate immunity .
In this study, we analyzed the sialotranscriptome of the salivary glands of unfed (unattached) and fed (3 or 4 days after attached) adult ticks (Additional file 1). There were 1179 up-regulated genes and 574 down-regulated genes detected from the differential expression databases. For functional annotation of the unique transcripts, we used BLASTx, comparing them against different databases and, finally, four up-regulated and two associated cysteine proteases, namely cathepsin B (CATB, KT194088), cathepsin L (CATL, KT194089), caspase–1 (CASP1, KT194090), caspase–8 (CASP8, KT194091), autophagy protease 4B (ATG4B, KT194092) and autophagy protease 4D (ATG4D, KT194093), were cloned successfully from the cDNA library of the salivary glands of R. haemaphysaloides. The confirmation of these genes will contribute to further research in vitro. The dynamic analysis of these target genes was shown in different developmental stages of R. haemaphysaloides by Q-PCR, and suggests their important roles during blood feeding. To our knowledge, this is the first analysis of the transcriptome of the salivary glands of female R. haemaphysaloides ticks. The characterization of the components of tick saliva, especially the proteases (cysteine proteases), is likely to be of value in the design of novel methods or drugs for the control of ticks and tick-borne diseases, as well as when searching for proteins that may have potential use in research on medical and veterinary diseases.
Collection of ticks and salivary glands
The R. haemaphysaloides colonies were maintained in the laboratory as described previously . For tissue collection, the salivary glands were dissected under a light microscope . The sample materials were stored at −80 °C until use.
cDNA library construction and sequencing
Total RNA was extracted from the unfed and fed salivary glands of female R. haemaphysaloides using TRIzol Reagent (Invitrogen, The Netherlands) according to the manufacturer’s protocol. The cDNA from two RNA-seq sequencing libraries was sheared to an average fragment size of about 300 values and was purified with Ampure beads. RNA-seq libraries were constructed according to the Illumina manufacturer’s instructions for 100-bp paired-ends, and sequenced. Raw reads were filtered to produce clean reads prior to assembly. Initially, the RNAs were extracted and constructed for two Illumina/Solexa cDNA libraries.
De novo assembly of transcriptome data
Low-quality regions in raw reads and adaptors were trimmed with ea-utils  prior to analyses. Read quality was then assessed, revealing that the mean quality scores of sequence reads was around 37 values, which suggested high sequencing quality. Transcriptome assembly was performed with Trinity assembles software  to obtain high-quality contigs. The Trinity assembly program first combines reads with a certain length of overlap to form longer fragments without N; these are called contigs. The sequencing reads are then mapped back to contigs; with paired-end reads it is able to detect contigs from the same transcript as well as measuring the distances between these contigs. Subsequently, Trinity connects the contigs using N to represent unknown sequences between each pair of contigs, following which Scaffolds are made.
Bioinformatic tools used and differential expression analysis
The BLASTn , CAP3 assembler  and ClustalW software  were used to compare, assemble, and align high quality expressed sequence tags, ESTs, respectively. For functional annotation of the transcripts we used BLASTx  to compare the nucleotide sequences with the non-redundant (NR) protein database of the National Center of Biological Information (NCBI) and to the Gene Ontology (GO) database . The gene expression profiles were compared by mapping RNA-seq reads using Bowtie 2 2.1.0 . Assembled sequences with high BLAST similarity to known protein sequences (E-value cut-off of 1E-6) from other species were annotated and the GO functional annotations were extracted using the Swiss-Pro (http://www.uniprot.org/) BLAST result by comparing with EMBL Uniprot eggNOG/GO Pathways databases. Analysis of GO terms was subsequently performed using a custom script. The GO terms belonging to cellular components, biological processes and molecular functions were listed. Bi-directional best hit (BBH) was used to search against the KEGG database  to obtain the KO (Reference pathway) number of the KEGG Annotation . The KO (Reference pathway) number of the transcriptome was also obtained, according to KEGG Annotation.
To detect the changes in global gene expression in the different tissues, we applied the MA-plot-based method with a random sampling model  to identify the differentially expressed genes by comparing the unfed library with the fed library. Genes with fold change >3 and P value <0.001 were regarded as differentially expressed genes.
Analysis of relative expression by quantitative real-time PCR
The Primers used in real-time PCR
Forward primer (5′-3′)
The data were normalized to the elongation factor-1 gene (EF-1) (accession number AB836665) . Relative gene expression data were analyzed using the 2—△Ct method [21, 22], and △Ct values were calculated by subtracting the average EF-1 Ct values from those for the average target gene.
Mean ± SE values for each group (n = 4) were calculated, and two-tailed t tests were used to compare differences between groups with Graphpad PRISM 5 software (GraphPad Software Inc., La Jolla, California).
Statistics sequencing amount for two fastq files
The summary of transcripts of salivary glands of R.haemaphysaloides
Total number of reference transcripts
The number of transcripts with BLASTX-hit
The number of transcripts with KEGG annotation
The number of transcripts detected with RPKM
The number of up-regulated genes
The number of down-regulated genes
GO from up-regulated expressed genes
GO cellular component /molecular function /biological process (up-regulated expressed genes)
GO from down-regulated expressed genes
GO cellular component /molecular function /biological process (down-regulated expressed genes)
To investigate the functions of these differentially expressed genes, all the differentially expressed genes were mapped to the GO database and compared with the whole transcriptome background. BiNGO  was used to identify enrichment GO terms by using custom annotation files from the transcriptome on the basis of the hypergeometric test (P-value <0.01) (Additional file 1). All the enrichment GO terms from the up-regulated expressed genes could be categorized into 135 GO terms (29 cellular component terms, 33 molecular function terms, and 73 biological process terms). The enrichment GO from down-regulated expressed genes could be categorized into 15 GO terms (5 cellular component terms, 9 molecular function terms, and 1 biological process term) (Table 3).
The summary of GO terms about proteases
(The number of related gene)
ubiquitin-specific protease activity
(The number of related gene)
apoptotic protease activator activity
(The number of related gene)
SUMO-specific protease activity
(The number of related gene)
small conjugating protein-specific protease activity
(The number of related gene)
NEDD8-specific protease activity
(The number of related gene)
ubiquitin-specific protease activator activity
(The number of related gene)
Bioinformatic evolutionary analysis
Dynamic profiling of cysteine proteases in salivary glands of ticks
Dynamic profiling of cysteine proteases in different developmental stages of ticks
Most of the proteins displayed in Fig. 3 have been functionally characterized in mammals, while few studies have described the exact function of caspases and autophagy proteases in ticks, although they possess only innate immunity. Almost all the cathepsins in ticks are involved in the process of blood digestion . The peptidases identified from different tick species with a proposed role in hemoglobin digestion are mainly cysteine peptidases, with some serine and metallopeptidases [25, 26]. As in humans, the caspases have the ability to regulate three alternative cell death pathways: apoptosis, pyroptosis, and necroptosis . Moreover, recent work has shown that cathepsin B and cathepsin D regulate the inflammasome-dependent and -independent macrophage responses induced by cytosolic flagellin . The research also revealed that cathepsin B contributes to NAIP5/NLRC4 inflammasome-induced pyroptosis and interleukin-1α (IL-1α) and IL-1β production in response to cytosolic flagellin . Unfortunately, there is little knowledge about the relationships among the cathepsins, caspases and autophagy proteases in ticks. It will be of great value to identify the function of these cysteine proteases because of their important roles, and this may contribute to the development of candidate vaccines or drugs for tick control by RNA interference.
When feeding on their hosts, ticks need to deal with host hemostasis, inflammation and immunity. Although recent progress in transcription research on hard ticks has shown that hundreds of different proteins are expressed in their salivary glands, many of them are known only as salivary proteins with unknown function .
Tick feeding is a slow and uninterruptible process, and blood digestion takes place in the gut epithelium, which is different with insect blood-feeders that feed and digest blood in the gut lumen rapidly with neutral pH [29, 30]. In this progress, the salivary glands act as a “pivot” in which anti-host defense molecules and digestion proteases are secreted and released to the host vasculature and to the tick’s midgut with the blood. Pathogens such as Babesia [31, 32], Anaplasma  and Borrelia [34, 35] infect tick salivary glands and are injected into the host during their meals . Ticks possess defense mechanisms that allow them to maintain pathogens and commensal microbes at a certain level that does not impair their fitness and further development. Whether any of the proteases also play a role in immune defense is unknown. However, evidence from insects suggests that metalloproteases may be important in cellular immune defense . In mosquitoes, serine proteases are reported to be up-regulated in response to invasion of the hemolymph by malarial parasites, contributing to the normal innate immune response .
Cysteine proteases are a component of the multi-enzyme hemoglobinolytic model for hard ticks; its numerous members play different roles in tick life cycles. To our knowledge, CATB and CATL are believed to be involved in the digestion of blood ; caspases seem also to participate in this process. In mammals, caspases lead to apoptosis and inflammation and may be associated with autophagy proteases, which are also involved in inflammation and cell death . Unfortunately, their functions in ticks are still unclear. After engorgement, the salivary glands become withered and apoptotic and may even vanish. This process is similar to apoptosis in mammals, but whether caspases and autophagy proteases are involved in this mechanism requires further confirmation.
Cathepsins, another component of the tick multiple enzyme system, are believed to be involved in the digestion of blood. The current knowledge of the molecular characteristics of tick digestive enzymes began to be assembled in the 1980s to 1990s by isolation and partial characterization of acidic aspartic peptidases of cathepsin D from soft and hard ticks [40, 41]. Later, Mendiola et al. reported that aspartic (cathepsin D-like) and cysteine (cathepsin L-like) peptidases are the major hemoglobinolytic enzymes in R. microplus. Two cathepsin L-type cysteine peptidases were partially characterized and cloned from the midgut of Haemaphysalis longicornis (H. longicornis) . Another cysteine peptidase gene homologous to cathepsin L (BmCL1) was shown to be expressed in the gut of partially engorged R. microplus females, and recombinant BmCL1 was optimally active against bovine hemoglobin at acidic pH [43, 44]. The research on Ixodes ricinus shows more detail about the cathepsins. There is a mechanistic model of the proteolytic pathway of hemoglobin degradation in the digestive vesicles of I. ricinus gut cells . Cathepsin D (CATD), supported by cathepsin L (CATL) and legumain (AE), is responsible for the primary events in the cleavage of hemoglobin. Subsequently, cathepsin B (CATB) and cathepsin L (CATL) participate in the secondary digestion, to generate smaller fragments. After the dipeptidase activities of CATB and cathepsin C (CATC), small fragments are degraded into dipeptides. Most CATB and CATL have been identified in the midgut in different tick species: for example, IrCB1 from I. ricinus [25, 26, 45], longipain from H. longicornis , and BmCL1 from R. microplus [43, 44, 47, 48]. In I. ricinus, IrCL1 was identified in the midgut, salivary glands, ovary and malpighian tubes [26, 49, 50]. Based on our report of the sialotranscriptome, we have identified for the first time that CATB and CATL existed in salivary glands from R. haemaphysaloides. These two peptidases may be involved in hemoglobin digestion, and therefore their mRNA expression reaches a peak in engorged ticks and their salivary glands. Given that digestion takes place in the gut epithelium, the midgut of R. haemaphysaloides may also contain CATB and CATL. When the blood flows into the gut lumen, it may stimulate and trigger the transcription of peptidases in gut epithelium. Under this circumstance, the mRNA expression in salivary glands may be affected and this may be the reason that there is always low transcription level in the fed state.
Recently, most of the knowledge about caspases has been derived from research on the human or mouse. There is little known about their function in insects, and much less in ticks. However, in-depth analysis on human caspases shows much more detail and offers suggestions for further research in other species. As a member of cysteine proteases family, Caspases have primary specificity for aspartic acid (Asp) residues; they cleave their substrates after tetrapeptide sequences containing Asp in the P1 position. All caspases are synthesized as inactive single-chain zymogens (procaspases) initially, and then processed into their active forms. Additional signals are required to the initiation of caspase activation pathways . Protein interaction domains are component of long prodomains in initiator caspases. For example, CARDs are contained in caspases−1, −2, −4, −5, −9, − 11, and −12, and death effector domains (DEDs) are in caspases-8 and −10 . Caspase−8 as an initiator caspase is activated via their DED-mediated interactions within the death-inducing signaling complex (DISC) in the extrinsic pathway. Ultimately the downstream effector caspases−3, −6, and −7 are activated and responsible for the classical phenotypic changes associated with apoptosis . Caspase-1 is regarded as the prototypical inflammatory caspase , and is responsible for the processing of proIL-1β and proIL-18 . Caspase−8 is a mediator of inflammation . Recent studies have revealed distinct roles for caspase−8, which are associated with the extrinsic apoptotic pathway. It is involved in the regulation of inflammation and is also proposed to have an anti-inflammatory role . In addition, caspase−8 is proposed to cleave proIL-1β into its active mature form [54–61], to be incorporated and activated within inflammasome complexes [62, 63], and to activate caspase-1 directly in an inflammasome-independent manner [64, 65].
As mentioned above, it is believed that there is only an innate immune system in ticks. This innate immunity plays really important roles in the tick life cycle, although it may not be as complex and complete as that in mammals. The sialotranscriptome of R. haemaphysaloides shows at least 3 or 4 caspases in this tick species, but only caspase−1 and caspase−8 were cloned successfully from the cDNA library. CASP8 reached its peak level after engorgement in our experiment; at this time, ticks are full of host blood, which contains cytokines and perhaps pathogens. As an initiator caspase, CASP8 may trigger and enhance the innate immunity, with CASP1, to defend against harmful effects as described above. The incomplete immune system of the tick may have restricted function when compared with that in mammals. It is still unclear how ticks deal with the cytokines ingested in blood, and whether caspases participate in the management of cytokines and in the mechanism of pathogen defense. Western blot showed that, in female adult R. haemaphysaloides, CASP1 may contain ~35 kDa monomers and ~70 kDa dimers (unpublished data), coinciding with that in humans. With reference to studies on human caspases, we will undertake further research on the function of caspases in R. haemaphysaloides using RNAi.
Autophagy related (ATG) genes are a complex and mysterious family containing numerous members. There were more than 30 of these genes characterized in yeast originally, and many orthologs have been identified as autophagy regulators in higher eukaryotes [66, 67]. Proteases participate in several stages of autophagy. In the initial steps of macroautophagy, the formation of autophagic vesicles requires the conjugation of phosphatidylethanolamine with ATG8 . Following translation, ATG8 is cleaved by ATG4 (the cysteine family protease) in yeast cells. Subsequently, the resulting ATG8G116 is involved in a ubiquitin-like conjugation reaction catalyzed by ATG7 and ATG3 . Similarly, ATG4 participates in processing of three mammalian homologs of ATG8 in mammalian cells, which is crucial for the autophagic pathway [39, 69, 70]. In the sialotranscriptome of R. haemaphysaloides, we found only two autophagy-related (ATG) genes; these were classified into ATG4B and ATG4D on the basis of their sequences in BLAST databases. After engorgement, the mRNA expression of ATG4B and ATG4D rises to a peak in adult ticks and their salivary glands, when the salivary glands are approaching apoptosis. It is thought that the ATG-protein may be involved in this process, and that the constitution of autophagosomes in ticks may be much simpler than that in yeast for the limited ATG genes, although its mechanism of action is still unknown.
It is generally suggested that the major function of lysosomal proteases is to maintain cellular homeostasis and differentiation by recycling cellular content . Cathepsin B and Cathepsin L, as the most abundant lysosomal proteases , are involved directly in the execution of autophagy [72, 73]. Cathepsin B has also been reported to regulate the activity of caspases . In addition, autophagy-related proteins are acknowledged to be involved in inflammation, infection and cancer , meaning that apoptosis and pyroptosis are closely associated with caspases. In view of these findings, it is thought that there must be a closely regulated relationship among the cysteine proteases.
Analysis of the sialotranscriptome of R. haemaphysaloides using two cDNA libraries, from unfed and fed ticks, identified many transcripts coding for different proteins. On the basis of the database of differentially expressed genes, we identified four up-regulated cysteine proteases and two associated genes for further study. Using Q-PCR, the up-regulated genes were found to be CATB, ATG4B, ATG4D and CASP8, which is slightly different from results of the sialotranscriptome, although mRNA expression of all six target genes reached a peak after engorgement. Moreover, the phylogenetic trees showed that the R. haemaphysaloides cysteine protease sequences are dispersed into different clades, which contain sequences from other species, suggesting an ancient origin for these genes. The phylogram of different tick species demonstrates the close relation between caspases and autophagy-related proteins. Furthermore, we observed dynamic changes of mRNA expression in different developmental stages of R. haemaphysaloides, suggesting their multiple functions during blood feeding.
Taken together, these results improve our knowledge of the salivary components and distribution of cysteine proteases of R. haemaphysaloides that will make a contribution to further study, and may help in the research to identify candidate antigens for anti-tick vaccines and to discover drugs to treat inflammation and cancer.
The protocols were approved by the Institutional Animal Care and Use Committee of the Shanghai Veterinary Research Institute, and followed the misconduct policy of BMC Genomics, and authorized by the Animal Ethical Committee of Shanghai Veterinary Research Institute.
Project support was provided by the “National Key Basic Research Program (973 program) of China” (Grant No.2015CB150300).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- De la Fuente J, Estrada-Pena A, Venzal JM, Kocan KM, Sonenshine DE. Overview: Ticks as vectors of pathogens that cause disease in humans and animals. Front Biosci. 2008;13:6938–46.View ArticlePubMedGoogle Scholar
- Nava S, Guglielmone AA, Mangold AJ. An overview of systematics and evolution of ticks. Front Biosci (Landmark Ed). 2009;14:2857–77.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Yin H, Lu W, Luo J. Babesiosis in China. Trop Anim Health Prod. 1997;29(4 Suppl):11S–5S.View ArticlePubMedGoogle Scholar
- Bhat HR, Naik SV, Ilkal MA, Banerjee K. Transmission of Kyasanur Forest disease virus by Rhipicephalus haemaphysaloides ticks. Acta Virol. 1978;22(3):241–4.PubMedGoogle Scholar
- Bowman AS, Coons LB, Needham GR, Sauer JR. Tick saliva: recent advances and implications for vector competence. Med Vet Entomol. 1997;11(3):277–85.View ArticlePubMedGoogle Scholar
- Oaks JF, McSwain JL, Bantle JA, Essenberg RC, Sauer JR. Putative new expression of genes in ixodid tick salivary gland development during feeding. J Parasitol. 1991;77(3):378–83.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Sojka D, Francischetti IM, Calvo E, Kotsyfakis M. Cysteine proteases from bloodfeeding arthropod ectoparasites. Adv Exp Med Biol. 2011;712:177–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhou J, Liao M, Hatta T, Tanaka M, Xuan X, Fujisaki K. Identification of a follistatin-related protein from the tick Haemaphysalis longicornis and its effect on tick oviposition. Gene. 2006;372:191–8.View ArticlePubMedGoogle Scholar
- Aronesty E. Comparison of Sequencing Utility Programs. Open Bioinform J. 2013;7:1–8.View ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.View ArticlePubMed CentralPubMedGoogle Scholar
- Huang X, Madan A. CAP3: A DNA sequence assembly program. Genome Res. 1999;9(9):868–77.View ArticlePubMed CentralPubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database issue):D480–484.PubMed CentralPubMedGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.View ArticlePubMedGoogle Scholar
- Tomizawa MNH. High mortality caused by high dose of dsRNA in the green rice leafhopper Nephotettix cincticeps (Hemiptera: Cicadellidae). Appl Entomol Zool. 2013;48:553–9.View ArticleGoogle Scholar
- Ginzinger DG. Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Exp Hematol. 2002;30(6):503–12.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods (San Diego, Calif). 2001;25(4):402–8.View ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.View ArticlePubMedGoogle Scholar
- Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21(16):3448–9.View ArticlePubMedGoogle Scholar
- Sojka D, Franta Z, Horn M, Caffrey CR, Mares M, Kopacek P. New insights into the machinery of blood digestion by ticks. Trends Parasitol. 2013;29(6):276–85.View ArticlePubMedGoogle Scholar
- Horn M, Nussbaumerova M, Sanda M, Kovarova Z, Srba J, Franta Z, et al. Hemoglobin digestion in blood-feeding ticks: mapping a multipeptidase pathway by functional proteomics. Chem Biol. 2009;16(10):1053–63.View ArticlePubMed CentralPubMedGoogle Scholar
- Creagh EM. Caspase crosstalk: integration of apoptotic and innate immune signalling pathways. Trends Immunol. 2014;35(12):631–40.View ArticlePubMedGoogle Scholar
- Lage SL, Buzzo CL, Amaral EP, Matteucci KC, Massis LM, Icimoto MY, et al. Cytosolic flagellin-induced lysosomal pathway regulates inflammasome-dependent and -independent macrophage responses. Proc Natl Acad Sci U S A. 2013;110(35):E3321–3330.View ArticlePubMed CentralPubMedGoogle Scholar
- Briegel H, Lea AO. Relationship between protein and proteolytic activity in the midgut of mosquitoes. J Insect Physiol. 1975;21(9):1597–604.View ArticlePubMedGoogle Scholar
- Grandjean O, Aeschlimann A. Contribution to the study of digestion in ticks: histology and fine structure of the midgut ephithelium of Ornithodorus moubata, Murray (Ixodoidea, Argasidae). Acta Trop. 1973;30(4):193–212.PubMedGoogle Scholar
- Chauvin A, Moreau E, Bonnet S, Plantard O, Malandrin L. Babesia and its hosts: adaptation to long-lasting interactions as a way to achieve efficient transmission. Vet Res. 2009;40(2):37.View ArticlePubMed CentralPubMedGoogle Scholar
- Florin-Christensen M, Schnittger L. Piroplasmids and ticks: a long-lasting intimate relationship. Front Biosci (Landmark Ed). 2009;14:3064–73.View ArticleGoogle Scholar
- Kocan KM, de la Fuente J, Blouin EF, Coetzee JF, Ewing SA. The natural history of Anaplasma marginale. Vet Parasitol. 2010;167(2–4):95–107.View ArticlePubMedGoogle Scholar
- Steere AC, Coburn J, Glickstein L. The emergence of Lyme disease. J Clin Invest. 2004;113(8):1093–101.View ArticlePubMed CentralPubMedGoogle Scholar
- A PG. Borrelia: Molecular Biology, Host Interaction and Pathogenesis. Clin Infect Dis. 2011;52:965–5.Google Scholar
- Hajdusek O, Sima R, Ayllon N, Jalovecka M, Perner J, de la Fuente J, et al. Interaction of the tick immune system with transmitted pathogens. Front Cell Infection Microbiol. 2013;3:26.View ArticleGoogle Scholar
- Willott E, Tran HQ. Zinc and Manduca sexta hemocyte functions. J Insect Sci. 2002;2:6.View ArticlePubMed CentralPubMedGoogle Scholar
- Xu W, Huang FS, Hao HX, Duan JH, Qiu ZW. Two serine proteases from Anopheles dirus haemocytes exhibit changes in transcript abundance after infection of an incompatible rodent malaria parasite, Plasmodium yoelii. Vet Parasitol. 2006;139(1–3):93–101.View ArticlePubMedGoogle Scholar
- Arroyo DS, Gaviglio EA, Peralta Ramos JM, Bussi C, Rodriguez-Galan MC, Iribarren P. Autophagy in inflammation, infection, neurodegeneration and cancer. Int Immunopharmacol. 2014;18(1):55–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Akov S, Samish M, Galun R. Protease activity in female Ornithodoros tholozani ticks. Acta Trop. 1976;33(1):37–52.PubMedGoogle Scholar
- Vundla WRMBM, Pearson DJ, Labongo VL. Characterization of aspartic proteinases from the gut of the tick, Rhipicephalus appendiculatus neuman. Insect Biochem Mol Biol. 1992;22:405–10.View ArticleGoogle Scholar
- Mulenga A, Sugimoto C, Ingram G, Ohashi K, Onuma M. Molecular cloning of two Haemaphysalis longicornis cathepsin L-like cysteine proteinase genes. J Vet Med Sci. 1999;61(5):497–502.View ArticlePubMedGoogle Scholar
- Renard G, Garcia JF, Cardoso FC, Richter MF, Sakanari JA, Ozaki LS, et al. Cloning and functional expression of a Boophilus microplus cathepsin L-like enzyme. Insect Biochem Mol Biol. 2000;30(11):1017–26.View ArticlePubMedGoogle Scholar
- Renard G, Lara FA, de Cardoso FC, Miguens FC, Dansa-Petretski M, Termignoni C, et al. Expression and immunolocalization of a Boophilus microplus cathepsin L-like enzyme. Insect Mol Biol. 2002;11(4):325–8.View ArticlePubMedGoogle Scholar
- Franta Z, Frantova H, Konvickova J, Horn M, Sojka D, Mares M, et al. Dynamics of digestive proteolytic system during blood feeding of the hard tick Ixodes ricinus. Parasites Vectors. 2010;3:119.View ArticlePubMed CentralPubMedGoogle Scholar
- Tsuji N, Miyoshi T, Battsetseg B, Matsuo T, Xuan X, Fujisaki K. A cysteine protease is critical for Babesia spp. transmission in Haemaphysalis ticks. PLoS Pathog. 2008;4(5):e1000062.View ArticlePubMed CentralPubMedGoogle Scholar
- Cruz CE, Fogaca AC, Nakayasu ES, Angeli CB, Belmonte R, Almeida IC, et al. Characterization of proteinases from the midgut of Rhipicephalus (Boophilus) microplus involved in the generation of antimicrobial peptides. Parasites Vectors. 2010;3:63.View ArticlePubMed CentralPubMedGoogle Scholar
- Clara RO, Soares TS, Torquato RJ, Lima CA, Watanabe RO, Barros NM, et al. Boophilus microplus cathepsin L-like (BmCL1) cysteine protease: specificity study using a peptide phage display library. Vet Parasitol. 2011;181(2–4):291–300.View ArticlePubMedGoogle Scholar
- Franta Z, Sojka D, Frantova H, Dvorak J, Horn M, Srba J, et al. IrCL1 - the haemoglobinolytic cathepsin L of the hard tick, Ixodes ricinus. Int J Parasitol. 2011;41(12):1253–62.View ArticlePubMedGoogle Scholar
- Sojka D, Franta Z, Horn M, Hajdusek O, Caffrey CR, Mares M, et al. Profiling of proteolytic enzymes in the gut of the tick Ixodes ricinus reveals an evolutionarily conserved network of aspartic and cysteine peptidases. Parasites Vectors. 2008;1(1):7.View ArticlePubMed CentralPubMedGoogle Scholar
- Salvesen GS, Walsh CM. Functions of caspase 8: the identified and the mysterious. Semin Immunol. 2014;26(3):246–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Garlanda C, Dinarello CA, Mantovani A. The interleukin-1 family: back to the future. Immunity. 2013;39(6):1003–18.View ArticlePubMed CentralPubMedGoogle Scholar
- Kaiser WJ, Upton JW, Long AB, Livingston-Rosanoff D, Daley-Bauer LP, Hakem R, et al. RIP3 mediates the embryonic lethality of caspase-8-deficient mice. Nature. 2011;471(7338):368–72.View ArticlePubMed CentralPubMedGoogle Scholar
- Hise AG, Tomalka J, Ganesan S, Patel K, Hall BA, Brown GD, et al. An essential role for the NLRP3 inflammasome in host defense against the human fungal pathogen Candida albicans. Cell Host Microbe. 2009;5(5):487–97.View ArticlePubMed CentralPubMedGoogle Scholar
- Gringhuis SI, Kaptein TM, Wevers BA, Theelen B, van der Vlist M, Boekhout T, et al. Dectin-1 is an extracellular pathogen sensor for the induction and processing of IL-1beta via a noncanonical caspase−8 inflammasome. Nat Immunol. 2012;13(3):246–54.Google Scholar
- Vince JE, Wong WW, Gentle I, Lawlor KE, Allam R, O'Reilly L, et al. Inhibitor of apoptosis proteins limit RIP3 kinase-dependent interleukin-1 activation. Immunity. 2012;36(2):215–27.View ArticlePubMedGoogle Scholar
- Gross O, Poeck H, Bscheider M, Dostert C, Hannesschlager N, Endres S, et al. Syk kinase signalling couples to the Nlrp3 inflammasome for anti-fungal host defence. Nature. 2009;459(7245):433–6.View ArticlePubMedGoogle Scholar
- Maelfait J, Vercammen E, Janssens S, Schotte P, Haegman M, Magez S, et al. Stimulation of Toll-like receptor 3 and 4 induces interleukin-1beta maturation by caspase–8. J Exp Med. 2008;205(9):1967–73.Google Scholar
- Bossaller L, Chiang PI, Schmidt-Lauber C, Ganesan S, Kaiser WJ, Rathinam VA, et al. Cutting edge: FAS (CD95) mediates noncanonical IL-1beta and IL-18 maturation via caspase−8 in an RIP3-independent manner. J Immunol. 2012;189(12):5508–12.Google Scholar
- Dickens LS, Powley IR, Hughes MA, MacFarlane M. The ‘complexities’ of life and death: death receptor signalling platforms. Exp Cell Res. 2012;318(11):1269–77.View ArticlePubMedGoogle Scholar
- Antonopoulos C, El Sanadi C, Kaiser WJ, Mocarski ES, Dubyak GR. Proapoptotic chemotherapeutic drugs induce noncanonical processing and release of IL-1beta via caspase−8 in dendritic cells. J Immunol. 2013;191(9):4789–803.Google Scholar
- Gurung P, Anand PK, Malireddi RK, Vande Walle L, Van Opdenbosch N, Dillon CP, et al. FADD and caspase−8 mediate priming and activation of the canonical and noncanonical Nlrp3 inflammasomes. J Immunol. 2014;192(4):1835–46.Google Scholar
- Man SM, Hopkins LJ, Nugent E, Cox S, Gluck IM, Tourlomousis P, et al. Inflammasome activation causes dual recruitment of NLRC4 and NLRP3 to the same macromolecular complex. Proc Natl Acad Sci U S A. 2014;111(20):7403–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Philip NH, Dillon CP, Snyder AG, Fitzgerald P, Wynosky-Dolfi MA, Zwack EE, et al. Caspase−8 mediates caspase-1 processing and innate immune defense in response to bacterial blockade of NF-kappaB and MAPK signaling. Proc Natl Acad Sci U S A. 2014;111(20):7385–90.Google Scholar
- Weng D, Marty-Roix R, Ganesan S, Proulx MK, Vladimer GI, Kaiser WJ, et al. Caspase-8 and RIP kinases regulate bacteria-induced innate immune responses and cell death. Proc Natl Acad Sci U S A. 2014;111(20):7391–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Virgin HW, Levine B. Autophagy genes in immunity. Nat Immunol. 2009;10(5):461–70.View ArticlePubMed CentralPubMedGoogle Scholar
- Mizushima N. The role of the Atg1/ULK1 complex in autophagy regulation. Curr Opin Cell Biol. 2010;22(2):132–9.View ArticlePubMedGoogle Scholar
- Ichimura Y, Kirisako T, Takao T, Satomi Y, Shimonishi Y, Ishihara N, et al. A ubiquitin-like system mediates protein lipidation. Nature. 2000;408(6811):488–92.View ArticlePubMedGoogle Scholar
- Tanida I, Komatsu M, Ueno T, Kominami E. GATE-16 and GABARAP are authentic modifiers mediated by Apg7 and Apg3. Biochem Biophys Res Commun. 2003;300(3):637–44.View ArticlePubMedGoogle Scholar
- Scherz-Shouval R, Sagiv Y, Shorer H, Elazar Z. The COOH terminus of GATE-16, an intra-Golgi transport modulator, is cleaved by the human cysteine protease HsApg4A. J Biol Chem. 2003;278(16):14053–8.View ArticlePubMedGoogle Scholar
- Turk B, Turk V. Lysosomes as “suicide bags” in cell death: myth or reality? J Biol Chem. 2009;284(33):21783–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Turk B, Stoka V, Rozman-Pungercar J, Cirman T, Droga-Mazovec G, Oresic K, et al. Apoptotic pathways: involvement of lysosomal proteases. Biol Chem. 2002;383(7–8):1035–44.PubMedGoogle Scholar
- Punnonen EL, Autio S, Marjomaki VS, Reunanen H. Autophagy, cathepsin L transport, and acidification in cultured rat fibroblasts. J Histochem Cytochem. 1992;40(10):1579–87.View ArticlePubMedGoogle Scholar