Genetic and morphometric divergence in threespine stickleback in the Chignik catchment, Alaska

Intergenerational effects of ocean temperature variation: Early life benefits are short-lived in threespine stickleback

Current climate change models predict an increase in temperature variability and extreme events such as heatwaves, and organisms need to cope with consequent changes to environmental variation. Non-genetic inheritance mechanisms can enable parental generations to prime their offspring's abilities to acclimate to environmental change-but they may also be deleterious. When parents are exposed to predictable environments, intergenerational plasticity can lead to better offspring trait performance in matching environments. Alternatively, parents exposed to variable or unpredictable environments may use plastic bet-hedging strategies to adjust the phenotypic variance among offspring. Here, we used a model species, the threespine stickleback (Gasterosteus aculeatus), to test whether putatively adaptive intergenerational effects can occur in response to shifts in environmental variation as well as to shifts in environmental mean, and whether parents employ plastic bet-hedging strategies in response to increasing environmental variation. We used a full-factorial, split-clutch experiment with parents and offspring exposed to three temperature regimes: constant, natural variation, and increased variation. We show that within-generation exposure to increased temperature variation reduces growth of offspring, but having parents that were exposed to natural temperature variation during gametogenesis may offset some early-life negative growth effects. However, these mitigating intergenerational effects do not appear to persist later in life. We found no indication that stickleback mothers plastically altered offspring phenotypic variance (egg size or clutch size) in response to temperature variation. However, lower inter-individual variance of juvenile fish morphology in offspring of increased variation parents may imply the presence of conservative bet-hedging strategies in natural populations. Overall, in our experiment, parental exposure to temperature variation had limited effects on offspring fitness-related traits. Natural levels of environmental variation promoted a potentially adaptive intergenerational response in early life development, but under more challenging conditions associated with increased environmental variation, the effect was lost.

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.

Salinity-induced transcriptome profiles in marine and freshwater threespine stickleback after an abrupt 6-hour exposure

Saltwater and freshwater environments have opposing physiological challenges, yet, there are fish species that are able to enter both habitats during short time spans, and as individuals they must therefore adjust quickly to osmoregulatory contrasts. In this study, we conducted an experiment to test for plastic responses to abrupt salinity changes in two populations of threespine stickleback, Gasterosteus aculeatus, representing two ecotypes (freshwater and ancestral saltwater). We exposed both ecotypes to abrupt native (control treatment) and non-native salinities (0‰ and 30‰) and sampled gill tissue for transcriptomic analyses after 6 h of exposure. To investigate genomic responses to salinity, we analyzed four different comparisons; one for each ecotype (in their control and exposure salinity; (1) and (2), one between ecotypes in their control salinity (3), and the fourth comparison included all transcripts identified in (3) that did not show any expressional changes within ecotype in either the control or the exposed salinity (4)). Abrupt salinity transfer affected the expression of 10 and 1530 transcripts for the saltwater and freshwater ecotype, respectively, and 1314 were differentially expressed between the controls, including 502 that were not affected by salinity within ecotype (fixed expression). In total, these results indicate that factors other than genomic expressional plasticity are important for osmoregulation in stickleback, due to the need for opposite physiological pathways to survive the abrupt change in salinity.

Parallel and non-parallel divergence within polymorphic populations of brook stickleback, Culaea inconstans (Actinopterygii: Gasterosteidae)

Studying parallel evolution allows us to draw conclusions about the repeatability of adaptive evolution. Whereas populations likely experience similar selective pressures in similar environments, it is not clear if this will always result in parallel divergence of ecologically relevant traits. Our study investigates the extent of parallelism associated with the evolution of pelvic spine reduction in brook stickleback populations. We find that populations with parallel divergence in pelvic spine morphology do not exhibit parallel divergence in head and body morphology but do exhibit parallel divergence in diet. In addition, we compare these patterns associated with pelvic reduction in brook stickleback to well-studied patterns of divergence between spined and unspined threespine stickleback. Whereas spine reduction is associated with littoral habitats and a benthic diet in threespine stickleback, spine reduction in brook stickleback is associated with a planktonic diet. Hence, we find that pelvic spine divergence is associated with largely non-parallel ecological consequences across species.

Allometric relationships in morphological traits associated with foraging, swimming ability, and predator defense reveal adaptations toward brackish and freshwater environments in the threespine stickleback

