Shotgun Proteomics Analysis Discards Alkali ... - ACS Publications

31 May 2017 - Immunology, Faculty of Biology, University of Vigo, E-36200 Vigo, Galicia, Spain. •S Supporting Information. ABSTRACT: Vitellogenin, t...
0 downloads 0 Views 1MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Shotgun proteomics analysis discards alkali labile phosphate (ALP) as a reliable method to assess vitellogenin levels in Mytilus galloprovincialis Paula Sanchez-Marin, Laura E. Fernández-González, Leonardo Mantilla-Aldana, Angel P. Diz, and Ricardo Beiras Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 31 May 2017 Downloaded from http://pubs.acs.org on June 2, 2017

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.

Environmental Science & Technology 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 34

Environmental Science & Technology

1

Shotgun proteomics analysis discards alkali

2

labile phosphate (ALP) as a reliable method to

3

assess vitellogenin levels in Mytilus

4

galloprovincialis

5

6

Paula Sánchez-Marín1,2,*, Laura E. Fernández-González1,2,3, Leonardo Mantilla-

7

Aldana1,2, Angel P. Diz1,3,# and Ricardo Beiras1,2, #.

8

9

1

Toralla Marine Science Station (ECIMAT), University of Vigo, Illa de Toralla, E-

10

36331, Vigo, Galicia, Spain.

11

2

12

of Vigo, E-36200, Vigo, Galicia, Spain.

13

3

14

University of Vigo, E-36200, Vigo, Galicia, Spain.

15

#

16

*Corresponding author phone and e-mail: (+34)986818780; [email protected].

Department of Ecology and Animal Biology, Faculty of Marine Sciences, University

Department of Biochemistry, Genetics and Immunology, Faculty of Biology,

Authors with equal contribution.

17

18

1 ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 34

19

ABSTRACT

20

Vitellogenin, the egg yolk precursor, is a well known biomarker of endocrine disruption

21

in oviparous vertebrates. In invertebrates, such as bivalves, it has been used in the last

22

10 years for the same purpose, despite the limited knowledge of invertebrate

23

endocrinology. In bivalves, vitellogenin levels are usually estimated using an indirect

24

technique, alkali labile phosphate (ALP), that assumes that vitellogenin is the most

25

abundant phosphorylated protein in the analyzed tissue. In this study, we applied

26

shotgun proteomics for the identification and quantification of vitellogenin in marine

27

mussel gonads and compared the results with those obtained with the ALP method. The

28

proteomic analysis revealed that vitellogenin is only detected in female gonads with

29

expression levels that are rather variable among female mussels at different stages of

30

gonad development. ALP analysis, on the contrary, detected similar amounts of

31

phosphorylated proteins regardless sex or gonad development stage. These results show

32

evidence that the ALP method is not providing reliable information about Vtg levels, at

33

least in marine mussel gonads. ALP is not a good proxy to assess Vtg levels in marine

34

mussels and careful verification of the adequacy of the procedure should be done before

35

ALP is further assumed as a proxy of Vtg in other bivalve mollusks.

36

37

KEYWORDS

38

Mytilus, endocrine disruption, biomarker, proteomics, mass spectrometry analysis,

39

label-free quantification

40

2 ACS Paragon Plus Environment

Page 3 of 34

Environmental Science & Technology

41

TOC/ABSTRACT ART

42 43

3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 34

44

INTRODUCTION

45

In the marine environment, the blue mussel (including species of the Mytilus edulis

46

complex) is widely used as a sentinel species in pollution biomonitoring because it is

47

sessile, ubiquitous, filters high volumes of water and accumulates pollutants in its

48

tissues. Its use in pollution monitoring studies includes not only the measurement of the

49

bioaccumulated pollutants,1-2 but also the analysis of several biochemical biomarkers of

50

exposure or effects, such as different antioxidant enzymatic activities,3 or physiological

51

energetic endpoints.4 More recently, a biomarker of endocrine disruption has also been

52

added to this list: this is alkali labile phosphate (ALP), an easy and inexpensive

53

alternative for the indirect quantification of vitellogenin (Vtg), the egg-yolk precursor

54

(reviewed in5-6). In oviparous vertebrates, Vtg is synthesized in the liver in response to

55

estrogens in maturing females, from where it is transported to the gonads through blood.

56

Vtg is normally undetectable in the plasma of males and juveniles, which makes its

57

anomalous presence in males an excellent biomarker of estrogenicity and feminization

58

in fishes.7-9 Vtg is usually measured in fish serum using an enzyme-linked

59

immunosorbent assay (ELISA), which has the limitation of being species specific.

60

When specific antibodies are not available for the test species, the ALP method has

61

been proved to be useful for the estimation of Vtg levels in fish serum,10 under the

62

assumption that Vtg is the most abundant phosphorylated protein in this fluid.11

63

The ALP assay was also modified for its use in mollusks, and applied in exposure

64

studies that showed increases in ALP levels in females of freshwater mussels and clams

65

that were injected 17β-estradiol or estrogen mimics such as nonylphenol. 12-13 Since

66

then, the ALP assay has been widely applied in several mollusk species, including blue

67

mussels.5-6 Concerning marine mussels, different studies using ALP as a proxy for

68

measuring Vtg levels have shown contrasting results. One study showed increases in 4 ACS Paragon Plus Environment

Page 5 of 34

Environmental Science & Technology

69

ALP levels in male and female mussels upon exposure to nonylphenol,14 while other

70

studies failed in detecting changes in ALP upon exposure to a related chemical,

71

bisphenol A.15-16 Other studies with the Mytilus family have used ALP in the evaluation

72

of the effects of xenobiotics with different modes of action, such as hydrocarbons or

73

pharmaceuticals. In these studies, increases or decreases in ALP as a response to those

74

chemicals were also observed, and the influence of sex in the response was very

75

variable, being in some cases only male or female specific, while in other cases both

76

sexes were affected. 17-19 Surprisingly, no attempt has been made to relate ALP levels

77

with exposure to natural or synthetic estrogens in marine mussels, so it is not confirmed

78

if changes in ALP levels can be attributed to estrogenic or androgenic effects of

79

contaminants. An alternative technique to study Vtg induction showed that Vtg mRNA

80

expression increased several orders of magnitude in both male and female mussels upon

81

exposure to natural and synthetic estrogens only when the mussels were at early

82

gametogenesis stage.20 In situ hybridization has shown that Vtg is only detected in the

83

female gonad, suggesting that this is the primary site of both synthesis and

84

accumulation. 21 However, since ALP can be measured in any tissue, ALP levels have

85

been reported in haemolymph, 22 digestive gland 14 and both unexposed male and

86

female gonads 23 in Mytilus. Recent studies have shown that the most abundant

87

phosphorylated protein present in mussel haemolymph –and reported in ALP

88

measurements– was in fact extrapalial protein precursor, and not Vtg.24 Therefore, the

89

usefulness of ALP for the detection of Vtg in mollusks has been questioned, and more

90

reliable methods are needed before this protein can be used as a biomarker of endocrine

91

disruption in marine mussels.

92

The present study tested whether results from ALP method are a good proxy for Vtg

93

levels in mussel gonads, the primary site of Vtg production and accumulation, using a 5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 34

