<title>Background</title> <p>Ticks are obligate hematophagous ectoparasites of wild and domestic animals and humans, and are important vectors of diseases to humans and animals worldwide [##REF##11298147##1##]. Ticks are classified in the subclass Acari, order Parasitiformes, suborder Ixodida and are distributed worldwide from Arctic to tropical regions [##REF##15938503##2##]. Despite efforts to control tick infestations, these ectoparasites remain a serious threat for human and animal health [##REF##16842264##3##,##REF##16497440##4##]. Recently, both vaccine studies using key tick antigens as well as characterization of tick gene function by RNA interference (RNAi) have provided new information on genes that impact tick life cycle and the tick-pathogen interface [##REF##16842264##3##, ####REF##16497440##4##, ##REF##17692140##5##, ##REF##17964755##6####17964755##6##].</p> <p>One of these genes, subolesin (also called 4D8), was discovered recently in <italic>Ixodes scapularis </italic>and was shown by both RNAi and immunization with recombinant proteins to protect against tick infestations, resulting in reduced tick survival, feeding and reproduction [##REF##12615446##7##, ####REF##16005748##8##, ##REF##16153760##9##, ##REF##16580098##10##, ##REF##16630571##11####16630571##11##]. The silencing of subolesin expression by RNAi in ticks resulted in degeneration of gut, salivary gland, reproductive and embryonic tissues as well as causing sterility in males [##REF##16580098##10##, ####REF##16630571##11##, ##REF##17196597##12##, ##REF##17372765##13####17372765##13##]. In addition, targeting subolesin by RNAi or vaccination decreased the ability of ticks to become infected with <italic>Anaplasma marginale </italic>or <italic>A. phagocytophilum </italic>[##REF##16816958##14##,##REF##17656154##15##].</p> <p>Evidence of the conservation of subolesin throughout evolution was provided by the high homology of amino acid sequences in higher eukaryotes, which suggests an essential conserved biological function for this protein [##REF##16580098##10##]. For example, the expression of subolesin orthologs has been detected in a variety of adult and immature tissues of several tick species [##REF##16005748##8##,##REF##16580098##10##,##REF##17196597##12##], <italic>Drosophila melanogaster </italic>[##REF##11901120##16##,##UREF##0##17##] and <italic>Caenorhabditis elegans </italic>[##REF##11231151##18##]. These studies have suggested that subolesin orthologs may be involved in the control of developmental processes in these organisms [##REF##16005748##8##,##REF##16580098##10##,##REF##17196597##12##,##REF##11901120##16##, ####UREF##0##17##, ##REF##11231151##18####11231151##18##]. However, despite the important role that subolesin plays in tick reproduction, development and pathogen infection, the biological function of subolesin and its orthologs has not been reported.</p> <p>Because of the profound effect of subolesin knockdown in ticks and other organisms we hypothesized that subolesin may play a role in gene expression, thus affecting multiple cellular processes. Herein, gene expression is understood as the process by which a gene gets turned on in a cell to make RNA and proteins and therefore may be affected at the transcriptional and/or translational levels. The objective of this study was to provide evidence for the role of subolesin in gene expression through a combination of methodological approaches that included characterization of subolesin-interacting proteins, the effect of gene knockdown on tick gene expression pattern and the prediction of subolesin post-translational modifications. Although the biological function of subolesin is unkown, the results provided evidence that this protein plays a role in gene expression in ticks and most likely other organisms.</p>
<title>Methods</title> <title>Ticks</title> <p><italic>I. scapularis </italic>female ticks were obtained from the laboratory colony maintained at the Oklahoma State University tick rearing facility. Off-host ticks were maintained in a 12 hr light: 12 hr dark photoperiod at 22–25°C and 95% relative humidity. <italic>R. microplus </italic>female ticks (Mozambique strain) were reared in Holstein-Friesian cattle at the Utrecht Centre for Tick-borne Diseases, Utrecht University. Animals were housed with the approval and supervision of the respective Institutional Animal Care and Use Committees.</p> <title>cDNA library construction</title> <p><italic>R. microplus </italic>eggs were used for RNA extraction and cDNA library construction. cDNA synthesis and amplification was performed using the Super SMART System™ principle (Clontech Laboratories, Mountain View, CA, USA) with adapted anchor primers containing unique <italic>Sfi</italic>I restriction sites for directional cloning (SMART IV: 5'-AAGCAGTGGTATCAACGCAGAGTGGCCATGGAGGCCGGG-3'and CDS III: 5'-ATTCTAGAGGCCTCCATGGCCGACATG(T)30VN-3'). Amplified cDNA was polished using the protocol described by the SMART PCR cDNA synthesis manual (Clontech Laboratories), purified using the DNA Extract II Kit (Macherey-Nagel, Düren, Germany) and subjected to directional cloning in pACT2 plasmid via the <italic>Sfi</italic>I sites. Repetitive electroporation of <italic>E. coli </italic>JM109 cells was performed to obtain a library with a titer exceeding 3 × 10<sup>7 </sup>cfu/ml.</p> <title>Yeast two-hybrid screen</title> <p>The yeast two-hybrid screen was performed as dictated in the MATCHMAKER Two-Hybrid user manual (Clontech Laboratories). Full-length subolesin cDNA from <italic>R. microplus </italic>was inserted in yeast expression vector pAS2_1 at <italic>Nde</italic>I and <italic>Pst</italic>I sites. Yeast strain AH109 was co-transformed with plasmid pAS2_1_subolesin (bait) and pACT2 <italic>R. microplus </italic>egg cDNA library (prey) in sequential transformation procedure. The co-transformants were selected after incubation at 30°C for 5 days on synthetic defined (SD) minimal medium lacking tryptophan and leucine. All colonies were pooled from the plates and replated on SD minimal medium lacking tryptophan, leucine, histidine and adenine. Positive clones were analyzed for β-galactosidase activity by colony-filter assays. Prey plasmids of positive clones were rescued using <italic>E. coli </italic>KC8 and subjected to DNA sequencing.</p> <title>Confirmation of protein-protein interaction by co-affinity purification</title> <p>The GI and GII cDNA fragments, encoding for putative prey proteins interacting with the subolesin bait, were amplified by PCR using oligonucleotide primers QEGI5: 5'-GGCCATGGAGGCCGGGATAGGA and QEGI3: 5'-GGAGATCTATGCTCATTAAGACAAATCTC, and QEGII5: 5'-GGCCATGGCCTCAACCAGGCCCACGGACAAACCCCTC and QEGII3: 5'-GGAGATCTGGCGACCGTTTGCCTCATGTC, respectively and the Access RT-PCR system (Promega, Madison, WI, USA). The PCR primers introduced <italic>Nco</italic>I and <italic>Bgl</italic>II restriction sites for cloning into the expression vector pQE-60 fused to a 6xHis tag (Qiagen Inc., Valencia, CA, USA). The full-length cDNA of <italic>I. scapularis </italic>subolesin was inserted into the expression plasmid pFLAG-CTC (Sigma, St. Louis, MO, USA) as described previously [##REF##16005748##8##]. Recombinant subolesin, GI and GII proteins were expressed in <italic>E. coli </italic>JM109 as described previously [##REF##16005748##8##]. For co-affinity purification of subolesin-interacting proteins, the Ni-NTA spin columns (Qiagen) were used following manufacturer's protocol. The supernatant of GI and GII cell lysates were loaded onto the columns. The columns were washed three times with washing buffer (Qiagen) before loading the purified subolesin. The columns were incubated at 4°C for 4 hours to allow protein-protein interaction and then washed again. The proteins were eluted with elution buffer (Qiagen) and the eluted proteins were subjected to 12% sodium dodecyl sulfate-polyacrylamide electrophoresis (SDS-PAGE) and immunoblotted with anti-subolesin antibodies [##REF##16005748##8##]. Controls included affinity purification of GI and GII in the absence of subolesin and loading subolesin onto Ni-NTA columns in the absence of His-tagged GI and GII proteins and in the presence of an unrelated His-tagged protein. GI and GII proteins were detected by Western blot using the HisDetector Western Blot Kit HRP (KPL, Gaithersburg, Maryland, USA).</p> <title>RNA interference in ticks</title> <p>Oligonucleotide primers homologous to <italic>I. scapularis </italic>and <italic>R. microplus </italic>subolesin containing T7 promoters were used for <italic>in vitro </italic>transcription and synthesis of subolesin dsRNA as described previously [##REF##16580098##10##,##REF##17196597##12##], using the Access RT-PCR system (Promega) and the Megascript RNAi kit (Ambion, Austin, TX, USA). For the synthesis of <italic>R. microplus </italic>GI and GII dsRNAs, oligonucleotide primers GI5: 5'-ATGGAGGCCGGGATAGGA and GI3: 5'-TGCTCATTAAGACAAATCTC, and GII5: 5'-GGCCCACGGACAAACCCCTC and GII3: 5'-CGACCGTTTGCCTCATGTC were used, respectively, with the same procedure described above. The dsRNA was purified and quantified by spectrophotometry.</p> <p>Unfed <italic>I. scapularis </italic>female ticks were injected with approximately 0.5 μl dsRNA in the lower right quadrant of the ventral surface of the exoskeleton of ticks [##REF##16580098##10##]. Freshly molted <italic>R. microplus </italic>females were placed on double-sided sticky tape with the ventral side upwards and injected with approximately 0.5 μl dsRNA into the anal aperture using a syringe mounted on a MM3301-M3 micromanipulator (World Precision Instruments (WPI), Berlin, Germany) and connected to an UMPII syringe pump (WPI) [##REF##17196597##12##]. Engorged <italic>R. microplus </italic>females were injected with 5 μl of dsRNA (5 × 10<sup>10</sup>–5 × 10<sup>11 </sup>molecules/μl) in the right spiracular plate within 6 hours after dropping off the host [##REF##17196597##12##]. The injections were done using a 10 μl syringe with a 1 inch, 33 gauge needle (Hamilton, Bonaduz, Switzerland). Control ticks were injected with injection buffer (10 mM Tris-HCl, pH 7, 1 mM EDTA) alone (saline negative control). We demonstrated previously that there is no difference between using an unrelated dsRNA or injection buffer alone for negative control in tick RNAi experiments [##REF##16580098##10##].</p> <p>Two RNAi experiments were conducted with <italic>I. scapularis</italic>. In the first experiment, designed to determine the time after dsRNA injection in which subolesin knockdown occurs, 20 female ticks per group were injected with subolesin dsRNA or injection buffer. The ticks were held in a humidity chamber for 1 day after which they were allowed to feed on a sheep for 10 days. Three ticks from each group were collected at 1, 3, 6, 9 and 11 days post injection (dpi). Only the ticks collected after 10 days of feeding (11 dpi) were fed together with male ticks to evaluate the effect of subolesin knockdown on tick weight. After collection, ticks were weighed and the tick weight was compared between subolesin dsRNA-injected and saline controls by Student's t-test (P = 0.05). Internal organs were then dissected and stored in RNALater (Ambion) for RNA extraction to determine subolesin mRNA levels by real-time RT-PCR. In the second experiment, 100 <italic>I. scapularis </italic>female ticks per group were injected with subolesin dsRNA or injection buffer alone. The ticks were held in a humidity chamber for 1 day and then allowed to feed on a sheep without males for 5 and 8 days (6 and 9 dpi) in groups of 50 dsRNA- and 50 saline-injected ticks each. Ten additional female ticks per group were fed with males for 10 days to compare the weight of replete dsRNA-injected and control ticks after RNAi by Student's t-test with unequal variance (P = 0.05). Unattached ticks were removed two days after infestation. After feeding, ticks were removed from the sheep and dissected for RNA extraction.</p> <p>For RNAi in <italic>R. microplus</italic>, unfed and replete female ticks were injected with GI dsRNA, GII or subolesin dsRNA or injection buffer alone as described previously [##REF##17196597##12##]. In the experiment with unfed ticks, 60–76 freshly molted female ticks were injected and fed on a calf for 15 days with untreated males to determine female tick mortality, weight, oviposition and hatching. Tick mortality was evaluated as the ratio of dead female ticks 15 dpi to the total number of female ticks placed on the calf and was compared between dsRNA and mock-injected control ticks by χ<sup>2</sup>-test as implemented in Mstat 4.01 (α = 0.01). Ticks completing feeding and egg masses oviposited by each tick were weighed individually and the average ± SD were calculated and compared between dsRNA and saline injected control ticks by Student's t-test with unequal variance (P = 0.05). In the experiment with replete female ticks, 6 ticks were used per group. Each tick was weighed before injection and stored individually in 2 ml Eppendorf tubes with pierced lids in an incubator at 27°C and 95% relative humidity to evaluate oviposition and hatching after RNAi. The weight of replete ticks before injection and the eggs mass per tick were compared with the control saline injected ticks by Student's t-test with unequal variance (P = 0.05). mRNA levels of RNAi targeted genes were determined by real-time RT-PCR. For real-time RT-PCR, internal organs were dissected and total RNA was extracted from 6 females of each dsRNA- or saline-injected group after 6 days of feeding and from 100 mg eggs oviposited by injected replete ticks.</p> <title>Suppression-subtractive hybridization (SSH)</title> <p>Total RNA was isolated from pooled guts and salivary glands of <italic>I. scapularis </italic>female ticks injected with subolesin dsRNA or injection buffer at 9 dpi (8 days of feeding) using TriReagent (Sigma) according to manufacturer's instructions. RNA quality was assessed by gel electrophoresis. SSH was performed at Evrogen JCS (Moscow, Russia) as described previously [##REF##17964755##6##,##REF##16487312##28##]. Tester and driver RNAs were subtracted in both directions to construct two SSH libraries enriched for differentially expressed cDNAs in subolesin dsRNA-injected (reverse-subtracted) and saline-injected (forward-subtracted) ticks. Approximately 100 clones from each library were randomly picked up and subjected to differential hybridization with subtracted and non-subtracted probes using the PCR-select differential screening kit (Clontech, Palo Alto, CA, USA), which resulted in &gt; 95% candidate differentially expressed cDNAs.</p> <title>Microarray construction and analysis</title> <p>Tick cDNA fragments (384 from each of the reverse and forward subtracted SSH libraries) were amplified by PCR using pAL-16 vector-specific primers, purified (MultiScreen PCR plates, Millipore, Billerica, MA, USA) and arrayed onto gamma amino propyl silane coated GAPS II slides (Corning, Lowell, MA, USA). Eight pools of 12 clones each from an unsubtracted <italic>I. scapularis </italic>cDNA library [##REF##16487312##28##] and subolesin cDNA were also arrayed and used to validate normalization. Total RNA was prepared from subolesin dsRNA- and saline-injected ticks at 6 and 9 dpi (5 and 8 days of feeding) using the RNeasy Mini Kit (Qiagen) including the on-column DNA digestion with the RNase-free DNase set following manufacturer's instructions. Total RNA (5 μg) were labeled using the 3DNA Array900 kit with Alexa Fluor dyes (Genisphere, Hatfield, PA, USA), Superscript II (Invitrogen, Carlsbad, CA, USA), the supplied formamide-based hybridization buffer and 24 × 60 mm LifterSlips (Erie Scientific, Portmouth, NH, USA) according to the manufacturer's (Genisphere) instructions. Hybridization signals were measured using a ScanArray Express (PerkinElmer, Boston, MA, USA) and the images were processed using GenePix Pro version 4.0 (Axon, Union City, CA, USA). Ratios were calculated as subolesin dsRNA-injected ticks versus saline-injected control ticks. Normalized ratio values obtained for each probe were averaged across 50 biological replicates and four technical replicates and significant differences were defined as displaying an expression fold change greater than 2-fold. All the microarray data were deposited at the NCBI Gene Expression Omnibus (GEO) under the platform accession number <ext-link ext-link-type="gen" xlink:href="GPL6394">GPL6394</ext-link> and the series number <ext-link ext-link-type="gen" xlink:href="GSE10222">GSE10222</ext-link>.</p> <title>Real-time RT-PCR analysis</title> <p>The RNA samples prepared as described above from <italic>I. scapularis </italic>female ticks and from <italic>R. microplus </italic>female ticks and eggs after RNAi experiments were used for real-time RT-PCR analysis. The RNA from <italic>I. scapularis </italic>female ticks injected with an unrelated 4A8 dsRNA [##REF##15824899##29##] was used as a control in the analysis of mRNA levels for selected differentially expressed genes after subolesin RNAi. Two primers were synthesized based on the sequences determined for subolesin, GI, GII and three differentially expressed genes after subolesin RNAi (identical to <italic>I. scapularis </italic>putative secreted salivary WC peptide [Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>] and putative secreted protein [Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>] and <italic>M. rosenbergii </italic>Cu-Zn SOD; [Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>]) and used for real-time RT-PCR analysis of mRNA levels in dsRNA- and control ticks. Real-time RT-PCR was done using the QuantiTec SYBR Green RT-PCR kit (Qiagen, Valencia, CA, USA) and a Bio-Rad iQ5 thermal cycler (Hercules, CA, USA) following manufacturer's recommendations and oligonucleotide primers and PCR conditions described in Table ##TAB##3##4##. mRNA levels were normalized against tick β-actin using the comparative Ct method. mRNA levels were compared between dsRNA-infected and control ticks by Student's t-Test (P = 0.05).</p> <title>Sequence analysis and database search</title> <p>The sequences of positive clones in the yeast two-hybrid screen were sequenced and searched for homology with the PSI-BLAST algorithm through NCBI web site. Proteins with similar domain architectures were characterized with CDART [##REF##12368255##30##,##UREF##1##31##].</p> <p>The tick cDNAs identified in the microarray analysis as down or up-regulated by more than two-fold after subolesin RNAi were sequenced and analyzed. Multiple sequence alignment were performed to exclude vector sequences and to identify redundant (not unique) sequences using the program AlignX (Vector NTI Suite V 8.0, InforMax, Invitrogen, Carisbad, CA, USA) with an engine based on the Clustal W algorithm. Searches for sequence similarity were performed with the BLASTX program against the non-redundant (nr) peptide sequence database and the database of tick specific sequences at NCBI and VectorBase. Protein ontology was determined using the protein reference database [##UREF##2##32##]. The possibility of subolesin RNAi off-target effects were analyzed by searching for exact complementarity between the seed region (bases 2–8) of all possible 20–22 bp subolesin siRNAs and cDNA sequences of tick differentially regulated genes identified in the microarray analysis, including the sequence of 3' UTRs when known [##REF##16489337##23##].</p> <p>The pattern search and predictions for tick GI, GII and subolesin post-translational modifications were done by PIR searching against the PROSITE database [##UREF##3##33##]. The search for human ortholog (C6orf166; Genbank accession number <ext-link ext-link-type="gen" xlink:href="NP_060534">NP_060534</ext-link>) protein-protein interactions was done by STRING [##REF##17098935##34##,##UREF##4##35##].</p> <title>Nucleotide sequence accession numbers</title> <p>The nucleotide sequences of GI, GII and the ESTs reported in this paper have been deposited in the GenBank database under accession numbers [Genbank:<ext-link ext-link-type="gen" xlink:href="EU436162">EU436162</ext-link>, Genbank:<ext-link ext-link-type="gen" xlink:href="EU436163">EU436163</ext-link> and Genbank:<ext-link ext-link-type="gen" xlink:href="FD482610">FD482610</ext-link>–<ext-link ext-link-type="gen" xlink:href="FD482874">FD482874</ext-link>], respectively.</p>
<title>Results</title> <title>Identification of subolesin interaction proteins by yeast two-hybrid screen and co-affinity purification</title> <p>For discovery of subolesin-interacting proteins, subolesin was used as a bait to search for preys in <italic>Rhipicephalus (Boophilus) microplus </italic>using the yeast two-hybrid system. Two sequences, GI and GII, were identified that encoded for candidate subolesin-interacting proteins. These sequences were represented in 50% and 10% of the positive clones, respectively. BLAST analysis of the GI sequences did not result in identity to known sequences, except for the EST910636 from a <italic>R. microplus </italic>cDNA library. However, a transduction/transcription domain was found within the GI open reading frame (ORF) (Fig. ##FIG##0##1##). The GII sequence was 99% identical (Expect = 4e-102) to <italic>Amblyomma </italic>sp. elongation factor-1 alpha (EF-1a; [Genbank:<ext-link ext-link-type="gen" xlink:href="AAK12647">AAK12647</ext-link>]) and to other proteins containing EF1_alpha_II and EF1_alpha_III domains (Fig. ##FIG##1##2##). Both GI and GII sequences contained multiple N-myristoylation and casein kinase II (CK2) and protein kinase C (PKC) phosphorylation sites (Figs. ##FIG##0##1## and ##FIG##1##2##).</p> <p>The proteins encoded by GI and GII were expressed in <italic>Escherichia coli </italic>fused to a 6xHis tag with molecular weights of 10 and 22 kDa, respectively and the interaction of these proteins with subolesin was confirmed by co-affinity purification on Ni-NTA columns (Fig. ##FIG##2##3##). In the absence of either GI or GII recombinant proteins or in the presence of an unrelated His-tagged protein, co-purification of subolesin did not occur, confirming the specificity of the interaction between subolesin and GI or GII (Fig. ##FIG##2##3## and data not shown).</p> <title>Characterization of the effect of GI and GII knockdown on tick survival, feeding, oviposition, hatching and egg development</title> <p>RNAi was used to study the effect of the knockdown of subolesin-interacting protein-coding genes, GI and GII, on female tick survival, feeding, oviposition, hatching and egg development. Subolesin, GI or GII dsRNAs or injection buffer (saline control) were injected into either unfed or replete <italic>R. microplus </italic>female ticks. The injection of GI dsRNA did not affect tick survival, feeding and oviposition, a fact that correlated with the absence of GI knockdown in ticks (Table ##TAB##0##1##). Although GI expression was reduced in tick eggs, it did not affect hatching (Table ##TAB##1##2##). Gene knockdown was demonstrated for GII and subolesin in both ticks and eggs (Tables ##TAB##0##1## and ##TAB##1##2##). GII knockdown produced a phenotype similar to that obtained with the silencing of subolesin expression. The injection of GII dsRNA in unfed ticks resulted in reduced tick survival, feeding, oviposition and hatching when compared to control ticks (Table ##TAB##0##1##). When GII dsRNA was injected into replete ticks, hatching and egg development were affected (Table ##TAB##1##2## and Fig. ##FIG##3##4##). Undifferentiated egg masses were observed in eggs oviposited by replete ticks injected with GII and subolesin dsRNA when compared to saline injected controls (Fig. ##FIG##3##4##).</p> <title>Characterization of the effect of subolesin knockdown on tick gene expression profile</title> <p>The results of the subolesin-protein interactions and RNAi experiments with GI and GII, encoding for subolesin-interacting proteins, suggested that subolesin may be involved in gene expression in ticks. Therefore, the subsequent experiments were directed towards analyzing the effect of subolesin knockdown on the tick gene expression profile.</p> <p>In the first experiment, the kinetics of subolesin silencing by RNAi was investigated in unfed <italic>I. scapularis </italic>female ticks. The results showed that over 90% silencing of subolesin transcription was obtained 6 days post injection (dpi) and continued until at least 11 dpi (Fig. ##FIG##4##5##). Significant differences (P = 0.02) in tick weight between subolesin dsRNA (78 ± 33 mg) and saline injected controls (142 ± 47 mg) were observed in ticks that completed feeding on the host in the presence of males at 11 dpi.</p> <p>In the second experiment, unfed <italic>I. scapularis </italic>female ticks were injected with subolesin dsRNA or injection buffer alone. Ticks were then placed on a sheep without males and collected at 6 and 9 dpi, by which time subolesin knockdown was regarded as significant based on the RNAi kinetics experiment (Fig. ##FIG##4##5##). Subolesin knockdown was corroborated by real-time RT-PCR as 81 ± 3% and 79 ± 5% silencing at 6 and 9 dpi, respectively. A group of ticks collected after 10 days of feeding with males was used to corroborate the effect of subolesin knockdown on tick weight, giving results similar to those obtained in the RNAi kinetics experiment described above. The ticks collected at 6 and 9 dpi were then dissected and used to extract RNA for suppression-subtractive hybridization (SSH) and microarray construction and analysis.</p> <p>The SSH libraries constructed with the RNA of subolesin dsRNA and saline injected ticks collected at 9 dpi were used for random clone amplification and microarray construction. This microarray was enriched for genes differentially expressed after subolesin knockdown. The microarray was hybridized with total RNA prepared from subolesin dsRNA- and saline-injected ticks at 6 and 9 dpi to evaluate how the expression pattern of candidate differentially expressed genes varied at two different time points after subolesin knockdown (Table ##TAB##2##3##). The results revealed 34 unique (not redundant) genes that were down or up-regulated by at least two-fold after subolesin knockdown at 6 or 9 dpi (Fig. ##FIG##5##6##). Of these genes, 28 were down-regulated and 6 were up-regulated after subolesin RNAi. Of the down-regulated genes, 4 were down-reglated at 6 dpi, 20 at 9 dpi and 4 at both 6 and 9 dpi (Fig. ##FIG##5##6##). Up-regulated genes were found at 9 dpi only (Fig. ##FIG##5##6##). In the microarray analysis, subolesin expression was silenced by approximately 35% at both 6 and 9 dpi. Although 62% of the unique differentially expressed genes did not have homology to available sequences or were identical to genes with unknown function, the remaining sequences indicated that various biological processes were affected by subolesin knockdown. These included regulation of protein and nucleic acid metabolism (54% of protein ontology assignments), energy pathways (15%), immunity (15%), cell communication and signal transduction (8%) and transport (8%).</p> <p>Three of the down-regulated genes after subolesin knockdown, identical to <italic>I. scapularis </italic>putative secreted salivary WC peptide [Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>] and putative secreted protein [Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>] and <italic>Macrobrachium rosenbergii </italic>copper/zinc superoxide dismutase (Cu-Zn SOD) [Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>], were selected to corroborate the results of the microarray analysis by real-time RT-PCR. The results demonstrated that the expression of these genes was silenced after subolesin RNAi by 98% and 100% (for [Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>]), 41% and 59% (for [Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>]) and 93% and 74% (for [Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>]) at 6 and 9 dpi, respectively. Ticks injected with an unrelated 4A8 dsRNA control did not show silencing of these genes after RNAi. The possibility of subolesin RNAi off-target effects was analyzed and complementary 7 bp regions were not found between all possible 20–22 bp subolesin siRNAs and tick cDNA sequences of differentially expressed genes identified in the microarray analysis.</p> <title>Prediction of subolesin conserved post-translational modifications</title> <p><italic>I. scapularis </italic>subolesin [Genbank:<ext-link ext-link-type="gen" xlink:href="AAV67031">AAV67031</ext-link>] and human ortholog C6orf166 [Genbank:<ext-link ext-link-type="gen" xlink:href="NP_060534">NP_060534</ext-link>] proteins were compared to predict conserved post-translational modifications. Three conserved PKC phosphorylation sites were found (Fig. ##FIG##6##7##), which were also present in all known tick subolesin protein sequences (data not shown).</p>
<title>Discussion</title> <p>Subolesin, discovered and characterized in <italic>I. scapularis </italic>as a tick protective antigen [##REF##12615446##7##, ####REF##16005748##8##, ##REF##16153760##9####16153760##9##], is an evolutionary conserved protein which is involved in modulation of tick blood digestion, reproduction and development [##REF##16580098##10##, ####REF##16630571##11##, ##REF##17196597##12##, ##REF##17372765##13####17372765##13##]. In other organisms, subolesin orthologs may be involved in the control of developmental processes [##REF##11901120##16##, ####UREF##0##17##, ##REF##11231151##18####11231151##18##]. Although the function of subolesin is unknown, these results suggest a conserved function for subolesin. Because of the profound effect of subolesin knockdown in ticks and other organisms [##REF##16580098##10##, ####REF##16630571##11##, ##REF##17196597##12##, ##REF##17372765##13####17372765##13##,##REF##11231151##18##], our hypothesis was that subolesin may have a role in gene expression, thus affecting multiple cellular processes. Therefore, the objective of this study was to provide evidence of the role of subolesin in gene expression. To test this hypothesis, three experiments were conducted. In the first series of experiments, subolesin-interacting proteins were identified and characterized in <italic>R. microplus</italic>, suggesting the interaction of subolesin with regulatory proteins. Therefore, in the second series of experiments, the effect of subolesin knockdown was analyzed in <italic>I. scapularis </italic>and showed the effect of subolesin on gene expression affecting different biological processes. Finally, post-translational modifications were predicted for tick subolesin. All together, the results of these experiments suggested a role for tick subolesin in gene expression.</p> <p>To identify proteins that interact with subolesin in yeast two-hybrid experiments, we used a cDNA library obtained from tick eggs because subolesin is expressed in tick embryos and gene knockdown affects egg development [##REF##16005748##8##,##REF##17196597##12##]. Two genes, GI and GII, were identified encoding for proteins that interact with tick subolesin. These proteins contained domains and post-translational modification sites found in proteins with regulatory functions. The transduction/transcription domain found in GI is present in phosphorylated proteins involved in transcriptional regulation and other cell functions related to gene expression [##REF##10373431##19##,##REF##15302935##20##]. The EF1_alpha_II and EF1_alpha_III domains present in GII are found in proteins with different functions such as protein biosynthesis, DNA binding, transcriptional regulation, RNA processing, structural constituent of cytoskeleton as well as ATP and GTP binding that are also involved in gene expression (see for example proteins with accession numbers [Genbank:<ext-link ext-link-type="gen" xlink:href="EAY95388">EAY95388</ext-link>, Genbank:<ext-link ext-link-type="gen" xlink:href="XP_657651">XP_657651</ext-link>, Genbank:<ext-link ext-link-type="gen" xlink:href="XP_404184">XP_404184</ext-link>, Genbank:<ext-link ext-link-type="gen" xlink:href="XP_312333">XP_312333</ext-link>, Genbank:<ext-link ext-link-type="gen" xlink:href="XP_001651261">XP_001651261</ext-link>]).</p> <p>The results reported herein suggested that subolesin may interact with regulatory proteins. In other organisms, subolesin orthologs interact with proteins with gene expression regulatory activities. The human subolesin ortholog interacts with LNXp80 [Genbank:<ext-link ext-link-type="gen" xlink:href="AK056823">AK056823</ext-link>], DIPA [Genbank:<ext-link ext-link-type="gen" xlink:href="NM_006848">NM_006848</ext-link>] and SPG21 [Genbank:<ext-link ext-link-type="gen" xlink:href="NM_016630">NM_016630</ext-link>] proteins. LNX is an E3 ubiquitin-protein ligase that mediates ubiquitination and subsequent proteasomal degradation of Numb, implicated in the control of cell fate decisions during development. DIPA interacts with the viral phosphoprotein hepatitis delta antigen (HDAG) and acts as a repressor of gene transcription [##REF##17014843##21##]. SPG21 binds to the hydrophobic C-terminal amino acids of CD4 which are involved in repression of T cell activation [##REF##11113139##22##]. In <italic>D. melanogaster</italic>, the subolesin ortholog <italic>bhringi </italic>(<italic>bhr</italic>; CG8580) may act to regulate Twist activity through recruitment of the chromatin remodeling Brahma complex [##UREF##0##17##]. Therefore, subolesin may exert its effect on gene expression through the interaction with GI, GII and possibly other regulatory proteins. Interestingly, the GII knockdown phenotype was similar to that obtained with subolesin, suggesting that these proteins may functionally interact in ticks.</p> <p>The results of the microarray analysis of gene expression profile in ticks after subolesin knockdown provided evidence for the role for subolesin in gene expression. Subolesin knockdown affected the expression of genes involved in multiple cellular pathways. The nonspecific effect of dsRNA injection on tick global gene expression was not addressed in these studies. However, the injection of an unrelated dsRNA did not affect the expression of selected genes differentially expressed after subolesin knockdown. Although we cannot rule out off-target effects of subolesin RNAi in ticks, evidence suggested that this was not a likely possibility to explain the effect of subolesin knockdown on tick gene expression pattern. Firstly, we did not find complementary sequences between subolesin and identified differentially expressed genes that could support off-target effects of subolesin RNAi. To search for complementary sequences between subolesin and identified differentially expressed genes, we used the approach proposed by Birmingham et al. [##REF##16489337##23##] who showed that although maximum complementarity by itself is an unsatisfactory predictor of off-target RNAi effects, a highly significant association exists between off-targeting and exact complementarity between the seed region (bases 2–8) of siRNA and their off-targeted gene 3' untranslated region (UTR). Secondly, the analysis of <italic>D. melanogaster </italic>subolesin ortholog RNAi off-target effects demonstrated the presence of a single off-targeted gene [##REF##17625558##24##], suggesting that off-target effects of subolesin RNAi may also be minimal in ticks.</p> <p>A common characteristic of many regulatory protein sequences is the presence of phosphorylation sites. Although we did not demonstrate phosphorylation of tick subolesin, there is evidence that the human ortholog protein undergoes phosphorylation at serine 21 (AS*PKRRR) [##REF##17081983##25##], a PKC phosphorylation site that is conserved in tick sequences. Therefore, as with other regulatory proteins, subolesin may be regulated by reversible phosphorylation by PKC.</p>
<title>Conclusion</title> <p>In summary, the results presented herein provide evidence that support a role for subolesin in gene expression in ticks and other organisms. The regulatory function of subolesin in gene expression may be exerted through interaction with other regulatory proteins at the transcriptional level (Fig. ##FIG##7##8##). In fact, a recent publication by Goto et al. [##REF##18066067##26##] that appeared after submission of our work renamed subolesin orthologs in insects and vertebrates as Akirins and proposed that they constitute transcription factors required for NF-kB-dependent gene expression in <italic>Drosophila </italic>and mice. Alternatively, subolesin may also interact with proteins involved in translational control in ticks (Fig. ##FIG##7##8##). The experimental approach used in this study may be important for the annotation of tick sequences that result from genome sequencing efforts [##REF##17656151##27##] and for the characterization of candidate tick protective antigens for the development of vaccines for the control of tick infestations and the transmission of tick-borne pathogens [##REF##16842264##3##].</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Subolesin is an evolutionary conserved protein that was discovered recently in <italic>Ixodes scapularis </italic>as a tick protective antigen and has a role in tick blood digestion, reproduction and development. In other organisms, subolesin orthologs may be involved in the control of developmental processes. Because of the profound effect of subolesin knockdown in ticks and other organisms, we hypothesized that subolesin plays a role in gene expression, and therefore affects multiple cellular processes. The objective of this study was to provide evidence for the role of subolesin in gene expression.</p> <title>Results</title> <p>Two subolesin-interacting proteins were identified and characterized by yeast two-hybrid screen, co-affinity purification and RNA interference (RNAi). The effect of subolesin knockdown on the tick gene expression pattern was characterized by microarray analysis and demonstrated that subolesin RNAi affects the expression of genes involved in multiple cellular pathways. The analysis of subolesin and interacting protein sequences identified regulatory motifs and predicted the presence of conserved protein kinase C (PKC) phosphorylation sites.</p> <title>Conclusion</title> <p>Collectively, these results provide evidence that subolesin plays a role in gene expression in ticks.</p>
<title>Authors' contributions</title> <p>JdlF conceived and coordinated the study, participated in its design and helped to draft the manuscript. CM-O carried out two hybrid screening. VN, PA, AMN, CA, JMPdlL, RCG, EFB and KMK carried out RNAi, microarray analysis and real-time RT-PCR experiments. MC carried out confirmation of protein-protein interactions. PA and JdlF carried out computational analyses. KMK, FJ, CG, PA and CM-O participated in study design and helped to draft the manuscript. All authors read and approved the final manuscript.</p>
<title>Acknowledgements</title> <p>We thank A. Taoufik (UCTD) for technical assistance. This research was supported by the Oklahoma Agricultural Experiment Station (project 1669), the Walter R. Sitlington Endowed Chair for Food Animal Research (K. M. Kocan, Oklahoma State University), Pfizer Animal Health, Kalamazoo, MI, USA, the Junta de Comunidades de Castilla-La Mancha, Spain (project 06036-00 ICS-JCCM), and the Wellcome Trust under the "Animal Health in the Developing World" initiative (project 0757990). V. Naranjo was founded by Consejería de Educación, JCCM, Spain. R.C. Galindo was funded by Ministerio de Educación y Ciencia, Spain.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Analysis of GI sequence encoding for subolesin-interacting protein.</bold> GI nucleotide and deduced amino acid sequences are shown. The GI sequence contains a transduction/transcription domain (boxed letters) and multiple N-myristoylation (solid line) and CK2 (dotted line) and PKC (dashed line) phosphorylation sites.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Analysis of GII sequence encoding for subolesin-interacting protein.</bold> GII nucleotide and deduced amino acid sequences are shown. The GII sequence contains EF1_alpha_II and EF1_alpha_III domains (boxed letters) and multiple N-myristoylation (solid line) and CK2 (dotted line) and PKC (dashed line) phosphorylation sites.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Co-affinity purification assays.</bold> The supernatant of GI- and GII-expressing <italic>E. coli </italic>cell lysates were loaded onto Ni-NTA columns, washed, then loaded with purified recombinant subolesin and washed again before protein elution for SDS-PAGE and Western blot. The top panel shows subolesin bait detection after affinity purification on Ni-NTA columns in the presence (+) or absence (-) of His-tagged GI/GII protein extracts, demonstrating binding to His-tagged prey. The bottom panel shows detection of GI/GII His-prey after affinity purification on Ni-NTA columns.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Representative <italic>R. microplus </italic>eggs from replete female ticks injected with injection buffer alone (saline control) showing normal development, or GII and subolesin dsRNA showing an undifferentiated egg mass.</bold> Eggs were incubated at 27°C and 95% relative humidity. Photographs were taken 22 dpi of replete females and 19 days after oviposition. Bar = 0.1 mm.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Kinetics of subolesin RNAi.</bold> Ticks were injected with subolesin dsRNA (RNAi group) or injection buffer alone (control group). Three ticks from each group were collected at 1, 3, 6, 9 and 11 days post injection to determine mean ± SD subolesin mRNA ratio by real-time RT-PCR.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Effect of subolesin knockdown on tick gene expression pattern.</bold> The expression fold change was determined by microarray hybridization at 6 and 9 days post injection (dpi). Clone ID (SSH library plate and well) are shown and correspond to entries in Table 3. The graph was constructed with the HCE software <ext-link ext-link-type="uri" xlink:href="http://www.cs.umd.edu/hcil/hce/hce3.html"/>.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Prediction of subolesin post-translational modifications.</bold> Sequence alignment of <italic>I. scapularis </italic>subolesin [Genbank:<ext-link ext-link-type="gen" xlink:href="AAV67031">AAV67031</ext-link>] and human ortholog [Genbank:<ext-link ext-link-type="gen" xlink:href="NP_060534">NP_060534</ext-link>] proteins using the one letter amino acid code. Identical amino acids are indicated with asterisks. Three conserved PKC phosphorylation sites were predicted by PIR searching against the PROSITE database.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>The regulatory function of subolesin on gene expression may be exerted through interactions with regulatory proteins that act at the transcriptional and/or translational levels.</bold> Tick subolesin interacts with GI and GII proteins with putative regulatory functions on gene expression. Other hypothetical interactions with regulatory proteins may exist based on data for the human ortholog. The search for human ortholog (C6orf166; [Genbank:<ext-link ext-link-type="gen" xlink:href="NP_060534">NP_060534</ext-link>]) protein-protein interactions was done by STRING.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p><italic>R. microplus </italic>tick survival, weight, oviposition and hatching after RNAi in unfed female ticks.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Experimental group</bold></td><td align="left"><bold>Number of ticks</bold></td><td align="left"><bold>Expression silencing (%)<sup>a</sup></bold></td><td align="left"><bold>Tick weight (mg)<sup>b</sup></bold></td><td align="left"><bold>Mortality (%)<sup>c</sup></bold></td><td align="left"><bold>Eggs per tick (mg)<sup>d</sup></bold></td><td align="left"><bold>Hatching rate (%)</bold></td></tr></thead><tbody><tr><td align="left">Saline control</td><td align="left">76</td><td align="left">---</td><td align="left">295 ± 85</td><td align="left">21</td><td align="left">111 ± 75</td><td align="left">&gt; 90</td></tr><tr><td align="left">GI</td><td align="left">60</td><td align="left">-44 ± 29</td><td align="left">266 ± 114</td><td align="left">17</td><td align="left">110 ± 72</td><td align="left">&gt; 90</td></tr><tr><td align="left">GII</td><td align="left">63</td><td align="left">98 ± 0.4*</td><td align="left">131 ± 161**</td><td align="left">60***</td><td align="left">1 ± 4****</td><td align="left">0*****</td></tr><tr><td align="left">Subolesin</td><td align="left">66</td><td align="left">91 ± 3*</td><td align="left">73 ± 94**</td><td align="left">80***</td><td align="left">0****</td><td align="left">ND</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Weight, oviposition and hatching of replete <italic>R. microplus </italic>female ticks after RNAi.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Experimental group</bold></td><td align="left"><bold>Tick weight (mg)<sup>a</sup></bold></td><td align="left"><bold>Eggs per tick (mg)<sup>b</sup></bold></td><td align="left"><bold>Hatching rate (%)<sup>c</sup></bold></td><td align="left"><bold>Expression silencing (%)<sup>d</sup></bold></td></tr></thead><tbody><tr><td align="left">Saline control</td><td align="left">318 (272–367)</td><td align="left">173 (128–229)</td><td align="left">&gt; 90</td><td align="left">---</td></tr><tr><td align="left">GI</td><td align="left">318 (260–369)</td><td align="left">123 (31–175)</td><td align="left">&gt; 90</td><td align="left">55 ± 8**</td></tr><tr><td align="left">GII</td><td align="left">315 (267–349)</td><td align="left">132 (95–182)</td><td align="left">0.2*</td><td align="left">84 ± 5**</td></tr><tr><td align="left">Subolesin</td><td align="left">318 (270–360)</td><td align="left">142 (101–175)</td><td align="left">0.4*</td><td align="left">46 ± 8**</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p><italic>I. scapularis </italic>differentially expressed genes after subolesin knockdown.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Clone ID<sup>a</sup></bold></td><td align="left"><bold>Accession number, name and species<sup>b</sup></bold></td><td align="left"><bold>Fold change (6 dpi)<sup>c</sup></bold></td><td align="left"><bold>SD (6 dpi)<sup>d</sup></bold></td><td align="left"><bold>Fold change (9 dpi)<sup>c</sup></bold></td><td align="left"><bold>SD (9 dpi)<sup>d</sup></bold></td></tr></thead><tbody><tr><td align="left">LibPlateC1_wellB12</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAT92189">AAT92189</ext-link>] KUN-6 (<italic>Ixodes pacificus</italic>)</td><td align="left">-1.1196</td><td align="left">0.2519</td><td align="left">-2.2057</td><td align="left">0.3109</td></tr><tr><td align="left">LibPlateC1_wellB6</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>] copper/zinc superoxide dismutase (<italic>Macrobrachium rosenbergii</italic>)</td><td align="left">-2.7645</td><td align="left">0.5455</td><td align="left">-2.5136</td><td align="left">0.2930</td></tr><tr><td align="left">LibPlateC1_wellD9</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66901">AAY66901</ext-link>] histone 2B (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.5922</td><td align="left">1.5332</td><td align="left">-2.5298</td><td align="left">0.9579</td></tr><tr><td align="left">LibPlateC1_wellE6</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>] putative secreted salivary WC peptide (<italic>Ixodes scapularis</italic>)</td><td align="left">1.2320</td><td align="left">1.4160</td><td align="left">-5.2598</td><td align="left">1.3686</td></tr><tr><td align="left">LibPlateC1_wellF10</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66817">AAY66817</ext-link>] anticoagulant Salp9-like (<italic>Ixodes scapularis</italic>)</td><td align="left">-2.0332</td><td align="left">0.2959</td><td align="left">-1.2037</td><td align="left">0.1124</td></tr><tr><td align="left">LibPlateC1_wellG12</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>] copper/zinc superoxide dismutase (<italic>Macrobrachium rosenbergii</italic>)</td><td align="left">-1.6426</td><td align="left">0.2527</td><td align="left">-2.2118</td><td align="left">0.4974</td></tr><tr><td align="left">LibPlateC1_wellH10</td><td align="left">No homology found</td><td align="left">-1.8346</td><td align="left">0.7216</td><td align="left">-2.0350</td><td align="left">1.0821</td></tr><tr><td align="left">LibPlateC1_wellH4</td><td align="left">No homology found</td><td align="left">-1.6161</td><td align="left">0.4640</td><td align="left">-2.0602</td><td align="left">0.5671</td></tr><tr><td align="left">LibPlateC2_wellB11</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>] putative secreted salivary WC peptide (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.8612</td><td align="left">0.1952</td><td align="left">-2.1981</td><td align="left">0.1664</td></tr><tr><td align="left">LibPlateC2_wellB12</td><td align="left">No homology found</td><td align="left">-1.1313</td><td align="left">0.1588</td><td align="left">-2.1133</td><td align="left">0.2023</td></tr><tr><td align="left">LibPlateC2_wellD5</td><td align="left">No homology found</td><td align="left">-1.6882</td><td align="left">0.0965</td><td align="left">-2.2689</td><td align="left">0.2158</td></tr><tr><td align="left">LibPlateC2_wellE3</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>] putative secreted salivary WC peptide (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.0794</td><td align="left">0.2033</td><td align="left">-2.2673</td><td align="left">0.2302</td></tr><tr><td align="left">LibPlateC2_wellG1</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>|<ext-link ext-link-type="gen" xlink:href="AF483711_1">AF483711_1</ext-link>] putative secreted protein (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.4409</td><td align="left">0.2000</td><td align="left">-3.5351</td><td align="left">0.3250</td></tr><tr><td align="left">LibPlateC3_wellB10</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAH67494">AAH67494</ext-link>] HIST1H3I protein (<italic>Homo sapiens</italic>)</td><td align="left">-2.0111</td><td align="left">0.2174</td><td align="left">-2.3299</td><td align="left">0.4323</td></tr><tr><td align="left">LibPlateC3_wellB8</td><td align="left">No homology found</td><td align="left">-2.0438</td><td align="left">1.5307</td><td align="left">-1.3772</td><td align="left">0.5537</td></tr><tr><td align="left">LibPlateC3_wellC7</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>] copper/zinc superoxide dismutase (<italic>Macrobrachium rosenbergii</italic>)</td><td align="left">-4.1699</td><td align="left">0.1770</td><td align="left">-3.1053</td><td align="left">0.0638</td></tr><tr><td align="left">LibPlateC3_wellD7</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>] copper/zinc superoxide dismutase (<italic>Macrobrachium rosenbergii</italic>)</td><td align="left">-2.5487</td><td align="left">0.1977</td><td align="left">-2.1210</td><td align="left">0.4374</td></tr><tr><td align="left">LibPlateC3_wellE2</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="XP_624446.2">XP_624446.2</ext-link>] Predicted similar to ruby CG11427-PA isoform 2 (<italic>Apis mellifera</italic>)</td><td align="left">-2.4606</td><td align="left">1.6665</td><td align="left">-1.0115</td><td align="left">0.3122</td></tr><tr><td align="left">LibPlateC3_wellE9</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>] copper/zinc superoxide dismutase (<italic>Macrobrachium rosenbergii</italic>)</td><td align="left">-3.4248</td><td align="left">0.3717</td><td align="left">-2.0792</td><td align="left">0.2988</td></tr><tr><td align="left">LibPlateC3_wellF10</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>|<ext-link ext-link-type="gen" xlink:href="AF483711_1">AF483711_1</ext-link>] putative secreted protein (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.3583</td><td align="left">0.2449</td><td align="left">-3.9731</td><td align="left">0.1679</td></tr><tr><td align="left">LibPlateC3_wellF11</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>|<ext-link ext-link-type="gen" xlink:href="AF483711_1">AF483711_1</ext-link>] putative secreted protein (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.0215</td><td align="left">0.1441</td><td align="left">-2.9246</td><td align="left">0.1172</td></tr><tr><td align="left">LibPlateC3_wellG5</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>|<ext-link ext-link-type="gen" xlink:href="AF483711_1">AF483711_1</ext-link>] putative secreted protein (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.1451</td><td align="left">0.2311</td><td align="left">-3.7711</td><td align="left">0.1890</td></tr><tr><td align="left">LibPlateC4_wellB1</td><td align="left">No homology found</td><td align="left">-1.7450</td><td align="left">0.6400</td><td align="left">-2.2268</td><td align="left">0.3409</td></tr><tr><td align="left">LibPlateC4_wellB10</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAQ01562">AAQ01562</ext-link>] von Willebrand factor (<italic>Ixodes ricinus</italic>)</td><td align="left">1.1202</td><td align="left">0.3737</td><td align="left">-3.9449</td><td align="left">0.2800</td></tr><tr><td align="left">LibPlateC4_wellB5</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="ABD83654">ABD83654</ext-link>] hemelipoglycoprotein precursor (<italic>Dermacentor variabilis</italic>)</td><td align="left">-1.9958</td><td align="left">0.5689</td><td align="left">-3.5813</td><td align="left">0.6757</td></tr><tr><td align="left">LibPlateC4_wellB7</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAP84098">AAP84098</ext-link>] ML domain-containing protein (<italic>Ixodes ricinus</italic>)</td><td align="left">-1.5703</td><td align="left">0.3164</td><td align="left">-4.5599</td><td align="left">0.4105</td></tr><tr><td align="left">LibPlateC4_wellB8</td><td align="left">No homology found</td><td align="left">-2.6551</td><td align="left">0.3564</td><td align="left">-3.5944</td><td align="left">0.3325</td></tr><tr><td align="left">LibPlateC4_wellC1</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66502">AAY66502</ext-link>] secreted salivary gland protein (<italic>Ixodes scapularis</italic>)</td><td align="left">1.0009</td><td align="left">0.2427</td><td align="left">-2.2922</td><td align="left">0.3518</td></tr><tr><td align="left">LibPlateC4_wellC2</td><td align="left">No homology found</td><td align="left">-1.0286</td><td align="left">0.0610</td><td align="left">-2.0853</td><td align="left">0.1338</td></tr><tr><td align="left">LibPlateC4_wellC6</td><td align="left">No homology found</td><td align="left">-1.1217</td><td align="left">0.1682</td><td align="left">-2.7104</td><td align="left">0.1235</td></tr><tr><td align="left">LibPlateC4_wellD3</td><td align="left">No homology found</td><td align="left">1.6859</td><td align="left">0.2075</td><td align="left">2.6468</td><td align="left">0.1954</td></tr><tr><td align="left">LibPlateC4_wellD4</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="XP_974124">XP_974124</ext-link>] PREDICTED: similar to CG2972-PA (<italic>Tribolium castaneum</italic>)</td><td align="left">-2.1335</td><td align="left">0.2315</td><td align="left">-2.6578</td><td align="left">0.3642</td></tr><tr><td align="left">LibPlateC4_wellF6</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="BAE53722">BAE53722</ext-link>] aspartic protease (<italic>Haemaphysalis longicornis</italic>)</td><td align="left">-1.4881</td><td align="left">0.2385</td><td align="left">-4.8576</td><td align="left">0.3690</td></tr><tr><td align="left">LibPlateC4_wellG3</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66629">AAY66629</ext-link>] putative secreted salivary protein (<italic>Ixodes scapularis</italic>)</td><td align="left">1.0601</td><td align="left">0.1625</td><td align="left">2.0049</td><td align="left">0.1481</td></tr><tr><td align="left">LibPlateC4_wellG4</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="XP_794044">XP_794044</ext-link>] Predicted similar to Endoplasmic reticulum-golgi intermediate</td><td align="left">-2.2835</td><td align="left">0.2324</td><td align="left">-1.7916</td><td align="left">0.5510</td></tr><tr><td align="left">LibPlateC4_wellH4</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66660">AAY66660</ext-link>] putative salivary secreted protein (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.5502</td><td align="left">0.2177</td><td align="left">-5.1239</td><td align="left">0.2056</td></tr><tr><td align="left">LibPlateC4_wellH6</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66982">AAY66982</ext-link>] cyclophilin A (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.0545</td><td align="left">0.0098</td><td align="left">2.0388</td><td align="left">0.2239</td></tr><tr><td align="left">LibPlateR1_wellG3</td><td align="left">No homology found</td><td align="left">1.2792</td><td align="left">0.2664</td><td align="left">2.2118</td><td align="left">0.1902</td></tr><tr><td align="left">LibPlateR1_wellH7</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66764">AAY66764</ext-link>] putative secreted salivary protein (<italic>Ixodes scapularis</italic>)</td><td align="left">1.2470</td><td align="left">0.0447</td><td align="left">-2.0598</td><td align="left">0.0708</td></tr><tr><td align="left">LibPlateR2_wellB2</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66713">AAY66713</ext-link>] putative secreted salivary protein (<italic>Ixodes scapularis</italic>)</td><td align="left">1.4186</td><td align="left">0.2490</td><td align="left">2.2141</td><td align="left">0.2145</td></tr><tr><td align="left">LibPlateR2_wellF10</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAK97818">AAK97818</ext-link>|<ext-link ext-link-type="gen" xlink:href="AF209915_1">AF209915_1</ext-link>] 16 kDa salivary gland protein A (<italic>Ixodes scapularis</italic>)</td><td align="left">-1.8100</td><td align="left">0.4667</td><td align="left">-2.4275</td><td align="left">0.3690</td></tr><tr><td align="left">LibPlateR2_wellG7</td><td align="left">[Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66713">AAY66713</ext-link>] putative secreted salivary protein (<italic>Ixodes scapularis</italic>)</td><td align="left">1.4101</td><td align="left">0.1056</td><td align="left">2.1495</td><td align="left">0.1347</td></tr><tr><td align="left">LibPlateR3_wellA7</td><td align="left">No homology found</td><td align="left">-1.0338</td><td align="left">0.1570</td><td align="left">2.0213</td><td align="left">0.1464</td></tr><tr><td align="left">LibPlateR4_wellA3</td><td align="left">No homology found</td><td align="left">-1.0990</td><td align="left">1.1277</td><td align="left">-3.4438</td><td align="left">1.6338</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>RT-PCR oligonucleotide primers and conditions.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Gene description</bold></td><td align="left"><bold>Upstream/downstream primer sequences (5'-3')</bold></td><td align="left"><bold>PCR annealing conditions</bold></td></tr></thead><tbody><tr><td align="left"><bold><italic>R. microplus</italic></bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">GI [Genbank:<ext-link ext-link-type="gen" xlink:href="EU436162">EU436162</ext-link>]</td><td align="left">CACCATCACTGAAAGCGGT<break/>GTGTTAAATAGTTATGCTCATTAAGAC</td><td align="left">50°C, 30 sec</td></tr><tr><td align="left">GII [Genbank:<ext-link ext-link-type="gen" xlink:href="EU436163">EU436163</ext-link>]</td><td align="left">GATCATTGACCTGGTACCTTCC<break/>GACTTGATGACACCGACGG</td><td align="left">50°C, 30 sec</td></tr><tr><td align="left">Subolesin [Genbank:<ext-link ext-link-type="gen" xlink:href="DQ923495">DQ923495</ext-link>]</td><td align="left">CACAGTCCGAGTGGCAGAT<break/>GATGCACTGGTGACGAGAGA</td><td align="left">50°C, 30 sec</td></tr><tr><td align="left">Beta actin [Genbank:<ext-link ext-link-type="gen" xlink:href="AY255624">AY255624</ext-link>]</td><td align="left">CAC GGT ATC GTC ACC AAC TG<break/>TGA TCT GCG TCA TCT TCT CG</td><td align="left">50°C, 30 sec</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"><bold><italic>I. scapularis</italic></bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">Putative secreted salivary WC peptide [Genbank:<ext-link ext-link-type="gen" xlink:href="AAY66498">AAY66498</ext-link>]</td><td align="left">GATATTGATCCAGCCGGAGA<break/>ATGTCGTCCCTCCATTGTGT</td><td align="left">60°C, 30 sec</td></tr><tr><td align="left">Putative secreted protein [Genbank:<ext-link ext-link-type="gen" xlink:href="AAM93633">AAM93633</ext-link>]</td><td align="left">TGAAGGCAACCATTGCAGTT<break/>ATTGATGGCAATCCTGTGGA</td><td align="left">56°C, 30 sec</td></tr><tr><td align="left">Identical to <italic>M. rosenbergii </italic>Cu-Zn SOD [Genbank:<ext-link ext-link-type="gen" xlink:href="AAZ29240">AAZ29240</ext-link>]</td><td align="left">TGACCTGGGCAACGTTGA<break/>ATGACGCAGCAGGCAATG</td><td align="left">54°C, 30 sec</td></tr><tr><td align="left">Subolesin [Genbank:<ext-link ext-link-type="gen" xlink:href="AY652654">AY652654</ext-link>]</td><td align="left">AGCAGCTCTGCTTCTCGTCT<break/>TCGTACTCGTCGCGTATCTG</td><td align="left">54°C, 30 sec</td></tr><tr><td align="left">Beta actin [Genbank:<ext-link ext-link-type="gen" xlink:href="AF426178">AF426178</ext-link>]</td><td align="left">GAGAAGATGACCCAGATCA<break/>GTTGCCGATGGTGATCACC</td><td align="left">50°C, 30 sec</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p><sup>a</sup>The expression silencing of target genes was determined by real-time RT-PCR and average ± SD mRNA levels calculated and compared between dsRNA and saline-injected control ticks by Student's t-test with unequal variance (*P &lt; 0.01). Amplification efficiencies were normalized against β-actin.</p><p><sup>b</sup>Ticks which completed feeding and those removed after 30 days (15 days after saline or dsRNA injection) were weighed individually and average ± S.D. calculated and compared between dsRNA and saline-injected control ticks by Student's t-test with unequal variance (**P &lt; 0.01).</p><p><sup>c</sup>Tick mortality was evaluated as the ratio of dead female ticks after day 30 (15 days after saline or dsRNA injection) to the total number of female ticks placed on the animal and was compared between dsRNA and saline-injected control ticks by χ<sup>2</sup>-test (***α &lt; 0.001).</p><p><sup>d</sup>The egg mass oviposited by each tick was weighed individually and average ± S.D. calculated and compared between dsRNA and saline-injected control ticks by Student's t-test with unequal variance (****P &lt; 0.01).</p><p><sup>e</sup>The hatching rate was determined 6 weeks after oviposition and compared with the saline-injected control ticks by Student's t-test with unequal variance (*****P &lt; 0.01).</p><p>Abbreviations: ND, not determined because ticks did not laid eggs.</p></table-wrap-foot> <table-wrap-foot><p><sup>a</sup>Replete <italic>R. microplus </italic>ticks were collected and weighed individually before injection with dsRNA, all within 6 hours after dropping of the host. The average weight and variation (between parentheses) of each group are shown. No significant statistical differences were observed (P &gt; 0.05).</p><p><sup>b</sup>The egg mass oviposited by each tick was weighed individually. The average egg mass and variation (between parentheses) was calculated and compared with the saline-injected control ticks by Student's t-test with unequal variance. No significant statistical differences were observed (P &gt; 0.05).</p><p><sup>c</sup>The hatching rate was determined 6 weeks post oviposition and compared with the saline-injected control ticks by Student's t-test with unequal variance (*P &lt; 0.01).</p><p><sup>d</sup>The expression silencing of target genes was determined in eggs by real-time RT-PCR and average ± S.D. mRNA levels calculated and compared between dsRNA and saline-injected control ticks by Student's t-test with unequal variance (**P &lt; 0.05). Amplification efficiencies were normalized against β-actin.</p></table-wrap-foot> <table-wrap-foot><p><sup>a</sup>Clone ID identifies the clone based on SSH Library Plate_well.</p><p><sup>b</sup>Name: This output only contains descriptions of the top blast hit (most significant alignment based on E value).</p><p><sup>c</sup>Fold change is of the global normalized ratio (log2(635/532)) of background-corrected means averaged between replicates (determined from valid spots only). Only genes down-regulated (negative values) and up-regulated (positive values) after subolesin knockdown by at least two fold at 6 or 9 dpi were considered.</p><p><sup>d</sup>SD is the standard deviation determined from the normalized average log2 ratio (determined from valid spots only).</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-372-1"/> <graphic xlink:href="1471-2164-9-372-2"/> <graphic xlink:href="1471-2164-9-372-3"/> <graphic xlink:href="1471-2164-9-372-4"/> <graphic xlink:href="1471-2164-9-372-5"/> <graphic xlink:href="1471-2164-9-372-6"/> <graphic xlink:href="1471-2164-9-372-7"/> <graphic xlink:href="1471-2164-9-372-8"/>
[{"surname": ["Gonzalez", "Baylies"], "given-names": ["K", "M"], "article-title": ["Bhringi: A novel Twist co-regulator"], "source": ["A Dros Res Conf"], "year": ["2005"], "volume": ["46"], "fpage": ["320B"]}, {"article-title": ["Conserved Domain Architecture Retrieval Tool"]}, {"article-title": ["Protein reference database"]}, {"article-title": ["Protein Information Resource"]}, {"article-title": ["Search Tool for the Retrieval of Interacting Genes/Proteins"]}]
BMC Genomics. 2008 Aug 2; 9:372