Dynamic Rearrangement in Snake Venom Gland Proteome: Insights

Aug 30, 2016 - ... Viperidae; Wied-Neuwied 1824) snake venom after up to 54 years of storage ... among specimens from Sri Lanka, India and Pakistan...
0 downloads 0 Views 2MB Size
Subscriber access provided by Northern Illinois University

Article

Dynamic rearrangement in snake venom gland proteome: insights into Bothrops jararaca intraspecific venom variation César Augusto de-Oliveira, Daniel R. Stuginski, Eduardo S. Kitano, Débora AndradeSilva, Tarcísio Liberato, Isabella Fukushima, Solange M.T. Serrano, and André Zelanis J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/acs.jproteome.6b00561 • Publication Date (Web): 30 Aug 2016 Downloaded from http://pubs.acs.org on September 1, 2016

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 free 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 accessible to all readers and 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.

Journal of Proteome Research 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 45

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

Dynamic rearrangement in snake venom gland proteome: insights into Bothrops jararaca intraspecific venom variation

César Augusto-de-Oliveira1, Daniel R. Stuginski2, Eduardo S. Kitano3, Débora Andrade-Silva3, Tarcísio Liberato1, Isabella Fukushima1, Solange. M. T. Serrano3, André Zelanis1*

1

Laboratório de Proteômica Funcional, Departamento de Ciência e Tecnologia,

Universidade Federal de São Paulo (ICT-UNIFESP), São José dos Campos, SP, Brazil. 2

Laboratório de Herpetologia, Instituto Butantan, São Paulo, Brazil.

3

Laboratório Especial de Toxinologia Aplicada, Center of Toxins, Immune-

Response and Cell Signaling (CeTICS), Instituto Butantan, São Paulo, Brazil.

*

Corresponding author: André Zelanis Department of Science and Technology, Federal University of São Paulo (ICT-UNIFESP), São José dos Campos, SP, Brazil, Rua Talim, 330, 12231-280, Tel.: +55 11 3385-4135 ext. 9727. E-mail: [email protected] Running title: Bothrops jararaca venom gland proteome 1 ACS Paragon Plus Environment

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 45

Abstract We carried out an analysis of the venom gland proteome of Bothrops jararaca taking into account two distinct phases of its ontogenetic development (i.e. newborn and adult) and the marked sexual dimorphism recently reported on its venom proteome. Proteomic data analysis showed a dynamic rearrangement in the proteome landscape of B. jararaca venom gland upon development and gender-related changes. Differentially expressed proteins covered a number of biological pathways related to protein synthesis, including proteins associated to transcription and translation, which were found to be significantly higher expressed in the newborn venom gland. Our results suggest that the variation in the expression levels of cellular proteins might give rise to an even higher variation in the levels of the expressed toxins. Upon ageing the venom gland proteome repertoire related to the protein synthesis together with ecological traits would have an impact on the toxin repertoire, which, in the case of B. jararaca species, would enable the species to deal with different prey types during its lifespan. Proteomic data are available via ProteomeXchange with identifier PXD004186.

Keywords: snake venom proteome, venom gland proteome, venom variability, Bothrops jararaca

2 ACS Paragon Plus Environment

Page 3 of 45

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

Abbreviations svVEGF: snake venom Vascular Endothelium Growth Factor; 5NUCL: Nucleotidase; AP: Aminopeptidase; BPP: Bradykynin Potentiating Peptides; CRISP: Cysteine-Rich Secretory Protein; CTL: C-Type lectin; CVF:Cobra Venom Factor; DIEST: Diesterase; DPP: Dipeptidyl Peptidase; NGF: Nerve Growth Factor; PLAI: phospholipase inhibitor gamma-like; PLB: phospholipase B; SVMP: Snake Venom Metalloproteinase; SVSP: Snake Venom Serine Proteinase; LAAO: L-Amino Acid Oxidase.

3 ACS Paragon Plus Environment

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 45

Introduction Investigators have long appreciated the composition and diversity of snake venoms and a number of venom proteomes have been reported so far.1-4 As a biological product derived from an intricate process of protein synthesis and secretion, the composition of snake venoms may vary according to several factors, such as habitat, gender, ontogenetic development and others.5-11 Owing to their highly specialized secretory function, snake venom glands represent exceptional models for the study of protein synthesis, secretion and long-term storage of potentially harmful components, such as proteases and myotoxins.12-18 Moreover, venom gland secretory cells encode rich information regarding toxin diversity as well as comprise the molecular machinery responsible for the venom production. As such, they are the main biological source of starting material for the study of the gene expression profiles. Nevertheless, with some exceptions, little attention has been given to the cellular transcripts/proteins composing snake venom glands or the biological processes related to snake venom synthesis. Even though a number of venom gland transcriptomes have been characterized so far, little is known regarding the composition of the venom gland proteome of any snake species. Bothrops jararaca is among the most studied South American snake species both in terms of ecological and toxinological points of view.10,11,19-26 According to recent epidemiological data from Brazilian Ministry of Health, this species is responsible to the majority of the accidents with humans in Southeastern Brazil, a highly populated region.27 Although this species is primarily generalist, it displays a significant change in feeding habits during lifespan with newborn specimens feeding on small ectothermic prey such as small lizards, amphibians 4 ACS Paragon Plus Environment

Page 5 of 45

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 arthropods and adults feeding mainly on mammals.19,20,28 Such changes correlate with the venom proteome rearrangement upon neonate to adult transition and may reflect different strategies to deal with distinct prey types.10,23 It has been shown that in B. jararaca noradrenergic innervation is a key activator of venom gland and that the stimulation of α1-adrenoreceptors triggers snake

venom

production

phosphatidylinositol

cycle

in

4,5-bisphosphate

secretory hydrolysis

cells, and

by

activating

ERK

signaling

pathway.16,29 Interestingly, Luna and coworkers29 reported sexual dimorphism in the activation of NFκB and AP-1 transcription factors in adult male and female B. jararaca venom gland. Furthermore, the authors also reported gender-related differences in the venom production cycle after 4 and 7 days after venom extraction in adult female and male snakes, respectively. Indeed, sexual dimorphism in adult B. jararaca venom was recently reported in which quantitative label-free proteomic analysis revealed the higher expression of Nerve Growth Factor in male venom whereas C-type lectins displayed more than 3 fold increase in expression in female venom.26 Moreover, peptidomic analysis confirmed the ubiquitous presence of the four Bradykinin Potentiating Peptides (BPPs) that lack the C-terminal Q-I-P-P sequence only in the female venom as gender molecular markers.26,30 Despite the number of reports on qualitative/quantitative variation in B. jararaca venom proteome as well as its venom gland transcriptomic profiles22,25 the molecular mechanisms underlying venom variability in this species remains poorly understood. Therefore, the aim of this work was to analyze the venom gland proteome of B. jararaca taking into account two distinct phases of its ontogenetic development (i.e. newborn and adult) and the sexual dimorphism 5 ACS Paragon Plus Environment

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 45

