Environmental adaptations as windows on molecular evolution

Phylogenomic analysis of the bowfin (Amia calva) reveals unrecognized species diversity in a living fossil lineage

The Bowfin (Amia calva), as currently recognized, represents the sole living member of the family Amiidae, which dates back to approximately 150 Ma. Prior to 1896, 13 species of extant Bowfins had been described, but these were all placed into a single species with no rationale or analysis given. This situation has persisted until the present day, with little attention given to re-evaluation of those previously described nominal forms. Here, we present a phylogenomic analysis based on over 21,000 single nucleotide polymorphisms (SNPs) from 94 individuals that unambiguously demonstrates the presence of at least two independent evolutionary lineages within extant Amia populations that merit species-level standing, as well as the possibility of two more. These findings not only expand the recognizable species diversity in an iconic, ancient lineage, but also demonstrate the utility of such methods in addressing previously intractable questions of molecular systematics and phylogeography in slowly evolving groups of ancient fishes.

Temperature effects on teleost immunity in the light of climate change

Temperature is an important environmental modulator of teleost immune activity. Susceptibility of teleosts to temperature variation depends on the species-specific adaptive temperature range, and the activity of the teleost immune system is generally temperature-dependent. Similar to many physiological and metabolic traits of ectotherms, temperature modulates the activity of immune traits. At low temperatures, acquired immunity of many teleost species is down-modulated, and their immuno-competence mainly depends on innate immunity. At intermediate temperatures, both innate and acquired immunity are fully active and provide optimal protection, including long-lasting immunological memory. When temperatures increase and reach the upper permissive range, teleost immunity is compromised. Moreover, temperature shifts may have negative effects on teleost immune functions, in particular if shifts occur rapidly with high amplitudes. On the contrary, short-term temperature increase may help teleost immunity to fight against pathogens transiently. A major challenge to teleosts therefore is to maintain immuno-competence throughout the temperature range they are exposed to. Climate change coincides with rising temperatures, and more frequent and more extreme temperature shifts. Both are likely to influence the immuno-competence of teleosts. Nonetheless, teleosts exist in habitats that differ substantially in temperature, ranging from below zero in the Arctic's to above 40°C in warm springs, illustrating their enormous potential to adapt to different temperature regimes. The present review seeks to discuss how changes in temperature variation, induced by climate change, might influence teleost immunity.

High migratory propensity constitutes a single stock of an exploited cutlassfish species in the Northwest Pacific: A microsatellite approach

Cutlassfishes, also known as hairtails, include multiple predatory fishes of the family Trichiuridae. They constitute a top marine fish commodity globally, yet the knowledge about their composition and intraspecific genetic structures is still limited. Trichiurus japonicus accounts for a major amount in the northwest Pacific fishery. Previous studies based on mitochondrial DNA markers reported incongruences in its population structure, hence prompting the need for high-resolution markers and avoiding possible shortcomings in its management. Here we genotyped ten novel de novo-assembled transcriptome-derived microsatellite markers on a total of 150 samples across five major fishing grounds (encompassing latitudes 22–39°N). These markers presented a high number of alleles and heterozygosity compared to other marine fishes, corresponding to the large effective population size of ~20,000 per location and cohort differentiation. Population structuring analyses suggested T. japonicus to be a homogenous well-mixed population. This configuration is likely attributed to the majority of its effective population migrates across locations, and the absence of oceanographic barriers at the continental shelves. Qingdao with reportedly high ocean productivity could be a genetic pseudosink based on the high heterozygosity and migratory preference. Moreover, the results of sign tests suggest that T. japonicus experienced a recent bottleneck likely concurrent with historical glaciation events. Further, we demonstrated satisfactory cross-amplifications of our markers on several congeners, indicating a great promise to use these markers to study the population genetics of trichiurids. Together, our findings will serve as an essential groundwork for enhancing resource conservation and management of cutlassfishes.

Ion Transporters and Osmoregulation in the Kidney of Teleost Fishes as a Function of Salinity

Euryhaline teleosts exhibit major changes in renal function as they move between freshwater (FW) and seawater (SW) environments, thus tolerating large fluctuations in salinity. In FW, the kidney excretes large volumes of water through high glomerular filtration rates (GFR) and low tubular reabsorption rates, while actively reabsorbing most ions at high rates. The excreted product has a high urine flow rate (UFR) with a dilute composition. In SW, GFR is greatly reduced, and the tubules reabsorb as much water as possible, while actively secreting divalent ions. The excreted product has a low UFR, and is almost isosmotic to the blood plasma, with Mg²⁺, SO₄^2–, and Cl^– as the major ionic components. Early studies at the organismal level have described these basic patterns, while in the last two decades, studies of regulation at the cell and molecular level have been implemented, though only in a few euryhaline groups (salmonids, eels, tilapias, and fugus). There have been few studies combining the two approaches. The aim of the review is to integrate known aspects of renal physiology (reabsorption and secretion) with more recent advances in molecular water and solute physiology (gene and protein function of transporters). The renal transporters addressed include the subunits of the Na⁺, K⁺- ATPase (NKA) enzyme, monovalent ion transporters for Na⁺, Cl^–, and K⁺ (NKCC1, NKCC2, CLC-K, NCC, ROMK2), water transport pathways [aquaporins (AQP), claudins (CLDN)], and divalent ion transporters for SO₄^2–, Mg²⁺, and Ca²⁺ (SLC26A6, SLC26A1, SLC13A1, SLC41A1, CNNM2, CNNM3, NCX1, NCX2, PMCA). For each transport category, we address the current understanding at the molecular level, try to synthesize it with classical knowledge of overall renal function, and highlight knowledge gaps. Future research on the kidney of euryhaline fishes should focus on integrating changes in kidney reabsorption and secretion of ions with changes in transporter function at the cellular and molecular level (gene and protein verification) in different regions of the nephrons. An increased focus on the kidney individually and its functional integration with the other osmoregulatory organs (gills, skin and intestine) in maintaining overall homeostasis will have applied relevance for aquaculture.

Transcriptomic evidences of local thermal adaptation for the native fish Colossoma macropomum (Cuvier, 1818)

Brazil has five climatically distinct regions, with an annual average temperature difference up to 14 ºC between the northern and southern extremes. Environmental variation of this magnitude can lead to new genetic patterns among farmed fish populations. Genetically differentiated populations of tambaqui (Colossoma macropomum Cuvier, 1818), an important freshwater fish for Brazilian continental aquaculture, may be associated with regional adaptation. In this study, we selected tambaquis raised in two thermally distinct regions, belonging to different latitudes, to test this hypothesis. De novo transcriptome analysis was performed to compare the significant differences of genes expressed in the liver of juvenile tambaqui from a northern population (Balbina) and a southeastern population (Brumado). In total, 2,410 genes were differentially expressed (1,196 in Balbina and 1,214 in Brumado). Many of the genes are involved in a multitude of biological functions such as biosynthetic processes, homeostasis, biorhythm, immunity, cell signaling, ribosome biogenesis, modification of proteins, intracellular transport, structure/cytoskeleton, and catalytic activity. Enrichment analysis based on biological networks showed a different protein interaction profile for each population, whose encoding genes may play potential functions in local thermal adaptation of fish to their respective farming environments.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

