Lifelong Exposure to Dioxin-Like PCBs Alters ... - ACS Publications

Aug 1, 2019 - offspring survival and could ultimately influence population dynamics and abundance. .... a nesting tile (102 mm long × 102 mm wide × ...
1 downloads 0 Views 4MB Size
Subscriber access provided by BUFFALO STATE

Ecotoxicology and Human Environmental Health

Lifelong exposure to dioxin-like PCBs alters paternal offspring care behavior and reduces male fish reproductive success David P. Coulter, Kara E. Huff Hartz, Maria S Sepulveda, Amy Godfrey, James E. Garvey, and Michael J. Lydy Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b03460 • Publication Date (Web): 01 Aug 2019 Downloaded from pubs.acs.org on August 2, 2019

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

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

Page 1 of 35

Environmental Science & Technology

1

Lifelong exposure to dioxin-like PCBs alters

2

paternal offspring care behavior and reduces male

3

fish reproductive success

4

David P. Coultera*, Kara E. Huff Hartza, Maria S. Sepúlvedab,c, Amy Godfreyb, James E.

5

Garveya, Michael J. Lydya

6

aCenter

7

Illinois University, Life Science II Room 251, 1125 Lincoln Drive, Carbondale, IL USA 62901

8

bDepartment

9

Lafayette, IN USA 47907

for Fisheries, Aquaculture, and Aquatic Sciences and Department of Zoology, Southern

of Forestry and Natural Resources, Purdue University, 195 Marsteller St., West

10

cLyles

11

USA 47907

12

KEYWORDS: dioxin-like PCBs; paternal offspring care; Pimephales promelas; endocrine

13

disruption; anti-androgenic; fathead minnow

School of Civil Engineering, Purdue University, 195 Marsteller St., West Lafayette, IN

14

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 35

15

ABSTRACT

16

Offspring survival, cohort performance, and ultimately population dynamics are strongly

17

influenced by maternal characteristics (e.g., fecundity), whereas paternal contribution is often

18

considered limited to genetic-driven fitness of males through sexual selection. However, male

19

contribution to reproductive success can be particularly influential in species exhibiting paternal

20

offspring care. Polychlorinated biphenyls (PCBs) are widespread, persistent contaminants that can

21

disrupt maternal reproductive processes and negatively affect offspring. In contrast, how PCBs

22

affect paternal reproductive success is largely unknown, but could ultimately affect population

23

dynamics. We examined the effects of lifelong PCB exposure on the reproductive processes of

24

male fathead minnows (Pimephales promelas), a species exhibiting sole paternal offspring care,

25

by examining endocrine-associated gene expression, testes histology, secondary sexual

26

characteristics, courtship ability, offspring care, and offspring survival. PCBs minimized male

27

secondary sexual characteristics, but did not affect gonadal endpoints or inhibit ability to court

28

females. Fathers exposed to high concentrations of dioxin-like PCBs had changes in gene

29

expression, reduced offspring care behavior, and higher embryo mortality, possibly due to fathers

30

spending less time within nests and less frequently tending to embryos. Through complex

31

interactions among gene expression, physical characteristics, and behavior, PCBs inhibit paternal

32

reproductive success and have the potential to suppress population size.

ACS Paragon Plus Environment

2

Page 3 of 35

33

Environmental Science & Technology

INTRODUCTION

34

Offspring survival and performance, and subsequent population demographics, are

35

strongly influenced by maternal characteristics and offspring investment.1,2 For example, female

36

size, age, and condition may affect the number, size, survival, and condition of offspring.1,3,4 In

37

comparison, less research has examined the role of paternal contribution to offspring and

38

subsequent population dynamics. When it is examined, male reproductive contribution is

39

oftentimes determined by the quantity and quality of sperm,5 which can be positively correlated

40

with male phenotype, including coloration6 and physical ornamentation.5 These characteristics

41

might improve offspring survival and could ultimately influence population dynamics and

42

abundance.

43

Male reproductive contributions may become particularly influential in species exhibiting

44

paternal offspring care. Male behaviors, including obtaining a nesting location7, defending

45

offspring from competitors and predators8, and tending to nests and offspring9, have been shown

46

to affect offspring fate. While a gradient of paternal care strategies exists, males have a stronger

47

effect on population dynamics in species where they care solely for offspring after birth.10

48

Parental reproductive effects are ultimately determined by a combination of parental

49

characteristics and environmental conditions.3,11 Extreme fluctuations in temperature or

50

precipitation, for example, can affect maternal nutritional and energetic provisioning from

51

somatic tissue to gonads and, therefore, influence biochemical composition and viability of

52

offspring.3,12 Natural variability in environmental conditions likely has relatively minor long-

53

term negative effects on many populations because such fluctuations oftentimes last for a

54

relatively short duration and affect only a few reproductive events (e.g., short-term drought13).

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 35

55

These effects are also mitigated by species evolving under natural variability in environmental

56

conditions and adjusting reproductive strategies accordingly.4 Anthropogenic alterations to

57

ecosystems, however, have the potential to cause long-term population effects. Chemical

58

contaminants (e.g., polychlorinated biphenyls; PCBs) are particularly problematic as they are

59

globally distributed oftentimes at high concentrations14, target numerous pathways including

60

reproductive endpoints15, bioaccumulate16, and remain in the environment for decades17, thereby

61

negatively affecting organisms throughout their lives and influencing numerous successive

62

reproductive events.

63

How PCBs affect male reproductive contribution is largely unknown. Similar to other

64

contaminants, some PCBs act as endocrine-disrupting compounds (e.g., dioxin-like or anti-

65

androgenic congeners) where they reduce testosterone levels and can produce intersex gonadal

66

tissue15,18,19. Observational studies suggest that PCBs may adversely affect aspects of paternal

67

offspring care. Field studies with birds revealed that males with increased PCB concentrations

68

spent less time defending nests containing eggs20 and were more frequently absent from nests for

69

longer periods21. Long-term effects on offspring survival have also occurred, where elevated

70

PCB concentrations in male fish spawned with uncontaminated females produced fewer

71

offspring.22 These studies support the assertion that PCB contamination reduces male

72

reproduction but provide sporadic information about mechanisms. A comprehensive

73

understanding of how lifelong PCB exposure affects the suite of male reproductive endpoints

74

would allow for more accurate predictions of how contaminants affect populations. Specifically,

75

population models are generally parameterized with data on maternal reproductive rates23, but

76

accuracy of model predictions would be improved if data on paternal and maternal reproductive

77

contributions were incorporated, particularly in species exhibiting paternal offspring care.

ACS Paragon Plus Environment

4

Page 5 of 35

78

Environmental Science & Technology

We conducted a laboratory experiment designed to identify the mechanisms by which

79

PCB contamination affects male reproductive contribution from individuals to populations. We

80

dosed male fathead minnows (Pimephales promelas), a widely distributed cavity-nesting species

81

exhibiting male secondary sexual characteristics and sole paternal offspring care, with one of two

82

PCB mixtures throughout their lives. Once mature, we bred males with uncontaminated females

83

and evaluated effects on 1) expression of genes from brain tissue involved in reproductive

