The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia.

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.

CircPrime: a web-based platform for design of specific circular RNA primers.

BACKGROUND: Circular RNAs (circRNAs) are covalently closed-loop RNAs with critical regulatory roles in cells. Tens of thousands of circRNAs have been unveiled due to the recent advances in high throughput RNA sequencing technologies and bioinformatic tools development. At the same time, polymerase chain reaction (PCR) cross-validation for circRNAs predicted by bioinformatic tools remains an essential part of any circRNA study before publication. RESULTS: Here, we present the CircPrime web-based platform, providing a user-friendly solution for DNA primer design and thermocycling conditions for circRNA identification with routine PCR methods. CONCLUSIONS: User-friendly CircPrime web platform ( http://circprime.elgene.net/ ) works with outputs of the most popular bioinformatic predictors of circRNAs to design specific circular RNA primers. CircPrime works with circRNA coordinates and any reference genome from the National Center for Biotechnology Information database).

Context-dependent parasite infection affects trophic niche in populations of sympatric stickleback species.

How parasites alter host feeding ecology remains elusive in natural populations. A powerful approach to investigate the link between infection and feeding ecology is quantifying unique and shared responses to parasite infection in related host species within a common environment. Here, 9 pairs of sympatric populations of the three-spined and nine-spined stickleback fishes were sampled across a range of freshwater and brackish habitats to investigate how parasites alter host feeding ecology: (i) biotic and abiotic determinants of parasite community composition, and (ii) to what extent parasite infection correlates with trophic niche specialization of the 2 species, using stable isotope analyses (δ15N and δ13C). It was determined that parasite community composition and host parasite load varied among sites and species and were correlated with dissolved oxygen. It was also observed that the digenean Cyathocotyle sp.'s abundance, a common directly infecting parasite with a complex life cycle, correlated with host δ13C in a fish species-specific manner. In 6 sites, correlations were found between parasite abundance and their hosts' feeding ecology. These effects were location-specific and occasionally host species or host size-specific. Overall, the results suggest a relationship between parasite infection and host trophic niche which may be an important and largely overlooked ecological factor. The population specificity and variation in parasite communities also suggest this effect is multifarious and context-dependent.

Sea ice loss increases genetic isolation in a high Arctic ungulate metapopulation.

Sea ice loss may have dramatic consequences for population connectivity, extinction-colonization dynamics, and even the persistence of Arctic species subject to climate change. This is of particular concern in face of additional anthropogenic stressors, such as overexploitation. In this study, we assess the population-genetic implications of diminishing sea ice cover in the endemic, high Arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) by analyzing the interactive effects of landscape barriers and reintroductions (following harvest-induced extirpations) on their metapopulation genetic structure. We genotyped 411 wild reindeer from 25 sampling sites throughout the entire subspecies' range at 19 microsatellite loci. Bayesian clustering analysis showed a genetic structure composed of eight populations, of which two were admixed. Overall population genetic differentiation was high (mean FST  = 0.21). Genetic diversity was low (allelic richness [AR] = 2.07-2.58; observed heterozygosity = 0.23-0.43) and declined toward the outer distribution range, where populations showed significant levels of inbreeding. Coalescent estimates of effective population sizes and migration rates revealed strong evolutionary source-sink dynamics with the central population as the main source. The population genetic structure was best explained by a landscape genetics model combining strong isolation by glaciers and open water, and high connectivity by dispersal across winter sea ice. However, the observed patterns of natural isolation were strongly modified by the signature of past harvest-induced extirpations, subsequent reintroductions, and recent lack of sea ice. These results suggest that past and current anthropogenic drivers of metapopulation dynamics may have interactive effects on large-scale ecological and evolutionary processes. Continued loss of sea ice as a dispersal corridor within and between island systems is expected to increase the genetic isolation of populations, and thus threaten the evolutionary potential and persistence of Arctic wildlife.

First genomic study on Lake Tanganyika sprat Stolothrissa tanganicae: a lack of population structure calls for integrated management of this important fisheries target species.