94

shotgun label-free proteomics approach. To this aim, several male and female mussel

95

gonads were analyzed for Vtg protein identification and quantification, and results

96

discussed in terms of the suitability of this protein as a valid biomarker for the detection

97

of effects caused by estrogenic endocrine disruptors in marine mussels.

98

99

EXPERIMENTAL

100

Sampling and histological analysis

101

Mussels (Mytilus galloprovincialis) between 40 and 50 mm length were collected in

102

unpolluted areas in the Galician Rías (NW Iberian Peninsula) in two different seasons,

103

representing different stages of the gametogenesis cycle: December 2014, when mussels

104

are mainly in the early gametogenesis stage, and in the pre-spawning season, April

105

2015, when most mussels are fully mature.25

106

Once in the laboratory, mussels were opened with a scalpel and the mantle-gonad was

107

excised and split in three (December) or two (April) fragments. Two of the gonad

108

fragments were introduced in microtubes, quick-frozen in liquid nitrogen and conserved

109

at -80 ºC until further processing for ALP and proteomic analysis. The third gonad

110

fragment was used for determination of sex and gonad gametogenesis stage after

111

standard histological processing using haematoxylin-eosin stain. For individuals

112

collected in April, sex and gametogenesis stage were determined by direct observation

113

(size of the gonad, color and texture) and confirmed by observation of a gonad smear

114

under the microscope.

115

ALP analysis

6 ACS Paragon Plus Environment

Page 7 of 34

Environmental Science & Technology

116

ALP extraction was performed following a method modified from Gagné et al. 26 Tissue

117

samples were homogenized on ice with micropestles in 1:5 volume of cold (4 ºC)

118

homogenization buffer (25 mM Hepes, 125 mM NaCl, 1 mM dithiothreitol and 1 mM

119

EDTA; pH 7.4), centrifuged (12,000 g, 30 min, 2 ºC) and the supernatant collected. An

120

aliquot was taken for total protein determination using the Bradford method27 adapted to

121

microplate, using bovine serum albumin (BSA) as standard. Another aliquot was

122

adjusted to 35% acetone, let on ice during 10 min, and centrifuged (10,000 g, 5 min, 2

123

ºC). The pellet was dissolved in 1 M NaOH at 60 ºC during 30 min to hydrolyze the

124

phosphate groups. Before the analysis of inorganic phosphate, a last step was

125

introduced as a modification of the Gagné et al. procedure in order to eliminate turbidity

126

interferences (see text S1 in supporting information). Trichloroacetic acid (TCA) 6.1 N

127

was added to the samples, in a dose of 25 µL per 100 µL of sample, and samples were

128

centrifuged (12,000 g, 15 min). Analysis of released inorganic phosphate was

129

performed in the supernatant. This last step prevents the formation of turbidity that

130

interferes with the spectrophotometric analysis of phosphate if TCA is directly added to

131

the sample in the spectrophotometric well (see text S1 and Figure S1).

132

Measurement of phosphate was performed using a spectrophotometric method. In

133

preliminary tests (see text S1), the phosphomolybdenum method 28 was used when the

134

concentration of phosphate in the sample was between 5-200 mg/L and for lower

135

concentrations, between 0.2 - 2 mg/L PO43- , the malachite green method29 was used

136

instead. The malachite green method was the preferred method afterwards due to the

137

low phosphate concentration of most samples.

138

Phosvitin (from egg yolk, Sigma) (8-10% P) was used as internal standard for quality

139

control. The recovery of PO43- from phosvitin was 96% assuming a P content of 9%.

140

Details on the preliminary tests, quality control and comparison with the results 7 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 34

141

obtained with the original method –without turbidity removal– are given in the

142

supporting information (Text S1).

143

Proteomic analysis

144

Protein extraction and quantification. Protein extraction was performed as previously

145

described,30 by sonication on ice, in 1:5 vol of lysis buffer (7 M urea, 2 M thiourea and

146

4% CHAPS), followed by centrifugation (20,000 g, 30 min, 4 ºC) and supernatant

147

storing at -80ºC until electrophoresis. Total protein concentration was measured with

148

the Bradford method adapted to microplate, using a similar matrix (10% lysis buffer) in

149

both BSA standards and samples.

150

1-D gel electrophoresis. In a first analysis, 20 µg of protein from a mature male and a

151

mature female gonad were separated according to their molecular weight by 12% SDS-

152

PAGE electrophoresis, and stained using OrioleTM fluorescent gel stain (Bio-Rad

153

Laboratories, USA). 1-D gels were visualized on a blue-light DarkReader (Clare

154

Chemical Research, CO, USA). EZ-Run Protein Marker (Fisher BioReagents, NJ, USA)

155

was used as molecular weight marker. The gel was fractionated in 7 sections according

156

to the protein's molecular weight (MW) (>116 kDa, 116-66.2 kDa, 66.2-45 kDa, 45-35

157

kDa, 35-25 kDa, 25-18.4 kDa, and 116 kDa fraction,

160

we proceed to Vtg quantification using a label-free MS analysis approach in different

161

samples focused only on the highest MW gel fraction.

162

Tryptic digestion. Gel sections were cut in small pieces, washed with ultrapure water

163

and acetonitrile, and proteins reduced by treatment with 10 mM dithiothreitol and

164

alkylated with 100 mM iodoacetamide. Samples were then dehydrated in acetonitrile,

8 ACS Paragon Plus Environment

Page 9 of 34

Environmental Science & Technology

165

rehydrated in 25mM ammonium bicarbonate solution followed by addition of

166

acetonitrile, and dried at room temperature. Forty µl of sequencing-grade modified

167

trypsin (12.5 ng µl-1; Promega, Madison, WI) in 25 mM ammonium bicarbonate

168

solution was added to each sample, kept on ice during 45 min, followed by an

169

incubation at 37 ºC overnight. Peptides were recovered and transferred to a clean

170

microtube using 0.1% formic acid/50% acetonitrile and dried under speed-vacuum

171

(Concentrator plus, Eppendorf, Hamburg, Germany).

172

Mass spectrometry (LC-MS/MS) analysis. Tryptic digests were acidified with 0.5%

173

formic acid and analysed in a LTQ-Orbitrap Elite mass spectrometer coupled to a

174

Proxeon EASY-nLC 1000 UHPLC system (Thermo Fisher, San Jose, CA). For protein

175

identification, MS/MS spectra were searched using PEAKS Studio v7.0 software

176

(Bioinformatics Solutions Inc., Waterloo, Canada) against two customized protein

177

databases: one containing 411,509 non-redundant protein sequences of mollusks

178

retrieved from NCBI in 11th May of 2016, and another one resulting from the six-frame

179

translation of 151,320 sequences from available M. galloprovincialis transcriptomes

180

(RNA-seq data) from different tissues, including the mantle-gonad. 31 This database

181

contains two sequences that were annotated as vitellogenin by sequence homology with

182

other species (CL1077.Contig3_All_5 and Unigene21_All_6). The second one

183

(Unigene21_All_6), hereafter referred to as "the M. galloprovincialis Vtg sequence",

184

was used for the identification and quantification of Vtg and included also in the

185

database retrieved from NCBI. A list of common contaminants was also added to both

