Hybrid zones and the speciation continuum in Heliconius butterflies

Tropical butterflies in the genus Heliconius have long been models in the study of the stages of speciation. Heliconius are unpalatable to predators, and many species are notable for multiple geographic populations with striking warning colour pattern differences associated with Müllerian mimicry. A speciation continuum is evident in Heliconius hybrid zones. Examples range from hybrid zones across which (a) there is little genetic differentiation other than at mimicry loci, but where hybrids are common, (b) to ‘bimodal‘ hybrid zones with strong genetic divergence and few hybrids, (c) through to ‘good’ sympatric species, with hybridization extremely rare or absent. Now, in this issue of Molecular Ecology, Arias et al. ( 2012 ) have found an intermediate case in Colombian Heliconius cydno showing evidence for assortative mating and molecular differences, but where hybrids are abundant.     

A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour‐intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra 2008 ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.     

Complex selection on a regulator of social cognition: Evidence of balancing selection,regulatory interactions and population differentiation in the prairie vole Avpr1a locus

Adaptive variation in social behaviour depends upon standing genetic variation, but we know little about how evolutionary forces shape genetic diversity relevant to brain and behaviour. In prairie voles (Microtus ochrogaster), variants at the Avpr1a locus predict expression of the vasopressin 1a receptor in the retrosplenial cortex (RSC), a brain region that mediates spatial and contextual memory; cortical V1aR abundance in turn predicts diversity in space use and sexual fidelity in the field. To examine the potential contributions of adaptive and neutral forces to variation at the Avpr1a locus, we explore sequence diversity at the Avpr1a locus and throughout the genome in two populations of wild prairie voles. First, we refine results demonstrating balancing selection at the locus by comparing the frequency spectrum of variants at the locus to a random sample of the genome. Next, we find that the four single nucleotide polymorphisms that predict high V1aR expression in the RSC are in stronger linkage disequilibrium than expected by chance despite high recombination among intervening variants, suggesting that epistatic selection maintains their association. Analysis of population structure and a haplotype network for two populations revealed that this excessive LD was unlikely to be due to admixture alone. Furthermore, the two populations differed considerably in the region shown to be a regulator of V1aR expression despite the extremely low levels of genomewide genetic differentiation. Together, our data suggest that complex selection on Avpr1a locus favours specific combinations of regulatory polymorphisms, maintains the resulting alleles at population‐specific frequencies, and may contribute to unique patterns of spatial cognition and sexual fidelity among populations.     

Pollen,wind and fire: how to investigate genetic effects of disturbance‐induced change in forest trees

Understanding the consequences of habitat disturbance on mating patterns although pollen and seed dispersal in forest trees has been a long‐standing theme of forest and conservation genetics. Forest ecosystems face global environmental pressures from timber exploitation to genetic pollution and climate change, and it is therefore essential to comprehend how disturbances may alter the dispersal of genes and their establishment in tree populations in order to formulate relevant recommendations for sustainable resource management practices and realistic predictions of potential adaptation to climate change by means of range shift or expansion (Kremer et al. 2012 ). However, obtaining reliable evidence of disturbance‐induced effects on gene dispersal processes from empirical evaluation of forest tree populations is difficult. Indeed, tree species share characteristics such as high longevity, long generation time and large reproductive population size, which may impede the experimenter's ability to assess parameters at the spatial and time scales at which any change may occur (Petit and Hampe 2006 ). It has been suggested that appropriate study designs should encompass comparison of populations before and after disturbance as well as account for demonstrated variation in conspecific density, that is, the spatial distribution of mates, and forest density, including all species and relating to alteration in landscape openness (Bacles & Jump 2011 ). However, more often than not, empirical studies aiming to assess the consequences of habitat disturbance on genetic processes in tree populations assume rather than quantify a change in tree densities in forests under disturbance and generally fail to account for population history, which may lead to inappropriate interpretation of a causal relationship between population genetic structure and habitat disturbance due to effects of unmonitored confounding variables (Gauzere et al. 2013). In this issue, Shohami and Nathan ( 2014 ) take advantage of the distinctive features of the fire‐adapted wind‐pollinated Aleppo pine Pinus halepensis (Fig. 1) to provide an elegant example of best practice. Thanks to long‐term monitoring of the study site, a natural stand in Israel, Shohami and Nathan witnessed the direct impact of habitat disturbance, here taking the shape of fire, on conspecific and forest densities and compared pre‐ and postdisturbance mating patterns estimated from cones of different ages sampled on the same surviving maternal individuals (Fig. 2). This excellent study design is all the more strong that Shohami and Nathan took further analytical steps to account for confounding variables, such as historical population genetic structure and possible interannual variation in wind conditions, thus giving high credibility to their findings of unequivocal fire‐induced alteration of mating patterns in P. halepensis. Most notably, the authors found, at the pollen pool level, a disruption of local genetic structure which, furthermore, they were able to attribute explicitly to enhanced pollen‐mediated gene immigration into the low‐density fire‐disturbed stand. This cleverly designed research provides a model approach to be followed if we are to advance our understanding of disturbance‐induced dispersal and genetic change in forest trees.     

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (F_ST ~0.1), over a short coastline (<60 km) in the North Sea‐Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post‐glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.     