84

endocrine pathways, 2) testicular cellular development, 3) secondary sexual characteristics, 4)

85

courtship and reproductive behavior, 5) offspring survival, and 6) long-term offspring

86

performance. We then used a simplified population model to illustrate the extent to which the

87

PCB effects on reproductive endpoints observed in our experiments could have on population

88

dynamics and overall abundance.

89 90

MATERIALS AND METHODS

91

Experimental design

92

The experiment consisted of three treatments: a control and two PCB treatments with

93

different congener mixtures and concentrations. Each treatment consisted of a recirculating

94

freshwater system of ten 20-L aquaria. Twenty-five, thirty-day-old fathead minnows produced

95

from a breeding colony (provided by US Environmental Protection Agency, Duluth, Minnesota,

96

USA) were placed into each aquarium and reared for six months to maturity. Throughout this

97

time, fish were dosed to treatment-specific PCB levels by feeding daily on freshwater

98

oligochaete worms, Lumbriculus variegatus. Control L. variegatus cultures were reared in

99

uncontaminated sediment, cultures for one PCB treatment (hereafter PCB-1) were reared in

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 35

100

control sediment mixed 30 d prior with Aroclor 1254, and cultures for the second PCB treatment

101

(hereafter PCB-2) were reared in sediment collected from a known, contaminated field site in a

102

local lake (see Supporting Information).24 Each day, 50 L. variegatus were removed from each

103

treatment culture and provided in equal amounts to fathead minnows (~28 mg L. variegatus·tank-

104

1·day-1).

105

similar across treatments (χ2 = 1.8, df = 2, P = 0.41). See the Supporting Information (SI) for

106

complete details on experimental design.

After growing fish in experimental aquaria to seven months old, final fish density was

107 108 109

Spawning trials At seven months old, one adult male fathead minnow from each replicate aquarium was

110

randomly selected to remain in the aquarium for spawning trials, with extra fish removed from

111

the study (see Supporting Information). Two uncontaminated females (seven months old from a

112

separately held colony) were placed into each aquarium at this time along with a nesting tile (102

113

mm long x 102 mm wide x 51 mm deep). After 24 h, breeding behavior was recorded (Go-Pro

114

Hero 3, GoPro, Inc., San Mateo, CA, USA) for 15 min, with the camera placed outside of

115

aquaria facing nest openings and camera lights covered. Recordings were used to quantify

116

breeding behavior (first five min discarded as acclimation) and included: duration each male

117

stayed inside the nest; duration a female stayed inside the nest; frequency a male contacted the

118

underside of the nest (i.e., nest maintenance); frequency a male circled inside the nest (i.e., nest

119

guarding); and, frequency a female entered the nest.25

120 121

Fish were allowed to breed for seven days. Nests were inspected daily to determine the number of viable embryos present. All adult fish were removed five days after first laying

ACS Paragon Plus Environment

6

Page 7 of 35

Environmental Science & Technology

122

embryos to prevent adults from eating hatched larvae. Larvae were fed uncontaminated, 24 h old

123

brine shrimp (Artemia) twice daily for 30 d post-hatch to assess possible paternal treatment

124

effects on offspring foraging ability and development. After 30 d post-hatch, larvae were

125

euthanized and measured for total body length and mass.

126

Male fathead minnows used in spawning trials were further examined for possible effects

127

of life-long PCB exposure. Males were visually inspected while spawning for the presence of a

128

dark colored body and dark vertical bands, which are both secondary sexual characteristics.

129

Euthanized males were photographed on their left side for geometric morphometric analysis to

130

quantify body shape. Body length and mass were measured and the brain and pituitary were

131

removed, flash-frozen in liquid nitrogen, and stored at -80°C for gene expression analyses.

132

Additional secondary sexual characteristics assessed included the fatpad (fleshy dorsal pad

133

anterior to the dorsal fin), which males use for cleaning nests, and nuptial tubercles (epidermal

134

protuberances on the face), used to stimulate females during breeding or in territoriality. Fatpads

135

were categorically scored based on the visual appearance on the body and were also assessed as

136

relative size by removing from the body and weighing to quantify fatpad index (fatpad mass /

137

body mass * 100). Nuptial tubercles were scored into one of three categories representing

138

increasing tubercle size, where higher scores reflect larger tubercle size. See the Supporting

139

Information for information on fatpad and tubercle scoring. Nuptial tubercles were also

140

quantified as total number. Testes were also weighed to determine gonadosomatic index (GSI;

141

gonad mass / body mass * 100) and preserved in 10% neutral-buffered formalin for histological

142

analyses. Fulton’s condition factor (K), an index of body condition, was calculated for offspring

143

as (M / L3)(100000), where M is total body mass (g) and L is total length (mm).

144

ACS Paragon Plus Environment

7

Environmental Science & Technology

145

Page 8 of 35

PCB quantification

146

Due to the cost of PCB analysis, three randomly selected males per treatment were

147

analyzed for PCB concentration following spawning trials. Concentration was quantified as total

148

PCBs (119 congeners), dioxin-like PCBs (12 congeners), and anti-androgenic PCBs (18

149

congeners). See the Supporting Information for complete details. Dioxin-like congeners were

150

examined because they have been found to alter maternal offspring care behavior in other