186

databases, and a decoy sequence database was included in the analysis to calculate the

187

false discovery rate (FDR). All cysteines were considered as fully

188

carbamidomethylated, while methionine oxidation was searched as a variable

189

modification. Positive protein identifications were only accepted when the number of 9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 34

190

matched peptide sequences > 2, unique peptides > 1, peptide spectrum matches (PSMs)

191

FDR < 0.1%, and protein identification PEAKS score > 20.32

192

Label-free quantification of vitellogenin. The relative abundance of Vtg in the different

193

samples was estimated using the normalized spectral abundance factor (NSAF),33 using

194

the spectral counts derived from protein identifications against the customized protein

195

database generated from RNA-seq data (see above). This was done using PatternLab for

196

proteomics 4.0.34 Positive protein identifications were only accepted when the number

197

of matched peptide sequences > 2, unique peptides > 1, PSMs FDR < 3%, peptide FDR

198

< 2% and protein FDR < 1 %. Calculation of fold change was performed by dividing the

199

mean NSAF value of Vtg in each treatment by the mean NSAF value of Vtg in the

200

treatment with the lowest value.

201

Statistics

202

Data normality and homoscedasticity was tested by the Shapiro-Wilk and Levene tests

203

respectively before carrying out an analysis of variance (ANOVA) followed by the post

204

hoc Tukey HSD when multiple treatments were compared. Statistical analysis were

205

performed with IBM SPSS Statistics 23.

206

207

RESULTS AND DISCUSSION

208

Gametogenesis stage

209

The individuals collected in April, hereafter referred to as "mature", presented very

210

thick, uniform gonads with no signs of spawning. The gonads were white or cream-

211

colored in the case of males and sperm poured easily when the gonads were lightly

212

pressed with a scalpel. Female gonads were orange-colored and with a granular texture.

10 ACS Paragon Plus Environment

Page 11 of 34

Environmental Science & Technology

213

The smear observation revealed spherical full grown oocytes in the case of females and

214

mobile sperm in the case of males.

215

Individuals collected in December, hereafter referred to as "immature", presented very

216

thin mantle tissues that were almost transparent. Histological observations revealed the

217

presence of very few follicles embedded in the connective tissue, and presented germ

218

cells at different stages of development. Gametogenesis stage was classified in the

219

histological stage I or II according to classification given in Martinez-Castro and

220

Vazquez.35 Histological sections of selected specimens (representative of these

221

gametogenesis stages) are shown in the supporting information (Table S2).

222

ALP measurements

223

Results of ALP analysis are presented in Table 1. ALP data normalized by total protein

224

was quite similar for all mussel groups, regardless sex or gametogenesis stage, with the

225

only exception of mature females, which showed significantly lower ALP levels than

226

males. These values are not expected under the assumptions that PO43- released in ALP

227

measurements would come mostly from Vtg proteins, since Vtg proteins are supposed

228

to be mainly, if not only, synthesized in female tissues, and should be higher in mature

229

individuals.36 Other studies have also found similar or even higher ALP values in males

230

compared to females samples of different bivalve mollusks,37-39 including blue

231

mussels15, 18, under control conditions. It has been proposed that these unexpected

232

results are explained by over-normalization due to the higher protein content of females

233

at advanced gametogenesis stages, and normalization by tissue weight was proposed to

234

solve this problem.36 ALP values normalized by tissue weight are indeed around twice

235

for mature females compared to mature males (Table 1), and higher for mature than for

236

immature individuals, but still there are not differences between males and females

237

when mussels are not mature. Despite this normalization gives values that are more 11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 34

238

similar to the expected if variation in ALP data would be mainly explained by

239

differences in Vtg levels, it is important to note that they follow exactly the same

240

pattern as total protein contents in mantle tissue, as confirmed by a statistically

241

significant high correlation coefficient among both variables (ALP normalized by tissue

242

weight vs total protein; r = 0.954; p = 0.046). It seems that the higher the amount of

243

protein extracted, the higher the amount of PO43- extracted from those proteins. This

244

phosphate may come from any phosphorylated proteins in the sample, indicating that

245

the basal levels of phosphorylation in gonads might be similar or higher than those of

246

Vtg, which compromises the use of ALP for Vtg determination in these tissues.

247

Other studies have shown higher ALP values (normalized by total protein) in females

248

than in males of blue mussels17 and other bivalves,13 as well as covariations between

249

ALP levels and gametogenesis cycle.40,12 These results are not concordant with others

250

provided in different studies,15, 18, 37-39 nor with those from the present study. This

251

important discrepancy among studies provide evidence that reported ALP levels show

252

low reproducibility due to high technical variation. In the present study, the

253

methodology was carefully examined to assure that ALP values are correct and not

254

overestimated due to sample turbidity (Text S1).

255

Identification of vitellogenin by shotgun proteomics

256

The protein sequence of M. galloprovincialis vitellogenin, consists on a 2828 amino

257

acids length sequence and presents homology with other vitellogenins described so far

258

for mollusks. It should be noted, however, that the last ~200 amino acids are not present

259

in other Vtg sequences, neither detected in LC-MS/MS analysis (see below), indicating

260

that they probably correspond to a 3' UTR (untranslated region). An NCBI-BLASTp

261

search41 reports the highest alignment scores (> 800) and lowest E-values (< 10-141 )

262

with Vtg sequences described for other marine bivalves, such as Pecten maximus, 12 ACS Paragon Plus Environment

Page 13 of 34

Environmental Science & Technology

263

Saccostrea gigas and Mizuhopecten yessoensis, with >55% of coverage and 30 to 40%

264

identity for complete sequences. An almost perfect alignment (99% identity) was also

265

found with the partial sequence of M. galloprovincialis Vtg available at NCBI

266

(accession number AHB63226.1), that only covers a 4% of the total protein length. It

267

also presents 94% identity with a Vtg sequence of 2593 amino acids length obtained

268

from a first release of the M. galloprovincialis genome.42 The alignment of the complete

269

M. galloprovincialis Vtg sequence with that of P. maximus (the one showing the lowest

270

E-value) is presented in Figure 1.

271

LC-MS/MS analysis of proteins extracted from the gonad of a mature female detected

272

Vtg with a coverage of 66% (280 peptides -all unique-, including post-translational

273

modifications, and 3100 peptide spectrum matches) (see Fig. 1 and supporting file 1).

274

Vtg was detected mainly in the >116 kDa fraction of the gel, which is in agreement with

275

the expected molecular weight of the protein, 326 kDa. On the contrary, none of the

276

peptides detected in the male gonad correspond to the M. galloprovincialis Vtg

277

sequence, neither against any other vitellogenin or vitelin-like sequence when the

278

analysis was repeated against the NCBI database for mollusks.

279

Subsequent proteomic analysis were only performed in the >116 KDa fraction of the

280

electrophoresis gel.

281

Label-free quantification of vitellogenin by shotgun proteomics

282

Results of the label-free LC-MS/MS quantification of Vtg are represented in Table 2.

283

We did not find any evidence of Vtg in male gonads. On the contrary, we provide

284

evidence of the presence of Vtg in female gonads in an amount that is largely dependent

285

