Jose Reck1*†
Anelise Webster1
Bruno Dall’Agnol1
Ronel Pienaar2
Minique H. de Castro3
Jonathan Featherston3
Ben J. Mans2,4,5*†- 1Instituto de Pesquisas Veterinárias Desidério Finamor (IPVDF), Eldorado do Sul, Brazil
- 2Epidemiology, Parasites and Vectors, Agricultural Research Council, Onderstepoort Veterinary Research, Pretoria, South Africa
- 3Agricultural Research Council, Biotechnology Platform, Pretoria, South Africa
- 4Department of Veterinary Tropical Diseases, Vector and Vector-borne Disease Research Programme, University of Pretoria, Pretoria, South Africa
- 5Department of Life and Consumer Sciences, University of South Africa, Pretoria, South Africa
Tick salivary glands produce and secrete a variety of compounds that modulate host responses and ensure a successful blood meal. Despite great progress made in the identification of ticks salivary compounds in recent years, there is still a paucity of information concerning salivary molecules of Neotropical argasid ticks. Among this group of ticks, considering the number of human cases of parasitism, including severe syndromes and hospitalization, Ornithodoros brasiliensis can be considered one of the major Neotropical argasid species with impact in public health. Here, we describe the transcriptome analysis of O. brasiliensis salivary glands (ObSG). The transcriptome yielded ~14,957 putative contigs. A total of 368 contigs were attributed to secreted proteins (SP), which represent approximately 2.5% of transcripts but ~53% expression coverage transcripts per million. Lipocalins are the major protein family among the most expressed SP, accounting for ~16% of the secretory transcripts and 51% of secretory protein abundance. The most expressed transcript is an ortholog of TSGP4 (tick salivary gland protein 4), a lipocalin first identified in Ornithodoros kalahariensis that functions as a leukotriene C4 scavenger. A total of 55 lipocalin transcripts were identified in ObSG. Other transcripts potentially involved in tick-host interaction included as: basic/acid tail secretory proteins (second most abundant expressed group), serine protease inhibitors (including Kunitz inhibitors), 5' nucleotidases (tick apyrases), phospholipase A2, 7 disulfide bond domain, cystatins, and tick antimicrobial peptides. Another abundant group of proteins in ObSG is metalloproteases. Analysis of these major protein groups suggests that several duplication events after speciation were responsible for the abundance of redundant compounds in tick salivary glands. A full mitochondrial genome could be assembled from the transcriptome data and confirmed the close genetic identity of the tick strain sampled in the current study, to a tick strain previously implicated in tick toxicoses. This study provides novel information on the molecular composition of ObSG, a Brazilian endemic tick associated with several human cases of parasitism. These results could be helpful in the understanding of clinical findings observed in bitten patients, and also, could provide more information on the evolution of Neotropical argasids.
Introduction
Ticks (Ixodida) are composed of three families, the hard (Ixodidae), soft (Argasidae), and monotypic Nuttalliellidae (Guglielmone et al., 2010). The Ixodida is obligate blood-feeding ecto-parasites that secrete a cocktail of bioactive salivary gland-derived components during feeding, to counteract the vertebrate host’s defense mechanisms, such as blood-clotting, platelet aggregation, and inflammation (Mans, 2019). Ixodids in all life stages feed for prolonged periods that may last for days to weeks during which hundreds to thousands of salivary gland proteins are differentially secreted, presumably to evade the host’s changing immune responses, but also to adapt to a changing feeding environment due to the host’s healing responses (Francischetti et al., 2009). Argasids show much more diverse feeding behavior in different life stages. In the larval stages, some species molt to nymphs without feeding, other species feed rapidly within minutes, while some species feed for prolonged periods of days-weeks resembling ixodids. In the nymphal stages, instars from some species do not feed, but molt to the next developmental stage, while females from some species also do not feed, although the majority of species requires a blood meal for successful oviposition (Hoogstraal, 1985).
In the case of Ornithodoros brasiliensis Aragão, 1923, the larvae molt to nymphs without feeding, while nymphal stages fed within 25–35 min similar to adults (Ramirez et al., 2016). Bites by this tick species result in toxicoses most probably due to injection of salivary gland components into the feeding site (Reck et al., 2011, 2013a, 2014; Dall’Agnol et al., 2019). Toxicoses by the “mouro tick” exhibit as erythemic and swollen lesions, hyperemia of oral/ocular mucosa pruritus, and tachypnea and slow wound-healing processes, and in recent years, cases of human hospitalization after tick bite have been reported (Reck et al., 2011, 2013a,b, 2014; Dall’Agnol et al., 2019). It has been shown that salivary gland homogenates can inhibit wound healing and endothelial cell proliferation in vitro (Reck et al., 2013b). To date, little is known about the salivary gland protein composition of this tick species that would help to explain these effects or symptoms.
Salivary gland transcriptome sequencing is a useful tool to generate catalogs of salivary gland-derived transcripts and has been important in the description of salivary gland protein sequence diversity (Mans, 2020). As such, characterization of the salivary gland transcriptome of O. brasiliensis would create an important resource to elucidate the molecular mechanisms behind mouro tick toxicoses. Argasid salivary gland transcriptomes have been characterized in detail using conventional cDNA library Sanger sequencing for Antricola delacruzi Estrada-Peña, Barros-Battesti, and Venzal, 2004 (Ribeiro et al., 2012), Argas monolakensis Schwan, Corwin, and Brown, 1992 (Mans et al., 2008a), Ornithodoros coriaceus Koch, 1844 (Francischetti et al., 2008a), and Ornithodoros parkeri Cooley, 1936 (Francischetti et al., 2008b). Argasid salivary gland transcriptomes have also been characterized using next-generation sequencing and assembly strategies for Ornithodoros moubata (Murray, 1877) (Pérez-Sánchez et al., 2021), Ornithodoros erraticus (Lucas, 1849) (Oleaga et al., 2021), and Ornithodoros rostratus Aragão, 1911 (Araujo et al., 2019). The salivary gland transcriptome of O. brasiliensis has not been described yet. The aim of the current study was to sequence the nymphal salivary gland transcriptome of O. brasiliensis ticks using next-generation sequencing.
Materials and Methods
Tick Collection and RNA Extraction
Nymphs of O. brasiliensis were collected in the field from a site previously implicated in tick parasitism of travelers in Brazil (Dall’Agnol et al., 2019) and maintained unfed in the laboratory for ~2 months before dissection. Salivary glands were dissected from 10 nymphs and placed in RNA later before storing at −70°C. Total RNA was extracted using the RNeasy Protect Mini Kit (QIAGEN Group). Briefly, glands were suspended in 500 μl RLT buffer and disrupted by 10X passage through an 18G needle followed by 10X passage using a 24G needle. Residual genomic DNA was removed with DNase I digestion. Total RNA quantification was performed using the Qubit fluorimeter 2.0 (Life Technologies, Carlsbad, CA).
Library Construction and Next-Generation Sequencing
For library preparation, 1.0 μg purified total RNA was used with the TruSeq stranded mRNA sample preparation kit (Illumina, San Diego, CA). Poly-A mRNA was isolated, fragmented (for 3 min), and converted to double-stranded cDNA, adapters ligated, and PCR amplified for 12 cycles. Amplified bands were size selected from 450 to 1,200 bp. Bands were excised, purified, and sequenced using the Illumina MiSeq system (300 bp×300 bp). Raw sequence reads were submitted to GenBank under BioProject PRJNA719007 with small read archive accession number: SRR14139641.
Transcriptome Assembly
Raw Illumina reads were quality trimmed (0.001 quality limit) and TruSeq adapters removed using CLC Genomics Workbench. Reads were imported as single or paired end reads. The paired end reads were merged to produce a merged dataset (Merged), while the single reads were used as unpaired (Single) and to produce a single-merged dataset (SM). Duplicates were also removed from these datasets to produce three duplicate removed datasets (Mddup, Sddup, and SMddup; Pienaar et al., 2021). These datasets were used to assemble the transcriptome using Trinity v2.4.0 and CLC Genomics Workbench v 20.0. Trinity was used with default parameter settings of kmer size 25. For CLC Genomics Workbench, kmer sizes were used in step sizes of 5 starting at 15 up to 60 and an additional assembly using kmer 64 (11 assemblies) with assembly parameters: mismatch cost-2, insertion cost-3, deletion cost-3, length fraction-0.9, similarity-0.9, minimum contig length-240, kmer size-variable, and bubble size-automatic. Given the different dataset structures used, a total of 72 assemblies were produced.
Extraction of the Mitochondrial Genome and Phylogenetic Analysis
The mitochondrial genome was identified in the assemblies by BLASTN analysis (Altschul et al., 1993), using the previously published mitochondrial genome for O. brasiliensis (Burger et al., 2014). The mitochondrial genome was annotated using the MITOS server to identify tRNA genes (Bernt et al., 2013). Protein coding and rRNA genes were identified using BLAST analysis (Altschul et al., 1993). The translated COI, CYTB, ND1, ND2, and ND4 proteins were used for phylogenetic analysis as previously described (Mans et al., 2015, 2019, 2021).
Extraction of CDS and Quality Assessment
Open reading frames (ORFs) were extracted using a Perl-script and chimeric and duplicate sequences removed by clustering at 90% identity using CD-HIT (Li and Godzik, 2006). The single dataset was mapped against the clustered ORFs using CLC Genomics Workbench and ORFs with TPM>1 (transcripts per million) were selected. BLASTX analysis against the ACARI database reduced the dataset further by selecting ORFs with E-values below 0.004 for further analysis. The transcriptome quality was measured for accuracy, completeness, contiguity, and chimerism using the Benchmarking Universal Single-Copy Orthologs (BUSCO; Simao et al., 2015). The final set of ORFs was submitted to GenBank under BioProject PRJNA719007.
Bioinformatic Analysis of the Transcriptome
To identify potential secretory peptides translated ORFs were submitted to SignalP (Petersen et al., 2011), while TMHMM and Phobius (Krogh et al., 2001; Kall et al., 2004) were used to identify membrane proteins. Potential housekeeping and secretory proteins were identified by BLASTP analysis against an ACARI database annotated using the KEGG database and GhostKOALA (Kanehisa et al., 2016a,b), the TSFAM database (Ribeiro and Mans, 2020), and an in-house annotation of secretory protein families (de Castro et al., 2016). To identify functional orthologs of proteins with experimentally verified functions, BLASTP analysis of secretory proteins was performed against the NCBI non-redundant database before phylogenetic analysis was performed to confirm clustering in functional clades with high bootstrap support. Protein families were aligned using MAFFT (Katoh and Standley, 2013), and maximum-likelihood analysis performed using IQ-Tree2 v 1.6.12 (Minh et al., 2020) with a standard 1,000,000 bootstraps.
Results and Discussion
Transcriptome Assembly
A total of 58,490,632 paired reads were generated that yielded 2,285,451 merged reads (Merged) after quality trimming and merging and 10,505,566 single reads (Single) after quality trimming, resulting in a combined 12,791,017 reads (SM). Removal of duplicate reads resulted in 1,301,316 merged reads (Mddup), 6,841,083 single reads (Sddup), and combined 8,142,399 reads (SMddup). Assembly resulted in 100,397 contigs after clustering with CDHIT that was further reduced to 44,052 contigs with TPM>1. This resulted in 16,908 contigs with E-values >0.004 and 14,957 contigs after manual curation. BUSCO analysis indicated 95.0% completeness with 92.6% as single genes, 2.4% as duplicated, 2.7% fragmented, and 2.3% missing from a set of 1,066 conserved genes. This compares well with other tick transcriptomes sequenced thus far (Figure 1A). Comparison to two other soft tick salivary gland transcriptomes (O. rostratus and Ornithodoros turicata which belong to the Neotropic and Nearctic Pavlovskyella), for which protein sequence data are available in the public repositories, indicates that O. brasiliensis generally has longer contigs although it also has higher numbers of short contigs, which may indicate that some contigs may be truncated (Figure 1B). Reciprocal best hit analysis indicated that O. brasiliensis shares 7,441 orthologs in total with these transcriptomes (Figure 1C). This can be considered a minimum, since this is limited by the numbers of contigs submitted for the other transcriptomes, i.e., O. rostratus (n = 6,602) and O. turicata (n = 7,544). As such, the percentage of orthologs shared for each transcriptome is 89% of O. rostratus and 82% of O. turicata. These measures were taken to indicate a well-represented high quality transcriptome. BLASTP analysis of ACARI database using the transcriptome retrieved as highest number of hits, proteins from related soft ticks, such as O. rostratus, O. turicata, O. erraticus, and O. moubata (Figure 2). This may be expected but is also a good measure of transcriptome quality.
Figure 1. Quality assessment of the transcriptome. (A) Comparison of the BUSCO analysis of all published salivary gland transcriptomes and Ornithodoros brasiliensis. The species are indicated on the horizontal axis, while the percentage of single complete, duplicated complete, missing, and fragmented is indicated on the vertical axis. Transcriptomes used in the comparison are as: H. dromedarii (Bensaoud et al., 2018), R. annulatus (Antunes et al., 2019), R. appendiculatus (de Castro et al., 2016), O. rostratus (Araujo et al., 2019), I. holocyclus (Rodriguez-Valle et al., 2018), O. turicata (Bourret et al., 2019), R. zambeziensis (de Castro et al., 2017), H. excavatum (Ribeiro et al., 2017), I. ricinus (Schwarz et al., 2013, 2014), A. triste (Garcia et al., 2014), R. pulchellus (Tan et al., 2015), A. sculptum (Moreira et al., 2017), A. maculatum (Karim et al., 2011), A. sculptum (Esteves et al., 2017), A. cajennense (Garcia et al., 2014), and A. americanum (Karim and Ribeiro, 2015). (B) Comparison of protein sequence length in amino acid residues presented up to 2,000 residues against the number of proteins. Proteins were binned in windows of 20 based on protein length. (C) Reciprocal best hit analysis of the transcriptome of O. brasiliensis (n = 14,957) against two closely related Neotropic and Nearctic tick species. Indicated are orthologs shared uniquely between species pairs or shared between all three species.
Figure 2. BLASTP analysis of 14,957 contigs against the ACARI database. Species with more than 10 hits are indicated.
Mitochondrial Genome Analysis
Previously, O. brasiliensis ticks shown to cause toxicoses and for which a mitochondrial genome was sequenced were collected at a site distant from the current study (Reck et al., 2013a; Burger et al., 2014). Since the collection sites were approximately 60 km apart, the question was raised regarding the genetic relationship of the ticks collected in this study and that of the mitochondrial genome previously sequenced for this species. The full-length mitochondrial genome was assembled in the transcriptome and was 99% identical to the previously sequenced genome and was deposited under accession number MW864544. Phylogenetic analysis indicated that both genomes cluster together with high support in a clade shared with O. rostratus in what has been described as the Neotropic Pavlovskyella clade (Figure 3). The sister-clade with high bootstrap support is that of the Afrotropic Ornithodoros sensu stricto group suggesting as previously indicated that these clades probably diverged during continental breakup of Gondwanaland ~127 MYA (Mans et al., 2019). This places the transcriptome within a phylogenetic context for comparative analysis and suggests that the transcriptomes of Afrotropic Ornithodoros and Neotropic Pavlovskyella should share extensive similarities with regard to orthologs and function. In addition, full-length 18S and 28S rRNA sequences were retrieved from the study and were deposited under GenBank accession numbers MW857182 and MW877711, respectively. These commonly used phylogenetic markers can therefore also be retrieved from transcriptome sequencing projects.
Figure 3. Phylogenetic analysis of argasid mitochondrial genes. Indicated is a maximum-likelihood tree of the protein coding genes COI, CYTB, ND1, ND2, and ND4 of the mitochondrial genome. Bootstrap values for 1,000,000 bootstraps are indicated. The tree was rooted with the Ixodidae. The mitochondrial genome obtained in the current study is underlined.
Composition of the Transcriptome
The transcriptome can be divided into housekeeping (n = 12,471), secretory (n = 368), and unknown categories (n = 2,118; Figure 4A). Housekeeping proteins are defined as all proteins not part of the secretory class that is involved in general cellular or organismal housekeeping functions, while secretory proteins are defined as those with secretory peptides, considered to be secreted into the feeding site during feeding, while unknown proteins refer to proteins with significant BLASTP hits in the database, but which has no annotation. These classifications have been used in all tick transcriptome studies dating back to Ribeiro et al. (2006). Housekeeping proteins account for ~83.3% of the transcriptome, secretory proteins for ~2.5%, and unknown proteins for ~14.1%. Reads mapped back to the transcriptome indicate that housekeeping proteins account for ~42% of the coverage, secretory proteins for ~53%, and unknowns for ~5%. This would suggest that secretory proteins are present at higher concentrations in the salivary glands relative to the housekeeping proteins. To corroborate this, the first 39 proteins with highest coverage are secretory proteins and comprise 43% of the total TPM coverage (80% of the secretory protein contribution; Figure 4B). The high abundance of secretory proteins was previously observed in soft tick salivary glands and may be explained by the fact that soft ticks synthesize secretory proteins and store them in large secretory granules that occupy most of the space in the salivary gland cells until secretion (Mans and Neitz, 2004; Mans et al., 2004). Abundance of secretory transcripts was also observed in conventional Sanger sequenced argasid salivary gland transcriptomes, as well as in those sequenced with next-generation sequencing technologies (Francischetti et al., 2008a,b; Mans et al., 2008a; Araujo et al., 2019).
Figure 4. Summary of transcriptome. (A) Comparison of the housekeeping, secretory, and unknown classes of contigs with regard to number of contigs contained in each class and the cumulative coverage of each class. (B) The first 500 contigs with the highest TPM coverage. Red dots indicate secretory proteins and black dots indicate housekeeping while white dots indicate unknown proteins.
Housekeeping Proteins
The housekeeping proteins can be divided into various functional classes as classified by the KEGG database (Figure 5). These include proteins involved in metabolism, protein synthesis (transcription and translation), protein folding, sorting, degradation and excretion (FSDE), environmental sensing processes, such as signal transduction systems, cellular processes, such as cell growth, death, and motility, and organismal processes, such as the circulatory, developmental, digestive, endocrine, excretory, immune, nervous, and sensory systems. Signal transduction processes are notably enriched with regard to number of proteins. The most abundant housekeeping classes include those involved in protein synthesis (transcription and translation), FSDE (Figure 6A), that cumulatively account for 54% of TPM coverage for housekeeping proteins. This may be expected for an organ, such as salivary glands whose primary function is the synthesis, folding, and sorting of secretory proteins.
Figure 5. Summary of housekeeping proteins as annotated by the GhostKOALA search of the KEGG database.
Figure 6. Abundance of housekeeping and secretory proteins expressed as percentage of total coverage (TPM). (A) For the housekeeping proteins, the classes with >1% total coverage (~90% of total coverage) are indicated. FSDE: folding, sorting, degradation, and excretion. (B) For the secretory proteins, the classes with >0.2% total coverage (~98% of total coverage) are indicated.
Some housekeeping functions with specific reference to tick biology may be highlighted. Of ~270 proteins identified in vertebrates to function as part of the protein secretory pathway (Gutierrez et al., 2020), 216 orthologs could be identified with confidence (Supplementary Table 1). The reason for many of those not identified may be due to multiple isoforms in the vertebrate secretory system which may be only represented by single proteins in ticks. This suggests that a large portion of the secretory system is present in the current transcriptome and is also conserved in ticks. This may be expected since the secretory pathway may be considered one of the Lineages of Life processes conserved in all Metazoa (Mans et al., 2016). Similarly, several orthologs of proteins previously implicated in neuronal control of salivary gland secretion in ticks (Šimo et al., 2014) were identified (Supplementary Table 2). This included hormones, such as bursicon, crustacean hyperglycemic hormone, eclosion hormone, elevenin, insulin, and orcokinin, as well as hormone receptors, such as allatostatin receptor, calcitonin gene-related peptide type 1 receptor, corazonin receptor 1, dopamine receptor 1, dopamine receptor 2-like, dopamine D2-like receptor, elevenin receptor 2, insulin-like growth factor 1 receptor, tachykinin-like peptides receptor, and periviscerokinin/Cap2b receptor.
It was previously indicated that the heme biosynthesis pathway is incomplete in ticks with ixodids lacking hemA-hemE, but possess hemF-hemH (Braz et al., 1999; Mans et al., 2016; Perner et al., 2016). Argasids possessed hemB and hemF-hemH (Mans et al., 2016). BLASTP analysis using hemA-hemH for Metaseiulus occidentalis indicated orthologs for hemB, hemF, and hemG, consistent with previous observations for argasids.
Secretory Proteins
Secretory protein families identified in the transcriptome included the majority of protein families generally found in tick salivary glands (Figures 6B, 7). The most abundant protein families in terms of sequence coverage (TPM) were the lipocalin, basic tail secretory family (BTSP), 7 disulfide bond domain (7DB), Kunitz-BPTI domains, 7 cysteine domain, and reprolysin, comprising ~98% of the total (Figure 6B). The lipocalin, basic tail, reprolysin, and Kunitz-BPTI families were also the most abundant in terms of the number of family members (Figure 7). This has previously been observed in soft tick salivary gland transcriptomes (Francischetti et al., 2008a,b; Mans et al., 2008a; Araujo et al., 2019). In addition, a number of peptides with secretory signals but no BLASTP hits were also identified as unknown proteins.
Figure 7. Summary of the secretory protein families found in the salivary gland transcriptome of O. brasiliensis. Indicated are numbers of proteins for each class.
Potential Proteins Functional at the Tick-Host Interface
Approximately 120 protein functions involved at the tick-host interface has been experimentally validated (Mans et al., 2019). BLASTP analysis using these proteins found 73 potential orthologs for 13 functions (Table 1). This includes the 5' nucleotidase family member apyrase that inhibits platelet aggregation by hydrolyzing ADP (Ribeiro et al., 1991; Mans et al., 1998; Stutzer et al., 2009).
Table 1. Orthologs for proteins with experimentally confirmed functions.
For the lipocalin family, a number of orthologs were found for proteins with known function, including leukotriene B4 (LTB4) scavengers (Mans and Ribeiro, 2008a). The LTB4 scavengers did not have the motifs conserved for complement C5 inhibition or thromboxane A2 scavenging and probably do not possess these functions (Nunn et al., 2005; Mans and Ribeiro, 2008a). Orthologs were also found for the biogenic amine scavengers and these orthologs possessed the biogenic amine-binding motif of the lower-binding site, while some also possessed the conserved amino acid residues involved in the upper-binding site of TSGP1 (Mans et al., 2008b). These orthologs therefore likely bind both histamine and serotonin. Orthologs for leukotriene C4 scavengers were also found (Mans and Ribeiro, 2008b).
In addition, an ortholog for the lipocalin (savicalin) was found. Savicalin was implicated in antimicrobial activity (Cheng et al., 2010). Other potential secretory antimicrobial orthologs present are defensins, microplusins, and trypsin inhibitor-like domains (Nakajima et al., 2001; Fogaça et al., 2004, 2006; Sasaki et al., 2008).
Four orthologs for adrenomedullin (vasodilation) were found in the transcriptome. Adrenomedullin was possibly acquired by soft ticks from the genus Ornithodoros by horizontal gene transfer from mammals and has thus far been found in O. coriaceus, O. moubata, O. parkeri, and O. rostratus (Iwanaga et al., 2014; Araujo et al., 2019). It would therefore seem as if the horizontal gene transfer occurred in the ancestral lineage to the Pavlovskyella and Ornithodoros (Figure 3). The presence of adrenomedullin in O. brasiliensis is therefore not surprising.
For members of the Kunitz-BPTI family, three orthologs of savignygrin were found. Single domain orthologs to the savignygrins were found that possessed the integrin RGD motif on the substrate binding presenting loop of the BPTI fold (Mans et al., 2002a). These inhibitors target the fibrinogen receptor integrin αIIbβ3 and inhibit platelet aggregation. Orthologs have also been found in A. monolakensis suggesting that these inhibitors were evolved in the ancestral argasid lineage (Mans et al., 2008c). A double-domain Kunitz-BPTI protein (Obras12157) with the same integrin RGD recognition motif located in the second Kunitz domain was also found. Such double Kunitz domain proteins were also observed in the transcriptomes of O. coriaceus and O. parkeri (Francischetti et al., 2008a,b). To date, no experimental evidence exists that these target integrin αIIbβ3, but it is likely that they target an integrin. The savignygrins are the Kunitz-BPTI proteins with the highest coverage for this family, as may be expected for an inhibitor that targets highly abundant platelet receptors in the host (Mans, 2019). Of interest is that no orthologs were found for the thrombin inhibitors like monobin, ornithodorin, or savignin (van de Locht et al., 1996; Nienaber et al., 1999; Mans et al., 2002b, 2008c), or the fXa inhibitors (Waxman et al., 1990; Gaspar et al., 1996; Joubert et al., 1998). This was also observed for Nearctic tick species, such as O. coriaceus and O. parkeri (Francischetti et al., 2008a,b). It was suggested that fXa inhibitors evolved exclusively within the genus Ornithodoros while thrombin inhibitors evolved in the ancestral argasid lineage (Mans et al., 2008c). Absence of these inhibitors in other Ornithodorinae lineages may indicate gene losses probably associated with host preferences. Possible orthologs to longistatin may be involved in blood coagulation modulation by activating plasminogen and in the immune system by targeting the receptor for advanced glycation end products (Anisuzzaman et al., 2011, 2014). Another ortholog associated with fibrinolysis is carboxypeptidase inhibitor that targets plasma carboxypeptidase B (thrombin-activatable fibrinolysis inhibitor) leading to fibrinolysis (Arolas et al., 2005). As such, a number of inhibitors that targets the blood coagulation cascade are present.
Orthologs of cystatin-2 from O. moubata were also detected. These cystatins inhibit cathepsin L and S allowing for modulation of the inflammatory responses in the host (Grunclová et al., 2006; Salát et al., 2010). Orthologs of serpins were also found that may inhibit elastase and cathepsin G acting as immunomodulators and platelet aggregation inhibitors (Prevot et al., 2006; Chmelar et al., 2011). It should be noted that the orthologs to which functions could be ascribed comprise ~19% of all the secretory proteins found in the transcriptome, suggesting that numerous undescribed functions still exist for this tick species.
Potential Proteins Involved in Tick Toxicoses of O. brasiliensis
The toxicoses and bite of O. brasiliensis are accompanied by pruritis, local edema and erythema, pain, and blisters, while histopathology of the feeding site indicates extensive subcutaneous edema and hemorrhage (Reck et al., 2013a, 2014). Possibly, salivary proteins that may contribute to this clinical outcome would include the clotting and platelet aggregation inhibitors that would induce hemorrhage. On the other hand, the LTB4 and LTC4 scavengers may inhibit edema and erythema (Mans and Ribeiro, 2008b). Additional contributors to edema and hemorrhage may be the secretory proteolytic enzymes found in the transcriptome that includes a cathepsin, two serine proteases, and 38 metalloproteases (astacin, gluzincin, and reprolysin). These may also impact in wound-healing processes leading to prolonged recovery times.
Evolutionary Perspectives on the Salivary Gland Transcriptome of O. brasiliensis
The salivary gland transcriptome of O. brasiliensis shows a number of evolutionary features previously found for argasids. This includes the presence of the major salivary gland secretory protein families conserved in all ticks (Mans et al., 2008a), as depicted in Figure 7. A number of orthologs for functions thus far conserved in argasids were also identified and include biogenic amine, leukotriene B4, and leukotriene C4 scavengers, as well as apyrase, savignygrin, and defensins (Mans and Ribeiro, 2008b). These functions are thus far conserved in all argasid species that feed on blood (Mans et al., 2016) and underscore their important role in blood-feeding. The only argasid salivary gland transcriptome sequenced to date that did not show these conserved features of protein families and conserved functions, were the salivary gland transcriptome for adult A. delacruzi (Ribeiro et al., 2012), which lacked all conserved features. This tick does not feed on a host in the adult phase and it is likely that this difference is due to differential expression in different life stages, and that these conserved features will be found in its blood-feeding stages. Since Antricola groups well within the Ornithodorinae, it may be expected that their lack of conserved protein families and functions in non-feeding phases is a derived trait. Conversely, the conserved protein families and functions as observed in O. brasiliensis and other argasids, indicate evolution of these traits in the last common ancestor to the Argasidae, or even the Ixodida (Mans et al., 2016).
Conclusion
The current study reported 14,957 unique transcripts for the nymphal salivary gland transcriptome of O. brasiliensis considered to be of high quality. The transcriptome is enriched with secretory proteins with high abundance that belongs to well-characterized secretory protein families. Several orthologs could be identified of experimentally verified functions, while the data indicated that numerous functions remain to be discovered. Some of the secretory proteins identified in the transcriptome could also be linked to the clinical disease syndrome caused by O. brasiliensis and opens new avenues to identify the proteins involved in toxicoses.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at https://www.ncbi.nlm.nih.gov/genbank/, SRR14139641.
Author Contributions
JR conceived the study, supplied the material, analyzed the data, and wrote the primary manuscript. AW supplied the material, analyzed the data, and critically reviewed the manuscript. BD supplied the material and critically reviewed the manuscript. RP, MC, and JF performed the experiments, analyzed the data, and critically reviewed the manuscript. BM conceived the study, performed the experiments, analyzed the data, and wrote the primary manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work is based on the research supported wholly/in part by the National Research Foundation of South Africa (grant numbers: 76499, 108998, 98862, and 118540) and the Economic Competitive Support Programme (30/01/V010).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fphys.2021.725635/full#supplementary-material
References
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1993). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1006/jmbi.1990.9999
Anisuzzaman, Islam, M. K., Alim, M. A., Miyoshi, T., Hatta, T., Yamaji, K., et al. (2011). Longistatin, a plasminogen activator, is key to the availability of blood-meals for ixodid ticks. PLoS Pathog. 7:e1001312. doi: 10.1371/journal.ppat.1001312
Anisuzzaman, H. T., Miyoshi, T., Matsubayashi, M., Islam, M. K., Alim, M. A., Anas, M. A., et al. (2014). Longistatin in tick saliva blocks advanced glycation end-product receptor activation. J. Clin. Invest. 124, 4429–4444. doi: 10.1172/jci74917
Antunes, S., Couto, J., Ferrolho, J., Sanches, G. S., Merino Charrez, J. O., De la Cruz Hernández, N., et al. (2019). Transcriptome and proteome response of Rhipicephalus annulatus tick vector to Babesia bigemina infection. Front. Physiol. 10:318. doi: 10.3389/fphys.2019.00318
Araujo, R. N., Silva, N. C. S., Mendes-Sousa, A., Paim, R., Costa, G. C. A., Dias, L. R., et al. (2019). RNA-seq analysis of the salivary glands and midgut of the Argasid tick Ornithodoros rostratus. Sci. Rep. 9:6764. doi: 10.1038/s41598-019-42899-z
Arolas, J. L., Lorenzo, J., Rovira, A., Castellà, J., Aviles, F. X., and Sommerhoff, C. P. (2005). A carboxypeptidase inhibitor from the tick Rhipicephalus bursa: isolation, cDNA cloning, recombinant expression, and characterization. J. Biol. Chem. 280, 3441–3448. doi: 10.1074/jbc.M411086200
Bensaoud, C., Nishiyama, M. Y. Jr., Ben Hamda, C., Lichtenstein, F., Castro de Oliveira, U., Faria, F., et al. (2018). De novo assembly and annotation of Hyalomma dromedarii tick (Acari: Ixodidae) sialotranscriptome with regard to gender differences in gene expression. Parasit. Vectors 11:314. doi: 10.1186/s13071-018-2874-9
Bernt, M., Donath, A., Jühling, F., Externbrink, F., Florentz, C., Fritzsch, G., et al. (2013). MITOS: improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 69, 313–319. doi: 10.1016/j.ympev.2012.08.023
Bourret, T. J., Boyle, W. K., Zalud, A. K., Valenzuela, J. G., Oliveira, F., and Lopez, J. E. (2019). The relapsing fever spirochete Borrelia turicatae persists in the highly oxidative environment of its soft-bodied tick vector. Cell. Microbiol. 21:e12987. doi: 10.1111/cmi.12987
Braz, G. R., Coelho, H. S., Masuda, H., and Oliveira, P. L. (1999). A missing metabolic pathway in the cattle tick Boophilus microplus. Curr. Biol. 9, 703–706. doi: 10.1016/s0960-9822(99)80312-1
Burger, T. D., Shao, R., Labruna, M. B., and Barker, S. C. (2014). Molecular phylogeny of soft ticks (Ixodida: Argasidae) inferred from mitochondrial genome and nuclear rRNA sequences. Ticks Tick Borne Dis. 5, 195–207. doi: 10.1016/j.ttbdis.2013.10.009
Cheng, P. H., Mans, B. J., Neitz, A. W., and Gaspar, A. R. (2010). Savicalin, a lipocalin from hemocytes of the soft tick, Ornithodoros savignyi. Exp. Appl. Acarol. 52, 313–326. doi: 10.1007/s10493-010-9368-6
Chmelar, J., Oliveira, C. J., Rezacova, P., Francischetti, I. M., Kovarova, Z., Pejler, G., et al. (2011). A tick salivary protein targets cathepsin G and chymase and inhibits host inflammation and platelet aggregation. Blood 117, 736–744. doi: 10.1182/blood-2010-06-293241
Dall’Agnol, B., Schott, D., Padilha, T., Antunes, P., Souza, U. A., Webster, A., et al. (2019). Clinical findings associated with Ornithodoros brasiliensis tick parasitism in travelers, southern Brazil. Wilderness Environ. Med. 30, 437–440. doi: 10.1016/j.wem.2019.06.012
de Castro, M. H., de Klerk, D., Pienaar, R., Latif, A. A., Rees, D. J., and Mans, B. J. (2016). De novo assembly and annotation of the salivary gland transcriptome of Rhipicephalus appendiculatus male and female ticks during blood feeding. Ticks Tick Borne Dis. 7, 536–548. doi: 10.1016/j.ttbdis.2016.01.014
de Castro, M. H., de Klerk, D., Pienaar, R., Rees, D. J. G., and Mans, B. J. (2017). Sialotranscriptomics of Rhipicephalus zambeziensis reveals intricate expression profiles of secretory proteins and suggests tight temporal transcriptional regulation during blood-feeding. Parasit. Vectors 10:384. doi: 10.1186/s13071-017-2312-4
Esteves, E., Maruyama, S. R., Kawahara, R., Fujita, A., Martins, L. A., Righi, A. A., et al. (2017). Analysis of the salivary gland transcriptome of unfed and partially fed Amblyomma sculptum ticks and descriptive proteome of the saliva. Front. Cell. Infect. Microbiol. 7:476. doi: 10.3389/fcimb.2017.00476
Fogaça, A. C., Almeida, I. C., Eberlin, M. N., Tanaka, A. S., Bulet, P., and Daffre, S. (2006). Ixodidin, a novel antimicrobial peptide from the hemocytes of the cattle tick Boophilus microplus with inhibitory activity against serine proteinases. Peptides 27, 667–674. doi: 10.1016/j.peptides.2005.07.013
Fogaça, A. C., Lorenzini, D. M., Kaku, L. M., Esteves, E., Bulet, P., and Daffre, S. (2004). Cysteine-rich antimicrobial peptides of the cattle tick Boophilus microplus: isolation, structural characterization and tissue expression profile. Dev. Comp. Immunol. 28, 191–200. doi: 10.1016/j.dci.2003.08.001
Francischetti, I. M., Mans, B. J., Meng, Z., Gudderra, N., Veenstra, T. D., Pham, V. M., et al. (2008b). An insight into the sialome of the soft tick, Ornithodorus parkeri. Insect Biochem. Mol. Biol. 38, 1–21. doi: 10.1016/j.ibmb.2007.09.009
Francischetti, I. M., Meng, Z., Mans, B. J., Gudderra, N., Hall, M., Veenstra, T. D., et al. (2008a). An insight into the salivary transcriptome and proteome of the soft tick and vector of epizootic bovine abortion, Ornithodoros coriaceus. J. Proteome 71, 493–512. doi: 10.1016/j.jprot.2008.07.006
Francischetti, I. M., Sa-Nunes, A., Mans, B. J., Santos, I. M., and Ribeiro, J. M. (2009). The role of saliva in tick feeding. Front. Biosci. (Landmark Ed) 14, 2051–2088. doi: 10.2741/3363
Garcia, G. R., Gardinassi, L. G., Ribeiro, J. M., Anatriello, E., Ferreira, B. R., Moreira, H. N., et al. (2014). The sialotranscriptome of Amblyomma triste, Amblyomma parvum and Amblyomma cajennense ticks, uncovered by 454-based RNA-seq. Parasit. Vectors 7:430. doi: 10.1186/1756-3305-7-430
Gaspar, A. R., Joubert, A. M., Crause, J. C., and Neitz, A. W. (1996). Isolation and characterization of an anticoagulant from the salivary glands of the tick, Ornithodoros savignyi (Acari: Argasidae). Exp. Appl. Acarol. 20, 583–598. doi: 10.1007/BF00052809
Grunclová, L., Horn, M., Vancová, M., Sojka, D., Franta, Z., Mares, M., et al. (2006). Two secreted cystatins of the soft tick Ornithodoros moubata: differential expression pattern and inhibitory specificity. Biol. Chem. 387, 1635–1644. doi: 10.1515/BC.2006.204
Guglielmone, A. A., Robbins, R. G., Apanaskevich, D. A., Petney, T. N., Estrada-Peña, A., Horak, I. G., et al. (2010). The Argasidae, Ixodidae and Nuttalliellidae (Acari: Ixodida) of the world: a list of valid species names. Zootaxa 2528, 1–28. doi: 10.11646/zootaxa.2528.1.1.1
Gutierrez, J. M., Feizi, A., Li, S., Kallehauge, T. B., Hefzi, H., Grav, L. M., et al. (2020). Genome-scale reconstructions of the mammalian secretory pathway predict metabolic costs and limitations of protein secretion. Nat. Commun. 11:68. doi: 10.1038/s41467-019-13867-y
Hoogstraal, H. (1985). Argasid and nuttalliellid ticks as parasites and vectors. Adv. Parasitol. 24, 135–238. doi: 10.1016/s0065-308x(08)60563-1
Iwanaga, S., Isawa, H., and Yuda, M. (2014). Horizontal gene transfer of a vertebrate vasodilatory hormone into ticks. Nat. Commun. 5;3373. doi: 10.1038/ncomms4373
Joubert, A. M., Louw, A. I., Joubert, F., and Neitz, A. W. (1998). Cloning, nucleotide sequence and expression of the gene encoding factor Xa inhibitor from the salivary glands of the tick, Ornithodoros savignyi. Exp. Appl. Acarol. 22, 603–619. doi: 10.1023/a:1006198713791
Kall, L., Krogh, A., and Sonnhammer, E. L. (2004). A combined transmembrane topology and signal peptide prediction method. J. Mol. Biol. 338, 1027–1036. doi: 10.1016/j.jmb.2004.03.016
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016a). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi: 10.1093/nar/gkv1070
Kanehisa, M., Sato, Y., and Morishima, K. (2016b). BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol. 428, 726–731. doi: 10.1016/j.jmb.2015.11.006
Karim, S., Singh, P., and Ribeiro, J. M. (2011). A deep insight into the sialotranscriptome of the gulf coast tick, Amblyomma maculatum. PLoS One 6:e28525. doi: 10.1371/journal.pone.0028525
Karim, S., and Ribeiro, J. M. (2015). An insight into the sialome of the Lone Star tick, Amblyomma americanum, with a glimpse on its time dependent gene expression. PLoS One 10:e0131292. doi: 10.1371/journal.pone.0131292
Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Krogh, A., Larsson, B., von Heijne, G., and Sonnhammer, E. L. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J. Mol. Biol. 305, 567–580. doi: 10.1006/jmbi.2000.4315
Li, W., and Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. doi: 10.1093/bioinformatics/btl158
Mans, B. J. (2019). Chemical equilibrium at the tick-host feeding interface: a critical examination of biological relevance in hematophagous behavior. Front. Physiol. 10:530. doi: 10.3389/fphys.2019.00530
Mans, B. J. (2020). Quantitative visions of reality at the tick-host interface: biochemistry, genomics, proteomics, and transcriptomics as measures of complete inventories of the tick sialoverse. Front. Cell. Infect. Microbiol. 10:574405. doi: 10.3389/fcimb.2020.574405
Mans, B. J., Andersen, J. F., Francischetti, I. M., Valenzuela, J. G., Schwan, T. G., Pham, V. M., et al. (2008a). Comparative sialomics between hard and soft ticks: implications for the evolution of blood-feeding behavior. Insect Biochem. Mol. Biol. 38, 42–58. doi: 10.1016/j.ibmb.2007.09.003
Mans, B. J., Andersen, J. F., Schwan, T. G., and Ribeiro, J. M. (2008c). Characterization of anti-hemostatic factors in the argasid, Argas monolakensis: implications for the evolution of blood-feeding in the soft tick family. Insect Biochem. Mol. Biol. 38, 22–41. doi: 10.1016/j.ibmb.2007.09.002
Mans, B. J., de Castro, M. H., Pienaar, R., de Klerk, D., Gaven, P., Genu, S., et al. (2016). Ancestral reconstruction of tick lineages. Ticks Tick Borne Dis. 7, 509–535. doi: 10.1016/j.ttbdis.2016.02.002
Mans, B. J., de Klerk, D., Pienaar, R., de Castro, M. H., and Latif, A. A. (2015). Next-generation sequencing as means to retrieve tick systematic markers, with the focus on Nuttalliella namaqua (Ixodoidea: Nuttalliellidae). Ticks Tick Borne Dis. 6, 450–462. doi: 10.1016/j.ttbdis.2015.03.013
Mans, B. J., Featherston, J., Kvas, M., Pillay, K. A., de Klerk, D. G., Pienaar, R., et al. (2019). Argasid and ixodid systematics: implications for soft tick evolution and systematics, with a new argasid species list. Ticks Tick Borne Dis. 10, 219–240. doi: 10.1016/j.ttbdis.2018.09.010
Mans, B. J., Gasper, A. R., Louw, A. I., and Neitz, A. W. (1998). Purification and characterization of apyrase from the tick, Ornithodoros savignyi. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 120, 617–624. doi: 10.1016/s0305-0491(98)10061-5
Mans, B. J., Kelava, S., Pienaar, R., Featherston, J., de Castro, M. H., Quetglas, J., et al. (2021). Nuclear (18S-28S rRNA) and mitochondrial genome markers of Carios (Carios) vespertilionis (Argasidae) support Carios Latreille, 1796 as a lineage embedded in the Ornithodorinae: re-classification of the Carios sensu Klompen and Oliver (1993) clade into its respective subgenera. Ticks Tick Borne Dis 12:101688. doi: 10.1016/j.ttbdis.2021.101688
Mans, B. J., Louw, A. I., and Neitz, A. W. (2002a). Savignygrin, a platelet aggregation inhibitor from the soft tick Ornithodoros savignyi, presents the RGD integrin recognition motif on the Kunitz-BPTI fold. J. Biol. Chem. 277, 21371–21378. doi: 10.1074/jbc.M112060200
Mans, B. J., Louw, A. I., and Neitz, A. W. (2002b). Amino acid sequence and structure modeling of savignin, a thrombin inhibitor from the tick, Ornithodoros savignyi. Insect Biochem. Mol. Biol. 32, 821–828. doi: 10.1016/s0965-1748(01)00169-2
Mans, B. J., and Neitz, A. W. (2004). Molecular crowding as a mechanism for tick secretory granule biogenesis. Insect Biochem. Mol. Biol. 34, 1187–1193. doi: 10.1016/j.ibmb.2004.07.007
Mans, B. J., and Ribeiro, J. M. (2008a). Function, mechanism and evolution of the moubatin-clade of soft tick lipocalins. Insect Biochem. Mol. Biol. 38, 841–852. doi: 10.1016/j.ibmb.2008.06.007
Mans, B. J., and Ribeiro, J. M. (2008b). A novel clade of cysteinyl leukotriene scavengers in soft ticks. Insect Biochem. Mol. Biol. 38, 862–870. doi: 10.1016/j.ibmb.2008.06.002
Mans, B. J., Ribeiro, J. M., and Andersen, J. F. (2008b). Structure, function, and evolution of biogenic amine-binding proteins in soft ticks. J. Biol. Chem. 283, 18721–18733. doi: 10.1074/jbc.M800188200
Mans, B. J., Venter, J. D., Coons, L. B., Louw, A. I., and Neitz, A. W. (2004). A reassessment of argasid tick salivary gland ultrastructure from an immuno-cytochemical perspective. Exp. Appl. Acarol. 33, 119–129. doi: 10.1023/b:appa.0000030012.47964.b3
Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015
Moreira, H. N. S., Barcelos, R. M., Vidigal, P. M. P., Klein, R. C., Montandon, C. E., Maciel, T. E. F., et al. (2017). A deep insight into the whole transcriptome of midguts, ovaries and salivary glands of the Amblyomma sculptum tick. Parasitol. Int. 66, 64–73. doi: 10.1016/j.parint.2016.10.011
Nakajima, Y., van der Goes van Naters-Yasui, A., Taylor, D., and Yamakawa, M. (2001). Two isoforms of a member of the arthropod defensin family from the soft tick, Ornithodoros moubata (Acari: Argasidae). Insect Biochem. Mol. Biol. 31, 747–751. doi: 10.1016/s0965-1748(01)00066-2
Nienaber, J., Gaspar, A. R., and Neitz, A. W. (1999). Savignin, a potent thrombin inhibitor isolated from the salivary glands of the tick Ornithodoros savignyi (Acari: Argasidae). Exp. Parasitol. 93, 82–91. doi: 10.1006/expr.1999.4448
Nunn, M. A., Sharma, A., Paesen, G. C., Adamson, S., Lissina, O., Willis, A. C., et al. (2005). Complement inhibitor of C5 activation from the soft tick Ornithodoros moubata. J. Immunol. 174, 2084–2091. doi: 10.4049/jimmunol.174.4.2084
Oleaga, A., Soriano, B., Llorens, C., and Pérez-Sánchez, R. (2021). Sialotranscriptomics of the argasid tick Ornithodoros moubata along the trophogonic cycle. PLoS Negl. Trop. Dis. 15:e0009105. doi: 10.1371/journal.pntd.0009105
Pérez-Sánchez, R., Carnero-Morán, Á., Soriano, B., Llorens, C., and Oleaga, A. (2021). RNA-seq analysis and gene expression dynamics in the salivary glands of the argasid tick Ornithodoros erraticus along the trophogonic cycle. Parasit. Vectors 14:170. doi: 10.1186/s13071-021-04671-z
Perner, J., Sobotka, R., Sima, R., Konvickova, J., Sojka, D., Oliveira, P. L., et al. (2016). Acquisition of exogenous haem is essential for tick reproduction. eLife 5:e12318. doi: 10.7554/eLife.12318
Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786. doi: 10.1038/nmeth.1701
Pienaar, R., de Klerk, D. G., de Castro, M. H., Featherston, J., and Mans, B. J. (2021). De novo assembled salivary gland transcriptome and expression pattern analyses for Rhipicephalus evertsi evertsi Neuman, 1897 male and female ticks. Sci. Rep. 11:1642. doi: 10.1038/s41598-020-80454-3
Prevot, P. P., Adam, B., Boudjeltia, K. Z., Brossard, M., Lins, L., Cauchie, P., et al. (2006). Anti-hemostatic effects of a serpin from the saliva of the tick Ixodes ricinus. J. Biol. Chem. 281, 26361–26369. doi: 10.1074/jbc.M604197200
Ramirez, D. G., Landulfo, G. A., Onofrio, V. C., Simons, S. M., Reck, J., Martins, J. R., et al. (2016). Laboratory life cycle of Ornithodoros brasiliensis (Acari: Argasidae): An endemic tick from southern Brazil. Ticks Tick Borne Dis. 7, 730–733. doi: 10.1016/j.ttbdis.2016.03.001
Reck, J., Bandarra, P., Pavarini, S., Termignoni, C., Driemeier, D., Martins, J. R., et al. (2014). Experimentally induced tick toxicosis in rats bitten by Ornithodoros brasiliensis (Chelicerata: Argasidae): a clinico-pathological characterization. Toxicon 88, 99–106. doi: 10.1016/j.toxicon.2014.06.017
Reck, J., Marks, F. S., Guimarães, J. A., Termignoni, C., and Martins, J. R. (2013a). Epidemiology of Ornithodoros brasiliensis (mouro tick) in the southern Brazilian highlands and the description of human and animal retrospective cases of tick parasitism. Ticks Tick Borne Dis. 4, 101–109. doi: 10.1016/j.ttbdis.2012.09.004
Reck, J., Marks, F. S., Termignoni, C., Guimarães, J. A., and Martins, J. R. (2013b). Ornithodoros brasiliensis (mouro tick) salivary gland homogenates inhibit in vivo wound healing and in vitro endothelial cell proliferation. Parasitol. Res. 112, 1749–1753. doi: 10.1007/s00436-013-3333-3
Reck, J., Soares, J. F., Termignoni, C., Labruna, M. B., and Martins, J. R. (2011). Tick toxicosis in a dog bitten by Ornithodoros brasiliensis. Vet. Clin. Pathol. 40, 356–360. doi: 10.1111/j.1939-165X.2011.00338.x
Ribeiro, J. M., Alarcon-Chaidez, F., Francischetti, I. M., Mans, B. J., Mather, T. N., Valenzuela, J. G., et al. (2006). An annotated catalog of salivary gland transcripts from Ixodes scapularis ticks. Insect Biochem. Mol. Biol. 36, 111–129. doi: 10.1016/j.ibmb.2005.11.005
Ribeiro, J. M., Endris, T. M., and Endris, R. (1991). Saliva of the soft tick, Ornithodoros moubata, contains anti-platelet and apyrase activities. Comp. Biochem. Physiol. A Comp. Physiol. 100, 109–112. doi: 10.1016/0300-9629(91)90190-n
Ribeiro, J. M., Labruna, M. B., Mans, B. J., Maruyama, S. R., Francischetti, I. M., Barizon, G. C., et al. (2012). The sialotranscriptome of Antricola delacruzi female ticks is compatible with non-hematophagous behavior and an alternative source of food. Insect Biochem. Mol. Biol. 42, 332–342. doi: 10.1016/j.ibmb.2012.01.003
Ribeiro, J. M., Slovák, M., and Francischetti, I. M. (2017). An insight into the sialome of Hyalomma excavatum. Ticks Tick Borne Dis. 8, 201–207. doi: 10.1016/j.ttbdis.2016.08.011
Ribeiro, J. M. C., and Mans, B. J. (2020). TickSialoFam (TSFam): A database that helps to classify tick salivary proteins, a review on tick salivary protein function and evolution, with considerations on the tick sialome switching phenomenon. Front. Cell. Infect. Microbiol. 10:374. doi: 10.3389/fcimb.2020.00374
Rodriguez-Valle, M., Moolhuijzen, P., Barrero, R. A., Ong, C. T., Busch, G., Karbanowicz, T., et al. (2018). Transcriptome and toxin family analysis of the paralysis tick, Ixodes holocyclus. Int. J. Parasitol. 48, 71–82. doi: 10.1016/j.ijpara.2017.07.007
Salát, J., Paesen, G. C., Rezácová, P., Kotsyfakis, M., Kovárová, Z., Sanda, M., et al. (2010). Crystal structure and functional characterization of an immunomodulatory salivary cystatin from the soft tick Ornithodoros moubata. Biochem. J. 429, 103–112. doi: 10.1042/BJ20100280
Sasaki, S. D., de Lima, C. A., Lovato, D. V., Juliano, M. A., Torquato, R. J., and Tanaka, A. S. (2008). BmSI-7, a novel subtilisin inhibitor from Boophilus microplus, with activity toward Pr1 proteases from the fungus Metarhizium anisopliae. Exp. Parasitol. 118, 214–220. doi: 10.1016/j.exppara.2007.08.003
Schwarz, A., Tenzer, S., Hackenberg, M., Erhart, J., Gerhold-Ay, A., Mazur, J., et al. (2014). A systems level analysis reveals transcriptomic and proteomic complexity in Ixodes ricinus midgut and salivary glands during early attachment and feeding. Mol. Cell Proteom. 13, 2725–2735. doi: 10.1074/mcp.M114.039289
Schwarz, A., von Reumont, B. M., Erhart, J., Chagas, A. C., Ribeiro, J. M., and Kotsyfakis, M. (2013). De novo Ixodes ricinus salivary gland transcriptome analysis using two next-generation sequencing methodologies. FASEB J. 27, 4745–4756. doi: 10.1096/fj.13-232140
Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Šimo, L., Sonenshine, D. E., Park, Y., and Žitňan, D. (2014). “Nervous and sensory systems: structure, function, genomics and proteomics,” in Biology of ticks, 2nd Edn. Vol. 1. eds. D. E. Sonenshine and R. M. Roe (New York: Oxford University Press), 309–367.
Stutzer, C., Mans, B. J., Gaspar, A. R., Neitz, A. W., and Maritz-Olivier, C. (2009). Ornithodoros savignyi: soft tick apyrase belongs to the 5'-nucleotidase family. Exp. Parasitol. 122, 318–327. doi: 10.1016/j.exppara.2009.04.007
Tan, A. W., Francischetti, I. M., Slovak, M., Kini, R. M., and Ribeiro, J. M. (2015). Sexual differences in the sialomes of the zebra tick. Rhipicephalus pulchellus. J. Proteomics 117, 120–144. doi: 10.1016/j.jprot.2014.12.014
van de Locht, A., Stubbs, M. T., Bode, W., Friedrich, T., Bollschweiler, C., Höffken, W., et al. (1996). The ornithodorin-thrombin crystal structure, a key to the TAP enigma? EMBO J. 15, 6011–6017. doi: 10.1002/j.1460-2075.1996.tb00989.x
Keywords: RNA, next-generation sequencing, argasid, nymph, secretory proteins
Citation: Reck J, Webster A, Dall’Agnol B, Pienaar R, de Castro MH, Featherston J and Mans BJ (2021) Transcriptomic Analysis of Salivary Glands of Ornithodoros brasiliensis Aragão, 1923, the Agent of a Neotropical Tick-Toxicosis Syndrome in Humans. Front. Physiol. 12:725635. doi: 10.3389/fphys.2021.725635
Edited by:
Tetsuya Tanaka, Kagoshima University, JapanReviewed by:
Michail Kotsyfakis, Academy of Sciences of the Czech Republic (ASCR), CzechiaAlejandro Cabezas-Cruz, Institut National de Recherche Pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), France
Copyright © 2021 Reck, Webster, Dall’Agnol, Pienaar, de Castro, Featherston and Mans. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jose Reck, am9zZS5yZWNrQGdtYWlsLmNvbQ==; Ben J. Mans, bWFuc2JAYXJjLmFncmljLnph
†These authors have contributed equally to this work