Comparative Tandem Mass Tag-Based Quantitative Proteomic Analysis of Tachaea chinensis Isopod During Parasitism

Parasitic isopods perforate and attach to the host integument via the mandibles and then feed on hemolymph and exudate from the wounds. Such isopods attack a variety of commercially important fish and crustacean hosts. Similar to other hematophagous parasites, isopods may also employ biomolecules that affect host blood conglutination and defense systems. In the present study, a tandem mass tag-based quantitative proteomic approach was used to identify differentially expressed proteins in Tachaea chinensis parasites of shrimp, by comparing parasitic (fed) and pre-parasitic (unfed) individuals. We identified 888 proteins from a total of 1,510 peptides, with a significant difference in 129 between the fed and unfed groups. Among these, 37 were upregulated and 92 were downregulated in unfed T. chinensis. This indicates that T. chinensis may require more energy before parasitism during its search for a host. In addition, as is the case for other blood-sucking parasites, it might secrete antihemostatic, anti-inflammatory, and immunomodulatory molecules to facilitate blood meal acquisition. To our knowledge, this study is the first to use a TMT-based proteomic approach to analyze the proteome of isopod parasites, and the results will facilitate our understanding of the molecular mechanisms of isopod parasitism on crustaceans.

Introduction

Hematophagous parasites such as ticks, mosquitoes, fleas, and leeches can bypass the clotting system of their vertebrate hosts by secreting anticoagulant compounds (Shu et al., 2010; Santiago et al., 2017). Isopod parasites mainly feed on blood (hemolymph) from hosts after perforating the integument with their mandibles, which can cause stress, tissue damage, secondary infection, and mortality (Wilson, 2008; Poore and Bruce, 2012; Williams and Boyko, 2012; Smit et al., 2014). There is evidence that isopod parasites may also inject anticoagulants or other compounds directly into the host so as to evade defensive responses (Nair and Nair, 1983; Manship et al., 2012). At present, only a few proteins for this purpose have been identified, including trypsin inhibitors and anticoagulants in Paragnathia formica. Identification of isopod anticoagulant agents may be useful in fully understanding parasite modulation of the host physiological system. This information may be valuable for the identification of novel target antigens and antithrombotic drug development.

Tachaea chinensis, one of the most common ectoparasites of farm-raised shrimps, is widely distributed in China and neighboring countries (Lu et al., 2007; Hua et al., 2018; Li et al., 2018). It is normally approximately 0.7 cm in size and can easily be observed with the naked eye following parasitization. After attacking its host, this isopod invariably remains attached until either the shrimp dies or another shrimp comes within close proximity. Over the past 2 years, T. chinensis has been shown to parasitize over 80% of aquacultured Palaemonetes sinensis (Li et al., 2018) and 90% of Exopalaemon carinicauda (Li et al., 2019). Additional research on isopod parasitization mechanisms is necessary to develop effective methods of control for use in shrimp aquaculture.

We performed tandem mass tag (TMT)-based quantitative proteomic profiling of T. chinensis to compare protein expression in fed and unfed parasites and gain insight into their integrated molecular mechanisms and host responses. This will provide an empirical basis for disease prevention and control efforts and support further research on the molecular biology of isopods.

Materials and Methods

Ethics Statement

Our study did not involve endangered or protected species. In China, the capture of isopod parasites and their host shrimp from rice fields does not require specific permits. All efforts were made to minimize animal suffering and discomfort. The experimental protocol was approved by the Animal Ethics Committee of Shenyang Agriculture University.

Animals

Tachaea chinensis isopod parasites (0.82 ± 0.17 cm) and their host shrimp P. sinensis (3.48 ± 0.35 g) used in this study were collected from a rice field in Panjin City, Liaoning Province, China, in November 2018, and transported to the aquaculture laboratory at Shenyang Agricultural University. They were acclimated in two 300 L², fiberglass recirculation tanks with a circular flow system. Water temperature was maintained at 24 ± 0.5°C, and the photoperiod was set at a 12:12 h light:dark cycle. After 2 weeks of acclimatization, 20 healthy T. chinensis were transferred to individual plastic tanks (15.8 cm diameter and 32.1 cm height), each containing 5 L water, with the same environmental conditions as during the acclimatization period.

Following this, one P. sinensis each was placed in 10 of the prepared tanks, forming the fed group, with the other tanks kept as unfed pre-parasitism controls. After 7 d, fed and unfed T. chinensis were removed to individual 2 mL RNAse-free tubes and immediately frozen in liquid nitrogen for storage until protein extraction.

Protein Extraction and SDS-PAGE Analysis