on the gametogenesis stage of mussel. Mature females presented 8.5 times more Vtg

286

than immature females (in stage II of gonad development). Immature females in stage I 13 ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 34

287

of gonad development (which do not present vitellogenic oocytes)35 did not present

288

enough amount of Vtg to be detected by the methodology used in this study, despite

289

there is evidence of its synthesis in adipogranular and follicle cells of female gonads at

290

early stages of development. 21 According to mean NSAF values, Vtg represents around

291

0.9% of the detected gonad proteins in the >116 KDa fraction for females in early

292

gametogenesis, and around 7.5% for mature females.

293

The fact that Vtg is only detected in female gonads at advanced stages of gametogenesis

294

is in concordance with a study using in situ hybridization of Vtg in Crassostrea gigas,

295

that shows the presence of this protein only in follicle cells of females in late

296

vitellogenesis or ripe state.43 In contrast, similar contents of Vtg were detected by HPLC

297

in males and females of the species Saccostrea glomerata, both unexposed and after

298

exposure to xenoestrogens.44-45 The technique used in the present study is highly

299

sensitive.46 In fact, from 108 to 242 proteins were detected only in the >116 KDa

300

fraction of female gonad samples, disregarding gametogenesis stage, and from 141 to

301

223 in male samples. Despite this, we were unable to detect Vtg neither in males nor in

302

females at the very early stage of gametogenesis (stage I, see Table S2 in supporting

303

information). Therefore, Vtg is either absent or expressed at extremely low levels in the

304

gonads of male individuals that were not exposed to estrogens. It would be interesting

305

in the future to test if Vtg synthesis is induced upon exposure to estrogens, and in case

306

this is confirmed in males, it would enable the use of Vtg as a valid biomarker for

307

endocrine disruption in this species. In fact, it has been already shown that Vtg mRNA

308

expression was induced both in males and females of the closely related species M.

309

edulis upon exposure to 17β-estradiol and synthetic estrogens.20 It is still unknown

310

whether or not this mRNA is effectively translated to proteins in males. A partial

311

sequence of Vtg was also found in an unpublished RNA-seq project from M. edulis and 14 ACS Paragon Plus Environment

Page 15 of 34

Environmental Science & Technology

312

M. galloprovincialis mature male gonads (Romero et al., manuscript in preparation), but

313

the Vtg protein was not found within highly expressed proteins in M. edulis sperm.47

314

Comparison of the indirect method (ALP) with direct Vtg quantification

315

The disagreement among the indirect ALP method and direct Vtg quantification by

316

shotgun proteomics is evident when comparing Fig. 2a and Fig. 2b. Since Vtg was so

317

low in females in stage I of gametogenesis that it could not be detected by LC-MS/MS,

318

only females in stage II were chosen for this comparison, and therefore ALP values

319

change slightly compared to those presented in Table 1. In any case, there were not

320

significant differences in ALP values of females in stage I (7 ± 3 mg PO43- / g of

321

protein; n = 4) or II (4 ± 1 mg PO43- / g of protein; n = 6).

322

The results of both methods are in high disagreement since ALP levels were as high in

323

males as in females, despite Vtg being undetectable in males by LC-MS/MS. Regarding

324

females, mature individuals presented 8.5 fold more Vtg than immature ones, while

325

ALP values were again similar for both gametogenesis stages. ALP normalization by

326

tissue weight does not represent a significant improvement in the results, since high

327

ALP levels are still observed in mature males, that exceed those observed in immature

328

females, despite LC-MS/MS analysis shows the opposite trend (Figure S3, supporting

329

information). This clearly indicates that ALP is not suitable as a proxy of Vtg in mussel

330

gonads.

331

In fish plasma, significant correlations have been found among ALP levels and Vtg

332

measured by direct techniques such as enzyme immunoassays,10, 48 but similar

333

correlations have not yet been shown for mollusks, despite this has been frequently

334

assumed. On the contrary, no correlation was observed between ALP levels and yolk

335

protein levels (measured by ELISA) in haemolymph of the freshwater mussel Unio

15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 34

336

tumidus, and, as occurred for blue mussels, ALP levels did not present differences

337

depending on the sex, while ELISA detected higher yolk protein levels in females.49 It

338

seems clear that although ALP was a suitable technique for fish plasma, where Vtg is

339

the most abundant phosphorylated protein,48, 50 this seems far from being the case for

340

mollusk tissues. Indeed, the ALP method was shown to be unsuitable for the

341

measurement of Vtg-like proteins in gastropods, masked by the background level of

342

phosphorylation of other proteins that fractionate in acetone.51 The origin of ALP in M.

343

galloprovincialis haemolymph was revealed to be extrapallial protein precursor, and not

344

Vtg,24 and the present results reveal that ALP in M. galloprovincialis gonads may well

345

come from different proteins, also extracted with acetone fractionation, with a

346

phosphorylation level that might be higher or similar to that of Vtg, as evidenced by the

347

high ALP values observed in males, and no evidence of Vtg expression on the basis of

348

LC-MS/MS analysis. Protein phosphorylation is indeed a common process initiated by

349

numerous factors,52 so this result is not surprising at all.

350

Not only there might be another phosphorylated proteins in mollusks gonads, but also, it

351

is reasonable to question if mollusks' Vtgs are also highly phosphorylated as those of

352

vertebrates,49 which is the main basis for the application of the ALP methodology.

353

Vertebrate Vtgs are highly phosphorylated especially because of the presence of high

354

serine contents in the phosvitin moiety, but invertebrate Vtgs seem to lack this phosvitin

355

domain.53 This is the case for the oysters Saccostrea glomerata and Crassostrea gigas,54

356

although their serine contents (12-14%) are comparable to that of vertebrate Vtgs.53-54

357

According to Gagné,36 Vtg from Mytilus edulis has about 14.8% of potential

358

phosphorylation sites (serine and threonine) , slightly higher than the mean serine and

359

threonine contents in eukariotic proteins (12.9%).53 Similarly, the M. galloprovincialis

360

Vtg sequence has 9.8% serine and 3.7% threonine in its composition, and the final 16 ACS Paragon Plus Environment

Page 17 of 34

Environmental Science & Technology

361

degree of phosphorylation of the protein –depending on the level of phosphorylation

362

that is actually achieved during post-translational modifications– is predicted to be 11%

363

on the basis of amino acid sequence. 55-56

364

The present results clearly showed that ALP was not suitable as a proxy for Vtg in this

365

species, probably due to a high ratio of background to Vtg phosphorylation. It is highly

366

probable that this is the case for other bivalves too, that may show a similar proteomic

367

composition in their gonads, so careful verification of the adequacy of the procedure

368

should be done before ALP is further assumed as a proxy of Vtg in bivalve mollusks.

369

As already noted by Morthorst et al.,49 in agreement with this study, the assumption that

370

ALP measurements reflected yolk protein levels in bivalves has been the foundation for

371

more than 50 publications, whose conclusions should be carefully revisited.

372

Environmental implications

373

There is a wide variety of chemicals with confirmed or suspected endocrine disruptor

374

ability,57 that are spilled in the environment and can reach coastal areas through

375

effluents of wastewater treatment plants, runoff or deposition. Deleterious effects of

376