BACKGROUND: Clupeid fisheries in Lake Tanganyika (East Africa) provide food for millions of people in one of the world's poorest regions. Due to climate change and overfishing, the clupeid stocks of Lake Tanganyika are declining. We investigate the population structure of the Lake Tanganyika sprat Stolothrissa tanganicae, using for the first time a genomic approach on this species. This is an important step towards knowing if the species should be managed separately or as a single stock. Population structure is important for fisheries management, yet understudied for many African freshwater species. We hypothesize that distinct stocks of S. tanganicae could be present due to the large size of the lake (isolation by distance), limnological variation (adaptive evolution), or past separation of the lake (historical subdivision). On the other hand, high mobility of the species and lack of obvious migration barriers might have resulted in a homogenous population. RESULTS: We performed a population genetic study on wild-caught S. tanganicae through a combination of mitochondrial genotyping (96 individuals) and RAD sequencing (83 individuals). Samples were collected at five locations along a north-south axis of Lake Tanganyika. The mtDNA data had low global FST and, visualised in a haplotype network, did not show phylogeographic structure. RAD sequencing yielded a panel of 3504 SNPs, with low genetic differentiation (FST = 0.0054; 95% CI: 0.0046-0.0066). PCoA, fineRADstructure and global FST suggest a near-panmictic population. Two distinct groups are apparent in these analyses (FST = 0.1338 95% CI: 0.1239,0.1445), which do not correspond to sampling locations. Autocorrelation analysis showed a slight increase in genetic difference with increasing distance. No outlier loci were detected in the RADseq data. CONCLUSION: Our results show at most very weak geographical structuring of the stock and do not provide evidence for genetic adaptation to historical or environmental differences over a north-south axis. Based on these results, we advise to manage the stock as one population, integrating one management strategy over the four riparian countries. These results are a first comprehensive study on the population structure of these important fisheries target species, and can guide fisheries management.

Temporally varying disruptive selection in the medium ground finch (Geospiza fortis).

Disruptive natural selection within populations exploiting different resources is considered to be a major driver of adaptive radiation and the production of biodiversity. Fitness functions, which describe the relationships between trait variation and fitness, can help to illuminate how this disruptive selection leads to population differentiation. However, a single fitness function represents only a particular selection regime over a single specified time period (often a single season or a year), and therefore might not capture longer-term dynamics. Here, we build a series of annual fitness functions that quantify the relationships between phenotype and apparent survival. These functions are based on a 9-year mark-recapture dataset of over 600 medium ground finches (Geospiza fortis) within a population bimodal for beak size. We then relate changes in the shape of these functions to climate variables. We find that disruptive selection between small and large beak morphotypes, as reported previously for 2 years, is present throughout the study period, but that the intensity of this selection varies in association with the harshness of environment. In particular, we find that disruptive selection was strongest when precipitation was high during the dry season of the previous year. Our results shed light on climatic factors associated with disruptive selection in Darwin's finches, and highlight the role of temporally varying fitness functions in modulating the extent of population differentiation.

Genomic variation at the tips of the adaptive radiation of Darwin's finches.

Adaptive radiation unfolds as selection acts on the genetic variation underlying functional traits. The nature of this variation can be revealed by studying the tips of an ongoing adaptive radiation. We studied genomic variation at the tips of the Darwin's finch radiation; specifically focusing on polymorphism within, and variation among, three sympatric species of the genus Geospiza. Using restriction site-associated DNA (RAD-seq), we characterized 32 569 single-nucleotide polymorphisms (SNPs), from which 11 outlier SNPs for beak and body size were uncovered by a genomewide association study (GWAS). Principal component analysis revealed that these 11 SNPs formed four statistically linked groups. Stepwise regression then revealed that the first PC score, which included 6 of the 11 top SNPs, explained over 80% of the variation in beak size, suggesting that selection on these traits influences multiple correlated loci. The two SNPs most strongly associated with beak size were near genes associated with beak morphology across deeper branches of the radiation: delta-like 1 homologue (DLK1) and high-mobility group AT-hook 2 (HMGA2). Our results suggest that (i) key adaptive traits are associated with a small fraction of the genome (11 of 32 569 SNPs), (ii) SNPs linked to the candidate genes are dispersed throughout the genome (on several chromosomes), and (iii) micro- and macro-evolutionary variation (roots and tips of the radiation) involve some shared and some unique genomic regions.

Recurrent adaptation in a low-dispersal trait.