151

organisms26-28 and to be estrogenic.29 Anti-androgenic congeners were examined due to their

152

endocrine-disrupting potential in males.30,31

153 154 155

Gene expression Expression of genes involved in reproductive function and known to be influenced by

156

PCBs was quantified in males using real-time polymerase chain reaction (see Supporting

157

Information). Genes examined from pituitary include luteinizing hormone β (lhb), follicle

158

stimulating hormone β (fshb), and thyroid stimulating hormone (tshb), and genes examined from

159

brain tissue included aromatase (cyp19b) and gonadotropin releasing hormone receptor (gnrhr).

160

Although brain and pituitary were initially collected on 30 total males, gene expression was only

161

quantified on samples where personnel could ensure that the pituitary was present in the tissue

162

sample after being removed from -80°C in order to guarantee comparable data across samples

163

(17 total males).

164 165

Histological analysis of testes

ACS Paragon Plus Environment

8

Page 9 of 35

Environmental Science & Technology

166

Testes were examined because PCBs have been found to affect gonadal tissue

167

development across taxa.18,19 One testis was collected from each male and fixed in 10% buffered

168

formalin, although five testes were ultimately not retained within tissue cassettes after

169

preservation and were not analyzed. Hematoxylin and eosin (H & E) longitudinal sections were

170

prepared following standard procedures and sections examined under a light microscope (10-

171

40X) for staging and evaluation of any abnormalities.32 The thickness of germinal epithelium and

172

the percent area occupied by different cell types (primary spermatocyte, secondary spermatocyte,

173

spermatid, spermatozoa or spermatogonia) were measured from each testes on five seminiferous

174

tubules per section, with three total sections (100 μm spacing between sections) examined per

175

male.

176 177 178

Population modelling The potential for long-term, population-level consequences to emerge from results

179

observed during spawning trials were assessed using a deterministic life-stage model (stagePop

180

package in R33). Embryo, larval, juvenile, and adult life stages were used in the model and, based

181

on spawning trial results, three simulations were performed. One simulated population growth

182

under baseline conditions where the probability of embryos surviving to larvae was the mean

183

embryo survival (98%) from the control treatment, whereas the other simulations used mean

184

embryo survival from the PCB-1 (86%) and PCB-2 (97%) treatments. In order to illustrate

185

effects of reduced embryo survival on a growing population, simulations started with 100

186

individuals, had a carrying capacity of 10,000 individuals, and were simulated for 100 years.

187

ACS Paragon Plus Environment

9

Environmental Science & Technology

188

Data analysis

189

See the Supporting Information for complete details of statistical analyses.

Page 10 of 35

190 191

RESULTS

192

Treatment PCB Concentrations

193

Male fathead minnows differed in total PCB concentration (χ2 = 86, df = 2, P < 0.001),

194

where control males had lower total PCB concentration than both PCB treatments, which were

195

not statistically different from each other (Figure 1). Mean (SE) estimated oral dose of all PCB

196

congeners was 4.9·10-5 (1.7·10-6) mg PCBs·kg fish-1·day-1 for control males; 0.006 (0.0003)

197

mg·kg-1·day-1 for PCB-1 males; and 0.011 (0.0004) mg·kg-1·day-1 for PCB-2 males. Total

198

concentration of dioxin-like congeners in males also differed across all treatments (χ2 = 637, df =

199

2, P < 0.001) and was highest for the PCB-1 treatment and intermediate for the PCB-2 treatment

200

(Figure 1; SI Table S6). Total concentration of anti-androgenic congeners was lower for the

201

control treatment compared to the two PCB treatments, which did not differ (χ2 = 107, df = 2, P

202

< 0.001; Figure 1; SI Table S6). Concentrations of total PCBs and dioxin-like PCBs in L.

203

variegatus fed to male fathead minnows remained stable throughout the study (See SI Table S5).

204 205 206

Gene expression and histology

ACS Paragon Plus Environment

10

Page 11 of 35

Environmental Science & Technology

207

Expression of lhb was upregulated in males exposed to the PCB-1 treatment compared to

208

control males but was similar to males in the PCB-2 treatment (χ2 = 6.8, df = 2, P = 0.03; Figure

209

2). Expression of fshb also differed among treatments (χ2 = 7.8, df = 2, P = 0.02), where males in

210

the PCB-1 treatment had increased fshb expression compared to males in the PCB-2 treatment,

211

and expression did not differ between the control and both PCB treatments (Figure 2).

212

Expression of gnrhr (χ2 = 0.9, df = 2, P = 0.63), cyp19b (χ2 = 0.7, df = 2, P = 0.78), and tshb (χ2

213

= 2.9, df = 2, P = 0.23) did not differ among treatments (Figure 2).

214 215

Histological sections of testes revealed no difference in thickness of germinal epithelium

216

among treatments (F2,22 = 1.3, P = 0.31). Composition of sperm cell type was marginally

217

significant across all treatments (F2,22 = 2.2, P = 0.06), but corrected pairwise comparisons

218

identified no significant differences across treatments.

219 220

Physical characteristics

221

Geometric morphometric analysis identified differences in male fathead minnow

222

morphology across treatments (F2,27 = 3.2, P < 0.001), where morphology of control males was

223

different from males in both PCB treatments (all F2,27 ≥ 3.7, P ≤ 0.004), and male morphology

224

was similar between PCB treatments (F2,27 = 1.4, P = 0.17). Control males had deeper bodies

225

than males in the two PCB treatments, particularly having a larger head and anterior half of the

226

body (Figure 3A). All control males had dark colored bodies which was statistically more than

227

the 64% (PCB-1) and 55% (PCB-2) of males with dark coloration in both PCB treatments

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 35

228

(Fisher’s exact test: P = 0.002; Figure 3B). In contrast, the number of males with vertical

229

banding on their bodies was similar for all treatments (overall mean = 25% of males; Fisher’s

230

exact test: P = 0.71). Fatpad index (relative fatpad mass) was similar among treatments (F2,27 =

231

0.3, P = 0.74), but fatpad score (categorical scoring of fatpad appearance on the body) was higher

232

for control males than males in both PCB treatments (χ2 = 10.6, df = 2, P = 0.005; Figure 3C).