xenoestrogens have been reported in recent years in European estuaries, and bivalve

377

molluscs exhibiting intersex, i.e. the simultaneous presence of both testicular and

378

ovarian tissue in the male gonad, are increasing in number and sites affected.58-59 Due to

379

the difficulty in analyzing and monitoring this myriad of chemicals, it is therefore

380

desirable to count with an early-warning biomarker of endocrine disruption in marine

381

biota.

382

Concerning mollusks, very little is known about their endocrinology, and despite it has

383

been frequently assumed that reproductive hormones play similar roles as those

384

described for vertebrates, the topic is under debate,60-61 and it has not been

17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 34

385

unequivocally demonstrated that vertebrate-like hormones such as estrogens play any

386

role in reproduction in mollusks, and it has been also proposed that these processes are

387

regulated by other mechanisms, such as by neuropeptides.62 Some recent studies

388

showed evidence supporting the theory that the estrogen receptor (ER) plays a role in

389

mollusks' reproduction and that ER and Vtg mRNA expression are induced in response

390

to estrogens.20, 54, 63-65 Also, direct Vtg analysis by HPLC showed that its synthesis is

391

induced in oysters in response to estrogens.44-45 More studies are needed to better

392

understand the mechanism of estrogen signaling in mollusks' endocrinology and its role

393

in Vtg induction. In particular, the analysis of Vtg using robust techniques (such as LC-

394

MS/MS) will be very useful in this research process, especially to confirm if the

395

upregulation of Vtg mRNA is reflected at the protein level.66

396

Also, Vtg analysis will be of interest in the quantitative determination of gametogenesis

397

stage in females, or for other applications, given the other non-reproductive functions

398

attributed to Vtg.67-68 If less expensive techniques are preferred for routine analysis of

399

Vtg in marine mussels, the use of other techniques such as immunoassays should be

400

explored, given that ALP is not a reliable alternative.

401 402

Supporting information

403

Details on the standardization and quality control of ALP analysis, including

404

comparison with the results obtained with the original method (without turbidity

405

removal) and full absorbance spectra. A table showing histological sections of gonads

406

representative of stage I and II of male and females mussels. An image file (supporting

407

file 1) obtained from PEAKS Studio v7.0 software showing the coverage detected by

408

LC-MS/MS of the Vtg sequence in a female gonad including PTMs.

18 ACS Paragon Plus Environment

Page 19 of 34

Environmental Science & Technology

409

410

ACKNOWLEDGEMENTS

411

LC-MS/MS analysis were performed in the CACTI and histological analysis were

412

performed in the ECIMAT, both centers belonging to University of Vigo, by Manuel

413

Marcos (LC-MS/MS) and Rosana Rodriguez (histology). Carlos Canchaya provided us

414

with a Vtg sequence obtained from the M. galloprovincialis genome project. Mónica

415

Rodriguez-Romero is also acknowledged for technical assistance. This study was

416

funded by the Spanish Government through projects CTM2013-48194-C3-3-R and

417

PCIN-2015-187-C03-03, and by "Consellería de Educación e Ordenación Universitaria"

418

Xunta de Galicia (Galician Regional Government), with cofunding from the European

419

Regional Development Fund (ERDF). P.S.-M. was supported by the People Programme

420

(Marie Curie Actions) of the European Union's Seventh Framework Programme

421

(FP7/2007-2013) under REA grant agreement n° 600391 and A.P.D. was supported by

422

the Spanish “Ministerio de Economía y Competitividad” (code AGL2014-52062-R),

423

Fondos Feder and Xunta de Galicia (“Grupos de Referencia Competitiva” ED431C

424

2016-037)".

19 ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 34

425

TABLES

426

Table 1. Levels of ALP measured in M. galloprovincialis gonad tissues, normalized

427

by protein or by tissue weight, and total protein extracted per gram of tissue. Mussel's group

ALP (µg PO43- / mg protein)

ALP (µg PO43- / g tissue)

Total protein (mg / g tissue)

Mature females

2.9 ± 0.8 a

221 ± 51 a

78 ± 7

Mature males

6.5 ± 2.9 b

122 ± 53 b

19 ± 2

Immature females

5.2 ± 2.6 ab

52 ± 29 c

11 ± 5

Immature males

5.5 ± 1.9 b

43 ± 14 c

8±2

428

Mean ± SD is represented (n = 10).

429

Values sharing the same superscript letter do not present significant differences among

430

them (HSD Tukey, α = 0.05).

431

20 ACS Paragon Plus Environment

Page 21 of 34

Environmental Science & Technology

432

Table 2. Proteins detected, Vtg spectral counts and NSAF values obtained from

433

LC-MS/MS analysis in the > 116 kDa fraction of proteins extracted from M.

434

galloprovincialis gonad tissues (three biological replicates per treatment).

435

Mussel's group

Nº Proteins detected

Vtg spectral counts

Vtg NSAF

Mature females

108 242 177

434 1182 801

0.0708 0.0781 0.0749

Mature males

162 141 195

0 0 0

n.a. n.a. n.a.

Immature females (stage II)

232 139 111

143 79 4

0.0141 0.0113 0.0006

Immature males

159 211 223

0 0 0

n.a. n.a. n.a.

n.a. = not applicable due to zero values in spectral counts.

436

21 ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 34

437

FIGURE CAPTIONS

438

Figure 1. Alignment of the M. galloprovincialis vitellogenin sequence (obtained from

439

the RNA_seq sequence Unigene21_All_6 )31 with that showing the best identity (lowest

440

E-value) in NCBI-Blast search, Pecten maximus (GenBank: CAQ06469.2). Red lines

441

represent the peptide coverage detected by LC-MS/MS in a M. galloprovincialis mature

442

female gonad sample.

443

444

Figure 2. Comparison of (a) indirect quantification of Vtg using ALP and (b) direct

445

quantification using LC-MS/MS in gonads of Mytilus galloprovincialis females (black

446

bars) or males (white bars) in different stages of gametogenesis: stage II of

447

gametogenesis (Immature) and fully mature individuals (Mature). The absent bars

448

corresponding to males in figure (b) correspond to 0 values (undetected Vtg).

449

22 ACS Paragon Plus Environment

Page 23 of 34

Environmental Science & Technology

450

REFERENCES

451

(1)

concentrations of selected hazardous substances in sediments and biota; 2009.

452 453

OSPAR CEMP assessment report: 2008/2009. Assessment of trends and

(2)

Kimbrough, K. L.; Johnson, W. E.; Lauenstein, G. G.; Christensen, J. D.; Apeti,

454

D. A. An Assessment of Two Decades of Contaminant Monitoring in the

455

Nation’s Coastal Zone; Silver Spring, MD. NOAA Technical Memorandum

456

NOS NCCOS 74. 105 pp.; 2008.

457

(3)

Vidal-Liñán, L.; Bellas, J.; Campillo, J. A.; Beiras, R. Integrated use of

458

antioxidant enzymes in mussels, Mytilus galloprovincialis, for monitoring

459

pollution in highly productive coastal areas of Galicia (NW Spain).

460

Chemosphere 2010, 78, 265-272.

461

(4)

Albentosa, M.; Viñas, L.; Besada, V.; Franco, A.; González-Quijano, A. First