The study of natural populations from contrasting environments has greatly enhanced our understanding of ecological-dependent selection, adaptation and speciation. Cases of parallel evolution in particular have facilitated the study of the molecular and genetic basis of adaptive variation. This includes the type and number of genes underlying adaptive traits, as well as the extent to which these genes are exchanged among populations and contribute repeatedly to parallel evolution. Yet, surprisingly few studies provide a comprehensive view on the evolutionary history of adaptive traits from mutation to widespread adaptation. When did key mutations arise, how did they increase in frequency, and how did they spread? In this issue of Molecular Ecology, Van Belleghem et al. (2015) reconstruct the evolutionary history of a gene associated with wing size in the salt marsh beetle Pogonus chalceus. Screening the entire distribution range of this species, they found a single origin for the allele associated with the short-winged ecotype. This allele seemingly evolved in an isolated population and rapidly introgressed into other populations. These findings suggest that the adaptive genetic variation found in sympatric short- and long-winged populations has an allopatric origin, confirming that allopatric phases may be important at early stages of speciation.

A gene with major phenotypic effects as a target for selection vs. homogenizing gene flow.

Genes with major phenotypic effects facilitate quantifying the contribution of genetic vs. plastic effects to adaptive divergence. A classical example is Ectodysplasin (Eda), the major gene controlling lateral plate phenotype in three-spined stickleback. Completely plated marine stickleback populations evolved repeatedly towards low-plated freshwater populations, representing a prime example of parallel evolution by natural selection. However, many populations remain polymorphic for lateral plate number. Possible explanations for this polymorphism include relaxation of selection, disruptive selection or a balance between divergent selection and gene flow. We investigated 15 polymorphic stickleback populations from brackish and freshwater habitats in coastal North-western Europe. At each site, we tracked changes in allele frequency at the Eda gene between subadults in fall, adults in spring and juveniles in summer. Eda genotypes were also compared for body size and reproductive investment. We observed a fitness advantage for the Eda allele for the low morph in freshwater and for the allele for the complete morph in brackish water. Despite these results, the differentiation at the Eda gene was poorly correlated with habitat characteristics. Neutral population structure was the best predictor of spatial variation in lateral plate number, suggestive of a substantial effect of gene flow. A meta-analysis revealed that the signature of selection at Eda was weak compared to similar studies in stickleback. We conclude that a balance between divergent selection and gene flow can maintain stickleback populations polymorphic for lateral plate number and that ecologically relevant genes may not always contribute much to local adaptation, even when targeted by selection.

Thyasirid species composition (Bivalvia: Thyasiridae) and genetic connectivity of Parathyasira equalis (A. E. Verrill & K. J. Bush, 1898) in deep basins of sub-Arctic fjords.

BACKGROUND: Thyasirid bivalves are often recorded as a dominant component of macrobenthic infaunal communities in depositional environments such as fjord basins. Fjord basins comprise patchy soft-bottom habitats bounded by steep walls and sills; however, little is known how this semi-isolated nature of fjords affects benthic populations. Accordingly, data on the composition and population connectivity of thyasirids can provide valuable information on the ecology of these ecosystems. RESULTS: The species composition of thyasirid bivalves has been studied in the basins of three sub-Arctic fjords (Nordland, Northern Norway). Overall, six thyasirid species were recorded: Parathyasira equalis, Parathyasira dunbari, Mendicula ferruginosa, Genaxinus eumyarius, Thyasira sarsii, and Thyasira obsoleta. The species composition remained stable within the basins during the sampling period (2013-2020) and suggested the importance of local reproduction over advection of individuals for population dynamics. Only one species, Parathyasira equalis, was common in all fjords. We have further investigated the population genetics of this species by combining two types of genetic markers: a 579 bp fragment of the cytochrome c oxidase subunit I (COI) gene and 4043 single-nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing. The latter provided a more in-depth resolution on the population genetics of this species and revealed a weak but significant differentiation of populations within fjords, further indicating limited connectivity between basins. CONCLUSION: Based on our findings, we conclude that limited dispersal between the basin communities results in weakly connected populations and might be an important structuring factor for macrobenthic communities.

Embryonic temperature has long-term effects on muscle circRNA expression and somatic growth in Nile tilapia.