233

Male body length (F2,27 = 0.1, P = 0.88), body mass (F2,27 = 0.2, P = 0.78), gonadosomatic index

234

(i.e., relative gonad mass; F2,27 = 0.7, P = 0.49), total number of nuptial tubercles (used to

235

stimulate females and in territoriality; χ2 = 6.1, df = 2, P = 0.05; no differences after pair-wise

236

comparisons), and nuptial tubercle score (F2,27 = 1.2, P = 0.31) were similar across treatments

237

(Table 1).

238

ACS Paragon Plus Environment

12

Page 13 of 35

Environmental Science & Technology

239

Table 1. Mean (SE) response variables of male fathead minnow physical characteristics after six

240

month exposure to PCB treatments that were not statistically different among treatments (n = 10

241

per treatment).

242 243 244 245

Response

Control

PCB-1

PCB-2

% with vertical bands

33

18

23

Body length (mm)

58 (2)

58 (2)

59 (1)

Body mass (g)

2.5 (0.5)

2.4 (0.7)

2.5 (0.5)

Fatpad index (%)a

4.6 (0.4)

4.2 (0.9)

4.0 (0.3)

Type 1 tuberclesb

9 (1)

10 (1)

13 (2)

Type 2 tuberclesc

8 (1)

5 (2)

7 (1)

Type 3 tuberclesd

1 (0.3)

1 (0.3)

1 (0.3)

GSIe

1.4 (0.1)

1.6 (0.2)

1.6 (0.2)

aFatpad

index: relative fatpad mass. bType 1 tubercles: small nuptial tubercles present. cType 2 tubercles: nuptial tubercles with a radial base, the presence of furrows, and a jagged vertical appearance. dType 3 tubercles: pronounced tubercles with a large base and deep furrows. eGSI: gonadosomatic index.

246 247 248

Spawning behavior and success Male fathead minnows in the PCB-1 treatment spent, on average, approximately 50% less

249

time inside spawning nests (χ2 = 6.5, df = 2, P = 0.04) than males in the other treatments (Figure

250

4A), and males in the PCB-1 treatment circled 69% less often within the spawning nest, on

251

average, than control males (χ2 = 6.4, df = 2, P = 0.04; Figure 4B). Males contacted the underside

252

of the nest (i.e., nest maintenance) equally across treatments (χ2 = 2.3, df = 2, P = 0.10; Table 2),

253

and the number of times females entered nests (χ2 = 0.1, df = 2, P = 0.71) and the duration

254

females stayed inside nests (χ2 = 3.1, df = 2, P = 0.21) was similar among treatments (Table 2).

ACS Paragon Plus Environment

13

Environmental Science & Technology

255

Page 14 of 35

The number of embryos laid on nesting tiles was variable and similar across treatments

256

(χ2 = 0.1, df = 2, P = 0.95), although a higher percentage of embryos in the PCB-1 treatment

257

became nonviable before hatch compared to the control (χ2 = 6.9, df = 2, P = 0.03; Figure 4C).

258

There were no differences in body length (χ2 = 2.5, df = 2, P = 0.29), body mass (χ2 = 1.1, df = 2,

259

P = 0.57), or body condition (Fulton’s K; χ2 = 3.7, df = 2, P = 0.47; Table 2) of larvae surviving

260

from hatch to 30 d post-hatch. However, the density of 30 d post-hatch larvae in each tank was a

261

significant covariate in analyses (P < 0.001).

262 263

Table 2. Mean (SE) response variables of male fathead minnow spawning behavior, spawning

264

success, and characteristics of 30 d post-hatch (dph) offspring across PCB treatments. No response

265

variables were statistically significant across treatments (n = 10 per treatment; 30 dph larvae:

266

control n = 18; PCB-1 n = 24; PCB-2 n = 49).

267

Response

Control

PCB-1

PCB-2

P-value

Frequency male touched nest

198 (48)

77 (39)

167 (41)

0.10

Frequency female entered nest

6 (2)

6 (1)

5 (2)

0.71

Seconds female in nest

226 (70)

222 (64)

105 (49)

0.21

Eggs laid

168 (39)

178 (48)

191 (53)

0.95

Larval survival to 30 dph (%)a

14 (12)

2 (1)

3 (1)

0.38

30 dph larval body length (mm)

17 (2)

17 (1)

16 (1)

0.29

30 dph larval body mass (g)

0.0858 (0.0328)

0.0537 (0.0156)

0.0442 (0.0044)

0.57

30 dph larval condition

1.15 (0.14)

0.86 (0.03)

0.87 (0.02)

0.47

aLarval

survival to 30 dph: percentage of hatched larvae surviving to 30 dph.

ACS Paragon Plus Environment

14

Page 15 of 35

268

Environmental Science & Technology

Male whole-body total PCB concentration and concentration of anti-androgenic

269

congeners were not related to the duration males stayed within nests (Figure 5A; anti-androgenic

270

congeners: R2 < 0.001, P = 0.99), frequency males circled inside of nests (Figure 5B; anti-

271

androgenic congeners: R2 < 0.001, P = 0.99), or the proportion of clutches becoming nonviable

272

(total PCBs: R2 = 0.49, P = 0.11; anti-androgenic congeners: R2 = 0.48, P = 0.11). In contrast,

273

male whole-body concentration of dioxin-like congeners was negatively related to duration in

274

nests (Figure 5C) and number of times males circled in nests (Figure 5D), but was not related to

275

proportion of clutches becoming nonviable (R2 = 0.03, P = 0.64).

276 277 278

Population model Model simulations indicated that reduced embryo survival representing that observed in

279

the spawning trials (98% control; 86% PCB-1; 97% PCB-1) can have a disproportionate effect

280

on population size and alter population trajectories (Figure 6). The simulated PCB-1 population

281

grew more slowly than the control and PCB-2 populations throughout the time simulated, and

282

had 39% lower adult abundance than the control at the largest difference between populations.

283

After 100 years, the PCB-1 population had 10% lower adult abundance than the control

284

population.

285 286 287 288

DISCUSSION Our results demonstrated that sublethal, lifelong PCB exposure affects a suite of male reproductive endpoints, including gene expression in the brain, male secondary sexual

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 35

