Integrative Proteomic Analysis of Digestive Tract ... - ACS Publications

Jul 19, 2019 - different biological features, including their resistance to environmental ... variety of food types, such as submerged and emergent ...
3 downloads 0 Views 1024KB Size
Subscriber access provided by BUFFALO STATE

Article

Integrative proteomic analysis of digestive tract glycosidases from the invasive golden apple snail, Pomacea canaliculata Sophia Escobar-Correas, Omar Mendoza-Porras, Federico Dellagnola, Michelle L Colgrave, and Israel A. Vega J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/acs.jproteome.9b00282 • Publication Date (Web): 19 Jul 2019 Downloaded from pubs.acs.org on July 22, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

INTEGRATIVE PROTEOMIC ANALYSIS OF DIGESTIVE TRACT GLYCOSIDASES FROM THE INVASIVE GOLDEN APPLE SNAIL, Pomacea canaliculata Sophia Escobar-Correasa,b, Omar Mendoza-Porrasc, Federico A. Dellagnolaa,b,d, Michelle L. Colgravec*, Israel A. Vegaa,b,d*

a IHEM,

CONICET, Universidad Nacional de Cuyo, Mendoza, Argentina.

b Universidad

Nacional de Cuyo, Facultad de Ciencias Médicas, Instituto de Fisiología,

Mendoza, Argentina. c CSIRO,

Agriculture & Food, 306 Carmody Road, St. Lucia, Queensland 4067, Australia

d Universidad

Nacional de Cuyo, Facultad de Ciencias Exactas y Naturales, Departamento de

Biología, Mendoza, Argentina.

Authors for correspondence: Michelle Colgrave E-mail: [email protected] Israel A. Vega E-mail: [email protected]

ACS Paragon Plus Environment

1

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 38

ABSTRACT The freshwater snail Pomacea canaliculata, an invasive species of global significance, possesses a well-developed digestive system and diverse feeding mechanisms enabling the intake of a wide variety of food. The identification of glycosidases in adult snails would increase the understanding of their digestive physiology and potentially generate new opportunities to eradicate and/or control this invasive species. In this study, liquid chromatography coupled to tandem mass spectrometry was applied to define the occurrence, diversity and origin of glycoside hydrolases along the digestive tract of P. canaliculata.. A range of cellulases, hemicellulases, amylases, maltases, fucosidases and galactosidases were identified across the digestive tract. The digestive gland and the contents of the crop and style sac yield a higher diversity of glycosidase-derived peptides. Subsequently, peptides derived from 81 glycosidases (46 proteins from the public database and 35 uniquely from the transcriptome database) that were distributed amongst 13 glycoside hydrolase families were selected and quantified using multiple reaction monitoring mass spectrometry. This study showed a high glycosidase abundance and diversity in the gut contents of P. canaliculata which participate in extracellular digestion of complex dietary carbohydrates. Salivary and digestive glands were the main tissues involved in their synthesis and secretion.

KEYWORDS Ampullariids, Liquid chromatography-mass spectrometry (LC-MS), Digestion, Cellulases, Hemicellulases, Amylases, Maltases, Fucosidases, Galactosidases.

ACS Paragon Plus Environment

2

Page 3 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

1. INTRODUCTION Apple snails (Caenogastropoda, Ampullariidae) are naturally distributed in humid tropical and subtropical ecosystems of Africa, America, and Asia1. The International Union for Conservative of Nature and the Invasive Species Specialist Group have included three American apple snails, Pomacea canaliculata, Pomacea maculata, and Marisa cornuarietis, in the list of the top “100 world’s worst invasive alien species”2. Pomacea species have also been introduced anthropogenically in South East Asia, China, Japan, USA, and Europe3-5 where they alter ecosystem dynamics4 affecting the abundance and richness of native species6, and producing economic losses in agricultural industries, in particular of rice7. The invasiveness of P. canaliculata has been associated with different biological features, including their resistance to environmental stressors8, pollutants9, 10, tissue regeneration in adult individuals11, and suitability to reproduce in different environments12,

13.

Remarkably, this

polyphagous species has diverse feeding mechanisms1, 14 including the intake of a wide variety of food types, such as submerged and emergent macrophytes, biofilms, periphyton and animal carrion15, 16. Additionally, it has a well-developed digestive system specialised for the digestion of complex polymers1, 17, 18. However, knowledge of the mechanisms of carbohydrate digestion by P. canaliculata remains incomplete. Cellulose-digestion enzymes have been associated with the gut microbiome and stomach of this snail19, 20. Endogenous cellulases from glycoside hydrolase families (GHFs) 9, 10 and 45 have been previously reported from Ampullaria crossean in Xiamen, China21-23, and from P. canaliculata collected in Bangkok, Thailand19; A. crossean should be a misidentification of P. canaliculata1. Also, various glycosidase activities (endo-1,4-β-D-xylanase, α- and β-mannosidase, α-fucosidase and cellulase) have been found in whole extracts of P. canaliculata24-26, but there is

ACS Paragon Plus Environment

3

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 38

no available information about the spatial and temporal framework where the digestion of complex dietary polymers occurs. The study of the set of digestive enzymes responsible for supplying macronutrients to this snail, may allow the identification of molecular targets for control strategies aiming to stop the spread of this invasive pest. In this study the occurrence, diversity and origin of glycosidases detected along the digestive tract of P. canaliculata were resolved by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). Subsequently, the abundance of each enzyme was followed by quantitative proteomics using multiple reaction monitoring mass spectrometry (MRM-MS).

2. MATERIALS AND METHODS

