Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.)

Polymorphism, selection and tandem duplication of transferrin genes in Atlantic cod (Gadus morhua)--conserved synteny between fish monolobal and tetrapod bilobal transferrin loci

BACKGROUND

The two homologous iron-binding lobes of transferrins are thought to have evolved by gene duplication of an ancestral monolobal form, but any conserved synteny between bilobal and monolobal transferrin loci remains unexplored. The important role played by transferrin in the resistance to invading pathogens makes this polymorphic gene a highly valuable candidate for studying adaptive divergence among local populations.

RESULTS

The Atlantic cod genome was shown to harbour two tandem duplicated serum transferrin genes (Tf1, Tf2), a melanotransferrin gene (MTf), and a monolobal transferrin gene (Omp). Whereas Tf1 and Tf2 were differentially expressed in liver and brain, the Omp transcript was restricted to the otoliths. Fish, chicken and mammals showed highly conserved syntenic regions in which monolobal and bilobal transferrins reside, but contrasting with tetrapods, the fish transferrin genes are positioned on three different linkage groups. Sequence alignment of cod Tf1 cDNAs from Northeast (NE) and Northwest (NW) Atlantic populations revealed 22 single nucleotide polymorphisms (SNP) causing the replacement of 16 amino acids, including eight surface residues revealed by the modelled 3D-structures, that might influence the binding of pathogens for removal of iron. SNP analysis of a total of 375 individuals from 14 trans-Atlantic populations showed that the Tf1-NE variant was almost fixed in the Baltic cod and predominated in the other NE Atlantic populations, whereas the NW Atlantic populations were more heterozygous and showed high frequencies of the Tf-NW SNP alleles.

CONCLUSIONS

The highly conserved synteny between fish and tetrapod transferrin loci infers that the fusion of tandem duplicated Omp-like genes gave rise to the modern transferrins. The multiple nonsynonymous substitutions in cod Tf1 with putative structural effects, together with highly divergent allele frequencies among different cod populations, strongly suggest evidence for positive selection and local adaptation in trans-Atlantic cod populations.

Strong gene flow and lack of stable population structure in the face of rapid adaptation to local temperature in a spring-spawning salmonid, the European grayling (Thymallus thymallus)

Gene flow has the potential to both constrain and facilitate adaptation to local environmental conditions. The early stages of population divergence can be unstable because of fluctuating levels of gene flow. Investigating temporal variation in gene flow during the initial stages of population divergence can therefore provide insights to the role of gene flow in adaptive evolution. Since the recent colonization of Lake Lesjaskogsvatnet in Norway by European grayling (Thymallus thymallus), local populations have been established in over 20 tributaries. Multiple founder events appear to have resulted in reduced neutral variation. Nevertheless, there is evidence for local adaptation in early life-history traits to different temperature regimes. In this study, microsatellite data from almost a decade of sampling were assessed to infer population structuring and its temporal stability. Several alternative analyses indicated that spatial variation explained 2-3 times more of the divergence in the system than temporal variation. Over all samples and years, there was a significant correlation between genetic and geographic distance. However, decomposed pairwise regression analysis revealed differing patterns of genetic structure among local populations and indicated that migration outweighs genetic drift in the majority of populations. In addition, isolation by distance was observable in only three of the six years, and signals of population bottlenecks were observed in the majority of samples. Combined, the results suggest that habitat-specific adaptation in this system has preceded the development of consistent population substructuring in the face of high levels of gene flow from divergent environments.

The future of Baltic Sea populations: local extinction or evolutionary rescue?

Environmental change challenges local and global survival of populations and species. In a species-poor environment like the Baltic Sea this is particularly critical as major ecosystem functions may be upheld by single species. A complex interplay between demographic and genetic characteristics of species and populations determines risks of local extinction, chances of re-establishment of lost populations, and tolerance to environmental changes by evolution of new adaptations. Recent studies show that Baltic populations of dominant marine species are locally adapted, have lost genetic variation and are relatively isolated. In addition, some have evolved unusually high degrees of clonality and others are representatives of endemic (unique) evolutionary lineages. We here suggest that a consequence of local adaptation, isolation and genetic endemism is an increased risk of failure in restoring extinct Baltic populations. Additionally, restricted availability of genetic variation owing to lost variation and isolation may negatively impact the potential for evolutionary rescue following environmental change.

Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.