Embryonic temperature has a lasting impact on muscle phenotype in vertebrates, involving complex molecular mechanisms that encompass both protein-coding and non-coding genes. Circular RNAs (circRNAs) are a class of regulatory RNAs that play important roles in various biological processes, but the effect of variable thermal conditions on the circRNA transcriptome and its long-term impact on muscle growth plasticity remains largely unexplored. To fill this knowledge gap, we performed a transcriptomic analysis of circRNAs in fast muscle of Nile tilapia (Oreochromis niloticus) subjected to different embryonic temperatures (24°C, 28°C and 32°C) and then reared at a common temperature (28°C) for 4 months. Nile tilapia embryos exhibited faster development and subsequently higher long-term growth at 32°C compared to those reared at 28°C and 24°C. Next-generation sequencing data revealed a total of 5,141 unique circRNAs across all temperature groups, of which 1,604, 1,531, and 1,169 circRNAs were exclusively found in the 24°C, 28°C and 32°C groups, respectively. Among them, circNexn exhibited a 1.7-fold (log2) upregulation in the 24°C group and a 1.3-fold (log2) upregulation in the 32°C group when compared to the 28°C group. Conversely, circTTN and circTTN_b were downregulated in the 24°C groups compared to their 28°C and 32°C counterparts. Furthermore, these differentially expressed circRNAs were found to have multiple interactions with myomiRs, highlighting their potential as promising candidates for further investigation in the context of muscle growth plasticity. Taken together, our findings provide new insights into the molecular mechanisms that may underlie muscle growth plasticity in response to thermal variation in fish, with important implications in the context of climate change, fisheries and aquaculture.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts.

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Contrasting parasite communities among allopatric colour morphs of the Lake Tanganyika cichlid Tropheus.

BACKGROUND: Adaptation to different ecological environments is thought to drive ecological speciation. This phenomenon culminates in the radiations of cichlid fishes in the African Great Lakes. Multiple characteristic traits of cichlids, targeted by natural or sexual selection, are considered among the driving factors of these radiations. Parasites and pathogens have been suggested to initiate or accelerate speciation by triggering both natural and sexual selection. Three prerequisites for parasite-driven speciation can be inferred from ecological speciation theory. The first prerequisite is that different populations experience divergent infection levels. The second prerequisite is that these infection levels cause divergent selection and facilitate adaptive divergence. The third prerequisite is that parasite-driven adaptive divergence facilitates the evolution of reproductive isolation. Here we investigate the first and the second prerequisite in allopatric chromatically differentiated lineages of the rock-dwelling cichlid Tropheus spp. from southern Lake Tanganyika (Central Africa). Macroparasite communities were screened in eight populations belonging to five different colour morphs. RESULTS: Parasite communities were mainly composed of acanthocephalans, nematodes, monogeneans, copepods, branchiurans, and digeneans. In two consecutive years (2011 and 2012), we observed significant variation across populations for infection with acanthocephalans, nematodes, monogeneans of the genera Gyrodactylus and Cichlidogyrus, and the copepod Ergasilus spp. Overall, parasite community composition differed significantly between populations of different colour morphs. Differences in parasite community composition were stable in time. The genetic structure of Tropheus populations was strong and showed a significant isolation-by-distance pattern, confirming that spatial isolation is limiting host dispersal. Correlations between parasite community composition and Tropheus genetic differentiation were not significant, suggesting that host dispersal does not influence parasite community diversification. CONCLUSIONS: Subject to alternating episodes of isolation and secondary contact because of lake level fluctuations, Tropheus colour morphs are believed to accumulate and maintain genetic differentiation through a combination of vicariance, philopatric behaviour and mate discrimination. Provided that the observed contrasts in parasitism facilitate adaptive divergence among populations in allopatry (which is the current situation), and promote the evolution of reproductive isolation during episodes of sympatry, parasites might facilitate speciation in this genus.

Damped long-term host-parasite Red Queen coevolutionary dynamics: a reflection of dilution effects?