Three fed and three unfed isopods were ground into powder and vortexed in 600 μL SDT buffer (pH 8.0, 4% SDS, 150 mM Tris-HCl, 1 mM DTT), respectively. The mixtures were heated at 100°C for 10 min, then sonicated at 35 W for 4 s, with 7 s intervals, for 10 min. These were centrifuged at 14,000 × g for 30 min, and supernatants were collected into 0.22 μm filter tubes. One microliter of the underlayer liquid of each sample was used for BCA quantitative analysis, and 1 μg of the protein sample from each group was subjected to SDS-PAGE (12.5% resolving gels) analysis (Figure S1).

Protein Digestion and TMT Labeling

Briefly, six protein concentrates (300 μg each) were mixed in ultrafiltration filtrate tubes (30 kDa cut-off, Sartorius, Gottingen, Germany) with 200 μL urea buffer (8 M urea, 150 mM Tris-HCl, pH 8.0), and the sample was centrifuged at 14,000 × g at 20°C for 30 min. The sample was washed twice by adding 200 μL UA and centrifuged at 14,000 × g at 20°C for 30 min. The flow-through from the collection tube was discarded. Next, 100 μL of indole-3-acetic acid (IAA) solution (50 mM IAA in UA buffer) was added to the filter tube and vortexed at 600 rpm in a Thermomixer comfort incubator (Eppendorf, Germany) for 1 min. Subsequently, the sample was incubated at room temperature for 30 min in the dark and centrifuged at 14,000 × g for 30 min at 20°C. Next, 100 μL UA was added to the filter unit, which was centrifuged at 14,000 × g for 20 min; this was carried out three times. The protein suspension in the filtrate tube was subjected to enzyme digestion with 52 μL of trypsin (Promega, Madison, WI, USA) buffer [6 μg trypsin (0.5 μg/μL) in 40 μL of dissolution buffer] for 16–18 h at 37°C. Finally, the filter unit was transferred to a new tube and centrifuged at 14,000 × g for 30 min. Peptides were collected in the filtrate and the peptide concentration was measured based on the optical density at 280 nm (OD280).

TMT labeling was performed using the TMT6plex™ Isobaric Label Reagent Set (Thermo Scientific) according to the manufacturer's instructions. The proteins in the unfed group were labeled with reagents 119, 127, and 128, and those in the infected group were labeled with reagents 129, 130, and 131. The labeling solution reaction was then incubated at room temperature for 1.5 h prior to further analysis.

SCX Fractionation and LC-MS/MS Analysis

SCX fractionation and LC-MS/MS analysis were conducted according to a published protocol (Xiao et al., 2019). The TMT-labeled samples were analyzed using an Easy-nLC nanoflow HPLC system connected to Orbitrap-Elite (Thermo Fisher Scientific, San Jose, CA, USA). One microgram of each sample was loaded onto each of two Thermo Scientific EASY columns using an autosampler at a flow rate of 200 nL/min. The sequential separation of peptides on Thermo Scientific EASY trap column (100 μm × 2 cm, 5 μm, 100 Å, C18) and analytical column (75 μm × 25 cm, 5 μm, 100 Å, C18) was achieved using a segmented 1 h gradient from 5 to 28% Solvent B (0.1% formic acid in 100% ACN) for 40 min, followed by 28–90% Solvent B for 2 min, and 90% Solvent B for 18 min. The column was re-equilibrated to its initial highly aqueous solvent composition before each analysis. The mass spectrometer was operated in positive ion mode, and MS spectra were acquired over a range of 350–2,000 m/z. The resolving powers of the MS scan and MS/MS scan at 100 m/z for the Orbitrap Elite were set as 60,000 and 15,000, respectively. The top 16 most intense signals in the acquired MS spectra were selected for further MS/MS analysis. The isolation window was 2 m/z, and ions were fragmented through high energy collisional dissociation with normalized collision energies of 35 eV. The maximum ion injection times were set at 10 ms for the survey scan and 100 ms for the MS/MS scans, and the automatic gain control target values for full scan modes was set to 1 e 6 and for MS/MS was 5 e 4. The dynamic exclusion duration was 30 s.

Data Analysis

Raw files were analyzed using the Proteome Discoverer 2.1 software (Thermo Fisher Scientific). A search for the fragmentation spectra was performed using the MASCOT search engine embedded in Proteome Discoverer against the NCBI_Peracarida_91190_20190313.fasta database. The following search parameters were used: monoisotopic mass, trypsin as the cleavage enzyme, two missed cleavages, TMT labeling, and carbamidomethylating of cysteine as fixed modifications. Peptide charges of 2⁺, 3⁺, and 4⁺ and the oxidation of methionine were specified as variable modifications. The mass tolerance was set to 20 ppm for precursor ions and to 0.1 Da for fragment ions. The results were filtered based on a false discovery rate of no more than 1%. The relative quantitative analysis of the proteins in the samples based on the ratios of TMT reporter ions from all unique peptides representing each protein was performed using Proteome Discoverer (version 2.1). The relative peak intensities of the TMT reporter ions released in each of the MS/MS spectra were used. Then, the final ratios obtained from the relative protein quantifications were normalized based on the median average protein quantification ratio. The fold change was set to >1.2 for protein upregulation and <0.85 for protein downregulation (P < 0.05). Protein functional annotation was conducted using the Universal Protein (UniProt) database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were also used.