289

characteristics, and offspring care behavior. However, embryo survival was reduced only in

290

males exposed to high concentrations of dioxin-like PCB congeners, where male offspring care

291

behaviors were reduced. Past studies have determined that PCBs can affect individual male

292

reproductive responses (e.g., gonadal effects18) or paternal behavior21, whereas we have

293

experimentally identified mechanisms by which male reproductive success is influenced when

294

exposed to contaminants throughout life. Our work also identified that paternal exposure to PCB

295

mixtures with greater dioxin-like PCBs led to reductions in embryo survival. In a broader

296

ecological context, we also demonstrated the role that males exhibiting paternal offspring care

297

can provide to female reproductive success22 and to subsequent population dynamics.

298

The reduced male reproductive success (elevated embryo mortality) observed in one PCB

299

treatment (PCB-1) may have been due to altered male offspring care behavior, with males

300

spending less time inside nests and less frequently tending to nests. These paternal behaviors are

301

essential to prevent embryo mortality due to fungal infection by providing embryos with

302

sufficient dissolved oxygen (tail fanning) and reducing fungal spread by removing embryos that

303

develop fungus.25 While we observed effects on other male reproductive endpoints, it is unclear

304

the degree to which these contributed to reduced male reproductive success. Contaminated males

305

had diminished secondary sexual characteristics compared to control males, but this did not alter

306

male ability to court females or affect the number of eggs laid in nests. Similarly, gonadal

307

endpoints examined suggest that fertility may not have been affected since these endpoints (GSI,

308

sperm cell development, germinal epithelial thickness) were similar to control males. It is

309

possible, though, that reduced male fertility occurred in the PCB-1 treatment, since we did not

310

directly quantify sperm production or quality. Along with reduced paternal offspring care

311

behavior, this may have contributed to decreased male reproductive success. While we attempted

ACS Paragon Plus Environment

16

Page 17 of 35

Environmental Science & Technology

312

to assess effects of paternal contaminant exposure on long-term (30 d post-hatch) offspring

313

performance, these results were inconclusive. Density of 30 d post-hatch offspring was a

314

significant covariate when comparing larval final length and mass across PCB treatments,

315

indicating larvae were inadvertently food-limited during the grow-out period which likely

316

masked treatment effects, including survival to 30 d post-hatch. However, the long-term effects

317

of paternal PCB exposure on offspring survival has previously been reported22, where negative

318

relationships were found between paternal PCB concentration and abundance of offspring alive

319

at the end of their first year of life.

320

One mechanism by which PCBs affect paternal reproductive ability is via endocrine

321

disruption. Several studies have reported that PCBs can interact with nuclear receptors, including

322

the estrogen (ERα/β), androgen (AR), and thyroid (TR) receptors, leading to estrogenic, anti-

323

estrogenic, and anti-androgenic effects, as well as thyroid disruption.34,35 Male fathead minnows

324

exposed to PCB mixtures in our study were altered via morphological changes and decreased

325

paternal offspring care behavioral. We also observed an upregulation in the expression of lhb in

326

the PCB-1 treatment, and while the average fsh expression for the PCB-1 group was nearly twice

327

that of the control, differences were not significant due to the overall high amount of variability

328

in the expression of this gene. An up-regulation of lhb and fsh is suggestive of an anti-androgenic

329

compensatory response by the pituitary driven by a decrease in androgen circulating levels,

330

which although not quantified in our study, has been reported as a common response after PCB

331

exposure in several male animal models, including humans.29,36,37 No changes in the expression

332

of gnrhr, cyp19b, and tshb were detected.

333

Differences in offspring care between PCB treatments illustrate how the specific

334

chemical mixture to which males are exposed affects male reproductive behavior and success.

ACS Paragon Plus Environment

17

Environmental Science & Technology

Page 18 of 35

335

Poorer paternal offspring care in males from the PCB-1 treatment was likely due to their elevated

336

concentration of dioxin-like congeners. Laboratory studies with female rats have shown altered

337

maternal care caused by dioxin-like congeners, although effects are complex and sometimes in

338

contrast to our findings with male fish. Female rats fed a dioxin-like PCB congener spent more

339

time building28 and staying in nests26-28, more time grooming28 and nursing offspring26,28, and

340

had upregulation of oxytocin receptor (OXTR) gene expression, which is positively associated

341

with maternal offspring care.28 Data from a variety of taxa demonstrate that dioxin-like PCB

342

congeners can affect parental behavior, and future work designed to quantify dose-response

343

relationships might be able to better identify and predict these trends. However, it is unclear why

344

we did not observe diminished paternal offspring care and embryo survival in our PCB-2

345

treatment, which had males with intermediate dioxin-like PCB concentrations. Specific dioxin-

346

like congeners may play a large role in affecting male reproductive behavior, since

347

concentrations of specific dioxin-like congeners varied between our PCB treatments.

348

Additionally, a threshold concentration of one or more specific congeners may be needed to

349

elicit a behavioral response.38

350

Anti-androgenic PCB congeners potentially contributed to diminished paternal

351

reproductive success. PCB-126 is considered one of the most potent anti-androgenic PCBs;

352

however, we did not detect this congener in our fish. Instead, we detected several other known

353

anti-androgenic PCBs at concentrations > 500 ng·g-1 in both mixtures, including PCBs 101, 138,

354

and 153. A major anti-androgenic PCB (PCB-118) was detected in PCB-1 at 1,225 ng·g-1, but

355

only at a concentration of 217 ng·g-1 in PCB-2. Additional anti-androgenic PCBs (PCBs- 47,

356

170, 180) were only detected in PCB-2. Follow-up studies are required to determine the role of

357

these congeners on diminished paternal reproductive success. However, regardless of the

ACS Paragon Plus Environment

18

Page 19 of 35

Environmental Science & Technology

358

mechanism, our results demonstrate how species exhibiting paternal offspring care can respond

359

quite differently to mixtures depending on the specific chemical composition.

360

Model simulations illustrated how slight reductions in embryo survival at levels observed

361

in our experiment can affect species at the population level. These model simulations were

362