An increase in biological diversity leads to a greater stability of ecosystem properties. For host-parasite interactions, this is illustrated by the 'dilution effect': a negative correlation between host biodiversity and disease risk. We show that a similar mechanism might stabilise host-parasite dynamics at a lower level of diversity, i.e. at the level of genetic diversity within host species. A long-term time shift experiment, based on a historical reconstruction of a Daphnia-parasite coevolution, reveals infectivity cycles with more stable amplitude in experienced than in naive hosts. Coevolutionary models incorporating an increase in host allelic diversity over time explain the detected asymmetry. The accumulation of resistance alleles creates an opportunity for the host to stabilise Red Queen dynamics. It leads to a larger arsenal enhancing the host performance in its coevolution with the parasite in which 'it takes all the running both antagonists can do to keep in the same place'.

Differential modes of MHC class IIB gene evolution in cichlid fishes.

Cichlid fishes are emblematic models for the study of adaptive radiation, driven by natural and sexual selection. Parasite mediated selection is an important component in these processes, and the evolution of their immune system therefore merits special attention. In this study, light is shed on the phylogeny of the b family of cichlid major histocompatibility complex (MHC) class IIB genes. Full-length coding sequences were used to reconstruct phylogenies using criteria of maximum parsimony, maximum likelihood and Bayesian inference. All analyses suggest monophyly of the b family of cichlid MHC class IIB genes, although sequences of the cichlid sister taxa are currently not available. Two evolutionary lineages of these genes, respectively encompassing the recently defined genomic regions DBB-DEB-DFB and DCB-DDB, show highly contrasting levels of differentiation. To explore putative causes for these differences, exon 2 sequences were screened for variation in recombination rate and strength of selection. The more diversified lineage of cichlid MHC class IIB b genes was found to have higher levels of both recombination and selection. This is consistent with the observation in other taxa that recombination facilitates the horizontal spread of positively selected sites across MHC loci and hence contributes to fast sequence evolution. In contrast, the lineage that showed low diversification might either be under stabilizing selection or is evolutionary constrained by its low recombination rate. We speculate whether this lineage might include MHC genes with non-classical functions.

Adding parasites to the guppy-predation story: insights from field surveys.

Studies of phenotypic variation in nature often consider only a single potential selective agent. In such cases, it remains an open question as to whether variation attributed to that single measured agent might be influenced by some other unmeasured agent. Previous research has shown that phenotypic variation in the Trinidadian guppy (Poecilia reticulata) is strongly influenced by predation regime, and we here ask whether parasitism might represent an additional important selective agent shaping this variation. We performed a field survey of 26 natural guppy populations of known predation regime in northern Trinidad. We quantified levels of parasitism of guppies by the monogenean ecotoparasite, Gyrodactylus, and examined whether this parasite was associated with guppy body size or male colour. Spatial variation in Gyrodactylus parasitism was consistent between years, and parasite prevalence was generally, but not always, higher at high-predation sites than at low-predation sites. Consistent with previous work, predation regime was related to guppy size and some aspects of male colour, whereas parasitism showed few and only minor associations with the same traits. Moreover, a consideration of parasitism did not alter any interpretations regarding associations between guppy traits and predation regimes. These results suggest that parasitism, at least as quantified in the present study, does not play a major role in shaping variation in guppy body size or colour. Nevertheless, considerable variation in these traits, even within a predation regime, suggests the likely importance of other selective agents beyond just predation regime.

Host-parasite 'Red Queen' dynamics archived in pond sediment.

Antagonistic interactions between hosts and parasites are a key structuring force in natural populations, driving coevolution. However, direct empirical evidence of long-term host-parasite coevolution, in particular 'Red Queen' dynamics--in which antagonistic biotic interactions such as host-parasite interactions can lead to reciprocal evolutionary dynamics--is rare, and current data, although consistent with theories of antagonistic coevolution, do not reveal the temporal dynamics of the process. Dormant stages of both the water flea Daphnia and its microparasites are conserved in lake sediments, providing an archive of past gene pools. Here we use this fact to reconstruct rapid coevolutionary dynamics in a natural setting and show that the parasite rapidly adapts to its host over a period of only a few years. A coevolutionary model based on negative frequency-dependent selection, and designed to mimic essential aspects of our host-parasite system, corroborated these experimental results. In line with the idea of continuing host-parasite coevolution, temporal variation in parasite infectivity changed little over time. In contrast, from the moment the parasite was first found in the sediments, we observed a steady increase in virulence over time, associated with higher fitness of the parasite.