The mass spectrometry proteomics data have been deposited in the ProteomeXchange Consortium via the PRIDE partner repository (identifier PXD015247). All analyzed data are available from the corresponding author upon reasonable request.

Verification of Protein Quantifications Using PRM Analysis

Parallel reaction monitoring (PRM) was used to verify the TMT-based quantitative proteomics results. Briefly, 2 μg of peptide from each sample was taken for LC-PRM/MS analysis. After sample loading, chromatographic separation was performed using a Thermo Scientific EASY-nLC nano-HPLC system, with two buffers. Solution A was 0.1% formic acid aqueous solution and solution B was a mixed solution of 0.1% formic acid, 95% acetonitrile, and water. The column was first equilibrated with 95% solution A. The sample was injected into a Trap column (100 μm × 20 mm, 5 μm C18, Dr. Maisch GmbH) and subjected to gradient separation through a chromatography column (75 μm × 150 mm, 3 μm C18, Dr. Maisch GmbH) at a flow rate of 250 nL/min. The liquid phase separation gradient was as follows: 0–25 min, linear gradient of B liquid from 5 to 18%; 25–45 min, linear gradient of B liquid from 18 to 50%; 45–48 min, linear gradient of B liquid from 50 to 95%; and 48–60 min, B liquid maintained at 95%. The peptides were separated and subjected to targeted PRM/MS using a Q-Exactive mass spectrometer (Thermo Scientific) for 60 min. The parameters were set as follows: detection mode, positive mode; parent ion scanning range, 350–1,500 m/z; capillary voltage, 1.8 kV; isolation width, 1.6 Th; first-order MS resolution, 70,000 at m/z 200; AGC target, 3 e 6; first-level maximum IT, 250 ms. Peptide secondary MS were obtained as follows: for each full scan, target peptides of precursor m/z were sequentially selected based on the inclusion list for second-order MS (MS2) scan. The parameters used were as follows: resolution, 35,000 at m/z 200; AGC target, 3 e 6; Level 2 maximum IT, 120 ms; MS2 activation type, HCD; peptide fragmentation, nitrogen; isolation window, 2.0 Th; normalized collision energy, 28 eV. Four proteins were randomly selected from the global proteomics analysis (re-labeled peptide GISNEGQNASIK^* as the reference standard mixture), including phosphoglucomutase (TFTTQETITNAATSAK), glucose-6-phosphate isomerase (LGAENFVFFHPR), xylose isomerase (YFGNLMDAGR, LSICGEESFGTGSDHIR), and malate dehydrogenase (IFGVTTLDIVR, IQDAGTEVVK). Skyline 3.5 was used to generate an initial PRM transition pair list for the four candidate DEPs.

Results

SDS-PAGE Analysis and Protein Profiling

BCA results indicated protein concentrations of 15.60, 20.94, and 21.55 μg/μL in fed isopods and 19.06, 25.07, and 22.03 μg/μL in unfed isopods. TMT analysis indicated 3,492 queries in the 133,946 spectra. Among them, a total of 888 unique proteins were identified across 1,510 peptides (Figure 1). There were 139 proteins between 0 and 20 kDa, and 282, 241, 92, 45, and 89 proteins of 20–40, 40–60, 60–80, and 80–100 kDa, respectively. Eighty-nine proteins had a mass of over 100 kDa (Figure 2).

FIGURE 1

Figure 1. Total numbers of tandem mass spectrometry spectra (MS/MS spectrum), peptide-spectrum matches (PSMs), peptides, unique peptides, and proteins in Tachaea chinensis.

FIGURE 2

Figure 2. Distribution of identified proteins with different molecular weights (the weight and pI values are for the target proteins from the NCBI_Peracarida_91190_20190313.fasta database).

Differentially Expressed Protein Analysis

Using a 1.2-fold increase or decrease in protein expression as a benchmark for a physiologically significant change, 129 differentially expressed proteins (DEPs, p < 0.05) were identified between the control and parasitism groups (Figure 3). Within these DEPs, 37 were upregulated and 92 were downregulated (Table S1, Figure S2) before parasitism.

FIGURE 3

Figure 3. Volcano plot of all proteins identified in the proteome database, for the fed and unfed parasite groups. Yellow points: differentially expressed proteins that were significantly downregulated after parasitism (fold change > 1.2; p < 0.05). Green points: differentially expressed proteins that were significantly upregulated after parasitism (fold change < 0.8; p < 0.05).

