Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

Phylogenetic and phylogeographic insights into Sri Lankan killifishes (Teleostei: Aplocheilidae)

Three nominal species of the killifish genus Aplocheilus are reported from the lowlands of Sri Lanka. Two of these, Aplocheilus dayi and Aplocheilus werneri, are considered endemic to the island, whereas Aplocheilus parvus is reported from both Sri Lanka and Peninsular India. Here, based on a collection from 28 locations in Sri Lanka, also including a dataset of Asian Aplocheilus downloaded from GenBank, we present a phylogeny constructed from the mitochondrial cytochrome b (cytb), mitochondrial cytochrome c oxidase subunit 1 (cox1), and nuclear recombination activating protein 1 (rag1), and investigate the interrelationships of the species of Aplocheilus in Sri Lanka. The endemic Sri Lankan aplocheilid clade comprising A. dayi and A. werneri is recovered as the sister group to the clade comprising A. parvus from Sri Lanka and Aplocheilus blockii from Peninsular India. The reciprocal monophyly of A. dayi and A. werneri is not supported in our molecular phylogeny. A. dayi and A. werneri display strong sexual dimorphism, but species-level differences are subtle, explained mostly by pigmentation patterns. Their phenotypes exhibit a parapatric distribution and may represent locally adapted forms of a single species. Alternatively, the present study does not rule out the possibility that A. dayi and A. werneri may represent an incipient species pair or that they have undergone introgression or hybridization in their contact zones. We provide evidence that the Nilwala-Gin region of southwestern Sri Lanka may have acted as a drought refugium for these fishes.

Population genetic differentiation of Daphnia sinensis in a lasting high-phosphorus Chinese lake, Lake Chaohu

Ecological shifts (e.g., eutrophication) can affect the genetic differentiation of zooplankton populations in lakes. However, the role of environmental change in a lasting high-phosphorus lake driving the genetic differentiation of zooplankton population structure over time is poorly understood. In this paper, the changes of the genetic diversity and differentiation of Daphnia sinensis population were studied by using the mitochondrial COI gene and microsatellite markers on modern groups (from January to June 2016) and historic groups (obtained from resting eggs in the sediments) in Lake Chaohu. Based on the microsatellite markers, six modern groups were clustered into two clusters (the WG cluster and SG cluster) during the seasonal dynamics, whereas the genetic differentiation of the five historic groups showed a wave-like pattern and had evolved into four clusters. Moreover, the haplotype network showed that six modern groups had one origin center whereas five historic groups had two origin centers based on the mitochondrial COI gene marker. Fu’s Fs neutral test and Tajima’s test indicated that the five historic groups deviated from neutral evolution and showed a bottleneck effect in the history process. Water temperature and total dissolved phosphorus were obviously associated with the seasonal genetic differentiation of D. sinensis, whereas nitrogen content of the sediments was significantly related to the long-term microevolution of D. sinensis in the high-phosphorus environment. Therefore, the changing pattern of D. sinensis population genetic structure was one of the environmental selections probably combined with co-evolutionary, where rapid-increasing nitrogen level had a large impact on D. sinensis population genetic structure in lasting high phosphorus environment in Lake Chaohu.

Development and characterization of 16 novel microsatellite markers by Transcriptome sequencing for Angelica dahurica and test for cross-species amplification

BACKGROUND

Angelica dahurica (Apiaceae) is an important herb in traditional Chinese medicine. Because of its important medicinal and economic values, its wild resources were over-exploited and increasingly reduced. Meanwhile, the diversity of cultivars of A. dahurica has decreased as a result of long-term artificial cultivation. However, there are no population genetics studies of natural A. dahurica reported yet, especially for using microsatellite markers (SSRs) to investigate population genetics of the species.

RESULTS

Sixteen polymorphic EST-SSR loci were isolated from A. dahurica with transcriptome sequencing technology (RNA-Seq). The number of alleles varied from 2 to 15 per polymorphic locus over populations with the observed and expected heterozygosities ranging from 0.000 to 1.000 and from 0.000 to 0.829, respectively. Significant deviations from Hardy-Weinberg equilibrium were observed at 8 loci. Tests of linkage disequilibrium showed 11 informative locus pairs were significant across all populations. Cross-species amplification showed that 14 out of 16 SSR loci have the transferability in cultivar-A. dahurica cv. 'Hangbaizhi' and A. decursiva.

CONCLUSIONS

The 16 newly developed loci microsatellite primers with RNA-Seq will be useful for further investigating population genetics of A. dahurica, cultivars and other members of this genus.

Genetic diversity of a widespread annual killifish from coastal Tanzania

Background

African annual killifishes (Nothobranchius spp.) are adapted to seasonally desiccating habitats (ephemeral pools), surviving dry periods as dormant eggs. Given their peculiar life history, geographic aspects of their diversity uniquely combine patterns typical for freshwater taxa (river basin structure and elevation gradient) and terrestrial animals (rivers acting as major dispersal barriers). However, our current knowledge on fine-scale inter-specific and intra-specific genetic diversity of African annual fish is limited to a single, particularly dry region of their distribution (subtropical Mozambique). Using a widespread annual killifish from coastal Tanzania and Kenya, we tested whether the same pattern of genetic divergence pertains to a wet equatorial region in the centre of Nothobranchius distribution.

Results