Morphological and physiological traits of Mediterranean sticklebacks living in the Camargue wetland (Rhone river delta)

Three-spined sticklebacks (Gasterosteus aculeatus L.) living at the southern limit of the species distribution range could possess specific morphological and physiological traits that enable these fish to live at the threshold of their physiological capacities. Morphological analysis was carried out on samples of sticklebacks living in different saline habitats of the Camargue area (Rhone delta, northern Mediterranean coast) obtained from 1993 to 2017. Salinity acclimation capacities were also investigated using individuals from freshwater-low salinity drainage canals and from mesohaline-euryhaline lagoons. Fish were maintained in laboratory conditions at salinity values close to those of their respective habitats: low salinity (LS, 5‰) or seawater (SW, 30‰). Fish obtained from a mesohaline brackish water lagoon (BW, 15‰) were acclimated to SW or LS. Oxygen consumption rates and branchial Na⁺ /K⁺ -ATPase (NKA) activity (indicator of fish osmoregulatory capacity) were measured in these LS or SW control fish and in individuals subjected to abrupt SW or LS transfers. At all the studied locations, only the low-plated "leiurus" morphotype showed no spatial or temporal variations in their body morphology. Gill rakers were only longer and denser in fish sampled from the LS-freshwater (FW) drainage canals. All fish presented similar physiological capacities. Oxygen consumption rates were not influenced by salinity challenge except in SW fish transferred to LS immediately and 1 h after transfer. However, and as expected, gill NKA activity was salinity dependent. Sticklebacks of the Camargue area sampled from habitats with contrasted saline conditions are homogenously euryhaline, have low oxygen consumption rates and do not appear to experience significantly greater metabolic costs when challenged with salinity. However, an observed difference in gill raker length and density is most probably related to the nutritional condition of their habitat, indicating that individuals can rapidly acclimatize to different diets.

Impact of pyrethroids and organochlorine pesticides residue on IGF-1 and CYP1A genes expression and muscle protein patterns of cultured Mugil capito

This study aimed to assess the levels of pyrethroids and organochlorine residues in the tissues of cultured Mugil capito and in water samples obtained from three different sites (Al-Hamol, Al-Riad and Sidi Salem; referred to as Area 1, Area 2, and Area 3, respectively) in the Delta region, Egypt. The study also assessed the biochemical markers in exposed mullet and evaluated the impact of these residues on the expression of insulin-like growth factor 1 (IGF-1) in muscle and cytochrome P4501A (CYP1A) in liver tissues using qRT-PCR and SDS-PAGE methods. The results revealed that pesticide residue levels in the water were variable, but were lower than detected levels in fish. Significant (P < 0.05) differences were found across the three study areas in terms of serum ALT, but the serum AST level was not significantly (P > 0.05) elevated in all study regions. Serum creatinine and urea levels were significantly (P < 0.05) elevated in area 3. Furthermore, glutathione and malondialdehyde concentrations significantly increased (P < 0.05) in liver tissues in area 3. Using the qRT-PCR technique, the results revealed that the expression level of IGF-1 was most significant in area 3, while the expression level of CYP1A was most significant in area 1. The protein profile showed some differences in band numbers and molecular weights of protein bands across different regions. Overall, the alteration in biochemical parameters revealed pesticide interference with the metabolic processes of fish. Furthermore, the pesticide pollution had an effect on the expression of IGF-1 and CYP1A genes and led to changes in the protein profile. Therefore, these markers can be used to monitor fish distress following exposure to the pollutant.

Thermal acclimation and subspecies-specific effects on heart and brain mitochondrial performance in a eurythermal teleost (Fundulus heteroclitus)

Mitochondrial performance may play a role in setting whole-animal thermal tolerance limits and their plasticity, but the relative roles of adjustments in mitochondrial performance across different highly aerobic tissues remain poorly understood. We compared heart and brain mitochondrial responses to acute thermal challenges and to thermal acclimation using high-resolution respirometry in two locally adapted subspecies of Atlantic killifish (Fundulus heteroclitus). We predicted that 5°C acclimation to would result in compensatory increases in mitochondrial performance, while 33°C acclimation would cause suppression of mitochondrial function to minimize the effects of high temperature on mitochondrial metabolism. In contrast, acclimation to both 33 and 5°C decreased mitochondrial performance compared to fish acclimated to 15°C. These adjustments could represent an energetic cost saving mechanism at temperature extremes. Acclimation responses were similar in both heart and brain; however, this effect was smaller in the heart which might indicate its importance in maintaining whole-animal thermal performance. Alternatively, larger acclimation effects in the brain might indicate greater thermal sensitivity compared to the heart. We detected only modest differences between subspecies that were dependent on the tissue assayed. These data demonstrate extensive plasticity in mitochondrial performance following thermal acclimation in killifish, and indicate that the extent of these responses differs between tissues, highlighting the importance and complexity of mitochondrial regulation in thermal acclimation in eurytherms.

Kidney morphology and candidate gene expression shows plasticity in sticklebacks adapted to divergent osmotic environments

Novel physiological challenges in different environments can promote the evolution of divergent phenotypes, either through plastic or genetic changes. Environmental salinity serves as a key barrier to the distribution of nearly all aquatic organisms, and species diversification is likely to be enabled by adaptation to alternative osmotic environments. The threespine stickleback (Gasterosteus aculeatus) is a euryhaline species with populations found both in marine and freshwater environments. It has evolved both highly plastic and locally adapted phenotypes due to salinity-derived selection, but the physiological and genetic basis of adaptation to salinity is not fully understood. We integrated comparative cellular morphology of the kidney, a key organ for osmoregulation, and candidate gene expression to explore the underpinnings of evolved variation in osmotic plasticity within two populations of sticklebacks from distinct salinity zones in the Baltic Sea: the high salinity Kattegat, representative of the ancestral marine habitat, and the low salinity Bay of Bothnia. A common-garden experiment revealed that kidney morphology in the ancestral high salinity population had a highly plastic response to salinity conditions, whereas this plastic response was reduced in the low salinity population. Candidate gene expression in kidney tissue revealed a similar pattern of population-specific differences, with a higher degree of plasticity in the native high salinity population. Together these results suggest that renal cellular morphology has become canalized to low salinity, and that these structural differences may have functional implications for osmoregulation.

Comparative enzymology-new insights from studies of an "old" enzyme, lactate dehydrogenase