According to GO enrichment analysis, 1,256, 206, and 274 proteins were enriched for the categories biological process, cell component, and molecular function, respectively, and 178, 71, and 83 of these were statistically significant (Figure S2). The top 10 significant pathways are shown in Figure 4. Fifty-seven KEGG proteins were enriched, and five were statistically significant. Their functions included ECM-receptor interaction, pentose phosphate pathways, protein processing in the endoplasmic reticulum, glutathione metabolism, and phototransduction-fly, with 3, 3, 7, 4, and 3 proteins mapped to them, respectively (Figure 5). According to the protein–protein interaction network analysis, the three major clusters of interaction were protein processing in endoplasmic reticulum, ribosomes, and glutathione metabolism (Figure 6).

FIGURE 4

Figure 4. GO functional classification of the identified proteins of the top 10 pathways identified in the proteome database, for the fed and unfed parasite groups, in the biological process, cell component, and molecular function categories.

FIGURE 5

Figure 5. KEGG annotation of the top five pathways from the proteome database of the fed and unfed parasite groups.

FIGURE 6

Figure 6. Protein-protein interaction network analysis of the three major network groups.

Candidate Proteins Related to Energy and Digestive Metabolism

Some differentially expressed proteins related to energy metabolism were significantly upregulated in unfed T. chinensis, including calcium-transporting ATPase sarcoplasmic/endoplasmic reticulum type-like proteins, NADH dehydrogenase [ubiquinone] flavoprotein 1, mitochondrial-like Na⁺/K⁺-ATPase, partial sodium potassium ATPase alpha subunit, partial smooth endoplasmic reticulum calcium ATPase, NADH dehydrogenase subunit 1 (mitochondrion), and cyclic AMP-responsive element-binding protein 1 (Table 1). Among the 92 significantly upregulated differentially expressed proteins, 20 were related to digestive metabolism, including 7, 11, and 2 in carbohydrate, protein, and lipid metabolism, respectively (Table 1).

TABLE 1

Table 1. Differentially expressed proteins involved in energy and digestive metabolism, comparing the fed and unfed Tachaea chinensis.

Candidate Proteins Related to Blood Sucking

Three hemocyanin proteins showed the greatest difference in expression between fed and unfed isopods (Table S1). Differentially expressed proteins related to blood sucking, including inorganic pyrophosphatase-like isoform X1, alpha-2-macroglobulin, bestrophin homolog, calreticulin, barrier-to-autointegration factor putative salivary alkaline phosphatase, inorganic pyrophosphatase-like isoform X1, neurocalcin homolog isoform X3, and insulin receptor substrate 2-B were significantly upregulated after parasitism (Table 2). Four heat shock proteins and two stress-activated proteins were also upregulated.

TABLE 2

Table 2. Differentially expressed proteins involved in blood sucking, comparing the fed and unfed Tachaea chinensis.

PRM Results

To validate the proteomic data, three proteins that were upregulated after parasitism (phosphoglucomutase, glucose-6-phosphate isomerase, and xylose isomerase) and one that was downregulated after parasitism (malate dehydrogenase) were selected for PRM analysis. In both the fed and unfed groups, the validated proteins showed expression trends similar to the proteomic expression trends, suggesting that the proteomic data were reliable (Figure 7).

FIGURE 7

Figure 7. Verification of the expression profiles of randomly selected significantly changed proteins comparing the fed and unfed T. chinensis groups, as determined by parallel reaction monitoring (PRM) and tandem mass tag (TMT) analysis. The data on the relative expression of proteins, for use in the TMT analysis, were obtained from Tables 1, 2. The data on the relative expression of proteins for use in the PRM analysis were calculated relative to the expression of the relabeled peptide GISNEGQNASIK*.

Discussion

Hematophagous parasites use antihemostatic, anti-inflammatory, and immunomodulatory compounds to facilitate blood meal acquisition during feeding (Shu et al., 2010; Hajdusek et al., 2013; Eliane et al., 2017; Santiago et al., 2017; Jalovecka et al., 2018). These promote the survival and establishment of pathogens in the host (Boonsriroj et al., 2015; Eliane et al., 2017). Identification and characterization of such compounds can elucidate molecular mechanisms of interaction among parasites, pathogens, and hosts, revealing new vaccine targets (Mcgowan et al., 1982; Slovák et al., 2013; Oleaga et al., 2015; Yu et al., 2015; Eliane et al., 2017; Tirloni et al., 2017). Although many researchers have assumed that isopod parasites are hematophagous parasites and feed mainly on fish blood or crustacean hemolymph, no prior studies have demonstrated this directly. Moreover, in contrast to ticks and mosquitoes, there is a lack of comprehensive genomic resources for isopod parasites to investigate infection and response mechanisms. We performed, to our knowledge, the first comparative proteomic analysis between fed and unfed T. chinensis to gain a better understanding of its feeding and antigen-related proteins. However, because very few protein sequences for isopods are available in the NCBI protein database and all proteins identified were based on homologous proteins, we were able to identify only 888 proteins. Moreover, to avoid potential contamination of the host protein samples and to screen out proteins associated with parasitic processing, we also searched the NCBI protein database for proteins of Peracarida and blood-sucking insects (ticks and mosquitos); we found that proteins from Armadillidium vulgare (442, 51.76%), Hyalella azteca (227, 26.58%), and Hirondellea gigas (143, 16.74%) were similar to those present in our samples (Table S1).