recently reported on its venom proteome26. Our results point to a striking rearrangement of B. jararaca venom gland proteome in two distinct phases of the ontogenetic development as well as gender-related differences and shed light into the remarkable intraspecific venom variability found in this species.

Material and Methods Bothrops jararaca venom and venom gland samples Venom from 694 two-week old newborns (359 male and 335 female specimens) and 110 adults (49 male and 61 female specimens older than 3 years) from São Paulo State (Brazil) was used in this study. The venom was extracted and centrifuged for 30 min at 2,000 x g, at 4 °C, to remove any scales or mucus, lyophilized and stored at -20°C until use. Venom protein concentrations were determined using the Bradford reagent (Sigma, USA) and bovine serum albumin (Sigma, USA) as a standard. B. jararaca specimens were obtained from the Herpetology Laboratory, Instituto Butantan, São Paulo, Brazil. Twelve animals were used to obtain venom glands, including 6 newborns (up to 45 days old; 3 male and 3 female specimens) and 6 wild caught adults (3 male and 3 female specimens), from different geographic regions of São Paulo State, Brazil (Supplementary table 1). The venom was extracted and 3 days later the animals were sacrificed by CO2 inhalation. The venom glands were carefully dissected, frozen in liquid nitrogen and kept at -80°C until use. All animal work has been conducted in agreement with the Ethical Principles in Animal Research, adopted by the Brazilian College of Animal Experimentation and was approved by the Ethical Committee for Animal Research of Butantan Institute (protocol # 1371/15). 6 ACS Paragon Plus Environment

Page 7 of 45

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

Protein extraction and quantitation Three milliliters of lysis buffer (50 mM HEPES, 200 mM NaCl, 2% CHAPS, pH 7.5) containing the HALTTM protease Inhibitor cocktail (Thermo Fisher, USA) and 5 mM EDTA (final concentration) were added to each venom gland pool and the samples were subjected to homogenization using a tissue homogenizer Polytron PT3100 (Kinematica AG, Switzerland) using 18,000 rpm cycles for 30 s in ice bath. Homogenates were clarified by centrifugation at 14,000 x g for 10 min, at 4°C, to remove debris and the clear supernatant was carefully removed. Due to the low weight of the individual newborn venom glands (data not shown), we chose to pool the 6 pairs of newborn venom glands (3 male + 3 female), whereas a total of 6 venom glands was used for the adult comparison: a male sample (a pool of 3 venom glands derived from 3 male individuals), a female sample (a pool of 3 venom glands derived from 3 female individuals) and an adult pool (a mixture where equal amounts of proteins were added from each adult pool to make up an 1 mg/mL protein solution) Protein concentrations were determined using the Bradford reagent (Sigma, USA) and bovine serum albumin (Sigma, USA) as a standard.

SDS-polyacrylamide gel electrophoresis and Western blot analysis SDS-PAGE was carried out according to Laemmli31 and proteins were stained with Colloidal Coomassie Blue G-250. Western blot analysis was carried out as described elsewhere32 using anti-bothropic antivenom from Instituto Butantan (batch # 0610188), generated by the immunization of horses with a mixture of adult venoms of B. jararaca (50%), B. alternatus (12.5%), B. jararacussu

7 ACS Paragon Plus Environment

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 8 of 45

(12.5%), B. moojeni (12.5%), B. neuwiedi (synonyms B. n. goyazensis, B. n. meridionalis, B. n. paranaensis and B. n. urutu)33 (6.25%) and B. pauloensis (6.25%). The secondary antibody used was anti-Horse IgG- peroxidase (Sigma, USA) with PierceTM ECL substrate (Thermo Scientific, USA). As a negative control, the membrane was incubated with non-immune horse serum as the primary antibody, followed by the incubation with the same secondary antibody.

In-solution trypsin digestion and reductive isotopic dimethylation labeling The in-solution trypsin digestion was performed according to the protocol described by Kleifeld et al.34 with slight modifications. Briefly, a solution of 6 M guanidine hydrochloride (GuHCl) was added to a sample of 100 µg of protein from each venom gland sample to a final concentration of 3 M GuHCl, followed by the addition of 5 mM dithiothreitol (DTT) (final concentration). The mixture was incubated at 65 °C for 60 min. Iodoacetamide (IAA) was then added to a final concentration of 15 mM and the samples were incubated for 60 min at room temperature, in the dark. To quench the excess of IAA, DTT was added to a final concentration of 15 mM. Clean-up of samples was performed by the addition of ice cold acetone (8 volumes) and methanol (1 volume), followed by the incubation of samples for 3 h at -80° C. After centrifugation at 14,000 x g for 10 min, protein pellets were washed twice with one volume of ice cold methanol and then resolubilized with NaOH solution (final concentration of 2.5 mM), followed by the addition of 50 mM HEPES buffer, pH 7.5, to a final volume of 100 µL. Trypsin (Proteomics grade; Sigma, USA) was added at 1:100 ratio (enzyme/substrate) and protein samples were incubated at 37 °C for 18 h.

8 ACS Paragon Plus Environment

Page 9 of 45

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

Tryptic peptides were differentially labeled via stable-isotope dimethyl labeling, as previously described35. In brief, tryptic peptides were submitted to reductive dimethylation with either light or heavy formaldehyde/cyanoborohydride solutions, as follows: newborn venom gland peptides (light) vs. adult venom gland peptides (heavy); adult female venom gland peptides (light) vs. adult male venom gland peptides (heavy). Tryptic peptides (pH 7.5) from each sample were incubated overnight at 37 °C with either light or heavy sodium cyanoborohydride

(NaBH3CN,

light,

or

NaBD3CN,

heavy)

concentration of 20 mM followed by the addition of formaldehyde or

13

to 12

a

final

CH2O (light)

CD2O (heavy) to a final concentration of 40 mM, resulting in mass

differences of +28.031300 Da and +36.075670 Da for the light and heavylabeled samples, respectively. The reaction was terminated by adding 1 M Tris (pH 6.8; to a final concentration of 200 mM) to each sample and the mixture was incubated for 2 h at 37 °C. Samples were then combined at 1:1 ratio into two pools: (i) newborn + adult and (ii) adult female + adult male. After desalting using C-18 cartridges (3M EmporeTM SPE Extraction disks, USA) peptide samples were dried in a SpeedVac and redissolved in 50 µL of 0.1% formic acid prior

to

nanoflow

liquid

chromatography/tandem

mass

spectrometry

(LC−MS/MS) analysis.

Mass spectrometric analysis An aliquot (5 µL) of the resulting peptide mixture was injected into a trap column packed with C18 (100 µm i.d. × 2 cm) for desalting with 100% solvent A (0.1% formic acid). Peptides were then eluted onto an analytical column (75 µm i.d. × 100 mm) packed in house with Aqua® C-18 5 µm beads (Phenomenex, 9 ACS Paragon Plus Environment

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 10 of 45

USA). Nanoflow liquid chromatography was performed on an Easy nanoLC system (Thermo Fisher Scientific, USA) coupled to an LTQ-Orbitrap Velos mass spectrometer (Thermo Fisher Scientific, USA). Peptides were loaded onto the column with solvent A (0.1% formic acid) and eluted with a 120 min linear gradient from 3 to 30% of solvent B (acetonitrile in 0.1% formic acid) at a flow rate of 200 nL/min. Spray voltage was set at 2.1kV, at 200 °C and the mass spectrometer was operated in data dependent mode, in which one full MS scan was acquired in the m/z range of 300-1650 followed by MS/MS acquisition using Collisional Induced Dissociation (CID) of the fifteen most intense ions from the MS scan. MS spectra were acquired in the Orbitrap analyzer at 60,000 resolution (at 400 m/z). Dynamic exclusion was defined by a list size of 500 features and exclusion duration of 60 s. For the survey (MS) scan AGC target value of 1,000,000 was set whereas the AGC target value for the fragment ion (MS/MS) spectra was set to 10,000 ions. The lower threshold for targeting precursor ions in the MS scans was 3,000 counts. Samples were analyzed as three technical replicates.

Proteomics data processing Mass spectrometric (RAW) data were analyzed with MaxQuant software36 (version 1.5.3.17). A False Discovery Rate (FDR) of 1% was required for both protein and peptide-to-spectrum match identifications. LTQ-Orbitrap Velos raw data were searched against a target database restricted to the taxonomy ‘Serpentes’ (UniProt release 09_2015; 58,895 sequences appended with 15,276 sequences derived from the translation of cDNAs encoding tissue proteins and toxins from B. jararaca venom gland.18 This database was also 10 ACS Paragon Plus Environment

Page 11 of 45

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

combined with the sequences of 245 common contaminants and concatenated with the reversed versions of all sequences. Enzyme specificity was set to trypsin

and

at

least

two

missed

cleavages

were

allowed;

cysteine

carbamidomethylation was selected as fixed modification whereas methionine oxidation, glutamine/asparagine deamidation and protein N-terminal acetylation were selected as variable modifications. Multiplicity was set to 2 to account for the two labeling states that were used (light and heavy dimethyl labeling at peptide N-terminus and lysine side chains). Peptide identification was based on a search with an initial mass deviation of the precursor ion of 7 ppm and the fragment mass tolerance was set to 20 ppm. For protein quantification, a minimum of two ratio counts was set and the ‘Re-quantify’ and ‘match between runs’ functions of MaxQuant software were enabled. As is observed from complex proteomes such as vertebrates, peptides can be shared between homologous proteins or splice variants, leading to “protein groups”. For each protein group in the MaxQuant’s ‘proteinGroups.txt’ file, the first protein entry was selected as representative. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE37 partner repository with the dataset identifier PXD004186.

Bioinformatic analysis All the light and heavy peptide intensities were log2-transformed and quantile-normalized using the ‘preprocessCore’ library in R scripting and statistical environment38,39 to correct for intra-experimental variation. Protein identifications in which light or heavy singletons were found were not considered for further analysis. Statistical analyses were performed using the 11 ACS Paragon Plus Environment

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 45

‘limma’ package in R/Bioconductor.40,41 After fitting a linear model to the data, an empirical Bayes moderated t-test was used for the comparisons (Female x Male and Newborn x Adult), and P values were adjusted for multiple testing with the Benjamini-Hochberg method. Proteins with an adjusted P-value < 0.05 and log2(fold change) > 1 and 170 kDa to

14 ACS Paragon Plus Environment

Page 15 of 45

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

~10 kDa. Since a small portion of MS/MS spectra matched to venom peptides, we decided to perform a Western blot analysis of venom gland extracts using a polyspecific antivenom (anti-bothropic antivenom, produced by Instituto Butantan, Brazil) in order to get a better picture of the residual (or newly synthesized) venom proteome. For this purpose, the best profile was acquired after analyzing venom gland extracts in 5-fold protein excess compared to venom proteins (i.e. 25 µg of venom gland protein extract and 5 µg of venom proteins). Immunostaining of B. jararaca venom gland extracts with antibothropic antivenom revelaed several protein bands ranging from 110 kDa to 10 kDa, with a similar pattern of recognition among all samples (Figure 2B). For several protein bands, the Western blot profile of venom gland extracts somehow mirrored the venom protein profiles, with some exceptions found for high molecular mass proteins. It is interesting to note that the protein profiles at the range of 20-15 kDa were considerably similar between venom gland extracts and venom proteins. On the other hand, a number of high molecular mass proteins (> 100 kDa) were detected in all the samples derived from venom gland extracts, which might correspond to still non-processed precursor protein forms.

Differentially expressed cellular proteins Reductive isotopic dimethylation allowed for the identification and quantification of proteins derived from distinct biological scenarios: (i) two stages of B. jararaca ontogenetic development (newborn and adult) and (ii) gender-related differences. All technical replicates presented an overall good

15 ACS Paragon Plus Environment

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 45

correlation as observed by their Pearson correlation coefficient ranging from 0.89 to 0.94 (Supplementary figures 1 and 2; Supplementary table 3). Almost 8% of the quantified proteins were found to be differentially expressed (51 proteins in total; 21 in newborn and 30 in adult; Supplementary table 4). A protein disulfide-isomerase (U3FBW7) was quantified with the highest expression among all of the identified proteins (Figure 3A).Protein disulfide-isomerases (PDIs) are Endoplasmic Reticulum (ER)-resident proteins that catalyze the formation, breakage and rearrangement of disulfide bonds and therefore are key enzymes responsible for protein folding and the maintenance of protein native structures.43 The identification of PDIs and molecular chaperones among differentially expressed proteins is a recurrent feature in quantitative proteomic studies of eukaryotic samples.44 However, this observation is often related to experimental/biological stress conditions (tumor tissues, tumor derived cell lines or experimentally induced toxic states). In fact, protein synthesis requires continuous activity of PDIs, which might explain the high levels of this protein in tissues where biosynthetic process are taking place, a feature already observed in the proteome of B. jararaca venom gland.17 In relation to gender-related differences the male x female comparison yielded almost 5% of differentially expressed proteins (31 proteins in total; 19 in male and 12 in female; Supplementary table 4). Structural proteins such as collagen alpha-2(I)chain (V8NPI7) and the calcium binding protein, (sorcin-like protein, P30626) were highly expressed in female’s venom gland whereas ribosomal and histone proteins displayed significantly higher expression values in male’s venom gland (Figure 3B). In addition, we also found higher expression of PDIs in male’s venom gland, a feature that might correlate with the higher abundance

16 ACS Paragon Plus Environment

Page 17 of 45

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

of high molecular mass disulfide-linked oligomeric toxins recently reported in the electrophoretic profiles of adult male B. jararaca venom.26 Clustering of newborn and adult differentially expressed proteins revealed specific patterns of expression (Figures 4A) in which proteins belonging to ‘General metabolism’ category were prominent and, in general, displayed a similar pattern of high expression between both samples. Further analysis revealed a striking higher expression of proteins related to biosynthetic process (transcription/translation category) among newborn’s differentially expressed

proteins,

followed

by

proteins

related

to

post-translational

modifications (Figure 4B). Inhibitors such as the antihemorrhagic factor BJ46a (Q9DGI0), small serum protein-1 (D9N567), murinoglobulin-2 (V8NK36) and Complement C3a protein (A0A081DUA7) were highly expressed in adult venom gland, followed by structural proteins such as tropomyosin and collagen. It is possible that these proteins might be present due to the contamination of venom gland extracts with muscle or blood proteins (in the case of structural and inhibitor proteins, respectively), a feature that was previously reported.17 Clustering of adult male and female differentially expressed proteins and the subsequent functional analysis revealed a number of proteins related to the process of transcription/translation as well as PTM/processing in male venom gland whereas ‘general metabolism’ and ‘inhibitor’ categories were prominent in female venom gland (Figure 4C and D).

Toxins Regardless of age and gender-related differences (i.e. when protein identification data were merged together) the venom proteome correlated well 17 ACS Paragon Plus Environment

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 45

with B. jararaca venom proteomes reported so far, both quantitatively (as observed from the number of proteins derived from each toxin family) and qualitatively (as illustrated by the diversity of the identified toxins).10,2225

Proteases comprised the majority of venom toxins (almost 60% of the overall

toxin composition) with SVMPs being the most abundant class (Figure 5A). Actually, among all peptide identifications, only a small portion matched to venom toxins (16%). Whether these toxins resulted from newly synthesized proteins or are remnants from an incomplete emptying of venom gland is difficult to address. We were able to identify peptides belonging to the precursor regions of SVMPs and SVSPs (Supplementary table 5). Sequence logos of the set of the identified peptides resulted in the recurrence of the canonical sequences ‘MCGVT’, and ‘SEL’, part of the cysteine-switch and pro-peptide regions in SVMPs and SVSPs, respectively (Figure 5B). These results suggest that newly synthesized (or remnant) proteinases are inactive either within the secreting cells or within the lumen of the venom gland, or both. In the case of SVMPs, our results are in line with the hypothesis that their processing starts within secretory vesicles and continues into the lumen of the venom gland just after enzyme secretion.45 Moreover, we also found peptides derived from the predicted mature Ntermini form of some C-type lectins and a peptide derived from the signal peptide from a phospholipase A2 enzyme (Supplementary figure 3 and Supplementary table 5). The latter observation illustrates a case where the protein had still not left the endoplasmic reticulum and, therefore, was inside the venom gland cells, a likely newly synthesized protein. These results corroborate with the previous report on B. jararaca venom gland activation in which C-type 18 ACS Paragon Plus Environment

Page 19 of 45

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

lectins were found to be present in the quiescent stage of secreting cells whereas PLA2s where more abundant in the activated stage17. Statistical

analysis

of

differentially

expressed

proteins

strongly

corroborated previous findings related to the venom proteome and venom gland transcriptomes of newborn and adult B. jararaca as well as adult male and female venoms/venom gland transcriptomes/proteomes (Figure 6 A and B) 8,10,22-26,30

. SVMPs comprised more than 60% of the differentially expressed

toxins within newborn venom glands, whereas SVSPs and CTLs accounted for more than 70% of the differentially expressed proteins in adult venom gland (Figure 6A-D). A previous work on the ontogenetic changes in B. jararaca venom showed that the main role of SVMPs in the venom of newborns is acting upon coagulation cascade components, including the activation of Factor X and prothrombin, a feature that likely enables newborn specimens to deal with distinct prey types at an early stage of their ontogeny.10,46 As the snake grows, there might be a diversification in the biological targets of venom toxins as illustrated by the qualitative and quantitative differences found in B. jararaca venom proteome and venom gland transcriptome.23,25 Changes in the dynamics of protein expression and proteome diversity may impact the composition of male and female venoms. In fact, we have recently identified gender molecular markers in adult B. jararaca venoms.26 The quantitative results reported herein also corroborate the higher expression of NGF and CTLs in male and female venoms, respectively (Figures 6E and F). In addition, our results suggest a higher expression of LAAO enzymes in male venom, a feature also observed in the two-dimensional venom protein profile of adult male venoms26. Indeed, changes in the protein composition of B. jararaca 19 ACS Paragon Plus Environment

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 45

venom gland during the venom production cycle were reported, with striking changes occurring 4 and 7 days after venom removal in adult female and male snakes, respectively.29 Therefore, the results reported herein regarding genderrelated differences in adult B. jararaca venom gland might also be related to our (venom) extraction protocol. As sexual maturation and ontogeny display a paramount role in numerous biological features, one might expect that differences in the venom phenotype would mirror the main developmental changes that occur in newborn and adult and male and female specimens. It was shown that the regulation of protein synthesis in B. jararaca venom gland is triggered by noradrenergic stimulation, after the activation of α1-and βadrenoreceptors16,29 and that there is a marked sexual dimorphism in the activation of the transcription factors NFκB and AP-1, whose expression peaks earlier in female B. jararaca venom glands.

Insights into the intraspecific variability in B. jararaca venom proteome We separately evaluated the quantification data from toxins and cellular proteins and, as a general observation, regardless of age or gender, the expression range of cellular proteins showed less variation than the toxin expression range (Figure 7). While cellular transcripts encode rich information that provide key features to understand the molecular basis of snake venom variation, their presence/ abundance does not necessary imply/correlate in the translation of a functional (protein) product. Therefore, a proteomic approach is desirable as a complementary approach. Additionally, several reports have shown a low to

20 ACS Paragon Plus Environment

Page 21 of 45

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

moderate correlation of mRNA and protein levels in eukaryotes including vertebrates47,48. Based on our results it is possible to hypothesize that at an early stage of B. jararaca development the protein biosynthetic process might be prominent in comparison to the adult phase. A number of reports showed that the expression of genes related to protein synthesis is diminished upon aging in the salivary glands of other vertebrates, such as mouse and humans.49-51 The biological meaning behind this hypothesis might be related to the snake metabolism and growing. Growing involves, besides metabolic adjustments, an increase in protein synthesis, which could explain the higher expression of proteins related to transcription and translation when compared to adults. In the same way, due to their higher metabolic rates and the necessity to grow, it is expected that juveniles need to eat more frequently than adults, in order to fuel this process, therefore an increase in the protein synthesis could also play an important role in food processing during the juvenile phase of these animals. The limiting step in protein synthesis is the translation initiation, when the ribosome subunits are assembled to begin the polypeptide formation52. Actually, in addition to a number of ribosomal proteins we found higher expression of several translation initiation factors among newborn venom gland proteins (Supplementary table 4). In any living organism, the maintenance of protein homeostasis

is

acquired

through

a

fine-tuned

balance

of

molecular

mechanisms, namely transcriptional and translational rates coupled with the stability of individual proteins which, in turn, results from the balance between the rate of protein synthesis and the rate of degradation.48,53-55 Therefore, an increase in the rate of protein synthesis should usually be followed by an 21 ACS Paragon Plus Environment

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 45

increase in the degradation rate, which ultimately results in a faster protein turnover and allows a rapid adaptation to a variation of nutritional and/or pathophysiological conditions55. Upon ageing, the cellular proteome repertoire related to the protein synthesis together with ecological traits, would have an impact on the toxin repertoire, which, in the case of B. jararaca species, would enable the species to deal with different prey types during its lifespan.

Conclusions The knowledge of the proteome of B. jararaca venom gland revealed specific profiles related to aging and sexual dimorphism and, in light of the understanding of the diversity of B. jararaca venom, it is possible to infer a dynamic proteome rearrangement of cellular proteins that likely results in the striking intraspecific variability of the venom in this species. It turns out that the variation in the expression levels of cellular proteins gives rise to an even higher variation in the dynamic range of the expressed toxins. Underlying mechanisms that modulate the expression of toxins may likely involve hormonal regulation, which, in turn, might account for the ontogenetic and gender-related differences found in the expression profiles of B. jararaca venom gland. Further analysis (i.e. metabolomics) would also assist in the integration of data related to venom production in this species and bridge the gaps towards the comprehensive knowledge of snake venom proteomes.

SUPPORTING INFORMATION

22 ACS Paragon Plus Environment

Page 23 of 45

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

following

files

are

available

free

of

charge

at

ACS

website http://pubs.acs.org:

Figures (S1-S3) and table S1.docx. Supplementary information on proteomic analysis and biological samples Supplementary tables(S2-S6).xls. Detailed analysis on proteomic data.

Acknowledgements We thank to Dr. Inácio Junqueira-de-Azevedo for providing the B. jararaca transcrpitome dataset used in this work and Dr. Sávio S. Sant’Anna for his assistance in providing access to B. jararaca specimens. This work was supported by grants from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (grant 441804/2014-4), Fundação de Amparo à Pesquisa do Estado de São Paulo (grants 2015/18096-0, 2014/06579-3 and 2013/074671) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

References 1. Fox,

J.W.;

Serrano,

S.M.

Exploring

snake

venom

proteomes:

multifaceted analyses for complex toxin mixtures. Proteomics 2008, 8(4),909-920. 2. Gutiérrez J.M.; Lomonte, B.; León, G.; Alape-Girón, A.; Flores-Díaz, M., Sanz, L.; Angulo, Y.; Calvete, J.J.Snake venomics and antivenomics: 23 ACS Paragon Plus Environment

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 45

Proteomic tools in the design and control of antivenoms for the treatment of snakebite envenoming. J Proteomics 2009, 72(2),165-182. 3. Calvete, J.J. Snake venomics: From the inventory of toxins to biology. Toxicon 2013, 75,44-62. 4. Sunagar, K.; Morgenstern, D.; Reitzel, A.M.; Moran, Y. Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom. J. Proteomics 2016, 135,6272. 5. Mackessy, S.P. Fibrinogenolytic proteases from the venoms of juvenile and adult northern pacific rattlesnakes (Crotalus viridis oreganus). Comp. Biochem.Physiol. 1993, 106 B (1), 181-189. 6. Daltry, J.C.; Wüster, W.; Thorpe, R.S. Diet and snake venom evolution. Nature 1996, 379, 537-540. 7. Mackessy, S.P. Venom ontogeny in the pacific rattlesnakes Crotalus viridis helleri and C.v. oreganus. Copeia 1998, 92-101. 8. Menezes, M.C.; Furtado, M.F.; Travaglia-Cardoso, S. R.; Camargo, A.C.; Serrano, S.M. Sex-based individual variation of snake venom proteome among eighteen Bothrops jararaca siblings. Toxicon 2006, 47, 304–312. 9. Alape-Giron, A.; Sanz, L.; Escolano, J.; Flores-Dias, M. et al. Snake venomics of the lancehead Bothrops asper. Geographical, individual and ontogenetic variations. J. Proteome Res. 2008, 7, 3556-3571. 10. Zelanis, A.; Tashima, A. K.; Rocha, M. M.; Furtado, M. F.; Camargo, A. C.; Ho, P. L.; Serrano, S. M. Analysis of the ontogenetic variation in the venom proteome/peptidome of Bothrops jararaca reveals different strategies to deal with prey. J. Proteome Res. 2010, 9 (5), 2278–2291.

24 ACS Paragon Plus Environment

Page 25 of 45

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

11. Dias, G.S.; Kitano, E.S.; Pagotto, A.H.; Sant'anna, S.S.; Rocha, M.M.;Zelanis, A.; Serrano, S.M. Individual variability in the venom proteome of juvenile Bothrops jararaca specimens. J. Proteome Res. 2013, 12(10),4585-4598. 12. De Lucca, F.L.; Haddad, A.; Kochva, E.; Rothschild, A.M.; Valeri, V. protein synthesis and morphological changes in the secretory epithelium of the venom gland of Crotalus durissus terrificus at different times after manual extraction of venom. Toxicon 1974,12(4),361-368. 13. Brown, R.S.; Brown, M.B.; Bdolah, A.; Kochva, E. Accumulation of some secretory enzymes in venom glands of Vipera palaestinae. Am. J. Physiol. 1975, 229(6),1675-1679. 14. Yamanouye, N.; Britto, L.R.G.; Carneiro, S.M.; Markus, R.P. Control of venom production and secretion by sympathetic outflow in the snake Bothrops jararaca. J. Exp. Biol. 1997, 200, 2547–2556. 15. Mackessy, S.P.; Baxter, L.M. Bioweapons synthesis and storage:the venom gland of front-fanged snakes. Zool. Anz. 2006, 245,147-159. 16. Kerchove, C.M.; Luna, M.S.; Zablith, M.B.; Lazari, M.F.M.;Smaili, S.S.; Yamanouye. N. α1-adrenoceptors trigger the snake venom production cycle

in

secretory

bisphosphate

cells

hydrolysis

by and

activating ERK

phosphatidylinositol

signaling

pathway.

4,5-

Comp.

Biochem.Physiol A 2008, 150,431-437. 17. Luna, M.S.; Valente, R.H.; Perales, J.; Vieira, M.L.; Yamanouye, N. Activation of Bothrops jararaca snake venom gland and venom production: A proteomic approach. J. Proteomics 2013, 94,460-472.

25 ACS Paragon Plus Environment

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

18. Junqueira-de-Azevedo,

I.L.;

Bastos,

C.M.;

Page 26 of 45

Ho,

P.L.;

Luna,

M.S.;Yamanouye, N.; Casewell, N.R. Venom-related transcripts from Bothrops jararaca tissues provide novel molecular insights into the production and evolution of snake venom. Mol. Biol. Evol. 2015, 32(3), 754-766. 19. Sazima, I. Natural history of the jararaca pitviper, Bothrops jararaca, in southeastern Brazil. In: Biology of Pitvipers. Campbell, J.A.; Brodie, E.D (Eds). Selva Tyler: Texas, 1992, pp. 199-216. 20. Martins, M.; Marques, O. A. V.; Sazima, I. Ecological and phylogenetic correlates of feeding habits in neotropical pitvipers of the genus Bothrops. In: Biology of the vipers. Schuett, G.W; Höggren, M.; Douglas, M. E.; Greene, H. W. (Eds), Eagle Mountain Publishing, LC: Utah, 2002 pp.307-328. 21. Hartmann, P. A.; Hartmann, M. T.; Giasson, l. O. M. Uso do habitat e alimentação em juvenis de Bothrops jararaca (Serpentes: Viperidae) na mata atlântica do sudeste do Brasil. Phyllomedusa 2003, 2(1), 35-41. 22. Cidade, D.A.; Simão, T.A.; Dávila, A.M.; Wagner, G.; Junqueira-deAzevedo, I.L.; Ho, P.L.; Bon, C.; Zingali, R.B.; Albano, R.M. Bothrops jararaca venom gland transcriptome: analysis of the gene expression pattern. Toxicon 2006, 48(4),437-461. 23. Zelanis, A.; Tashima, A. K.; Pinto, A. F.; Leme, A. F.; Stuginski, D. R.; Furtado, M. F.; Sherman, N. E.; Ho, P. L.; Fox J. W.; Serrano S. M. Bothrops jararaca venom proteome rearrangement upon neonate to adult transition. Proteomics 2011, 11 (21), 4218–4228.

26 ACS Paragon Plus Environment

Page 27 of 45

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

24. Zelanis, A.; Serrano, S. M.; Reinhold, V. N. N-glycome profiling of Bothrops jararaca newborn and adult venoms. J. Proteomics 2012, 75(3), 774–782. 25. Zelanis, A.; Andrade-Silva, D.; Rocha, M. M.; Furtado, M. F.; Serrano, S. M.; Junqueira-de-Azevedo, I. L.; Ho, P. L. A transcriptomic view of the proteome variability of newborn and adult Bothrops jararaca snake venoms. PLoS Negl. Trop. Dis. 2012, 6 (3), e1554. 26. Zelanis, A.;Menezes, M.C.; Kitano, E.S.; Liberato, T.; Tashima, A.K.; Pinto, A.F.M.; Sherman, N.E.;Ho, P.L.; Fox, J.W.; Serrano, S.M.T. Proteomic identification of gender molecular markers in Bothrops jararaca venom. J. Proteomics 2016,139,26-37. 27. SINAN/SVS,

Brazilian

Ministry

of

Health,

http://tabnet.datasus.gov.br/cgi/deftohtm.exe?sinannet/cnv/animaisbr.def ;access June, 19th 2016. 28. Sazima, I. Caudal luring in two neotropical pitvipers, Bothrops jararaca and B. jararacussu. Copeia 1991,1,245-248. 29. Luna,

M.S.;

Hortencio,

T.M.A.;

Ferreira,

Z.S.;

Yamanouye,

N.

Sympathetic outflow activates the venom gland of the snake Bothrops jararaca Bothrops jararaca by regulating the activation of transcription factors and the synthesis of vbenom gland proteins. J. Exp. Biol. 2009, 212,1535-1543. 30. Pimenta, D. C.; Prezoto, B. C.; Konno, K.; Melo, R. L.; Furtado, M. F.; Camargo A. C.; Serrano, S. M. Mass spectrometric analysis of the individual variability of Bothrops jararaca venom peptide fraction.

27 ACS Paragon Plus Environment

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 45

Evidence for sex-based variation among the bradykinin-potentiating peptides. Rapid Commun. Mass Spectrom. 2007, 21(6), 1034–1042. 31. Laemmli, U.K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970, 227, 680-685. 32. Paes-Leme, A.F.; Kitano, E.S.; Furtado, M.F.; Valente, R.H.; Camargo, A.C.; Fox, J.W.; Serrano, S.M.T. Analysis of the subproteomes of proteinases and heparin-binding toxins of eight Bothrops venoms. Proteomics 2009 9(3), 733-745. 33. Silva, V.X.; Rodrigues, M.T. Taxonomic revision of the Bothrops neuwiedi complex (Serpentes, Viperidae) with description of a new species. Phyllomedusa 2008, 7(1):45-90. 34. Kleifeld, O.; Doucet, A.; Prudova, A.; auf dem Keller,U.;

Gioia, M.;

Kizhakkedathu, J.N.; Overall, C.M. Identifying and quantifying proteolytic events and the natural N-terminome by terminal amine isotopic labeling of substrates, Nat. Protoc. 2011, 6 (10), 1578–1611. 35. Boersema, P.J.; Raijmakers, R.; Lemeer, S.; Mohammed, S.; Heck, A.J. Multiplex peptide stable isotope dimethyl labeling for quantitative proteomics. Nat. Protoc. 2009, 4(4),484-494. 36. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26 (12),1367–1372. 37. Vizcaíno J.A.; Csordas, A.; del-Toro, N.; Dianes, J.A.; Griss, J.; Lavidas, I.; Mayer, G.; Perez-Riverol, Y.; Reisinger, F.; Ternent, T.; Xu, Q.W.; Wang, R.; Hermjakob, H. 2016 update of the PRIDE database and related tools. Nucleic Acids Res. 2016, 44(D1), D447-D456. 28 ACS Paragon Plus Environment

Page 29 of 45

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

38. Bolstad, B.M.; Irizarry, R.A.; Astrand, M.; Speed, T.P. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003,19(2), 185-193. 39. Ihaka, R.; Gentleman, R. R: a language for data analysis and graphics, J. Comput. Graph. Stat. 1996, 5(3), 299–314. 40. Gentleman, R.C. et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5(10), R80. 41. Smyth, G.K. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 2004, 3, 3. 42. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674-3676. 43. Wilkinson, B.; Gilbert, H.F. Protein disulfide isomerase. Biochim. Biophys. Acta. 2004,1699(1-2),35-44. 44. Petrak J.; Ivanek, R.; Toman, O.; Cmejla, R.; Cmejlova, J.; Vyoral, D.; Zivny, J.; Vulpe, C.D. Déjà vu in proteomics. Hit parade of repeatedly identified differentially expressed proteins. Proteomics 2008, 9,17441749. 45. Portes-Junior, J.A.; Yamanouye, P.S.; Sant'Anna,

S.S.; Nogueira,

N.; Carneiro, F.C.; Junqueira,

S.M.; Knittel, M.; Magalhães

G.S.; Domont, G.B.; Moura-da-Silva, A.M. Unraveling the processing and activation of snake venom metalloproteinases. J. Proteome Res. 2014,13(7),3338-3348.

29 ACS Paragon Plus Environment

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 45

46. Serrano, S.M.; Oliveira, A.K.; Menezes, M.C.; Zelanis, A. The proteinaserich proteome of Bothrops jararaca venom. J.Toxicol. Toxin. Rev., 2014, 1-16. 47. Maier, T.; Güel, M.; Serrano, L. Correlation of mRNA and protein in complex biological samples. FEBS Lett ., 2009,584(24),3966-3973. 48. Schwannhäusser, B.; Busse, D.; Li, N.; Dittmar, G.; Schuchhardt, J.; Wolf, J.; Chen, W.; Selbach, M. Global quantification of mammalian gene expression control. Nature 2011, 473,337-342. 49. Vissink, A.; Spijkevert, F.K.;Van Nieuw Amerongen, A. Aging and saliva: a review of the literature. Spec. Care Dentist.1996,16(3),95-103. 50. Hiratsuka, K.; Kamino, Y.; Nagata, T.; Takahashi, Y.; Asai, S.; Ishikawa, K.; Abiko, Y. Microarray analysis of gene expression changes in aging in mouse submandibular gland. J. Dent. Res.2002, 81(10),679-682. 51. Srivastava, A.; Wang, J.; Zhou, H.; Melvin, J.E.; Wong, D.T. Age and gender related differences in human parotid gland gene expression. Arch. Oral Biol. 2008,53(11),1058,1070. 52. Ma, X.M.; Blenis, J. Molecular mechanisms of mTOR-mediated translational control. Nat. Rev. Mol. Cell. Biol. 2009,10(5), 307-318. 53. Yen,

H.C.;

Xu,

Q.;

Chou,

D.M.;

Zhao,

Global protein stability profiling in mammalian cells.

Z.;

Elledge,

Science

S.J. 2008,

322(5903),918-923. 54. Toyama, B.H.; Hetzer, M.W. Protein homeostasis: live long, won't prosper. Nat. Rev. Mol. Cell Biol. 2013, 14, 55-61. 55. Bourgoin-Voillard, S.; Goron, A.; Seve, M.; Moinard, C. Regulation of the proteome by amino acids. Proteomics 2016,16,831-846.

30 ACS Paragon Plus Environment

Page 31 of 45

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

56. Martens, L.; Vandekerckhove, J.; Gevaert, K. DBToolkit: processing protein databases for peptide-centric proteomics. Bioinformatics 2005, 21,3584-3585. 57. Colaert, N.; Helsens, K.; Martens, L.; Vandekerckhove, J.; Gevaert, K. Improved visualization of protein consensus sequences by iceLogo. Nat Methods 2009,6(11),786-7.

Figure legends

Scheme 1. Schematic overview of the analysis of B. jararaca venom gland proteome.

Figure 1. Overall composition of B. jararaca venom gland proteome (A) and biological annotation of cellular proteins into functional categories (B). The diversity of Post Translational Modification (PTM)/processing category (C).

Figure 2. B. jararaca venom and venom gland protein profiles. (A) SDS-PAGE (12% SDS-polyacrylamide gel) profile of B. jararaca venoms and venom gland extracts (25 µg of proteins) under reducing conditions. Mobility of molecular mass markers is indicated on the left. (B) Western blot analysis of B. jararaca venoms (5 µg) and venom gland extracts (25 µg). After SDS-PAGE separation on 12% SDS-polyacrylamide gels proteins were transferred to a nitrocellulose membrane and immunostained with anti-bothropic antivenom.

31 ACS Paragon Plus Environment

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 45

Figure 3. Volcano plots showing the statistically differentially expressed cellular proteins (t-test; adjusted p-value ≤ 0.05 and log2 fold change 1) in adult and newborn venom glands (A) and adult male and female venom glands (B). Colored dots refer to statistically differentially expressed proteins in adult (blue) or newborn (red) and adult male (blue) and female (red) cellular proteins.

Figure 4. The heatmaps show the normalized (z-score) expression values for differentially expressed proteins (left) together with their functional annotation presented as the percentage of identifications matching specific functional categories (right). Newborn vs. adult venom gland proteins (A and B) and adult male vs. female cellular proteins (C and D). The expression values were measured across three technical replicates. Median expression values were clustered as described in Materials and Methods.

Figure 5. The overall composition of B. jararaca venom proteome as inferred from the analysis of the venom gland proteome (A). (B) Sequence logos of the set of the identified peptides showing the canonical sequences ‘MCGVT’, and ‘SEL’, part of the cysteine-switch and propeptide regions in SVMPs and SVSPs, respectively. Sequence logos were generated with IceLogo57 .

Figure 6. Volcano plot showing the statistically differentially expressed toxins (ttest; adjusted p-value ≤ 0.05 and log2 fold change 1) in adult and newborn venom glands (A) and adult male and female venom glands (B). Differentially expressed proteins in newborn (C) and adult (D) and male (E) and female (F) are presented as the percentages matching to specific toxin classes.

32 ACS Paragon Plus Environment

Page 33 of 45

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

Colored dots refer to statistically differentially expressed proteins in adult (red) or newborn (green) and adult male (blue) and female (pink) venom proteins.

Figure 7. Variation in the log2 median of newborn/adult ratio (A) and male/female ratio found for cellular proteins and toxins (C and D), respectively. The dashed red lines correspond to the standard deviation related to the measures.

33 ACS Paragon Plus Environment

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 45

For TOC only

34 ACS Paragon Plus Environment

Page 35 of 45

Journal of Proteome Research

Dynamic proteome rearrangement in Bothrops jararaca venom gland 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Venom gland 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 Bothrops jararaca 41 (Photo courtesy of André Zelanis, Copyright 2016) 42

Ontogenetic changes

Biological annotation

Sexual dimorphism Proteome diversity

Quantitative analysis

Intraspecific venom variation Venom gland proteome

ACS Paragon Plus Environment

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

A

B

Page 36 of 45

Functional categories 3%

4%

4% Degradation

General metabolism

16%

Inhibitor

18% PTM/processing

47%

Structural

13%

84%

Transcription/translation Transport/sorting

1% 9% Toxins

Secretion

Cellular proteins ACS Paragon Plus Figure 1 Environment

Vesicular trafficking

1%

Page 37 of 45

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

Journal of Proteome Research

PTM/processing

C

C-terminal processing

2% 1% 2%

Folding

6%

Glycation 17%

Glycosylation

2%

70%

N-terminal modification N-terminal processing (secretion) Protein processing

Figure 1

ACS Paragon Plus Environment

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

A

Page 38 of 45

B

kDa

kDa

170 130 100 70

102 76 52

55

38 31

40 35

24

25

17 15

12

Venoms

Venom gland extracts

Venoms

ACS Paragon Plus Environment

Figure 2

Venom gland extracts

Page 39 of 45 Up Adult

A

B

V8NPI7

T1E4R4

U3FBW7

V8NPI7

3 2 V8NSS7

-log10(adj.Pvalue)

V8PA09

-log10(adj.Pvalue)

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

Up Female

Up Newborn Journal of Proteome Research Up Male

2

1.5

1

1

0.5

0

0 -4

-2

0

log2(fold change)

2

4

-2

ACS Paragon Plus Environment

Figure 3

-1

0

log2(fold change)

1

2

3

Journal of Proteome Research

Page 40 of 45

B

UniProt #

Adult

Newborn

Inhibitor

2 1 0 -1

-2

Functional category

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

Adult

Expression value (z-score)

A

Newborn

Transcription/translation

PTM/processing

Structural

General metabolism

0

ACS Paragon Plus Environment

Figure 4

20

40 60 Percentage

80

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

Male

Journal of Proteome Research

D

UniProt #

Female

Male

Inhibitor

2 1 0 -1

Functional category

C

Female

Expression value (z-score)

Page 41 of 45

PTM/processing

Structural

-2

Transcription/translation

General metabolism

0 ACS Paragon Plus Environment

Figure 4

20

40 60 Percentage

80

SVMPs Cysteine-switch

Journal of Proteome Research

B

A

4.0

svVEGF SVSP Bits

SVMP PLB

2.0

PLAI PLA2 NGF LAAO

1

2

3

DPP

4

5

6

7

8

Amino acid residues

DIEST

SVSPs Propeptide region

4.0

CVF

CTL CRISP BPP

Bits

Toxin class

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

Page 42 of 45

5NUCL

2.0

AP 0

10

20 30 Percentage

40

50

ACS Paragon Plus Environment

Figure 5

1

2

3

4

5

6

Amino acid residues

7

8

Page 43 of 45

Up Adult

Up Newborn

B

Up Male

Up Female

-log10(adj.Pvalue)

-log10(adj.Pvalue)

A

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

Journal of Proteome Research

log2(fold change)

log2(fold change) ACS Paragon Plus Environment

Figure 6

Newborn

C 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

Journal of Proteome Research

SVSP 3%

AP CRISP 3% 3%

Page 44 of 45

Male

E CTL 9% DPP 3%

DIEST 3%

SVSP 18%

PLA2 9%

SVMP 64%

NGF 9%

CTL 14%

LAAO 9% PLA2 9%

SVMP 41%

PLB 3%

Adult

D SVSP 36%

Female

F

VF 8%

SVSP 20%

CTL 36%

SVMP 10%

CTL 36%

SVMP 32% LAAO 5% PLAI 5%

PLA2 8%

PLAI 4%

ACS Paragon Plus Environment

Figure 6

Page 45 of 45

Journal of Proteome Research

A

log2(ratio adult/newborn)

log2(ratio adult/newborn)

8 4 0 0

200

-4

400

Toxins

600

Quantified proteins

-8

8 4 0 0

-8

5

0 0

50

100

150

-5

Quantified proteins

10

D

-4

10

-10

100

200

300

Quantified proteins

400

500

600

log2 (ratio male/female)

B

C

Cellular proteins

log2 ratio (male/female)

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

5 0 0

40

60

-5

Quantified proteins -10

Figure 7

20

ACS Paragon Plus Environment

80

100

120