Comparative enzymology explores the molecular mechanisms that alter the properties of enzymes to best fit and adapt them to the biotic demands and abiotic stresses that affect the cellular environment in which these protein catalysts function. For many years, comparative enzymology was primarily concerned with analyzing enzyme functional properties (e.g. substrate affinities, allosteric effectors, responses to temperature or pH, stabilizers, denaturants, etc.) in order to determine how enzyme properties were optimized to function under changing conditions. More recently it became apparent that posttranslational modifications of enzymes play a huge role in metabolic regulation. At first, such modifications appeared to target just crucial regulatory enzymes but recent work is showing that many dehydrogenases are also targets of posttranslational modification leading to substantial changes in enzyme properties. The present article focuses in particular on lactate dehydrogenase (LDH) showing that stress-induced changes in enzyme properties can be linked with reversible posttranslational modifications; e.g. changes in the phosphorylation state of LDH occur in response to dehydration stress in frogs and anoxia exposure of turtles and snails. Furthermore, these studies show that LDH is also a target of other posttranslational modifications including acetylation, methylation and ubiquitination that change in response to anoxia or dehydration stress. Selected new methods for exploring posttranslational modifications of dehydrogenases are discussed and new challenges for the future of comparative enzymology are presented that will help to achieve a deeper understanding of biochemical adaptation through enzyme regulation.

Transcriptional response of stress-regulated genes to cadmium exposure in the cockle Cerastoderma glaucum from the gulf of Gabès area (Tunisia)

This study investigates cadmium effects on key messenger RNA (mRNA) expression (MT, MnSOD, CuZnSOD, CAT, ABCB1, HSP70, and CO1) by qPCR in the cockle Cerastoderma glaucum after chronic exposure to two high but environmentally relevant concentrations of CdCl2 (50 μg/L and 5 mg/L) for 12 h to 18 days. Cd accumulation measured in cockles' tissues is significantly higher in both treatment conditions compared to controls and in a dose-dependent manner. Stress on stress tests performed at different times of the experiment clearly demonstrated that exposure to both concentrations of Cd significantly affects cockle survival time in air. Important changes in gene transcription were also highlighted. In particular, MT, HSP70, CAT, and CuZnSOD seem to be relevant biomarkers of Cd exposure because (1) their mRNA levels increase upon exposure and (2) they are highly correlated to Cd accumulation in tissues. Results may be useful for control strategies and for the use of cockles as sentinel organisms.

The Elizabeth River Story: A Case Study in Evolutionary Toxicology

The Elizabeth River system is an estuary in southeastern Virginia, surrounded by the towns of Chesapeake, Norfolk, Portsmouth, and Virginia Beach. The river has played important roles in U.S. history and has been the location of various military and industrial activities. These activities have been the source of chemical contamination in this aquatic system. Important industries, until the 1990s, included wood treatment plants that used creosote, an oil-derived product that is rich in polycyclic aromatic hydrocarbons (PAH). These plants left a legacy of PAH pollution in the river, and in particular Atlantic Wood Industries is a designated Superfund site now undergoing remediation. Numerous studies examined the distribution of PAH in the river and impacts on resident fauna. This review focuses on how a small estuarine fish with a limited home range, Fundulus heteroclitus (Atlantic killifish or mummichog), has responded to this pollution. While in certain areas of the river this species has clearly been impacted, as evidenced by elevated rates of liver cancer, some subpopulations, notably the one associated with the Atlantic Wood Industries site, displayed a remarkable ability to resist the marked effects PAH have on the embryonic development of fish. This review provides evidence of how pollutants have acted as evolutionary agents, causing changes in ecosystems potentially lasting longer than the pollutants themselves. Mechanisms underlying this evolved resistance, as well as mechanisms underlying the effects of PAH on embryonic development, are also described. The review concludes with a description of ongoing and promising efforts to restore this historic American river.

Synaptosomal lactate dehydrogenase isoenzyme composition is shifted toward aerobic forms in primate brain evolution

With the evolution of a relatively large brain size in haplorhine primates (i.e. tarsiers, monkeys, apes, and humans), there have been associated changes in the molecular machinery that delivers energy to the neocortex. Here we investigated variation in lactate dehydrogenase (LDH) expression and isoenzyme composition of the neocortex and striatum in primates using quantitative Western blotting and isoenzyme analysis of total homogenates and synaptosomal fractions. Analysis of isoform expression revealed that LDH in synaptosomal fractions from both forebrain regions shifted towards a predominance of the heart-type, aerobic isoform LDH-B among haplorhines as compared to strepsirrhines (i.e. lorises and lemurs), while in the total homogenate of the neocortex and striatum there was no significant difference in LDH isoenzyme composition between the primate suborders. The largest increase occurred in synapse-associated LDH-B expression in the neocortex, with an especially remarkable elevation in the ratio of LDH-B/LDH-A in humans. The phylogenetic variation in the ratio of LDH-B/LDH-A was correlated with species-typical brain mass but not the encephalization quotient. A significant LDH-B increase in the subneuronal fraction from haplorhine neocortex and striatum suggests a relatively higher rate of aerobic glycolysis that is linked to synaptosomal mitochondrial metabolism. Our results indicate that there is a differential composition of LDH isoenzymes and metabolism in synaptic terminals that evolved in primates to meet increased energy requirements in association with brain enlargement.

The Adaptive Cline at LDH (Lactate Dehydrogenase) in Killifish Fundulus heteroclitus Remains Stationary After 40 Years of Warming Estuaries

Since the 1970s, water temperatures along the Atlantic seaboard of the United States have risen by an average of 0.5 °C in summer months and 2.2 °C in winter months. In response, the distribution and abundance of several nearshore species have changed dramatically, but no study has attempted to document whether estuarine populations have evolved greater thermal tolerance. Here, we re-examine the classic latitudinal cline at lactate dehydrogenase (LDH) in the killifish Fundulus heteroclitus that was originally described by Dennis Powers and associates from samples collected between 1970 and 1972. Laboratory and field evidences indicated that northern and southern isozymes at muscle LDH are locally adapted to cold and warm temperatures, respectively. Despite the potential for evolutionary response at this adaptive locus, we detected no significant shift of the LDH cline from 20 to 30 F. heteroclitus collected at each of 13 locations between the early 1970s and 2010. We conclude that the microevolution of LDH-mediated thermal tolerance has not occurred, that shifts in alleles are too incremental to be distinguished from random processes, or that F. heteroclitus uses phenotypic and genetic mechanisms besides LDH to respond to warmer waters.

Genetic variation at aryl hydrocarbon receptor (AHR) loci in populations of Atlantic killifish (Fundulus heteroclitus) inhabiting polluted and reference habitats

Background