In populations of Nothobranchius melanospilus species group across its range, we genotyped a part of mitochondrial cytochrome oxidase subunit 1 (COI) gene (83 individuals from 22 populations) and 10 nuclear microsatellite markers (251 individuals from 16 populations). We found five lineages with a clear phylogeographic structure but frequent secondary contact. Mitochondrial lineages were largely congruent with main population genetic clusters identified on microsatellite markers. In the upper Wami basin, populations are isolated as a putative Nothobranchius prognathus, but include also a population from a periphery of the middle Ruvu basin. Other four lineages (including putative Nothobranchius kwalensis) coexisted in secondary contact zones, but possessed clear spatial pattern. Main river channels did not form apparent barriers to dispersal. The most widespread lineage had strong signal of recent population expansion.

Conclusions

We conclude that dispersal of a Nothobranchius species from a wet part of the genus distribution (tropical lowland) is not constrained by main river channels and closely related lineages frequently coexist in secondary contact zones. We also demonstrate contemporary connection between the Ruvu and Rufiji river basins. Our data do not provide genetic support for existence of recently described cryptic species from N. melanospilus complex, but cannot resolve this issue.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

Genetic data reveals a complex history of multiple admixture events in presently allopatric wild gingers (Asarum spp.) showing intertaxonomic clinal variation in calyx lobe length

Clinal variation is a major pattern of observed phenotypic diversity and identifying underlying demographic processes is a necessary step to understand the establishment of clinal variation. The wild ginger series Sakawanum (genus Asarum) comprises four taxa, which exhibit intertaxonomic clinal variation in calyx lobe length across two continental islands isolated by a sea strait. To test alternative hypotheses of the evolutionary history and to determine the implications for the formation of clinal variation, we conducted approximate Bayesian computation (ABC) analysis and ecological niche modeling (ENM). ABC analysis indicated that the scenario assuming multiple admixture events was strongly supported. This scenario assumed two admixture events occurred between morphologically distinct taxa, likely leading to the generation of intermediate taxa. One of the admixture events was estimated to have occurred during the last glacial maximum (LGM), during which the taxa were estimated to have formed a common refugia in southern areas by ENM analysis. Although four taxa are currently distributed allopatrically on different islands and trans-oceanic dispersal appears unlikely, the formation of a land bridge and the geographic range shift to refugia would have allowed secondary contact between previously isolated taxa. This study suggests that clinal variation can be shaped by demographic history including multiple admixtures due to climatic oscillations.

Population genomics of rapid evolution in natural populations: polygenic selection in response to power station thermal effluents

BACKGROUND

Examples of rapid evolution are common in nature but difficult to account for with the standard population genetic model of adaptation. Instead, selection from the standing genetic variation permits rapid adaptation via soft sweeps or polygenic adaptation. Empirical evidence of this process in nature is currently limited but accumulating.

RESULTS

We provide genome-wide analyses of rapid evolution in Fundulus heteroclitus populations subjected to recently elevated temperatures due to coastal power station thermal effluents using 5449 SNPs across two effluent-affected and four reference populations. Bayesian and multivariate analyses of population genomic structure reveal a substantial portion of genetic variation that is most parsimoniously explained by selection at the site of thermal effluents. An FST outlier approach in conjunction with additional conservative requirements identify significant allele frequency differentiation that exceeds neutral expectations among exposed and closely related reference populations. Genomic variation patterns near these candidate loci reveal that individuals living near thermal effluents have rapidly evolved from the standing genetic variation through small allele frequency changes at many loci in a pattern consistent with polygenic selection on the standing genetic variation.

CONCLUSIONS

While the ultimate trajectory of selection in these populations is unknown and we survey only a minority of genomic loci, our findings suggest that polygenic models of adaptation may play important roles in large, natural populations experiencing recent selection due to environmental changes that cause broad physiological impacts.

Invasion genetics of the mummichog (Fundulus heteroclitus): recent anthropogenic introduction in Iberia

Human activities such as trade and transport have increased considerably in the last decades, greatly facilitating the introduction and spread of non-native species at a global level. In the Iberian Peninsula, Fundulus heteroclitus, a small euryhaline coastal fish with short dispersal, was found for the first time in the mid-1970s. Since then, F. heteroclitus has undergone range expansions, colonizing the southern region of Portugal, southwestern coast of Spain and the Ebro Delta in the Mediterranean Sea. Cytochrome b sequences were used to elucidate the species invasion pathway in Iberia. Three Iberian locations (Faro, Cádiz and Ebro Delta) and 13 other locations along the native range of F. heteroclitus in North America were sampled. Results revealed a single haplotype, common to all invasive populations, which can be traced to the northern region of the species’ native range. We posit that the origin of the founder individuals is between New York and Nova Scotia. Additionally, the lack of genetic structure within Iberia is consistent with a recent invasion scenario and a strong founder effect. We suggest the most probable introduction vector is associated with the aquarium trade. We further discuss the hypothesis of a second human-mediated introduction responsible for the establishment of individuals in the Ebro Delta supported by the absence of adequate muddy habitats linking Cádiz and the Ebro Delta. Although the species has a high tolerance to salinity and temperature, ecological niche modelling indicates that benthic habitat constraints prevent along-shore colonisation suggesting that such expansions would need to be aided by human release.

Can secondary contact following range expansion be distinguished from barriers to gene flow?