Host habitat rather than evolutionary history explains gut microbiome diversity in sympatric stickleback species.

Host-associated microbiota can influence host phenotypic variation, fitness and potential to adapt to local environmental conditions. In turn, both host evolutionary history and the abiotic and biotic environment can influence the diversity and composition of microbiota. Yet, to what extent environmental and host-specific factors drive microbial diversity remains largely unknown, limiting our understanding of host-microbiome interactions in natural populations. Here, we compared the intestinal microbiota between two phylogenetically related fishes, the three-spined stickleback (Gasterosteus aculeatus) and the nine-spined stickleback (Pungitius pungitius) in a common landscape. Using amplicon sequencing of the V3-V4 region of the bacterial 16S rRNA gene, we characterised the α and ß diversity of the microbial communities in these two fish species from both brackish water and freshwater habitats. Across eight locations, α diversity was higher in the nine-spined stickleback, suggesting a broader niche use in this host species. Habitat was a strong determinant of ß diversity in both host species, while host species only explained a small fraction of the variation in gut microbial composition. Strong habitat-specific effects overruled effects of geographic distance and historical freshwater colonisation, suggesting that the gut microbiome correlates primarily with local environmental conditions. Interestingly, the effect of habitat divergence on gut microbial communities was stronger in three-spined stickleback than in nine-spined stickleback, possibly mirroring the stronger level of adaptive divergence in this host species. Overall, our results show that microbial communities reflect habitat divergence rather than colonisation history or dispersal limitation of host species.

The fitness landscape of a community of Darwin's finches.

Divergent natural selection should lead to adaptive radiation-that is, the rapid evolution of phenotypic and ecological diversity originating from a single clade. The drivers of adaptive radiation have often been conceptualized through the concept of "adaptive landscapes," yet formal empirical estimates of adaptive landscapes for natural adaptive radiations have proven elusive. Here, we use a 17-year dataset of Darwin's ground finches (Geospiza spp.) at an intensively studied site on Santa Cruz (Galápagos) to estimate individual apparent lifespan in relation to beak traits. We use these estimates to model a multi-species fitness landscape, which we also convert to a formal adaptive landscape. We then assess the correspondence between estimated fitness peaks and observed phenotypes for each of five phenotypic modes (G. fuliginosa, G. fortis [small and large morphotypes], G. magnirostris, and G. scandens). The fitness and adaptive landscapes show 5 and 4 peaks, respectively, and, as expected, the adaptive landscape was smoother than the fitness landscape. Each of the five phenotypic modes appeared reasonably close to the corresponding fitness peak, yet interesting deviations were also documented and examined. By estimating adaptive landscapes in an ongoing adaptive radiation, our study demonstrates their utility as a quantitative tool for exploring and predicting adaptive radiation.

Quantifying population structure on short timescales.

Quantifying the contribution of the various processes that influence population genetic structure is important, but difficult. One of the reasons is that no single measure appropriately quantifies all aspects of genetic structure. An increasing number of studies is analysing population structure using the statistic D, which measures genetic differentiation, next to G(ST) , which quantifies the standardized variance in allele frequencies among populations. Few studies have evaluated which statistic is most appropriate in particular situations. In this study, we evaluated which index is more suitable in quantifying postglacial divergence between three-spined stickleback (Gasterosteus aculeatus) populations from Western Europe. Population structure on this short timescale (10 000 generations) is probably shaped by colonization history, followed by migration and drift. Using microsatellite markers and anticipating that D and G(ST) might have different capacities to reveal these processes, we evaluated population structure at two levels: (i) between lowland and upland populations, aiming to infer historical processes; and (ii) among upland populations, aiming to quantify contemporary processes. In the first case, only D revealed clear clusters of populations, putatively indicative of population ancestry. In the second case, only G(ST) was indicative for the balance between migration and drift. Simulations of colonization and subsequent divergence in a hierarchical stepping stone model confirmed this discrepancy, which becomes particularly strong for markers with moderate to high mutation rates. We conclude that on short timescales, and across strong clines in population size and connectivity, D is useful to infer colonization history, whereas G(ST) is sensitive to more recent demographic events.