Peter Foged Larsen

Population genomics of marine fishes: identifying adaptive variation in space and time.

Studies of adaptive evolution have experienced a recent revival in population genetics of natural populations and there is currently much focus on identifying genomic signatures of selection in space and time. Insights into local adaptation, adaptive response to global change and evolutionary consequences of selective harvesting can be generated through population genomics studies, allowing the separation of the effects invoked by neutral processes (drift‐migration) from those due to selection. Such knowledge is important not only for improving our basic understanding of natural as well as human‐induced evolutionary processes, but also for predicting future trajectories of biodiversity and for setting conservation priorities. Marine fishes possess a number of features rendering them well suited for providing general insights into adaptive genomic evolution in natural populations. These include well‐described population structures, substantial and rapidly developing genomic resources and abundant archived samples enabling temporal studies. Furthermore, superior possibilities for conducting large‐scale experiments under controlled conditions, due to the economic resources provided by the large and growing aquaculture industry, hold great promise for utilizing recent technological developments. Here, we review achievements in marine fish genomics to date and highlight potential avenues for future research, which will provide both general insights into evolution in high gene flow species, as well as specific knowledge which can lead to improved management of marine organisms.

Adaptive differences in gene expression in European flounder (Platichthys flesus)

Population structure was previously believed to be very limited or absent in classical marine fishes, but recently, evidence of weakly differentiated local populations has been accumulating using noncoding microsatellite markers. However, the evolutionary significance of such minute genetic differences remains unknown. Therefore, in order to elucidate the relationship between genetic markers and adaptive divergence among populations of marine fishes, we combined cDNA microarray and microsatellite analysis in European flounders (Platichthys flesus). We demonstrate that despite extremely low levels of neutral genetic divergence, a high number of genes were significantly differentially expressed between North Sea and Baltic Sea flounders maintained in a long‐term reciprocal transplantation experiment mimicking natural salinities. Several of the differentially regulated genes could be directly linked to fitness traits. These findings demonstrate that flounders, despite little neutral genetic divergence between populations, are differently adapted to local environmental conditions and imply that adaptation in gene expression could be common in other marine organisms with similar low levels of population subdivision.

Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered european eel