Don't throw the baby out with the bathwater: identifying and mapping paralogs in salmonids

Many eukaryotic genomes contain a large fraction of gene duplicates (or paralogs) as a result of ancient or recent whole‐genome duplications (Ohno 1970 ; Jaillon et al. 2004 ; Kellis et al. 2004 ). Identifying paralogs with NGS data is a pervasive problem in both ancient polyploids and neopolyploids. Likewise, paralogs are often treated as a nuisance that has to be detected and removed (Everett et al. 2012 ). In this issue of Molecular Ecology Resources, Waples et al. ( 2015 ) show that exclusion might not be necessary and how we may miss out on important genomic information in doing so. They present a novel statistical approach to detect paralogs based on the segregation of RAD loci in haploid offspring and test their method by constructing linkage maps with and without these duplicated loci in chum salmon, Oncorhynchus keta (Fig.  1 ). Their linkage map including the resolved paralogs shows that these are mostly located in the distal regions of several linkage groups. Particularly intriguing is their finding that these homoeologous regions appear impoverished in transposable elements (TE). Given the role that TE play in genome remodelling, it is noteworthy that these elements are of low abundance in regions showing residual tetrasomic inheritance. This raises the question whether re‐diploidization is constrained in these regions and whether they might have a role to play in salmonid speciation. This study provides an original approach to identifying duplicated loci in species with a pedigree, as well as providing a dense linkage map for chum salmon, and interesting insights into the retention of gene duplicates in an ancient polyploid.     

Molecular population genetics of the melanic plumage polymorphism in Arctic skuas (Stercorarius parasiticus): evidence for divergent selection on plumage colour

The Arctic skua (Stercorarius parasiticus) is a classic example of an avian plumage polymorphism, with variation in melanin‐based ventral plumage coloration defining pale, intermediate and dark morphs in adults of both sexes. However, despite several decades of field research, there is an incomplete understanding of how the polymorphism in ventral plumage colour is maintained and the selective forces involved. Here, we investigate selection on a locus (MC1R) that is strongly associated with plumage colour variation in Arctic skuas using patterns of nucleotide variation and comparison to neutral loci (nuclear introns and mtDNA). We find that three linked nonsynonymous mutations in MC1R, including the single mutation described previously, are associated with plumage colour in the Arctic skua. The position of nonsynonymous mutations on a MC1R haplotype network implies that divergent selection drove the initial evolution of the colour morphs. Comparisons of F_STs of MC1R vs. nuclear introns among five skua populations differing in proportion of dark morphs along an approximate north–south cline reveal a signature of divergent selection on MC1R. In contrast, we find limited evidence for balancing selection on MC1R within populations, although the power is low. Our results provide strong evidence for both past and ongoing selection on MC1R, and, by implication, plumage colour in Arctic skuas. The results suggest that a fruitful avenue for future ecological studies will be analysis of selection on morphs in colonies at the extremes along the morph ratio cline.     

Climate change and the spread of vector‐borne diseases: using approximate Bayesian computation to compare invasion scenarios for the bluetongue virus vector Culicoides imicola in Italy

Bluetongue (BT) is a commonly cited example of a disease with a distribution believed to have recently expanded in response to global warming. The BT virus is transmitted to ruminants by biting midges of the genus Culicoides, and it has been hypothesized that the emergence of BT in Mediterranean Europe during the last two decades is a consequence of the recent colonization of the region by Culicoides imicola and linked to climate change. To better understand the mechanism responsible for the northward spread of BT, we tested the hypothesis of a recent colonization of Italy by C. imicola, by obtaining samples from more than 60 localities across Italy, Corsica, Southern France, and Northern Africa (the hypothesized source point for the recent invasion of C. imicola), and by genotyping them with 10 newly identified microsatellite loci. The patterns of genetic variation within and among the sampled populations were characterized and used in a rigorous approximate Bayesian computation framework to compare three competing historical hypotheses related to the arrival and establishment of C. imicola in Italy. The hypothesis of an ancient presence of the insect vector was strongly favoured by this analysis, with an associated P ≥ 99%, suggesting that causes other than the northward range expansion of C. imicola may have supported the emergence of BT in southern Europe. Overall, this study illustrates the potential of molecular genetic markers for exploring the assumed link between climate change and the spread of diseases.     