Comparative proteomic analysis of unfed and fed T. chinensis indicated 129 proteins expressed significantly, with 37 upregulated and 92 downregulated before parasitism. Among the 37 upregulated proteins, seven were identified as contributing to ATP synthesis (Table 1). These results suggest that high energy metabolic proteins are expressed in pre-parasitism isopods while searching for hosts. Meanwhile, previous research has shown that genes related to ATP generation in host shrimp P. sinensis are upregulated under parasitism, which might be induced by a rapid depletion of ATP content in the parasite itself after parasitism, during its development and reproduction (Li et al., 2018). A similar profile has been demonstrated in fleas, where energy costs for feeding are lower than those in unfed individuals (Sarfati, 2005). Moreover, in crustaceans, because proteins are major contributors to hemolymph density (Simonetta et al., 2011), in their parasites, the major digestive enzymes are thought to be proteases. Eleven proteolytic enzymes were identified among the 129 DEPs, and all were upregulated after parasitism (Table 1). Seven carbohydrate metabolism-related proteins were also upregulated, suggesting that hemolymph glucose may also be an important nutrient for T. chinensis. Trehalose, which regulates energy metabolism and glucose generation via trehalose catabolism, was also upregulated after parasitism, suggesting that access to host trehalose is an important feature of isopod parasites.

Many researchers have assumed that isopod parasites can feed on fish blood or crustacean hemolymph, comparable to an “aquatic mosquito” (Wilson, 2008; Manship et al., 2012; Smit et al., 2014; Nagler et al., 2017). We previously found that hemocyanin in P. sinensis is significantly downregulated following parasitism by T. chinensis (Li et al., 2019). The present results indicate that hemocyanin in T. chinensis was significantly upregulated after parasitism, suggesting a role in obtaining hemolymph from P. sinensis. In hematophagous arthropods such as ticks, mosquitos, fleas, and flies, hematophagy evolved independently over millions of years, leading to various morphological adaptations and diverse strategies to overcome the barriers imposed by hosts (Ribeiro, 1995; Garcia et al., 2015). However, one example of convergent evolution and adaptation for hematophagy is the development of highly functional salivary glands, which produce a large diversity of anticlotting, antiplatelet, vasodilatory, and pathogen-transporting substances (Kim et al., 2016; Tirloni et al., 2016). Although it is known that isopod parasite morphology has been strongly modified for parasitism over 168 million years (Nagler et al., 2017), the type and composition of blood-sucking-related substances remain to be understood.

During blood sucking, the salivary gland of hematophagous parasites produces a cocktail of anti-hemostatic, anti-inflammatory, and immunomodulatory molecules that facilitate blood-meal acquisition. Saravanan et al. (2003) found that α₂-macroglobulin, which they consider to have an anticoagulant function, is highly expressed in tick salivary glands. In the present study, α₂-macroglobulin was significantly upregulated in isopods that had recently fed, which indicates that T. chinensis may also secrete anticoagulatory molecules during feeding. Moreover, Saurabh et al. (2011) found that fish macroglobulin was downregulated during Argulus siamensis parasitism. Similarly, we previously found that shrimp macroglobulin was downregulated during T. chinensis parasitism. Moreover, calreticulin, which functions in hematophagy in ticks through host immunosuppression or antihemostasis (Parizi et al., 2009), was also significantly upregulated after parasitism in the present study. Furthermore, calcineurin subunit B type 2 isoform X1, putative salivary alkaline phosphatase, insulin receptor substrate 2-B, PI-actitoxin-Aeq3a, and neurocalcin homolog isoform X3, which are similar to substances in the salivary glands of hematophagous parasites, were also upregulated in T. chinensis during feeding (Kim et al., 2015). However, the interaction of those proteins between parasites and hosts and the mechanisms of paralyzation and inhibition of inflammation remain unclear and require further study.

Anticlotting-related proteins are considered potential targets for the development of new drugs and vaccines to control and prevent ectoparasites such as ticks (Oleaga et al., 2017), mosquitos (Popova-Butler and Donald, 2009; Tetreau et al., 2012), Schistosoma (Delcroix et al., 2006), and sea louse (Carpio et al., 2011, 2016). However, the hosts of these parasites are vertebrates, in which the composition of the coagulation system differs from that of invertebrates. In shrimps, the hemolymph clotting system comprises transglutaminase and clotting proteins and plays an important role in the innate immune response and prevention of blood loss during injury and wound healing (Maningas et al., 2013). Isopod parasites may inject anticoagulants or other compounds directly into the blood to obtain a blood meal (Nair and Nair, 1983; Manship et al., 2012), but the constituents of these compounds are still unclear. Here, we identified two protein disulfide-isomerase (PDI) proteins that participate in parasite–host cell interactions. PDI-specific antibodies may constitute part of the mucosal antibody repertoire, which is possibly involved in defense against parasites. The potential roles of those proteins for isopod parasite control will be investigated in our further study.