Freshwater colonization by threespine stickleback has led to divergence in morphology between ancestral marine and derived freshwater populations, making them ideal for studying natural selection on phenotypes. In an open brackish-freshwater system, we previously discovered two genetically distinct stickleback populations that also differ in geometric shape: one mainly found in the brackish water lagoon and one throughout the freshwater system. As shape and size are not perfectly correlated, the aim of this study was to identify the morphological trait(s) that separated the populations in geometric shape. We measured 23 phenotypes likely to be important for foraging, swimming capacity, and defense against predation. The lateral plate morphs in freshwater displayed few significant changes in trait sizes, but the low plated expressed feeding traits more associated with benthic habitats. When comparing the completely plated genetically assigned populations, the freshwater, the hybrids, the migrants and the lagoon fish, many of the linear traits had different slopes and intercepts in trait-size regressions, precluding our ability to directly compare all traits simultaneously, which most likely results from low variation in body length for the lagoon and migrant population. We found the lagoon stickleback population to be more specialized toward the littoral zone, displaying benthic traits such as large, deep bodies with smaller eyes compared to the freshwater completely plated morph. Further, the lagoon and migrant fish had an overall higher body coverage of lateral plates compared to freshwater fish, and the dorsal and pelvic spines were longer. Evolutionary constraints due to allometric scaling relationships could explain the observed, overall restricted, differences in morphology between the sticklebacks in this study, as most traits have diversified in common allometric trajectories. The observed differences in foraging and antipredation traits between the fish with a lagoon and freshwater genetic signature are likely a result of genetic or plastic adaptations toward brackish and freshwater environments.

Tempo and mode of allopatric divergence in the weakly electric fish Sternopygus dariensis in the Isthmus of Panama

Spatial isolation is one of the main drivers of allopatric speciation, but the extent to which spatially-segregated populations accumulate genetic differences relevant to speciation is not always clear. We used data from ultraconserved elements (UCEs) and whole mitochondrial genomes (i.e., mitogenomes) to explore genetic variation among allopatric populations of the weakly electric fish Sternopygus dariensis across the Isthmus of Panama. We found strong genetic divergence between eastern and western populations of S. dariensis. Over 77% of the UCE loci examined were differentially fixed between populations, and these loci appear to be distributed across the species' genome. Population divergence occurred within the last 1.1 million years, perhaps due to global glaciation oscillations during the Pleistocene. Our results are consistent with a pattern of genetic differentiation under strict geographic isolation, and suggest the presence of incipient allopatric species within S. dariensis. Genetic divergence in S. dariensis likely occurred in situ, long after the closure of the Isthmus of Panama. Our study highlights the contribution of spatial isolation and vicariance to promoting rapid diversification in Neotropical freshwater fishes. The study of spatially-segregated populations within the Isthmus of Panama could reveal how genetic differences accumulate as allopatric speciation proceeds.

Study of morphological variation of northern Neotropical Ariidae reveals conservatism despite macrohabitat transitions

BACKGROUND

Morphological convergence triggered by trophic adaptations is a common pattern in adaptive radiations. The study of shape variation in an evolutionary context is usually restricted to well-studied fish models. We take advantage of the recently revised systematics of New World Ariidae and investigate skull shape evolution in six genera of northern Neotropical Ariidae. They constitute a lineage that diversified in the marine habitat but repeatedly adapted to freshwater habitats. 3D geometric morphometrics was applied for the first time in catfish skulls and phylogenetically informed statistical analyses were performed to test for the impact of habitat on skull diversification after habitat transition in this lineage.

RESULTS

We found that skull shape is conserved throughout phylogeny. A morphospace analysis revealed that freshwater and marine species occupy extreme ends of the first principal component axis and that they exhibit similar Procrustes variances. Yet freshwater species occupy the smallest shape space compared to marine and brackish species (based on partial disparity), and marine and freshwater species have the largest Procrustes distance to each other. We observed a single case of shape convergence as derived from 'C-metrics', which cannot be explained by the occupation of the same habitat.

CONCLUSIONS

Although Ariidae occupy such a broad spectrum of different habitats from sea to freshwater, the morphospace analysis and analyses of shape and co-variation with habitat in a phylogenetic context shows that conservatism dominates skull shape evolution among ariid genera.

Phenotypic plasticity in gene expression and physiological response in red drum Sciaenops ocellatus exposed to a long-term freshwater environment

Red drum (Sciaenops ocellatus) is a euryhaline fish commonly found in the Gulf of Mexico and along the Atlantic coast of North America. Because of high commercial demand and its euryhaline characteristics, aquaculture of this species has diversified from marine to low-salinity aquaculture systems. In recent years, interest in the feasibility of producing red drum in inland freshwater systems has grown and this prompted us to investigate its osmoregulatory capacity after rearing for 8 months in a freshwater aquaculture system. We compared the activities of several genes and enzymes involved in the osmoregulatory process in freshwater-acclimatized (FW) and seawater (SW) red drum. The gene expression profiles were variable: the expression of genes encoding Na⁺/K⁺-ATPase (NKA) and the cystic fibrosis transmembrane regulator (CFTR) was slightly higher in SW than FW fish, while phosphoenolpyruvate carboxykinase (PEPCK) and the glucocorticoid receptor messenger RNA (mRNA) levels were higher in FW red drum. The total plasma K concentration was 60.3% lower, and gill NKA activity was 63.5% lower in FW than in SW fish. PEPCK activity was twofold higher in FW than in SW red drum. Similarly, liver glycogen was 60% higher in FW fish. In summary, both gene expression and the enzyme activity data support the phenotypic plasticity of red drum and suggest that the limited capacity for ion homeostasis observed, in particular the low plasma K concentration, was due to the composition of freshwater and does not necessarily reflect a physiological inability to osmoregulate.