2.1. Snail and culture conditions Adult snails (shell length = 33–46 mm) from a cultured strain of P. canaliculata were used. Room temperature was regulated (24-26 °C) and artificial lighting was provided 14 h per day. Room relative humidity was 80%. The animals were maintained in aquaria containing 6 L of tap water. Aquarium water was changed thrice per week. Unless otherwise indicated, animals were fed ad libitum with lettuce from Monday through Friday. This diet was supplemented with fish food pellets with a 40% of total protein (Peishe Car Shulet, Argentina) on Thursday and with toilet paper on Friday. The original stock of snails was collected at the Rosedal Lake (Palermo, Buenos Aires, Argentina; neotype locality). Animals from the original population and from the laboratory strain have been deposited at the collection of Museo Argentino de Ciencias Naturales (Buenos Aires, Argentina; lots MACN-In 35707 and MACN-In 36046, respectively). Sacrifice procedures and

ACS Paragon Plus Environment

4

Page 5 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

sampling were carried out in agreement with the policies of the Institutional Animal Care and Use Committee of the Facultad de Ciencias Médicas of the Universidad Nacional de Cuyo (Approval Protocol N° 55/2015).

2.2. Animal acclimation, sampling, and protein extraction A first experiment was aimed at determining the glycoside hydrolases involved in carbohydrate digestion at different sites along the digestive tract of P. canaliculata. Five sampling sites were collected from four replicate snails: (a) the medial portion of the oesophagus i.e. the crop (Cr), which retains food and salivary gland juices that are originated in the buccal cavity; (b) the final portion of the stomach, i.e. the style sac (SS); (c) the coiled gut (CG), which receives contents of the small intestine; (d) both salivary glands (SG); and (e) the digestive gland (DG). The animals were acclimated to a feeding scheme, which has been described elsewhere18. Briefly, animals were fed ad libitum and then fasted for 24 h. After that, each animal was isolated in a beaker containing 70 mL water with three fish food pellets. Once the first pellet was swallowed, fresh lettuce was added into the beaker. Animals consumed all the food that was given. After 90 min, each animal was immersed in an ice-water bath for 10 min to minimise pain and then the shell was cracked and removed. Autostatic forceps were fixed on the posterior oesophagus to prevent the passage of contents between the crop and the stomach during sampling. The crop and the style sac contents were collected using a 1 mL syringe with an 18-gauge needle and were thoroughly dispersed in Protein Extraction Buffer (PEB = 100 mM Tris-HCl pH 7.4, 10 mM NaCl, 0.25%, Triton X-100). Samples from salivary glands, digestive gland, and coiled gut were obtained and homogenised on ice (five cycles of 15 s) with an Ultraturrax homogeniser in 500 μL of PEB. All samples were centrifuged at 12,350 xg for 15 min at 4°C; the pellets were discarded, and the

ACS Paragon Plus Environment

5

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 38

protein concentration of the supernatant was quantified according to Bradford protein assay. Aliquots of 330 μg protein were fixed in 70% ethanol (for shipment to Australia) and stored at −80°C until use.

2.3. Protein preparation and tryptic digestion for mass spectrometry After removing ethanol, the proteins were lyophilised and reconstituted in 165 μL urea buffer (8M Urea, 0.1M Tris-HCl pH 8.5). The samples were vortexed (30 s, five times), and sonicated (2 min; once); these procedures were repeated twice. Proteins were quantified by Bradford protein assay using the manufacturers’ instructions yielding protein concentrations: salivary gland = 1.11.3 μg of protein/μL; crop content = 3.2-9.2 μg of protein/μL; digestive gland = 3.3-9.6 μg of protein/μL; style sac content = 1.6-10.3 μg of protein/μL; coiled gut = 6.4-12.2 μg of protein/μL. Soluble proteins (100 μg) from each sample extract were transferred to a 10 kDa MWCO filter unit (Millipore) and centrifuged at 20,800 xg for 10 min at room temperature. Three technical replicates were prepared from each sample. The retained proteins (>10 kDa) were subsequently reduced, alkylated and digested using a modified version of the filter-assisted sample preparation (FASP) workflow27 that is compatible with reagents (detergents, chaotropes) and downstream analysis by LC-MS/MS. The filters were washed three times with urea buffer and then centrifuged (20,800 xg, 10 min). Cysteine residues were reduced with 100 mM dithiothreitol (DTT) for 60 min at room temperature. Excess DTT was removed by centrifugation (20,800 xg, 10 min). Free sulfhydryl groups were alkylated by incubation with 50 mM iodoacetamide for 20 min at room temperature followed by centrifugation (20,800 xg, 10 min). After reduction and alkylation, each filter unit was washed twice with 200 μL of urea buffer and then twice with 200 μL of 50 mM ammonium

ACS Paragon Plus Environment

6

Page 7 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

