Cost of Tolerance: Physiological Consequences of Evolved

Jul 6, 2017 - School of Marine Sciences, University of Maine, Orono, Maine United States. ‡ Nicholas School of the Environment, Duke University, Dur...
1 downloads 10 Views 5MB Size
Subscriber access provided by UNIV OF NEWCASTLE

Article

Cost of tolerance: physiological consequences of evolved resistance to inhabit a polluted environment in teleost fish Fundulus heteroclitus Nishad Jayasundara, Pani W. Fernando, Joshua S. Osterberg, Kristina M Cammen, Thomas F. Schultz, and Richard T. Di Giulio Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b01913 • Publication Date (Web): 06 Jul 2017 Downloaded from http://pubs.acs.org on July 17, 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 38

Environmental Science & Technology

1

Cost of tolerance: physiological consequences of evolved resistance to inhabit a polluted

2

environment in teleost fish Fundulus heteroclitus.

3 4 5

Nishad Jayasundara1,2,*, Pani W. Fernando3, Joshua S. Osterberg4, Kristina M.

6

Cammen2,4, Thomas F. Schultz4, and Richard T. Di Giulio2

7 8

1

School of Marine Sciences, University of Maine, Orono, ME, USA.

9

2

Nicholas School of the Environment, Duke University, Durham, NC, USA.

10

3

Department of Mathematics and Information Technology, University of Leoben,

11

Leoben, Austria.

12

4

13

USA.

Duke Marine Lab, Nicholas School of the Environment, Duke University, Beaufort, NC,

14 15

*

Corresponding author : [email protected]

16 17

Key words: Adaptation, fitness costs, polycyclic aromatic hydrocarbons, pollution,

18

metabolic rate, thermal plasticity

19 20 21 22 23

ACS Paragon Plus Environment

Environmental Science & Technology

24

Abstract

25

Anthropogenic stressors, including pollutants, are key evolutionary drivers. It is

26

hypothesized that rapid evolution to anthropogenic changes may alter fundamental

27

physiological processes (e.g., energy metabolism), compromising an organism’s capacity

28

to respond to additional stressors. The Elizabeth River (ER) Superfund site represents a

29

“natural-experiment” to explore this hypothesis in several subpopulations of Atlantic

30

killifish that have evolved a gradation of resistance to a ubiquitous pollutant—polycyclic

31

aromatic hydrocarbons (PAH). We examined bioenergetic shifts and associated

32

consequences in PAH-resistant killifish by integrating genomic, physiological, and

33

modeling approaches. Population genomics data revealed that genomic regions encoding

34

bioenergetic processes are under selection in PAH-adapted fish from the most

35

contaminated ER site and ex vivo studies confirmed altered mitochondrial function in

36

these fish. Further analyses extending to differentially PAH-resistant subpopulations

37

showed organismal level bioenergetic shifts in ER fish that are associated with increased

38

cost of living, decreased performance, and altered metabolic response to temperature

39

stress—an indication of reduced thermal plasticity. A movement model predicted a

40

higher energetic cost for PAH-resistant subpopulations when seeking an optimum habitat.

41

Collectively, we demonstrate that pollution adaption and inhabiting contaminated

42

environments may result in physiological shifts leading to compromised organismal

43

capacity to respond to additional stressors.

44 45 46

ACS Paragon Plus Environment

Page 2 of 38

Page 3 of 38

Environmental Science & Technology

47

Introduction

48

Anthropogenic activities have contributed to rapid environmental change, altering the

49

trajectory and rate of biological evolution.1 Adaptive changes to abiotic stressors are

50

likely to alter physiological and biochemical processes to maintain homeostasis, but may

51

have tradeoffs. Notably, abiotic stressors often affect energy metabolism via effects on

52

reaction rates and macromolecular structures, and due to the increased energy demand to

53

mount an adequate response. Thus, stress responses over evolutionary time may pose a

54

selection pressure on organismal metabolic phenotype, potentially altering the capacity of

55

organisms to respond to additional stressors.2 Here, we examined this hypothesis in

56

subpopulations of the estuarine fish, Atlantic killifish (Fundulus heteroclitus), inhabiting

57

the Elizabeth River (ER), Virginia, USA that have evolved resistance to a highly toxic

58

class of ubiquitous pollutants, polycyclic aromatic hydrocarbons (PAHs).3

59 60

Evolutionary adaptation to anthropogenic contaminants is a key element of biological

61

pollution response. Pollution adaptation has occurred across taxa4-6 including several

62

teleost species such as Atlantic killifish that show evolved resistance to complex

63

hydrocarbon compounds and mixtures.3, 7-10 Despite a number of studies examining the

64

mechanisms of resistance in these fish, potential tradeoffs, including the capacity to

65

mount an adequate response to additional stressors, are poorly understood.7, 11 In fact,

66

consequences of contemporary evolutionary shifts in physiological processes across

67

vertebrates are largely unknown, but remain an important consideration in determining

68

“winners” and “losers” of global environmental change.12

69

ACS Paragon Plus Environment

Environmental Science & Technology

70

Killifish are remarkably eurytolerant, withstanding fluctuations in temperature, salinity,