Secondary contact is the reestablishment of gene flow between sister populations that have diverged. For instance, at the end of the Quaternary glaciations in Europe, secondary contact occurred during the northward expansion of the populations which had found refugia in the southern peninsulas. With the advent of multi-locus markers, secondary contact can be investigated using various molecular signatures including gradients of allele frequency, admixture clines, and local increase of genetic differentiation. We use coalescent simulations to investigate if molecular data provide enough information to distinguish between secondary contact following range expansion and an alternative evolutionary scenario consisting of a barrier to gene flow in an isolation-by-distance model. We find that an excess of linkage disequilibrium and of genetic diversity at the suture zone is a unique signature of secondary contact. We also find that the directionality index ψ, which was proposed to study range expansion, is informative to distinguish between the two hypotheses. However, although evidence for secondary contact is usually conveyed by statistics related to admixture coefficients, we find that they can be confounded by isolation-by-distance. We recommend to account for the spatial repartition of individuals when investigating secondary contact in order to better reflect the complex spatio-temporal evolution of populations and species.

Sea-level driven glacial-age refugia and post-glacial mixing on subtropical coasts, a palaeohabitat and genetic study

Using a novel combination of palaeohabitat modelling and genetic mixture analyses, we identify and assess a sea-level-driven recolonization process following the Last Glacial Maximum (LGM). Our palaeohabitat modelling reveals dramatic changes in estuarine habitat distribution along the coast of California (USA) and Baja California (Mexico). At the LGM (approx. 20 kya), when sea level was approximately 130 m lower, the palaeo-shoreline was too steep for tidal estuarine habitat formation, eliminating this habitat type from regions where it is currently most abundant, and limiting such estuaries to a northern and a southern refugium separated by 1000 km. We assess the recolonization of estuaries formed during post-LGM sea-level rise through examination of refugium-associated alleles and approximate Bayesian computation in three species of estuarine fishes. Results reveal sourcing of modern populations from both refugia, which admix in the newly formed habitat between the refuges. We infer a dramatic peak in habitat area between 15 and 10 kya with subsequent decline. Overall, this approach revealed a previously undocumented dynamic and integrated relationship between sea-level change, coastal processes and population genetics. These results extend glacial refugial dynamics to unglaciated subtropical coasts and have significant implications for biotic response to predicted sea-level rise.

New insight into hybridization and unidirectional introgression between Ammodytes japonicus and Ammodytes heian (Trachiniformes, Ammodytidae)

Based on northern (NOL) and southern (SOL) mitochondrial lineages, recently, it proposed the new species Ammodytes heian and revived the species name Ammodytes japonicus to describe sand lances from the northwestern Pacific Ocean. This study used molecular methods to investigate genetic relationships between the two sand lance species in Korea and Japan. In total, 154 specimens were collected from four locations in Korea (Baengnyeongdo in the Yellow Sea, Tongyeong in the Korean Strait, and Jumunjin and Gijang in the East Sea), and 50 specimens were collected from a single location in Japan (Wakkanai in the Okhotsk Sea). Mitochondrial DNA analysis demonstrated that the individuals from Baengnyeongdo and Tongyeong all belonged to the SOL, whereas those from Gijang, Jumunjin, and Wakkanai included individuals from both the NOL and SOL (over 75% NOL). Population structure analyses were performed on the same individuals using seven microsatellite DNA markers. The population structure analysis based on 201 specimens identified two clusters (named as northern group and southern group), with the admixture proportion (q) of < 0.1 for the northern group in the Backyeongdo and Tongyeong sand lances and < 0.1 for the southern group in the Wakkanai sand lances. The high heterogeneity indicated that the former was probably A. japonicus and the latter probably A. heian. However, the admixture proportion in the Jumunjin and Gijang sand lances was 0.71–0.75 for the southern group, indicating that hybridization and unidirectional introgression from SOL to NOL occurs in southwestern margin of the East Sea. Our findings illustrate the speciation process based on different patterns of gene flow between Korean and Japanese sand lance, which is strongly influenced by both the paleo-climatic change and the contemporary local oceanic current pattern.

Evolved genetic and phenotypic differences due to mitochondrial-nuclear interactions

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only pathway with both nuclear and mitochondrial encoded proteins. The importance of the interactions between these two genomes has recently received more attention because of their potential evolutionary effects and how they may affect human health and disease. In many different organisms, healthy nuclear and mitochondrial genome hybrids between species or among distant populations within a species affect fitness and OxPhos functions. However, what is less understood is whether these interactions impact individuals within a single natural population. The significance of this impact depends on the strength of selection for mito-nuclear interactions. We examined whether mito-nuclear interactions alter allele frequencies for ~11,000 nuclear SNPs within a single, natural Fundulus heteroclitus population containing two divergent mitochondrial haplotypes (mt-haplotypes). Between the two mt-haplotypes, there are significant nuclear allele frequency differences for 349 SNPs with a p-value of 1% (236 with 10% FDR). Unlike the rest of the genome, these 349 outlier SNPs form two groups associated with each mt-haplotype, with a minority of individuals having mixed ancestry. We use this mixed ancestry in combination with mt-haplotype as a polygenic factor to explain a significant fraction of the individual OxPhos variation. These data suggest that mito-nuclear interactions affect cardiac OxPhos function. The 349 outlier SNPs occur in genes involved in regulating metabolic processes but are not directly associated with the 79 nuclear OxPhos proteins. Therefore, we postulate that the evolution of mito-nuclear interactions affects OxPhos function by acting upstream of OxPhos.

Fine-scale genetic structure due to adaptive divergence among microhabitats