The non-migratory killifish Fundulus heteroclitus inhabits clean and polluted environments interspersed throughout its range along the Atlantic coast of North America. Several populations of this species have successfully adapted to environments contaminated with toxic aromatic hydrocarbon pollutants such as polychlorinated biphenyls (PCBs). Previous studies suggest that the mechanism of resistance to these and other “dioxin-like compounds” (DLCs) may involve reduced signaling through the aryl hydrocarbon receptor (AHR) pathway. Here we investigated gene diversity and evidence for positive selection at three AHR-related loci (AHR1, AHR2, AHRR) in F. heteroclitus by comparing alleles from seven locations ranging over 600 km along the northeastern US, including extremely polluted and reference estuaries, with a focus on New Bedford Harbor (MA, USA), a PCB Superfund site, and nearby reference sites.

Results

We identified 98 single nucleotide polymorphisms within three AHR-related loci among all populations, including synonymous and nonsynonymous substitutions. Haplotype distributions were spatially segregated and F-statistics suggested strong population genetic structure at these loci, consistent with previous studies showing strong population genetic structure at other F. heteroclitus loci. Genetic diversity at these three loci was not significantly different in contaminated sites as compared to reference sites. However, for AHR2 the New Bedford Harbor population had significant F_ST values in comparison to the nearest reference populations. Tests for positive selection revealed ten nonsynonymous polymorphisms in AHR1 and four in AHR2. Four nonsynonymous SNPs in AHR1 and three in AHR2 showed large differences in base frequency between New Bedford Harbor and its reference site. Tests for isolation-by-distance revealed evidence for non-neutral change at the AHR2 locus.

Conclusion

Together, these data suggest that F. heteroclitus populations in reference and polluted sites have similar genetic diversity, providing no evidence for strong genetic bottlenecks for populations in polluted locations. However, the data provide evidence for genetic differentiation among sites, selection at specific nucleotides in AHR1 and AHR2, and specific AHR2 SNPs and haplotypes that are associated with the PCB-resistant phenotype in the New Bedford Harbor population. The results suggest that AHRs, and especially AHR2, may be important, recurring targets for selection in local adaptation to dioxin-like aromatic hydrocarbon contaminants.

Microsatellite length variation in candidate genes correlates with habitat in the gilthead sea bream Sparus aurata

The genetic basis and evolutionary implications of local adaptation in high gene flow marine organisms are still poorly understood. In several Mediterranean fish species, alternative migration patterns exist between individuals entering coastal lagoons that offer favourable conditions for growth and those staying in the sea where environmental conditions are less subject to rapid and stressful change. Whether these coexisting strategies are phenotypically plastic or include a role for local adaptation through differential survival needs to be determined. Here, we explore the genetic basis of alternate habitat use in western Mediterranean populations of the gilthead sea bream (Sparus aurata). Samples from lagoonal and open-sea habitats were typed for three candidate gene microsatellite loci, seven anonymous microsatellites and 44 amplified fragment length polymorphism markers to test for genotype-environment associations. While anonymous markers globally indicated high levels of gene flow across geographic locations and habitats, non-neutral differentiation patterns correlated with habitat type were found at two candidate microsatellite loci located in the promoter region of the growth hormone and prolactin genes. Further analysis of these two genes revealed that a mechanism based on habitat choice alone could not explain the distribution of genotype frequencies at a regional scale, thus implying a role for differential survival between habitats. We also found an association between allele size and habitat type, which, in the light of previous studies, suggests that polymorphisms in the proximal promoter region could influence gene expression by modulating transcription factor binding, thus providing a potential explanatory link between genotype and growth phenotype in nature.

Can diversifying selection be distinguished from history in geographic clines? A population genomic study of killifish (Fundulus heteroclitus)

A common geographical pattern of genetic variation is the one-dimensional cline. Clines may be maintained by diversifying selection across a geographical gradient but can also reflect historical processes such as allopatry followed by secondary contact. To identify loci that may be undergoing diversifying selection, we examined the distribution of geographical variation patterns across the range of the killifish (Fundulus heteroclitus) in 310 loci, including microsatellites, allozymes, and single nucleotide polymorphisms. We employed two approaches to detect loci under strong diversifying selection. First, we developed an automated method to identify clinal variation on a per-locus basis and examined the distribution of clines to detect those that exhibited signifcantly steeper slopes. Second, we employed a classic [Formula: see text]-outlier method as a complementary approach. We also assessed performance of these techniques using simulations. Overall, latitudinal clines were detected in nearly half of all loci genotyped (i.e., all eight microsatellite loci, 12 of 16 allozyme loci and 44% of the 285 SNPs). With the exception of few outlier loci (notably mtDNA and malate dehydrogenase), the positions and slopes of Fundulus clines were statistically indistinguishable. The high frequency of latitudinal clines across the genome indicates that secondary contact plays a central role in the historical demography of this species. Our simulation results indicate that accurately detecting diversifying selection using genome scans is extremely difficult in species with a strong signal of secondary contact; neutral evolution under this history produces clines as steep as those expected under selection. Based on these results, we propose that demographic history can explain all clinal patterns observed in F. heteroclitus without invoking natural selection to either establish or maintain the pattern we observe today.

Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus

BACKGROUND

Geographic variation in the thermal environment impacts a broad range of biochemical and physiological processes and can be a major selective force leading to local population adaptation. In the intertidal copepod Tigriopus californicus, populations along the coast of California show differences in thermal tolerance that are consistent with adaptation, i.e., southern populations withstand thermal stresses that are lethal to northern populations. To understand the genetic basis of these physiological differences, we use an RNA-seq approach to compare genome-wide patterns of gene expression in two populations known to differ in thermal tolerance.

RESULTS

Observed differences in gene expression between the southern (San Diego) and the northern (Santa Cruz) populations included both the number of affected loci as well as the identity of these loci. However, the most pronounced differences concerned the amplitude of up-regulation of genes producing heat shock proteins (Hsps) and genes involved in ubiquitination and proteolysis. Among the hsp genes, orthologous pairs show markedly different thermal responses as the amplitude of hsp response was greatly elevated in the San Diego population, most notably in members of the hsp70 gene family. There was no evidence of accelerated evolution at the sequence level for hsp genes. Among other sets of genes, cuticle genes were up-regulated in SD but down-regulated in SC, and mitochondrial genes were down-regulated in both populations.

CONCLUSIONS

Marked changes in gene expression were observed in response to acute sub-lethal thermal stress in the copepod T. californicus. Although some qualitative differences were observed between populations, the most pronounced differences involved the magnitude of induction of numerous hsp and ubiquitin genes. These differences in gene expression suggest that evolutionary divergence in the regulatory pathway(s) involved in acute temperature stress may offer at least a partial explanation of population differences in thermal tolerance observed in Tigriopus.

Phylogeography study of Ammodytes personatus in Northwestern Pacific: Pleistocene isolation, temperature and current conducted secondary contact