bicarbonate (AmBic, pH 8.5); each washed was followed by centrifugation (20,800 xg, 10 min). Proteins were digested overnight at 37°C with sequencing grade trypsin (V5111, Promega; protein/trypsin ratio = 50/1) in AmBic (pH 8.5). After digestion, the resulting tryptic peptides (99% confidence, with one or more unique peptides (>95% confidence) using the false discovery rate (FDR) analysis provided within the ProteinPilot software.

ACS Paragon Plus Environment

8

Page 9 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Peptide summaries (Supporting Table 1) generated by ProteinPilot were analysed. Peptides with low peak intensity ( DG > SG > CG). The fucosidases were mainly found in tissues (SG, DG and CG, range 13-15), with the style sac content yielding only two representatives and none

ACS Paragon Plus Environment

10

Page 11 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

detected in the crop. The number of galactosidases (range 11-22) were variable between contents and tissues (Cr = DG > SS > CG > SG). A comparison of the tissue-specific expression was undertaken examining unique and shared peptides of the identified glycosidases (Figure 3) in order to avoid issues with protein inference arising from the use of a non-redundant protein database. The crop content showed approximately four times more unique peptides than the style sac content (410 versus 102). The digestive gland showed the highest number of unique peptides compared to salivary gland and coiled gut (165 > 95 > 30). Thirty-nine peptides were shared among digestive contents and tissues of the digestive tract of this snail. Crop and style sac contents had the highest number of peptides uniquely detected (438 peptides) while coiled gut, salivary and digestive glands shared 96 peptides.

3.2. Peptide quantification and glycoside hydrolase families A total of 250 peptides belonging to 81 glycosidases (46 from the public databases and 35 uniquely from the transcriptome database) were selected (Supporting Table 3) for LC-MRM-MS relative quantitation. Figure 4 shows an analysis of the relative percent of each peptide in relation to glycosidase type and their distribution along the digestive tract. Peptides of cellulases and amylases (Figure 4A, C) showed higher values for both contents (crop and style sac) as opposed to the digestive gland. Also, peptides from cellulases and amylases showed lower or null percent values in samples taken from salivary glands and coiled gut. Some peptides of these enzymes (cellulases Pc125214 and A0A2T7PT81; amylases A0A2T7NPM3, A0A2T7PLX3 and Pc109575) showed an idiosyncratic behavior for which lower values in crop and higher values in style sac were found (Supporting Figure 1).

ACS Paragon Plus Environment

11

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 38

The distribution of peptides from hemicellulases (Figure 4B), maltases (Figure 4D) and galactosidases (Figure 4F) were more variable across the different samples taken from the P. canaliculata digestive tract. Unexpected findings were the different profiles (Supporting Figure 1) found for peptides of maltases A0A2T7PY79 (VTASSLTK and QNYPVTLDTVR) and A0A2T7NPM8 (QGGVDLPTADSCQK). On the other hand, fucosidases were exclusively found in tissues (Figure 4E), with a significantly higher proportion in the digestive gland compared to the salivary gland and coiled gut. An analysis of amino acid sequence similarities in the carbohydrate-active enzymes database (CAZy)37 showed 13 glycoside hydrolase families. A total of thirty cellulases in three families: sixteen GHF9 (Figure 5A), thirteen GHF10 (Figure 5B) and one GHF45 (Figure 5C), were highly expressed in crop and style sac contents. Nineteen hemicellulases were distributed into six families of glycoside hydrolases: two GHF2 (Figure 5D), five GHF3 (Figure 5E), four GHF5 (Figure 5F), three GHF10 (Figure 5G), three GHF16 (Figure 5H), and two GHF38 (Figure 5I). These GHFs showed different patterns of distribution along the digestive tract of P. canaliculata. One peptide with a different pattern of expression was found in hemicellulase A0A2T7NWR9 of GHF10 (Supporting Figure 1). GHF2 and GHF3 showed higher values in digestive gland and in crop and style sac contents. Specifically, GHF5, GHF10, and GHF16 hemicellulases showed higher relative abundance in the crop and style sac with a pattern similar to that observed for the GHF9, GHF10 and GHF45 cellulases. In contrast, GHF38 hemicellulases were exclusively found in coiled gut and both salivary and digestive glands. GHF13 comprised six amylases (Figure 5J) that were highly expressed in the crop and style sac contents. Three maltases of GHF13 (Figure 5K) were detected in salivary glands and had low percent values in samples taken from digestive contents. GHF31 maltases (Figure 5J) were highly expressed in crop and style sac contents and lower in

ACS Paragon Plus Environment

12

Page 13 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

digestive gland, with the exception being the maltase A0A2T7NWB7 expressed at tissue level (SG > DG = CG). Eight galactosidases were grouped within either the GH1 (Figure 5N) and GH35 (Figure 5O) families. The five GHF1 proteins were highly expressed in crop and style sac contents, while the three GHF35 proteins were exclusively expressed in digestive tissues, wherein the expression levels followed the pattern DG > SG > CG. The six monitored fucosidases of the GHF29 (Figure 5M) showed an expression pattern similar to that of the GH35 galactosidases. An analysis of variance for glycoside hydrolase peptides across the digestive contents and tissues from P. canaliculata was undertaken (GLMMs - LSD of Fischer; Supporting Table 6). Statistically significant differences in the variance were often found when GLMMs was run using the individual as the random factor. For some peptides (grey; Supporting Table 6) the differences may also be explained by the individual factor. Generally, all peptides of the same hydrolase showed similar abundance and distribution patterns. Figure 6 shows some examples of each type of glycosidase. For maltase A0A2T7NWB7 and hemicellulase Pc127144, the salivary glands showed a significantly higher abundance than digestive gland and coiled gut (Figure 6A and B, respectively). For fucosidase Pc125661, hemicellulase A0A2T7P4T1, and galactosidase A0A2T7PBX1 (Figure 6C-E), the digestive gland showed often significantly higher peptides abundance than salivary glands and coiled gut, but the fucosidase Pc125661 was not found in either of digestive contents, crop or style sac. Peptides of maltase Pc10364, cellulase A0A2T7PPN2, hemicellulase Pc123712, amylase Pc109575 and galactosidase A0A2T7PCL9 (Figure 6F-J) were mainly found in digestive contents, with the crop content significantly higher than digestive gland.

4. DISCUSSION

ACS Paragon Plus Environment

13

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 38

This work reports the diversity, occurrence, and origin of intracellular and extracellular enzymes involved in the digestion of carbohydrates along the digestive tract of the apple snail P. canaliculata, using a proteomic approach and database searching employing a custom database, comprising the public protein database of Molluscan origin (UniProt) and the transcript-derived database of P. canaliculata29. This proteomic approach revealed six different types of glycosidases that are involved in the digestive processes of this invasive snail, i.e. the cellulases, hemicellulases, amylases, maltases, fucosidases and galactosidases.. The digestive gland and contents of the crop and style sac revealed a higher diversity of glycosidase-derived peptides (Figures 2 and 3). Peptides representing 81 glycosidases (46 from the public databases and 35 uniquely from the transcriptome database) distributed across 13 GHFs were selected and quantified in digestive contents and glands of this snail (Figures 4-6). Conserved expression patterns of each type of glycosidase and their families were found in the digestive tract of the apple snail P. canaliculata, despite the high biological variability of some peptides (Supplementary Table 2).

4.1. Diversity, occurrence and origin of glycosidases 4.1.1. Cellulases Cellulose is an energetically-valuable carbohydrate biopolymer found in algae and plant cell walls38. In this work, thirty cellulases from families GHF9, 10 and 45 were found in crop and style sac contents and to a lesser extent within the digestive gland. They were comparatively lower or null in samples taken from salivary glands and coiled gut. These findings indicate that cellulases are mainly originated in the digestive gland and released through their two main ducts into the stomach vestibule and in consequence they are found downstream in the style sac. Cellulases that were found upstream in the crop content can be explained by the absence of a sphincter separating

ACS Paragon Plus Environment

14

Page 15 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

the stomach and posterior oesophagus and the backward flux driven by pumping action of the muscular gizzard18. The genome of P.canaliculata has revealed expanded cellulase gene families39, which may give rise to increased gut enzyme versatility found in this study. It has also been shown that cellulose expression is influenced by the feeding state (i.e. satiated or starved) of apple snails40. In addition, cellulases A0A2T7PPN6 (GHF9), A0A2T7NYY1 (GHF10), A0A2T7NZR0 (GHF10) and Pc89752 (GHF10) were uniquely found in the digestive contents indicating a possible origin in unicellular glands of the gut or in commensal symbiotic organisms41. To date, a bacterial cellulase secreted by Bacillus sp. from A. crossean has been characterised at molecular and biochemical levels41, 42. Together, these findings could explain why P. canaliculata is able to survive on an exclusive cellulose diet for 60 days43.

4.1.2. Hemicellulases Hemicelluloses are complex polysaccharides present in the plant diet of apple snails and they are putative sources of diverse monosaccharides (glucose, xylose, mannose, or galactose) bound by different linkages44. In this work, hemicellulases with different sugar-specificities were detected along the digestive tract of P. canaliculata. Almost all hemicellulases (except two mannosidases from GHF38) originated in the digestive gland fulfilling extracellular functions at the crop and stomach levels. Two different quantitative profiles were found in digestive gland: representatives from GHF2 (glucuronidases and mannosidases) and GHF3 (xylosidases) showed higher tissue values than contents while representatives from GHF5 (mannosidases), 10 (xylanases) and 16 (β1,3-glucanases) families showed lower tissues values than contents. It is possible that the digestive gland expressed different hemicellulose isoforms which are associated with the ingested food types or the feeding state of the snails40; however, future research will have to clarify this point.

ACS Paragon Plus Environment

15

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 38

Mannosidases from GHF38 were detected in the salivary glands, but these peptides might be overestimated since this tissue showed about 7x less soluble proteins compared to other tissues and the samples were normalized by protein load prior to digestion and analysis. Even though α- or β-mannosidases, and β-xylosidases activities have been reported from viscera of P. canaliculata24-26 and β-1,3-glucanases from digestive fluids of marine gastropods45-47, it is still possible that these enzymes originate in the digestive gland because it is the most voluminous visceral organ43.

4.1.3. Amylases and maltases As expected, amylases and maltases were distributed along the digestive tract of P. canaliculata. These enzymes work together to hydrolyse starch and glycogen arising from the consumption of plants, algae, periphyton and animal carrion, respectively16, 48. In this study, six amylases of the GHF13 (Figure 5J) were highly expressed in the crop and style sac contents indicating that these enzymes are working within the gut extracellular medium. The presence of six amylases may be due to the structure of the amylase gene and their differential tissue expression as was reported in marine bivalves Crassostrea gigas49 and Pecten maximus50. Although the tissue origin of amylases cannot be assigned in this study (Figures 4C-5J), previous research has reported the occurrence of a powerful amylase in the saliva of P. canaliculata17. It is possible that the disaccharides hydrolysis in P. canaliculata begin in the snail mouth since the salivary glands showed a high abundance from peptides of three maltases (A0A2T7NZL2 and A0A2T7PNS7 from GHF13; A0A2T7NWB7 from GHF31). In addition, the maltose released from dietary starch is presumably hydrolysed downstream of the digestive tract since nine maltases (one GHF13 member; eight GHF31 members) were often found in the crop and style sac contents,

ACS Paragon Plus Environment

16

Page 17 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

and digestive gland. Amylase, maltase and other glycosidases activities have been reported in the salivary gland from adult individuals of the freshwater gastropod Biomphalaria straminea51. Also, an α-D-glycosidase, with high maltase but null amylase enzymatic activity, was isolated and characterised in the visceral mass from the marine gastropod Aplysia fasciata52; however, the tissue origin was unknown.

4.1.4. Fucosidases Fucosidases are known to occur in various gastropods (Haliotis rufescens53; Charonia lampas54; Turbo cornutus55, Aplysia kurodai56), in the bivalve Chamelea gallina57, and in the cephalopod Octopus vulgaris58. The present study identified six fucosidases of the GHF29 which showed tissue-specific expression that were practically null in the crop and style sac contents (Figure 5M). These results agree with previous reports of fucosidase isoforms from haemocytes, digestive gland (hepatropancreas), and viscera of P. canaliculata24, 59-61. Although the physiological role of these enzymes is still unknown, some hypotheses suggest: (a) that fucosidades were associated with lysosomal compartments, where post-translational modifications of glycoproteins occur61, 62; or (b) that dietary algal polysaccharides with high fucose content (as fucoidan) were absorbed by endocytosis from digestive epithelial cells and then hydrolysed in granular and vacuolar compartments63, 64.

4.1.5. Galactosidases In this study, eight galactosidases from GHF1 (Figure 5N) and GHF35 (Figure 5O) were monitored. Five GHF1 were highly expressed in crop and style sac contents and to a lesser extent in digestive gland. This finding indicates that galactosidases GHF1 would participate in the

ACS Paragon Plus Environment

17

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 38

extracellular hydrolysis of galactose-rich dietary polysaccharides deriving from plant cell walls (i.e. galactomannan and galactoglucomannan)65. On the other hand, digestive tissues showed exclusively GHF35 galactosidases (Figure 5O) indicating that they may participate in the metabolism of storage carbohydrates. This is the first report of galactosidases in digestive tissues of an adult ampullariid, even though α- and β-galactosidase activities are known to occur during embryo development of Pomacea sp.66. Also, β-galactosidase activity has been reported in juveniles of the freshwater snail, Biomphalaria glabrata67. This enzymatic distribution appears to be associated with the presence of galactan, a storage carbohydrate composed mainly of galactose, which is produced in the albumen glands and accumulates around the eggs of the spawn of various snails68, 69.

4.2. Conclusions and future work The glycosidase diversity found in the luminal contents along the digestive tract indicates that extracellular digestion of complex carbohydrates from a diet is physiologically important in the apple snail P. canaliculata. Also, the consecutive occurrence of glycosidases in contents and glandular tissues ensure complete digestion of varied carbohydrate sources by this polyphagous snail. The glycosidases were notably lower in both diversity and abundance in the coiled gut of P. canaliculata. They likely originate in the salivary and digestive glands. The understanding of the digestive physiology of P. canaliculata has increased in the last decades18, 19, 70, and recent advances in ‘omics’ technologies can contribute to further elucidate key digestive mechanisms of P. canaliculata20, 29, 39. It has been postulated that biological invasions may be prevented with possible applications of physiological knowledge71. The comprehensive investigation of the enzymes involved in the digestive physiology of P. canaliculata provides a

ACS Paragon Plus Environment

18

Page 19 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

foundation for the development of strategies aiming to eradicate and/or control this invasive species.

ACS Paragon Plus Environment

19

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 38

FIGURES

Figure 1. A schematic view of workflow. Snail acclimation and sampling have been reported previously (18). After protein extraction of samples, tryptic digestion was performed using a filteraided sample preparation. The peptides obtained were analysed by LC-MS/MS and then identified by searching against a custom database, i.e. public protein database of Mollusca (Uniprot) and transcript database of P. canaliculata (29). The unannotated proteins were searched in BLASTp. After manual removal of redundant proteins, a final protein list was obtained.

ACS Paragon Plus Environment

20

Page 21 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Figure 2. Classification and distribution of glycosidases. Proteins from the global proteomic profiling were classified taking account their functional groups. The distribution along the digestive tract is shown: 1, salivary gland; 2, crop content; 3, digestive gland; 4, style sac content; 5, coiled gut.

ACS Paragon Plus Environment

21

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 38

Figure 3. Unique and shared peptides of glycosidases from digestive tract of P.canaliculata. Peptides derived from glycosidases that were identified by LC-MS/MS and then compared among them. Each organ can be differentiated by a unique peptide profile. The value in parentheses indicates the number of peptides found in each sampled tissue/content: SG, salivary gland; Cr, crop content; DG, digestive gland; SS, style sac content; CG, coiled gut.

ACS Paragon Plus Environment

22

Page 23 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Figure 4. Peptide distribution and abundance along the digestive tract. A total of 250 peptides from 81 glycosidases were selected for LC-MRM-MS quantitation. The abundance of each peptide was calculated as the sum of MRM peak areas. This value was then compared to the average of all measurements (i.e. 100% represents the mean area for each peptide along the digestive tract). Those peptides that showed a different pattern within each family are shown in grey.

ACS Paragon Plus Environment

23

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 38

ACS Paragon Plus Environment

24

Page 25 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Figure 5. Peptide distribution and abundance according to glycoside hydrolase families. Peptides were classified in 13 glycoside hydrolase families (GHFs) based on the amino acid sequences similarities in the carbohydrate-active enzymes database (CAZy; http://www.cazy.org/). The abundance of each peptide was calculated as the sum of MRM peak areas. This value was then compared to the average of all measurements (i.e. 100% represents the mean area for each peptide along the digestive tract). Those peptides that showed a different pattern within each family are shown in grey. The lines and error bars represent the mean ± SD for each tissue/content type.

ACS Paragon Plus Environment

25

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 38

Figure 6. Peptide expression patterns of selected glycosidases along the gastrointestinal tract of P. canaliculata. The MRM peak area (mean ±SD) for two peptides from each protein across the digestive tract was evaluated statistically using multiple comparisons. Different letters

ACS Paragon Plus Environment

26

Page 27 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

indicate statistically significant differences among groups. The significance level was fixed at p < 0.05.

ASSOCIATED CONTENT Supporting Table 1. Peptide identification metrics. Supporting Table 2. BLASTp results for proteins identified as glycosidases. Supporting Table 3. MRM transitions for peptides monitored by LC-MRM-MS. Supporting Table 4. Co-efficient of variance of peptides monitored in MRM-MS from samples of the digestive tract of P. canaliculata. Supporting Table 5. Occurrence of glycosidase peptides in the digestive tract of P. canaliculata. Supporting Table 6. Quantitation and statistical analysis of 81 glycosidases from P. canaliculata. Supporting Figure 1: Quantitation of peptides showing a different behaviour (grey) along the digestive tract.

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes The authors declare no competing financial interest.

ACKNOWLEDGMENT This work was supported by grants from Universidad Nacional de Cuyo; Consejo Nacional de Investigaciones Científicas y Técnicas and Fondo Nacional de Ciencia y Técnica of Argentina.

ACS Paragon Plus Environment

27

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 38

ABBREVIATIONS AmBic, Ammonium Bicarbonate; CG, Coiled Gut; Cr, Crop; DG, Digestive Gland; DTT, Dithiothreitol; FDR, false discovery rate; GHF, Glycoside Hydrolase Family; GLMMs, Generalized Linear Mixed Models; LSD, Least Significant Difference; MRM, multiple reaction monitoring; PEB, Protein Extraction Buffer; SG, Salivary Gland; SS, Style Sac.

REFERENCES 1.

Hayes KA, Burks RL, Castro-Vazquez A, Darby PC, Heras H, Martin PR, et al. Insights

from an Integrated View of the Biology of Apple Snails (Caenogastropoda: Ampullariidae). Malacologia. 2015; 58(1-2):245-302. 2.

Lowe S, Browne M, Boudjelas S, De Poorter M. 100 of the world's worst invasive alien

species: A selection from the global invasive species database. Published by The Invasive Species Specialist Group (ISSG) a specialist group of the Species Survival Commission (SSC) of the World Conservation Union (IUCN), 12pp. Aliens: The Invasive Species Bulletin. 2000; 12. 3.

Rawlings TA, Hayes KA, Cowie RH, Collins TM. The identity, distribution, and impacts

of non-native apple snails in the continental United States. BMC Evol Biol. 2007; 7(1):97. 4.

Nentwig W, Bacher S, Kumschick S, Pysek P, Vila M. More than "100 worst" alien species

in Europe. Biological Invasions. 2018; 20(6):1611-21. 5.

Yusa Y, Wada T. Impact of the introduction of apple snails and their control in Japan.

Naga, The ICLARM Quarterly. 1999; 22(3):9-13. 6.

Horgan FG, Stuart AM, Kudavidanage EP. Impact of invasive apple snails on the

functioning and services of natural and managed wetlands. Acta Oecologica-International Journal of Ecology. 2014; 54:90-100.

ACS Paragon Plus Environment

28

Page 29 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

7.

Naylor R. Invasions in agriculture: Assessing the cost of the golden apple snail in Asia.

Ambio. 1996; 25(7):443-8. 8.

Giraud-Billoud M, Castro-Vazquez A, Campoy-Diaz AD, Giuffrida PM, Vega IA.

Tolerance to hypometabolism and arousal induced by hibernation in the apple snail Pomacea canaliculata (Caenogastropoda, Ampullariidae). Comparative Biochemistry and Physiology Part B, Biochemistry & Molecular Biology. 2018; 224:129-37. 9.

Vega IA, Arribére MA, Almonacid AV, Guevara SR, Castro-Vazquez A. Apple snails and

their endosymbionts bioconcentrate heavy metals and uranium from contaminated drinking water. Environmental Science and Pollution Research. 2012; 19(8):3307-16. 10. Campoy-Diaz AD, Arribere MA, Guevara SR, Vega IA. Bioindication of mercury, arsenic and uranium in the apple snail Pomacea canaliculata (Caenogastropoda, Ampullariidae): Bioconcentration and depuration in tissues and symbiotic corpuscles. Chemosphere. 2018; 196:196-205. 11. Accorsi A, Ottaviani E, Malagoli D. Effects of repeated hemolymph withdrawals on the hemocyte populations and hematopoiesis in Pomacea canaliculata. Fish & Shellfish Immunology. 2014; 38(1):56-64. 12. Albrecht E, Carreno N, Castro-Vazquez A. A quantitative study of copulation and spawning in the South American apple-snail. Veliger. 1996; 39(2):142-7. 13. Estebenet AL, Martín PR. Pomacea canaliculata (Gastropoda: Ampullariidae): life-history traits and their plasticity. Biocell. 2002; 26(1):83-9. 14. Saveanu L, Martín PR. Pedal surface collecting as an alternative feeding mechanism of the invasive apple snail Pomacea canaliculata (Caenogastropoda: Ampullariidae). Journal of Molluscan Studies. 2012; 79(1):11-8.

ACS Paragon Plus Environment

29

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 38

15. Robert HC. Apple snails (Ampullariidae) as agricultural pests: their biology, impacts and management. Barke Centre for Agriculture and Bioscience International. 2002:145-92. 16. Cazzaniga N, Estebenet A. Revisión y notas sobre los hábitos alimentarios de los Ampullariidae (Gastropoda). Historia Natural. 1984; 4(22):213-24. 17. Andrews EB. The functional anatomy of the gut of the prosobranch gastropod Pomacea canaliculata and of some other pilids. Proceedings of the Zoological Society of London; 1965; 145(1):19-36. 18. Godoy MS, Castro-Vasquez A, Vega IA. Endosymbiotic and host proteases in the digestive tract of the invasive snail Pomacea canaliculata: diversity, origin and characterization. PLoS One. 2013; 8(6):e66689. 19. Imjongjirak C, Amparyup P, Sittipraneed S. Cloning, genomic organization and expression of two glycosyl hydrolase family 10 (GHF10) genes from golden apple snail (Pomacea canaliculata). DNA Sequence. 2008; 19(3):224-36. 20. Liu C, Zhang Y, Ren Y, Wang H, Li S, Jiang F, et al. The genome of the golden apple snail Pomacea canaliculata provides insight into stress tolerance and invasive adaptation. GigaScience. 2018; 7(9):giy101. 21. Li YH, Guo R, Yin QY, Ding M, Zhang SL, Xu GJ, et al. Purification and Characterization of Two Endo‐β‐1, 4‐glucanases from Mollusca, Ampullaria crossean. Acta Biochimica et Biophysica Sinica. 2005; 37(10):702-8. 22. Guo R, Ding M, Zhang S-L, Xu G-j, Zhao F-k. Molecular cloning and characterization of two novel cellulase genes from the mollusc Ampullaria crossean. Journal of Comparative Physiology B. 2008; 178(2):209-15.

ACS Paragon Plus Environment

30

Page 31 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

23. Wang J, Ding M, Li Y-H, Chen Q-X, Xu G-J, Zhao F-K. Isolation of a multi-functional endogenous cellulase gene from mollusc, Ampullaria crossean. Acta Biochimica et Biophysica Sinica. 2003; 35(10):941-6. 24. Hirata K, Nakahara Y, Kimura Y, Funatsu G. Purification and some properties of a betaxylosidase and an alpha-fucosidase from apple snails (Pomacea canaliculata). Bioscience Biotechnology and Biochemistry. 1996; 60(2):249-54. 25. Hirata K, Aso Y, Ishiguro M. Purification and some properties of beta-mannosidase, betaN-acetylglucosaminidase, and beta-galactosidase from apple snails (Pomacea canaliculata). Journal of the Faculty of Agriculture Kyushu University. 1998; 42(3-4):463-72. 26. Hirata K, Aso Y, Ishiguro M. Properties of alpha-mannosidase partially purified from the apple snail, Pomacea canaliculata. Bioscience Biotechnology and Biochemistry. 1998; 62(11):2242-5. 27. Wiśniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nature methods. 2009; 6(5), 359. 28. Colgrave ML, Byrne K, Howitt CA. Liquid chromatography-mass spectrometry analysis reveals hydrolyzed gluten in beers crafted to remove gluten. Journal of Agricultural and Food Chemistry. 2017; 65(44):9715-25. 29. Sun J, Wang M, Wang H, Zhang H, Zhang X, Thiyagarajan V, et al. De novo assembly of the transcriptome of an invasive snail and its multiple ecological applications. Molecular Ecology Resources. 2012; 12(6):1133-44. 30. Shilov IV, Seymour SL, Patel AA, Loboda A, Tang WH, Keating SP, et al. The Paragon Algorithm, a next generation search engine that uses sequence temperature values and feature

ACS Paragon Plus Environment

31

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 38

probabilities to identify peptides from tandem mass spectra. Mol Cell Proteomics. 2007; 6(9):1638-55. 31. Deutsch EW, Csordas A, Sun Z, Jarnuczak A, Perez-Riverol Y, Ternent T, et al. The ProteomeXchange consortium in 2017: supporting the cultural change in proteomics public data deposition. Nucleic Acids Research. 2016; 54(D1):D1100-D1106. 32. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Research. 2018; 47(D1):D442-D50. 33. MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, et al. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010; 26(7):966-8. 34. Soni H, Kango N. Hemicellulases in lignocellulose biotechnology: recent patents. Recent Patents on Biotechnology. 2013; 7(3):207-18. 35. Juturu V, Wu JC. Insight into microbial hemicellulases other than xylanases: a review. Journal of Chemical Technology & Biotechnology. 2013; 88(3):353-63. 36. Schomburg D, Salzmann M. Enzyme handbook 4, Class 3: Hydrolases. Springer-Verlag Berlin Heidelberg. 1991; 1. 37. Lombard V, Golaconda Ramulu H, Drula E, Coutinho P M, Henrissat B. The carbohydrateactive enzymes database (CAZy) in 2013. Nucleic acids research.2013; 42(D1), D490-D495. 38. Tanimura A, Liu W, Yamada K, Kishida T, Toyohara H. Animal cellulases with a focus on aquatic invertebrates. Fisheries Science. 2013; 79(1):1-13. 39. Sun J, Mu H, Ip JC, Li R, Xu T, Accorsi A, Sánchez Alvarado A, Ross E, Lan Y, Sun Y, Castro-Vazquez A, Vega IA, Heras H, Ituarte S, Van Bocxlaer B, Hayes KA, Cowie RH, Zhao Z,

ACS Paragon Plus Environment

32

Page 33 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Zhang Y, Qian P, Qiu J. Signatures of Divergence, Invasiveness and Terrestralization Revealed by Four Apple Snail Genomes. Molecular Biology and Evolution. 2019, msz084. 40. Yang L, Cheng T-y, Zhao F-y. Comparative profiling of hepatopancreas transcriptomes in satiated and starving Pomacea canaliculata. BMC Genetics. 2017; 18(1):18. 41. Zhang S, Yin Q-y, Li Y-h, Ding M, Xu G-j, Zhao F-k. Molecular and biochemical characterization of Ba-EGA, a cellulase secreted by Bacillus sp. AC-1 from Ampullaria crosseans. Applied Microbiology and Biotechnology. 2007; 75(6):1327-34. 42. Li Y-H, Ding M, Wang J, Xu G-j, Zhao F. A novel thermoacidophilic endoglucanase, BaEGA, from a new cellulose-degrading bacterium, Bacillus sp. AC-1. Applied Microbiology and Biotechnology. 2006; 70(4):430-6. 43. Vega IA. Pigmented corpuscles in the midgut gland of Pomacea canaliculata (Caenogastropoda, Ampullariidae): A cyanobacterial endocytobiosis. Universidad Nacional de Cuyo, Mendoza, Argentina. 2006:103. 44. Cosgrove DJ. Growth of the plant cell wall. Nature Reviews Molecular Cell Biology. 2005; 6(11):850. 45. Kumagai Y, Ojima T. Isolation and characterization of two types of β-1, 3-glucanases from the common sea hare Aplysia kurodai. Comparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology. 2010; 155(2):138-44. 46. Kumagai Y, Ojima T. Enzymatic properties and the primary structure of a β-1, 3-glucanase from the digestive fluid of the Pacific abalone Haliotis discus hannai. Comparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology. 2009; 154(1):113-20.

ACS Paragon Plus Environment

33

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 38

47. Pesentseva MS, Kusaykin MI, Anastyuk SD, Sova VV, Zvyagintseva TN. Catalytic properties and mode of action of endo-(1→ 3)-β-d-glucanase and β-d-glucosidase from the marine mollusk Littorina kurila. Carbohydrate Research. 2008; 343(14):2393-400. 48. Saveanu L, Manara E, Martín PR. Carrion consumption and its importance in a freshwater trophic generalist: the invasive apple snail Pomacea canaliculata. Marine and Freshwater Research. 2017; 68(4):752-9. 49. Sellos DY, Van Wormhoudt A. Structure of the of alpha-amylase genes in crustaceans and molluscs: evolution of the exon/intron organization. Biologia-Bratislava-. 2002; 57(2):191-6. 50. Le Moine S, Sellos D, Moal J, Daniel J, San Juan Serrano F, Samian J, et al. Amylase in Pecten maximus (mollusca, bivalves) protein and cDNA characterization; quantification of the expression in the digestive gland. Molecular Marine Biology and Biotechnology. 1997; 6(3):22837. 51. Moura KR, Terra WR, Ribeiro AF. The functional organization of the salivary gland of Biomphalaria straminea (Gastropoda: Planorbidae): secretory mechanisms and enzymatic determinations. Journal of Molluscan Studies. 2004; 70(1):21-9. 52. Andreotti G, Giordano A, Tramice A, Mollo E, Trincone A. Hydrolyses and transglycosylations performed by purified α-D-glucosidase of the marine mollusc Aplysia fasciata. Journal of Biotechnology. 2006; 122(2):274-84. 53. Garcia-Esquivel Z, Felbeck H. Activity of digestive enzymes along the gut of juvenile red abalone, Haliotis rufescens, fed natural and balanced diets. Aquaculture. 2006; 261(2):615-25. 54. Iijima Y, EGAMI F. Purification of α-L-Fucosidase from the Liver of a Marine Gastropod, Charonia lampas. The Journal of Biochemistry. 1971; 70(1):75-8.

ACS Paragon Plus Environment

34

Page 35 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

55. Iijima Y, Muramatsu T, Egami F. Purification of α-l-fucosidase from the liver of a marine gastropod, Turbo cornutus, and its action on blood group substances. Archives of Biochemistry and Biophysics. 1971; 145(1):50-4. 56. Hiraiwa M, Uda Y. Purification and characterization of α-L-fucosidase from the liver of seahare, Aplysia kurodai. Chemical and Pharmaceutical Bulletin. 1985; 33(4):1563-70. 57. Reglero A, Cabezas JA. Glycosidases of Molluscs: Purification and Properties of α‐l‐Fucosidase from Chamelea gallina. European Journal of Biochemistry. 1976; 66(2):379-87. 58. D'aniello A, Hakimi J, Cacace GM, Ceccarini C. The purification and characterization of α-l-fucosidase from the hepatopancreas of Octopus vulgaris. The Journal of Biochemistry. 1982; 91(3):1073-80. 59. Boraldi F, Lofaro FD, Accorsi A, Ross E, Malagoli D. Towards the molecular deciphering of Pomacea canaliculata immunity: first proteomic analysis of circulating hemocytes. Proteomics. 2019; 19(4):e1800314. 60. Endo T, Tsukada T, Hiraiwa M, Uda Y, Kobata A. Purification and characterization of an α-L-fucosidase from Pomacea canaliculata. Archives of Biochemistry and Biophysics. 1993; 302(1):152-60. 61. Hirata K, Aso Y, Yasuda S, Ishiguro M. Purification of Pomacea canaliculata α-fucosidase isoforms with different thermostabilities. Bioscience, Biotechnology, and Biochemistry. 1998; 62(4):807-10. 62. Perrella NN, Cantinha RS, Nakano E, Lopes AR. Characterization of α-l-fucosidase and other digestive hydrolases from Biomphalaria glabrata. Acta Tropica. 2015; 141:118-27.

ACS Paragon Plus Environment

35

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 38

63. Daniel R, Berteau O, Jozefonvicz J, Goasdoue N. Degradation of algal (Ascophyllum nodosum) fucoidan by an enzymatic activity contained in digestive glands of the marine mollusc Pecten maximus. Carbohydrate Research. 1999; 322(3-4):291-7. 64. Berteau O, McCort I, Goasdoué N, Tissot B, Daniel R. Characterization of a new α-Lfucosidase isolated from the marine mollusk Pecten maximus that catalyzes the hydrolysis of α-Lfucose from algal fucoidan (Ascophyllum nodosum). Glycobiology. 2002; 12(4):273-82. 65. Nishinari K, Takemasa M, Osaka JHZ, Takahashi R. 2.19 Storage Plant Polysaccharides: Xyloglucans, Galactomannans, Glucomannans. Comprehensive Glycoscience-From Chemistry to Systems Biology. 2007; 2:613-652. 66. Oliveira FW, Chavante SF, Santos EA, Dietrich CP, Nader HB. Appearance and fate of a β-galactanase, α, β-galactosidases, heparan sulfate and chondroitin sulfate degarding enzymes during embryoding development of the mollusc Pomacea sp. Biochimica et Biophysica Acta (BBA)-General Subjects. 1994; 1200(3):241-8. 67. Bretting H, Jacobs G, Müller S, Breitfeld O. The galactan-degrading enzymes in the snail Biomphalaria glabrata. Comparative Biochemistry and Physiology B, Comparative Biochemistry. 1991; 99(3):629-36. 68. Bretting H, Messer M, Bornaghi L, Kröger L, Mischnick P, Thiem J. Galactan biosynthesis in snails: a comparative study of β-(1→ 6) galactosyltransferases from Helix pomatia and Biomphalaria glabrata. Journal of Comparative Physiology B. 2000; 170(8):601-13. 69. Bretting H, Jacobs G, Benecke I, König WA, Thiem J. The occurrence of L-galactose in snail galactans. Carbohydrate Research. 1985; 139:225-36.

ACS Paragon Plus Environment

36

Page 37 of 38 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

70. Li YH, Yin QY, Ding M, Zhao FK. Purification, characterization and molecular cloning of a novel endo-beta-1,4-glucanase AC-EG65 from the mollusc Ampullaria crossean. Comparative Biochemistry and Physiology B-Biochemistry & Molecular Biology. 2009; 153(2):149-56. 71. Lennox R, Choi K, Harrison PM, Paterson JE, Peat TB, Ward TD, et al. Improving sciencebased invasive species management with physiological knowledge, concepts, and tools. Biological Invasions. 2015; 17(8):2213-27.

ACS Paragon Plus Environment

37

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 38

“For TOC Only”

ACS Paragon Plus Environment

38