- Open Access
The tick endosymbiont Candidatus Midichloria mitochondrii and selenoproteins are essential for the growth of Rickettsia parkeri in the Gulf Coast tick vector
Microbiomevolume 6, Article number: 141 (2018)
Pathogen colonization inside tick tissues is a significant aspect of the overall competence of a vector. Amblyomma maculatum is a competent vector of the spotted fever group rickettsiae, Rickettsia parkeri. When R. parkeri colonizes its tick host, it has the opportunity to dynamically interact with not just its host but with the endosymbionts living within it, and this enables it to modulate the tick’s defenses by regulating tick gene expression. The microbiome in A. maculatum is dominated by two endosymbiont microbes: a Francisella-like endosymbiont (FLE) and Candidatus Midichloria mitochondrii (CMM). A range of selenium-containing proteins (selenoproteins) in A. maculatum ticks protects them from oxidative stress during blood feeding and pathogen infections. Here, we investigated rickettsial multiplication in the presence of tick endosymbionts and characterized the functional significance of selenoproteins during R. parkeri replication in the tick.
FLE and CMM were quantified throughout the tick life stages by quantitative PCR in R. parkeri-infected and uninfected ticks. R. parkeri infection was found to decrease the FLE numbers but CMM thrived across the tick life cycle. Our qRT-PCR analysis indicated that the transcripts of genes with functions related to redox (selenogenes) were upregulated in ticks infected with R. parkeri. Three differentially expressed proteins, selenoprotein M, selenoprotein O, and selenoprotein S were silenced to examine their functional significance during rickettsial replication within the tick tissues. Gene silencing of the target genes was found to impair R. parkeri colonization in the tick vector. Knockdown of the selenogenes triggered a compensatory response from other selenogenes, as observed by changes in gene expression, but oxidative stress levels and endoplasmic reticulum stress inside the ticks were also found to have heightened.
This study illustrates the potential of this new research model for augmenting our understanding of the pathogen interactions occurring within tick hosts and the important roles that symbionts and various tick factors play in regulating pathogen growth.
Ticks are blood-feeding ectoparasites of both humans and animals and are important from a public health perspective because they serve as competent vectors of various disease-causing infectious agents. Many tick-borne pathogens are equipped to infect various tick organs where they can multiply. Infection of the salivary glands enables tick pathogens to readily infect vertebrate hosts upon tick feeding. The spotted fever group rickettsial (SFGR) agent, Rickettsia parkeri, is maintained in tick populations through transstadial (between life stage molts) and transovarial or vertical transmission (deposition into eggs for next-generation pathogen development) . A. maculatum, the Gulf Coast tick, is an arthropod vector with increasing public health significance because of its role as the primary vector of R. parkeri in the USA . Rickettsial diseases are caused by obligate intracellular Gram-negative bacteria, and these organisms infect humans on all continents except Antarctica [3,4,5]. In modern times, the rate and ease of global movement has increased the risk of transporting ticks and tick-borne diseases that may have previously been restricted to one region.
The microorganisms that occupy a tick vector are collectively called the tick microbiome; however, the collection of commensal, symbiotic, and pathogenic microbes associated with ticks are more specifically termed the “pathobiome.” Although the inclusion of all three microbial types appears to be counterintuitive at first sight, it is possible that microbes living in association with pathogens within ticks can positively influence pathogen transmission. For instance, rickettsial endosymbionts are thought to alter the transmission of other rickettsial pathogens, as shown by the inverse relationship between the infection prevalence of Rickettsia rickettsii (pathogen) and R. peacockii (symbiont) in the Rocky Mountain wood tick Dermacentor andersoni [6, 7]. Likewise, the presence of Coxiella-related symbionts in the salivary glands of A. americanum ticks has been proposed to impair the transmission of Ehrlichia chaffeensis . In addition to symbionts, ticks maintain a natural bacterial flora predominantly composed of Proteobacteria, Firmicutes, and Bacteroides phyla [9,10,11,12,13,14], which have also been implicated in pathogen maintenance interference in the tick. For example, when Ixodes scapularis ticks are hatched and raised in a sterile environment, their microbiota is altered such that they experience impaired gut integrity and reduced colonization ability towards Borrelia burgdorferi . As seen with other arthropod vectors, altering the tick microbiome may also result in a modulated type of immune response that can interfere with pathogen survival and infection .
A dynamic interaction happens between tick vectors and their associated disease-causing agents, and this has been referred to as a continuous “bellum omnium contra omes” or war of all against all . An unavoidable interaction between a pathogen and the obligate symbiont(s) in a vector occurs during colonization and transmission. However, understanding about the interactions between rickettsial endosymbionts and pathogenic bacteria in ticks and how they influence each other is limited. There are few published reports on the roles played by symbionts in ticks and whether these bacteria have an impact on tick proliferation or transmission [18, 19]. The symbionts commonly associated with hard ticks belong to Rickettsia, Francisella, Coxiella, Wolbachia, and Candidatus Midichloria genera . Francisella-like endosymbionts (FLE), which have been detected in many ticks [21,22,23], are γ-proteobacterial symbionts and are related to the bacterium that causes tularemia, Francisella tularensis . Genetically distinct FLEs have been reported in D. variabilis and D. andersoni  and across the tick’s developmental stages . Gerhart et al.  hypothesized that pathogenic F. tularensis was capable of transforming into symbiotic FLE in ticks. Candidatus Midichloria mitochondrii (CMM), an α-proteobacterial symbiont first detected in I. ricinus, has a unique intramitochondrial lifestyle . Based on phylogenetic and statistical studies of the 16S rRNA sequences from Midichloria and “similar organisms,” CMM is proposed to belong to a novel family known as “Candidatus Midichloriaceae”  and is widespread in various ixodid ticks . However, our understanding of the interactions between endosymbionts (FLE, CMM) and pathogenic bacteria (R. parkeri) in tick tissues and how they influence each other remains limited.
In the absence of preventive measures, the increasing number of tick-borne diseases poses a significant threat to public health. To survive, ticks must maintain homeostasis (stable equilibrium maintained by physiological processes) and obtain gigantic blood meals of up to 100 times their unfed weight. Selenium (Se) is an essential trace element that is incorporated as selenocysteine (Sec) into selenoproteins (SELENO), many of which form an essential line of defense against oxidative stress damage . These proteins are also responsible for myriad other functions including Se transport, protein folding, and endoplasmic reticulum-associated degradation (ERAD). The endoplasmic reticulum (ER) is involved in intracellular signaling, protein synthesis and protein folding, glycosylation, and secretion of saliva via the exocytotic pathway. Tick saliva composition, as revealed by our sialotranscriptome (from the Greek, sialo means saliva), indicated the presence of over 5000 putative secreted peptides representative of dozens of protein families [32, 33]. Protein folding is dependent on the oxidation of disulfide bridges via reactive oxygen species (ROS). A heightened oxidative environment can impair protein folding, leading to the accumulation of unfolded or misfolded proteins and ultimately ER stress. ER homeostasis can be disrupted by a variety of insults such as the accumulation of misfolded proteins, elevated levels of ROS, pathogen infections, and abnormalities in Ca+ 2 signaling. These disturbances are able to trigger the unfolded-protein response (UPR), a protective counter-measure that acts to reestablish homeostatic balance and promote survival by increasing the production of the chaperones involved in protein folding or by inhibiting global translation and eliminating chronically misfolded proteins. Proteins that fail to properly fold are eliminated via ERAD. ER-resident selenoproteins play a critical role in modulating oxidative and ER stress during prolonged tick feeding on the host. We have discovered multiple factors involved in the synthesis of the tick selenoproteome (a full set of novel selenoproteins in ticks), including a novel eukaryotic elongation factor (eEFSec), a novel SECIS-binding protein (SBP2), and Sec-tRNASec [31,32,33,34]. The serendipitous RNA-Seq findings from our experimental gene silencing of eEFSec indicated that dramatic changes in the expression patterns of the transcripts encoding secreted salivary proteins had occurred . Our further studies revealed that selenoproteins and antioxidants participate in SFGR colonization within the tick vector and in their vertical transmission to the next generation [31, 34,35,36,37,38]. Interestingly, I. scapularis Salp25D (a glutathione peroxidase) and D. variabilis SELENOM confer a survival advantage on B. burgdorferi  and Anaplasma marginale .
In the present study, we used two approaches to gain a better understanding of the replication and physiology of A. maculatum R. parkeri infected (Rp+) and R. parkeri-free (Rp−) colonies isolated from field collection and continuously propagated at The University of Southern Mississippi, USA. First, we examined the potential interplay between pathogenic Rickettsia (R. parkeri) and the dominant non-pathogenic tick symbionts by quantifying the FLE and CMM symbiont loads with or without rickettsial infections. Second, we investigated the differential gene expression of specific tick selenogenes in R. parkeri-infected tissues. Third, we utilized an RNA interference approach to deplete the expression of differentially regulated SELENOM, SELENOO, and SELENOS selenogenes and assess their functional significance in pathobiome maintenance in the tick vector. Overall, we have shown that R. parkeri replication success is correlated with the quantity of CMM present in the tick at the expense of FLE in the tick, and that selenogenes play important roles in tick–pathogen interactions.
Quantitation of R. parkeri, FLE, and CMM across the tick life cycle
R. parkeri-infected (Rp+) and R. parkeri-free (Rp−) A. maculatum tick colonies were established and maintained in the laboratory for studying the dynamics of symbiont–Rickettsia interactions. Our previous microbiome analysis of A. maculatum identified FLE and R. parkeri in A. maculatum ticks . Recently, our Illumina-sequenced sialotranscriptome work has also detected significant numbers of reads from CMM in A. maculatum (Shahid Karim, unpublished results).
We observed that R. parkeri, together with the FLE and CMM intracellular symbionts, are transovarially and transstadially transmitted in Rp+ A. maculatum ticks (Figs. 1, 2, and 3). We also found that the relative concentrations of R. parkeri increased across the tick developmental stages from eggs to adults (Fig. 1). The infected eggs that hatched into infected larvae contained similar concentrations of R. parkeri, but the R. parkeri load significantly increased (> 3.5 fold) with the blood meal in the fed larvae (p = 0.0041). R. parkeri was transstadially transmitted to the nymphal stages and after the nymphal blood meal the R. parkeri load significantly increased (> 2 fold) (p = 0.01). The concentrations of R. parkeri in molted adult males and females were only slightly greater than in the fed nymphs. In the adult ticks, the R. parkeri load decreased in the midgut tissues and the salivary glands upon prolonged blood feeding on the host (Fig. 1b).
The relative concentrations of FLE and CMM in the total bacterial load were estimated in the Rp+ and Rp− ticks to assess the potential interplay that may occur during bacterial replication among R. parkeri, FLE, and CMM (Figs. 2 and 3) across the immature and mature developmental stages of the ticks. The total bacterial concentration was significantly reduced in the Rp− eggs compared with the Rp+ eggs (Fig. 2a–c). In the unfed Rp+ larvae, the total bacterial concentration was lower than in the unfed Rp− ticks, but a significant effect was not observed for FLE or CMM (Fig. 2d–f). However, the larval blood meal in the Rp+ ticks greatly enhanced the total bacterial load and the FLE and CMM concentrations relative to those in the Rp− ticks (Fig. 3d–f). The outcome of this experiment in the nymphs did not parallel those for the larval stage (Fig. 3g–i). Indeed, the Rp+ ticks had reduced FLE loads in the fed and unfed nymphs as compared with the Rp− ticks. Similarly, the total bacterial and CMM loads were seen to increase in the Rp+ nymphs compared with the Rp− nymphs regardless of feeding status (Fig. 3g–i).
The presence of an R. parkeri infection increased the total bacterial load in the Rp+-infected female gut tissues, salivary glands, and ovarian tissues (Fig. 2). The total bacterial load  in naïve ticks (from the Oklahoma State Tick Rearing Facility), along with the FLE levels, were found to decrease rapidly in both tick midgut and salivary gland tissues after the blood meal in the Rp− female adults (Additional file 1: Figure S1), whereas the CMM level remained fairly constant in the midgut tissues over the course of the blood meal but decreased rapidly in the salivary glands after the blood meal (Additional file 1: Figure S1).
Rickettsia parkeri infection differentially regulates tick selenogene expression
We assessed the differential expression of ten tick selenoprotein genes in Rp+ and Rp− adult female ticks at 5 days post-infestation (dpi) of the ticks by quantitative reverse transcriptase (qRT)-PCR (Fig. 4). Six genes (eEFSec, SELENOK, SELENOS, SELENOO, SELENON, and SELENOX) were significantly upregulated (p < 0.05) in the midgut tissues, and five selenogenes (eEFSec, SELENOM, SELENOK, SELENOS, TrxR) were significantly upregulated (p < 0.05) in the salivary glands, while the levels of SELENOO, SELENON, and SELENOX remained unchanged (Fig. 4). Interestingly, only SELENOM and SELENOO genes were significantly upregulated in the Rp+-infected ovarian tissues, but SELENOT was not differentially regulated in any of the tested tissues (Fig. 4). Based on these results, the following three differentially expressed selenogenes encoding the following proteins were selected to determine their functional significance in tick blood feeding and pathogen colonization within the tick vector: endoplasmic reticulum-resident Selenoprotein M (gene, SELENOM), Selenoprotein S (gene, SELENOS), and mitochondrial resident Selenoprotein O (gene, SELENOO).
Temporal- and tissue-specific SELENOO and SELENOS transcript levels were assessed in the unfed and partially-fed midguts and salivary glands in from the Rp− female adults (Additional file 2: Figure S2). Interestingly, the gene expression level of SELENOO in the midgut tissues was upregulated 2–2.5-fold, while the salivary glands showed decreased transcriptional expression from 2 to 8 dpi (Additional file 2: Figure S2a). In contrast, the SELENOS transcript level was fourfold upregulated during the early phase of tick feeding compared with the tissue levels in the unfed ticks (Additional file 2: Figure S2). SELENOS expression remained unchanged in the midgut tissues from 2 to 8 dpi (Additional file 2: Figure S2b). The transcriptional activity of SELENOM has been reported to gradually increase up to 2 dpi but diminish thereafter in both tissue types .
Impact of selenoprotein gene silencing on tick physiology
We employed the RNAi approach to silence SELENOM, SELENOO, and SELENOS tick selenoproteins in Rp− and Rp+ A. maculatum ticks. The results obtained for the silencing of selenoprotein genes in Rp− ticks has been published for SELENOM  and is shown here as an Additional file 3: Figure S3. To examine the roles of each of the three selenogenes in R. parkeri-infected ticks, transcripts from the three selenogenes were individually depleted using RNAi (Fig. 5), and the effects of the depletion on R. parkeri replication were determined (Fig. 6). Silencing of these selenogenes did not negatively impact blood feeding in the ticks (data not shown). We achieved > 90% depletion of the transcriptional expression of all selenogenes under investigation in the uninfected ticks (except for SELENOO in Rp− ticks) (Additional file 3: Figure S3). Despite the effectiveness of the SELENOM knockdowns, the transcriptional expression levels of TrxR, Mn-SOD, Duox, Cat, Salp25D, SELENOK, SELENOO, and SELENOS were 2–3-fold upregulated in the partially blood-fed salivary glands of the Rp+ ticks (Fig. 5a). Moreover, TrxR, Mn-SOD, and Salp25D showed 2–3-fold increases in their expression levels in the midgut tissues (Fig. 5a). SELENOO depletion resulted in a 2–5-fold upregulation of Catalase, Cu/Zn-SOD, and Duox in the midgut tissues, whereas the other selenogenes were down-regulated in the salivary gland tissues (Fig. 5b). The SELENOS knockdown produced a 2–16-fold upregulation of Cu/Zn-SOD and Duox in the midgut tissues, while the transcript levels of the other tested genes remained unchanged or were downregulated in the salivary glands tissues (Fig. 5c). Tick SELENOM, SELENOO, and SELENOS silencing had no significant impact on the total oxidative stress levels, as estimated by the MDA lipid peroxidation method (Fig. 5d). The impact of knocking down these selenoproteins on ER stress was estimated by measuring the transcriptional expression of the unfolded protein response sensor genes, IRE1 and AFT6 (Table 1). The SELENOO and SELENOS knockdowns resulted in upregulation of both ATF6 and IRE1 sensor genes in the Rp+ tick midgut but not in the Rp− ticks (Additional file 4: Figure S4a), nonetheless opposite in SELENOM (Additional file 4: Figure S4b).
Impact of selenoprotein silencing on total bacterial load and R. parkeri replication
Knocking-down SELENOM resulted in a significant decrease in the R. parkeri concentration and the total bacterial load in the salivary gland tissues, but not in the midgut tissues after 5 dpi in the female ticks (Fig. 6a, b). The total bacterial load decreased in the midguts upon SELENOO knockdown, but the result was not statistically significant (p = 0.07), and the bacterial load in the salivary glands remained unaffected (Fig. 6c). The SELENOO knockdown depleted R. parkeri (p = 0.0062) in the midgut tissues, but the R. parkeri levels remained unchanged in the salivary glands (p = 0.06) (Fig. 6d). Finally, the SELENOS gene silencing depleted the total bacterial load in the midgut but not in the salivary gland tissues (Fig. 6e). Similarly, the SELENOS knockdown did not alter the R. parkeri load in the midgut tissues (p = 0.97), unlike in the salivary gland tissues where it was reduced, but not significantly so (p = 0.06) (Fig. 6f).
Quantification of tick symbionts upon selenogene silencing
Quantifying how the symbiont load changes in ticks may provide insight into potential pathogen–symbiont interactions inside ticks where selenogenes have been silenced. Therefore, FLE and CMM were quantified after selenogene knockdown in the tick tissues to assess how their levels might be altered by co-infection with R. parkeri (Fig. 7). In the SELENOM knockdown, the FLE levels were not significantly altered in the gut tissues (p = 0.35), but they were significantly reduced in the salivary glands (p = 0.0113). CMM loads in both midgut and salivary gland tissues were not affected by the SELENOM knockdown (p = 0.59 and p = 0.118, respectively). In contrast, the CMM levels were significantly depleted in the SELENOO knocked-down midgut tissues (p = 0.0042) (Fig. 7d) and in the salivary glands (p = 0.025) with the SELENOS knockdown (Fig. 7f).
Pathogen proliferation within tick tissues is an important aspect of the overall vector competence of ticks. This study investigated possible interactions between pathogenic R. parkeri and two non-pathogenic symbionts found in A. maculatum and whether selenoproteins might influence replication in R. parkeri and in the two tick symbionts, FLE and CMM (Fig. 8 and Table 1). In ticks, blood meal uptake can adversely affect bacterial replication because of the oxidative stress that is associated with it ; however, R. parkeri is able to defy oxidative stress and can multiply in the tick despite blood meal acquisition by the immature tick stages (Fig. 1a). The total bacterial load in Rp+ ticks compared with that in Rp− ticks is possibly higher because of the additional rickettsial and symbiont load in the Rp+ ticks (Fig. 2) which, like with Borrelia burgdorferi in I. scapularis, is able to increase in the tick tissues even in the presence of the blood meal . Both of the aforementioned studies lack detail about how many bacteria a tick inoculates into the host during feeding. In fact, they only considered the bacterial load residing in the tick tissues, particularly in the unfed and blood-fed stages (Figs. 1 and 2). Nevertheless, our results confirm that transovarial and transstadial maintenance of R. parkeri, FLE, and CMM occurs in A. maculatum, a finding reported previously for other tick species [42,43,44]. The genomic sequence of FLE from A. maculatum has been published , but no direct experiments have confirmed the specific role of this agent in the tick. But interestingly, an elegant study by Duron et al.  identified a Francisella type called F-Om, which they discovered was a “nutritional mutualist” involved in the synthesis of B vitamins in the in the African soft tick, Ornithodoros moubata. The role played by CMM in its host has not elucidated, but its presence may be related to the supply of important nutrients that are deficient in tick blood meals . Microbial symbionts play crucial roles in arthropod physiology and pathogen colonization, and innovative use of these symbionts could offer a novel method for controlling vector-borne disease transmission . The CMM microbes identified in one study have a unique intra-mitochondrial life cycle  and have been reported to occur in other blood feeding arthropods . Interestingly, large CMM loads have been found in I. holocyclus , and these microbes are known to flourish in the presence of a blood meal in I. ricinus , but this was not the case here with A. maculatum (Additional file 1: Figure S1). Rather, CMM in A. maculatum can be viewed as a successful colonizing partner of R. parkeri. Except in eggs and unfed larva, CMM loads were consistently higher in Rp+ ticks than in Rp− tick tissues including the midguts, salivary glands, and ovaries from fed or unfed female adult ticks. Unlike CMM, FLE levels became reduced in R. parkeri-infected female tick midguts but not in the salivary glands, ovaries, eggs, nymphs, or larvae (Figs. 2 and 3). These results suggest a possible synergistic relationship between R. parkeri and CMM during their transstadial transmission. Cafiso et al.  detected and quantified CMM bacteria in seven tick species and hypothesized diverse roles of this bacteria in variety of different tick species. Clearly, more work is needed in this area, but our results provide the first evidence that FLE could be an important symbiont in ticks, which like Wolbachia or Chromobacterium (Csp_P) in mosquitoes are able to modulate the replication of the malaria parasite and dengue virus [51, 52]. The refractoriness of pathogen transmission to humans from the arthropod vector has been linked to the mosquito symbiont, Wolbachia . The presence of FLE possibly interferes with R. parkeri replication in A. maculatum, which bears similarity to defensive symbiosis where the presence of one or more symbionts interferes with pathogenic bacteria colonization of the vector gut  but differs from nutritional symbiosis where the bacterial symbionts supplement the nutrients in the host [54, 55]. Here, based on the higher concentration of CMM that we found in the infected tick tissues and the reduced R. parkeri load, we believe that the presence of CMM favors the replication of R. parkeri replication (Figs. 2, 3, and 6). We also believe that CMM might have a synergistic effect on R. parkeri, and displacing this symbiont may also displace R. parkeri-infected ticks in the field, as was recently shown in a study on Candidatus Rickettsia andeanae .
We further investigated the dynamic interplay occurring among R. parkeri, CMM, and FLE microbes in ticks depleted of R. parkeri via selenoprotein silencing (Figs. 6 and 7). We found that infection with R. parkeri differentially regulated a battery of tick selenoproteins in the midgut, salivary glands, and ovarian tissues (Fig. 4). We showed the differential expression of selenogenes within and among the tick organs with respect to the pathogen infection. The selenoproteins in the salivary glands were highly expressed during infection with R. parkeri compared with those in the midguts and ovarian tissues (Fig. 4). Pathogen development and the secretory functions of the tick salivary glands might have resulted in the higher levels of selenoprotein expression in tick salivary glands that were observed. Furthermore, during tick feeding, the tick salivary glands probably remained under stress because these glands participate in the constant supply of anti-hemostatic, anti-inflammatory, and analgesic compounds during the continuous flow of the blood meal [57, 58]. ROS generation is one of the first lines of defense against invading microbes [59, 60]. However, despite minimal investigation, evidence is now accumulating that the tick selenoproteome and antioxidant enzymes may play critical roles in detoxifying ROS and in maintaining both vector microbiota and R. parkeri colonization [10, 31, 34, 36,37,38].
The compensatory actions of the redox genes following selenogene transcript depletion via RNAi differed between the Rp− and Rp+ ticks (Fig. 5, Additional file 3: Figure S3). SELENOM depletion in the Rp+ ticks showed evidence of compensation by overexpression of TrxR, Mn-SOD, and Salp25D in the midgut and ovary tissues, whereas CAT, SELENOK, and SELENOS were only upregulated in the salivary glands (Fig. 5a). Cu/Zn-SOD and Duox were significantly upregulated after SELENOO and SELENOS silencing; these are involved in known defense mechanisms against invading pathogens and are also involved in repairing the tissue damage from Rickettsia-dependent superoxide generation . Superoxide generation is associated with rickettsial infections , and tick extracellular Cu/Zn-SOD is the main quencher for dismutation of the superoxides generated during rickettsial infections . Upregulated Cu/Zn-SOD probably provides the redox balance required to offset the superoxide radicals generated by the rickettsial infections in the ticks after SELENOO and SELENOS were knocked down, but Cu/Zn-SOD was not upregulated in the SELENOM knockdowns (Fig. 5). Further investigation of the unfolded protein response sensor genes (ATF6 and IRE1) provided evidence of the altered protein folding homeostasis inside the ER, the organelle necessary for the proper folding of all secretory and transmembrane proteins (Additional file 4: Figure S4). The sensor genes for the unfolded protein response (ATF6 and IRE1) and tick selenoprotein silencing potentially cause ER stress . Knockdown of the mitochondrial resident selenogene, SELENOO, likely induced high oxidative stress in the gut tissues (Fig. 5, Additional file 4: Figure S4), which in turn induced ER stress. Studies have shown that there is a physical and biochemical interaction between the ER and mitochondria , and mitochondrial ROS can induce ER stress . We proposed a model to summarize the important points arising from our study (Fig. 8). In this model, we suggest that successful R. parkeri replication within the tick vector is enhanced by the presence of CMM, probably by displacing FLE. The selenogenes responding to R. parkeri infection by transcriptional upregulation favor R. parkeri replication, and this in turn enhances the overall vectorial competence of A. maculatum for R. parkeri.
The successful growth of a human pathogenic spotted fever group rickettsia, R. parkeri, inside its competent vector, A. maculatum, offers it a chance to dynamically interact with tick symbionts and modulate its host’s defenses by upregulating tick selenoproteins. This study illustrates the potential of a new research model aimed at providing better understanding of tick–pathogen interactions and the important roles played by symbionts and various tick factors in regulating pathogen growth.
All the animal experiments were performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health, USA. The protocols for tick blood feeding were approved by the Institutional Animal Care and Use Committee of the University of Southern Mississippi (protocol#15101501 and 15011402).
Ticks and tissue preparations
A. maculatum ticks were maintained at the University of Southern Mississippi according to established methods . The Rp+ tick colonies and Rp− tick colonies from Mississippi field collections were established and maintained in the laboratory. A. maculatum colonies containing individual Rp+ and Rp− ticks were established in our laboratory in 2013. Questing unfed adult ticks were collected from Mississippi Sandhill Crane, National Wildlife Refuge, Gautier, Mississippi (www.fws.gov/refuge/mississippi_sandhill_crane/) using the drag cloth method during July and August 2013. The hundreds of ticks collected from the field were blood-fed on sheep and allowed to fully engorge and drop off. Fully engorged female adult ticks were kept in snap vials for egg laying. Typically, a single fully engorged female A. maculatum lays on average a 450–950 mg egg mass containing 15,000–18,000 eggs. To determine the presence of R. parkeri in an egg mass, 10–20 mg of it from each gravid female was subjected to genomic DNA extraction using the DNeasy extraction kit as described by the manufacturers (Qiagen, CA) followed by PCR amplification of the extracted DNA to identify SFGR species . Individual Rp+ and Rp− eggs were selected from individual gravid females and allowed to hatch into unfed larva. The unfed larval ticks were blood-fed by allowing them to infest golden Syrian hamsters until they reached repletion. Fully engorged larvae were allowed to molt into nymphs and then blood-fed on hamsters. Fully engorged nymphs molted as male or female ticks. Each developmental stage was routinely tested for the presence of R. parkeri infection. Closed colonies from the third and fourth generation of the original wild-caught ticks were used in this study after confirming the presence or absence of infection in the colonies. Freshly laid eggs, freshly molted larvae and blood-fed larvae, freshly hatched nymphs, and blood-fed nymphs were collected from Rp+ and Rp− colonies separately. Eggs, unfed larvae (from three individual ticks, 20 mg each), three fed and pooled larval batches, unfed nymphs (20 mg), and fed nymphs were stored in RNAlater (Invitrogen, Carlsbad, CA). At least three biological replicates were used in all the experiments. Tick tissues from the unfed and partially blood-fed female adult ticks were dissected and stored immediately in RNAlater (Invitrogen) prior to extracting the mRNA using our previously described laboratory method .
RNA preparation, cDNA synthesis, and qRT-PCR
Total RNA extraction, cDNA synthesis, and qRT-PCR were conducted as previously described . The gene-specific primer sequences, which were designed to amplify specific cDNA fragments from A. maculatum tissues, are listed in Additional file 5: Table S1. Transcriptional gene expression of the tick genes in Rp− ticks was normalized against β-actin gene expression, while GAPDH gene expression was used to normalize tick gene expression in the Rp+ tick tissues because it is stably expressed irrespective of the infection status . The synthesized cDNA was used to measure mRNA levels by qRT-PCR using the CFX96 Real Time System (Bio-Rad Inc., Hercules, CA) as described previously [10, 38].
Double-stranded RNA (dsRNA) synthesis, tick injections, and hematophagy
dsRNA was synthesized to allow for the in vivo analysis of SELENOM, SELENOO, and SELENOS in the ticks. Tick manipulations were performed according to the methods described previously [36, 68]. The dsRNAs for each selenoprotein gene (dsSELENOM, dsSELENOO, dsSELENOS) were diluted to working concentrations of 1 μg/μL in nuclease-free water. The same protocol was used to synthesize dsLacZ, which was used as an irrelevant dsRNA control. Twenty-five unfed adult female ticks were each microinjected with 1 μl of dsRNA or dsLacZ using a 27-gauge needle, then kept overnight at 37 °C to alleviate needle trauma and promote their survival, after which they were blood-fed using routine laboratory procedures .
Quantification of total bacterial loads
The bacterial load in each tick tissue was estimated as described previously [10, 15]. The bacterial copy numbers were normalized against A. maculatum actin expression in the uninfected ticks and GAPDH expression in the Rp+ ticks.
Quantification of FLE in tick tissues
The FLE from A. maculatum tick tissues was quantified using the primers described elsewhere . The serially diluted copies (108 to 101) of each gene were PCR-amplified using predetermined thermal cycling conditions, and the Ct values for known dilutions were used to construct a standard curve from which the copy number of each gene was calculated. The 25 μL qRT-PCRs comprised 125 nM of each primer, SYBR Green Master Mix (Bio-Rad, Inc. USA), and the serially diluted PCR products prepared for each standard curve. The reaction mixtures were subjected to the thermal cycle parameters of 95 °C for 5 min followed by 29 cycles of 95 °C for 30s, 52 °C for 30s, and 72 °C for 30s with a final extension of 72 °C for 5 min in a CFX96 Real Time System (Bio-Rad Inc.). The FLE copy numbers were normalized against the A. maculatum GAPDH gene. As with the other qRT-PCRs, all the samples were run in triplicate.
Quantification of CMM in tick tissues
We followed the protocol described previously for quantifying CMM . The CMM-specific GyrB gene and tick GAPDH gene were PCR-amplified from A. maculatum ticks using the primers shown in Table 1. The amplified PCR products serially diluted tenfold (108 to 101 copies) were used to generate a standard curve. The qRT-PCRs comprised 400 nM of each primer and 25 ng of the cDNA samples. The reaction mixture containing SYBR Green (Bio-Rad Inc. USA) was subjected to thermal cycling at 95 °C for 2 min, 40 cycles at 95 °C for 15 s, and at 60 °C for 30s, and a melting curve from 55 °C to 95 °C with increasing increments of 0.5 °C per cycle was prepared using the CFX96 Real Time System (Bio-Rad Inc.). The standard curves generated were used to calculate the copy numbers of the CMM GyrB gene and the tick GAPDH gene. The CMM copy numbers were normalized against the A. maculatum GAPDH gene. As with the other qRT-PCRs, all the samples were run in triplicate.
Quantification of the R. parkeri load in tick tissues
The level of infection with R. parkeri within the tick tissues across the developmental stages was quantified using a slightly modified version of a previously published method [9, 70]. The R. parkeri load was estimated as the ratio of R. parkeri-specific rompB gene copies to tick GAPDH copies. GAPDH and rompB genes were amplified using 250 nM of each specific primer (Table 1) in a reaction containing SYBR Green Master Mix (Bio-Rad Inc.) in the CFX96 Real Time System (Bio-Rad Inc.) with thermal cycling conditions of 95 °C for 10 min, followed by 35 cycles of 95 °C for 15 s, 60 °C for 30s, and 72 °C for 30s. The standard curves for tick GAPDH and rompB were prepared based on the amplification profiles of known concentrations of purified GAPDH and rompB PCR products. The standard curves generated were used to estimate the copy numbers of each gene in the tick samples.
Quantification of total oxidative stress levels
The malondialdehyde lipid peroxidation assay kit (Sigma-Aldrich, St. Louis, MO, USA) was used to quantify lipid degradation as a result of oxidative damage . All the procedures followed the manufacturer’s recommendations, and all the samples were balanced by weight.
All data are expressed as mean values ± SEM unless otherwise indicated. Statistical significance between two experimental groups or their respective controls was determined by a t test (p value, 0.05). Comparative differences among multiple experimental groups were determined by analysis of variance with statically significant p values of < 0.05 (GraphPad Prism 6.05, La Jolla, CA). Transcriptional expression levels were determined using Bio-Rad software (Bio-Rad CFX MANAGER v.3.1), and the gene expression values obtained were considered statistically significant if a p value of 0.05 was obtained when compared with the control.
Raoult D, Roux V. Rickettsioses as paradigms of new or emerging infectious diseases. Clin Microbiol Rev. 1997;10:694–719.
Paddock CD, Goddard J. The evolving medical and veterinary importance of the Gulf Coast tick (Acari: Ixodidae). J Med Entomol. 2015;52:230–52.
Walker DH, Ismail N. Emerging and re-emerging rickettsioses: endothelial cell infection and early disease events. Nat Rev Microbiol. 2008;6:375–86.
Parola P, Paddock CD, Raoult D. Tick-borne rickettsioses around the world: emerging diseases challenging old concepts. Clin Microbiol Rev. 2005;18:719–56.
Rikihisa Y. Anaplasma phagocytophilum and Ehrlichia chaffeensis: subversive manipulators of host cells. Nat Rev Microbiol. 2010;8:328–39.
Burgdorfer W, Brinton LP, Hughes LE. Isolation and characterization of symbiotes from the Rocky Mountain wood tick, Dermacentor andersoni. J Invertebr Pathol. 1973;22:424–34.
Telford SR. Status of the “East Side hypothesis” (transovarial interference) 25 years later. Ann N Y Acad Sci. 2009;1166:144–50.
Klyachko O, Stein BD, Grindle N, Clay K, Fuqua C. Localization and visualization of a coxiella-type symbiont within the lone star tick, Amblyomma americanum. Appl Environ Microbiol. 2007;73:6584–94. https://doi.org/10.1128/AEM.00537-07.
Budachetri K, Browning RE, Adamson SW, Dowd SE, Chao C-C, Ching W-M, et al. An insight into the microbiome of the Amblyomma maculatum (Acari: Ixodidae). J Med Entomol. 2014;51:119–29.
Budachetri K, Karim S. An insight into the functional role of thioredoxin reductase, a selenoprotein, in maintaining normal native microbiota in the Gulf Coast tick (Amblyomma maculatum). Insect Mol Biol. 2015;24:570–81.
Budachetri K, Gaillard D, Williams J, Mukherjee N, Karim S. A snapshot of the microbiome of Amblyomma tuberculatum ticks infesting the gopher tortoise, an endangered species. Tick Tick Borne Dis. 2016;7:1225–9.
Budachetri K, Williams J, Mukherjee N, Sellers M, Moore F, Karim S. The microbiome of neotropical ticks parasitizing on passerine migratory birds. Ticks Tick Borne Dis. 2017;8(1):170–3.
Karim S, Budachetri K, Mukherjee N, Williams J, Kausar A, Hassan JM, Adamson SW, et al. A study of ticks and tick-borne livestock pathogens in Pakistan. PLoS Negl Trop Dis. 2017;11(6). https://doi.org/10.1371/journal.pntd.0005681.
Narasimhan S, Fikrig E. Tick microbiome: the force within. Trends Parasitol. 2015;31:315–23.
Narasimhan S, Rajeevan N, Liu L, Zhao YO, Heisig J, Pan J, et al. Gut microbiota of the tick vector Ixodes scapularis modulate colonization of the Lyme disease spirochete. Cell Host Microbe. 2014;15:58–71.
Cirimotich CM, Ramirez JL, Dimopoulos G. Native microbiota shape insect vector competence for human pathogens. Cell Host Microbe. 2011;10:307–10.
Chmelař J, Kotál J, Kopecký J, Pedra JHF, Kotsyfakis M. All for one and one for all on the tick-host battlefield. Trends Parasitol. 2016; https://doi.org/10.1016/j.pt.2016.01.004.
Niebylski ML, Peacock MG, Fischer ER, Porcella SF, Schwan TG. Characterization of an endosymbiont infecting wood ticks, Dermacentor andersoni, as a member of the genus Francisella. Appl Environ Microbiol. 1997a;63:3933–40.
Niebylski ML, Schrumpf ME, Burgdorfer W, Fischer ER, Gage KL, Schwan TG. Rickettsia peacockii sp. nov., a new species infecting wood ticks, Dermacentor andersoni, in Western Montana. Int J Syst Bacteriol. 1997b;47:446–52.
Ahantarig A, Trinachartvanit W, Baimai V, Grubhoffer L. Hard ticks and their bacterial endosymbionts (or would be pathogens). Folia Microbiol (Praha). 2013;58:419–28.
Sun LV, Scoles GA, Fish D, O’Neill SL. Francisella-like endosymbionts of ticks. J Invertebr Pathol. 2000;76:301–3.
Scoles GA. Phylogenetic analysis of the Francisella-like endosymbionts of Dermacentor ticks. J Med Entomol. 2004;41:277–86.
Vanov IN, Mitkova N, Reye AL, Hubschen JM, Vatcheva-Dobrevska RS, Dobreva EG, Kantardjiev TV, Muller CP. Detection of new Francisella-like tick endosymbionts in Hyalomma spp. and Rhipicephalus spp. (Acari: Ixodidae) from Bulgaria. Appl Env Microbiol. 2011:77:5562–65. https://doi.org/10.1128/AEM.02934-10AEM.02934-10.
Foley JE, Nieto NC. Tularemia. Vet Microbiol. 2010;140:332–8.
Dergousoff SJ, Chilton NB. Association of different genetic types of Francisella-like organisms with the Rocky Mountain wood tick (Dermacentor andersoni) and the American dog tick (Dermacentor variabilis) in localities near their northern distributional limits. Appl Environ Microbiol. 2012;78:965–71.
Liu J-N, Yu Z-J, Liu L-M, Li N-X, Wang R-R, Zhang C-M, et al. Identification, distribution and population dynamics of Francisella-like endosymbiont in Haemaphysalis doenitzi (Acari: Ixodidae). Sci Rep. 2016;6:35178.
Gerhart JG, Moses AS, Raghavan R. A Francisella-like endosymbiont in the Gulf Coast tick evolved from a mammalian pathogen. Sci Rep. 2016;6:33670.
Epis S, Mandrioli M, Genchi M, Montagna M, Sacchi L, Pistone D, Sassera D. Localization of the bacterial symbiont Candidatus Midichloria mitochondrii within the hard tick Ixodes ricinus by whole-mount FISH staining. Ticks Tick Borne Dis. 2013:4:39–45. https://doi.org/10.1016/j.ttbdis.2012.06.005.
Montagna M, Sassera D, Epis S, Bazzocchi C, Vannini C, Lo N, et al. ‘Candidatus Midichloriaceae’ fam. nov. (Rickettsiales), an ecologically widespread clade of intracellular alphaproteobacteria. Appl Environ Microbiol. 2013;79:3241–8.
Epis S, Sassera D, Beninati T, Lo N, Beati L, Piesman J, et al. Midichloria mitochondrii is widespread in hard ticks (Ixodidae) and resides in the mitochondria of phylogenetically diverse species. Parasitology. 2008;135:485–94.
Budachetri K, Crispell G, Karim S. Amblyomma maculatum SECIS binding protein 2 and putative selenoprotein P are indispensable for pathogen replication and tick fecundity. Insect Biochem Mol Biol. 2017. https://doi.org/10.1016/j.ibmb.2017.07.006.
Karim S, Singh P, Ribeiro JMC. A deep insight into the sialotranscriptome of the Gulf Coast tick, Amblyomma maculatum. PLoS One. 2011; https://doi.org/10.1371/journal.pone.0028525.
Karim S, Ribeiro JM. An insight into the sialome of the lone star tick, Amblyomma americanum, with a glimpse on its time dependent gene expression. PLoS One. 2015. https://doi.org/10.1371/journal.pone.0131292.
Adamson SW, Browning RE, Budachetri K, Ribeiro JMC, Karim S. Knockdown of selenocysteine-specific elongation factor in Amblyomma maculatum alters the pathogen burden of Rickettsia parkeri with epigenetic control by the Sin3 histone deacetylase corepressor complex. PLoS One. 2013. https://doi.org/10.1371/journal.pone.0082012.
Adamson S, Browning R, Singh P, Nobles S, Villarreal A, Karim S. Transcriptional activation of antioxidants may compensate for selenoprotein deficiencies in Amblyomma maculatum (Acari: Ixodidae) injected with selK- or selM-dsRNA. Insect Mol Biol. 2014;23:497–510.
Kumar D, Budachetri K, Meyers VC, Karim S. Assessment of tick antioxidant responses to exogenous oxidative stressors and insight into the role of catalase in the reproductive fitness of the Gulf Coast tick, Amblyomma maculatum. Insect Mol Biol. 2016;25:283–94.
Budachetri K, Kumar D, Karim S. Catalase is a determinant of the colonization and transovarial transmission of Rickettsia parkeri in the Gulf Coast tick Amblyomma maculatum. Insect Mol Biol. 2017;26(4):414–8.
Crispell G, Budachetri K, Karim S. Rickettsia parkeri colonization in Amblyomma maculatum: the role of superoxide dismutases. Parasit Vectors. 2016;9:291.
Narasimhan S, Sukumaran B, Bozdogan U, Thomas V, Lian X, DePonte K, Marcantonio N, Anderson FJ, Kantor F, Fikrig E. A tick antioxidant facilitates the Lyme disease agent’s successful migration from mammalian host to the arthropod. Cell Host Microbe. 2007;2:7–18.
Kocan KM, Zivkovic Z, Blouin EF, Naranjo V, Almazán C, Mitra R, et al. Silencing of genes involved in Anaplasma marginale-tick interactions affects the pathogen developmental cycle in Dermacentor variabilis. BMC Dev Biol. 2009;9:42.
Piesman J, Schneider BS, Zeidner NS. Use of quantitative PCR to measure density of Borrelia burgdorferi in the midgut and salivary glands of feeding tick vectors. J Clin Microbiol. 2001;39:4145–8.
Baldridge GD, Scoles GA, Burkhardt NY, Schloeder B, Kurtti TJ, Munderloh UG. Transovarial transmission of Francisella-like endosymbionts and Anaplasma phagocytophilum variants in Dermacentor albipictus (Acari: Ixodidae). J Med Entomol. 2009;46:625–32.
Sassera D, Beninati T, Bandi C, Bouman EAP, Sacchi L, Fabbi M, et al. ‘Candidatus Midichloria mitochondrii’, an endosymbiont of the tick Ixodes ricinus with a unique intramitochondrial lifestyle. Int J Syst Evol Microbiol. 2006;56:2535–40.
Wright CL, Gaff HD, Sonenshine DE, Hynes WL. Experimental vertical transmission of Rickettsia parkeri in the Gulf Coast tick, Amblyomma maculatum. Ticks Tick Borne Dis. 2015;6:568–73.
Duron O, Morel O, Noel V, Buysse M, Binetury F, Lancelot R, et al. Tick-bacteria mutualism depends on B vitamin synthesis pathways. Curr Biol. 2018; https://doi.org/10.1016/j.cub.2018.04.038.
Rio RVM, Attardo GM, Weiss BL. Grandeur alliances: symbiont metabolic integration and obligate arthropod hematophagy. Trends Parasitol. 2016;32:739–49.
Richard S, Seng P, Parola P, Raoult D, Davoust B, Brouqui P. Detection of a new bacterium related to ‘Candidatus Midichloria mitochondrii’ in bed bugs. Clin Microbiol Infect. 2009;15:84–5.
Gofton AW, Oskam CL, Lo N, Beninati T, Wei H, McCarl V, et al. Inhibition of the endosymbiont ‘Candidatus Midichloria mitochondrii’ during 16S rRNA gene profiling reveals potential pathogens in Ixodes ticks from Australia. Parasit Vectors. 2015;8:345.
Sassera D, Lo N, Bouman EAP, Epis S, Mortarino M, Bandi C. ‘Candidatus Midichloria’ endosymbionts bloom after the blood meal of the host, the hard tick Ixodes ricinus. Appl Environ Microbiol. 2008;74:6138–40.
Cafiso A, Bazzocchi C, Marco DL, Opara MN, Sassera D, Plantard O. Molecular screening for Midichloria in hard and soft ticks reveals variable prevalence levels and bacterial loads in different tick species. Ticks Tick Borne Dis. 2016;7:1186–92.
Hughes GL, Rivero A, Rasgon JL. Wolbachia can enhance Plasmodium infection in mosquitoes: implications for malaria control? PLoS Pathog. 2014; https://doi.org/10.1371/journal.ppat.1004182.
Ramirez JL, Short SM, Bahia AC, Saraiva RG, Dong Y, Kang S, et al. Chromobacterium Csp_P reduces malaria and dengue infection in vector mosquitoes and has entomopathogenic and in vitro anti-pathogen activities. PLoS Pathog. 2014; https://doi.org/10.1371/journal.ppat.1004398.
Clay K. Defensive symbiosis: a microbial perspective. Funct Ecol. 2014;28:293–8.
Baumann P. Biology of bacteriocyte-associated endosymbionts of plant sap-sucking insects. Annu Rev Microbiol. 2005;59:155–89.
Smith TA, Driscoll T, Gillespie JJ, Raghavan R. A Coxiella-like endosymbiont is a potential vitamin source for the lone star tick. Genome Biol Evol. 2015; https://doi.org/10.1093/gbe/evv016.
Paddock CD, Denison AM, Dryden MW, Noden BH, Lash RR, Abdelghani SS, et al. High prevalence of “Candidatus Rickettsia andeanae” and apparent exclusion of Rickettsia parkeri in adult Amblyomma maculatum (Acari: Ixodidae) from Kansas and Oklahoma. Ticks Tick Borne Dis. 2015;6:297–302.
Binnington KC. Sequential changes in salivary gland structure during attachment and feeding of the cattle tick, Boophilus microplus. Int J Parasitol. 1978;8:97–115.
Francischetti IMB, Sa-Nunes A, Mans BJ, Santos IM, Ribeiro JMC. The role of saliva in tick feeding. Front Biosci. 2009;14:2051–88.
Hoffman JS. The immune responses of Drosophila. Nature. 2003;426:33–8.
Ha EM, Oh CT, Bae YS, Lee WJ. A direct role for dual oxidase in Drosophila gut immunity. Science. 2005;310:847–50.
Hong JE, Santucci LA, Tian X, Silverman DJ. Superoxide dismutase-dependent, catalase-sensitive peroxides in human endothelial cells infected by Rickettsia rickettsii. Infect Immun. 1998;66:1293–8.
Santucci LA, Gutierrez PL, Silverman DJ. Rickettsia rickettsii induces superoxide radical and superoxide dismutase in human endothelial cells. Infect Immun. 1992;60:5113–8.
Gardner BM, Pincus D, Gotthardt K, Gallagher CM, Walter P. Endoplasmic reticulum stress sensing in the unfolded protein response. Cold Spring Harb Perspect Biol. 2013;5:a013169.
Csordás G, Várnai P, Golenár T, Roy S, Purkins G, Schneider TG, et al. Imaging interorganelle contacts and local calcium dynamics at the ER-mitochondrial interface. Mol Cell. 2010;39:121–32.
Pierre N, Barbé C, Gilson H, Deldicque L, Raymackers J-M, Francaux M. Activation of ER stress by hydrogen peroxide in C2C12 myotubes. 2014; https://doi.org/10.1016/j.bbrc.2014.05.143.
Patrick CD, Hair JA. Laboratory rearing procedures and equipment for multi-host ticks (Acarina: Ixodidae). J Med Entomol. 1975;12:389–90.
Browning R, Adamson SW, Karim S. Choice of a stable set of reference genes for qRT-PCR analysis in Amblyomma maculatum (Acari: Ixodidae). J Med Entomol. 2012;49:1339–46.
Bullard RL, Williams J, Karim S. Temporal gene expression analysis and RNA silencing of single and multiple members of gene family in the lone star tick Amblyomma americanum. PLoS One. 2016; https://doi.org/10.1371/journal.pone.0147966.
Karim S, Adamson SW, Simpson SJ, Casas J. RNA interference in ticks: a functional genomics tool for the study of physiology. In: Jockusch E, editor. Advances in insect physiology small RNAs: their diversity, roles and practical uses, vol. 42. San Diego: Academic Press; 2012. p. 119–54.
Jiang J, Stromdahl EY, Richards AL. Detection of Rickettsia parkeri and Candidatus Rickettsia andeanae in Amblyomma maculatum Gulf Coast ticks collected from humans in the United States. Vector Borne Zoonotic Dis. 2012;12:175–82.
This work was supported by grants from the National Institute of Allergy and Infectious Diseases (award #AI099919), United States Department of Agriculture, National Institute of Food and Agriculture (awards #2017-67016-26864), and the National Institutes of General Medical Sciences (award # P20RR016476). The funders played no role in the study design, data collection, analysis, decision to publish, or manuscript preparation.
Ethics approval and consent to participate
Use of animals for tick blood-feeding was approved by the IACUC of the University of Southern Mississippi.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. FLE and CMM loads across the blood meal in naïve tick tissues. FLE loads in midguts (a) and salivary glands (b), and CMM loads in midguts (c) and salivary glands (d) at different time points in the tick tissues. (DOCX 117 kb)
Figure S2. Time-dependent SELENOO (a) and SELENOS (b) transcriptional expression levels in uninfected (naïve) tick midguts and salivary gland tissues during the adult female blood meal. The change in transcriptional activity of (a) SELENOO and (b) SELENOS in A. maculatum midgut and salivary gland tissues was normalized to that of the unfed tick using β-actin as a reference gene. (DOCX 135 kb)
Figure S3. Knockdown of tick selenoproteins in naïve (uninfected) A. maculatum. (a) dsSELENOO and (b) dsSELENOS. Gene expression in naïve ticks was normalized against tick β-actin as the reference gene. Compensatory antioxidant expression levels were measured for eukaryotic elongation factor (eEFSec), selenoproteins (SELENOM, SELENOK, SELENOS, SELENOO, TrxR, SELENON, SELENOT), mitochondrial superoxide dismutase (Mn-SOD), cytosolic superoxide dismutase (Cu/Zn-SOD), catalase (Cat), glutathione reductase (GSHR), and glutathione peroxidase (Salp25D). (DOCX 231 kb)
Figure S4. The unfolded protein response. The unfolded protein response estimation for the salivary glands (SG) and midgut (MG) based on transcriptional gene expression of the ATF6 and IRE1 sensor genes when selenogenes SELENOO and SELENOS were silenced in Rp− ticks (a) and Rp+ ticks (b). The effects of the SELENOM silenced tissues were also measured in the Rp+ ticks (b). (DOCX 177 kb)
Table S1. Gene-specific PCR primers and probes used for qRT-PCR and RNAi studies. (DOCX 42 kb)