To assess the role of historical process and contemporary factors in shaping population structures in Northwestern Pacific, mitochondrial control region sequences were analyzed to characterize the phylogeography and population structure of the Japanese sand lance Ammodytes personatus. A total of 429 individuals sampled from 17 populations through the species' range are sequenced. Two distinct lineages are detected, which might have been divergent in the Sea of Japan and Pacific costal waters of Japanese Island, during the low sea level. Significant genetic structure is revealed between the Kuroshio and Oyashio Currents. However, significant genetic structure is also detected in the Sea of Japan, contracting expected homogenization hypothesis in Tsushima Current. The haplotype frequency of lineages in both sides of Japanese Island and significant genetic structure between north and south groups revealed that the distribution of lineage B and north group were highly limited by the annual sea temperature. The lack of lineage B in Qingdao population with low sea temperature reflects the sea temperature barrier. Lack of genetic structure in the south group and north group populations indicated that ocean currents within groups facilitated the dispersal of A. personatus.

Intraspecific variation in the thermal plasticity of mitochondria in killifish

Populations of the Atlantic killifish (Fundulus heteroclitus) inhabit salt marshes and estuaries along the eastern coast of North America from Newfoundland to northern Florida, and are thus exposed to a large range of temperatures. Previous studies have shown higher whole-organism metabolic rates in the northern subspecies (F. h. macrolepidotus) compared with the southern subspecies (F. h. heteroclitus) of these fish. Here, we examine phenotypic plasticity in the response to cold temperatures between the two subspecies by acclimating fish to 5, 15 and 25°C and comparing several mitochondrial and muscle properties. The relative area of oxidative muscle versus glycolytic muscle fibers was greater in the northern subspecies at the 5 and 15°C acclimation temperatures. However, there were no differences in capillary density between the two subspecies or at different temperatures. Mitochondrial volume and surface densities increased in response to cold temperature acclimation in red and white muscle, but only in the northern killifish. Citrate synthase activities also increased in the northern killifish at 5 and 15°C. The ratio of calculated [free ADP] to [ATP] increased in the 5°C acclimated southern killifish but not in the northern killifish at 5°C when compared with the 15°C acclimation group, suggesting that there are differences in adenylate signaling for mitochondrial respiration between subspecies at low temperature. Taken together, our data indicate that the northern subspecies have a greater ability to increase mitochondrial capacity at colder temperatures compared with the southern subspecies, providing one of the few examples of intraspecific variation in phenotypic plasticity in mitochondrial amount in response to cold temperatures.

Differences in salinity tolerance and gene expression between two populations of Atlantic cod (Gadus morhua) in response to salinity stress

Populations of marine fish, even from contrasting habitats, generally show low genetic differentiation at neutral genetic markers. Nevertheless, there is increasing evidence for differences in gene expression among populations that may be ascribed to adaptive divergence. Studying variation in salinity tolerance and gene expression among Atlantic cod (Gadus morhua) from two populations distributed across a steep salinity gradient, we observed high mortality (45% North Sea cod and 80% Baltic Sea cod) in a reciprocal common garden setup. Quantitative RT-PCR assays for expression of hsp70 and Na/K-ATPase α genes demonstrated significant differences in gene regulation within and between populations and treatment groups despite low sample sizes. Most interesting are the significant differences observed in expression of the Na/K-ATPase α gene in gill tissue between North Sea and Baltic cod. The findings strongly suggest that Atlantic cod are adapted to local saline conditions, despite relatively low levels of neutral genetic divergence between populations.

Evolutionary factors affecting Lactate dehydrogenase A and B variation in the Daphnia pulex species complex

BACKGROUND

Evidence for historical, demographic and selective factors affecting enzyme evolution can be obtained by examining nucleotide sequence variation in candidate genes such as Lactate dehydrogenase (Ldh). Two closely related Daphnia species can be distinguished by their electrophoretic Ldh genotype and habitat. Daphnia pulex populations are fixed for the S allele and inhabit temporary ponds, while D. pulicaria populations are fixed for the F allele and inhabit large stratified lakes. One locus is detected in most allozyme surveys, but genome sequencing has revealed two genes, LdhA and LdhB.

RESULTS

We sequenced both Ldh genes from 70 isolates of these two species from North America to determine if the association between Ldh genotype and habitat shows evidence for selection, and to elucidate the evolutionary history of the two genes. We found that alleles in the pond-dwelling D. pulex and in the lake-dwelling D. pulicaria form distinct groups at both loci, and the substitution of Glutamine (S) for Glutamic acid (F) at amino acid 229 likely causes the electrophoretic mobility shift in the LDHA protein. Nucleotide diversity in both Ldh genes is much lower in D. pulicaria than in D. pulex. Moreover, the lack of spatial structuring of the variation in both genes over a wide geographic area is consistent with a recent demographic expansion of lake populations. Neutrality tests indicate that both genes are under purifying selection, but the intensity is much stronger on LdhA.

CONCLUSIONS

Although lake-dwelling D. pulicaria hybridizes with the other lineages in the pulex species complex, it remains distinct ecologically and genetically. This ecological divergence, coupled with the intensity of purifying selection on LdhA and the strong association between its genotype and habitat, suggests that experimental studies would be useful to determine if variation in molecular function provides evidence that LDHA variants are adaptive.

Physiological adaptations of small mammals to desert ecosystems

Adaptations of animals to the xeric environment have been studied in various taxonomic groups and across several deserts. Despite the impressive data that have been accumulated, the focus in most of these studies is mainly on the significance of one variable at a time. Here, we attempt to integrate between responses of several physiological systems, challenged by increasing diet and water salinity and extreme temperatures, acquired in different studies of thermo and osmo-regulatory adaptations, of small rodents, to the xeric environment. Studies have shown differential thermoregulatory responses to increased dietary salinity, which were attributed to habitat and habits of the relevant species. In the thermoregulatory studies, a potential adaptive significance of low metabolic rate was demonstrated. From an evolutionary point of view, the most important adaptation is in the timing of reproduction, as it enables the transfer of genetic properties to the next generation in an unpredictable ecosystem, where reproduction might not occur every year. Results in this aspect show that increased dietary salinity, through an increase in vasopressin plasma levels, plays an important role as a regulator of the reproductive system. We assume that the amount of food existing in the habitat and the amount of reserves in the animal in the form of white adipose tissue are important for reproduction. Photoperiod affects all studied physiological responses, emphasizing the importance of pre-acclimation to seasonal characteristics. We summarize the existing data and suggest neuro-endocrine pathways, which have a central role in these adaptations by affecting thermoregulation, osmoregulation and reproduction to create the optimal response to xeric conditions. These hypotheses can be used as the basis for future studies.

The biochemical basis of life history adaptation: molecular and enzymological causes of NADP(+)-isocitrate dehydrogenase activity differences between morphs of Gryllus firmus that differ in lipid biosynthesis and life history