simplistic demonstrations of potential effects of reduced embryo survival and do not incorporate

363

complex interactions that would affect wild populations, although recent work using a large-

364

scale model also demonstrated strong population-level consequences due to PCB exposure in

365

whales.23 Some of the interactions not incorporated into our model could act to magnify

366

contaminant effects on population dynamics. For example, males leaving nests unattended for

367

long periods of time could not only have increased offspring mortality from lack of offspring

368

care (e.g., fungal infections), but could also have additional mortality of offspring from nest

369

predation, which alone can be a substantial source of offspring mortality.9 Also, environmental

370

stressors including altered temperature39 or additional contaminants40 could further reduce the

371

number of offspring reaching sexual maturity. In contrast, wild population abundances may be

372

minimally affected by inhibited male reproductive success due to intraspecific competition

373

among breeding males. Contaminant loads vary among individuals within a system41, so highly

374

contaminated males could be less competitive than less contaminated males for nest sites or

375

females25 and, therefore, have fewer breeding attempts. The ability of contaminated males to

376

compete for nesting sites and females is unknown, though, and needs further investigation.

377

Finally, quantifying the combined effects of PCBs on maternal and paternal reproductive success

378

will improve population models, since both sexes ultimately become dosed in the wild.

379 380

The population-level consequences of inhibited male reproductive contribution resulting from PCB exposure to populations are obviously complex. For example, PCBs can not only

ACS Paragon Plus Environment

19

Environmental Science & Technology

Page 20 of 35

381

affect reproductive endpoints, but can influence many other physiological pathways (e.g.,

382

metabolic42) that have the potential to also cause population-level effects. Nevertheless, evidence

383

from our study, similar work finding long-term paternal effects on offspring survival22, and

384

whole ecosystem studies43 demonstrate the potential for contaminant mixtures to alter population

385

dynamics. Future research should focus on bridging information gaps among individual-level

386

consequences of contaminants on reproductive success, species population dynamics, and

387

ecosystem interactions to understand the broader ecological consequences of specific

388

contaminant mixtures.

389 390

ACKNOWLEDGEMENTS

391

Amy Turlington, Federico Sinche, Andy Coursey, L. Zoe Almeida, and Tim Sesterhenn provided

392

laboratory and data analysis assistance. Funding was provided by the Center for Fisheries,

393

Aquaculture, and Aquatic Sciences at Southern Illinois University-Carbondale. The authors

394

declare no competing financial interests.

395 396

ASSOCIATED CONTENT

397

SUPPORTING INFORMATION

398

The Supporting Information is available free of charge on the ACS Publications website. Further

399

details about experimental design; materials and methods; sample processing for PCB

400

concentrations, gene expression, and histology; data and statistical analyses; results (PDF) 

ACS Paragon Plus Environment

20

Page 21 of 35

Environmental Science & Technology

401 402

REFERENCES

403

(1) Badyaev, A. V.; Hamstra, T. L.; Oh, K. P.; Seaman, D. A. A. Sex-biased maternal effects

404

reduce ectoparasite-induced mortality in a passerine bird. P. Natl. Acad Sci USA. 2006, 103,

405

14406−14411.

406

(2) Marshall, D. J.; Heppell, S. S.; Munch, S. B.; Warner, R. R. The relationship between

407

maternal phenotype and offspring quality: do older mothers really produce the best offspring?

408

Ecology 2010, 91, 2862−2873.

409 410

(3) Benton, T. G; St Clair, J. J. H.; Plaistow, S. J. Maternal effects mediated by maternal age: from life histories to population dynamics. J. Anim. Ecol. 2008, 77, 1038−1046.

411

(4) Feiner, Z. S.; Wang, H. Y.; Einhouse, D. W.; Jackson, J. R.; Rutherford, E. S.; Schelb, C.;

412

Vandergoot, C. S.; Zorn, T. G.; Höök, T. O. Thermal environment and maternal effects shape egg

413

size in a freshwater fish. Ecosphere 2016, 7, e01304.

414 415

(5) Kortet, R.; Vainikka, A.; Rantala, M. J.; Taskinen, J. Sperm quality, secondary sexual characters and parasitism in roach (Rutilus rutilus L.). Biol. J. Linn. Soc. 2004, 81, 111−117.

416

(6) Sheldon, B. C.; Merilö, J.; Qvarnström, A.; Gustafsson, L.; Ellegren, H. Paternal genetic

417

contribution to offspring condition predicted by size of male secondary sexual character. P. R. Soc.

418

B 1997, 264, 297−302.

419 420

(7) Townsend, D. S.; Stewart, M. M.; Pough, F. H. Male parental care and its adaptive significance in a neotropical frog. Anim. Behav. 1984, 32, 421−431.

ACS Paragon Plus Environment

21

Environmental Science & Technology

421 422 423 424 425 426 427 428

Page 22 of 35

(8) Hakkarainen, H.; Korpimäki, E. Nest defence of Tengmalm's owls reflects offspring survival prospects under fluctuating food conditions. Anim. Behav. 1994, 48, 843−849. (9) Bain, M. B.; Helfrich, L. A. Role of male parental care in survival of larval bluegills. T. Am. Fish. Soc. 1983, 112, 47−52. (10) Møller, A. P.; Thornhill, R. Male parental care, differential parental investment by females and sexual selection. Anim. Behav. 1998, 55, 1507−1515. (11) Bestion, E.; Teyssier, A.; Aubret, F.; Clobert, J.; Cote, J. Maternal exposure to predator scents: offspring phenotypic adjustment and dispersal. P. R. Soc. B 2014, 281, 20140701.

429

(12) Feiner, Z. S.; Coulter, D. P.; Guffey, S. C.; Höök, T. O. Does overwinter temperature affect

430

maternal body composition and egg traits in yellow perch Perca flavescens? J. Fish Biol. 2016 88,

431

1524−1543.

432 433

(13) Matthews, W. J.; Marsh-Matthews, E. Effects of drought on fish across axes of space, time and ecological complexity. Freshwater Biol. 2003, 48, 1232−1253.

434

(14) Meijer, S. N.; Ockenden, W. A.; Sweetman, A.; Breivik, K.; Grimalt, J. O.; Jones, K. C.

435