In conclusion, our results represent the first comparative proteomic study to detect key proteins from the whole body of fed and unfed T. chinensis isopod parasites. Owing to the lack of genomic data, we were able to identify only 888 proteins, and all were found in both stages. Differentially expressed proteins related to energy metabolism were upregulated in unfed individuals, particularly those involved in ATP generation. This may indicate that T. chinensis may require considerable energy during its search for a host. We also found that similar to other hematophagous parasites, T. chinensis secretes antihemostatic, anti-inflammatory, and immunomodulatory molecules to facilitate blood meal acquisition. Our study provides valuable empirical data that will support future molecular research on isopod parasitization of crustaceans.

Data Availability Statement

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

YL and QC designed the experiments. XinL and ZH generated biological samples. XinL, ZH, and WX performed the experiments. YL, XinL, and QC analyzed data. YL and WX performed statistic data analysis. YL and XiaL contributed reagents, materials and analysis tools. YL, XinL, ZH, and WX wrote the paper. All authors read and approved the final manuscript.

Funding

This work was supported by the earmarked fund for Modern Agro-industry Technology Research System (No. CARS-48), the National Natural Science Foundation of China (No. 31502190), and the Talent Introduction Program of Shenyang Agricultural University (Nos. 880416005 and 880417024).

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.

Acknowledgments

We thank Panjin Guanghe Fisheries Co., Ltd. for providing the shrimps.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2019.00350/full#supplementary-material

Figure S1. Separation of proteins by SDS-PAGE from six unfed (A1, A2, and A3) and fed (B1, B2, and B3) Tachaea chinensis.

Figure S2. Hierarchical clustering analysis of differentially expressed proteins in unfed (A1, A2, and A3) and fed (B1, B2, and B3) Tachaea chinensis.

Table S1. Lists of proteins indented from Tachaea chinensis.

References

Boonsriroj, H., Manalo, D. L., Kimitsuki, K., Shimatsu, T., Shiwa, N., Shinozaki, H., et al. (2015). A pathological study of the salivary glands of rabid dogs in the Philippines. J. Vet. Med. Sci. 78, 35–42. doi: 10.1292/jvms.15-0308

PubMed Abstract | CrossRef Full Text | Google Scholar

Carpio, Y., Basabe, L., Acosta, J., Rodríguez, A., Mendoza, A., Lisperger, A., et al. (2011). Novel gene isolated from Caligus rogercresseyi: a promising target for vaccine development against sea lice. Vaccine 29, 2810–2820. doi: 10.1016/j.vaccine.2011.01.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Carpio, Y., Velazquez, J., Leal, Y., Herrera, N., García, C., Acosta, J., et al. (2016). Sea lice vaccine development based on an immunogenic peptide derived from the ribosomal protein P0. Fish Shellfish Immunol. 53, 64–64. doi: 10.1016/j.fsi.2016.03.061

CrossRef Full Text | Google Scholar

Delcroix, M., Sajid, M., Caffrey, C. R., Lim, K.-C., Dvorak, J., Hsieh, I., et al. (2006). A multienzyme network functions in intestinal protein digestion by a platyhelminth parasite. J. Biol. Chem. 281, 39316–39329. doi: 10.1074/jbc.M607128200

PubMed Abstract | CrossRef Full Text | Google Scholar

Eliane, 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

CrossRef Full Text | Google Scholar

Garcia, G. R., Maruyama, S. R., Malardo, T., Zangirolamo, A. F., and Gardinassi, L. G. (2015). The biology of hematophagous arthropods addressed by molecular high-throughput approaches. Austin J. Trop. Med. Hygiene 1:1004.

Google Scholar

Hajdusek, O., Síma, R., Ayllón, N., Jalovecká, M., Perner, J., De La Fuente, J., et al. (2013). Interaction of the tick immune system with transmitted pathogens. Front. Cell. Infect. Microbiol. 3:26. doi: 10.3389/fcimb.2013.00026

PubMed Abstract | CrossRef Full Text | Google Scholar