71

and oxygen in estuaries along the North American Atlantic Coast.13 This eurytolerance to

72

physical stressors and evolved resistance to PAHs, enables the use of ER killifish as a

73

unique model to investigate tradeoffs of pollution adaptation. Several sites of the ER are

74

differentially PAH contaminated and the degree of resistance to developmental toxicity

75

of PAHs in fish inhabiting these sites directly corresponds to the level of contaminants

76

(Figure 1)3, 14, and are genetically different compared to fish from clean sites.8, 15, 16

77

Therefore, these differentially resistant ER subpopulations provide a “natural

78

experiment” to investigate consequences of evolved resistance and inhabiting a polluted

79

environment, and the impact of adaptive or ancillary shifts in physiological and

80

biochemical processes on capacity of organisms to maintain optimum performance in a

81

multi-stressor context.

82 83

To this end, a comparative physiological approach and an analysis of an existing

84

population genomic dataset was taken to examine bioenergetic consequences of PAH-

85

resistance and inhabiting a chronically polluted environment in ER fish,14 relative to a

86

clean reference subpopulation, Kings Creek (KC) (Figure 1). We focused on energy

87

metabolic processes since they may be under significant selection pressure due to (i)

88

increased ATP demand for PAH detoxification and to maintain cellular homeostasis; (ii)

89

effects of PAHs on mitochondrial integrity17 as well as cardiac development and

90

function18-20 potentially altering oxygen and nutrient circulation; and (iii) presumed role

91

of the aryl hydrocarbon receptor, a key protein involved in PAH metabolism, as a

92

mitochondrial regulator.21

ACS Paragon Plus Environment

Page 4 of 38

Page 5 of 38

Environmental Science & Technology

93 94

First, we examined a restriction site-associated DNA (RAD) genomic dataset comparing

95

fish from a highly contaminated site (Atlantic Wood (AW) Superfund site; Fig 1) with

96

Mains Creek (MC, very low PAH site in the ER) and KC to test if genes associated with

97

energy metabolism are under selection. We then examined cardiac mitochondrial

98

function utilizing an ex vivo oxygen consumption rate (OCR) assay in AW and KC fish.

99

Considering the role of cardiac bioenergetics in defining thermal limits of teleosts,22-24 we

100

evaluated cardiac mitochondrial function at 24°C and 34°C as an index for thermal

101

plasticity. Subsequently, we extended our analysis to include four ER subpopulations

102

(Figure 1) and characterized shifts in organismal aerobic metabolism at different

103

temperatures and also measured thermal tolerance. To characterize potential ecological

104

significance of altered bioenergetics, we measured swimming performance and the costs

105

of finding an optimum environment based on a novel probabilistic movement model.

106

Collectively, the current study provides a comprehensive analysis of bioenergetic

107

consequences and potential ecological significances of rapid evolution to anthropogenic

108

contaminants in inhabiting a chronically polluted environment in a vertebrate.

109 110

Materials and methods

111

Fish collection and care. Fish were collected from the reference (KC) and ER sites

112

(AW-highly contaminated; Money Point (MP)-contaminated and under remediation;

113

Jones Creek (JC)-moderately contaminated; MC-low levels of contaminants) (Figure 1,

114

SI). MP remediation started in 2009 and fish were collected from a fully remediated

115

section. Upon capture, fish were maintained in static glass tanks (30 cm × 30 cm × 75

ACS Paragon Plus Environment

Environmental Science & Technology

116

cm, at 24–25°C, 14:10 h light-dark cycle, and13–15 ppt artificial seawater (ASW; Instant

117

Ocean, Foster & Smith, Rhinelander, WI, USA) and were used within 2-3 months of

118

acclimation. Fish were fed ad libitum a mix of Cyclop-eeze (Argent Chemical

119

Laboratories, Inc., WA, USA) and Zeigler’s Adult Zebrafish Complete Diet (Zeigler

120

Bros., Inc., Gardners, PA, USA). Water changes were conducted once a week. Fish were

121

carefully chosen to represent a similar size range; weights and lengths are in Table S1.

122

Only males were chosen, but gender identification was difficult at times in small

123

juveniles. We specifically focused on wild-caught lab acclimated fish, as opposed to lab

124

reared F1 individuals, to characterize the metabolic phenotype of each subpopulation that

125

is most representative of their physiological function at each ER site. Detailed methods

126

for each experiment are available in SI.

127 128

RAD analysis. A set of genomic scaffolds that were previously characterized as

129

exhibiting genetic differentiation across multiple comparisons of fish from polluted and

130

unpolluted sites (Osterberg and colleagues, unpublished data) were queried for genes

131

related to energy metabolism and maintaining mitochondrial function and integrity based

132

on user defined categories determined by Gene Ontology information and previous

133

literature (Table S2). Data were drawn from a double digest RAD sequencing25, 26

134

analysis of 83,144 loci and 12,071 single nucleotide polymorphisms (SNPs) that

135

compared genome-wide variation among AW, KC, and MC fish (n=32 per population)

136

(SI). RAD loci with smoothed pairwise FST values27 that were significantly greater than

137

the average across the genome (P