It has been suggested that adaptive evolution on ecological timescales shapes communities. However, adaptation among environments relies on isolation or large selection coefficients that exceed migration effects. This reliance is tempered if adaptation is polygenic-does not depend on one allele completely replacing another but instead requires small allele frequency changes at many loci. Thus, whether individuals can evolve adaptation to fine-scale habitat variation (for example, microhabitats) is not resolved. Here we analyze the genetic divergence of the teleost fish, Fundulus heteroclitus, among microhabitats that are <200 m apart in three separate saltmarshes using 4741 single-nucleotide polymorphisms (SNPs). Among these SNPs, 1.3-2.3% have large and highly significant differences among microhabitats (mean F_ST=0.15; false discovery rate ⩽1%). The divergence among microhabitats for these outlier SNPs is larger than that among populations, exceeds neutral expectation and indicates surprising population structure among microhabitats. Thus, we suggest that polygenic selection is surprisingly effective in altering allele frequencies among many different SNPs that share similar biological functions in response to environmental and ecological differences over very small geographic distances. We acknowledge the evolutionary difficulty of large genetic divergence among well-connected habitats. Therefore, these studies are only the first step to discern whether natural selection is responsible and capable of effecting genetic divergence on such a fine scale.

Steep, coincident, and concordant clines in mitochondrial and nuclear-encoded genes in a hybrid zone between subspecies of Atlantic killifish, Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus-specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non-neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non-neutral nuclear clines, SNPs in two nuclear-encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin-binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non-unimodal, suggesting selection against advanced-generation hybrids, possibly due to effects on processes involved in oxidative metabolism.

Secondary contact and local adaptation contribute to genome-wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome-wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome-wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an FST outlier approach. Our results provide novel insight into the well-studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.

Evidence for a bimodal distribution of hybrid indices in a hybrid zone with high admixture

The genetic structure of a hybrid zone can provide insights into the relative roles of the various factors that maintain the zone. Here, we use a multilocus approach to characterize a hybrid zone between two subspecies of killifish (Fundulus heteroclitus, Walbaum 1792) found along the Atlantic coast of North America. We first analysed clinal variation along the Atlantic coast using a single-nucleotide polymorphism in the mitochondrial DNA (mtDNA) displacement loop (D-loop) and a panel of nine nuclear microsatellite markers. A model constraining all clines to the same width and centre was not significantly different from a model in which the clines were allowed to vary independently. Locus-by-locus analysis indicated that the majority of nuclear clines shared the same centre as the mtDNA cline, and the widths of these clines were also narrower than that predicted by a neutral model, suggesting that selection is operating to maintain the hybrid zone. However, two of the nuclear clines had widths greater than the neutral prediction and had centres that were displaced relative to the mtDNA cline centre. We also found that a marsh located near the centre of the mtDNA cline demonstrated a bimodal distribution of nuclear hybrid index values, suggesting a deficit of first-generation hybrids and backcrossed genotypes. Thus, selection against hybrid genotypes may be playing a role in maintaining this hybrid zone and the associated steep nuclear and mtDNA clines.

Comparative pattern of genetic structure in two Mediterranean killifishes Aphanius fasciatus and Aphanius iberus inferred from both mitochondrial and nuclear data

In this study, genetic variation was assessed in Aphanius fasciatus and Aphanius iberus characterized by similar ecological traits but with very different distribution ranges in the Mediterranean area. Five populations of A. iberus and five of A. fasciatus were analysed using five polymorphic microsatellite loci and partial mitochondrial control region (D-loop) sequences. Congruent results were found with both nuclear and mitochondrial molecular markers. The results showed that similar levels of genetic divergence, based on mitochondrial control region sequences, are present among populations of A. iberus and among populations of A. fasciatus despite the very different geographic distance existing among the examined populations of the two species (low geographic distance in A. iberus and high in A. fasciatus). A possible explanation could be that the populations of A. iberus were isolated for a longer time than the populations of A. fasciatus supporting the hypothesis that the split in the lineage leading to A. iberus is older than the split in the lineage leading to A. fasciatus. The possibility that the wide circum-Mediterranean distribution of A. fasciatus ensures the high connectivity of its populations, preventing, in some cases, local differentiation, however, cannot be ruled out.

Pronounced population genetic differentiation in the rock bream Oplegnathus fasciatus inferred from mitochondrial DNA sequences

The population genetic structure of the rock bream (Oplegnathus fasciatus) along the coastal waters of China was estimated based on three mtDNA fragments (D-loop, COI, and Cytb). A total of 112 polymorphic sites were checked, which defined 63 haplotypes. A pattern with high levels of haplotype diversity (hCOI = 0.886 ± 0.034, hCytb = 0.874 ± 0.023) and low levels of nucleotide diversity (лCOI = 0.009 ± 0.005, лCytb = 0.006 ± 0.003) was detected based on the COI and Cytb fragments, and high levels of genetic diversity (hD-loop = 0.995 ± 0.007, лD-loop = 0.021 ± 0.011) were detected from the mtDNA D-loop. The population genetic diversity of O. fasciatus in south China was significantly higher than those of north China. Three genealogical clades were checked in the O. fasciatus populations based on the NJ and MST analyses of mtDNA COI gene sequence, and the genetic distances among the clades ranged from 0.018 to 0.025. Significant population genetic differentiation was also checked based on the Fst (0.331, p = 0.000) and exact p (0.000) test analyses. No significant population differentiations were checked based on mtDNA D-loop and Cytb fragments. Using a variety of phylogenetic methods, coalescent reasoning, and molecular dating interpreted in conjunction with paleoclimatic and physiographic evidences, we inferred that the genetic make-up of extant populations of O. fasciatus was shaped by Pleistocene environmental impacts on the historical demography of this species. Coalescent analyses (neutrality tests, mismatch distribution analysis, and Bayesian skyline analyses) showed that the species along coastline of China has experienced population expansions originated in its most recent history at about 169-175 kya before present.