Although whole-organism aspects of life-history physiology are well studied and molecular information (e.g., transcript abundance) on life-history variation is accumulating rapidly, much less information is available on the biochemical (enzymological) basis of life-history adaptation. The present study investigated the biochemical and molecular causes of specific activity differences of the lipogenic enzyme, NADP(+)-isocitrate dehydrogenase, between genetic lines of the wing-polymorphic cricket, Gryllus firmus, which differ in lipid biosynthesis and life history. With one exception, variation among 21 Nadp(+)-Idh genomic sequences, which spanned the entire coding sequence of the gene, was restricted to a few synonymous substitutions within and among replicate flight-capable or flightless lines. No NADP(+)-IDH electromorph variation was observed among individuals within or among lines as determined by polyacrylamide gel electrophoresis. Nor did any NADP(+)-IDH kinetic or stability parameter, such as K(M) for substrate or cofactor, k(cat), or thermal denaturation, differ between flight-capable and flightless lines. By contrast, line differences in NADP(+)-IDH-specific activity strongly covaried with transcript abundance and enzyme protein concentration. These results demonstrate that NADP(+)-IDH-specific activity differences between artificially selected lines of G. firmus are due primarily, if not exclusively, to genetic variation in regulators of NADP(+)-IDH gene expression, with no observed contribution from altered catalytic efficiency of the enzyme due to changes in amino acid sequence or posttranslational modification. Kinetic analyses indicate that in vitro differences in enzyme-specific activity between flight-capable and flightless lines likely occur in vivo. This study constitutes the most comprehensive analysis to date of the biochemical and molecular causes of naturally occurring genetic variation in enzyme activity that covaries strongly with life history.

Transcriptomic assessment of resistance to effects of an aryl hydrocarbon receptor (AHR) agonist in embryos of Atlantic killifish (Fundulus heteroclitus) from a marine Superfund site

Background

Populations of Atlantic killifish (Fundulus heteroclitus) have evolved resistance to the embryotoxic effects of polychlorinated biphenyls (PCBs) and other halogenated and nonhalogenated aromatic hydrocarbons that act through an aryl hydrocarbon receptor (AHR)-dependent signaling pathway. The resistance is accompanied by reduced sensitivity to induction of cytochrome P450 1A (CYP1A), a widely used biomarker of aromatic hydrocarbon exposure and effect, but whether the reduced sensitivity is specific to CYP1A or reflects a genome-wide reduction in responsiveness to all AHR-mediated changes in gene expression is unknown. We compared gene expression profiles and the response to 3,3',4,4',5-pentachlorobiphenyl (PCB-126) exposure in embryos (5 and 10 dpf) and larvae (15 dpf) from F. heteroclitus populations inhabiting the New Bedford Harbor, Massachusetts (NBH) Superfund site (PCB-resistant) and a reference site, Scorton Creek, Massachusetts (SC; PCB-sensitive).

Results

Analysis using a 7,000-gene cDNA array revealed striking differences in responsiveness to PCB-126 between the populations; the differences occur at all three stages examined. There was a sizeable set of PCB-responsive genes in the sensitive SC population, a much smaller set of PCB-responsive genes in NBH fish, and few similarities in PCB-responsive genes between the two populations. Most of the array results were confirmed, and additional PCB-regulated genes identified, by RNA-Seq (deep pyrosequencing).

Conclusions

The results suggest that NBH fish possess a gene regulatory defect that is not specific to one target gene such as CYP1A but rather lies in a regulatory pathway that controls the transcriptional response of multiple genes to PCB exposure. The results are consistent with genome-wide disruption of AHR-dependent signaling in NBH fish.

Integrating evolutionary and functional approaches to infer adaptation at specific loci

Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally,population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation.

Gene expression analysis for the identification of selection and local adaptation in fishes

In recent years, variation in gene expression has been recognized as an important component of environmental adaptation in multiple model species, including a few fish species. There is, however, still little known about the genetic basis of adaptation in gene expression resulting from variation in the aquatic environment (e.g. temperature, salinity and oxygen) and the physiological effect and costs of such differences in gene expression. This review presents and discusses progress and pitfalls of applying gene expression analyses to fishes and suggests simple frameworks to get started with gene expression analysis. It is emphasized that well-planned gene expression studies can serve as an important tool for the identification of selection in local populations of fishes, even for non-traditional model species where limited genomic information is available. Recent studies focusing on gene expression variation among natural fish populations are reviewed, highlighting the latest applications that combine genetic evidence from neutral markers and gene expression data.

Population genetic structure of the point-head flounder, Cleisthenes herzensteini, in the Northwestern Pacific

Intraspecific phylogenies can provide useful insights into how populations have been shaped by historical and contemporary processes. To determine the population genetic structure and the demographic and colonization history of Cleisthenes herzensteini in the Northwestern Pacific, one hundred and twenty-one individuals were sampled from six localities along the coastal regions of Japan and the Yellow Sea of China. Mitochondrial DNA variation was analyzed using DNA sequence data from the 5' end of control region. High levels of haplotype diversity (>0.96) were found for all populations, indicating a high level of genetic diversity. No pattern of isolation by distance was detected among the population differentiation throughout the examined range. Analyses of molecular variance (AMOVA) and the conventional population statistic Fst revealed no significant population genetic structure among populations. According to the exact test of differentiation among populations, the null hypothesis that C. herzensteini within the examined range constituted a non-differential mtDNA gene pool was accepted. The demographic history of C. herzensteini was examined using neutrality test and mismatch distribution analyses and results indicated Pleistocene population expansion (about 94-376 kya) in the species, which was consistent with the inference result of nested clade phylogeographical analysis (NCPA) showing contiguous range expansion for C. herzensteini. The lack of phylogeographical structure for the species may reflect a recent range expansion after the glacial maximum and insufficient time to attain migration-drift equilibrium.

Parallel adaptive evolution of Atlantic cod on both sides of the Atlantic Ocean in response to temperature

Despite the enormous economic and ecological importance of marine organisms, the spatial scales of adaptation and biocomplexity remain largely unknown. Yet, the preservation of local stocks that possess adaptive diversity is critical to the long-term maintenance of productive stable fisheries and ecosystems. Here, we document genomic evidence of range-wide adaptive differentiation in a broadcast spawning marine fish, Atlantic cod (Gadus morhua), using a genome survey of single nucleotide polymorphisms. Of 1641 gene-associated polymorphisms examined, 70 (4.2%) tested positive for signatures of selection using a Bayesian approach. We identify a subset of these loci (n=40) for which allele frequencies show parallel temperature-associated clines (p<0.001, r2=0.89) in the eastern and western north Atlantic. Temperature associations were robust to the statistical removal of geographic distance or latitude effects, and contrasted 'neutral' loci, which displayed no temperature association. Allele frequencies at temperature-associated loci were significantly correlated, spanned three linkage groups and several were successfully annotated supporting the involvement of multiple independent genes. Our results are consistent with the evolution and/or selective sweep of multiple genes in response to ocean temperature, and support the possibility of a new conservation paradigm for non-model marine organisms based on genomic approaches to resolving functional and adaptive diversity.