Genetic polymorphism and its potential relation to environmental stress in five populations of the European flounder Platichthys flesus, along the French Atlantic coast

In this study, new DNA markers were explored for the flounder Platichthys flesus. cDNA and genomic sequences of the genes encoding the glyceraldehyde-3-phosphate-deshydrogenase (GAPDH), the cytosolic creatine kinase (CK), the prostaglandin D synthase (PGDS) and the betaine homocysteine methyltransferase (BHMT) were characterized. The tumour suppressor p53 gene structure was already described. A PCR-SSCP (Single Strand Conformation Polymorphism) analysis was finally conducted to study the genetic polymorphism of different populations of flounders collected along the French Atlantic coast. Four highly contaminated French estuaries (Seine, Vilaine, Loire and Gironde) were sampled and compared to a reference estuary (Ster) to explore possible selective effect of the environment on specific allelic frequencies. Our results showed that two loci p53 and PGDS, could be potential markers of chemical stress: p53A allele frequency increased in contaminated systems compared to the reference system. In the Vilaine estuary, PGDS polymorphism could be related to pesticide stress.

Differential expression of the heat shock protein Hsp70 in natural populations of the tilapia, Sarotherodon melanotheron, acclimatised to a range of environmental salinities

Background

The relationship between environmental variation and induction of heat shock proteins (Hsps) has been much documented under experimental conditions. However, very little is known about such induction in natural populations acclimatised to prevailing environmental conditions. Furthermore, while induction of stress proteins has been well documented in response to environmental contaminants and thermal stressors, little is known about whether factors, such as extreme salinity, are also potential inductors. The black-chinned tilapia Sarotherodon melanotheron is unusual for its ability to colonise estuarine environments in West Africa that are characterised by extremely high salinities. The relationships between mRNA levels of the 70 kDa heat shock protein (Hsp70) and Na⁺, K⁺-ATPase1α (Naka) in the gills, environmental salinity, and a life-history trait (condition factor) were investigated in wild populations of this species sampled from three locations in the Saloum estuary, at salinities ranging from 40 to 100 psu.

Results

The highest Hsp70 and Naka mRNA levels, and the poorest condition factors were recorded in the most saline sampling site (100 psu). The Hsp70 and Naka mRNA were correlated amongst themselves and showed a direct positive correlation with environmental salinity, and a negative correlation with fish condition factor. Thus, the Hsp70 is constitutively overexpressed by S. melanotheron acclimatised to extreme hypersalinity.

Conclusions

These results indicate that, although S. melanotheron can colonise extremely saline environments, the overexpression of Hsp70 combined with the higher Naka mRNA expression reveals that this represents a chronic stress. The induction of Hsp70 was, therefore, a biomarker of chronic hyper-osmotic stress which presumably can be linked to the impaired growth performance and precocious reproduction that have been demonstrated in the populations at the extremely saline sites.

Local selection and population structure in a deep-sea fish, the roundnose grenadier (Coryphaenoides rupestris)

Local populations within a species can become isolated by stochastic or adaptive processes, though it is most commonly the former that we quantify. Using presumably neutral markers we can assess the time-dependent process of genetic drift, and thereby quantify patterns of differentiation in support of the effective management of diversity. However, adaptive differences can be overlooked in these studies, and these are the very characteristics that we hope to conserve by managing neutral diversity. In this study, we used 16 hypothetically neutral microsatellite markers to investigate the genetic structure of the roundnose grenadier in the North Atlantic. We found that one locus was a clear outlier under directional selection, with F(ST) values much greater than at the remaining loci. Differentiation between populations at this locus was related to depth, suggesting directional selection, presumably acting on a linked locus. Considering only the loci identified as neutral, there remained significant population structure over the region of the North Atlantic studied. In addition to a weak pattern of isolation by distance, we identified a putative barrier to gene flow between sample sites either side of the Charlie-Gibbs Fracture Zone, which marks the location where the sub-polar front crosses the Mid-Atlantic Ridge. This may reflect a boundary across which larvae are differentially distributed in separate current systems to some extent, promoting differentiation by drift. Structure due to both drift and apparent selection should be considered in management policy.