462

measurements of the scope for growth (SFG) in mussels from a large scale

463

survey in the North-Atlantic Spanish coast. Sci. Tot. Environ. 2012, 435-436,

464

430-445.

465

(5)

Porte, C.; Janer, G.; Lorusso, L. C.; Ortiz-Zarragoitia, M.; Cajaraville, M. P.;

466

Fossi, M. C.; Canesi, L. Endocrine disruptors in marine organisms: Approaches

467

and perspectives. Comp. Biochem. Phys. C 2006, 143, 303-315.

468

(6)

Matozzo, V.; Gagné, F.; Marin, M. G.; Ricciardi, F.; Blaise, C. Vitellogenin as a

469

biomarker of exposure to estrogenic compounds in aquatic invertebrates: A

470

review. Environ. Int. 2008, 34, 531-545.

471 472

(7)

Jones, P. D.; De Coen, W. M.; Tremblay, L.; Giesy, J. P. Vitellogenin as a biomarker for environmental estrogens. Water Sci. Technol. 2000, 42, 1-14.

23 ACS Paragon Plus Environment

Environmental Science & Technology

473

(8)

Page 24 of 34

Diniz, M. S.; Peres, I.; Magalhães-Antoine, I.; Falla, J.; Pihan, J. C. Estrogenic

474

effects in crucian carp (Carassius carassius) exposed to treated sewage effluent.

475

Ecotox. Environ. Safe. 2005, 62, 427-435.

476

(9)

Sumpter, J. P.; Jobling, S. Vitellogenesis as a biomarker for estrogenic

477

contamination of the aquatic environment. Environ. Health Persp. 1995, 103,

478

173-178.

479

(10)

Nagler, J. J.; Ruby, S. M.; Idler, D. R.; So, Y. P. Serum phosphoprotein

480

phosphorus and calcium levels as reproductive indicators of vitellogenin in

481

highly vitellogenic mature female and estradiol-injected immature rainbow trout

482

(Salmo gairdneri). Can. J. Zool. 1987, 65, 2421-2425.

483

(11)

Hallgren, P.; Mårtensson, L.; Mathiasson, L. Improved spectrophotometric

484

vitellogenin determination via alkali-labile phosphate in fish plasma–a cost

485

effective approach for assessment of endocrine disruption. Int. J. Environ. An.

486

Ch. 2009, 89, 1023-1042.

487

(12)

Blaise, C.; Gagné, F.; Pellerin, J.; Hansen, P. D. Determination of vitellogenin-

488

like properties in Mya arenaria hemolymph (Saguenay Fjord, Canada): A

489

potential biomarker for endocrine disruption. Environ. Toxicol. 1999, 14, 455-

490

465.

491

(13)

Gagné, F.; Blaise, C.; Salazar, M.; Salazar, S.; Hansen, P. D. Evaluation of

492

estrogenic effects of municipal effluents to the freshwater mussel Elliptio

493

complanata. Comp. Biochem. Phys. C 2001, 128, 213-225.

494

(14)

Ricciardi, F.; Matozzo, V.; Marin, M. G. Effects of 4-nonylphenol exposure in

495

mussels (Mytilus galloprovincialis) and crabs (Carcinus aestuarii) with

496

particular emphasis on vitellogenin induction. Mar. Pollut. Bull. 2008, 57, 365-

497

372.

24 ACS Paragon Plus Environment

Page 25 of 34

Environmental Science & Technology

498

(15)

Aarab, N.; Lemaire-Gony, S.; Unruh, E.; Hansen, P. D.; Larsen, B. K.;

499

Andersen, O. K.; Narbonne, J. F. Preliminary study of responses in mussel

500

(Mytilus

501

tetrabromodiphenyl ether. Aquat. Toxicol. 2006, 78, Supplement, S86-S92.

502

(16)

edilus)

exposed

to

bisphenol

A,

diallyl

phthalate

and

Ortiz-Zarragoitia, M.; Cajaraville, M. P. Biomarkers of exposure and

503

reproduction-related effects in mussels exposed to endocrine disrupters. Arch.

504

Environ. Cont. Tox. 2006, 50, 361-369.

505

(17)

Baussant, T.; Ortiz-Zarragoitia, M.; Cajaraville, M. P.; Bechmann, R. K.; Taban,

506

I. C.; Sanni, S. Effects of chronic exposure to dispersed oil on selected

507

reproductive processes in adult blue mussels (Mytilus edulis) and the

508

consequences for the early life stages of their larvae. Mar. Pollut. Bull. 2011, 62,

509

1437-1445.

510

(18)

Gonzalez-Rey, M.; Bebianno, M. J. Does selective serotonin reuptake inhibitor

511

(SSRI) fluoxetine affects mussel Mytilus galloprovincialis? Environ. Pollut.

512

2013, 173, 200-209.

513

(19)

Gonzalez-Rey, M.; Bebianno, M. J. Does non-steroidal anti-inflammatory

514

(NSAID) ibuprofen induce antioxidant stress and endocrine disruption in mussel

515

Mytilus galloprovincialis? Environ. Toxicol. Phar. 2012, 33, 361-371.

516

(20)

Ciocan, C. M.; Cubero-Leon, E.; Puinean, A. M.; Hill, E. M.; Minier, C.; Osada,

517

M.; Fenlon, K.; Rotchell, J. M. Effects of estrogen exposure in mussels, Mytilus

518

edulis, at different stages of gametogenesis. Environ. Pollut. 2010, 158, 2977-

519

2984.

520 521

(21)

Agnese, M.; Verderame, M.; De Meo, E.; Prisco, M.; Rosati, L.; Limatola, E.; del Gaudio, R.; Aceto, S.; Andreuccetti, P. A network system for vitellogenin

25 ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 34

522

synthesis in the mussel Mytilus galloprovincialis (L.). J. Cell. Physiol. 2013,

523

228, 547-555.

524

(22)

Pampanin, D. M.; Marangon, I.; Volpato, E.; Campesan, G.; Nasci, C. Stress

525

biomarkers

526

galloprovincialis) collected in the urban area of Venice (Venice Lagoon, Italy).

527

Environ. Pollut. 2005, 136, 103-107.

528

(23)

and

alkali-labile

phosphate

level

in

mussels