Responses of the European flounder Platichthys flesus to the chemical stress in estuaries: load of contaminants, gene expression, cellular impact and growth rate

European flounder responses to the chemical stress were assessed by a comparative approach on four estuaries displaying contrasted patterns of contamination. The contamination typology of the estuaries was investigated by individual measurements of contaminants in fish. Molecular and physiological responses were studied by gene expression, genotoxicity, neurotoxicity and growth rate. Fishes in contaminated estuaries were characterized by high levels of bioaccumulated contaminants, slow energetic metabolism and reduced growth rate, in contrast to the fish responses in the reference site. A seasonal effect was highlighted for contaminated flounder populations, with high PCB levels, high genotoxicity and elevated detoxification rate in summer compared with winter.

Demographic history and population structure of blackfin flounder (Glyptocephalus stelleri) in Japan revealed by mitochondrial control region sequences

The demographic history and population genetic structure of the blackfin flounder (Glyptocephalus stelleri) along coastal regions of Japan were investigated. Genetic variation in DNA sequences was examined from the first hypervariable region of the mitochondrial DNA control region. A high level of haplotypic diversity (h = 0.99 +/- 0.004) was detected, indicating a high level of intrapopulation genetic diversity. The starburst structure of the minimum spanning tree suggested a very recent origin for most haplotypes. The demographic history of G. stelleri was examined using neutrality tests and mismatch distribution analysis, which also indicated a Pleistocene population expansion at about 124,100-413,400 years ago. Hierarchical molecular variance analysis and conventional population Fst comparisons revealed no significant genetic differentiation throughout the range examined.

Classic approach revitalizes genomics: Complete characterization of a candidate gene for thermal adaptation in two coral reef fishes

Lactate dehydrogenase-B (ldh-b) encodes a metabolic enzyme (LDH-B) which plays an important role in maintaining aerobic performance and in thermal acclimation and/or adaptation of fish. As the first step in understanding the effect this enzyme has on the ability of tropical coral reef fishes to cope with thermal stress, we characterized both coding and non-coding regions of ldh-b in two congeneric perciformes, Plectropomus leopardus and Plectropomus laevis. Ldh-b was 4666 and 4539bp in length in P. leopardus and P. laevis, respectively, with coding regions comprising 1005bp in both species. We report a high level of sequence homology between the coding regions of ldh-b in these two species, with 98.1% identity of nucleotides corresponding to 100% amino acid identity between the deduced protein sequences. Comparison between non-coding (intron) regions of both species revealed the presence of several indels, despite the high level of homology observed (95.9% identity of intron nucleotides). Potential regulatory motifs and elements, including twenty-six simple sequence repeat motifs (mono-, di-, tri- and tetranucleotide) and twenty-three putative microRNA elements are identified within the introns of both species, further supporting recent demonstrations that such short motifs and elements exhibit widespread positioning throughout non-coding regions of the genome. This novel characterization of ldh-b in these two coral reef fishes allows for a wide range of future studies (e.g. analytical comparisons of ldh-b and LDH-B among different fish genera from different thermal environments and habitats).

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

Do mitochondrial properties explain intraspecific variation in thermal tolerance?

As global temperatures rise, there is a growing need to understand the physiological mechanisms that determine an organism's thermal niche. Here, we test the hypothesis that increases in mitochondrial capacity with cold acclimation and adaptation are associated with decreases in thermal tolerance using two subspecies of killifish (Fundulus heteroclitus) that differ in thermal niche. We assessed whole-organism metabolic rate, mitochondrial amount and mitochondrial function in killifish acclimated to several temperatures. Mitochondrial enzyme activities and mRNA levels were greater in fish from the northern subspecies, particularly in cold-acclimated fish, suggesting that the putatively cold-adapted northern subspecies has a greater capacity for increases in mitochondrial amount in response to cold acclimation. When tested at the fish's acclimation temperature, maximum ADP-stimulated (State III) rates of mitochondrial oxygen consumption in vitro were greater in cold-acclimated northern fish than in southern fish but did not differ between subspecies at higher acclimation temperatures. Whole-organism metabolic rate was greater in fish of the northern subspecies at all acclimation temperatures. Cold acclimation also changed the response of mitochondrial respiration to acute temperature challenge. Mitochondrial oxygen consumption was greater in cold-acclimated northern fish than in southern fish at low test temperatures, but the opposite was true at high test temperatures. These differences were reflected in whole-organism oxygen consumption. Our data indicate that the plasticity of mitochondrial function and amount differs between killifish subspecies, with the less high-temperature tolerant, and putatively cold adapted, northern subspecies having greater ability to increase mitochondrial capacity in the cold. However, there were few differences in mitochondrial properties between subspecies at warm acclimation temperatures, despite differences in both whole-organism oxygen consumption and thermal tolerance at these temperatures.

Swimming performance and energetics as a function of temperature in killifish Fundulus heteroclitus

Populations of the common killifish Fundulus heteroclitus are found along a latitudinal temperature gradient in habitats with high thermal variability. The objectives of this study were to assess the effects of temperature and population of origin on killifish swimming performance (assessed as critical swimming speed, U(crit)). Acclimated fish from northern and southern killifish populations demonstrated a wide zone (from 7 degrees to 33 degrees C) over which U(crit) showed little change with temperature, with performance declining significantly only at lower temperatures. Although we observed significant differences in swimming performance between a northern and a southern population of killifish in one experiment, with northern fish having an approximately 1.5-fold-greater U(crit) than southern fish across all acclimation temperatures, we were unable to replicate this finding in other populations or collection years, and performance was consistently high across all populations and at both low (7 degrees C) and high (23 degrees C) acclimation temperatures. The poor swimming performance of southern killifish from a single collection year was correlated with low muscle [glycogen] rather than with other indicators of fuel stores or body condition. Killifish acclimated to 18 degrees C and acutely challenged at temperatures of 5 degrees , 18 degrees , 25 degrees , or 34 degrees C showed modest thermal sensitivity of U(crit) between 18 degrees and 34 degrees C, with performance declining substantially at 5 degrees C. Thus, much of the zone of relative thermal insensitivity of swimming performance is intrinsic in this species rather than acquired as a result of acclimation. These data suggest that killifish are broadly tolerant of changing temperatures, whether acute or chronic, and demonstrate little evidence of local adaptation in endurance swimming performance in populations from different thermal habitats.

Interpopulation differences in expression of candidate genes for salinity tolerance in winter migrating anadromous brown trout (Salmo trutta L.)

Background

Winter migration of immature brown trout (Salmo trutta) into freshwater rivers has been hypothesized to result from physiologically stressful combinations of high salinity and low temperature in the sea.