Next‐generation phylogenetics takes root

It has been a tumultuous 5 years in phylogeography and phylogenetics during which both fields have struggled to harness the power of next‐generation sequencing (NGS) (Ekblom & Galindo 2010 ; McCormack et al. 2012a ). Fortunately, several methodological approaches appear to be taking root. In this issue of Molecular Ecology, O'Neill et al. 2013 ) employ one such method – parallel tagged sequencing (PTS) – to elucidate the phylogeography of a tiger salamander (Ambystoma tigrinum) species complex. This study demonstrates a practical application of NGS on a scale appropriate (and not overkill) for most biologists interested in phylogeography (~100 loci for ~100 individuals), and their results highlight several analytical challenges that lie ahead for researchers employing NGS techniques.     

The genomics of adaptation,divergence and speciation: a congealing theory

In this issue, Flaxman et al. ( 2014 ) report the results of sophisticated whole‐genome simulations of speciation with gene flow, enhancing our understanding of the process by building on previous single‐locus, multilocus and analytical works. Their findings provide us with new insights about how genomes can diverge and the importance of statistical and chromosomal linkage in facilitating reproductive isolation. The authors characterize the conditions under which, even with high gene flow and weak divergent selection, reproductive isolation between populations can occur due to the emergent stochastic process of genomewide congealing, where numerous statistically or physically linked loci of small effect allow selection to limit effective migration rates. The initial congealing event can occur within a broad range conditions, and once initiated, the self‐reinforcing process leads to rapid divergence and ultimately two reproductively isolated populations. Flaxman et al.'s ( 2014 ) work is a valuable contribution to our understanding of speciation with gene flow and in making a more predictive field of evolutionary genomics and speciation.     

Population genomic analyses reveal possible drivers of population divergence

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. 2011 ), Saccharomyces uvarum (Almeida et al. 2014 ), Coccidioides spp. (Neafsey et al. 2010 ) and Cryptococcus gatti (Engelthaler et al. 2014 ). In this issue of Molecular Ecology, Branco et al. ( 2015 ) apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.     

Characterizing restriction enzyme‐associated loci in historic ragweed (Ambrosia artemisiifolia) voucher specimens using custom‐designed RNA probes

Population genetic studies of nonmodel organisms frequently employ reduced representation library (RRL) methodologies, many of which rely on protocols in which genomic DNA is digested by one or more restriction enzymes. However, because high molecular weight DNA is recommended for these protocols, samples with degraded DNA are generally unsuitable for RRL methods. Given that ancient and historic specimens can provide key temporal perspectives to evolutionary questions, we explored how custom‐designed RNA probes could enrich for RRL loci (Restriction Enzyme‐Associated Loci baits, or REALbaits). Starting with genotyping‐by‐sequencing (GBS) data generated on modern common ragweed (Ambrosia artemisiifolia L.) specimens, we designed 20 000 RNA probes to target well‐characterized genomic loci in herbarium voucher specimens dating from 1835 to 1913. Compared to shotgun sequencing, we observed enrichment of the targeted loci at 19‐ to 151‐fold. Using our GBS capture pipeline on a data set of 38 herbarium samples, we discovered 22 813 SNPs, providing sufficient genomic resolution to distinguish geographic populations. For these samples, we found that dilution of REALbaits to 10% of their original concentration still yielded sufficient data for downstream analyses and that a sequencing depth of ~7m reads was sufficient to characterize most loci without wasting sequencing capacity. In addition, we observed that targeted loci had highly variable rates of success, which we primarily attribute to similarity between loci, a trait that ultimately interferes with unambiguous read mapping. Our findings can help researchers design capture experiments for RRL loci, thereby providing an efficient means to integrate samples with degraded DNA into existing RRL data sets.     

Herbarium specimens reveal a historical shift in phylogeographic structure of common ragweed during native range disturbance

Invasive plants provide ample opportunity to study evolutionary shifts that occur after introduction to novel environments. However, although genetic characters pre‐dating introduction can be important determinants of later success, large‐scale investigations of historical genetic structure have not been feasible. Common ragweed (Ambrosia artemisiifolia L.) is an invasive weed native to North America that is known for its allergenic pollen. Palynological records from sediment cores indicate that this species was uncommon before European colonization of North America, and ragweed populations expanded rapidly as settlers deforested the landscape on a massive scale, later becoming an aggressive invasive with populations established globally. Towards a direct comparison of genetic structure now and during intense anthropogenic disturbance of the late 19th century, we sampled 45 natural populations of common ragweed across its native range as well as historical herbarium specimens collected up to 140 years ago. Bayesian clustering analyses of 453 modern and 473 historical samples genotyped at three chloroplast spacer regions and six nuclear microsatellite loci reveal that historical ragweed's spatial genetic structure mirrors both the palaeo‐record of Ambrosia pollen deposition and the historical pattern of agricultural density across the landscape. Furthermore, for unknown reasons, this spatial genetic pattern has changed substantially in the intervening years. Following on previous work relating morphology and genetic expression between plants collected from eastern North America and Western Europe, we speculate that the cluster associated with humans’ rapid transformation of the landscape is a likely source of these aggressive invasive populations.     