Hua, C. J., Li, W. X., Zhang, D., Zou, H., Li, M., Jakovlić, I., et al. (2018). Basal position of two new complete mitochondrial genomes of parasitic Cymothoida (Crustacea: Isopoda) challenges the monophyly of the suborder and phylogeny of the entire order. Parasit. Vectors 11:628. doi: 10.1186/s13071-018-3162-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jalovecka, M., Hajdusek, O., Sojka, D., Kopacek, P., and Malandrin, L. (2018). The complexity of piroplasms life cycles. Front. Cell. Infect. Microbiol. 8:248. doi: 10.3389/fcimb.2018.00248

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, T. K., Tirloni, L., Pinto, A. F., Moresco, J., Da Silva Vaz Junior, I., and Mulenga, A. (2016). Ixodes scapularis tick saliva proteins sequentially secreted every 24 h during blood feeding. PLoS Negl. Trop. Dis. 10:e0004323. doi: 10.1371/journal.pntd.0004323

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y. H., Slam, M. S., and You, M. J. (2015). Proteomic screening of atigenic proteins from the hard tick, Haemaphysalis longicornis (Acari: Ixodidae). Korean J. Parasitol. 53:85. doi: 10.3347/kjp.2015.53.1.85

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y. D., Li, X., Xu, W. B., Han, Z. B., Zhao, Y. Y., Dong, J., et al. (2019). Comparative iTRAQ-based quantitative proteomic analysis of the Chinese grass shrimp (Palaemonetes sinensis) infected with the isopod parasite Tachaea chinensis. Parasit. Vectors 12:415. doi: 10.1186/s13071-019-3675-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y. D., Xu, W. B., Li, X., Jiang, H. B., She, Q. X., Han, Z. B., et al. (2018). Comparative transcriptome analysis of Chinese grass shrimp (Palaemonetes sinensis) infected with isopod parasite Tachaea chinensis. Fish Shellfish Immunol. 82, 153–161. doi: 10.1016/j.fsi.2018.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, X., Zou, W., Peng, H., Shu, Q., Chen, W., Xie, J., et al. (2007). New records of two parasitic isopods in Peal River. J. Fish. Sci. China 14, 1023–1026.

Google Scholar

Maningas, M. B. B., Hidehiro, K., and Ikuo, H. (2013). Molecular mechanisms of the shrimp clotting system. Fish Shellfish Immunol. 34, 968–972. doi: 10.1016/j.fsi.2012.09.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Manship, B. M., Walker, A. J., Jones, L. A., and Davies, A. J. (2012). Blood feeding in juvenile Paragnathia formica (Isopoda: Gnathiidae): biochemical characterization of trypsin inhibitors, detection of anticoagulants, and molecular identification of fish hosts. Parasitology 139, 744–754. doi: 10.1017/S0031182011002320

PubMed Abstract | CrossRef Full Text | Google Scholar

Mcgowan, M. J., Homer, J. T., and Barker, R. W. (1982). Polyacrylamide gel electrophoresis of unfed, partially fed and gut tissues of Amblyomma americanum (L.). Florida Entomol. 65, 286–288. doi: 10.2307/3494293

CrossRef Full Text | Google Scholar

Nagler, C., HyŽný, M., and Haug, J. T. (2017). 168 million years old “marine lice” and the evolution of parasitism within isopods. BMC Evol. Biol. 17:76. doi: 10.1186/s12862-017-0915-1

CrossRef Full Text | Google Scholar

Nair, G. A., and Nair, N. B. (1983). Effect of infestation with the isopod, Alitropus typus M. Edwards (Crustacea: Flabellifera: Aegidae) on the haematological parameters of the host fish, Channa striatus (Bloch). Aquaculture 30, 11–19. doi: 10.1016/0044-8486(83)90147-3

CrossRef Full Text | Google Scholar

Oleaga, A., Obolo-Mvoulouga, P., Manzano-Román, R., and Pérez-Sánchez, R. (2015). Midgut proteome of an argasid tick, Ornithodoros erraticus: a comparison between unfed and engorged females. Parasit. Vectors 8:525. doi: 10.1186/s13071-015-1148-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Oleaga, A., Obolo-Mvoulouga, P., Manzano-Román, R., and Pérez-Sánchez, R. (2017). A proteomic insight into the midgut proteome of Ornithodoros moubata females reveals novel information on blood digestion in argasid ticks. Parasit. Vectors 10:366. doi: 10.1186/s13071-017-2300-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Parizi, L. F., Rech, H., Ferreira, C. A., Imamura, S., Ohashi, K., Onuma, M., et al. (2009). Comparative immunogenicity of Haemaphysalis longicornis and Rhipicephalus (Boophilus) microplus calreticulins. Vet. Parasitol. 164, 282–290. doi: 10.1016/j.vetpar.2009.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Poore, G. C., and Bruce, N. L. (2012). Global diversity of marine isopods (except Asellota and crustacean symbionts). PLoS ONE 7:e43529. doi: 10.1371/annotation/3260cd00-89cf-4e08-ab25-07e0be598ab4

PubMed Abstract | CrossRef Full Text | Google Scholar

Popova-Butler, A., and Donald, H. D. (2009). Proteomic analysis of the mosquito Aedes aegypti midgut brush border membrane vesicles. J. Insect Physiol. 55, 264–272. doi: 10.1016/j.jinsphys.2008.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribeiro, J. M. (1995). Blood-feeding arthropods: live syringes or invertebrate pharmacologists? Infect. Agents Dis. 4, 143–152.