What is environmental stress? Insights from fish living in a variable environment

Although the term environmental stress is used across multiple fields in biology, the inherent ambiguity associated with its definition has caused confusion when attempting to understand organismal responses to environmental change. Here I provide a brief summary of existing definitions of the term stress, and the related concepts of homeostasis and allostasis, and attempt to unify them to develop a general framework for understanding how organisms respond to environmental stressors. I suggest that viewing stressors as environmental changes that cause reductions in performance or fitness provides the broadest and most useful conception of the phenomenon of stress. I examine this framework in the context of animals that have evolved in highly variable environments, using the Atlantic killifish, Fundulus heteroclitus, as a case study. Consistent with the extreme environmental variation that they experience in their salt marsh habitats, killifish have substantial capacity for both short-term resistance and long-term plasticity in the face of changing temperature, salinity and oxygenation. There is inter-population variation in the sensitivity of killifish to environmental stressors, and in their ability to acclimate, suggesting that local adaptation can shape the stress response even in organisms that are broadly tolerant and highly plastic. Whole-organism differences between populations in stressor sensitivity and phenotypic plasticity are reflected at the biochemical and molecular levels in killifish, emphasizing the integrative nature of the response to environmental stressors. Examination of this empirical example highlights the utility of using an evolutionary perspective on stressors, stress and stress responses.

Understanding population structure and historical demography in a conservation context: population genetics of the endangered Kirengeshoma palmata (Hydrangeaceae)

PREMISE OF THE STUDY

Both historical and contemporary microevolutionary processes greatly influence the genetic patterns of East Asian plant endemics, but the spatial and temporal contexts of these processes remain poorly understood. Here, we investigate the relative influences of historical and contemporary gene flow and drift on the population genetic structure of Kirengeshoma palmata, a perennial herb from East China and South Japan.

METHODS

We used data from nine polymorphic microsatellite loci to assess the levels of genetic diversity, effective population size, and contemporary and historical gene flow for six of the seven known populations.

KEY RESULTS

We found high levels of inbreeding and allelic diversity within populations. Both contemporary and historical migration rates among populations were low, and a test of alternate models of population history strongly favored a model of long-term drift-migration equilibrium. We inferred declines in population size ca. 10,000-100,000 yr ago, but failed to detect recent declines. Bayesian clustering divided K. palmata populations into three genetic clusters, two of which were consistent with a glacial refugium hypothesis for two mountain ranges in East China.

CONCLUSIONS

These results suggest that anthropogenic fragmentation has had little effect on the genetic characteristics of Chinese K. palmata. Rather, past decline in population size due to Late Pleistocene climate change as well as restricted pollen and seed dispersal may have contributed to low levels of both historical and contemporary gene flow, resulting in high genetic differentiation between adjacent mountain ranges due to genetic drift and inbreeding.

Present-day genetic structure of Atlantic salmon (Salmo salar) in Icelandic rivers and ice-cap retreat models

Due to an improved understanding of past climatological conditions, it has now become possible to study the potential concordance between former climatological models and present-day genetic structure. Genetic variability was assessed in 26 samples from different rivers of Atlantic salmon in Iceland (total of 2,352 individuals), using 15 microsatellite loci. F-statistics revealed significant differences between the majority of the populations that were sampled. Bayesian cluster analyses using both prior information and no prior information on sampling location revealed the presence of two distinguishable genetic pools - namely, the Northern (Group 1) and Southern (Group 2) regions of Iceland. Furthermore, the random permutation of different allele sizes among allelic states revealed a significant mutational component to the genetic differentiation at four microsatellite loci (SsaD144, Ssa171, SSsp2201 and SsaF3), and supported the proposition of a historical origin behind the observed variation. The estimated time of divergence, using two different ABC methods, suggested that the observed genetic pattern originated from between the Last Glacial Maximum to the Younger Dryas, which serves as additional evidence of the relative immaturity of Icelandic fish populations, on account of the re-colonisation of this young environment following the Last Glacial Maximum. Additional analyses suggested the presence of several genetic entities which were likely to originate from the original groups detected.

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.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus)

Background

The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.

Results

Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.

Conclusions

The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Significant population genetic structure detected in the small yellow croaker Larimichthys polyactis inferred from mitochondrial control region