Global distribution and budget of PCBs and HCB in background surface soils: implications for

436

sources and environmental processes. Environ. Sci. Technol. 2003, 37, 667−672.

437 438

(15) Oberdörster, E.; Cheek, A. O. Gender benders at the beach: endocrine disruption in marine and estuarine organisms. Environ. Toxicol. Chem. 2001, 20, 23−36.

439

(16) Sobek, A.; McLachlan, M. S.; Borgå, K.; Asplund, L.; Lundstedt-Enkel, K.; Polder, A.;

440

Gustafsson, Ö. A comparison of PCB bioaccumulation factors between an arctic and a temperate

441

marine food web. Sci. Total Environ. 2010, 408, 2753−2760.

ACS Paragon Plus Environment

22

Page 23 of 35

442 443

Environmental Science & Technology

(17) Sinkkonen, S.; Paasivirta, J. Degradation half-life times of PCDDs, PCDFs and PCBs for environmental fate modeling. Chemosphere 2000, 40, 943−949.

444

(18) Ahmad, S. U.; Tariq, S.; Jalali, S.; Ahmad, M. M. Environmental pollutant Aroclor 1242

445

(PCB) disrupts reproduction in adult male rhesus monkeys (Macaca mulatta). Environ. Res. 2003,

446

93, 272−278.

447

(19) Jofré, M. B.; Karasov, W. H. Effect of mono-ortho and di-ortho substituted polychlorinated

448

biphenyl (PCB) congeners on leopard frog survival and sexual development. Chemosphere 2008,

449

70, 1609−1619.

450 451

(20) Fyfe, R. W.; Risebrough, R. W.; Walker, W. Pollutant effects on the reproduction of the prairie falcons and merlins of the Canadian prairies. Can. Field Nat. 1976, 90, 346−355.

452

(21) Bustnes, J. O.; Bakken, V.; Erikstad, K. E.; Mehlum, F.; Skaare, J. U. Patterns of incubation

453

and nest-site attentiveness in relation to organochlorine (PCB) contamination in glaucous gulls. J.

454

Appl. Ecol. 2001, 38, 791−801.

455 456

(22) Hillis, J. J.; Garvey, J. E.; Lydy, M. J. Contaminants reduce male contribution to reproduction at the population scale. Ecosphere 2015, 6, 1−12.

457

(23) Desforges, J.; Hall, A.; McConnell, B.; Rosing-Asvid, A.; Barber, J. L.; Brownlow, A.; De

458

Guise, S.; Eulaers, I.; Jepson, P. D.; Letcher, R. J.; Levin, M.; Ross, P. S.; Samarra, F.; Víkingson,

459

G.; Sonne, C.; Dietz, R. Predicting global killer whale population collapse from PCB pollution.

460

Science 2018, 361, 1373−1376.

ACS Paragon Plus Environment

23

Environmental Science & Technology

Page 24 of 35

461

(24) Archer, M. C.; Harwood, A. D.; Nutile, S. A.; Huff Hartz, K. E.; Mills, M. A.; Garvey, J.

462

E.; Lydy, M. J. The value of using multiple metrics to evaluate PCB exposure. Arch. Environ. Con.

463

Tox. 2018, 74, 361−371.

464

(25) Salierno, J. D.; Kane, A. S. 17α-Ethinylestradiol alters reproductive behaviors, circulating

465

hormones, and sexual morphology in male fathead minnows (Pimephales promelas). Environ.

466

Toxicol. Chem. 2009, 28, 953−961.

467 468 469 470

(26) Cummings, J. A.; Nunez, A. A.; Clemens, L. G. A cross-fostering analysis of the effects of PCB 77 on the maternal behavior of rats. Physiol. Behav. 2005, 85, 83−91. (27) Simmons, S. L.; Cummings, J. A.; Clemens, L. G.; Nunez, A. A. Exposure to PCB 77 affects the maternal behavior of rats. Physiol. Behav. 2005, 84, 81−86.

471

(28) Dover, E. N.; Mankin, D. E.; Cromwell, H. C.; Phuntumart, V.; Meserve, L. A.

472

Polychlorinated biphenyl exposure alters oxytocin receptor gene expression and maternal behavior

473

in rat model. Endocrine Disruptors 2015, 3, e979681.

474

(29) Plíšková, M.; Vondráček, J.; Canton, R. F.; Nera, J.; Kočan, A.; Petrík, J.; Trnovec, T.;

475

Sanderson, T.; van den Berg, M.; Machala, M. Impact of polychlorinated biphenyls contamination

476

on estrogenic activity in human male serum. Environ. Health Persp. 2005, 113, 1277−1284.

477

(30) Goncharov, A.; Rej, R.; Negoita, S.; Schymura, M.; Santiago-Rivera, A.; Morse, G.;

478

Carpenter, D. O. Lower serum testosterone associated with elevated polychlorinated biphenyl

479

concentrations in Native American men. Environ. Health Persp. 2009, 117, 1454−1460.

ACS Paragon Plus Environment

24

Page 25 of 35

Environmental Science & Technology

480

(31) Takeuchi, S.; Anezaki, K.; Kojima, H. Effects of unintentional PCBs in pigments and

481

chemical products on transcriptional activity via aryl hydrocarbon and nuclear hormone receptors.

482

Environ. Pollut. 2017, 227, 306−313.

483 484 485 486

(32) Johnson, R.; Wolff, J.; Braunbeck, T. OECD guidance document for the diagnosis of endocrine-related histopathology of fish gonads. Histopathology 2009, 1−96. (33) Kettle, H.; Nutter, D. stagePop: modelling stage‐structured populations in r. Methods in Ecol. Evol. 2015, 6, 1484−1490.

487

(34) Bonefeld-Jørgensen, E. C.; Andersen, H. R.; Rasmussen, T. H.; Vinggaard, A. M. Effect

488

of highly bioaccumulated polychlorinated biphenyl congeners on estrogen and androgen receptor

489

activity. Toxicology 2001, 158, 141−153.

490

(35) Zhang, Q.; Lu, M.; Wang, C.; Du, J.; Zhou, P.; Zhao, M. Characterization of estrogen

491