(Mytilus

Gonzalez-Rey, M.; Bebianno, M. J. Effects of non-steroidal anti-inflammatory

529

drug (NSAID) diclofenac exposure in mussel Mytilus galloprovincialis. Aquat.

530

Toxicol. 2014, 148, 818-831.

531

(24)

Oliveri, C.; Peric, L.; Sforzini, S.; Banni, M.; Viarengo, A.; Cavaletto, M.;

532

Marsano, F. Biochemical and proteomic characterisation of haemolymph serum

533

reveals the origin of the alkali-labile phosphate (ALP) in mussel (Mytilus

534

galloprovincialis). Comp. Biochem. Phys. D 2014, 11, 29-36.

535

(25)

the bays of Galicia (N.W. Spain). Aquaculture 1995, 130, 269-277.

536 537

Villalba, A. Gametogenic cycle of cultured mussel, Mytilus galloprovincialis, in

(26)

Gagné, F.; Blaise, C.; Pellerin, J.; Pelletier, E.; Douville, M.; Gauthier-Clerc, S.;

538

Viglino, L. Sex alteration in soft-shell clams (Mya arenaria) in an intertidal

539

zone of the Saint Lawrence River (Quebec, Canada). Comp. Biochem. Phys. C

540

2003, 134, 189-198.

541

(27)

Bradford, M. M. A rapid and sensitive method for the quantitation of microgram

542

quantities of protein utilizing the principle of protein-dye binding. Anal.

543

Biochem. 1976, 72, 248-254.

544

(28)

Stanton, M. G. Colorimetric determination of inorganic phosphate in the

545

presence of biological material and adenosine triphosphate. Anal. Biochem.

546

1968, 22, 27-34.

26 ACS Paragon Plus Environment

Page 27 of 34

Environmental Science & Technology

547

(29)

Hallgren, P.; Mårtensson, L.; Mathiasson, L. A new spectrophotometric method

548

for improved indirect measurement of low levels of vitellogenin using malachite

549

green. Int. J. Environ. An. Ch. 2011, 92, 894-908.

550

(30)

Diz, A. P.; Dudley, E.; Cogswell, A.; Macdonald, B. W.; Kenchington, E. L. R.;

551

Zouros, E.; Skibinski, D. O. F. Proteomic analysis of eggs from mytilus edulis

552

females differing in mitochondrial dna transmission mode. Mol. Cell.

553

Proteomics 2013, 12, 3068-3080.

554

(31)

Moreira, R.; Pereiro, P.; Canchaya, C.; Posada, D.; Figueras, A.; Novoa, B.

555

RNA-Seq in Mytilus galloprovincialis: Comparative transcriptomics and

556

expression profiles among different tissues. BMC Genomics 2015, 16, doi:

557

10.1186/s12864-015-1817-5.

558

(32)

Atasaral-Şahin, Ş.; Romero, M. R.; Cueto, R.; González-Lavín, N.; Marcos, M.;

559

Diz, A. P. Subtle tissue and sex-dependent proteome variation in mussel

560

(Mytilus galloprovincialis) populations of the Galician coast (NW Spain) raised

561

in a common environment. Proteomics 2015, 15, 3993-4006.

562

(33)

Zybailov, B.; Mosley, A. L.; Sardiu, M. E.; Coleman, M. K.; Florens, L.;

563

Washburn, M. P. Statistical analysis of membrane proteome expression changes

564

in Saccharomyces cerevisiae. J. Proteome Res. 2006, 5, 2339-2347.

565

(34)

Carvalho, P. C.; Lima, D. B.; Leprevost, F. V.; Santos, M. D. M.; Fischer, J. S.

566

G.; Aquino, P. F.; Moresco, J. J.; Yates, J. R.; Barbosa, V. C. Integrated analysis

567

of shotgun proteomic data with PatternLab for proteomics 4.0. Nature Protocols

568

2016, 11, 102-117.

569

(35)

Martínez-Castro, C.; Vázquez, E. Reproductive cycle of the cockle

570

Cerastoderma edule (Linnaeus 1758) in the Ría de Vigo (Galicia, Northwest

571

Spain). J. Shellfish Res. 2012, 31, 757-767.

27 ACS Paragon Plus Environment

Environmental Science & Technology

572

(36)

Gagné, F. Neuroendocrine Disruption. In Biochemical Ecotoxicology: Principles and Methods; Academic Press: London, U.K. 2014; pp 145-170.

573 574

Page 28 of 34

(37)

Quinn, B.; Gagné, F.; Costello, M.; McKenzie, C.; Wilson, J.; Mothersill, C.

575

The endocrine disrupting effect of municipal effluent on the zebra mussel

576

(Dreissena polymorpha). Aquat. Toxicol. 2004, 66, 279-292.

577

(38)

Saavedra, L.; Leonardi, M.; Morin, V.; Quiñones, R. A. Induction of

578

vitellogenin-like lipoproteins in the mussel aulacomya ater under exposure to

579

17β-estradiol. Rev. Biol. Mar. Oceanog. 2012, 47, 429-438.

580

(39)

in the clam Tapes philippinarum? Environ. Res. 2005, 97, 43-49.

581 582

Matozzo, V.; Marin, M. G. Can 4-nonylphenol induce vitellogenin-like proteins

(40)

Ortiz-Zarragoitia, M.; Cajaraville, M. P. Intersex and oocyte atresia in a mussel

583

population from the Biosphere's Reserve of Urdaibai (Bay of Biscay). Ecotox.

584

Environ. Safe. 2010, 73, 693-701.

585

(41)

Altschul, S. F.; Madden, T. L.; Schäffer, A. A.; Zhang, J.; Zhang, Z.; Miller, W.;

586

Lipman, D. J. Gapped BLAST and PSI-BLAST: A new generation of protein

587

database search programs. Nucleic Acids Res. 1997, 25, 3389-3402.

588

(42)

Murgarella, M.; Puiu, D.; Novoa, B.; Figueras, A.; Posada, D.; Canchaya, C. A

589

first

590

galloprovincialis. PLoS ONE 2016, 11, doi: 10.1371/journal.pone.0151561.

591

(43)

insight

into

the

genome

of

the

filter-feeder

mussel

Mytilus

Llera-Herrera, R.; Vázquez-Boucard, C.; García-Gasca, A.; Huvet, A. Co-

592

expression and regulation of ovarian vitellogenins in the Pacific oyster

593

Crassostrea gigas. Aquac. Res. 2014, 45, 448-459.

594 595

(44)

Andrew, M. N.; Dunstan, R. H.; O'Connor, W. A.; Van Zwieten, L.; Nixon, B.; MacFarlane, G. R. Effects of 4-nonylphenol and 17α-ethynylestradiol exposure

28 ACS Paragon Plus Environment

Page 29 of 34

Environmental Science & Technology

596

in the Sydney rock oyster, Saccostrea glomerata: Vitellogenin induction and

597

gonadal development. Aquat. Toxicol. 2008, 88, 39-47.

598

(45)

Andrew, M. N.; O'Connor, W. A.; Dunstan, R. H.; MacFarlane, G. R. Exposure

599

to 17α-ethynylestradiol causes dose and temporally dependent changes in

600

intersex, females and vitellogenin production in the sydney rock oyster.

601

Ecotoxicology 2010, 19, 1440-1451.

602

(46)

Michalski, A.; Damoc, E.; Lange, O.; Denisov, E.; Nolting, D.; Müller, M.;

603

Viner, R.; Schwartz, J.; Remes, P.; Belford, M.; Dunyach, J. J.; Cox, J.;

604

Horning, S.; Mann, M.; Makarov, A. Ultra high resolution linear ion trap

605

orbitrap mass spectrometer (orbitrap elite) facilitates top down LC MS/MS and

606

versatile peptide fragmentation modes. Mol. Cell. Proteomics 2012, 11, doi:

607

10.1074/mcp.O111.013698.

608

(47)

Diz, A. P.; Dudley, E.; Skibinski, D. O. F. Identification and characterization of

609

highly expressed proteins in sperm cells of the marine mussel Mytilus edulis.

610

Proteomics 2012, 12, 1949-1956.

611

(48)

Verslycke, T.; Vandenbergh, G. F.; Versonnen, B.; Arijs, K.; Janssen, C. R.

612

Induction of vitellogenesis in 17α-ethinylestradiol-exposed rainbow trout

613

(Oncorhynchus mykiss): a method comparison. Comp. Biochem. Phys. C 2002,

614

132, 483-492.

615

(49)

Morthorst, J. E.; Holbech, H.; Jeppesen, M.; Kinnberg, K. L.; Pedersen, K. L.;

616

Bjerregaard, P. Evaluation of yolk protein levels as estrogenic biomarker in

617

bivalves; Comparison of the alkali-labile phosphate method (ALP) and a

618

species-specific immunoassay (ELISA). Comp. Biochem. Phys. C 2014, 166, 88-

619

95.

29 ACS Paragon Plus Environment

Environmental Science & Technology

620

(50)

Page 30 of 34

Emmersen, B. K.; Petersen, I. M. Natural occurrence, and experimental

621

induction by estradiol-17-β, of a lipophosphoprotein (vitellogenin) in flounder

622

(Platichtys flesus, L.). Comp. Biochem. Phys. B 1976, 54, 443-446.

623

(51)

Gagnaire, B.; Gagné, F.; André, C.; Blaise, C.; Abbaci, K.; Budzinski, H.;

624

Dévier, M. H.; Garric, J. Development of biomarkers of stress related to

625

endocrine disruption in gastropods: Alkali-labile phosphates, protein-bound

626

lipids and vitellogenin-like proteins. Aquat. Toxicol. 2009, 92, 155-167.

627

(52)

2009, 37, 627-641.

628 629

(53)

Byrne, B. M.; Gruber, M.; Ab, G. The evolution of egg yolk proteins. Prog. Biophys. Mol. Bio. 1989, 53, 33-69.

630 631

Johnson, L. N. The regulation of protein phosphorylation. Biochem. Soc. T.

(54)

Tran, T. K. A.; MacFarlane, G. R.; Kong, R. Y. C.; O'Connor, W. A.; Yu, R. M.

632

K. Mechanistic insights into induction of vitellogenin gene expression by

633

estrogens in Sydney rock oysters, Saccostrea glomerata. Aquat. Toxicol. 2016,

634

174, 146-158.

635

(55)

Blom, N.; Gammeltoft, S.; Brunak, S. Sequence and structure-based prediction

636

of eukaryotic protein phosphorylation sites. Journal of Molecular Biology 1999,

637

294, 1351-1362.

638

(56)

Blom, N.; Sicheritz-Pontén, T.; Gupta, R.; Gammeltoft, S.; Brunak, S.

639

Prediction of post-translational glycosylation and phosphorylation of proteins

640

from the amino acid sequence. PROTEOMICS 2004, 4, 1633-1649.

641

(57)

IEH. Chemicals purported to be endocrine disrupters: a compilation of published

642

lists (Web Report W20). Leicester, UK, MRC Institute for Environment and

643

Health, retrieved from http://le.ac.uk/ieh.

30 ACS Paragon Plus Environment

Page 31 of 34

Environmental Science & Technology

644

(58)

Langston, W. J.; Burt, G. R.; Chesman, B. S. Feminisation of male clams

645

Scrobicularia plana from estuaries in Southwest UK and its induction by

646

endocrine-disrupting chemicals. Mar. Ecol. Prog. Ser. 2007, 333, 173-184.

647

(59)

Tankoua, O. F.; Amiard-Triquet, C.; Denis, F.; Minier, C.; Mouneyrac, C.;

648

Berthet, B. Physiological status and intersex in the endobenthic bivalve

649

Scrobicularia plana from thirteen estuaries in northwest France. Environ. Pollut.

650

2012, 167, 70-77.

651

(60)

Scott, A. P. Do mollusks use vertebrate sex steroids as reproductive hormones?

652

Part I: Critical appraisal of the evidence for the presence, biosynthesis and

653

uptake of steroids. Steroids 2012, 77, 1450-1468.

654

(61)

Scott, A. P. Do mollusks use vertebrate sex steroids as reproductive hormones?

655

II. Critical review of the evidence that steroids have biological effects. Steroids

656

2013, 78, 268-281.

657

(62)

Morishita, F.; Furukawa, Y.; Matsushima, O.; Minakata, H. Regulatory actions

658

of neuropeptides and peptide hormones on the reproduction of molluscs. Can. J.

659

Zool. 2010, 88, 825-845.

660

(63)

Tran, T. K. A.; MacFarlane, G. R.; Kong, R. Y. C.; O'Connor, W. A.; Yu, R. M.

661

K. Potential mechanisms underlying estrogen-induced expression of the

662

molluscan estrogen receptor (ER) gene. Aquat. Toxicol. 2016, 179, 82-94.

663

(64)

Ni, J.; Zeng, Z.; Ke, C. Sex steroid levels and expression patterns of estrogen

664

receptor gene in the oyster Crassostrea angulata during reproductive cycle.

665

Aquaculture 2013, 376-379, 105-116.

666

(65)

Ni, J.; Zeng, Z.; Kong, D.; Hou, L.; Huang, H.; Ke, C. Vitellogenin of Fujian

667

oyster, Crassostrea angulata: Synthesized in the ovary and controlled by

668

estradiol-17β. Gen. Comp. Endocr. 2014, 202, 35-43.

31 ACS Paragon Plus Environment

Environmental Science & Technology

669

(66)

(67)

Zhang, S.; Wang, S.; Li, H.; Li, L. Vitellogenin, a multivalent sensor and an antimicrobial effector. Int. J. Biochem. Cell B. 2011, 43, 303-305.

672 673

Liu, Y.; Beyer, A.; Aebersold, R. On the Dependency of Cellular Protein Levels on mRNA Abundance. Cell 2016, 165, 535-550.

670 671

Page 32 of 34

(68)

Jubeaux, G.; Simon, R.; Salvador, A.; Lopes, C.; Lacaze, E.; Quéau, H.;

674

Chaumot, A.; Geffard, O. Vitellogenin-like protein measurement in caged

675

Gammarus fossarum males as a biomarker of endocrine disruptor exposure:

676

Inconclusive experience. Aquat. Toxicol. 2012, 122-123, 9-18.

677

678

32 ACS Paragon Plus Environment

Page 33 of 34

Environmental Science & Technology

Figure 1. Alignment of the M. galloprovincialis vitellogenin sequence (obtained from the RNA_seq sequence Unigene21_All_6 31) with that showing the best identity (lowest E-value) in NCBI-Blast search, Pecten maximus (GenBank: CAQ06469.2). Red lines represent the peptide coverage detected by LC-MS/MS in a M. galloprovincialis mature female gonad sample. 664x1436mm (72 x 72 DPI)

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 2. Comparison of (a) indirect quantification of Vtg using ALP and (b) direct quantification using LCMS/MS in gonads of Mytilus galloprovincialis females (black bars) or males (white bars) in different stages of gametogenesis: stage II of gametogenesis (Immature) and fully mature individuals (Mature). The absent bars corresponding to males in figure (b) correspond to 0 values (undetected Vtg). 331x209mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 34 of 34