The population genetic structure of the small yellow croaker (Larimichthys polyactis) between China and Korea was further estimated by broad-scale sampling locations (Gulf of Bohai, Yellow Sea including Korea). One hundred and seventeen individuals from eight localities from coastal waters of China and Korea were analyzed based on mtDNA control region sequences (5' mtDNA CR). A total of 97 polymorphic sites were checked, which defined 136 haplotypes. A pattern with high levels of haplotype diversity (h=0.994 ± 0.002) and nucleotide diversity (л=0.020 ± 0.010) was detected in the examined range, and the genetic diversity of Korea populations was higher than that of China populations. Population genetic structure analyses (MDS, AMOVA, Fst, Barrier) showed that significant genetic differentiation existed between China and Korea populations. The migration analysis indicated asymmetry migration also existed among populations, which was consistent with the result of population genetic structure. Using a variety of phylogenetic methods, coalescent reasoning, and molecular dating interpreted in conjunction with paleoclimateic and physiographic evidence, we inferred that the genetic make-up of extant populations of L. polyactis was shaped by Pleistocene environmental impacts on the historical demography of this species. Coalescent analyses (Neutrality tests, Mismatch distribution analysis, Bayesian skyline analyses) showed that the species along coastline of China and Korea has experienced population expansions originated in its most recent history at about 32-196 kya and 166-662 kya before present, respectively.

Going where traditional markers have not gone before: utility of and promise for RAD sequencing in marine invertebrate phylogeography and population genomics

Characterization of large numbers of single-nucleotide polymorphisms (SNPs) throughout a genome has the power to refine the understanding of population demographic history and to identify genomic regions under selection in natural populations. To this end, population genomic approaches that harness the power of next-generation sequencing to understand the ecology and evolution of marine invertebrates represent a boon to test long-standing questions in marine biology and conservation. We employed restriction-site-associated DNA sequencing (RAD-seq) to identify SNPs in natural populations of the sea anemone Nematostella vectensis, an emerging cnidarian model with a broad geographic range in estuarine habitats in North and South America, and portions of England. We identified hundreds of SNP-containing tags in thousands of RAD loci from 30 barcoded individuals inhabiting four locations from Nova Scotia to South Carolina. Population genomic analyses using high-confidence SNPs resulted in a highly-resolved phylogeography, a result not achieved in previous studies using traditional markers. Plots of locus-specific FST against heterozygosity suggest that a majority of polymorphic sites are neutral, with a smaller proportion suggesting evidence for balancing selection. Loci inferred to be under balancing selection were mapped to the genome, where 90% were located in gene bodies, indicating potential targets of selection. The results from analyses with and without a reference genome supported similar conclusions, further highlighting RAD-seq as a method that can be efficiently applied to species lacking existing genomic resources. We discuss the utility of RAD-seq approaches in burgeoning Nematostella research as well as in other cnidarian species, particularly corals and jellyfishes, to determine phylogeographic relationships of populations and identify regions of the genome undergoing selection.

Continental phylogeography of an ecologically and morphologically diverse Neotropical songbird, Zonotrichia capensis

BACKGROUND

The Neotropics are exceptionally diverse, containing roughly one third of all extant bird species on Earth. This remarkable species richness is thought to be a consequence of processes associated with both Andean orogenesis throughout the Tertiary, and climatic fluctuations during the Quaternary. Phylogeographic studies allow insights into how such events might have influenced evolutionary trajectories of species and ultimately contribute to a better understanding of speciation. Studies on continentally distributed species are of particular interest because different populations of such taxa may show genetic signatures of events that impacted the continent-wide biota. Here we evaluate the genealogical history of one of the world's most broadly-distributed and polytypic passerines, the rufous-collared sparrow (Zonotrichia capensis).

RESULTS

We obtained control region DNA sequences from 92 Zonotrichia capensis individuals sampled across the species' range (Central and South America). Six additional molecular markers, both nuclear and mitochondrial, were sequenced for a subset of individuals with divergent control region haplotypes. Median-joining network analysis, and Bayesian and maximum parsimony phylogenetic analyses all recovered three lineages: one spanning Middle America, the Dominican Republic, and north-western South America; one encompassing the Dominican Republic, Roraima (Venezuela) and La Paz (Bolivia) south to Tierra del Fuego, Argentina; and a third, including eastern Argentina and Brazil. Phylogenetic analyses suggest that the Middle American/north-western South American clade is sister to the remaining two. Bayesian and maximum likelihood coalescent simulations used to study lineage demographic history, diversification times, migration rates and population expansion together suggested that diversification of the three lineages occurred rapidly during the Pleistocene, with negligible gene flow, leaving genetic signatures of population expansions.

CONCLUSIONS

The Pleistocene history of the rufous-collared sparrow involved extensive range expansion from a probable Central American origin. Its remarkable morphological and behavioral diversity probably represents recent responses to local conditions overlying deeper patterns of lineage diversity, which are themselves produced by isolation and the history of colonization of South America.

Interindividual variation in Complex I activity in Fundulus heteroclitus along a steep thermocline

The first enzyme in the oxidative phosphorylation pathway is Complex I (EC 1.6.5.3). Complex I is a large heteromeric enzyme complex with 45 protein subunits that translocates H(+) ions across the mitochondrial inner membrane. Among northern and southern populations of the teleost fish Fundulus heteroclitus, Complex I subunits have fixed amino acid substitutions. Additionally, there are differences in oxidative phosphorylation activity among populations of F. heteroclitus. To investigate whether these differences are related to Complex I, enzyme activity was measured in 121 individuals from five populations of F. heteroclitus and its sister species Fundulus grandis acclimated to a constant 20°C temperature. Within each population, Complex I activity is highly variable among individuals of F. heteroclitus (coefficient of variation percentage among individuals has a mean of 90% in the five F. heteroclitus populations), and the mean Complex I activity among populations is significantly different at the latitudinal extremes of the range. Importantly, Complex I activity is more similar between F. heteroclitus from the southernmost population and its sister species F. grandis than to the northern populations of F. heteroclitus, suggesting important evolutionary differences. Unexpectedly, the activity is nearly fourfold higher in southern populations than northern populations. Mitochondrial density appears to compensate partially for decreased activity in northern individuals; activity per wet weight is only twofold higher in southern populations. We suggest that some of the variation in Complex I activity is genetically based and thus is being influenced by directional selection. However, this conclusion presents a conundrum: there should not be so much variation in Complex I activity within a population if this variation is biologically important.

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.

Microsatellite and mtDNA analysis of lake trout, Salvelinus namaycush, from Great Bear Lake, Northwest Territories: impacts of historical and contemporary evolutionary forces on Arctic ecosystems

Resolving the genetic population structure of species inhabiting pristine, high latitude ecosystems can provide novel insights into the post-glacial, evolutionary processes shaping the distribution of contemporary genetic variation. In this study, we assayed genetic variation in lake trout (Salvelinus namaycush) from Great Bear Lake (GBL), NT and one population outside of this lake (Sandy Lake, NT) at 11 microsatellite loci and the mtDNA control region (d-loop). Overall, population subdivision was low, but significant (global F(ST) θ = 0.025), and pairwise comparisons indicated that significance was heavily influenced by comparisons between GBL localities and Sandy Lake. Our data indicate that there is no obvious genetic structure among the various basins within GBL (global F(ST) = 0.002) despite the large geographic distances between sampling areas. We found evidence of low levels of contemporary gene flow among arms within GBL, but not between Sandy Lake and GBL. Coalescent analyses suggested that some historical gene flow occurred among arms within GBL and between GBL and Sandy Lake. It appears, therefore, that contemporary (ongoing dispersal and gene flow) and historical (historical gene flow and large founding and present-day effective population sizes) factors contribute to the lack of neutral genetic structure in GBL. Overall, our results illustrate the importance of history (e.g., post-glacial colonization) and contemporary dispersal ecology in shaping genetic population structure of Arctic faunas and provide a better understanding of the evolutionary ecology of long-lived salmonids in pristine, interconnected habitats.

Evolutionary and functional analyses of cytochrome P4501A promoter polymorphisms in natural populations

The functional importance of variable, transcriptional regulatory sequences within and among natural populations is largely unexplored. We analysed the cytochrome P4501A (CYP1A) promoter in three populations of the minnow, Fundulus heteroclitus, because two SNPs in the promoter and first intron of CYP1A are under selection in populations adapted to pollutants. To define the importance of these SNPs, 1630 bp of the CYP1A promoter and first intron and exon were sequenced in eight individuals from three populations: a population from a polluted environment resistant to some aromatic pollutants and two flanking reference populations. CYP1A is induced by many aromatic pollutants, but in populations adapted to pollutants, CYP1A has been shown to be refractory to induction. We were interested in understanding whether variation in the CYP1A promoter explains mechanism(s) of adaptation to these aromatic pollutants. The CYP1A promoter was extremely variable (an average of 9.3% of the promoter nucleotides varied among all populations) and exhibited no fixed differences between populations. As CYP1A is poorly inducible in adapted fish, we hypothesized that CYP1A promoter regions might vary functionally between populations. Unexpectedly, in vitro analysis showed significantly greater transcription from CYP1A promoters found in the population from the polluted environment relative to promoters found in both reference populations. Thus, despite extensive variation among populations and lack of fixed differences between populations, individuals from a polluted environment have significantly enhanced promoter activity. These data demonstrate that intraspecific variation, which provides the raw material for natural selection to act on, can occur while maintaining promoter function.

The effects of mitochondrial genotype on hypoxic survival and gene expression in a hybrid population of the killifish, Fundulus heteroclitus

The physiological link between oxygen availability and mitochondrial function is well established. However, whether or not fitness variation is associated with mitochondrial genotypes in the field remains a contested topic in evolutionary biology. In this study, we draw on a population of the teleost fish, Fundulus heteroclitus, where functionally distinct subspecies hybridize, likely as a result of past glacial events. We had two specific aims: (i) to determine the effect of mtDNA genotype on survivorship of male and female fish under hypoxic stress and (ii) to determine the effect of hypoxic stress, sex and mtDNA genotype on gene expression. We found an unexpected and highly significant effect of sex on survivorship under hypoxic conditions, but no significant effect of mtDNA genotype. Gene expression analyses revealed hundreds of transcripts differentially regulated by sex and hypoxia. Mitochondrial transcripts and other predicted pathways were among those influenced by hypoxic stress, and a transcript corresponding to the mtDNA control region was the most highly suppressed transcript under the conditions of hypoxia. An RT-PCR experiment on the control region was consistent with microarray results. Effects of mtDNA sequence variation on genome expression were limited; however, a potentially important epistasis between mtDNA sequence and expression of a nuclear-encoded mitochondrial translation protein was discovered. Overall, these results confirm that mitochondrial regulation is a major component of hypoxia tolerance and further suggest that purifying selection has been the predominant selective force on mitochondrial genomes in these two subspecies.

Common mechanism underlies repeated evolution of extreme pollution tolerance

Human alterations to the environment can exert strong evolutionary pressures, yet contemporary adaptation to human-mediated stressors is rarely documented in wildlife populations. A common-garden experimental design was coupled with comparative transcriptomics to discover evolved mechanisms enabling three populations of killifish resident in urban estuaries to survive normally lethal pollution exposure during development, and to test whether mechanisms are unique or common across populations. We show that killifish populations from these polluted sites have independently converged on a common adaptive mechanism, despite variation in contaminant profiles among sites. These populations are united by a similarly profound desensitization of aryl-hydrocarbon receptor-mediated transcriptional activation, which is associated with extreme tolerance to the lethal effects of toxic dioxin-like pollutants. The rapid, repeated, heritable and convergent nature of evolved tolerance suggests that ancestral killifish populations harboured genotypes that enabled adaptation to twentieth-century industrial pollutants.

Repeated sampling detects gene flow in a flightless ground beetle in a fragmented landscape

Secondary clines level down in the course of time if the gene flow is not interrupted. Temporally repeated sampling of populations in a cline allows the investigation not only of its occurrence but also of the estimation of the amount of ongoing gene flow. We reinvestigated an allozyme gradient in Carabus auronitens populations in the Westphalian Lowlands (northwestern Germany) 15 to 20 years after it had originally been recorded. A total of 977 individuals of this flightless woodland species from 29 sample sites were genotyped at the diallelic Est-1 locus in 2005-2006 and compared to former findings, collected in 1985-1994 from the same populations. Both data sets showed clinal variation. Pairwise differences between the samples of both data sets indicated significant decrease in the steepness of the cline during the past 15 to 20 years. The estimated average gene flow per generation is 0.6% of each beetle population. Ongoing gene flow in the flightless ground beetle C. auronitens led to a less pronounced cline despite a stable degree of fragmentation (and connectivity) of the landscape. Migration and gene flow were obviously enabled by the numerous hedgerows. The corridors are seen to be a prerequisite for migration between populations and for possible future range shifts of forest insect species.

Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.

Niche dimensions in fishes: an integrative view

Current shifts in ecosystem composition and function emphasize the need for an understanding of the links between environmental factors and organism fitness and tolerance. The examples discussed here illustrate how recent progress in the field of comparative physiology may provide a better mechanistic understanding of the ecological concepts of the fundamental and realized niches and thus provide insights into the impacts of anthropogenic disturbance. Here we argue that, as a link between physiological and ecological indicators of organismal performance, the mechanisms shaping aerobic scope and passive tolerance set the dimensions of an animal's niche, here defined as its capacity to survive, grow, behave, and interact with other species. We demonstrate how comparative studies of cod or killifish populations in a latitudinal cline have unraveled mitochondrial mechanisms involved in establishing a species' niche, performance, and energy budget. Riverine fish exemplify how the performance windows of various developmental stages follow the dynamic regimes of both seasonal temperatures and river hydrodynamics, as synergistic challenges. Finally, studies of species in extreme environments, such as the tilapia of Lake Magadi, illustrate how on evolutionary timescales functional and morphological shifts can occur, associated with new specializations. We conclude that research on the processes and time course of adaptations suitable to overcome current niche limits is urgently needed to assess the resilience of species and ecosystems to human impact, including the challenges of global climate change.

Regional differentiation and post-glacial expansion of the Atlantic silverside, Menidia menidia, an annual fish with high dispersal potential

The coastal marine environment of the Northwest Atlantic contains strong environmental gradients that create distinct marine biogeographic provinces by limiting dispersal, recruitment, and survival. This region has also been subjected to numerous Pleistocene glacial cycles, resulting in repeated extirpations and recolonizations in northern populations of marine organisms. In this study, we examined patterns of genetic structure and historical demography in the Atlantic silverside, Menidia menidia, an annual marine fish with high dispersal potential but with well-documented patterns of clinal phenotypic adaptation along the environmental gradients of the Northwest Atlantic. Contrary to previous studies indicating genetic homogeneity that should preclude regional adaptation, results demonstrate subtle but significant (FST = 0.07; P < 0.0001) genetic structure among three phylogeographic regions that partially correspond with biogeographic provinces, suggesting regional limits to gene flow. Tests for non-equilibrium population dynamics and latitudinal patterns in genetic diversity indicate northward population expansion from a single southern refugium following the last glacial maximum, suggesting that phylogeographic and phenotypic patterns have relatively recent origins. The recovery of phylogeographic structure and the partial correspondence of these regions to recognized biogeographic provinces suggest that the environmental gradients that shape biogeographic patterns in the Northwest Atlantic may also limit gene flow in M. menidia, creating phylogeographic structure and contributing to the creation of latitudinal phenotypic clines in this species.

Embryonic gene expression among pollutant resistant and sensitive Fundulus heteroclitus populations

Changes in gene expression, coupled with biochemical, physiological, and behavioral alterations, play a critical role in adaptation to environmental stress. Our goal was to explore ways natural populations may have adapted to local, polluted environments. We took advantage of natural populations of Fundulus heteroclitus, one of the few studied fish species in North America that has established resistant populations in highly contaminated urban estuaries. We analyzed morphology, physiology, and gene expression of developing F. heteroclitus embryos during late organogenesis (stage 31); these embryos were from both resistant and sensitive populations and were raised in a common, unpolluted environment. While cardiac heart rates show significant differences between embryos of parents from clean and heavily contaminated Superfund sites, time-to-stage, embryo morphology, and gene expression profile analyses do not differ significantly between untreated embryos from resistant and sensitive populations. Further evaluation that includes tissue-specific approaches in gene expression analysis and larger sample sizes may be necessary to highlight important phenotypes associated with mechanisms of sensitivity and resistance among natural F. heteroclitus embryo populations. Alternatively, population differences may be masked by developmental canalization, and biologically important differences between sensitive and resistant embryos may only manifest with exposure (e.g., be dependent on gene by environment interactions).