QST–FST comparisons with unbalanced half‐sib designs

Q_ST, a measure of quantitative genetic differentiation among populations, is an index that can suggest local adaptation if Q_ST for a trait is sufficiently larger than the mean F_ST of neutral genetic markers. A previous method by Whitlock and Guillaume derived a simulation resampling approach to statistically test for a difference between Q_ST and F_ST, but that method is limited to balanced data sets with offspring related as half‐sibs through shared fathers. We extend this approach (i) to allow for a model more suitable for some plant populations or breeding designs in which offspring are related through mothers (assuming independent fathers for each offspring; half‐sibs by dam); and (ii) by explicitly allowing for unbalanced data sets. The resulting approach is made available through the R package QstFstComp.     

High genetic diversity and population differentiation in Boechera fecunda, a rare relative of Arabidopsis

Conservation of endangered species becomes a critical issue with the increasing rates of extinction. In this study, we use 13 microsatellite loci and 27 single-copy nuclear loci to investigate the population genetics of Boechera fecunda, a rare relative of Arabidopsis thaliana, known from only 21 populations in Montana. We investigated levels of genetic diversity and population structure in comparison to its widespread congener, Boechera stricta, which shares similar life history and mating system. Despite its rarity, B. fecunda had levels of genetic diversity similar to B. stricta for both microsatellites and nucleotide polymorphism. Populations of B. fecunda are highly differentiated, with a majority of genetic diversity existing among populations (F(ST) = 0.57). Differences in molecular diversity and allele frequencies between western and eastern population groups suggest they experienced very different evolutionary histories.     

Evolutionary dynamics of the leaf phenological cycle in an oak metapopulation along an elevation gradient

It has been predicted that environmental changes will radically alter the selective pressures on phenological traits. Long‐lived species, such as trees, will be particularly affected, as they may need to undergo major adaptive change over only one or a few generations. The traits describing the annual life cycle of trees are generally highly evolvable, but nothing is known about the strength of their genetic correlations. Tight correlations can impose strong evolutionary constraints, potentially hampering the adaptation of multivariate phenological phenotypes. In this study, we investigated the evolutionary, genetic and environmental components of the timing of leaf unfolding and senescence within an oak metapopulation along an elevation gradient. Population divergence, estimated from in situ and common‐garden data, was compared to expectations under neutral evolution, based on microsatellite markers. This approach made it possible (1) to evaluate the influence of genetic correlation on multivariate local adaptation to elevation and (2) to identify traits probably exposed to past selective pressures due to the colder climate at high elevation. The genetic correlation was positive but very weak, indicating that genetic constraints did not shape the local adaptation pattern for leaf phenology. Both spring and fall (leaf unfolding and senescence, respectively) phenology timings were involved in local adaptation, but leaf unfolding was probably the trait most exposed to climate change‐induced selection. Our data indicated that genetic variation makes a much smaller contribution to adaptation than the considerable plastic variation displayed by a tree during its lifetime. The evolutionary potential of leaf phenology is, therefore, probably not the most critical aspect for short‐term population survival in a changing climate.     

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. F_ST outlier analysis revealed that promoter, 5′‐UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within‐population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.     

Breaking RAD: an evaluation of the utility of restriction site‐associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site‐associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome‐scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq‐based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD‐tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best‐case modelling scenario, we found that RADseq using six‐ or eight‐base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD‐tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.     

Fine‐scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Geographic variation is commonly observed in plant resistance traits, where plant species might experience different selection pressure across a heterogeneous landscape. Arabidopsis halleri subsp. gemmifera is dimorphic for trichome production, generating two morphs, trichome‐producing (hairy) and trichomeless (glabrous) plants. Trichomes of A. halleri are known to confer resistance against the white butterfly, cabbage sawfly, and brassica leaf beetle, but not against flea beetles. We combined leaf damage, microclimate, and microsatellite loci data of 26 A. halleri populations in central Japan, to explore factors responsible for fine‐scale geographic variation in the morph frequency. We found that hairy plants were less damaged than glabrous plants within populations, but the among‐site variation was the most significant source of variation in the individual‐level damage. Fixation index () of a putative trichome locus exhibited a significant divergence along population‐level damage with an exception of an outlier population, inferring the local adaptation to herbivory. Notably, this outlier was a population wherein our previous study reported a balancing role of the brassica leaf beetle Phaedon brassicae on the morph frequency. This differentiation of the trichome locus was unrelated to neutral genetic differentiation (evaluated by of microsatellite loci) and meteorological factors (including temperature and solar radiation). The present findings, combined with those of our previous work, provide suggestive evidence that herbivore‐driven divergence and occasional outbreak of a specific herbivore have jointly contributed to the ecogeographic pattern in the frequency of two morphs.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Landscape genomic signatures indicate reduced gene flow and forest‐associated adaptive divergence in an endangered neotropical turtle