Allopatric divergence and speciation in coral reef fish: the three-spot dascyllus, Dascyllus trimaculatus, species complex

Long pelagic larval phases and the absence of physical barriers impede rapid speciation and contrast the high diversity observed in marine ecosystems such as coral reefs. In this study, we used the three-spot dascyllus (Dascyllus trimaculatus) species complex to evaluate speciation modes at the spatial scale of the Indo-Pacific. The complex includes four recognized species and four main color morphs that differ in distribution. Previous studies of the group using mitochondrial DNA revealed a noncongruence between color morphs and genetic groupings; with two of the color morphs grouped together and one color morph separated into three clades. Using extensive geographic sampling of 563 individuals and a combination of mitochondrial DNA sequences and 13 nuclear microsatellites, we defined population/species boundaries and inferred different speciation modes. The complex is composed of seven genetically distinct entities, some of which are distinct morphologically. Despite extensive dispersal abilities and an apparent lack of barriers, observed genetic partitions are consistent with allopatric speciation. However, ecological pressure, assortative mating, and sexual selection, were likely important during periods of geographical isolation. This study therefore suggests that primarily historical factors later followed by ecological factors caused divergence and speciation in this group of coral reef fish.

Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua)

BACKGROUND

Marine fishes have been shown to display low levels of genetic structuring and associated high levels of gene flow, suggesting shallow evolutionary trajectories and, possibly, limited or lacking adaptive divergence among local populations. We investigated variation in 98 gene-associated single nucleotide polymorphisms (SNPs) for evidence of selection in local populations of Atlantic cod (Gadus morhua L.) across the species distribution.

RESULTS

Our global genome scan analysis identified eight outlier gene loci with very high statistical support, likely to be subject to directional selection in local demes, or closely linked to loci under selection. Likewise, on a regional south/north transect of central and eastern Atlantic populations, seven loci displayed strongly elevated levels of genetic differentiation. Selection patterns among populations appeared to be relatively widespread and complex, i.e. outlier loci were generally not only associated with one of a few divergent local populations. Even on a limited geographical scale between the proximate North Sea and Baltic Sea populations four loci displayed evidence of adaptive evolution. Temporal genome scan analysis applied to DNA from archived otoliths from a Faeroese population demonstrated stability of the intra-population variation over 24 years. An exploratory landscape genetic analysis was used to elucidate potential effects of the most likely environmental factors responsible for the signatures of local adaptation. We found that genetic variation at several of the outlier loci was better correlated with temperature and/or salinity conditions at spawning grounds at spawning time than with geographic distance per se.

CONCLUSION

These findings illustrate that adaptive population divergence may indeed be prevalent despite seemingly high levels of gene flow, as found in most marine fishes. Thus, results have important implications for our understanding of the interplay of evolutionary forces in general, and for the conservation of marine biodiversity under rapidly increasing evolutionary pressure from climate and fisheries induced changes in local environments.

Genetic composition of Atlantic and Mediterranean recruits of European eel Anguilla anguilla based on EST-linked microsatellite loci