PubMed Abstract | Google Scholar

Santiago, P. B., Araújo, C. N. D., Motta, F. N., Praça, Y. R., Charneau, S., Bastos, I. M. D., et al. (2017). Proteases of haematophagous arthropod vectors are involved in blood-feeding, yolk formation and immunity - a review. Parasit. Vectors 10:79. doi: 10.1186/s13071-017-2005-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Saravanan, T., Weise, C., Sojka, D., and Kopáček, P. (2003). Molecular cloning, structure and bait region splice variants of α-macroglobulin from the soft tick. Insect Biochem. Mol. Biol. 33, 841–851. doi: 10.1016/S0965-1748(03)00083-3

CrossRef Full Text | Google Scholar

Sarfati, M. (2005). Energy costs of blood digestion in a host-specific haematophagous parasite. J. Exp. Biol. 208, 2489–2496. doi: 10.1242/jeb.01676

PubMed Abstract | CrossRef Full Text | Google Scholar

Saurabh, S., Mohanty, B. R., and Sahoo, P. K. (2011). Expression of immune-related genes in rohu Labeo rohita (Hamilton) by experimental freshwater lice Argulus siamensis (Wilson) infection. Vet. Parasitol. 175, 119–128. doi: 10.1016/j.vetpar.2010.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, Q. T., Glanfield, A., Gobert, G. N., and Jones, M. K. (2010). Heme and blood-feeding parasites: friends or foes? Parasit. Vectors 3:108. doi: 10.1186/1756-3305-3-108

CrossRef Full Text | Google Scholar

Simonetta, L., Martinis, M., and Ferrero, E. A. (2011). Ecological relevance of hemolymph total protein concentration in seven unrelated crustacean species from different habitats measured predictively by a density-salinity refractometer. J. Mar. Biol. 2011, 1–7. doi: 10.1155/2011/153654

CrossRef Full Text | Google Scholar

Slovák, M., Hajnická, V., Labuda, M., and Fuchsberger, N. (2013). Comparison of the protein profiles of salivary gland extracts derived from three species of unfed and partially fed ixodid ticks analysed by SDS-PAGE. Folia Parasitol. 47, 67–71. doi: 10.14411/fp.2000.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Smit, N. J., Bruce, N. L., and Hadfield, K. A. (2014). Global diversity of fish parasitic isopod crustaceans of the family Cymothoidae. Int. J. Parasitol. Parasites Wildl. 3, 188–197. doi: 10.1016/j.ijppaw.2014.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tetreau, G., Bayyareddy, K., Jones, C. M., Stalinski, R., Riaz, M. A., Paris, M., et al. (2012). Larval midgut modifications associated with Btiresistance in the yellow fever mosquito using proteomic and transcriptomic approaches. BMC Genomics 13:248. doi: 10.1186/1471-2164-13-248

PubMed Abstract | CrossRef Full Text | Google Scholar

Tirloni, L., Kim, T. K., Coutinho, M. L., Ali, A., Seixas, A., Termignoni, C., et al. (2016). The putative role of Rhipicephalus microplus salivary serpins in the tick-host relationship. Insect Biochem. Mol. Biol. 71, 12–28. doi: 10.1016/j.ibmb.2016.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tirloni, L., Kim, T. K., Pinto, A. F. M., Yates, J. R., Da Silva Vaz Junior, I., and Mulenga, A. (2017). Tick-host range adaptation: changes in protein profiles in unfed adult Ixodes scapularis and Amblyomma americanum saliva stimulated to feed on different hosts. Front. Cell. Infect. Microbiol. 7:517. doi: 10.3389/fcimb.2017.00517

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, J. D., and Boyko, C. B. (2012). The global diversity of parasitic isopods associated with crustacean hosts (Isopoda: Bopyroidea and Cryptoniscoidea). PLoS ONE 7:e35350. doi: 10.1371/journal.pone.0035350

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, G. D. F. (2008). Global diversity of Isopod crustaceans (Crustacea; Isopoda) in freshwater. Hydrobiologia 595, 231–240. doi: 10.1007/s10750-007-9019-z

CrossRef Full Text | Google Scholar

Xiao, C. F., Gao, J. H., Zhang, Y. X., Wang, Z., and Zhang, D. H. (2019). Quantitative proteomics of potato leaves infected with Phytophthora infestans provides insights into coordinated and altered protein expression during early and late disease stages. Int. J. Mol. Sci. 20:E136. doi: 10.3390/ijms20010136

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X., Gong, H., Zhou, Y., Zhang, H., Cao, J., and Zhou, J. (2015). Differential sialotranscriptomes of unfed and fed Rhipicephalus haemaphysaloides, with particular regard to differentially expressed genes of cysteine proteases. Parasit. Vectors 8:597. doi: 10.1186/s13071-015-1213-7

PubMed Abstract | CrossRef Full Text | Google Scholar