Low genetic and phenotypic divergence in a contact zone between freshwater and marine sticklebacks: gene flow constrains adaptation

BACKGROUND

Distinct hybrid zones and phenotypic and genomic divergence is often observed between marine and freshwater threespine sticklebacks (Gasterosteus aculeatus). Nevertheless, cases also exist where marine-freshwater divergence is diffuse despite seemingly similar environmental settings. In order to assess what characterizes these highly different outcomes, we focused on the latter kind of system in the Odder River, Denmark. Here, a previous study based on RAD (Restriction site Associated DNA) sequencing found non-significant genome-wide differentiation between marine and freshwater sticklebacks. In the present study, we analyzed samples on a finer geographical scale. We assessed if the system should be regarded as panmictic, or if fine-scale genetic structure and local selection was present but dominated by strong migration. We also asked if specific population components, that is the two sexes and different lateral plate morphs, contributed disproportionally more to dispersal.

RESULTS

We assessed variation at 96 SNPs and the Eda gene that affects lateral plate number, conducted molecular sex identification, and analyzed morphological traits. Genetic differentiation estimated by F_ST was non-significant throughout the system. Nevertheless, spatial autocorrelation analysis suggested fine scale genetic structure with a genetic patch size of 770 m. There was no evidence for sex-biased dispersal, but full-plated individuals showed higher dispersal than low- and partial-plated individuals. The system was dominated by full-plated morphs characteristic of marine sticklebacks, but in the upstream part of the river body shape and frequency of low-plated morphs changed in the direction expected for freshwater sticklebacks. Five markers including Eda were under possible diversifying selection. However, only subtle clinal patterns were observed for traits and markers.

CONCLUSIONS

We suggest that gene flow from marine sticklebacks overwhelms adaptation to freshwater conditions, but the short genetic patch size means that the effect of gene flow on the most upstream region must be indirect and occurs over generations. The occurrence of both weak unimodal and strong bimodal hybrid zones within the same species is striking. We suggest environmental and demographic factors that could determine these outcomes, but also highlight the possibility that long-term population history and the presence or absence of genomic incompatibilities could be a contributing factor.

Small changes in gene expression of targeted osmoregulatory genes when exposing marine and freshwater threespine stickleback (Gasterosteus aculeatus) to abrupt salinity transfers

Salinity is one of the key factors that affects metabolism, survival and distribution of fish species, as all fish osmoregulate and euryhaline fish maintain osmotic differences between their extracellular fluid and either freshwater or seawater. The threespine stickleback (Gasterosteus aculeatus) is a euryhaline species with populations in both marine and freshwater environments, where the physiological and genomic basis for salinity tolerance adaptation is not fully understood. Therefore, our main objective in this study was to investigate gene expression of three targeted osmoregulatory genes (Na⁺/K⁺-ATPase (ATPA13), cystic fibrosis transmembrane regulator (CFTR) and a voltage gated potassium channel gene (KCNH4) and one stress related heat shock protein gene (HSP70)) in gill tissue from marine and freshwater populations when exposed to non-native salinity for periods ranging from five minutes to three weeks. Overall, the targeted genes showed highly plastic expression profiles, in addition the expression of ATP1A3 was slightly higher in saltwater adapted fish and KCNH4 and HSP70 had slightly higher expression in freshwater. As no pronounced changes were observed in the expression profiles of the targeted genes, this indicates that the osmoregulatory apparatuses of both the marine and landlocked freshwater stickleback population have not been environmentally canalized, but are able to respond plastically to abrupt salinity challenges.

Is isolation by adaptation driving genetic divergence among proximate Dolly Varden char populations?

Numerous studies of population genetics in salmonids and other anadromous fishes have revealed that population structure is generally organized into geographic hierarchies (isolation by distance), but significant structure can exist in proximate populations due to varying selective pressures (isolation by adaptation). In Chignik Lakes, Alaska, anadromous Dolly Varden char (Salvelinus malma) spawn in nearly all accessible streams throughout the watershed, including those draining directly to an estuary, Chignik Lagoon, into larger rivers, and into lakes. Collections of Dolly Varden fry from 13 streams throughout the system revealed low levels of population structure among streams emptying into freshwater. However, much stronger genetic differentiation was detected between streams emptying into freshwater and streams flowing directly into estuarine environments. This fine-scale reproductive isolation without any physical barriers to migration is likely driven by differences in selection pressures across freshwater and estuarine environments. Estuary tributaries had fewer larger, older juveniles, suggesting an alternative life history of smolting and migration to the marine environment at a much smaller size than occurs in the other populations. Therefore, genetic data were consistent with a scenario where isolation by adaptation occurs between populations of Dolly Varden in the study system, and ecological data suggest that this isolation may partially be a result of a novel Dolly Varden life history of seawater tolerance at a smaller size than previously recognized.