Results

We sampled brown trout from two Danish populations entering different saline conditions and quantified expression of the hsp70 and Na/K-ATPases α 1b genes following acclimation to freshwater and full-strength seawater at 2°C and 10°C. An interaction effect of low temperature and high salinity on expression of both hsp70 and Na/K-ATPase α 1b was found in trout from the river entering high saline conditions, while a temperature independent up-regulation of both genes in full-strength seawater was found for trout entering marine conditions with lower salinities.

Conclusion

Overall our results support the hypothesis that physiologically stressful conditions in the sea drive sea-run brown trout into freshwater rivers in winter. However, our results also demonstrate intra-specific differences in expression of important stress and osmoregulative genes most likely reflecting adaptive differences between trout populations on a regional scale, thus strongly suggesting local adaptations driven by the local marine environment.

The influence of gene-environment interactions on GHR and IGF-1 expression and their association with growth in brook charr, Salvelinus fontinalis (Mitchill)

Background

Quantitative reaction norm theory proposes that genotype-by-environment interaction (GxE) results from inter-individual differences of expression in adaptive suites of genes in distinct environments. However, environmental norms for actual gene suites are poorly documented. In this study, we investigated the effects of GxE interactions on levels of gene transcription and growth by documenting the impact of rearing environment (freshwater vs. saltwater), sex and genotypic (low vs. high estimated breeding value EBV) effects on the transcription level of insulin-like growth factor (IGF-1) and growth hormone receptor (GHR) in brook charr (Salvelinus fontinalis).

Results

Males grew faster than females (μ_♀ = 1.20 ± 0.07 g·d^-1, μ_♂ = 1.46 ± 0.06 g·d^-1) and high-EBV fish faster than low-EBV fish (μ_LOW = 0.97 ± 0.05 g·d^-1, μ_HIGH = 1.58 ± 0.07 g·d^-1; p < 0.05). However, growth was markedly lower in saltwater-reared fish than freshwater sibs (μ_FW = 1.52 ± 0.07 g·d^-1, μ_SW = 1.15 ± 0.06 g·d^-1), yet GHR mRNA transcription level was significantly higher in saltwater than in freshwater (μ_SW = 0.85 ± 0.05, μ_FW = 0.61 ± 0.05). The ratio of actual growth to units in assayed mRNA ('individual transcript efficiency', iTE; g·d^-1·u^-1) also differed among EBV groups (μ_LOW = 2.0 ± 0.24 g·d^-1·u^-1; μ_HIGH = 3.7 ± 0.24 g·d^-1·u^-1) and environments (μ_SW = 2.0 ± 0.25 g·d^-1·u^-1; μ_FW = 3.7 ± 0.25 g·d^-1·u^-1) for GHR. Males had a lower iTE for GHR than females (μ_♂ = 2.4 ± 0.29 g·d^-1·u^-1; μ_♀ = 3.1 ± 0.23 g·d^-1·u^-1). There was no difference in IGF-1 transcription level between environments (p > 0.7) or EBV groups (p > 0.15) but the level of IGF-1 was four times higher in males than females (μ_♂ = 2.4 ± 0.11, μ_♀ = 0.58 ± 0.09; p < 0.0001). We detected significant sexual differences in iTE (μ_♂ = 1.3 ± 0.59 g·d^-1·u^-1; μ_♀ = 3.9 ± 0.47 g·d^-1·u^-1), salinities (μ_SW = 2.3 ± 0.52 g·d^-1·u^-1; μ_FW = 3.7 ± 0.53 g·d^-1·u^-1) and EBV-groups (μ_LOW = 2.4 ± 0.49 g·d^-1·u^-1; μ_HIGH = 3.8 ± 0.49 g·d^-1·u^-1). Interaction between EBV-group and environment was detected for both GHR (p = 0.027) and IGF-1 (p = 0.019), and for iTE in the two genes (p < 0.0001; p < 0.05, respectively), where increased divergence in levels of GHR and IGF-1 transcription occurred among EBV-groups in the saltwater environment.

Conclusion

Our results show that both environment and sex have major impacts on the expression of mRNA for two key genes involved in the physiological pathway for growth. We also demonstrate for the first time, at least in fish, genotype-by-environment interaction at the level of individual gene transcription. This work contributes significantly to ongoing efforts towards documenting environmentally and sexually induced variance of gene activity and understanding the resulting phenotypes.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

Marine macrophysiology: Studying physiological variation across large spatial scales in marine systems

A new approach toward understanding marine ecosystems has emerged through the integration of ecological physiology and macroecology. This multidisciplinary approach, titled here marine macrophysiology, facilitates unique insight into the foundation of macro-scale ecological patterns, such as biogeographic distributions, via examination of functional attributes of marine organisms across large spatial scales. For example, these broad-scale physiological inquiries confer the ability to directly assess the abundant-center hypothesis (aka Brown's principle) which proposes that species have decreased performance toward their range edges. By extension, the marine macrophysiological perspective also stands to clarify our understanding of more complex macro-scale phenomena such as biological invasions, the design of marine protected areas, and species' responses to global climate change. In this article, we review recent marine macrophysiology research and offer insights into future directions for this emerging field.

Recombinational micro-evolution of functionally different metallothionein promoter alleles from Orchesella cincta

Background

Metallothionein (mt) transcription is elevated in heavy metal tolerant field populations of Orchesella cincta (Collembola). This suggests that natural selection acts on transcriptional regulation of mt in springtails at sites where cadmium (Cd) levels in soil reach toxic values This study investigates the nature and the evolutionary origin of polymorphisms in the metallothionein promoter (pmt) and their functional significance for mt expression.

Results

We sequenced approximately 1600 bp upstream the mt coding region by genome walking. Nine pmt alleles were discovered in NW-European populations. They differ in the number of some indels, consensus transcription factor binding sites and core promoter elements. Extensive recombination events between some of the alleles can be inferred from the alignment. A deviation from neutral expectations was detected in a cadmium tolerant population, pointing towards balancing selection on some promoter stretches. Luciferase constructs were made from the most abundant alleles, and responses to Cd, paraquat (oxidative stress inducer) and moulting hormone were studied in cell lines. By using paraquat we were able to dissect the effect of oxidative stress from the Cd specific effect, and extensive differences in mt induction levels between these two stressors were observed.

Conclusion

The pmt alleles evolved by a number of recombination events, and exhibited differential inducibilities by Cd, paraquat and molting hormone. In a tolerant population from a metal contaminated site, promoter allele frequencies differed significantly from a reference site and nucleotide polymorphisms in some promoter stretches deviated from neutral expectations, revealing a signature of balancing selection. Our results suggest that the structural differences in the Orchesella cincta metallothionein promoter alleles contribute to the metallothionein -over-expresser phenotype in cadmium tolerant populations.