receptor α activities in polychlorinated biphenyls by in vitro dual-luciferase reporter gene assay.

492

Environ. Pollut. 2014, 189, 169−175.

493

(36) Murugesan, P.; Muthusamy, T.; Balasubramanian, K.; Arunakaran, J. Studies on the

494

protective role of vitamin C and E against polychlorinated biphenyl (Aroclor 1254)—induced

495

oxidative damage in Leydig cells. Free Radical Res. 2005, 39, 1259−1272.

496

(37) Selvakumar, K.; Banu, L. S.; Krishnamoorthy, G.; Venkataraman, P.; Elumalai, P.;

497

Arunakaran, J. Differential expression of androgen and estrogen receptors in PCB (Aroclor 1254)-

498

exposed rat ventral prostate: impact of alpha-tocopherol. Exp. Toxicol. Pathol. 2011, 63, 105−112.

499

(38) Kuzyk, Z. A.; Hodson, P. V.; Solomon, S. M.; Reimer, K. J. Biological responses to PCB

500

exposure in shorthorn sculpin from Saglek Bay, Labrador. Sci. Total Environ. 2005, 351, 285−300.

ACS Paragon Plus Environment

25

Environmental Science & Technology

Page 26 of 35

501

(39) Coulter, D. P.; Höök, T. O.; Mahapatra, C. T.; Guffey, S. C.; Sepúlveda, M. S. Fluctuating

502

water temperatures affect development, physiological responses and cause sex reversal in fathead

503

minnows. Environ. Sci. Technol. 2015, 49, 1921−1928.

504

(40) Boone, M. D.; Semlitsch, R. D.; Little, E. E.; Doyle, M. C. Multiple stressors in amphibian

505

communities: effects of chemical contamination, bullfrogs, and fish. Ecol. Appl. 2007, 17,

506

291−301.

507

(41) Rasmussen, J. B.; Rowan, D. J.; Lean, D. R. S.; Carey, J. H. Food chain structure in Ontario

508

lakes determines PCB levels in lake trout (Salvelinus namaycush) and other pelagic fish. Can. J.

509

Fish. Aquat. Sci. 1990, 47, 2030−2038.

510

(42) Schmidt, K.; Steinberg, C. E.; Pflugmacher, S.; Staaks, G. B. Xenobiotic substances such

511

as PCB mixtures (Aroclor 1254) and TBT can influence swimming behavior and

512

biotransformation activity (GST) of carp (Cyprinus carpio). Environ. Toxicol. 2004, 19, 460−470.

513

(43) Kidd, K. A.; Blanchfield, P. J.; Mills, K. H.; Palace, V. P.; Evans, R. E.; Lazorchak, J. M.;

514

Flick, R. W. Collapse of a fish population after exposure to a synthetic estrogen. P. Natl. Acad.

515

USA 2007, 104, 8897−8901.

ACS Paragon Plus Environment

26

Page 27 of 35

Environmental Science & Technology

516

Figure captions

517

Figure 1. Mean (SE) whole-body concentration of all PCB congeners (total), anti-androgenic

518

congeners, and dioxin-like congeners from male fathead minnows after six-month exposure to

519

PCBs (n = 3 males per treatment). Different letters denote differences across treatments within

520

congener categories.

521

Figure 2. Mean (SE) gene expression from brain (gnrhr, cyp19b) and pituitary (lhb, fshb, tshb)

522

samples of male fathead minnows exposed to PCB treatments for six months during development

523

(control n = 5; PCB-1 n=6; PCB-2 n = 6). Different letters indicate differences in expression across

524

treatments for each gene. Reference gene was ribosomal protein L8 (rpl8) and expression was

525

normalized to the control.

526

Figure 3. Differences in secondary sexual characteristics of male fathead minnows exposed to one

527

of three PCB treatments for six months during development. (A) Warped outline drawing (scale

528

factor: 3) illustrating mean morphology of males from the control and males from both PCB

529

treatments. (B) Percent of mature males from each treatment having dark body coloration and (C)

530

mean (SE) fatpad score. Different letters indicate statistical differences among treatments (n = 10

531

per treatment).

532

Figure 4. Differences in spawning behavior and success in male fathead minnows exposed to one

533

of three PCB treatments for six months during development. Mean (SE) (A) duration males

534

remained in the spawning nest, (B) frequency males circled inside of nests (indicator of nest

535

guarding), and (C) percent clutch mortality. Different letters indicate differences across treatments

536

(n = 10 per treatment).

ACS Paragon Plus Environment

27

Environmental Science & Technology

Page 28 of 35

537

Figure 5. Relationships between male fathead minnow whole-body PCB concentration and

538

spawning behavior (points represent individual males; n = 3 per treatment). Total PCB

539

concentration was not related to (A) duration males spent inside spawning nests or (B) frequency

540

males circled inside nests (nest guarding). However, concentration of dioxin-like congeners was

541

negatively related to (C) time in nests and (D) nest guarding. Several points represent multiple

542

individuals that had similar PCB concentrations or spawning behavioral responses. PCB

543

concentrations below reporting limit were set to 0 ng·g-1. Similar analyses with anti-androgenic

544

PCBs were not significant.

545

Figure 6. Modeled fathead minnow population growth under baseline embryo survival rates

546

(Control: mean 98% embryo survival), similar embryo survival rates from the PCB-2 treatment

547

(mean 97% survival), and reduced embryo survival rates from the PCB-1 treatment (mean 86%

548

survival). The arrow indicates the percent difference in adult abundance at the point of maximum

549

difference in population trajectories.

550

ACS Paragon Plus Environment

28

Page 29 of 35

Environmental Science & Technology

551

AUTHOR INFORMATION

552

Corresponding Author

553

*E-mail: [email protected]. Phone: (618) 453-2608

554 555 556

TOC ART

557

ACS Paragon Plus Environment

29

Environmental Science & Technology Page 30 of 35

ACS Paragon Plus Environment

Page 31 of 35

Environmental Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology Page 32 of 35

ACS Paragon Plus Environment

Page 33 of Environmental 35 Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology

ACS Paragon Plus Environment

Page 34 of 35

Page 35 ofEnvironmental 35 Science & Technology

ACS Paragon Plus Environment