Anguilla anguilla glass eels arriving at two Mediterranean and two Atlantic sites were tested for differences in genetic composition between regions using a total of 23 microsatellite loci developed from an expressed sequence tag (EST) library. Hierarchical analysis of molecular variance indicated a non-significant difference between regions (Mediterranean v. Atlantic), which contrasted with the significant differences observed between samples within regions. The existence of a single spawning site for all A. anguilla individuals and extensive migration loop with great opportunity for mixing of individuals might explain the homogeneity in genetic composition found between regions. The observation of a (small-scale) pattern of genetic patchiness among intra-annual samples (arrival waves) within geographic regions does not conflict with the lack of (large-scale) geographic sub-structuring found between the Mediterranean and Atlantic regions, but most likely is a consequence of the strong dependence of A. anguilla on oceanic conditions in the Sargasso Sea that might result in a limited parental contribution to each spawning event. The comparison of Atlantic and Mediterranean A. anguilla glass eel recruits based on EST-linked microsatellite loci provides evidence supporting the hypothesis of panmixia A. anguilla across Europe.

Spatial and temporal population genetic structure of four northeastern Pacific littorinid gastropods: the effect of mode of larval development on variation at one mitochondrial and two nuclear DNA markers

We investigated the effect of development mode on the spatial and temporal population genetic structure of four littorinid gastropod species. Snails were collected from the same three sites on the west coast of Vancouver Island, Canada in 1997 and again in 2007. DNA sequences were obtained for one mitochondrial gene, cytochrome b (Cyt b), and for up to two nuclear genes, heat shock cognate 70 (HSC70) and aminopeptidase N intron (APN54). We found that the mean level of genetic diversity and long-term effective population sizes (N(e)) were significantly greater for two species, Littorina scutulata and L. plena, that had a planktotrophic larval stage than for two species, Littorina sitkana and L. subrotundata, that laid benthic egg masses which hatched directly into crawl-away juveniles. Predictably, two poorly dispersing species, L. sitkana and L. subrotundata, showed significant spatial genetic structure at an 11- to 65-km geographical scale that was not observed in the two planktotrophic species. Conversely, the two planktotrophic species had more temporal genetic structure over a 10-year interval than did the two direct-developing species and showed highly significant temporal structure for spatially pooled samples. The greater temporal genetic variation of the two planktotrophic species may have been caused by their high fecundity, high larval dispersal, and low but spatially correlated early survivorship. The sweepstakes-like reproductive success of the planktotrophic species could allow a few related females to populate hundreds of kilometres of coastline and may explain their substantially larger temporal genetic variance but lower spatial genetic variance relative to the direct-developing species.

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.

Intraspecific variation in expression of candidate genes for osmoregulation, heme biosynthesis and stress resistance suggests local adaptation in European flounder (Platichthys flesus)

Despite the recent discovery of significant genetic structuring in a large number of marine organisms, the evolutionary significance of these often minute genetic differences are still poorly understood. To elucidate the adaptive relevance of low genetic differentiation among marine fish populations, we studied expression differences of osmoregulatory and stress genes in genetically weakly differentiated populations of the European flounder (Platichthys flesus), distributed across a natural salinity gradient. Flounders were maintained in a long-term reciprocal transplantation experiment mimicking natural salinities in the North Sea and the Baltic Sea. Applying real-time quantitative PCR and microarray analysis we studied expression of four candidate genes (hsp70, angiotensinogen, Na/K-ATPase-alpha and 5-aminolevulinic acid synthase (ALAS)) in gill, kidney and liver tissues. Genes involved in osmoregulative processes (Na/K-ATPases-alpha and angiotensinogen) showed highly plastic but similar expression in the two populations dependent on environmental salinity. However, we observed a unique sixfold up-regulation of hsp70 in kidney tissue of flounder from the North Sea following long-term acclimation to Baltic salinities. Similarly, significant differences between North Sea and Baltic flounders in expression of ALAS in response to different salinities were found in gill and liver tissue. These findings strongly suggest that gene expression in flounders is shaped by adaptation to local environmental conditions. This identification of adaptive differences in high gene flow marine organisms adds a new dimension to our current understanding of evolutionary processes in the sea and is of paramount importance for identification, protection and sustainable management of marine biodiversity.

Landscape genetics goes to sea

A recent study revealing geographical and environmental barriers to gene flow in the harbour porpoise shows the great potential of 'landscape genetics' when applied to marine organisms.