BackgroundOnce highly abundant, the European eel (Anguilla anguilla L.; Anguillidae; Teleostei) is considered to be critically endangered and on the verge of extinction, as the stock has declined by 90-99% since the 1980s. Yet, the species is poorly characterized at molecular level with little sequence information available in public databases.ResultsThe first European eel transcriptome was obtained by 454 FLX Titanium sequencing of a normalized cDNA library, produced from a pool of 18 glass eels (juveniles) from the French Atlantic coast and two sites in the Mediterranean coast. Over 310,000 reads were assembled in a total of 19,631 transcribed contigs, with an average length of 531 nucleotides. Overall 36% of the contigs were annotated to known protein/nucleotide sequences and 35 putative miRNA identified.ConclusionsThis study represents the first transcriptome analysis for a critically endangered species. EeelBase, a dedicated database of annotated transcriptome sequences of the European eel is freely available at http://compgen.bio.unipd.it/eeelbase. Considering the multiple factors potentially involved in the decline of the European eel, including anthropogenic factors such as pollution and human-introduced diseases, our results will provide a rich source of data to discover and identify new genes, characterize gene expression, as well as for identification of genetic markers scattered across the genome to be used in various applications.

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-α and 5-aminolevulinic acid synthase (ALAS)) in gill, kidney and liver tissues. Genes involved in osmoregulative processes (Na/K-ATPases-α 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.

A resource of genome-wide single-nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel

Reduced representation genome sequencing such as restriction‐site‐associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single‐nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome‐wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome‐wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19 703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long‐term effective population size was estimated to range between 132 000 and 1 320 000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82 425 loci and 376 918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.

The Effects of Vinclozolin, an Anti-Androgenic Fungicide, on Male Guppy Secondary Sex Characters and Reproductive Success

Abstract Despite the enormous volume of research concerning the various effects of chemicals with endocrine-disrupting properties in fish, there is still very little evidence that endocrine disruption can adversely affect individual fertility and, hence, pose problems for the population. In the present study, guppies (Poecilia reticulata) were fed with the anti-androgenic fungicide vinclozolin at concentrations ranging from 1.8 to 180 mg/kg from 8–14 wk of age. Male sperm count and the intensity of his sexual display behavior were significantly reduced in treatment groups, which was in line with the results of previous studies. Here, we show further that these impairments translate into reduced fertility, measured as the size of the females first clutch. Also, this reduced fertility was correlated to the male sperm count, but not to the intensity of the male sexual display. Finally, by crossing exposed with unexposed animals, we show that the adverse effect of vinclozolin on reproduction is mediated through the male alone.

Stocking impact and temporal stability of genetic composition in a brackish northern pike population ( Esox lucius L.), assessed using microsatellite DNA analysis of historical and contemporary samples

During the last decade, brackish northern pike populations in Denmark have been subject to stocking programmes, using nonindigenous pike from freshwater lakes, in order to compensate for drastic population declines. The present study was designed to investigate the genetic impact of stocking freshwater pike into a brackish pike population in Stege Nor, Denmark. We analysed polymorphism at eight microsatellite loci in samples representing the indigenous Stege Nor population prior to stocking (ie from 1956 to 1957), along with a sample of the contemporary Stege Nor population and samples from the three populations used for stocking. Despite large numbers of stocked fry, the results from both individual and population level admixture analyses demonstrated extremely poor performance and <1% introgression of stocked freshwater pike into the brackish pike population. Furthermore, pairwise FST estimates between samples demonstrated close genetic relationship among temporal samples from Stege Nor, indicating temporal stability over the last 45 years. We also estimated the effective population size (Ne) of pike in Stege Nor and applied a test for recent population bottlenecks. The harmonic mean of Ne was relatively high (>250), but there were indications of bottlenecks in all samples and populations. We ascribe this finding to historical rather than recent bottlenecks, possibly dating back to founder events associated with postglacial recolonisation.

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

BackgroundWinter 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.ResultsWe 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.ConclusionOverall 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.

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.

Local Adaptation at the Transcriptome Level in Brown Trout: Evidence from Early Life History Temperature Genomic Reaction Norms

Local adaptation and its underlying molecular basis has long been a key focus in evolutionary biology. There has recently been increased interest in the evolutionary role of plasticity and the molecular mechanisms underlying local adaptation. Using transcriptome analysis, we assessed differences in gene expression profiles for three brown trout (Salmo trutta) populations, one resident and two anadromous, experiencing different temperature regimes in the wild. The study was based on an F2 generation raised in a common garden setting. A previous study of the F1 generation revealed different reaction norms and significantly higher QST than FST among populations for two early life-history traits. In the present study we investigated if genomic reaction norm patterns were also present at the transcriptome level. Eggs from the three populations were incubated at two temperatures (5 and 8 degrees C) representing conditions encountered in the local environments. Global gene expression for fry at the stage of first feeding was analysed using a 32k cDNA microarray. The results revealed differences in gene expression between populations and temperatures and population × temperature interactions, the latter indicating locally adapted reaction norms. Moreover, the reaction norms paralleled those observed previously at early life-history traits. We identified 90 cDNA clones among the genes with an interaction effect that were differently expressed between the ecologically divergent populations. These included genes involved in immune- and stress response. We observed less plasticity in the resident as compared to the anadromous populations, possibly reflecting that the degree of environmental heterogeneity encountered by individuals throughout their life cycle will select for variable level of phenotypic plasticity at the transcriptome level. Our study demonstrates the usefulness of transcriptome approaches to identify genes with different temperature reaction norms. The responses observed suggest that populations may vary in their susceptibility to climate change.