Ixodes scapularis and Ixodes ricinus tick cell lines respond to infection with tick-borne encephalitis virus: transcriptomic and proteomic analysis
© Weisheit et al. 2015
Received: 8 August 2015
Accepted: 11 November 2015
Published: 18 November 2015
Ixodid ticks are important vectors of a wide variety of viral, bacterial and protozoan pathogens of medical and veterinary importance. Although several studies have elucidated tick responses to bacteria, little is known about the tick response to viruses. To gain insight into the response of tick cells to flavivirus infection, the transcriptomes and proteomes of two Ixodes spp cell lines infected with the flavivirus tick-borne encephalitis virus (TBEV) were analysed.
RNA and proteins were isolated from the Ixodes scapularis-derived cell line IDE8 and the Ixodes ricinus-derived cell line IRE/CTVM19, mock-infected or infected with TBEV, on day 2 post-infection (p.i.) when virus production was increasing, and on day 6 p.i. when virus production was decreasing. RNA-Seq and mass spectrometric technologies were used to identify changes in abundance of, respectively, transcripts and proteins. Functional analyses were conducted on selected transcripts using RNA interference (RNAi) for gene knockdown in tick cells infected with the closely-related but less pathogenic flavivirus Langat virus (LGTV).
Differential expression analysis using DESeq resulted in totals of 43 and 83 statistically significantly differentially-expressed transcripts in IDE8 and IRE/CTVM19 cells, respectively. Mass spectrometry detected 76 and 129 statistically significantly differentially-represented proteins in IDE8 and IRE/CTVM19 cells, respectively. Differentially-expressed transcripts and differentially-represented proteins included some that may be involved in innate immune and cell stress responses. Knockdown of the heat-shock proteins HSP90, HSP70 and gp96, the complement-associated protein Factor H and the protease trypsin resulted in increased LGTV replication and production in at least one tick cell line, indicating a possible antiviral role for these proteins. Knockdown of RNAi-associated proteins Argonaute and Dicer, which were included as positive controls, also resulted in increased LGTV replication and production in both cell lines, confirming their role in the antiviral RNAi pathway.
This systems biology approach identified several molecules that may be involved in the tick cell innate immune response against flaviviruses and highlighted that ticks, in common with other invertebrate species, have other antiviral responses in addition to RNAi.
KeywordsTick Ixodes Flavivirus Tick-borne encephalitis virus Tick cell line Innate immunity Antiviral response
Ticks are hematophagous ectoparasites that are second only to mosquitoes in their importance as vectors of viral, bacterial and protozoan pathogens that cause human disease, and they are probably the most important vectors of livestock disease worldwide .
Tick-borne encephalitis virus (TBEV), a member of the family Flaviviridae, is a medically-important virus transmitted by ticks. TBEV is one of the most important tick-borne viruses in Europe, Russia and many parts of Asia, causing tick-borne encephalitis in humans with an estimated annual number of disease cases >10,000 [2, 3]. The Western European subtype of TBEV is transmitted by Ixodes ricinus ticks in Central, Eastern and Northern Europe, whereas the Siberian and Far-Eastern subtypes are transmitted by Ixodes persulcatus ticks in Siberia, parts of Russia, Latvia and Finland and the latter subtype additionally in Central and Eastern Asia including China and Japan [4, 5]. Other tick species may also transmit TBEV under certain ecological conditions ; however it is not known if Ixodes scapularis ticks found in the United States, where TBEV does not occur, are capable of transmitting the virus.
Langat virus (LGTV), a close relative of TBEV, was isolated from Ixodes granulatus ticks in Malaysia . Although the virus is antigenically closely related to TBEV, there are no reports of naturally-acquired cases of human disease caused by LGTV. The attenuated LGTV strain E5 was tested as a candidate live vaccine against TBEV in animals and human volunteers. It resulted in high levels of neutralising antibodies which cross-reacted with TBEV, Powassan virus and Kyasanur Forest disease virus [7, 8]. Due to its close antigenic relationship with TBEV, low pathogenicity and lack of naturally-occurring cases of disease in humans and animals, LGTV is a useful experimental model for more virulent tick-borne flavivirus infections.
Most knowledge of the response of arthropods to microorganisms has been obtained from studies in insects. These have revealed the involvement in the antiviral response of several signaling pathways including RNA interference (RNAi) [9, 10], Toll, Immune deficiency (IMD), and Janus kinase-signal transducers and activators of transcription (JAK/STAT), as well as melanisation, autophagy and possibly heat shock proteins (HSPs) (reviewed by [11–14]). RNAi, Toll, IMD and JAK/STAT pathway components have been identified in the genome of the tick I. scapularis [15, 16], but in comparison to insects there is only limited knowledge on tick innate immune responses to pathogen infection [15, 17–19]. A recent study reported a role for the JAK/STAT pathway in I. scapularis ticks during Anaplasma phagocytophilum infection . This study showed that silencing of STAT or JAK, but not Toll-1, TAK1 or TAB1, which are components of the Toll and IMD pathways, resulted in an increase in A. phagocytophilum in infected ticks and that the JAK/STAT pathway controls bacterial infection by regulating the expression of antimicrobial peptides of the 5.3 kD gene family. Other important regulatory molecules with a possible role in tick innate immune responses include RNA-dependent RNA polymerase, subolesin and ubiquitin-related molecules [21–24].
The only antiviral innate immune response described to date in ticks is RNAi [25, 26]. RNAi has been efficiently used for gene knockdown in ticks and tick cell lines [27–29]. Tick cell lines have been used as tools to understand LGTV and TBEV interactions with their vectors [30–38]. Recently, Dicer (Dcr) and several orthologues of Argonaute (Ago) 2, a key member of the exogenous siRNA pathway in insects, were identified in ticks and Dcr 90, Ago 16 and Ago 30 were shown to mediate an antiviral response .
The present study was carried out with the aim of identifying transcripts and proteins with a possible role in tick innate antiviral responses. We first characterised TBEV infection in the tick cell lines IDE8 derived from the only tick species with a sequenced genome, I. scapularis, and IRE/CTVM19 derived from I. ricinus, a natural vector of TBEV. We then investigated differences in transcript and protein abundance between TBEV-infected and mock-infected tick cells using the Illumina HiSeq2000 platform and LC-MS/MS, respectively. Statistically significantly differentially-expressed transcripts and differentially-represented proteins were identified, annotated and grouped according to their biological function. Finally, using LGTV which could be handled at a lower level of biosafety containment than TBEV, we silenced selected transcripts and proteins by RNAi, to elucidate their effect on virus replication and production.
This study was carried out in strict accordance with the Czech national law and guidelines on the use of experimental animals and protection of animals against cruelty (the Animal Welfare Act Number 246/1992 Coll.). The protocol was approved by the Committee on the Ethics of Animal Experiments of the Institute of Parasitology and of the Departmental Expert Committee for the Approval of Projects of Experiments on Animals of the Czech Academy of Sciences (Permit Number: 165/2010).
Tick and mammalian cell lines
The I. scapularis-derived cell line IDE8  was maintained in ambient air at 32 °C in L-15B medium  supplemented with 10 % tryptose phosphate broth (TPB), 5 % fetal calf serum (FCS), 0.1 % bovine lipoprotein (MP Biomedicals), 2 mM L-glutamine and 100 units/ml penicillin and 100 μg/ml streptomycin (pen/strep). The I. ricinus-derived cell line IRE/CTVM19  was maintained in ambient air at 28 °C in L-15 (Leibovitz) medium supplemented with 10 % TPB, 20 % FCS, 2 mM L-glutamine and pen/strep . Baby hamster kidney (BHK-21) cells (C-13, ATCC, cat: CCL-10) and African green monkey kidney epithelial (Vero) cells (ECACC, cat: 84113001) were grown at 37 °C in a humidified atmosphere of 5 % CO2 in air. Porcine kidney stable (PS) cells were grown at 37 °C in ambient air . BHK-21 cells were maintained in Glasgow Minimal Essential Medium (GMEM) supplemented with 5 % newborn calf serum (NBCS), 10 % TPB, 2 mM L-glutamine and pen/strep (GMEM/5 % NBCS). Vero cells were grown in Dulbecco’s Minimal Essential Medium containing 10 % FCS and pen/strep. PS cells were maintained in L-15 (Leibovitz) medium supplemented with 3 % NBCS, 2 mM L-glutamine, pen/strep and 0.25 μg/ml amphotericin B (L-15/3 % NBCS).
Virus strains, propagation and virus titration
The TBEV strain Neudoerfl was kindly provided by Professor F.X. Heinz, Institute of Virology, Medical University of Vienna, Austria, and had been passaged five times by intracranial infection of suckling mice before use in the present study. Suckling CD1 mice were intracranially infected with 1 μl of TBEV-infected mouse brain suspension corresponding to 100 plaque-forming units (PFU) per mouse, or mock-infected with the same volume of uninfected mouse brain suspension. After the onset of symptoms, 4 to 5 days post infection (p.i.), the TBEV-infected mice were euthanised and the brains removed. The mock-infected mice were euthanised 2 days later, to prevent any possibility of cross contamination while handling the samples. The brains were homogenised in L-15/3 % NBCS to obtain a 20 % mouse brain suspension (w/v) using a Tissue Lyser II (Retsch) at 30 Hz (30/s) for 2 min. The homogenate was then centrifuged for 10 min at 16,000 x g at 4 °C and the clarified supernatant was used for infection of tick cell lines.
The LGTV strain TP-21 was kindly provided by Dr Sonja Best, Laboratory of Virology, Rocky Mountain Laboratories, NIAID, NIH, Hamilton, Montana, USA and was propagated in Vero cells prior to being used for infection. TBEV was titrated on PS cells as described previously [44, 45]; the titre of the stock used in the experiments was 6x107 PFU/ml. LGTV was titrated on BHK-21 cells in 12-well plates using Avicel (RC-581, FMC Biopolymer) as an overlay. In brief, cells were seeded at a density of 1.5x105 cells per well in GMEM/5 % NBCS and incubated overnight. When the cells were 80 % confluent, medium was removed and replaced with supernatant of test samples which had been 10-fold serially diluted in GMEM/2 % NBCS. After incubation on a shaker for 1 h, cells were overlaid with 1 ml of Avicel suspension (1.2 g Avicel in 100 ml PBS) mixed in a 1:1 ratio with 2x Minimal Essential Medium (Gibco) supplemented with 5 % FCS. Cells were incubated for 4 days, fixed in 10 % neutral buffered formaldehyde (Leica), stained with 0.1 % aqueous toluidine blue for 30 min and plaques were counted. The titre of the LGTV stock used in the experiments was 2x106 PFU/ml.
Infection of tick cell lines
For establishing a TBEV growth curve, both tick cell lines were seeded at a density of 5x105 cells per ml in 2 ml total medium volume in flat-sided tubes (Nunc) and 24 h later were infected with TBEV diluted in the respective tick cell growth medium to a multiplicity of infection (MOI) of 5. After infection, supernatant was collected over a 10 day time-period for virus titration. Tick cells used for transcriptomic and proteomic analysis were seeded at a density of 1.5x106 (IDE8) or 1x106 (IRE/CTVM19) cells per ml in flat-sided tubes and were infected 24 h later with TBEV at MOI 5. Preliminary experiments (data not shown) revealed that RNA and protein yields from equivalent numbers of cells were higher from IRE/CTVM19 cells than from IDE8 cells, presumably because IRE/CTVM19 cells are larger than IDE8 cells (authors’ observations), and established that the cell densities used were the minimum required to produce acceptable RNA and protein yields. After incubation for the required time cells were harvested by pipetting for analysis as indicated below. To validate differential expression of transcripts in IDE8 and IRE/CTVM19 cells upon LGTV infection, cells were seeded at a density of 5x105 per ml in flat-sided tubes, mock-infected or infected 24 h later with LGTV at MOI 5 and the medium was changed 2 h p.i.. At 2 and 6 days p.i., cells were harvested and RNA was isolated for qRT-PCR analysis as described below.
Tick cells were seeded at a density of 5x105 cells per ml on 12 mm diameter glass coverslips in 24-well plates, incubated overnight and two replicate wells were infected with TBEV at each of MOI 0.1, 1 and 5. The immunostaining procedure was carried out in the 24-well plates. At day 2 p.i. the cells on the coverslips were washed in PBS, fixed in 10 % neutral buffered formaldehyde for 1 h and washed for 5 min in PBS. The cells were then permeabilised with 300 μl of 0.3 % TritonX-100 for 30 min and subsequently with 0.1 % SDS for 10 min. After permeabilisation, the cells were washed in PBS and then blocked with 300 μl of 1 % bovine serum albumin (BSA) in PBS for 60 min. The blocking solution was removed and the primary antibody, Anti-Flavivirus Group antigen antibody (clone D1-4G2-4-15, Millipore, recognising the E protein of flaviviruses) diluted 1:100 in 1 % BSA, was added and incubated for 2 h at room temperature. The cells were washed three times for 5–10 min in PBS and the secondary antibody, Goat anti-Mouse IgG (H + L) DyLight 488 conjugate (Pierce Thermo Scientific) diluted 1:1000 in 1 % BSA, was added and incubated for 1 h. After three washes of 5 min each in PBS, the coverslips were mounted on glass microscope slides with Vectashield HardSet mounting medium containing DAPI (Vector Laboratories). Images were taken of randomly-selected fields using an Olympus Fluoview FV10 confocal microscope and the percentage of green cells determined by visual counting of DyLight 488-positive and negative cells.
RNA and protein isolation
IDE8 and IRE/CTVM19 cells were seeded at densities of, respectively, 1x106 and 1.5x106 cells per ml into 24 tubes per cell line and, 24 h later, 12 tubes per cell line were infected with TBEV at MOI 5 and 12 tubes were mock-infected with the same volume (300 μl) of diluted uninfected mouse brain suspension. On each of days 2 and 6 p.i., cells from six tubes per treatment were harvested by pipetting and the cell suspension from each replicate tube was split into two aliquots of 1 ml, which were both centrifuged at 500 x g for 5 min and the supernatants discarded. One aliquot from each replicate tube was used for RNA isolation using 1 ml of TriReagent (Sigma) according to the manufacturer’s instructions. To improve RNA purity, RNA samples were further purified using the RNeasy Mini kit (Qiagen) according to the manufacturer’s instructions. RNA samples were stored at −80 °C. The second aliquot from each replicate tube was used for protein isolation as follows. The cell pellet was washed twice with ice-cold PBS and resuspended in 350 μl ice-cold PBS supplemented with 1 % Triton X-100, 50 μl cOmplete, Mini, EDTA-free Protease Inhibitor Cocktail (Roche) and 3.5 μl Halt Phosphatase Inhibitor Cocktail (Thermo Scientific Pierce). After incubation for 1 h on ice, the cell suspension was homogenised at 4 °C using a micro pestle and centrifuged at 200 x g for 5 min to remove cell debris. Supernatants were collected and protein concentration was determined with the BCA protein assay kit (Thermo Scientific Pierce) using BSA as standard. Samples were stored at −80 °C until use.
Protein quality was tested by SDS-PAGE with subsequent Coomassie staining as follows. Protein samples were mixed 1:1 (v/v) with 2x Laemmli buffer (Biorad) supplemented with 5 % β-mercaptoethanol (Sigma), heated at 96 °C for 10 min and then loaded onto discontinuous SDS-PAGE (0.75 mm thick, 4 % stacking and 12 % resolving) gels. The gels were run at 40 V for 30 min followed by 120 V for 30 min. For staining gels, a solution consisting of 0.25 g Coomassie Brilliant Blue R-250 (Thermo Scientific Pierce) dissolved in 40 % water, 50 % methanol and 10 % glacial acetic acid was prepared. The gels containing protein were immersed in staining solution for 3 h prior to de-staining in a solution containing 50 % methanol, 40 % water and 10 % glacial acetic acid. The de-staining solution was changed several times until the protein bands were clearly visible. Samples showing good protein quality with clear, distinct bands and widely-distributed molecular masses were considered suitable for proteomic analysis.
RNA sequencing and assembly
RNA integrity was assessed using the RNA 6000 Nano Kit (Agilent) according to the manufacturer’s instructions and tested for RNA integrity using a 2100 Bioanalyser (Agilent) according to the manufacturer’s instructions. Before sequencing, infection levels were measured by qRT-PCR using primers targeting the TBEV NS5 protein (Additional file 1). Aliquots of only those samples showing satisfactory RNA quality, and presence or absence of TBEV infection in the case of infected cells and mock-infected controls respectively, were pooled according to time-point, cell line and condition. The pooled RNA samples containing total RNA were processed by ARK-Genomics (http://www.ark-genomics.org/) according to the Truseq RNA sample guide 1500813 (Illumina Inc). In brief, mRNA molecules containing poly(A) tails were purified from total RNA using poly-T oligo‐attached magnetic beads. The resulting mRNA was fragmented, first and second strand cDNAs were synthesised, ends repaired and adapters ligated. After PCR amplification, the cDNA library was quantified, multiplexed and sequenced on the HiSeq2000 platform, generating paired end reads of approximately 2 x 100 bp in length. The reads were sorted into samples according to cell line, time-point and treatment using the software CASAVA 1.8 (Illumina, https://support.illumina.com/sequencing/sequencing_software/casava.ilmn). Reads obtained from the I. scapularis-derived cell line IDE8 were mapped with TopHat 2.0.3  against the I. scapularis reference genome (iscapularis.SUPERCONTIGS-Wikel.IscaW1.fa). Counts of reads mapping to the genome were generated with HTSeq count 0.5.3p9 (http://www-huber.embl.de/users/anders/HTSeq/doc/count.html). The unmapped reads were de novo assembled with CLC genomic workbench 5.1 (http://www.clcbio.com/products/clc-genomics-workbench/) and mapped with BWA 0.6.1  against the mapped, filtered (5x 400b) reads for generating counts using a Perl script. The reads obtained from the I. ricinus cell line IRE/CTVM19 were de novo assembled as described for the unmapped reads from IDE8. Only reads mapping unambiguously to contigs were counted.
Differential gene expression analysis and annotation
Each assembled contig was assumed to represent a transcript and, since the majority of reads generated during sequencing mapped unambiguously, it was assumed that the count data reflected the expression of each transcript. As reported in previous studies [48–51], we did not use biological replicates for RNA-seq but used pooled RNA isolated from replicate samples; the algorithm used to quantitate transcriptomics data allows the use of non-replicated samples [52, 53]. Differential gene expression was analysed using DESeq in R following the script for working without replicates . DESeq uses a very conservative approach in calling statistical significance in samples without biological replicates. This results in fewer transcripts being called statistically significant; thus some important transcripts might have been missed, whereas the transcripts that were included were strongly supported. Transcripts that were greater than log2 2-fold differentially expressed, and those statistically significantly differentially expressed, were annotated first using Blast2GO  with a Blastx algorithm against the NCBI nr database using a threshold of E-value < 10−6 as cut-off. Those sequences which did not result in any blast hits with Blast2GO were blasted manually using Blastx and Blastn algorithms against the nr and nt NCBI databases and were included when they showed more than 50 % coverage and more than 25 % sequence similarity. All sequences obtained by either of the two approaches were additionally blasted against the UniProt/Swiss-Prot and VectorBase databases to retrieve ontology information, including ontology information for conserved domains provided by NCBI and UniProt. For the statistically significantly differentially-expressed transcripts, literature research was performed in addition to database information retrieval to assign biological process groups.
For those samples which passed both the RNA and protein quality checks in each experimental group, protein extracts equivalent to 100 μg for each group, obtained by pooling equal aliquots from the replicates, were suspended in 100 μl of Laemmli buffer supplemented with 5 % β-mercaptoethanol and applied to 1.2 cm-wide wells of a conventional discontinuous SDS-PAGE gel (0.75 mm thick, 4 % stacking, and 12 % resolving). The electrophoretic run was stopped as soon as the dye front entered 3 mm into the resolving gel. The whole proteome, concentrated within the stacking/resolving gel interface, was visualised using Bio-Safe Coomassie Stain G-250 (BioRad), excised and cut into cubes of 2 x 2 mm. The gel pieces were destained in a 1:1 mixture of acetonitrile and water and digested overnight at 37 °C with 60 ng/ml sequencing grade trypsin (Promega, Madison, WI) as described previously . Trifluoroacetic acid was added to a final concentration of 1 % to stop digestion, and peptides were desalted onto OMIX Pipette tips C18 (Agilent Technologies, Santa Clara, CA, USA) as described previously , dried down and stored at −20 °C until required for mass spectrometry analysis. The desalted protein digests were resuspended in 0.1 % formic acid and analysed by reversed phase liquid chromatography coupled to mass spectrometry (RP-LC-MS/MS) using an Easy-nLC II system coupled to an ion trap LTQ-Orbitrap-Velos-Pro mass spectrometer (Thermo Scientific, San Jose, CA, USA). The peptides were concentrated (on-line) by reverse phase chromatography using a 0.1 mm × 20 mm C18 RP precolumn (Thermo Scientific), and separated using a 0.075 mm x 100 mm C18 RP column (Thermo Scientific) operating at 0.3 μl/min. Peptides were eluted using a 180-min gradient from 5 % to 40 % solvent B in solvent A (Solvent A: 0.1 % formic acid in water, solvent B: 0.1 % formic acid, 80 % acetonitrile in water). ESI ionisation was carried out using a nano-bore emitters stainless steel ID 30 μm (Thermo Scientific) interface. Peptides were detected in survey scans from 400 to 1600 atomic mass units (amu, 1 μscan), followed by fifteen data-dependent MS/MS scans (Top 15), using an isolation width of 2 mass-to-charge ratio units, normalised collision energy of 35 %, and dynamic exclusion applied during 30 s periods.
Proteomic data analysis and annotation
Peptide identification from the MS/MS raw data was carried out using the SEQUEST algorithm (Proteome Discoverer 1.3, Thermo Scientific). Database searches were performed against UniProt-Arthropoda.fasta and UniProt-Flaviviridae.fasta. The following constraints were used for the searches: tryptic cleavage after Arg and Lys, up to two missed cleavage sites, and tolerances of 10 ppm for precursor ions and 0.8 Da for MS/MS fragment ions, and the searches were performed allowing optional methionine oxidation and cysteine carbamidomethylation. Searches were performed against a decoy database in an integrated decoy approach. A false discovery rate (FDR) < 0.01 was considered as a condition for successful peptide assignments and at least 2 peptides per protein in at least one of the samples analysed was the condition for subsequent protein identification (Additional file 2). The total number of peptide assignments for each protein were normalised against the total number of peptide assignments in each sample (control and infected tick cell lines at days 2 and 6 p.i.) and differential representation of individual proteins between different samples was determined using Chi-square test statistics with Bonferroni correction in the IDEG6 software (http://telethon.bio.unipd.it/bioinfo/IDEG6 form/) (p < 0.05) as published . Samples with a p-value equal to or lower than 0.05 were considered statistically significant. Statistically significantly differentially-represented proteins were allocated to biological process groups using ontology information available on the UniProt/Swiss-Prot and Panther databases, including information for conserved domains. Information was curated manually through literature search.
For verification of infection status of TBEV-infected and mock-infected cells, 1 μg of total RNA was reverse-transcribed using Superscript III (Invitrogen) and random hexamers according to the manufacturer’s instructions. For verification of RNA-Seq data and knockdowns followed by LGTV infection, 1 μg of total RNA was reverse-transcribed using the High-Capacity cDNA Reverse Transcription kit (Applied Biosystems) according to the manufacturer’s instructions.
Verification of TBEV infection and RNA-Seq data by qRT-PCR
The gene coding for TBEV NS5 protein was cloned into the pJET vector using the CloneJET PCR Cloning Kit (Fermentas) according to the manufacturer’s instructions. In brief, the plasmid pTND/ΔME  was linearised, purified and used as DNA template for amplification of TBEV NS5 using KOD polymerase (Novagen). An aliquot of the 187 bp PCR product was visualised by gel electrophoresis and, since only one band with the correct size was visible, the non-purified product was directly used for ligation. For ligation, 10 μl 2x reaction buffer, 2 μl non-purified PCR product, 1 μl pJET1.2 blunt cloning vector, 6 μl nuclease-free water and 1 μl T4 DNA ligase were mixed by vortexing. The ligation mixture was incubated for 5 min at room temperature before using directly for transformation of DH5α. To check if the correct insert was cloned into the vector, the plasmid was linearised and sent for sequencing to GATC Biotech (London, UK).
For verification of RNA-Seq data, primers for 12 transcripts (Additional file 1) were designed, using as template species-specific sequences or identical regions from sequences common to both I. scapularis and I. ricinus obtained by HiSeq2000. The same samples from which aliquots had been pooled for the transcriptome profiling were used individually for qRT-PCR analysis. Beta actin and ribosomal protein L13A were used as housekeeping genes for normalisation. Primer efficiencies were calculated for each primer and the quantity of gene transcripts in infected samples relative to controls was calculated using the 2-∆∆CT method [59, 60].
dsRNA production and silencing of tick transcripts
Long dsRNA transcripts (407 - 615 bp in length) specific to genes from each cell line were produced from PCR products flanked by T7 promoter sequences using the MegaScript RNAi kit (Ambion). In brief, cDNA generated by reverse transcription from total RNA of tick cells or from the plasmid pIRES2-eGFP (Clontech) was used as template to generate specific PCR products using T7 primers (Additional file 1) by PCR. PCR products were gel-purified and sequenced. The gel-purified PCR products were subjected to an additional PCR amplification and were then transcribed using the MegaScript kit according to the manufacturer’s instructions.
For knockdown experiments, LGTV was used at a low MOI to ensure that not all cells would be infected initially, thereby allowing virus to spread from cell to cell which might enhance detection of any effect of transcript knockdown on virus replication and/or production. Cell lines IDE8 and IRE/CTVM19 were seeded at a density of 5x105 cells per ml in 24-well plates and were incubated in humidified self-sealing polythene bags. For IDE8 cells, 300 ng of dsRNA was added to the supernatant twice, at 8 h and 48 h post-seeding. Approximately 72 h post-seeding, cells were infected with LGTV at MOI 0.01; 48 h later supernatant was collected for plaque assay and cells were harvested for RNA extraction. To achieve a good knockdown in IRE/CTVM19 cells, cells were transfected 24 h post-seeding with 400 ng of dsRNA using Lipofectamine® 2000 (Invitrogen) as previously described  and, after incubation for a further 48 h, were infected with LGTV at MOI 0.5. At 24 h p.i supernatant was collected for plaque assay and RNA was extracted using TriReagent as above. Non-specific dsRNA encoding eGFP was used as a negative control, to provide a baseline level of activation above which the effect of the specific exogenous dsRNA was measured. Additional controls, in which samples were not treated with dsRNA prior to infection, were included to test whether or not addition of any non-specific dsRNA triggers an innate immune response in tick cells and to provide a baseline for virus replication and virus titres in untreated cells. For the detection of Ago and Dcr knockdowns, PCR was carried out (95 °C for 2 min, 95 °C for 30 s, primer set specific annealing temperature (Additional file 1) for 30s, 72 °C for 50 s, 72 °C for 7 min) using GoTaq reaction mix (Promega), according to the manufacturer’s instructions, together with 2 μl of the cDNA reaction and the corresponding primers (Additional file 1).
PCR products were run on a 1 % agarose gel and gel images were taken and used to quantify mRNA knockdown with Image Lab software (BioRad) normalised to beta actin. Gene knockdowns and LGTV RNA levels were measured by qRT-PCR with, respectively, target gene-specific primers or LGTV NS5-specific primers (Additional file 1), using FastStart Universal SYBR Green Master (Rox) (Roche) with a temperature profile of 95 °C for 10 min, 95 °C for 15 or 20s, specific annealing temperature for 20 or 30s, 72 °C for 15 or 30s and 95 °C for 15 s. All qRT-PCR reactions were followed by melting curve generation (60-95 °C) to confirm amplification specificity. Primer efficiencies were calculated for each primer and the quantity of gene transcripts in infected samples relative to controls was calculated using the 2-∆∆CT method [59, 60].
Statistical analysis of gene knockdown experiments
Gene knockdowns were done in quadruplicate in three independent experiments. Only those samples in which a knockdown occurred were included in subsequent statistical analysis. Analysis was done in GraphPad Prism (http://www.graphpad.com/scientific-software/prism/). Statistical significance of the three independent experiments was analysed using the two-way Analysis of Variance Fisher’s LSD test (P = 0.05).
Results and discussion
Characterisation of TBEV growth in tick cells
To measure newly-produced virus within defined time-periods, tick cells were infected with TBEV at MOI 5. The cells were then washed and fresh medium was added at 2 h p.i. for time-points 12, 18 and 24 h, and at 24 h prior to each sampling for the subsequent daily time-points up to day 10 p.i. (Fig. 1b). Supernatants were collected at each time-point and the virus titre measured by plaque assay on PS cells. The pattern of TBEV infection was similar in both cell lines with the highest level of virus production between 2 and 4 days p.i. (Fig. 1b). Interestingly, the maximum titre in the cell line IRE/CTVM19 derived from a known TBEV vector was approximately 2 log10 higher than in the IDE8 cell line. The higher level of virus production in the I. ricinus cell line compared to the cell line derived from I. scapularis, which is not known to be a natural vector of TBEV, confirms previous studies [31, 34]. The lower virus titres in IDE8 cells might be an indicator of reduced vector capacity of I. scapularis for this virus, which could be due to cells being less efficiently infectable or having a more rapid and efficient antiviral response.
In order to examine how these two cell lines react to virus infection, two time-points were chosen for transcriptomic and proteomic analysis: one early in infection at day 2 when virus production was increasing, and one later in infection at day 6 when virus production was decreasing.
Generation of samples for transcriptomic and proteomic analysis
Six replicate tubes per time-point per cell line were either infected with TBEV at MOI 5 or mock-infected. At days 2 and 6 p.i., both RNA and protein were isolated from each tube by dividing the cell suspension in half, thus ensuring that both transcriptomic and proteomic analyses were carried out on samples derived from each tube. Because of the biosafety restrictions on handling TBEV, it was only possible to generate a single set of samples, which yielded sufficient material for a single sequencing run and proteomic analysis for each treatment. Therefore the results will be considered in this context as a baseline for future studies.
Prior to RNA sequencing, RNA samples were tested for the presence of TBEV by qRT-PCR (Additional file 3). Only mock-infected samples negative for TBEV with NS5 RNA levels below the detection limit of the assay, and infected samples positive for TBEV with NS5 levels above 10,000 copies, were included in the subsequent analysis. Furthermore, only high quality samples showing no signs of RNA degradation were used. The soluble proteins extracted from mock-infected and TBEV-infected IDE8 and IRE/CTVM19 cells at days 2 and 6 p.i. were tested for protein quality prior to pooling. Only those RNA samples and protein samples which passed both the RNA and protein quality checks (Additional file 3) were pooled, guaranteeing that pooled RNA samples and protein samples contained RNA and protein, respectively, from the same original samples.
Illumina sequencing, assembly and differential gene expression analysis
Sequencing depth, assembly of RNA-Seq data and total number of proteins identified by MS from TBEV-infected and mock-infected (control) IDE8 and IRE/CTVM19 cells on days 2 (2d) and 6 (6d) p.i
Mean reads per lane
Total no. of contigs assembled
Mean contig length
Total no. of proteins identified by RP-LC-MS/MS
IDE8 control 2d
3.71E + 07
IDE8 infected 2d
3.25E + 07
IDE8 control 6d
3.18E + 07
IDE8 infected 6d
2.96E + 07
IRE/CTVM19 control 2d
3.03E + 07
IRE/CTVM19 infected 2d
2.70E + 07
IRE/CTVM19 control 6d
4.44E + 07
IRE/CTVM19 infected 6d
2.61E + 07
Number of statistically significantly differentially expressed transcripts that were up- or down-regulated upon TBEV infection of IDE8 and IRE/CTVM19 cells on days 2 (2d) and 6 (6d) p.i
Protein identification and differential protein representation
Number of statistically significantly differentially-represented proteins upon TBEV infection of IDE8 and IRE/CTVM19 cells on days 2 (2d) and 6 (6d) p.i
Annotation and ontology of tick cell transcripts and proteins
In IDE8 cells, all 42 statistically significantly differentially-expressed tick cell transcripts (Fig. 2a) were annotated; 32 were most closely related to transcripts from I. scapularis, 2 to transcripts from other Ixodes spp, and 8 to transcripts from rodent species. In IRE/CTVM19 cells (Fig. 2b), only 56 of the 81 statistically significantly differentially-expressed tick cell transcripts could be annotated; 54 corresponded to I. scapularis and one each to Rattus norvegicus and Harpegnathos saltator. The other 25 transcripts did not return any blast hits and were excluded from further analysis. This lack of homology has been reported in other tick studies [66–68] and might be attributed to factors such as low sequence and/or assembly quality, lack of homology to the I. scapularis genome due to its fragmented state or that these sequences represent novel species-specific transcripts .
With respect to the statistically significantly differentially-represented proteins, the majority in both IDE8 (Fig. 2c) and IRE/CTVM19 (Fig. 2d) corresponded to I. scapularis followed by other tick species including Amblyomma spp., Hyalomma marginatum rufipes and Haemaphysalis qinghaiensis.
The I. scapularis genome is currently not fully annotated, and annotation of transcripts or proteins, including the inference of their functional role, relies in the majority of cases on sequence similarity to evolutionarily quite distant species, including mammals and insects with well-annotated genomes. Thus although homology is observed for a specific transcript or protein, it might have evolved functions within the tick different from those within other species. This comparative approach might therefore be misleading and makes it difficult to infer true biological role by sequence similarity, conserved domains and/or literature search. Currently the only method to identify possible target genes within large datasets of ticks is to infer biological function from other better-annotated organisms or from sequence similarity to other model or non-model organisms.
In both cell lines, high levels of virus (Fig. 1b) were associated with up-regulation of transcripts annotated as immunity or metabolism and down-regulation of transcripts annotated as cell stress and protein folding. At the protein level, proteins annotated as cell stress or protein folding were, as with transcripts, amongst those consistently down-regulated. Equivalent biological process groups, as well as specific transcripts and proteins, have been observed to be differentially expressed in studies of mosquitoes and mosquito cells upon virus infection [65, 69, 70]. Although some of these studies showed different directions of representation, this might be attributed to different sampling times, different host and/or virus species and/or in vivo versus cell line usage. When comparing the transcript and protein profiles for each cell line at each time-point individually, there is little correlation at the biological process group level (Fig. 3a and b compared to Fig. 4a and b) between transcripts and proteins, apart from those involved in protein folding and cell stress, both of which are generally down-regulated in both TBEV-infected cell lines at both time-points. This lack of correlation between statistically significantly differentially-expressed transcripts and differentially-represented proteins was also observed in studies on tick cell responses to infection with intracellular bacteria [67, 71, 72], and probably reflects the different half-lives of mRNA and proteins and differential regulation of systems at the transcriptional, post-transcriptional, translational or post-translational levels. Novel approaches could increase the correlation between these datasets, for example proteomics informed by transcriptomics  which has recently been applied to ticks in vivo [51, 74].
To validate the differential gene expression observed during RNA-Seq, twelve transcripts differentially expressed in the transcriptomics data, and/or coding for differentially represented proteins in the proteomics data, were selected for qRT-PCR analysis. Preference was given to those putatively involved in immunity or cell stress, and transcripts/proteins with a range of different expression levels were chosen. Transcripts for the housekeeping genes ribosomal protein L13A and beta actin were used for normalisation since neither was differentially expressed in either of the cell lines at either time-point in the transcriptomic data. Most of the twelve transcripts showed similar fold changes by qRT-PCR and RNA-Seq, or at least confirmed the trend seen in the sequencing data (Additional file 4). The differences observed between the two techniques qRT-PCR and RNA-Seq for complement factor H in both cell lines and coagulation factor in IDE8 cells at day 6 p.i. (Additional file 4, A) might result from splice variants of these transcripts, different RNA processing techniques, primer design, the choice of reference genes and/or the different normalisation methods – normalisation across the whole transcriptome for sequencing depth with RNA-Seq versus normalisation against specific transcripts with qRT-PCR . Similar observations, in which the differential gene expression data by RNA-Seq was, for most transcripts, in good agreement with qRT-PCR data in showing at least a similar trend in differential expression, have also been reported in other transcriptomic studies [77–79].
Functional role of selected cell stress and immunity genes and proteins in tick cells during LGTV infection
Selected transcripts were then silenced using sequence-specific dsRNA in IDE8 and/or IRE/CTVM19 cells, the cells were then infected with LGTV and virus replication and infectious virus production were measured by qRT-PCR and plaque assay respectively. Genes possibly involved in immunity such as those encoding complement factor H or trypsin or those possibly involved in cell stress such as those encoding calreticulin, HSP90, gp96 and HSP70 were silenced in both tick cell lines. Transcripts encoding Argonaute (Ago 30) and Dicer (Dcr 90) were included as positive controls , while cells treated with non-specific dsRNA against eGFP were used as baseline controls. Three independent experiments with quadruplicate samples were conducted per cell line; thus 12 samples were analysed in total per cell line. Only those replicates in which silencing was confirmed were included in the statistical analysis.
RNAi is probably the most important antiviral response in insects . The importance of the RNAi response in the antiviral defence response in tick cells was recently confirmed by detecting specific viRNAs in tick cells infected with LGTV and observing a proviral effect upon silencing of orthologues of key members of the RNAi pathway (Ago 30 and Dcr 90) . Although Ago 30 and Dcr 90 were not differentially expressed upon TBEV infection in the present study, they were included as positive controls. Both proteins are known to be involved in RNAi  and knockdown would be expected to result in an increase in levels of viral RNA as well as infectious virus. In IDE8 cells, silencing of Dcr 90 resulted in a significant increase in LGTV RNA levels and silencing of Ago 30 resulted in a significant increase in levels of both LGTV RNA and infectious virus (Fig. 7d and e). In IRE/CTVM19 cells, significant increases in both LGTV RNA levels (Fig. 8d) and infectious virus (Fig. 8e) were observed following knockdown of both Ago 30 and Dcr 90. This confirms the role of RNAi as an antiviral response in tick cells.
Knockdown of ISCW022021, annotated as complement factor H, resulted in an increase in LGTV RNA and infectious virus in IDE8 cells (Fig. 7) but not in IRE/CTVM19 cells (Fig. 8). Complement Factor H, up-regulated on day 6 p.i. in IRE/CTVM19 cells (Fig. 6), functions in vertebrates as a negative regulator of the alternative pathway of the complement system and as a pattern recognition molecule binding with high efficiency to host-specific molecular signatures, such as heparin and sialic acid, thereby protecting uninfected cells from the complement system . In vertebrates, the complement system is an important innate immune response against different families of viruses [85–87]. However some viruses, such as West Nile virus (WNV), are able to evade the complement system by binding and recruiting Factor H, resulting in decreased complement activation and reduced targeting of virus-infected cells . Although ticks have been shown to have a complement system, with all α2 -macroglobulin family proteins, insect thioester-containing and macroglobulin-related proteins [89, 90], that functions against different types of bacteria [89, 91], nothing is known about the antiviral effect of the tick complement system.
If complement factor H interacts with LGTV in tick cells similarly to WNV in mammalian cells, a decrease in virus replication and/or production would be expected but the opposite was seen in the present study. The increase in viral RNA levels and infectious virus in tick cells upon silencing of complement factor H might have resulted in exhaustion of the complement system and lack of any complement antiviral-activity. The fact that complement factor H was up-regulated in response to virus infection suggests an antiviral role, supported by the observed increase in viral RNA levels and infectious virus upon silencing of complement factor H. The designation of ISCW022021 as complement Factor H by the Ixodes scapularis Genome Project  might be misleading and, instead of being involved in the complement system, it might be involved in another antiviral innate immune response. Further sequence and functional analyses are required to determine the role of ISCW022021 within the tick cell innate immune response.
Silencing of HSP90 and HSP70 resulted in an increase in LGTV RNA levels in IDE8 cells (Fig. 7). This suggests that HSP90 and HSP70 might be involved in loading of siRNA duplexes into Ago 2, as observed in Drosophila ; thus knockdown of either protein would lead to an impairment of RNAi, which would result in reduction of degradation of viral RNA, as suggested by the higher viral RNA levels seen in the present study in cultures in which HSP90 and HSP70 were silenced, compared to unsilenced controls. It would be interesting to test whether simultaneous knockdown of both HSP70 and HSP90 would augment the increase in viral RNA levels. Knockdown of trypsin also resulted in a significant increase in LGTV production accompanied by a slight non-significant increase in viral RNA levels. The putative antiviral effect of trypsin might be due to its serine protease activity, since serine proteases are involved in the modulation of several immune signaling pathways [94–96] and, of these, one or more might mediate an antiviral role in tick cells.
In contrast to the results with IDE8 cells, in IRE/CTVM19 cells only silencing of Ago 30 and Dcr 90 resulted in significantly increased virus replication and production (Fig. 8). Silencing of HSP90 and gp96 in IRE/CTVM19 cells resulted in a significant increase in virus production without affecting LGTV RNA levels, suggesting an antiviral role for these proteins at the post-transcriptional level in this cell line. HSP90 and gp96 are both heat-shock proteins which are involved in folding of different client proteins. Inhibition of HSP90 in mammalian cells has been shown to block viral replication [97, 98] and this protein has been proposed to be an important factor in the replication of a wide spectrum of RNA viruses . In the present study in tick cells, however, HSP90 seemed to be involved in the antiviral response with an inhibitory influence on virus RNA levels in IDE8 and at the post-translational level in IRE/CTVM19. The ER-based heat-shock protein gp96 is important for the folding of Toll-like receptors (TLRs) and integrins in mammals and Drosophila . The putative antiviral role of gp96 observed in the present study might be due to its capacity for folding TLRs or other client proteins involved in the antiviral response, which upon silencing would cause an increase in virus production.
Silencing of complement factor H, which resulted in increased LGTV replication and production in IDE8 cells, did not show any effect in IRE/CTVM19 cells. This could be due to the less efficient and more variable silencing in the latter cell line, compared to IDE8 cells. Additionally the different responses of the two cell lines to LGTV infection could represent a cell line-specific response towards flavivirus infection; however most of the transcripts tested in silencing experiments were differentially expressed upon TBEV infection in IRE/CTVM19 but not in IDE8. The different responses could indicate a species-specific response since the two cell lines were derived from different tick species, or could be due to the heterogeneity of the cell lines  or presence of endogenous viruses. Both IDE8 and IRE/CTVM19 cells are persistently infected with endogenous viruses, St Croix River virus and unidentified reovirus-like particles respectively [100, 101], which could affect the innate immune response towards infection with another virus. The presence of an endogenous virus could either suppress or persistently activate certain immune responses thereby affecting silencing of genes and the effect on virus replication and production. Furthermore, each cell line might have a different timing in the response to virus infection, with IDE8 cells possibly activating a response faster than IRE/CTVM19, which could explain the higher virus titres observed for IRE/CTVM19 in the TBEV growth curve experiment.
Interestingly, LGTV production in samples of both cell lines not treated with dsRNA prior to LGTV infection was significantly higher in comparison to samples treated with control dsRNA against eGFP, suggesting that dsRNA treatment alone triggers an antiviral immune response. This is in contrast to studies on mosquitoes and Drosophila in which an antiviral response, possibly RNAi, was shown to be triggered by virus-specific dsRNA but not by non-specific dsRNA [102, 103]. A possible explanation for this difference could be the presence of RNA-dependent RNA polymerase in ticks  that could be involved in boosting the non-specific antiviral immune response seen in tick cells against non-specific dsRNA. However, studies in other arthropod systems including sandfly cells , shrimp  and honey bees  also showed that non-specific dsRNA can trigger an antiviral state affecting virus infection. Interestingly, dsRNA encoding eGFP resulted in an increase in Dcr 2 levels in Bombyx mori , suggesting that Dcr 2 recognises dsRNA as a pathogen-associated molecular pattern, which might result in the expression of Vago, an interferon-like molecule inducing an antiviral state in neighbouring cells, as described in Drosophila and mosquitoes [108, 109]. Possible homologues of insect Dicers were recently identified in the I. scapularis genome and phylogenetic analysis suggests that Dcr 89 is a possible homologue of insect Dcr 2, whereas Dcr 90 is a possible homologue of Dcr 1 . If Dcr 2 and Vago are also induced in tick cells upon addition of non-specific dsRNA, thereby causing an antiviral state in neighbouring uninfected cells and resulting in reduced virus infection levels in the culture overall as observed in the present study, this could suggest a non-specific antiviral response which recognises dsRNA as foreign in tick cells. Vago has been shown to be present within the I. scapularis genome ; it would therefore be interesting to investigate whether it has a similar function in ticks as it has in insects.
The observation that, in both tick cell lines, silencing of Ago 30, Dcr 90 and HSP90 resulted in increased LGTV NS5 expression and/or virus production strengthens the hypothesis that these proteins are involved in the antiviral response in tick cells. Additionally, the finding that silencing of complement factor H, trypsin and HSP70 in IDE8 and gp96 in IRE/CTVM19 resulted in a proviral effect is encouraging and warrants further experiments to elucidate their roles in the antiviral response in tick cells.
Other transcripts and proteins that may be involved in innate immunity
In addition to complement factor H and trypsin, which were grouped into the biological process group innate immunity and were functionally analysed in the previous section, several other transcripts and proteins with a possible role in the innate immune response of tick cells to virus infection were differentially expressed/represented. Those most likely to be involved in innate immunity are discussed below; others are discussed in Additional file 6.
Of the differentially-expressed transcripts in IDE8 cells, two were inferred by GO ontology descriptors to be involved in innate immunity: the class B scavenger receptor CD36 and the cysteine protease longipain (Fig. 6). CD36 is a surface receptor on tick haemocytes that is up-regulated upon bacterial infection in H. longicornis, where it is involved in granulocyte-mediated phagocytosis of Escherichia coli . CD36 has also been suggested to be involved in the RNAi pathway . RNAi is currently the only antiviral pathway known to be effective in ticks . Up-regulation of CD36 in IDE8 cells upon TBEV infection is therefore intriguing and could be an indicator of up-regulation of the RNAi response. Longipain is present in and on the surface of lysosomes in the midgut of the tick Haemaphysalis longicornis and, in addition to playing a role in blood digestion, is involved in dose-dependent killing of Babesia parasites . Longipain has a high homology to cathepsin B , which in mosquitoes is indirectly involved in defence responses against viruses by triggering apoptosis as observed during dengue virus (DENV) infection . It is not known whether longipain has a similar function in ticks.
In IRE/CTVM19 cells, more transcripts with a potential link to an antiviral response were up-regulated than in IDE8 (Fig. 6). One of these, peroxinectin, up-regulated on day 6 p.i. (Fig. 6), is a cell adhesive peroxidase which is stored in haemocyte granules in crustaceans. Upon an immune stimulus, peroxinectin is released from cells by degranulation and is activated by serine proteases to stimulate cell adhesion, encapsulation, phagocytosis and peroxidase activity [114–117]. In the mud crab, white spot syndrome virus (WSSV) infection results in increased expression of peroxinectin within the first 48 h; this increase is associated with a latent period of WSSV infection, suggesting that peroxinectin is involved in the early defence response against this virus . Furthermore, peroxinectin in crustaceans is associated with the prophenoloxidase (proPO) system, as both require the same activating enzyme, a trypsin-like serine protease [119, 120]. However in contrast to other arthropods, ticks are assumed to lack the proPO system, since no proPO-related gene has so far been identified [15, 18], although controversial reports of the existence of this innate immune response in ticks exist [121–123]. Interestingly, the serine protease trypsin was also up-regulated on day 6 in IRE/CTVM19 cells (Fig. 6).
At the protein level, in IDE8 on day 2 and IRE/CTVM19 on day 6, the 4SNc-Tudor (Tudor-SN) domain protein was under-represented. Tudor-SN is part of the RISC complex in Drosophila, Caenorhabditis elegans and mammals  and is involved in binding and possibly cleaving hyper-edited dsRNAs and miRNAs in Xenopus laevis and humans respectively [125, 126]. Tudor-SN is expressed, and suspected to be part of the RISC complex, in ticks ; it plays a role in tick feeding and the RNAi pathway but does not appear to be involved in innate responses to TBEV or LGTV .
In IRE/CTVM19 at day 6 p.i., Cniwi, a PIWI 1 protein, was under-represented. PIWI proteins are part of the piRNA pathway which was initially thought to protect germline cells from transposable elements in Drosophila. However, virus-specific piRNA molecules were found to be expressed for a variety of different viruses in Drosophila  and mosquitoes [129–132], suggesting a possible antiviral role for the piRNA pathway. In a recent study, the antiviral role of PIWI proteins was confirmed when knockdown of PIWI proteins resulted in an increase in Semliki Forest virus replication and production in mosquito cells . It would be tempting to speculate that this pathway might also be important for the antiviral response against TBEV in tick cells. Under-representation of these two proteins in infected cells would be consistent with suppression of RNAi or perhaps customisation of the RNAi system to virus infection.
Other transcripts and proteins that may be involved in cell stress
In addition to HSP70, HSP90, gp96 and calreticulin, which were grouped into the biological process group cell stress and were functionally analysed above, several other transcripts and proteins involved in cell stress (Fig. 6) might also have an important role in the response of tick cells to virus infection. Those most likely to be involved in cell stress are discussed below; others are discussed in Additional file 6.
Heat-shock proteins are the most abundant and ubiquitous soluble proteins in all forms of life and are involved in a multitude of housekeeping functions essential for cell survival . Studies on tick cell responses to bacterial infection revealed pathogen- and species-specific differences in the expression of HSPs [71, 134]. HSP20 and the small HSP alpha-crystallin B chain were under-represented in both IDE8 and IRE/CTVM19 cells on day 2 p.i. at the transcript level (Fig. 6). In vertebrate cells the generation of large amounts of viral proteins leads, through the unfolded protein response of the endoplasmic reticulum (ER), to cell stress and an increase in HSPs [135, 136]. However, there is evidence that in vertebrates some HSPs may be controlled to disrupt virus replication.
Other proteins involved in cell stress include the three ER chaperones calreticulin, calnexin and gp96. In mammalian cells, calreticulin, a soluble lectin-like chaperone, is involved in Ca2+ homeostasis and is important for the processing and maturation of viral glycoproteins [137, 138]. Knockdown of calreticulin in Vero cells reduces the yield of infectious DENV particles . Similarly, knockdown of calreticulin in Babesia bigemina-infected Rhipicephalus microplus ticks reduced the level of this protozoan parasite . Calnexin, a membrane-bound chaperone of the ER similar to calreticulin, that was up-regulated at the transcript level in IRE/CTVM19 at day 2 p.i., has been shown to be important for viral glycoprotein processing and maturation . In Vero cells, both calnexin and calreticulin have been shown to be important for the production of infectious DENV particles by interacting with the glycosylated DENV E protein, facilitating proper folding and assembly of DENV proteins . Under-representation of calnexin in IRE/CTVM19 cells at day 6 could be interpreted as an antiviral response selected to curtail virus production.
The tumour-rejection antigen gp96 which was up-regulated in IRE/CTVM19 cells on day 2 p.i. and under-represented at the protein level on day 6 p.i. (Fig. 6), is important in mammalian cells for chaperoning TLRs and integrins. Up-regulation of ER chaperones such as gp96 and calreticulin upon virus infection could be a sign of ER stress, which in mammals can lead to triggering of apoptosis or the unfolded protein response leading to inhibition of translation or apoptosis . There is currently no published information on translational inhibition or regulation of the unfolded protein response in ticks.
Transcripts and proteins that may be involved in nucleic acid processing
Several transcripts and proteins involved in nucleic acid processing functions, such as replication, transcription, processing of nucleic acid or translation, were differently expressed upon TBEV infection in tick cells. This is not surprising since viruses require the nucleic acid processing machinery of the host to amplify their genome and many viruses perturb these processes in cells or manipulate them for their own advantage. Differential regulation of this group of transcripts and proteins was also observed in several other transcriptomic and proteomic studies of arthropods upon virus infection (e.g. [65, 66, 70, 83, 142]). Many of these transcripts and proteins might be involved in replication and translation of TBEV in tick cells and might be interesting targets for future research to understand virus infections in tick cells and ticks. Histones and elongation factor (EF)-1 alpha were differentially represented at both the transcript and protein levels. Several viral proteins have been shown to target histone proteins and host chromatin to interfere with host gene expression by various mechanisms and for different purposes . The C protein of DENV for example targets core histones during infection to disrupt the host cell genetic machinery in favour of viral replication . EF-1 alpha was up-regulated and shown to be important for virus replication in mammalian and mosquito cells during DENV and WNV infection [145–147]. Interestingly, the eukaryotic translation initiation factor eIF3 was under-represented at the protein level in IDE8 cells at day 6 p.i. and in IRE/CTVM19 cells at day 2 p.i. as observed during WNV infection in Vero cells . This is surprising since eIF3 together with the 40S ribosomal subunit, also under-represented in the present study, have been shown to be important during the initial phase of protein synthesis, and flaviviruses are thought to prevent host cell protein shutoff, at least in mammalian cells . However, a recent study using yellow fever virus (YFV) in mammalian cells found that NS5 interacts with eIF3L, a subunit of eIF3, and that overexpression of this subunit facilitates YFV translation but does not affect global protein synthesis . This suggests that eIF3 is also important in tick cells for replication of flaviviruses but that down-regulation of this initiation factor might have an antiviral effect. Furthermore EF-2, t-RNA synthetasen and several other participants in the translation of RNA were also under-represented in the present study.
Another interesting observation is that in both cell lines at both time-points DEAD-box RNA helicase was up-regulated at the protein level which was also seen at the transcriptional level for DENV in Aedes aegypti cells . This is interesting, since Dcr 2, a DExD/H-box helicase, was shown to be capable of sensing viral dsRNA in Drosophila leading to the production of possibly antiviral molecules .
This is the first study that combines transcriptomic and proteomic analysis to investigate the response of tick cells to infection with a medically important virus. Despite the limitations imposed on the study by biosafety considerations, the findings represent a valuable baseline for future research. Tick cells responded to TBEV infection by changing the expression and/or representation of cellular genes and proteins involved in a variety of biological processes, including metabolism, transport, protein folding, nucleic acid processing, signaling, cell stress and immunity, revealing a complex response of tick cells to virus infection at both the transcriptome and proteome levels, as observed in other arthropods upon virus infection [65, 66, 69, 70, 83, 142, 148, 151–153]. Some of these transcripts and/or proteins, such as those involved in nucleic acid processing, transport, metabolism, protein folding and cell stress, have also been identified in other species as important host cell factors exploited by viruses to support their life cycle in processes including endocytosis, trafficking, viral RNA transcription and translation and virus maturation. Further analyses of these transcripts/proteins using techniques such as Western blotting, although limited by the lack of specific antibodies to tick proteins, might reveal specific factors required for successful infection, replication and production of viruses in ticks and tick cell lines. The dataset created in the present study represents an important starting point for elucidating the viral life cycle and virus-vector relationships.
Several of the identified transcripts and/or proteins have a possible role in immune-related pathways such as the ubiquitin-proteasome pathway, phagocytosis, the complement system, the piRNA pathway and the unfolded protein response. It is obvious from the present study that RNAi is not the only mechanism involved in the antiviral defence response of ticks and that further research is required to elucidate the cellular mechanisms behind virus infection in tick cell lines and ticks. While some of the changes observed upon virus infection could be a response to any microorganism or to cell stress, others appear to be more specific to virus infection, including a Cniwi protein and the 4SNc Tudor domain protein, which are components of the RNAi pathway. Therefore tick cells seem to be able to respond differently to viral and bacterial infection. This is not surprising since tick cells have been shown to be able to raise a specific response against certain bacterial infections [89, 91].
An unexpected but interesting observation was the down-regulation/under-representation of the heat-shock proteins upon TBEV infection in both tick cell lines. This is surprising since vertebrate and invertebrate HSPs are usually up-regulated upon virus infection and some viruses exploit the presence of HSPs to support virus infection [135, 136]. However, HSPs might also have an antiviral effect due to their implication in loading siRNAs into the RISC complex in Drosophila. Thus the down-regulation of HSPs might either be a cellular innate immune response protecting the cell by possibly preventing correct folding of the large number of viral proteins, thereby reducing the production of viral particles, or a response induced by the virus to prevent an efficient RNAi response. Knockdown of the RNAi components Ago 30, Dcr 90, the putative complement component complement factor H, the serine protease trypsin and the HSPs HSP90, HSP70 and gp96, suggests a role for these in the antiviral response of tick cells.
Overall the two cell lines showed a complex expression pattern upon TBEV infection, with differences in expression when compared to each other at both the transcript and protein levels (Fig. 6). These differences could be caused by the necessity of using different methods used to assemble the transcriptomic data, or could represent cell line-specific responses to TBEV infection. One cell line might react more slowly in response to virus infection than the other. Alternatively the response might be species-specific, since the two cell lines are derived from different tick species. Furthermore, both cell lines are heterogeneous, with a range of different cell types present within each culture which could be responding differently to virus infection.
This study enhances the understanding of viral infection of tick cells by identifying transcripts and proteins which may have a role in the innate antiviral defence response of ticks by augmenting or limiting virus production. This preliminary knowledge can be used in future studies to identify important host cell factors required for viral infection, as well as elucidating the innate immune response of tick cells to virus infection.
The authors would like to thank the Tick Cell Biobank at The Pirbright Institute (http://www.pirbright.ac.uk/research/tickcell/Default.aspx) and Dr Ulrike Munderloh of the University of Minnesota for provision of tick cell lines, Dr Sonja Best of the Rocky Mountain Laboratories, NIAID, NIH, Hamilton, Montana and Dr Esther Schnettler of the University of Glasgow for provision of the TP21 strain of Langat virus, and Dr Christian Mandl and Prof. Franz Heinz of the Medical University of Vienna for provision of the Neudoerfl strain of tick-borne encephalitis virus and the pTND/∆ME plasmid. SW and MP were Early Stage Researchers supported by the POSTICK ITN (Post-graduate training network for capacity building to control ticks and tick-borne diseases) within the FP7- PEOPLE – ITN programme (EU Grant No. 238511). The research was supported by: grants BFU2011-23896 from the Ministerio de Economía y Competitividad, Spain and the European Union FP7 ANTIGONE project number 278976 (MV, JF); the Czech Science Foundation (GACR) (15-03044S) and by Institutional support RVO: 60077344 from Biology Centre ASCR, Institute of Parasitology (HT, LG); the Czech Science Foundation (projects nos P502/11/2116 and GA14-29256S) and the MEYS of the Czech Republic (project LO1218) under the NPU I programme (DR); the United Kingdom Biotechnology and Biological Sciences Research Council’s National Capability Grant to the Pirbright Institute (LBS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
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.
- Jongejan F, Uilenberg G. The global importance of ticks. Parasitology. 2004;129:S3–14.PubMedView ArticleGoogle Scholar
- Vaccines against tick-borne encephalitis. WHO position paper – recommendations. Vaccine. 2011;29:8769–70.View ArticleGoogle Scholar
- Süss J. Tick-borne encephalitis 2010: Epidemiology, risk areas, and virus strains in Europe and Asia — An overview. Ticks Tick Borne Dis. 2011;2:2–15.PubMedView ArticleGoogle Scholar
- Lindquist L, Vapalahti O. Tick-borne encephalitis. Lancet. 2008;371:1861–71.PubMedView ArticleGoogle Scholar
- Dobler G, Gniel D, Petermann R, Pfeffer M. Epidemiology and distribution of tick-borne encephalitis. Wien Med Wochenschr. 2012;162:230–8.PubMedView ArticleGoogle Scholar
- Smith CE. A virus resembling Russian spring–summer encephalitis virus from an Ixodid tick in Malaya. Nature. 1956;178:581–2.PubMedView ArticleGoogle Scholar
- Price WH, Thind IS, Teasdall RD, O’Leary W. Vaccination of human volunteers against Russian spring-summer (RSS) virus complex with attenuated Langat E5 virus. Bull World Hlth Org. 1970;42:89–94.Google Scholar
- Price WH, Thind IS. Immunization of mice against Russian spring-summer virus complex and monkeys against Powassan virus with attenuated Langat E5 virus. Am J Trop Med Hyg. 1973;22:100–8.PubMedGoogle Scholar
- Blair CD. Mosquito RNAi is the major innate immune pathway controlling arbovirus infection and transmission. Future Microbiol. 2011;6:265–77.PubMed CentralPubMedView ArticleGoogle Scholar
- Donald CL, Kohl A, Schnettler E. New insights into control of arbovirus replication and spread by insect RNA interference pathways. Insects. 2012;3:511–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Fragkoudis R, Attarzadeh-Yazdi G, Nash AA, Fazakerley JK, Kohl A. Advances in dissecting mosquito innate immune responses to arbovirus infection. J Gen Virol. 2009;90:2061–72.PubMedView ArticleGoogle Scholar
- Kingsolver MB, Huang Z, Hardy RW. Insect antiviral innate immunity: pathways, effectors, and connections. J Mol Biol. 2013;425:4921–36.PubMed CentralPubMedView ArticleGoogle Scholar
- Merkling SH, van Rij RP. Beyond RNAi: antiviral defense strategies in Drosophila and mosquito. J Insect Physiol. 2013;59:159–70.PubMedView ArticleGoogle Scholar
- Rückert C, Bell-Sakyi L, Fazakerley J, Fragkoudis R. Antiviral responses of arthropod vectors: an update on recent advances. VirusDisease. 2014;25:249–60.PubMed CentralPubMedView ArticleGoogle Scholar
- Kopácek P, Hajdusek O, Buresová V, Daffre S. Tick innate immunity. Adv Exp Med Biol. 2010;708:137–62.PubMedView ArticleGoogle Scholar
- Smith AA, Pal U. Immunity-related genes in Ixodes scapularis – perspectives from genome information. Front Cell Infect Microbiol. 2014;4:116. doi:https://doi.org/10.3389/fcimb.2014.00116.PubMed CentralPubMedView ArticleGoogle Scholar
- Taylor D. Innate immunity in ticks : a review. Acarol Soc Japan. 2006;15:109–27.View ArticleGoogle Scholar
- Hajdušek O, Síma R, Ayllón N, Jalovecká M, Perner J, de la Fuente J, et al. Interaction of the tick immune system with transmitted pathogens. Front Cell Infect Microbiol. 2013;3:26. doi:https://doi.org/10.3389/fcimb.2013.00026.PubMed CentralPubMedView ArticleGoogle Scholar
- Ayllón N, Villar M, Galindo RC, Kocan KM, Šíma R, López 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, e1005120. Available: http://dx.plos.org/https://doi.org/10.1371/journal.pgen.1005120.
- Liu L, Dai J, Zhao YO, Narasimhan S, Yang Y, Zhang L, et al. Ixodes scapularis JAK-STAT pathway regulates tick antimicrobial peptides, thereby controlling the agent of human granulocytic anaplasmosis. J Infect Dis. 2012;206:1233–41.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurscheid S, Lew-Tabor AE, Rodriguez Valle M, Bruyeres AG, Doogan VJ, Munderloh UG, et al. Evidence of a tick RNAi pathway by comparative genomics and reverse genetics screen of targets with known loss-of-function phenotypes in Drosophila. BMC Mol Biol. 2009;10:1–21.View ArticleGoogle Scholar
- Zivkovic Z, Torina A, Mitra R, Alongi A, Scimeca S, Kocan KM, et al. Subolesin expression in response to pathogen infection in ticks. BMC Immunol. 2010;11:7. doi:https://doi.org/10.1186/1471-2172-11-7.PubMed CentralPubMedView ArticleGoogle Scholar
- De la Fuente J, Kocan KM, Blouin EF, Zivkovic Z, Naranjo V, Almazán C, et al. Functional genomics and evolution of tick-Anaplasma interactions and vaccine development. Vet Parasitol. 2010;167:175–86.PubMedView ArticleGoogle Scholar
- Severo MS, Choy A, Stephens KD, Sakhon OS, Chen G, Chung D-WD, et al. The E3 ubiquitin ligase XIAP restricts Anaplasma phagocytophilum colonization of Ixodes scapularis ticks. J Infect Dis. 2013;208:1830–40.PubMed CentralPubMedView ArticleGoogle Scholar
- Garcia S, Billecocq A, Crance J-M, Munderloh UG, Garin D, Bouloy M. Nairovirus RNA sequences expressed by a Semliki Forest virus replicon induce RNA interference in tick cells. J Virol. 2005;79:8942–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Garcia S, Billecocq A, Crance J-M, Prins M, Garin D, Bouloy M. Viral suppressors of RNA interference impair RNA silencing induced by a Semliki Forest virus replicon in tick cells. J Gen Virol. 2006;87:1985–9.PubMedView ArticleGoogle Scholar
- Blouin EF, Manzano-Roman R, de la Fuente J, Kocan KM. Defining the role of subolesin in tick cell culture by use of RNA interference. Ann N Y Acad Sci. 2008;1149:41–4.PubMedView ArticleGoogle Scholar
- De la Fuente J, Kocan KM, Almazán C, Blouin EF. RNA interference for the study and genetic manipulation of ticks. Trends Parasitol. 2007;23:427–33.PubMedView ArticleGoogle Scholar
- Barry G, Alberdi P, Schnettler E, Weisheit S, Kohl A, Fazakerley JK, et al. Gene silencing in tick cell lines using small interfering or long double-stranded RNA. Exp Appl Acarol. 2013;59:319–38.PubMed CentralPubMedView ArticleGoogle Scholar
- Kopecký J, Stanková I. Interaction of virulent and attenuated tick-borne encephalitis virus strains in ticks and a tick cell line. Folia Parasitol (Praha). 1998;45:245–50.Google Scholar
- Lawrie CH, Uzcátegui NY, Armesto M, Bell-Sakyi L, Gould EA. Susceptibility of mosquito and tick cell lines to infection with various flaviviruses. Med Vet Entomol. 2004;18:268–74.PubMedView ArticleGoogle Scholar
- Senigl F, Kopecký J, Grubhoffer L. Distribution of E and NS1 proteins of TBE virus in mammalian and tick cells. Folia Microbiol (Praha). 2004;49:213–6.View ArticleGoogle Scholar
- Senigl F, Grubhoffer L, Kopecký J. Differences in maturation of tick-borne encephalitis virus in mammalian and tick cell line. Intervirology. 2006;49:239–48.PubMedView ArticleGoogle Scholar
- Růzek D, Bell-Sakyi L, Kopecký J, Grubhoffer L. Growth of tick-borne encephalitis virus (European subtype) in cell lines from vector and non-vector ticks. Virus Res. 2008;137:142–6.PubMedView ArticleGoogle Scholar
- Schrauf S, Mandl CW, Bell-Sakyi L, Skern T. Extension of flavivirus protein C differentially affects early RNA synthesis and growth in mammalian and arthropod host cells. J Virol. 2009;83:11201–10.PubMed CentralPubMedView ArticleGoogle Scholar
- Offerdahl DK, Dorward DW, Hansen BT, Bloom ME. A three-dimensional comparison of tick-borne flavivirus infection in mammalian and tick cell lines. PLoS One. 2012;7, e47912. doi:https://doi.org/10.1371/journal.pone.0047912.PubMed CentralPubMedView ArticleGoogle Scholar
- Yoshii K, Yanagihara N, Ishizuka M, Sakai M, Kariwa H. N-linked glycan in tick-borne encephalitis virus envelope protein affects viral secretion in mammalian cells, but not in tick cells. J Gen Virol. 2013;94:2249–58. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23824303.PubMedView ArticleGoogle Scholar
- Schnettler E, Tykalová H, Watson M, Sharma M, Sterken MG, Obbard DJ, et al. Induction and suppression of tick cell antiviral RNAi responses by tick-borne flaviviruses. Nucleic Acids Res. 2014;42:9436–46.PubMed CentralPubMedView ArticleGoogle Scholar
- Munderloh UG, Liu Y, Wang M, Chen C, Kurtti TJ. Establishment, maintenance and description of cell lines from the tick Ixodes scapularis. J Parasitol. 1994;80:533–43.PubMedView ArticleGoogle Scholar
- Munderloh UG, Kurtti TJ. Formulation of medium for tick cell culture. Exp Appl Acarol. 1989;7:219–29.PubMedView ArticleGoogle Scholar
- Bell-Sakyi L, Zweygarth E, Blouin EF, Gould EA, Jongejan F. Tick cell lines: tools for tick and tick-borne disease research. Trends Parasitol. 2007;230:450–7.View ArticleGoogle Scholar
- Bell-Sakyi L. Ehrlichia ruminantium grows in cell lines from four ixodid tick genera. J Comp Pathol. 2004;130:285–93.PubMedView ArticleGoogle Scholar
- Kozuch O, Mayer V. Pig kidney epithelial (PS) cells: a perfect tool for the study of flaviviruses and some other arboviruses. Acta Virol. 1975;19:498.PubMedGoogle Scholar
- De Madrid AT, Porterfield JS. A simple micro-culture method for the study of group B arboviruses. Bull World Hlth Org. 1969;40:113–21.Google Scholar
- Růzek D, Vancová M, Tesarová M, Ahantarig A, Kopecký J, Grubhoffer L. Morphological changes in human neural cells following tick-borne encephalitis virus infection. J Gen Virol. 2009;90:1649–58.PubMedView ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36. doi:https://doi.org/10.1186/gb-2013-14-4-r36.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.PubMed CentralPubMedView ArticleGoogle Scholar
- Sonenshine DE, Bissinger BW, Egekwu N, Donohue KV, Khalil SM, Roe RM. First transcriptome of the testis-vas deferens-male accessory gland and proteome of the spermatophore from Dermacentor variabilis (Acari: Ixodidae). PLoS One. 2011;6, e24711. doi:https://doi.org/10.1371/journal.pone.0024711.PubMed CentralPubMedView ArticleGoogle Scholar
- Francischetti IM, Anderson JM, Manoukis N, Pham VM, Ribeiro JM. An insight into the sialotranscriptome and proteome of the coarse bontlegged tick, Hyalomma marginatum rufipes. J Proteomics. 2011;74:2892–908.PubMed CentralPubMedView ArticleGoogle Scholar
- Schicht S, Qi W, Poveda L, Strube C. Whole transcriptome analysis of the poultry red mite Dermanyssus gallinae (De Geer, 1778). Parasitology. 2014;141:336–46.PubMedView ArticleGoogle Scholar
- Villar M, Popara M, Ayllón N, de Mera IG F, Mateos-Hernández L, Galindo RC, et al. A systems biology approach to the characterization of stress response in Dermacentor reticulatus tick unfed larvae. PLoS One. 2014;9, e89564. doi:https://doi.org/10.1371/journal.pone.0089564.
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106. doi:https://doi.org/10.1186/gb-2010-11-10-r106.PubMed CentralPubMedView ArticleGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq - A Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.PubMedView ArticleGoogle Scholar
- Shevchenko A, Tomas H, Havlis J, Olsen JV, Mann M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat Protoc. 2006;1:2856–60.PubMedView ArticleGoogle Scholar
- Popara M, Villar M, Mateos-Hernández L, de Mera IGF, Marina A, del Valle M, et al. Lesser protein degradation machinery correlates with higher BM86 tick vaccine efficacy in Rhipicephalus annulatus when compared to Rhipicephalus microplus. Vaccine. 2013;31:4728–35.PubMedView ArticleGoogle Scholar
- Whelan JA, Russel NB, Whelan MA. A method for the absolute quantification of cDNA using real-time PCR. J Immunol Methods. 2003;278:261–9.PubMedView ArticleGoogle Scholar
- Gehrke R, Ecker M, Aberle SW, Allison SL, Heinz FX, Mandl CW. Incorporation of tick-borne encephalitis virus replicons into virus-like particles by a packaging cell line. J Virol. 2003;77:8924–33.PubMed CentralPubMedView ArticleGoogle 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. 2001;25:402–8.PubMedView ArticleGoogle Scholar
- Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.PubMedView ArticleGoogle Scholar
- Bhat BKM, Yunker CE. Susceptibility of a tick cell line (Dermacentor parumapertus Neumann) to infection with arboviruses. In: Kurstak E, editor. Arctic and Tropical arboviruses. New York, NY: Academic Press; 1979. p. 263–75.Google Scholar
- Leake CJ, Pudney M, Varma M. Studies on arboviruses in established tick cell lines. In: Kurstak E, Maramorosch K, Dubendorfer A, editors. Invertebrate systems in vitro. Amsterdam: Elsevier/North Holland Biomedical Press; 1980. p. 327–35.Google Scholar
- Rehacek J. Arthropod cell cultures in studies of tick-borne togaviruses and orbiviruses in Central Europe. In: Yunker CE, editor. Arboviruses in arthropod cells in vitro, vol. I. Boca Raton: CRC Press; 1987. p. 115–32.Google Scholar
- Wallner G, Mandl CW, Kunz C, Heinz FX. The flavivirus 3′-noncoding region: extensive size heterogeneity independent of evolutionary relationships among strains of tick-borne encephalitis virus. Virology. 1995;213:169–78.PubMedView ArticleGoogle Scholar
- Bonizzoni M, Dunn WA, Campbell CL, Olson KE, Marinotti O, James AA. Complex modulation of the Aedes aegypti transcriptome in response to dengue virus infection. PLoS One. 2012;7, e50512. doi:https://doi.org/10.1371/journal.pone.0050512.PubMed CentralPubMedView ArticleGoogle Scholar
- McNally KL, Mitzel DN, Anderson JM, Ribeiro JMC, Valenzuela JG, Myers TG, et al. Differential salivary gland transcript expression profile in Ixodes scapularis nymphs upon feeding or flavivirus infection. Ticks Tick Borne Dis. 2012;3:18–26.PubMed CentralPubMedView ArticleGoogle Scholar
- De la Fuente J, Blouin EF, Manzano-Roman R, Naranjo V, Almazán C, Pérez de la Lastra JM, et al. Functional genomic studies of tick cells in response to infection with the cattle pathogen, Anaplasma marginale. Genomics. 2007;90:712–22.PubMedView ArticleGoogle Scholar
- Gibson AK, Smith Z, Fuqua C, Clay K, Colbourne JK. Why so many unknown genes? Partitioning orphans from a representative transcriptome of the lone star tick Amblyomma americanum. BMC Genomics. 2013;14:135. doi:https://doi.org/10.1186/1471-2164-14-135.
- Sim S, Dimopoulos G. Dengue virus inhibits immune responses in Aedes aegypti cells. PLoS One. 2010;5, e10678. doi:https://doi.org/10.1371/journal.pone.0010678.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang M, Zheng X, Wu Y, Gan M, He A, Li Z, et al. Differential proteomics of Aedes albopictus salivary gland, midgut and C6/36 cell induced by dengue virus infection. Virology. 2013;444:109–18.PubMedView ArticleGoogle Scholar
- Villar M, Ayllón N, Busby AT, Galindo RC, Blouin EF, Kocan KM, et al. 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 Proteomics. 2010;2010:657261. doi:https://doi.org/10.1155/2010/657261.
- Villar M, Torina A, Nuñez Y, Zivkovic Z, Marina A, Alongi A, et al. Application of highly sensitive saturation labeling to the analysis of differential protein expression in infected ticks from limited samples. Proteome Sci. 2010;8:43. doi:https://doi.org/10.1186/1477-5956-8-43.PubMed CentralPubMedView ArticleGoogle Scholar
- Evans V, Barker G, Heesom K, Fan J, Bessant C, Matthews D. De novo derivation of proteomes from transcriptomes for transcript and protein identification. Nat Methods. 2012;9:1207–11.PubMed CentralPubMedView ArticleGoogle Scholar
- Mudenda L, Pierlé SA, Turse JE, Scoles GA, Purvine SO, Nicora CD, et al. Proteomics informed by transcriptomics identifies novel secreted proteins in Dermacentor andersoni saliva. Int J Parasitol. 2014;44:1029–37.PubMedView ArticleGoogle Scholar
- Nookaew I, Papini M, Pornputtapong N, Scalcinati G, Fagerberg L, Uhlén M, et al. A comprehensive comparison of RNA-Seq-based transcriptome analysis from reads to differential gene expression and cross-comparison with microarrays: a case study in Saccharomyces cerevisiae. Nucleic Acids Res. 2012;40:10084–97.PubMed CentralPubMedView ArticleGoogle Scholar
- Devonshire AS, Sanders R, Wilkes TM, Taylor MS, Foy CA, Huggett JF. Application of next generation qPCR and sequencing platforms to mRNA biomarker analysis. Methods. 2013;59:89–100.PubMedView ArticleGoogle Scholar
- Hegedus Z, Zakrzewska A, Agoston VC, Ordas A, Rácz P, Mink M, et al. Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol. 2009;46:2918–30.PubMedView ArticleGoogle Scholar
- Zeng D, Chen X, Xie D, Zhao Y, Yang C, Li Y, et al. Transcriptome analysis of Pacific white shrimp (Litopenaeus vannamei) hepatopancreas in response to Taura syndrome Virus (TSV) experimental infection. PLoS One. 2013;8, e57515. doi:https://doi.org/10.1371/journal.pone.0057515.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhu J-Y, Yang P, Zhang Z, Wu G-X, Yang B. Transcriptomic immune response of Tenebrio molitor pupae to parasitization by Scleroderma guani. PLoS One. 2013;8, e54411. doi:https://doi.org/10.1371/journal.pone.0054411.PubMed CentralPubMedView ArticleGoogle Scholar
- Price WH, Parks JJ, Ganaway J, O’Leary W, Lee R. The ability of an attenuated isolate of Langat virus to protect primates and mice against other members of the Russian spring-summer virus complex. Am J Trop Med Hyg. 1963;12:787–99.PubMedGoogle Scholar
- Bellés X. Beyond Drosophila: RNAi in vivo and functional genomics in insects. Annu Rev Entomol. 2010;55:111–28.PubMedView ArticleGoogle Scholar
- Terenius O, Papanicolaou A, Garbutt JS, Eleftherianos I, Huvenne H, Kanginakudru S, et al. RNA interference in Lepidoptera: an overview of successful and unsuccessful studies and implications for experimental design. J Insect Physiol. 2011;57:231–45.PubMedView ArticleGoogle Scholar
- Sim S, Ramirez JL, Dimopoulos G. Dengue virus infection of the Aedes aegypti salivary gland and chemosensory apparatus induces genes that modulate infection and blood-feeding behavior. PLoS Pathog. 2012;8, e1002631. doi:https://doi.org/10.1371/journal.ppat.1002631.PubMed CentralPubMedView ArticleGoogle Scholar
- Makou E, Herbert AP, Barlow PN. Functional anatomy of complement factor H. Biochemistry. 2013;52:3949–62.PubMedView ArticleGoogle Scholar
- Blue CE, Spiller OB, Blackbourn DJ. The relevance of complement to virus biology. Virology. 2004;319:176–84.PubMedView ArticleGoogle Scholar
- Favoreel HW. Virus complement evasion strategies. J Gen Virol. 2003;84:1–15.PubMedView ArticleGoogle Scholar
- Lachmann PJ, Davies A. Complement and immunity to viruses. Immunol Rev. 1997;159:69–77.PubMedView ArticleGoogle Scholar
- Chung KM, Liszewski MK, Nybakken G, Davis AE, Townsend RR, Fremont DH, et al. West Nile virus nonstructural protein NS1 inhibits complement activation by binding the regulatory protein factor H. Proc Natl Acad Sci U S A. 2006;103:19111–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Buresova V, Hajdusek O, Franta Z, Loosova G, Grunclova L, Levashina EA, et al. Functional genomics of tick thioester-containing proteins reveal the ancient origin of the complement system. J Innate Immun. 2011;3:623–30.PubMedView ArticleGoogle Scholar
- Kopacek P, Hajdusek O, Buresova V. Tick as a model for the study of a primitive complement system. In: Mylonakis E, Ausubel FM, Gilmore M, Casadevall A, editors. Recent advances on model hosts. New York, NY: Springer New York; 2012. p. 83–93.View ArticleGoogle Scholar
- Buresova V, Hajdusek O, Franta Z, Sojka D, Kopacek P. IrAM—An alpha2-macroglobulin from the hard tick Ixodes ricinus: characterization and function in phagocytosis of a potential pathogen Chryseobacterium indologenes. Dev Comp Immunol. 2009;33:489–98.PubMedView ArticleGoogle Scholar
- Hill CA, Wikel SK. The Ixodes scapularis Genome Project: an opportunity for advancing tick research. Trends Parasitol. 2005;21:151–3.PubMedView ArticleGoogle Scholar
- Iwasaki S, Kobayashi M, Yoda M, Sakaguchi Y, Katsuma S, Suzuki T, et al. Hsc70/Hsp90 chaperone machinery mediates ATP-dependent RISC loading of small RNA duplexes. Mol Cell. 2010;39:292–9.PubMedView ArticleGoogle Scholar
- De Gregorio E, Spellman PT, Rubin GM, Lemaitre B. Genome-wide analysis of the Drosophila immune response by using oligonucleotide microarrays. Proc Natl Acad Sci U S A. 2001;98:12590–5.PubMed CentralPubMedView ArticleGoogle Scholar
- De Gregorio E, Han S-J, Lee W-J, Baek M-J, Osaki T, Kawabata S-I, et al. An immune-responsive Serpin regulates the melanization cascade in Drosophila. Dev Cell. 2002;3:581–92.PubMedView ArticleGoogle Scholar
- Ligoxygakis P, Pelte N, Hoffmann JA, Reichhart J-M. Activation of Drosophila Toll during fungal infection by a blood serine protease. Science. 2002;297:114–6.PubMedView ArticleGoogle Scholar
- Hung J-J, Chung C-S, Chang W. Molecular chaperone Hsp90 is important for vaccinia virus growth in cells. J Virol. 2002;76:1379–90.PubMed CentralPubMedView ArticleGoogle Scholar
- Connor JH, McKenzie MO, Parks GD, Lyles DS. Antiviral activity and RNA polymerase degradation following Hsp90 inhibition in a range of negative strand viruses. Virology. 2007;362:109–19.PubMed CentralPubMedView ArticleGoogle Scholar
- Morales C, Wu S, Yang Y, Hao B, Li Z. Drosophila glycoprotein 93 is an ortholog of mammalian heat shock protein gp96 (grp94, HSP90b1, HSPC4) and retains disulfide bond-independent chaperone function for TLRs and integrins. J Immunol. 2009;183:5121–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Attoui H, Stirling JM, Munderloh UG, Billoir F, Brookes SM, Burroughs JN, et al. Complete sequence characterization of the genome of the St Croix River virus, a new orbivirus isolated from cells of Ixodes scapularis. J Gen Virol. 2001;82:795–804.PubMedView ArticleGoogle Scholar
- Alberdi MP, Dalby MJ, Rodriguez-Andres J, Fazakerley JK, Kohl A, Bell-Sakyi L. Detection and identification of putative bacterial endosymbionts and endogenous viruses in tick cell lines. Ticks Tick Borne Dis. 2012;3:137–46.PubMed CentralPubMedView ArticleGoogle Scholar
- Keene KM, Foy BD, Sanchez-Vargas I, Beaty BJ, Blair CD, Olson KE. RNA interference acts as a natural antiviral response to O’nyong-nyong virus (Alphavirus; Togaviridae) infection of Anopheles gambiae. Proc Natl Acad Sci U S A. 2004;101:17240–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Saleh M-C, Tassetto M, van Rij RP, Goic B, Gausson V, Berry B, et al. Antiviral immunity in Drosophila requires systemic RNA interference spread. Nature. 2009;458:346–50.PubMed CentralPubMedView ArticleGoogle Scholar
- Pitaluga AN, Mason PW, Traub-Cseko YM. Non-specific antiviral response detected in RNA-treated cultured cells of the sandfly, Lutzomyia longipalpis. Dev Comp Immunol. 2008;32:191–7.PubMedView ArticleGoogle Scholar
- Robalino J, Bartlett TC, Chapman RW, Gross PS, Browdy CL, Warr GW. Double-stranded RNA and antiviral immunity in marine shrimp: inducible host mechanisms and evidence for the evolution of viral counter-responses. Dev Comp Immunol. 2007;31:539–47.PubMedView ArticleGoogle Scholar
- Flenniken ML, Andino R. Non-specific dsRNA-mediated antiviral response in the honey bee. PLoS One. 2013;8, e77263. doi:https://doi.org/10.1371/journal.pone.0077263.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu J, Smagghe G, Swevers L. Transcriptional response of BmToll9-1 and RNAi machinery genes to exogenous dsRNA in the midgut of Bombyx mori. J Insect Physiol. 2013;59:646–54.PubMedView ArticleGoogle Scholar
- Deddouche S, Matt N, Budd A, Mueller S, Kemp C, Galiana-Arnoux D, et al. The DExD/H-box helicase Dicer-2 mediates the induction of antiviral activity in Drosophila. Nat Immunol. 2008;9:1425–32.PubMedView ArticleGoogle Scholar
- Paradkar PN, Trinidad L, Voysey R, Duchemin J-B, Walker PJ. Secreted Vago restricts West Nile virus infection in Culex mosquito cells by activating the Jak-STAT pathway. Proc Natl Acad Sci U S A. 2012;109:18915–20.PubMed CentralPubMedView ArticleGoogle Scholar
- Rückert C. Alphavirus and flavivirus infection of Ixodes tick cell lines: an insight into tick antiviral immunity. PhD thesis, University of Edinburgh. 2014. Available: http://hdl.handle.net/1842/10063
- Aung KM, Boldbaatar D, Umemiya-Shirafuji R, Liao M, Tsuji N, Xuenan X, et al. HlSRB, a Class B scavenger receptor, is key to the granulocyte-mediated microbial phagocytosis in ticks. PLoS One. 2012;7, e33504. doi:https://doi.org/10.1371/journal.pone.0033504.PubMed CentralPubMedView ArticleGoogle Scholar
- Aung KM, Boldbaatar D, Umemiya-Shirafuji R, Liao M, Xuenan X, Suzuki H, et al. Scavenger receptor mediates systemic RNA interference in ticks. PLoS One. 2011;6, e28407. doi:https://doi.org/10.1371/journal.pone.0028407.PubMed CentralPubMedView ArticleGoogle 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, e1000062. doi:https://doi.org/10.1371/journal.ppat.1000062.
- Johansson MW, Söderhäll K. Isolation and purification of a cell adhesion factor from crayfish blood cells. J Cell Biol. 1988;106:1795–803.PubMedView ArticleGoogle Scholar
- Kobayashi M, Johansson M, Söderhäll K. The 76 kDa cell-adhesion factor from crayfish haemocytes promotes encapsulation in vitro. Cell Tissue Res. 1990;260:113–8.View ArticleGoogle Scholar
- Thörnqvist P, Johansson MW, Söderhäll K. Opsonic activity of cell adhesion proteins and β-1,3-glucan binding proteins from two crustaceans. Dev Comp Immunol. 1994;18:3–12.PubMedView ArticleGoogle Scholar
- Johansson MW, Lind MI, Holmblad T, Thörnqvist P-O, Söderhäll K. Peroxinectin, a novel cell adhesion protein from crayfish blood. Biochem Biophys Res Commun. 1995;216:1079–89.PubMedView ArticleGoogle Scholar
- Du Z-Q, Ren Q, Huang A-M, Fang W-H, Zhou J-F, Gao L-J, et al. A novel peroxinectin involved in antiviral and antibacterial immunity of mud crab, Scylla paramamosain. Mol Biol Rep. 2013;40:6873–81.Google Scholar
- Lin X, Cerenius L, Lee BL, Söderhäll K. Purification of properoxinectin, a myeloperoxidase homologue and its activation to a cell adhesion molecule. Biochim Biophys Acta. 2007;1770:87–93.PubMedView ArticleGoogle Scholar
- Sritunyalucksana K, Wongsuebsantati K, Johansson MW, Söderhäll K. Peroxinectin, a cell adhesive protein associated with the proPO system from the black tiger shrimp, Penaeus monodon. Dev Comp Immunol. 2001;25:353–63.Google Scholar
- Zhioua E, Yeh MT, LeBrun A. Assay for phenoloxidase activity in Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis. J Parasitol. 1997;83:553–4.PubMedView ArticleGoogle Scholar
- Kadota K, Satoh E, Ochiai M, Inoue N, Tsuji N, Igarashi I, et al. Existence of phenol oxidase in the argasid tick Ornithodoros moubata. Parasitol Res. 2002;88:781–4.PubMedView ArticleGoogle Scholar
- Simser JA, Mulenga A, Macaluso KR, Azad AF. An immune responsive factor D-like serine proteinase homologue identified from the American dog tick, Dermacentor variabilis. Insect Mol Biol. 2004;13:25–35.PubMedView ArticleGoogle Scholar
- Caudy AA, Ketting RF, Hammond SM, Denli AM, Bathoorn AMP, Tops BBJ, et al. A micrococcal nuclease homologue in RNAi effector complexes. Nature. 2003;425:411–4.PubMedView ArticleGoogle Scholar
- Scadden ADJ. The RISC subunit Tudor-SN binds to hyper-edited double-stranded RNA and promotes its cleavage. Nat Struct Mol Biol. 2005;12:489–96.PubMedView ArticleGoogle Scholar
- Yang W, Chendrimada TP, Wang Q, Higuchi M, Seeburg PH, Shiekhattar R, et al. Modulation of microRNA processing and expression through RNA editing by ADAR deaminases. Nat Struct Mol Biol. 2006;13:13–21.PubMed CentralPubMedView ArticleGoogle Scholar
- Ayllón N, Naranjo V, Hajdušek O, Villar M, Galindo RC, Kocan KM, et al. Nuclease Tudor-SN is involved in tick dsRNA-mediated RNA interference and feeding but not in defense against flaviviral or Anaplasma phagocytophilum rickettsial infection. PLoS One. 2015;10, e0133038. doi:https://doi.org/10.1371/journal.pone.0133038.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu Q, Luo Y, Lu R, Lau N, Lai EC, Li W-X, et al. Virus discovery by deep sequencing and assembly of virus-derived small silencing RNAs. Proc Natl Acad Sci U S A. 2010;107:1606–11.PubMed CentralPubMedView ArticleGoogle Scholar
- Vodovar N, Bronkhorst AW, van Cleef KWR, Miesen P, Blanc H, van Rij RP, et al. Arbovirus-derived piRNAs exhibit a ping-pong signature in mosquito cells. PLoS One. 2012;7, e30861. doi:https://doi.org/10.1371/journal.pone.0030861.PubMed CentralPubMedView ArticleGoogle Scholar
- Morazzani EM, Wiley MR, Murreddu MG, Adelman ZN, Myles KM. Production of virus-derived ping-pong-dependent piRNA-like small RNAs in the mosquito soma. PLoS Pathog. 2012;8, e1002470. doi:https://doi.org/10.1371/journal.ppat.1002470.PubMed CentralPubMedView ArticleGoogle Scholar
- Léger P, Lara E, Jagla B, Sismeiro O, Mansuroglu Z, Coppée JY, et al. Dicer-2- and Piwi-mediated RNA interference in Rift Valley fever virus-infected mosquito cells. J Virol. 2013;87:1631–48.PubMed CentralPubMedView ArticleGoogle Scholar
- Schnettler E, Donald CL, Human S, Watson M, Siu RWC, McFarlane M, et al. Knockdown of piRNA pathway proteins results in enhanced Semliki Forest virus production in mosquito cells. J Gen Virol. 2013;94:1680–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Srivastava P. Roles of heat-shock proteins in innate and adaptive immunity. Nat Rev Immunol. 2002;2:185–94.PubMedView ArticleGoogle Scholar
- Busby AT, Ayllón N, Kocan KM, Blouin EF, de la Fuente G, Galindo RC, et al. Expression of heat shock proteins and subolesin affects stress responses, Anaplasma phagocytophilum infection and questing behaviour in the tick, Ixodes scapularis. Med Vet Entomol. 2012;26:92–102.PubMedView ArticleGoogle Scholar
- Nagy PD, Wang RY, Pogany J, Hafren A, Makinen K. Emerging picture of host chaperone and cyclophilin roles in RNA virus replication. Virology. 2011;411:374–82.PubMedView ArticleGoogle Scholar
- Zhao L, Jones W. Expression of heat shock protein genes in insect stress responses. Invertebr Surviv J. 2012;9:93–101.Google Scholar
- Michalak M, Corbett EF, Mesaeli N, Nakamura K, Opas M. Calreticulin: one protein, one gene, many functions. Biochem J. 1999;344:281–92.PubMed CentralPubMedView ArticleGoogle Scholar
- Pieren M, Galli C, Denzel A, Molinari M. The use of calnexin and calreticulin by cellular and viral glycoproteins. J Biol Chem. 2005;280:28265–71.PubMedView ArticleGoogle Scholar
- Limjindaporn T, Wongwiwat W, Noisakran S, Srisawat C, Netsawang J, Puttikhunt C, et al. Interaction of dengue virus envelope protein with endoplasmic reticulum-resident chaperones facilitates dengue virus production. Biochem Biophys Res Commun. 2009;379:196–200.PubMedView ArticleGoogle Scholar
- Antunes S, Galindo RC, Almazán C, Rudenko N, Golovchenko M, Grubhoffer L, et al. Functional genomics studies of Rhipicephalus (Boophilus) annulatus ticks in response to infection with the cattle protozoan parasite, Babesia bigemina. Int J Parasitol. 2012;42:187–95.PubMedView ArticleGoogle Scholar
- Pastorino B, Nougairède A, Wurtz N, Gould E, de Lamballerie X. Role of host cell factors in flavivirus infection: Implications for pathogenesis and development of antiviral drugs. Antiviral Res. 2010;87:281–94.PubMedView ArticleGoogle Scholar
- Colpitts TM, Cox J, Vanlandingham DL, Feitosa FM, Cheng G, Kurscheid S, et al. Alterations in the Aedes aegypti transcriptome during infection with West Nile, dengue and yellow fever viruses. PLoS Pathog. 2011;7, e1002189. doi:https://doi.org/10.1371/journal.ppat.1002189.PubMed CentralPubMedView ArticleGoogle Scholar
- Wei H, Zhou M-M. Viral-encoded enzymes that target host chromatin functions. Biochim Biophys Acta. 2010;1799:296–301.PubMed CentralPubMedView ArticleGoogle Scholar
- Colpitts TM, Barthel S, Wang P, Fikrig E. Dengue virus capsid protein binds core histones and inhibits nucleosome formation in human liver cells. PLoS One. 2011;6, e24365. doi:https://doi.org/10.1371/journal.pone.0024365.PubMed CentralPubMedView ArticleGoogle Scholar
- Blackwell J, Brinton M. Translation elongation factor-1 alpha interacts with the 3’stem-loop region of West Nile virus genomic RNA. J Virol. 1997;71:6433–44.PubMed CentralPubMedGoogle Scholar
- Davis WG, Blackwell JL, Shi P-Y, Brinton MA. Interaction between the cellular protein eEF1A and the 3’-terminal stem-loop of West Nile virus genomic RNA facilitates viral minus-strand RNA synthesis. J Virol. 2007;81:10172–87.PubMed CentralPubMedView ArticleGoogle Scholar
- Pattanakitsakul S, Rungrojcharoenkit K, Kanlaya R, Sinchaikul S, Noisakran S, Chen S-T, et al. Proteomic analysis of host responses in HepG2 cells during dengue virus infection. J Proteome Res. 2007;6:4592–600.PubMedView ArticleGoogle Scholar
- Pastorino B, Boucomont-Chapeaublanc E, Peyrefitte CN, Belghazi M, Fusaï T, Rogier C, et al. Identification of cellular proteome modifications in response to West Nile virus infection. Mol Cell Proteomics. 2009;8:1623–37.PubMed CentralPubMedView ArticleGoogle Scholar
- Emara MM, Brinton MA. Interaction of TIA-1/TIAR with West Nile and dengue virus products in infected cells interferes with stress granule formation and processing body assembly. Proc Natl Acad Sci U S A. 2007;104:9041–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Morais AT, Terzian AC, Duarte DV, Bronzoni RV, Madrid MC, Gavioli AF, et al. The eukaryotic translation initiation factor 3 subunit L protein interacts with Flavivirus NS5 and may modulate yellow fever virus replication. Virol J. 2013;10:205. doi:https://doi.org/10.1186/1743-422X-10-205.PubMed CentralPubMedView ArticleGoogle Scholar
- Bartholomay L, Cho W, Rocheleau T, Boyle JP, Beck ET, Fuchs JF, et al. Description of the transcriptomes of immune response-activated hemocytes from the mosquito vectors Aedes aegypti and Armigeres subalbatus. Infect Immun. 2004;72:4114–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Waldock J, Olson KE, Christophides GK. Anopheles gambiae antiviral immune response to systemic O’nyong-nyong infection. PLoS Negl Trop Dis. 2012;6, e1565. doi:https://doi.org/10.1371/journal.pntd.0001565.PubMed CentralPubMedView ArticleGoogle Scholar
- Xi Z, Ramirez JL, Dimopoulos G. The Aedes aegypti toll pathway controls dengue virus infection. PLoS Pathog. 2008;4, e1000098. doi:https://doi.org/10.1371/journal.ppat.1000098.PubMed CentralPubMedView ArticleGoogle Scholar
- Lai R, Takeuchi H, Lomas LO, Jonczy J, Rigden DJ, Rees HH, et al. A new type of antimicrobial protein with multiple histidines from the hard tick, Amblyomma hebraeum. FASEB J. 2004;18:1447–9.PubMedGoogle Scholar
- Silva FD, Rezende CA, Rossi DCP, Esteves E, Dyszy FH, Schreier S, et al. Structure and mode of action of microplusin, a copper II-chelating antimicrobial peptide from the cattle tick Rhipicephalus (Boophilus) microplus. J Biol Chem. 2009;284:34735–46.PubMed CentralPubMedView ArticleGoogle Scholar
- Mulenga A, Macaluso KR, Simser JA, Azad AF. Dynamics of Rickettsia-tick interactions: identification and characterization of differentially expressed mRNAs in uninfected and infected Dermacentor variabilis. Insect Mol Biol. 2003;12:185–93.PubMedView ArticleGoogle Scholar
- Mulenga A, Simser JA, Macaluso KR, Azad A. Stress and transcriptional regulation of tick ferritin HC. Insect Mol Biol. 2004;13:423–33.PubMedView ArticleGoogle Scholar
- Ong ST, Ho JZS, Ho B, Ding JL. Iron-withholding strategy in innate immunity. Immunobiology. 2006;211:295–314.PubMedView ArticleGoogle Scholar
- Fric J, Zelante T, Wong AYW, Mertes A, Yu H-B, Ricciardi-Castagnoli P. NFAT control of innate immunity. Blood. 2012;120:1380–9.PubMedView ArticleGoogle Scholar
- Rao A, Luo C, Hogan PG. Transcription factors of the NFAT family: regulation and function. Annu Rev Immunol. 1997;15:707–47.PubMedView ArticleGoogle Scholar
- Wu H, Peisley A, Graef IA, Crabtree GR. NFAT signaling and the invention of vertebrates. Trends Cell Biol. 2007;17:251–60.PubMedView ArticleGoogle Scholar
- Zanoni I, Granucci F. Regulation and dysregulation of innate immunity by NFAT signaling downstream of pattern recognition receptors (PRRs). Eur J Immunol. 2012;42:1924–31.PubMedView ArticleGoogle Scholar
- Keyser P, Borge-Renberg K, Hultmark D. The Drosophila NFAT homolog is involved in salt stress tolerance. Insect Biochem Mol Biol. 2007;37:356–62.PubMedView ArticleGoogle Scholar
- Song X, Hu J, Jin P, Chen L, Ma F. Identification and evolution of an NFAT gene involving Branchiostoma belcheri innate immunity. Genomics. 2013;102:355–62.PubMedView ArticleGoogle Scholar
- Stafforini DM, McIntyre TM, Zimmerman GA, Prescott SM. Platelet-activating factor acetylhydrolases. J Biol Chem. 1997;272(29):17895–8.PubMedView ArticleGoogle Scholar
- McIntyre TM, Prescott SM, Stafforini DM. The emerging roles of PAF acetylhydrolase. J Lipid Res. 2009;50(Suppl):S255–9. doi:https://doi.org/10.1194/jlr.R800024-JLR200.PubMed CentralPubMedGoogle Scholar
- Prescott SM, Zimmerman GA, Stafforini DM, McIntyre TM. Platelet-activating factor and related lipid mediators. Annu Rev Bio. 2000;69:419–45.View ArticleGoogle Scholar
- Cheeseman MT, Bates PA, Crampton JM. Preliminary characterisation of esterase and platelet-activating factor (PAF)-acetylhydrolase activities from cat flea (Ctenocephalides felis) salivary glands. Insect Biochem Mol Biol. 2001;31:157–64.PubMedView ArticleGoogle Scholar
- Figueiredo MB, Genta FA, Garcia ES, Azambuja P. Lipid mediators and vector infection: Trypanosoma rangeli inhibits Rhodnius prolixus hemocyte phagocytosis by modulation of phospholipase A2 and PAF-acetylhydrolase activities. J Insect Physiol. 2008;54:1528–37.PubMedView ArticleGoogle Scholar
- Garcia ES, Castro DP, Figueiredo MB, Genta FA, Azambuja P. Trypanosoma rangeli: a new perspective for studying the modulation of immune reactions of Rhodnius prolixus. Parasit Vectors. 2009;2:33. doi:https://doi.org/10.1186/1756-3305-2-33.PubMed CentralPubMedView ArticleGoogle Scholar
- Rachinsky A, Guerrero FD, Scoles GA. Differential protein expression in ovaries of uninfected and Babesia-infected southern cattle ticks, Rhipicephalus (Boophilus) microplus. Insect Biochem Mol Biol. 2007;37:1291–308.PubMedView ArticleGoogle Scholar
- Behnke J, Hendershot LM. The large Hsp70 Grp170 binds to unfolded protein substrates in vivo with a regulation distinct from conventional Hsp70s. J Biol Chem. 2013;289:2899–907.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu S-W, Wang H, Poitras MF, Coombs C, Bowers WJ, Federoff HJ, et al. Mediation of poly(ADP-ribose) polymerase-1-dependent cell death by apoptosis-inducing factor. Science. 2002;297:259–63.PubMedView ArticleGoogle Scholar
- Odorizzi G. The multiple personalities of Alix. J Cell Sci. 2006;119:3025–32.PubMedView ArticleGoogle Scholar
- Le Blanc I, Luyet P-P, Pons V, Ferguson C, Emans N, Petiot A, et al. Endosome-to-cytosol transport of viral nucleocapsids. Nat Cell Biol. 2005;7:653–64.PubMed CentralPubMedView ArticleGoogle Scholar
- Sangsuriya P, Rojtinnakorn J, Senapin S, Flegel TW. Identification and characterization of Alix/AIP1 interacting proteins from the black tiger shrimp, Penaeus monodon. J Fish Dis. 2010;33:571–81.PubMedView ArticleGoogle Scholar
- Anderson P, Kedersha N. RNA granules. J Cell Biol. 2006;172:803–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Anderson P, Kedersha N. Stress granules: the Tao of RNA triage. Trends Biochem Sci. 2008;33:141–50.PubMedView ArticleGoogle Scholar
- Valiente-Echeverría F, Melnychuk L, Mouland AJ. Viral modulation of stress granules. Virus Res. 2012;169:430–7.PubMedView ArticleGoogle Scholar
- Lloyd RE. Regulation of stress granules and P-bodies during RNA virus infection. WIREs RNA. 2013;4:317–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Li W, Li Y, Kedersha N, Anderson P, Emara M, Swiderek KMK, et al. Cell proteins TIA-1 and TIAR interact with the 3’ stem-loop of the West Nile virus complementary minus-strand RNA and facilitate virus replication. J Virol. 2002;76:11989–2000.PubMed CentralPubMedView ArticleGoogle Scholar