Demystifying the RAD fad

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart ( 2014 ) entitled ‘Recent novel approaches for population genomics data analysis’. Restriction‐site‐associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. ( 2007 ) and Baird et al. ( 2008 ) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. 2011 ), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart 2014 ). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. 2007 ; Baird et al. 2008 ), double digest RAD (ddRAD, Peterson et al. 2012 ), ezRAD (Toonen et al. 2013 ) and 2bRAD (Wang et al. 2012 ). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

Ecological genomics in full colour

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning et al. ( 2014 ) detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson et al. ( 2014 ) identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson et al. ( 2014 ) go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics – efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.     

Genetic evidence for the origin of Aedes aegypti,the yellow fever mosquito,in the southwestern Indian Ocean

Aedes aegypti is among the best‐studied mosquitoes due to its critical role as a vector of human pathogens and ease of laboratory rearing. Until now, this species was thought to have originated in continental Africa, and subsequently colonized much of the world following the establishment of global trade routes. However, populations of this mosquito on the islands in the southwestern Indian Ocean (SWIO), where the species occurs with its nearest relatives referred to as the Aegypti Group, have received little study. We re‐evaluated the evolutionary history of Ae. aegypti and these relatives, using three data sets: nucleotide sequence data, 18,489 SNPs and 12 microsatellites. We found that: (a) the Aegypti Group diverged 16 MYA (95% HPD: 7–28 MYA) from its nearest African/Asian ancestor; (b) SWIO populations of Ae. aegypti are basal to continental African populations; (c) after diverging 7 MYA (95% HPD: 4–15 MYA) from its nearest formally described relative (Ae. mascarensis), Ae. aegypti moved to continental Africa less than 85,000 years ago, where it recently (<1,000 years ago) split into two recognized subspecies Ae. aegypti formosus and a human commensal, Ae. aegypti aegypti; (d) the Madagascar samples form a clade more distant from all other Ae. aegypti than the named species Ae. mascarensis, implying that Madagascar may harbour a new cryptic species; and (e) there is evidence of introgression between Ae. mascarensis and Ae. aegypti on Réunion, and between the two subspecies elsewhere in the SWIO, a likely consequence of recent introductions of domestic Ae. aegypti aegypti from Asia.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

Two are better than one: combining landscape genomics and common gardens for detecting local adaptation in forest trees

Predicting likely species responses to an alteration of their local environment is key to decision‐making in resource management, ecosystem restoration and biodiversity conservation practice in the face of global human‐induced habitat disturbance. This is especially true for forest trees which are a dominant life form on Earth and play a central role in supporting diverse communities and structuring a wide range of ecosystems. In Europe, it is expected that most forest tree species will not be able to migrate North fast enough to follow the estimated temperature isocline shift given current predictions for rapid climate warming. In this context, a topical question for forest genetics research is to quantify the ability for tree species to adapt locally to strongly altered environmental conditions (Kremer et al. 2012 ). Identifying environmental factors driving local adaptation is, however, a major challenge for evolutionary biology and ecology in general but is particularly difficult in trees given their large individual and population size and long generation time. Empirical evaluation of local adaptation in trees has traditionally relied on fastidious long‐term common garden experiments (provenance trials) now supplemented by reference genome sequence analysis for a handful of economically valuable species. However, such resources have been lacking for most tree species despite their ecological importance in supporting whole ecosystems. In this issue of Molecular Ecology, De Kort et al. ( 2014 ) provide original and convincing empirical evidence of local adaptation to temperature in black alder, Alnus glutinosa L. Gaertn, a surprisingly understudied keystone species supporting riparian ecosystems. Here, De Kort et al. ( 2014 ) use an innovative empirical approach complementing state‐of‐the‐art landscape genomics analysis of A. glutinosa populations sampled in natura across a regional climate gradient with phenotypic trait assessment in a common garden experiment (Fig. 1 ). By combining the two methods, De Kort et al. ( 2014 ) were able to detect unequivocal association between temperature and phenotypic traits such as leaf size as well as with genetic loci putatively under divergent selection for temperature. The research by De Kort et al. ( 2014 ) provides valuable insight into adaptive response to temperature variation for an ecologically important species and demonstrates the usefulness of an integrated approach for empirical evaluation of local adaptation in nonmodel species (Sork et al. 2013 ).