Human‐induced transformations of ecosystems usually result in fragmented populations subject to increased extinction risk. Fragmentation is also often associated with novel environmental heterogeneity, which in combination with restricted gene flow may increase the opportunity for local adaptation. To manage at‐risk populations in these landscapes, it is important to understand how gene flow is changing, and how populations respond to habitat loss. We conducted a landscape genomics analysis using Restriction‐site Associated DNA sequencing to investigate the evolutionary response of the critically endangered Dahl's Toad‐headed turtle (Mesoclemmys dahli) to severe habitat modification. The species has lost almost all of its natural habitat in the southwestern part of its range and about 70% in the northeast. Based on least cost path analysis across different resistance surfaces for 3,211 SNPs, we found that the landscape matrix is restricting gene flow, causing the fragmentation of the species into at least six populations. Genome scans and allele‐environment association analyses indicate that the population fragments in the deforested grasslands of the southwest are adaptively different from those in the more forested northeast. Populations in areas with no forest had low levels of adaptive genetic diversity and the fixation of ancestrally‐polymorphic SNPs, consistent with directional selection in this novel environment. Our results suggest that this forest‐stream specialist is adapting to pond‐grassland conditions, but it is also suffering from negative consequences of habitat loss, including genetic erosion, isolation, small effective population sizes, and inbreeding. We recommend gene flow restoration via genetic rescue to counteract these threats, and provide guidance for this strategy.     

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.     

StAMPP: an R package for calculation of genetic differentiation and structure of mixed‐ploidy level populations

Statistical Analysis of Mixed‐Ploidy Populations (StAMPP) is a freely available R package for calculation of population structure and differentiation based on single nucleotide polymorphism (SNP) genotype data from populations of any ploidy level, and/or mixed‐ploidy levels. StAMPP provides an advance on previous similar software packages, due to an ability to calculate pairwise F_ST values along with confidence intervals, Nei's genetic distance and genomic relationship matrixes from data sets of mixed‐ploidy level. The software code is designed to efficiently handle analysis of large genotypic data sets that are typically generated by high‐throughput genotyping platforms. Population differentiation studies using StAMPP are broadly applicable to studies of molecular ecology and conservation genetics, as well as animal and plant breeding.     

Investigation of the geographical scale of adaptive phenological variation and its underlying genetics in Arabidopsis thaliana

