http://orcid.org/0000-0003-2262-6808
Figures
Abstract
The three-spined stickleback Gasterosteus aculeatus is an important model for studying microevolution and parallel adaptation to freshwater environments. Marine and freshwater forms differ markedly in their phenotype, especially in the number of lateral plates, which are serially repeated elements of the exoskeleton. In fishes, thyroid hormones are involved in adaptation to salinity, as well as the developmental regulation of serially repeated elements. To study how thyroid hormones influence lateral plate development, we manipulated levels of triiodothyronine and thiourea during early ontogeny in a marine and freshwater population with complete and low plate phenotypes, respectively. The development of lateral plates along the body and keel was heterochronic among experimental groups. Fish with a low dosage of exogenous triiodothyronine and those treated with thiourea exhibited retarded development of bony plates compared to both control fish and those treated with higher a triiodothyronine dosage. Several triiodothyronine-treated individuals of the marine form expressed the partial lateral plate phenotype. Some individuals with delayed development of lateral plates manifested 1–2 extra bony plates located above the main row of lateral plates.
Citation: Bolotovskiy AA, Levina MA, DeFaveri J, Merilä J, Levin BA (2018) Heterochronic development of lateral plates in the three-spined stickleback induced by thyroid hormone level alterations. PLoS ONE 13(3): e0194040. https://doi.org/10.1371/journal.pone.0194040
Editor: Paul Eckhard Witten, Universiteit Gent, BELGIUM
Received: October 17, 2017; Accepted: February 25, 2018; Published: March 9, 2018
Copyright: © 2018 Bolotovskiy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work received support from №286889 Suomen Akatemia (Academy of Finland) - JD; №218343 Suomen Akatemia (Academy of Finland) - JM; №15-04-03586 Russian Foundation for Basic Research (RFBR) - BAL; and №15-34-20416 Russian Foundation for Basic Research (RFBR) - BAL.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The three-spined stickleback Gasterosteus aculeatus L. is a circumpolar species distributed across the Northern hemisphere, and it has emerged as an important model organism for studying many contemporary problems in ecology and evolution [1–3]. The species is phenotypically very diverse, but there are two basic forms recognized as marine and freshwater ecotypes [4]. The marine form (MF) is ancestral, and its morphology is thought to be conserved since the Miocene [5, 6]. Multiple colonizations of freshwater habitats have resulted in significant phenotypic changes [7, 8]. One of the most striking morphological differences between the two main forms is in the number of bony lateral plates protecting the body. These bony plates, as well as fin rays, are part of the exoskeleton, which is largely independent from the endoskeleton in terms of function, development and evolution [9–11]. The marine form (MF) of the three-spined stickleback has numerous (30–36) lateral bony plates covering most or all of the body, while the freshwater form (FF) usually has less than 10 (typically 5–7) lateral plates along the anterior part of the body [12]. An intermediate partially plated form occurs mainly in hybrids between low and fully plated phenotypes, but also as an independently evolved form in some populations (e.g. [13, 14]).
Apart from plate numbers, there are also many other phenotypic differences between marine and freshwater forms [15–17], some of which have a genetic basis that has been identified [18–21]. The fast evolution of the MF to FF occurs via rapid assembly of freshwater genotypes from standing genetic variation in the ancestral MF, driven by strong selection acting simultaneously on multiple loci favoring life in freshwater [22–24]. However, environmental conditions have also been shown to affect the expression of traits [15, 25]. An important mediator of how genes respond to the environment is hormones [26]. In particular, thyroid hormones (THs) are known to regulate development of morphological, physiological, biochemical, and behavioral traits of many fish species, including three-spined sticklebacks [27–31]. Specifically, THs signaling pathway is involved in complex adaptations to the freshwater environment in stickleback [31]. For example, the FF exhibits a lower plasma concentration of thyroid hormone (thyroxine), as well as a lower metabolic rate, the latter of which is likely adaptive for permanent residency in small water bodies. Kitano et al. [31] suggested that divergence in hormone concentrations underlay evolutionary changes in the coordinated physiological and behavioral traits that comprise an organism's integrated phenotype.
Several recent studies have demonstrated the strong influence of THs on fish phenotype, including skeletal characters. Serially repeated elements (SREs), such as the number of lateral line scales, fin rays, and pharyngeal teeth, have been shown to change significantly in experiments with thyroid status regulation in cyprinid fishes [32–36]. In these studies, adding exogenous hormones during early development resulted in an accelerated developmental rate and fewer numbers of SREs (oligomerization), while exposure to goitrogen caused delayed development and increased SRE numbers (polymerization).
The aim of the present study was to examine the influence of thyroid hormones on the development of meristic traits, with special attention towards the lateral bony plates in the marine and freshwater forms of the three-spined stickleback. In particular, we were interested in whether the developmental and phenotypic responses of stickleback to hormonal treatments were similar to those previously observed in cyprinids.
Results
T3 concentrations in wild and experimental fish
Among the wild fish, plasma T3 concentrations differed significantly between the MF and FF groups (Mann-Whitney U-test, U = 1.000, df = 8, p = 0.039), with the MF fish having concentrations twice as high as the FF fish (Median = 0.40 ng/ml and 0.20 ng/ml, respectively; Fig 1A).
(A) Plasma T3 concentrations of wild caught three-spined sticklebacks. (B) Whole-body T3 concentrations in laboratory-raised three-spined sticklebacks (61 dpf). Middle point is median, box and whisker are quartiles and range respectively. Samples are combined and each contains between 4 and 13 individuals; the sample size for each group is indicated above the bar. Different lowercase letters above the bars indicate significant differences between forms (p < 0.05, Mann-Whitney U-test for wild fish) or between experimental groups (p < 0.05, Kruskal-Wallis test with Dunn’s post hoc test for laboratory-raised fish). Differences between groups were tested independently within each panel and each cross.
Whole body T3 concentrations differed between some of the experimental groups. The concentrations were significantly higher in TH-fish compared to the Thio-treated fish in both FF and MF crosses (Kruskal-Wallis tests, p < 0.05, Fig 1B). T3 concentration in TH-fish was slightly higher, though not significantly, than the control group in both MF and FF crosses (Fig 1B). A dose-dependent trend towards higher T3 concentration between the two TH groups in the MF cross was noted. There was no difference in T3 concentration between the Thio and control fish in the MF crosses, but the Thio-fish in the FF crosses had significantly lower T3 concentrations compared to the controls (Fig 1B).
Experimental effects on meristic characters
All non-treated fish (controls) of the MF and FF crosses demonstrated either complete or low-plated phenotypes, respectively (Fig 2A). Control fish from the MF and FF crosses differed significantly in LP (MF > FF), Va (MF > FF), Vc (FF > MF) and C (FF > MF) (Mann-Whitney, p < 0.05; Table 1). Among the MF treated groups, Thio-fish had significantly more LPs than TH-0.5 fish (Kruskal-Wallis, p < 0.05, Table 1). Moreover, these groups are significantly different, even upon exclusion of the TH-0.5 individuals with a partial lateral plate phenotype (see below).
(A) Alizarin red stained control fish of the freshwater (FF, 130 dpf age) and marine (MF, 101 dpf) forms. (B) Sites of sampling of parental fish at the White Sea basin. Maps were obtained from http://www.sasgis.org.
Four of the TH-0.5 treated individuals of the MF cross lacked several lateral plates on one (n = 3) or both sides (n = 1) of the body (Fig 3). Apart from lateral plates, there were no significant differences in mean values of all other meristic traits between treatment groups, except for a non-significant tendency for the control fish to have higher D count as compared to the Thio-fish in the FF crosses (Kruskall-Wallis, p = 0.07).
Ellipses delineate areas with no lateral plates.
Heterochronies in LP development in the marine form
Fish from the different experimental groups demonstrated heterochronic development of LP number (Fig 4A and 4B). The first anterior plates (three plates) appeared in the control, TH-0.5, and TH-1 groups almost at the same age (21–22 dpf). Thio-treated fish had no plates at 22 dpf; the first anterior plates in the Thio-fish appeared by 25 dpf (Fig 4A).
(A) Heterochronies in development (appearance) of lateral (left y-axis, four upper lines) and keel plate (right y-axis, four lower lines) in different experimental groups of the MF fish. Development of keel plates is given until forming a single continuous row with anterior plates. Plotted values are treatment-specific means. (B) Differences in degree of development of the lateral plates in MF cross at 42 dpf. Representative specimens of three sizes (largest, intermediate, and smallest) are shown from each experimental group.
The keel plates appeared first in the TH-1 fish (28 dpf), and keel plates in the control, TH-0.5 and Thio-fish appeared between 29–34 dpf (Fig 4A).
The fully plated phenotype (i.e. anterior and keel plates form a single continuous row, as seen in the typical adult marine form) was attained earliest (42 dpf) in the TH-1 individuals (Fig 4B). By 48 dpf, full plate development was complete in most of the control and TH-0.5 fish, whereas Thio-treated fish demonstrated suppression of LP development; half of the individuals took 54 dpf to express the fully plated phenotype (Fig 4A).
In contrast to lateral plates, all other meristic characters did not demonstrate obvious developmental shifts in response to either of the hormone treatments.
Induction of an additional row of bony plates by thyroid regulation
Seven of the treated individuals with a fully plated phenotype (MF) developed an additional small-sized bony plate(s) on the body, behind the head and above the main row of lateral plates (Fig 5). This extra plating was detected in four individuals from the TH-0.5 group and three individuals from the Thio group.
Additional plates are located above the main row of lateral plates, behind the head as exemplified from a TH-0.5 fish of marine form (fully plated phenotype, right side), age 101 dpf, TL = 34.7 mm (alizarin red S stained). Scale bar = 1 mm.
Discussion
When experienced early in ontogeny, experimental changes in thyroid hormone concentrations can result in developmental heterochronies of skeletal structures in fishes [27, 32, 34, 37]. In general, hormonally-induced heterochronies frequently affect the number of serially repeated skeletal elements like scales, fin rays and pharyngeal teeth [36, 38, 39]. Specifically, thyroid treatment usually induces accelerated development of calcification and squamation, which often translates into reduced numbers of skeletal elements and scales, whereas thyroid disruptors tend to delay development and can lead to increased numbers of skeletal elements [27, 32, 34, 37].
In our experiment, thyroid treatment did not have a significant effect on whole body T3 concentration. Moreover, among our lab-reared samples, T3 levels were similar between marine and freshwater offspring, both in the control and TH treatments. An earlier study of lab-reared stickleback reported contrasting results to ours–lab-reared anadromous stickleback had higher levels of T3 than resident freshwater individuals [31]. However, it is important to note the difference in age/developmental stage of the samples used in both studies. Specifically, the lab-reared stickleback in the previous study were sampled after 8–12 months of treatment and likely closer to reproductive condition, whereas our study commenced after 3 months when the fish were still juveniles. Hence, considering the correlation between thyroid status and reproductive status [40], it is possible that we did not observe differences in T3 levels between experimental fish because they were not in reproductive condition. Alternatively, the discrepancy between our results could be methodological–T3 concentration in the Kitano et al. [31] study was measured directly from plasma, whereas whole bodies were analyzed in the juveniles in our study. The same is true for the wild fish in our study, in which T3 levels were higher in the marine form, as determined from plasma samples. Hence, it is possible that T3 levels do differ between marine and freshwater lab-reared fish in our study, but these differences were diffused by other substances and lipids in the tissue that might have interfered with the assay. In fact, earlier stickleback studies measuring T3 in response to perchlorate exposure have also failed to detect differences in thyroid levels between treatment groups, in spite of significant morphological developmental differences [41, 42]. There too, T3 concentration was measured using whole-bodies, possibly suggesting that the use of whole bodies is not sensitive enough to detect changes in T3 levels in stickleback.
Regardless of the lack of detectable T3 differences between TH-treated and control lab-reared fish, thyroid treatment nevertheless had a significant effect on the development of lateral plate number. As expected, the TH-treated fish demonstrated an accelerated developmental rate of lateral plates. The final number of plates did not differ from those in the other groups, however there was a tendency for the TH-0.5-treated fish to express the partially plated phenotype (i.e. fewer plates). On the other hand, exposure to thiourea significantly retarded the development of lateral plates–the full range of LPs was registered later than in the other treatment groups. This prolonged developmental period resulted in an increased total number of lateral plates as compared with other groups (although this was only significant when compared with the TH-0.5 fish).
Interestingly, this heterochronic development only affected the marine form, most notably in the development of the keel and/or the posterior non-structural plates. Plate development is known to be largely controlled by the ectodysplasin A (Eda) gene, where different mutations lead to the expression of the different plate morphs [22]. Accordingly, most of the fully plated marine samples in our study had divergent alleles at the Eda locus compared to the low plated freshwater samples. In light of the fact that thyroid hormones can suppress or activate expression of many genes–including those involved in ossification, skeletal development and lateral line/seismosensory organs in fishes [43, 44]–it is possible that thyroid hormones also regulate the Eda gene in stickleback. Moreover, since the different forms/plate morphs responded differently to the thyroid treatments (i.e. heterochronic development only observed in the marine form), it is possible that thyroid hormones are regulated differently by the different Eda genotypes. In fact, an earlier study has demonstrated divergence in thyroid hormone physiology between marine (fully plated) and freshwater (low plated) stickleback, where both thyroid hormone levels and mRNA levels in thyroid stimulating hormone (TSHβ2) were significantly higher in marine populations [31]. In addition, only the marine form displayed changes in expression of TSHβ2 in response to photoperiod, suggesting differences in thyroid hormone signaling pathways between forms [31].
Some individuals of the fully plated phenotype (marine form) manifested additional one or two small-sized bony plates on the body, posterior to the head and above the main row of lateral plates (Fig 5). This extra row of bony plates was observed in some of the experimental fish exposed to thiourea, which suppressed/delayed development of lateral bony plates. Previous experiments investigating the effect of thyroid regulation on developmental rate in cyprinid fishes show that artificially retarded development can result in increased (or polymerized) SRE numbers, and scales in particular [32, 39, 45, 46]. Hence, the appearance of extra plates in the thiourea-treated fish are in line with earlier findings from cyprinids, particularly since both the number of scales and number of anteroposterior rows of scales were increased in earlier studies as they were in ours. However, extra plates were also observed in some of the low dosage triiodothyronine-treated fish, which was unexpected. It is noteworthy that some of these fish did express an initial delayed development until 42 dpf (Fig 4). It is possible that this was a result of nonmonotonic (non-lineare) dose-response, where the effects of low doses cannot be predicted by the effects observed at high doses [47]. We nevertheless suggest that delayed lateral plate development in three-spined stickleback is achieved through thyroid regulation, which promoted the appearance of the phenotype with additional bony scutes on the body.
Although we are not aware of any other published examples of such additional bony plates in wild or lab-reared three-spined sticklebacks, we have noted the presence of additional plates in several wild sticklebacks from the Baltic Sea. Unlike the lab-reared individuals that expressed the extra plates above the main row, these wild caught specimens expressed extra plating below the main row of plates (S1 Fig). Interestingly, lateral plate morphology has been exhaustively screened throughout the entire Baltic Sea and connecting Danish straits [48]. However, the only sampling location where stickleback with additional plates were found was along the Russian coast of the Gulf of Finland. A recent survey of dioxin concentrations in the Baltic Sea revealed that sediments from this particular region consistently had the highest levels of dioxins [49], which are known endocrine disruptors with specific effects on thyroid metabolizing enzymes [50]. For example, exposure of dioxins has been shown to reduce the concentration of thyroid hormone [51], hence it is possible that the extra plating observed in the wild stickleback from the Gulf of Finland was a result of thyroid disruption induced by the elevated dioxins.
Conclusions
Manipulation of thyroid level in developing sticklebacks resulted in heterochronies of development of lateral plates. Thiourea-treated fish had a significantly higher number of lateral plates compared to triiodothyronine-treated fish. Some sticklebacks with delayed development displayed an extra row of small-sized bony plates located above the main lateral plate row. Further studies are needed for understanding the mechanisms of thyroid effect on phenotypic features of sticklebacks, particularly those that are regulated by genes controlling ectodermal elements. In addition, thyroid hormones are a useful experimental tool for studying the mechanisms of loss and gain of serially repeated skeletal elements in fishes.
Methods
All used methods were approved by the ethics committee of Papanin Institute for Biology of Inland Waters Russian Academy of Sciences (Borok, Yaroslavl Prov., Russia), and all experiments were performed in accordance with approved guidelines and regulations.
Experimental design
Parents of the MF were caught in the Kandalaksha Bay of the White Sea near the Pertsov White Sea Biological Station of Lomonosov Moscow State University (WSBS MSU; 66°33’12"N, 33°06’17"E), Republic of Karelia, Russia. Parents of the FF were caught in an isolated freshwater lake, Mashinnoe (66°17'47"N, 33°22'03"E), Republic of Karelia, Russia (Fig 2B). Females of each form were stripped; eggs from different females within the same form (5 FF and 4 MF) were mixed in separate Petri dishes. Fish were killed with an overdose of anesthetic MS-222. Sperm was collected by macerating the testes of several males within each form (3 FF and 5 MF), and used to fertilize the mixed clutches of eggs from females of the same form. All parents were then fixed in 10% formalin and stained with Alizarin red S (Reachem, Russia); see below for details on staining protocol. Fertilized eggs were incubated in 5 L plastic containers using seawater from the White Sea or freshwater from a spring for incubation of marine and freshwater eggs, respectively. The progeny were kept at the following treatment regimes, which started just after fertilization (within 30 min) and continued until the end of experiment (101 dpf and 130 dpf for MF and FF, respectively): (i) alkaline solution of 3,5,3’-triiodothyronine T3 (Sigma Aldrich) at a concentration of 0.5 and 1 ng/ml for the MF, and 0.5 ng/ml for the FF (referred to as TH-fish); (ii) thiourea (Chimmed, Russia) solution, a goitrogen that blocks the synthesizing activity of the thyroid gland and usually leads to decreased endogenous TH-level (0.015% for both forms; referred to as Thio-fish); (iii) clean water (control). Previous studies on zebrafish D. rerio [52–54] and rainbow trout Oncorhynchus mykiss [55, 56] demonstrate that exogenous thyroid hormones are effectively absorbed by embryos and young fish from water (i.e. chorion is not a barrier for T3).
During the incubation period, 100% water changes were performed twice daily, and hormone/goitrogen were added to maintain the designated concentrations. Several days before hatching, experimental groups were transferred from the WSBS MSU to the laboratory in IBIW RAS, where seawater (18–20‰) was prepared by mixing sea salt (Tetra Marine SeaSalt, Germany) with distilled water. Water temperatures in all treatments was 10–11°C during the first days of egg incubation. After hatching, the larvae were transferred to 40 L tanks and temperature was gradually increased to 20–22°C. Tank bottoms were cleaned and 30% water changes were performed daily; designated concentrations of hormone/goitrogen were maintained through addition of each respective solution. The water in the tanks was gently aerated but not filtered. Larvae and juveniles were exposed to a 12L:12D photoperiod. Larvae were first fed with plankton food (Sera Micron, Germany) and Artemia nauplii. Frozen bloodworms and Daphnia were introduced to their diet at the juvenile stage, and provided ad libitum. The experimental and control fish were kept under identical conditions (viz. stocking density, temperature, oxygen conditions, etc.).
Genotyping of parental fish
The parents were genotyped to verify that their lateral plate phenotype were concordant with their genotypes at the Ectodisplasin A (Eda) locus, known to contribute to lateral plate phenotype [22]. DNA was extracted from ethanol-preserved pectoral fin tissue using a standard salt method [57]. An indel polymorphism in intron 1 of the Eda gene (Stn382F 5’-CCC TTA GAG AAT TTC CTA GCA G-3’, and Stn382R 5’-CTT GTC CCG GAT CAT ACG C-3’ from Colosimo et al. [22]) was amplified, generating characteristic band sizes of 150 bp for the low-plated allele, and a 218 bp fragment for the fully plated allele [22]. PCR conditions followed the PFC8 protocol from Colosimo et al. [22], except that a total of 30 cycles of PCR and an annealing temperature of 56°C were used. PCR products were visualized on 2% agarose gels. All FF parents were homozygous for the low-plated allele, and most MF parents were homozygous for the completely plated allele, except for one male that was heterozygous (S2 Fig).
Staining and morphological analyses
At each sampling point during development, three to eight individuals from each group were collected and preserved in 10% neutral-buffered formalin before further staining with alizarin red S [58]. Fish were cleared in 1% KOH, followed by gradual transfer to 100% glycerol for the final clearing and storage [59].
In addition to analysis of lateral plates development, we also studied definitive (adult) number of LP as well as other meristic traits, including number of anal (A), caudal (C), dorsal (D), and pectoral (P) fin rays; total number of lateral plates (LP); number of gill rakers on the lower and upper branches of the first branchial arch (GR); number of total (Vt), abdominal (Va), and caudal (Vc) vertebrae (Fig 6.). Number of pectoral rays, lateral plates and gill rakers were counted on both sides and given as averages. Counts of caudal vertebrae included all vertebrae having a normally developed haemal spine [25]. All counts were done by one person, at least twice, under a binocular microscope (Micromed MC-2-ZOOM) in order to minimize errors. Total length (TL) varied between 19.6–54.0 mm and 24.8–42.2 mm for FF and MF fish, respectively. Although some of the examined individuals were smaller than 28 mm in TL (the size at which three-spined stickleback attain an adult phenotype [60]), their meristic measurements did not differ from those of the other fish within the same group, and were hence included in the analyses.
A–number of anal fin rays, C–caudal fin rays, D–dorsal fin rays, P–pectoral fin rays (clipped fin shown), LP–total number of lateral plates, GR–number of gill rakers on the lower and upper branches of the first branchial arch, Vt–number of total vertebrae, Va–abdominal vertebrae, and Vc–caudal vertebrae.
Thyroid hormone assay
Thyroid hormone (T3) levels were analyzed both in wild adult fish (only females with eggs) and experimental juveniles. Blood samples were collected from wild-caught females sampled at the same locations as the parents used in the crosses described above. Blood from several females (4–13 individuals) was pooled for each T3-analysis. In total, three pooled blood samples of the FF (from a total of 19 individuals) and six pooled blood samples of the MF (from a total of 37 individuals) were analyzed. Blood was collected into sterile tubes. Plasma was centrifuged (5000 g for 30 min) and frozen at -20°C until assayed. The whole bodies of experimental fish (juveniles, 61 days post fertilization–dpf) were homogenized with preliminary freezing (-20°C). The hormone was extracted from whole-body homogenates in buffer (0.1 M phosphate buffered saline pH 7.4, 0.1% Triton X-100, 1 mM propylthiouracil–Sigma Aldrich). Samples were diluted 1:1 with buffer. Homogenates were centrifuged (10000 g for 30 min) and supernatants were analyzed for T3 concentration. Plasma T3 concentrations and whole body T3 concentrations were assessed by enzyme-linked immunosorbent assay (ELISA) using commercial kits (tT3, Monobind Inc., CA, USA) for total T3. The sensitivity of the T3 assay was 0.04 ng/ml according to the manufacturer’s instructions. Results were read on a photometer Stat Fax 303 Plus (Awareness Technology Inc., FL, USA) at 450 nm wavelength against 630 nm for the blank solution. All samples were analyzed in duplicate.
Statistical analysis
Statistical analyses of the data was performed using Statistica® software (Statsoft Inc., ver. 8.0). All data were tested for normality of variance using the Kolmogorov-Smirnov one-sample test and the Shapiro-Wilk W test. Due to small sample sizes and non-normal distributions, non-parametric Mann-Whitney U-test and Kruskal-Wallis tests were used to compare differences among wild populations and experimental groups (Kruskal-Wallis test was followed by the Dunn’s post hoc test). Results were considered significant when p < 0.05.
Supporting information
S1 Fig. Additional bony plates detected in wild-caught sticklebacks from the Russian coast of the Gulf of Finland (population code PET.
See "Sample collection" in DeFaveri & Merilä [48] for further details on collection).
https://doi.org/10.1371/journal.pone.0194040.s001
(DOCX)
S2 Fig. PCR genotyping of MF and FF parents using endogenous Eda locus polymorphism shows that MF parents are either homozygous for the 218-bp fully-plated allele (lanes 2–9) or heterozygous for the 218-bp allele and the 150-bp low-plated allele (lane 10).
In contrast, all FF parents were homozygous for the 150-bp low-plated allele (lanes 11–14, 16–19).
https://doi.org/10.1371/journal.pone.0194040.s002
(DOCX)
S1 Table. T3 concentrations in the wild caught three-spined sticklebacks and in the laboratory-raised three-spined sticklebacks.
https://doi.org/10.1371/journal.pone.0194040.s003
(XLS)
S2 Table. Values of meristic traits in different experimental groups for MF fish.
https://doi.org/10.1371/journal.pone.0194040.s004
(XLS)
S3 Table. Values of meristic traits in different experimental groups for FF fish.
https://doi.org/10.1371/journal.pone.0194040.s005
(XLS)
S4 Table. Development of lateral pates in different experimental groups for MF fish.
https://doi.org/10.1371/journal.pone.0194040.s006
(XLS)
Acknowledgments
We very acknowledge N.S. Mugue and A.E. Barmintseva for their help in the field and preliminary discussion of the obtained results. We are also very thankful to anonymous reviewers as well as to Prof. Brian Hall for their valuable critics and advices during manuscript preparation. The initial part of this study was realized at the White Sea Biological Station of Lomonosov Moscow State University (WSBS MSU).
References
- 1.
Ziuganov VV. The Family Gasterosteidae of World Fish Fauna. Fauna of USSR. Fishes 5. Leningrad: “Nauka” Publishing house; 1991.
- 2. Schluter D. Ecological causes of adaptive radiation. Am Nat. 1996; 148: 40–64.
- 3. McKinnon JS, Rundle HD. Speciation in nature: the threespine stickleback model systems. Trends Ecol Evol. 2002; 17: 480–488.
- 4.
Bell MA, Foster SA. The evolutionary biology of the threespine stickleback. Oxford: 1994.
- 5.
Bell MA. Paleobiology and evolution of threespine stickleback. In: Bell MA, Foster SA, editors. The evolutionary biology of the threespine stickleback. Oxford: 1994. pp. 438–471.
- 6. Walker JA, Bell MA. Net evolutionary trajectories of body shape evolution within a microgeographic radiation of threespine sticklebacks (Gasterosteus aculeatus). J Zool. 2000; 252: 293–302.
- 7. Leinonen T, Cano JM, Mäkinen H, Merilä J. Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks. J Evolution Biol. 2006; 19: 1803–1812.
- 8. Reid DT, Peichel CL. Perspectives on the genetic architecture of divergence in body shape in sticklebacks. Integr Comp Biol. 2010; 50: 1057–1066. pmid:21082067
- 9. Shimada A, Kawanishi T, Kaneko T, Yoshihara H, Yano T, Inohaya K, et al. Trunk exoskeleton in teleosts is mesodermal in origin. Nat Commun. 2013; 4: 1639. pmid:23535660
- 10. Hall BK. Endoskeleton/exo (dermal) skeleton—mesoderm/neural crest: two pair of problems and a shifting paradigm. J Appl Ichthyol. 2014; 30: 608–615.
- 11. Witten PE, Hall BK. Teleost skeletal plasticity: modulation, adaptation, and remodelling. Copeia. 2015; 103: 727–739.
- 12.
Bell MA. Lateral plate evolution in the threespine stickleback: getting nowhere fast. In: Hendry AP, Kinnison MT, editors. Microevolution Rate, Pattern, Process. Springer Netherlands: 2001. pp. 445–461.
- 13. Klepaker T. Lateral plate polymorphism in marine and estuarine populations of the threespine stickleback (Gasterosteus aculeatus) along the coast of Norway. Copeia. 1996; 4: 832–838.
- 14.
Mattern MY. Phylogeny, systematics and taxonomy of sticklebacks. In Östlund-Nilsson S, Mayer I, Huntingford FA, editors. Biology of the three-spined stickleback. New York: 2007. pp. 1–40.
- 15. Heuts MJ. Racial divergence in fin ray variation patterns in Gasterosteus aculeatus. J Genet. 1949; 49: 183–191. pmid:24536667
- 16. Hagen DW, Gilbertson LG. Geographic variation and environmental selection in Gasterosteus aculeatus L. in the Pacific northwest, America. Evolution. 1972; 26: 32–51. pmid:28555771
- 17. Reimchen TE, Stinson EM, Nelson JS. Multivariate differentiation of parapatric and allopatric populations of threespine stickleback in the Sangan River watershed, Queen Charlotte Islands. Can J Zool. 1985; 63: 2944–2951.
- 18. Peichel CL, Nereng KS, Ohgi KA, Cole BL, Colosimo PF, Buerkle CA, et al. The genetic architecture of divergence between threespine stickleback species. Nature. 2001; 414: 901–905. pmid:11780061
- 19. Schluter D, Clifford EA, Nemethy M, McKinnon JS. Parallel evolution and inheritance of quantitative traits. Am Nat. 2004; 163: 809–822. pmid:15266380
- 20. Berner D, Kaeuffer R, Grandchamp AC, Raeymaekers JAM, Räsänen K, Hendry AP. Quantitative genetic inheritance of morphological divergence in a lake–stream stickleback ecotype pair: implications for reproductive isolation. J Evolution Biol. 2011; 24: 1975–1983.
- 21. Leinonen T, Herczeg G, Cano JM, Merilä J. Predation-imposed selection on threespine stickleback (Gasterosteus aculeatus) morphology: a test of the refuge use hypothesis. Evolution. 2011; 65: 2916–2926. pmid:21967432
- 22. Colosimo PF, Hosemann KE, Balabhadra S, Villarreal G, Dickson M, Grimwood J, et al. Widespread parallel evolution in sticklebacks by repeated fixation of ectodysplasin alleles. Science. 2005; 307: 1928–1933. pmid:15790847
- 23. Jones FC, Grabherr MG, Chan YF, Russell P, Mauceli E, Johnson J, et al. The genomic basis of adaptive evolution in threespine sticklebacks. Nature. 2012; 484: 55–61. pmid:22481358
- 24. Terekhanova NV, Logacheva MD, Penin AA, Neretina TV, Barmintseva AE, Bazykin GA, et al. Fast evolution from precast bricks: genomics of young freshwater populations of threespine stickleback Gasterosteus aculeatus. PLoS Genet. 2014; 10.
- 25. Lindsey CC. Experimental study of meristic variation in a population of threespine sticklebacks, Gasterosteus aculeatus. Can J Zool. 1962; 40: 271–312.
- 26. Dufty AM Jr, Clobert J, Møller AP. Hormones, developmental plasticity and adaptation. Trends Ecol Evol. 2002; 17: 190–196.
- 27. Brown DD. The role of thyroid hormone in zebrafish and axolotl development. Proc Nat Acad Sci USA. 1997; 94: 13011–13016. pmid:9371791
- 28. de Jesus EGT, Toledo JD, Simpas MS. Thyroid hormones promote early metamorphosis in grouper (Epinephelus coioides) larvae. Gen Comp Endocrinol. 1998; 11: 10–16.
- 29. Hseu J-R, Chang H-F, Ting Y-Y. Inducing effects of hormones in metamorphosis-arrested grouper larvae. J Fish Soc Taiwan. 2002; 29: 79–83.
- 30. Shiao JC, Hwang PP. Thyroid hormones are necessary for the metamorphosis of tarpon Megalops cyprinoides leptocephali. J Exp Mar Biol Ecol. 2006; 331: 121–132.
- 31. Kitano J, Lema SC, Luckenbach JA, Mori S, Kawagishi Y, Kusakabe M, et al. Adaptive divergence in the thyroid hormone signaling pathway in the stickleback radiation. Curr Biol. 2010; 20: 2124–2130 pmid:21093265
- 32. Smirnov SV, Dzerzhinskii KF, Levin BA. On the relationship between scale number in the lateral line in the African barbel Barbus intermedius and the rate of ontogeny (by experimental data). J Ichthyol. 2006; 46: 129–132.
- 33. Smirnov SV, Levin BA. Reduced number of serial elements in the African barb Barbus intermedius with accelerated ontogeny: a new type of paedomorphosis. Dokl Biol Sci. 2007; 413: 140–142.
- 34. Levin BA. Drastic shift in the number of lateral line scales in the common roach Rutilus rutilus as a result of heterochronies: experimental data. J Appl Ichthyol. 2010; 26: 303–306.
- 35. Shkil FN, Borisov VB, Abdissa B, Smirnov SV. Role of thyroid hormone in the ontogeny and morphological diversification of Barbus intermedius sensu Banister, 1973 of Lake Tana in Ethiopia. Russ J Dev Biol. 2010; 41: 369–380.
- 36. Bolotovskiy AA, Levin BA. Effect of thyroid hormones on the development of pharyngeal dentition in roach Rutilus rutilus (Cyprinidae, Cypriniformes). J Ichthyol. 2011; 51: 474–483.
- 37. Shkil FN, Kapitanova DV, Borisov VB, Abdissa B, Smirnov SV. Thyroid hormone in skeletal development of cyprinids: effects and morphological consequences. J Appl Ichthyol. 2012; 28: 398–405.
- 38. Kapitanova DV, Shkil FN. Effects of thyroid hormone level alterations on the development of supraneural series in zebrafish, Danio rerio. J Appl Ichthyol. 2014; 30: 821–824.
- 39. Levin BA, Levina MA. Poly- and oligomerization of scales in the blue bream Ballerus ballerus (Cyprinidae) as a consequence of thyroid status regulation. J Appl Ichthyol. 2014; 30: 809–813.
- 40. Cyr JG, Eales JG. Interrelationships between thyroidal and reproductive endocrine systems in fish. Rev Fish Biol Fish. 1996; 6: 165–200.
- 41. Petersen AM, Dillon D, Bernhardt RR, Torunsky R, Postlethwait JH, von Hippel FA, et al. Perchlorate disrupts embryonic androgen synthesis and reproductive development in threespine stickleback without changing whole-body levels of thyroid hormone. Gen Comp Endocrinol. 2015; 210: 130–144. pmid:25448260
- 42. Gardell AM, Dillon DM, Smayda LC, von Hippel FA, Cresko WA, Postlethwait JH, et al. Perchlorate exposure does not modulate temporal variation of whole-body thyroid and androgen hormone content in threespine stickleback. Gen Comp Endocrinol. 2015; 219: 45–52. pmid:25733204
- 43. Rohner N, Bercsényi M, Orbán L, Kolanczyk ME, Linke D, Brand M, et al. Duplication of fgfr1 permits Fgf signaling to serve as a target for selection during domestication. Curr Biol. 2009; 19: 1642–1647. pmid:19733072
- 44. Hulbert AJ. Thyroid hormones and their effects: a new perspective. Biol Rev. 2000; 75: 519–631. pmid:11117200
- 45. Levin BA, Bolotovskiy AA, Levina MA. Body size determines the number of scales in cyprinid fishes as inferred from hormonal manipulation of developmental rate. J Appl Ichthyol. 2012; 28: 393–397.
- 46. Kondo S, Kuwahara Y, Kondo M, Naruse K, Mitani H, Wakamatsu Y, et al. The medaka rs-3 locus required for scale development encodes ectodysplasin-A receptor. Curr Biol. 2001; 11: 1202–1206. pmid:11516953
- 47. Vandenberg LN, Colborn T, Hayes TB, Heindel JJ, Jacobs DR Jr, Lee DH, et al. Hormones and endocrine-disrupting chemicals: low-dose effects and nonmonotonic dose responses. Endocrine reviews. 2012; 33(3): 378–455. pmid:22419778
- 48. DeFaveri J, Merilä J. Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks. J Evolution Biol. 2013; 26: 1700–1715.
- 49. Verta M, Salo S, Korhonen M, Assmuth T, Kiviranta H, Koistinen J, et al. Dioxin concentrations in sediments of the Baltic Sea–A survey of existing data. Chemosphere. 2007; 67: 1762–1775. pmid:17223167
- 50. Aoki Y. Polychlorinated biphenyls, polychloronated dibenzo-p-dioxins, and polychlorinated dibenzofurans as endocrine disrupters–what we have learned from Yusho disease. Environ Res. 2001; 86: 2–11. pmid:11386736
- 51. Sewall CH, Flagler N, Vandenheuvel JP, Clark GC, Tritscher AM, Maronpot RM, et al. Alterations in thyroid function in female Sprague-Dawley rats following chronic treatment with 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin. Toxicol Appl Pharm. 1995; 132: 237–244.
- 52. Liu YW, Lo LJ, Chan WK. Temporal expression and T3 induction of thyroid hormone receptors α1 and β1 during early embryonic and larval development in zebrafish, Danio rerio. Mol Cell Endocrinol. 2000; 159: 187–195. pmid:10687864
- 53. Liu YW, Chan WK. Thyroid hormones are important for embryonic to larval transitory phase in zebrafish. Differentiation. 2002; 70: 36–45. pmid:11963654
- 54. Walpita CN, Grommen SV, Darras VM, Van der Geyten S. The influence of stress on thyroid hormone production and peripheral deiodination in the Nile tilapia (Oreochromis niloticus). Gen Comp Endocrinol. 2007; 150: 18–25. pmid:16914151
- 55. Raine JC, Leatherland JF. Trafficking of L-triiodothyronine between ovarian fluid and oocytes of rainbow trout (Oncorhynchus mykiss). Comp Biochem Phys B. 2003; 136: 267–274.
- 56. Raine JC, Cameron C, Vijayan MM, Lamarre J, Leatherland JF. The effect of elevated oocyte triiodothyronine content on development of rainbow trout embryos and expression of mRNA encoding for thyroid hormone receptors. J Fish Biol. 2004; 65: 206–226.
- 57. Aljanabi SM, Martinez I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 1997; 25: 4692–4693. pmid:9358185
- 58. Taylor WR, van Dyke GG. Revised procedures for staining and clearing small fishes and other vertebrates for bone and cartilage study. Cybium. 1985; 9: 107–119.
- 59. Walker MB, Kimmel CB. A two-color acid-free cartilage and bone stain for zebrafish larvae. Biotech Histochem. 2007; 82: 23–28. pmid:17510811
- 60.
Lindsey CC. Factors Controlling Meristic Variation. In Hoar WS, Randall DJ, editors. Fish physiology 11. New York: 1988. pp. 197–274.