Despite the increasing number of genomic tools, identifying the genetics underlying adaptive complex traits remains challenging in the model species Arabidopsis thaliana. This is due, at least in part, to the lack of data on the geographical scale of adaptive phenotypic variation. The aims of this study were (i) to tease apart the historical roles of adaptive and nonselective processes in shaping phenological variation in A. thaliana in France and (ii) to gain insights into the spatial scale of adaptive variation by identifying the putative selective agents responsible for this selection. Forty‐nine natural stands from four climatically contrasted French regions were characterized (i) phenologically for six traits, (ii) genetically using 135 SNP markers and (iii) ecologically for 42 variables. Up to 63% of phenological variation could be explained by neutral genetic diversity. The remaining phenological variation displayed stronger associations with ecological variation within regions than among regions, suggesting the importance of local selective agents in shaping adaptive phenological variation. Although climatic conditions have often been suggested as the main selective agents acting on phenology in A. thaliana, both edaphic conditions and interspecific competition appear to be strong selective agents in some regions. In a first attempt to identify the genetics of phenological variation at different geographical scales, we phenotyped worldwide accessions and local polymorphic populations from the French RegMap in a genome‐wide association (GWA) mapping study. The genomic regions associated with phenological variation depended upon the geographical scale considered, stressing the need to account for the scale of adaptive phenotypic variation when choosing accession panels for GWAS.     

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST–FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST) and quantitative genetic (Q_ST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.     

Fungal evolutionary genomics provides insight into the mechanisms of adaptive divergence in eukaryotes

Fungi are ideal model organisms for dissecting the genomic bases of adaptive divergence in eukaryotes. They have simple morphologies and small genomes, occupy contrasting, well‐identified ecological niches and tend to have short generation times, and many are amenable to experimental approaches. Fungi also display diverse lifestyles, from saprotrophs to pathogens or mutualists, and they play extremely important roles in both ecosystems and human activities, as wood decayers, mycorrhizal fungi, lichens, endophytes, plant and animal pathogens, and in fermentation or drug production. We review here recent insights into the patterns and mechanisms of adaptive divergence in fungi, including sources of divergence, genomic variation and, ultimately, speciation. We outline the various ecological sources of divergent selection and genomic changes, showing that gene loss and changes in gene expression and in genomic architecture are important adaptation processes, in addition to the more widely recognized processes of amino acid substitution and gene duplication. We also review recent findings regarding the interspecific acquisition of genomic variation and suggesting an important role for introgression, hybridization and horizontal gene transfers (HGTs). We show that transposable elements can mediate several of these genomic changes, thus constituting important factors for adaptation. Finally, we review the consequences of divergent selection in terms of speciation, arguing that genetic incompatibilities may not be as widespread as generally thought and that pleiotropy between adaptation and reproductive isolation is an important route of speciation in fungal pathogens.     

High‐throughput sequencing of transposable element insertions suggests adaptive evolution of the invasive Asian tiger mosquito towards temperate environments

Invasive species represent unique opportunities to evaluate the role of local adaptation during colonization of new environments. Among these species, the Asian tiger mosquito, Aedes albopictus, is a threatening vector of several human viral diseases, including dengue and chikungunya, and raises concerns about the Zika fever. Its broad presence in both temperate and tropical environments has been considered the reflection of great “ecological plasticity.” However, no study has been conducted to assess the role of adaptive evolution in the ecological success of Ae. albopictus at the molecular level. In the present study, we performed a genomic scan to search for potential signatures of selection leading to local adaptation in one‐hundred‐forty field‐collected mosquitoes from native populations of Vietnam and temperate invasive populations of Europe. High‐throughput genotyping of transposable element insertions led to the discovery of more than 120,000 polymorphic loci, which, in their great majority, revealed a virtual absence of structure between the biogeographic areas. Nevertheless, 92 outlier loci showed a high level of differentiation between temperate and tropical populations. The majority of these loci segregate at high insertion frequencies among European populations, indicating that this pattern could have been caused by recent adaptive evolution events in temperate areas. An analysis of the overlapping and neighbouring genes highlighted several candidates, including diapause, lipid and juvenile hormone pathways.     

Large‐scale screening and characterisation of Lemna aequinoctialis and Spirodela polyrhiza strains for starch production

• Duckweed is considered a promising feedstock for bioethanol production due to its high biomass and starch production. Selection of duckweed strains with high starch accumulation is essential for application of duckweeds to bioethanol production. Geographic differentiation had a large influence on genetic diversity of duckweeds.
• Biomass production, starch content and starch amount in geographically isolated strains of 20 Lemna aequinoctialis and Spirodela polyrhiza were calculated to evaluate their potential for bioethanol production. The influence of different collection time, culture medium and NaCl concentration on starch accumulation of the best strains were analysed.
• The results showed that biomass production, starch content and starch production of duckweeds demonstrated clonal dependency. The best strain was L. aequinoctialis 6000, with biomass production of 15.38 ± 1.47 g m^?2, starch content of 28.68 ± 1.10% and starch production of 4.39 ± 0.25 g m^?2. Furthermore, starch content of L. aequinoctialis 6000 was highest after 8 h of light, tap water was the best medium for starch induction, and NaCl did not induce starch accumulation.
• This study suggests duckweed biomass production and starch production demonstrate clonal dependency, indicating that extensive clonal comparisons will be required to identify the most suitable isolates for duckweed selective breeding for bioethanol.

     

Using occupancy modelling to compare environmental DNA to traditional field methods for regional‐scale monitoring of an endangered aquatic species

Environmental DNA (eDNA) monitoring approaches promise to greatly improve detection of rare, endangered and invasive species in comparison with traditional field approaches. Herein, eDNA approaches and traditional seining methods were applied at 29 research locations to compare method‐specific estimates of detection and occupancy probabilities for endangered tidewater goby (Eucyclogobius newberryi). At each location, multiple paired seine hauls and water samples for eDNA analysis were taken, ranging from two to 23 samples per site, depending upon habitat size. Analysis using a multimethod occupancy modelling framework indicated that the probability of detection using eDNA was nearly double (0.74) the rate of detection for seining (0.39). The higher detection rates afforded by eDNA allowed determination of tidewater goby occupancy at two locations where they have not been previously detected and at one location considered to be locally extirpated. Additionally, eDNA concentration was positively related to tidewater goby catch per unit effort, suggesting eDNA could potentially be used as a proxy for local tidewater goby abundance. Compared to traditional field sampling, eDNA provided improved occupancy parameter estimates and can be applied to increase management efficiency across a broad spatial range and within a diversity of habitats.     

A trait‐based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life‐history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species‐level life‐history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G′_ST, and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G′_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.     

Formalizing the definition of meta‐analysis in Molecular Ecology

Meta‐analysis, the statistical synthesis of pertinent literature to develop evidence‐based conclusions, is relatively new to the field of molecular ecology, with the first meta‐analysis published in the journal Molecular Ecology in 2003 (Slate & Phua 2003). The goal of this article is to formalize the definition of meta‐analysis for the authors, editors, reviewers and readers of Molecular Ecology by completing a review of the meta‐analyses previously published in this journal. We also provide a brief overview of the many components required for meta‐analysis with a more specific discussion of the issues related to the field of molecular ecology, including the use and statistical considerations of Wright's F_ST and its related analogues as effect sizes in meta‐analysis. We performed a literature review to identify articles published as ‘meta‐analyses’ in Molecular Ecology, which were then evaluated by at least two reviewers. We specifically targeted Molecular Ecology publications because as a flagship journal in this field, meta‐analyses published in Molecular Ecology have the potential to set the standard for meta‐analyses in other journals. We found that while many of these reviewed articles were strong meta‐analyses, others failed to follow standard meta‐analytical techniques. One of these unsatisfactory meta‐analyses was in fact a secondary analysis. Other studies attempted meta‐analyses but lacked the fundamental statistics that are considered necessary for an effective and powerful meta‐analysis. By drawing attention to the inconsistency of studies labelled as meta‐analyses, we emphasize the importance of understanding the components of traditional meta‐analyses to fully embrace the strengths of quantitative data synthesis in the field of molecular ecology.