STRONG SELECTION GENOME‐WIDE ENHANCES FITNESS TRADE‐OFFS ACROSS ENVIRONMENTS AND EPISODES OF SELECTION

Fitness trade‐offs across episodes of selection and environments influence life‐history evolution and adaptive population divergence. Documenting these trade‐offs remains challenging as selection can vary in magnitude and direction through time and space. Here, we evaluate fitness trade‐offs at the levels of the whole organism and the quantitative trait locus (QTL) in a multiyear field study of Boechera stricta (Brassicaceae), a genetically tractable mustard native to the Rocky Mountains. Reciprocal local adaptation was pronounced for viability, but not for reproductive components of fitness. Instead, local genomes had a fecundity advantage only in the high latitude garden. By estimating realized selection coefficients from individual‐level data on viability and reproductive success and permuting the data to infer significance, we examined the genetic basis of fitness trade‐offs. This analytical approach (Conditional Neutrality‐Antagonistic Pleiotropy, CNAP) identified genetic trade‐offs at a flowering phenology QTL (costs of adaptation) and revealed genetic trade‐offs across fitness components (costs of reproduction). These patterns would not have emerged from traditional ANOVA‐based QTL mapping. Our analytical framework can be applied to other systems to investigate fitness trade‐offs. This task is becoming increasingly important as climate change may alter fitness landscapes, potentially disrupting fitness trade‐offs that took many generations to evolve.     

Long‐term in situ persistence of biodiversity in tropical sky islands revealed by landscape genomics

Tropical mountains are areas of high species richness and endemism. Two historical phenomena may have contributed to this: (i) fragmentation and isolation of habitats may have promoted the genetic differentiation of populations and increased the possibility of allopatric divergence and speciation and (ii) the mountain areas may have allowed long‐term population persistence during global climate fluctuations. These two phenomena have been studied using either species occurrence data or estimating species divergence times. However, only few studies have used intraspecific genetic data to analyse the mechanisms by which endemism may emerge at the microevolutionary scale. Here, we use landscape analysis of genomic SNP data sampled from two high‐elevation plant species from an archipelago of tropical sky islands (the Trans‐Mexican Volcanic Belt) to test for population genetic differentiation, synchronous demographic changes and habitat persistence. We show that genetic differentiation can be explained by the degree of glacial habitat connectivity among mountains and that mountains have facilitated the persistence of populations throughout glacial/interglacial cycles. Our results support the ongoing role of tropical mountains as cradles for biodiversity by uncovering cryptic differentiation and limits to gene flow.     

Geographical patterns of turnover and nestedness-resultant components of allelic diversity among populations

The analysis of geographical patterns in population divergence has always been a powerful way to infer microevolutionary processes involved in population differentiation, and several approaches have been used to investigate such patterns. Most frequently, multivariate spatial patterns of population differentiation are analyzed by computing pairwise genetic distances or F_ST (or related statistics, such as ?_ST from AMOVA), which are then correlated with geographical distances or landscape features. However, when calculating distances, especially based on presence-absence of alleles in local populations, there would be a confounding effect of allelic richness differences in the population differentiation. Moreover, the relative magnitude of these components and their spatial patterns can help identifying microevolutionary processes driving population differentiation. Here we show how recent methodological advances in ecological community analyses that allows partitioning dissimilarity into turnover (turnover) and richness differences, or nestedness-resultant dissimilarity, can be applied to allelic variation data, using an endemic Cerrado tree (Dipteryx alata) as a case study. Individuals from 15 local populations were genotyped for eight microsatellite loci, and pairwise dissimilarities were computed based on presence-absence of alleles. The turnover of alleles among populations represented 69?% of variation in dissimilarity, but only the richness difference component shows a clear spatial structure, appearing as a westward decrease of allelic richness. We show that decoupling richness difference and turnover components of allelic variation reveals more clearly how similarity among populations reflects geographical patterns in allelic diversity that can be interpreted in respect to historical range expansion in the species.     

Microevolution in island forms: the roles of drift and directional selection in morphological divergence of a passerine bird

Abstract.— Theory predicts that in small isolated populations random genetic drift can lead to phenotypic divergence; however this prediction has rarely been tested quantitatively in natural populations. Here we utilize natural repeated island colonization events by members of the avian species complex, Zosterops lateralis , to assess whether or not genetic drift alone is an adequate explanation for the observed patterns of microevolutionary divergence in morphology. Morphological and molecular genetic characteristics of island and mainland populations are compared to test three predictions of drift theory: (1) that the pattern of morphological change is idiosyncratic to each island; (2) that there is concordance between morphological and neutral genetic shifts across island populations; and (3) for populations whose time of colonization is known, that the rate of morphological change is sufficiently slow to be accounted for solely by genetic drift. Our results are not consistent with these predictions. First, the direction of size shifts was consistently towards larger size, suggesting the action of a nonrandom process. Second, patterns of morphological divergence among recently colonized populations showed little concordance with divergence in neutral genetic characters. Third, rate tests of morphological change showed that effective population sizes were not small enough for random processes alone to account for the magnitude of microevolutionary change. Altogether, these three lines of evidence suggest that drift alone is not an adequate explanation of morphological differentiation in recently colonized island Zosterops and therefore we suggest that the observed microevolutionary changes are largely a result of directional natural selection.     

DNA turnover and the molecular clock

Summary Many detailed studies on the mechanisms by which different components of eukaryotic nuclear genomes have diverged reveal that the majority of sequences are seemingly not passively accumulating base substitutions in a clocklike manner solely determined by laws of diffusion at the population level. It appears that variation in the rates, units, biases, and gradients of several DNA turnover mechanisms are contributing to the course of DNA divergence. Turnover mechanisms have the potential to retard, maintain, or accelerate the rate of DNA differentiation between populations. Furthermore, examples are known of coding and noncoding DNA subject to the simultaneous operation of several turnover mechanisms leading to complex patterns of fine-scale restructuring and divergence, generally uninterpretable using selection and/or neutral drift arguments in isolation. Constancy in the rate of divergence, where observed over defined periods of time, could be a reflection of constancy in the rates and units of turnover. However, a consideration of the generally large disparity between rates of turnover and mutation reveals that DNA clocks, which would be independently driven by turnover in separate genomic components, would tend to be episodic. The utility of any given DNA sequence for measuring time and species relationships, like individual proteins, is proportional to the extent to which all contributing forces to the evolution of the sequence, internal and external, are understood.     

Disentangling the effects of historic vs. contemporary landscape structure on population genetic divergence

Increasing habitat fragmentation poses an immediate threat to population viability, as gene flow patterns are changed in these altered landscapes. Patterns of genetic divergence can potentially reveal the impact of these shifts in landscape connectivity. However, divergence patterns not only carry the signature of altered contemporary landscapes, but also historical ones. When considered separately, both recent and historical landscape structure appear to significantly affect connectivity among 51 wood frog ( Rana sylvatica ) populations. However, by controlling for correlations among landscape structure from multiple time periods, we show that patterns of genetic divergence reflect recent landscape structure as opposed to landscape structure prior to European settlement of the region (before 1850s). At the same time, within-population genetic diversities remain high and a genetic signature of population bottlenecks is lacking. Together, these results suggest that metapopulation processes – not drift-induced divergence associated with strong demographic bottlenecks following habitat loss – underlie the strikingly rapid consequences of temporally shifting landscape structure on these amphibians. We discuss the implications of these results in the context of understanding the role of population demography in the adaptive variation observed in wood frog populations.     

Genome‐wide epigenetic isolation by environment in a widespread Anolis lizard

Epigenetic changes can provide a pathway for organisms to respond to local environmental conditions by influencing gene expression. However, we still know little about the spatial distribution of epigenetic variation in natural systems, how it relates to the distribution of genetic variation and the environmental structure of the landscape, and the processes that generate and maintain it. Studies examining spatial patterns of genetic and epigenetic variation can provide valuable insights into how ecological and population processes contribute to epigenetic divergence across heterogeneous landscapes. Here, we perform a comparative analysis of spatial genetic and epigenetic variation based on 8,459 single nucleotide polymorphisms (SNPs) and 8,580 single methylation variants (SMVs) from eight populations of the Puerto Rican crested anole, Anolis cristatellus, an abundant lizard in the adaptive radiations of anoles on the Greater Antilles that occupies a diverse range of habitats. Using generalized dissimilarity modelling and multiple matrix regression, we found that genome‐wide epigenetic differentiation is strongly correlated with environmental divergence, even after controlling for the underlying genetic structure. We also detected significant associations between key environmental variables and 96 SMVs, including 42 located in promoter regions or gene bodies. Our results suggest an environmental basis for population‐level epigenetic differentiation in this system and contribute to better understanding how environmental gradients structure epigenetic variation in nature.     

Cross‐sex genetic correlations and the evolution of sex‐specific local adaptation: Insights from classical trait clines in Drosophila melanogaster

Natural selection varies widely among locations of a species’ range, favoring population divergence and adaptation to local environmental conditions. Selection also differs between females and males, favoring the evolution of sexual dimorphism. Both forms of within‐species evolutionary diversification are widely studied, though largely in isolation, and it remains unclear whether environmental variability typically generates similar or distinct patterns of selection on each sex. Studies of sex‐specific local adaptation are also challenging because they must account for genetic correlations between female and male traits, which may lead to correlated patterns of trait divergence between sexes, whether or not local selection patterns are aligned or differ between the sexes. We quantified sex‐specific divergence in five clinally variable traits in Drosophila melanogaster that individually vary in their magnitude of cross‐sex genetic correlation (i.e., from moderate to strongly positive). In all five traits, we observed parallel male and female clines, regardless of the magnitude of their genetic correlation. These patterns imply that parallel spatial divergence of female and male traits is a reflection of sexually concordant directional selection imposed by local environmental conditions. In such contexts, genetic correlations between the sexes promote, rather than constrain, local adaptation to a spatially variable environment.     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Population genetics of the olive‐winged bulbul (Pycnonotus plumosus) in a tropical urban‐fragmented landscape

With increasing urbanization, urban‐fragmented landscapes are becoming more and more prevalent worldwide. Such fragmentation may lead to small, isolated populations that face great threats from genetic factors that affect even avian species with high dispersal propensities. Yet few studies have investigated the population genetics of species living within urban‐fragmented landscapes in the Old World tropics, in spite of the high levels of deforestation and fragmentation within this region. We investigated the evolutionary history and population genetics of the olive‐winged bulbul (Pycnonotus plumosus) in Singapore, a highly urbanized island which retains <5% of its original forest cover in fragments. Combining our own collected and sequenced samples with those from the literature, we conducted phylogenetic and population genetic analyses. We revealed high genetic diversity, evidence for population expansion, and potential presence of pronounced gene flow across the population in Singapore. This suggests increased chances of long‐term persistence for the olive‐winged bulbul and the ecosystem services it provides within this landscape.     

Speciation in parasites: a population genetics approach

Parasite speciation and host-parasite coevolution should be studied at both macroevolutionary and microevolutionary levels. Studies on a macroevolutionary scale provide an essential framework for understanding the origins of parasite lineages and the patterns of diversification. However, because coevolutionary interactions can be highly divergent across time and space, it is important to quantify and compare the phylogeographic variation in both the host and the parasite throughout their geographical range. Furthermore, to evaluate demographic parameters that are relevant to population genetics structure, such as effective population size and parasite transmission, parasite populations must be studied using neutral genetic markers. Previous emphasis on larger-scale studies means that the connection between microevolutionary and macroevolutionary events is poorly explored. In this article, we focus on the spatial fragmentation of parasites and the population genetics processes behind their diversification in an effort to bridge the micro- and macro-scales.     

Genetic and epigenetic control of the efficiency and fidelity of cross‐species prion transmission

Self‐perpetuating amyloid‐based protein isoforms (prions) transmit neurodegenerative diseases in mammals and phenotypic traits in yeast. Although mechanisms that control species specificity of prion transmission are poorly understood, studies of closely related orthologues of yeast prion protein Sup35 demonstrate that cross‐species prion transmission is modulated by both genetic (specific sequence elements) and epigenetic (prion variants, or ‘strains’) factors. Depending on the prion variant, the species barrier could be controlled at the level of either heterologous co‐aggregation or conversion of the aggregate‐associated heterologous protein into a prion polymer. Sequence divergence influences cross‐species transmission of different prion variants in opposing ways. The ability of a heterologous prion domain to either faithfully reproduce or irreversibly switch the variant‐specific prion patterns depends on both sequence divergence and the prion variant. Sequence variations within different modules of prion domains contribute to transmission barriers in different cross‐species combinations. Individual amino acid substitutions within short amyloidogenic stretches drastically alter patterns of cross‐species prion conversion, implicating these stretches as major determinants of species specificity.     

Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

Population genomics of eusocial insects: the costs of a vertebrate‐like effective population size

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life‐history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome‐wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA‐seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes – a marker of the effective population size – in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.     

Synergistic effects of climate and land‐use change influence broad‐scale avian population declines

Climate and land‐use changes are expected to be the primary drivers of future global biodiversity loss. Although theory suggests that these factors impact species synergistically, past studies have either focused on only one in isolation or have substituted space for time, which often results in confounding between drivers. Tests of synergistic effects require congruent time series on animal populations, climate change and land‐use change replicated across landscapes that span the gradient of correlations between the drivers of change. Using a unique time series of high‐resolution climate (measured as temperature and precipitation) and land‐use change (measured as forest change) data, we show that these drivers of global change act synergistically to influence forest bird population declines over 29 years in the Pacific Northwest of the United States. Nearly half of the species examined had declined over this time. Populations declined most in response to loss of early seral and mature forest, with responses to loss of early seral forest amplified in landscapes that had warmed over time. In addition, birds declined more in response to loss of mature forest in areas that had dried over time. Climate change did not appear to impact populations in landscapes with limited habitat loss, except when those landscapes were initially warmer than the average landscape. Our results provide some of the first empirical evidence of synergistic effects of climate and land‐use change on animal population dynamics, suggesting accelerated loss of biodiversity in areas under pressure from multiple global change drivers. Furthermore, our findings suggest strong spatial variability in the impacts of climate change and highlight the need for future studies to evaluate multiple drivers simultaneously to avoid potential misattribution of effects.     

Niche specificity influences gene flow across fine‐scale habitat mosaics in Succulent Karoo plants

While the tempo of diversification in biodiversity hotspots has received much attention, the spatial scale of diversification has often been overlooked. Addressing this deficiency requires understanding the drivers of population divergence and the spatial scales at which they operate in species‐rich clades and ecosystems. South Africa's Succulent Karoo (SK) hotspot provides an excellent system for such research, being both compact (ca. 110,000 km²) and home to spectacular in‐situ radiations, such as the ruschioid Aizoaceae. Here we use GBS to document genetic structure in two co‐occurring ruschioid species, at both coarse (>10 km) and fine (<500 m) spatial scales. Where Ruschia burtoniae shows strong between‐population genetic differentiation and no gene flow, Conophytum calculus shows weak differentiation, with high levels of admixture suggesting recent or ongoing gene flow. Community analysis and transplant experiments reveal that R. burtoniae occupies a narrow, low‐pH edaphic niche, and at scales of a few hundred metres, areas of elevated genetic turnover correspond to patches of edaphically unsuitable habitat. In contrast, C. calculus occupies a broader niche and exhibits isolation‐by‐distance without a habitat effect. We suggest that edaphic specialisation, coupled with highly restricted seed and pollen dispersal in heterogeneous landscapes, has played a major role in driving rapid diversification at small spatial scales in this system. However, the contrasting patterns in our study species show that these factors do not influence all organisms uniformly, being strongly modulated by lineage‐specific traits that influence both the spatial scale of gene flow and habitat specificity.     

Effects of historical forest contraction on the phylogeographic structure of Australo‐Papuan populations of the red‐legged pademelon (Macropodidae: Thylogale stigmatica)

Geographic patterns of genetic variation are strongly influenced by historical changes in species habitats. Whether such patterns are common to co‐distributed taxa may depend on the extent to which species vary in ecology and vagility. We investigated whether broad‐scale phylogeographic patterns common to a number of small‐bodied vertebrate and invertebrate species in eastern Australian forests were reflected in the population genetic structure of an Australo‐Papuan forest marsupial, the red‐legged pademelon (Macropodidae: Thylogale stigmatica). Strong genetic structuring of mtDNA haplotypes indicated the persistence of T. stigmatica populations across eastern Australia and southern New Guinea in Pleistocene refugial areas consistent with those inferred from studies of smaller, poorly dispersing species. However, there was limited divergence of haplotypes across two known historical barriers in the northeastern Wet Tropics (Black Mountain Barrier) and coastal mideastern Queensland (Burdekin Gap) regions. Lack of divergence across these barriers may reflect post‐glacial recolonization of forests from a large, central refugium in the Wet Tropics. Additionally, genetic structure is not consistent with the present delimitation of subspecies T. s. wilcoxi and T. s. stigmatica across the Burdekin Gap. Instead, the genetic division occurs further to the south in mideastern Queensland. Thus, while larger‐bodied marsupials such as T. stigmatica did persist in Pleistocene refugia common to a number of other forest‐restricted species, species‐specific local extinction and recolonization events have resulted in cryptic patterns of genetic variation. Our study demonstrates the importance of understanding individualistic responses to historical climate change in order to adequately conserve genetic diversity and the evolutionary potential of species.     

Genetic effect of transportation infrastructure on roe deer populations (Capreolus capreolus)

Anthropogenic transportation infrastructure is a major factor of habitat fragmentation leading to genetic population fragmentation in wildlife. Assessing and understanding the impact of this deterministic factor on genetic diversity and divergence of populations is crucial to appraise the viability of wildlife populations in fragmented landscapes. In this study, the roe deer is used as an example species for the assessment of genetic differentiation of populations separated by an anthropogenic barrier. In order to detect genetic discontinuities, we screened 12 polymorphic microsatellites on 222 individuals out of 11 roe deer populations that were sampled on the east and the westside of a fenced motorway in Central Switzerland. The interaction between landscape structure and microevolutionary processes such as gene flow and drift were assessed and evaluated by different population genetic methods like F-statistics, Mantel test, spatial autocorrelation analyses, Monmonier algorithm, and principal component analysis in conjunction with geographic information system data (synthesis map). We revealed an influence of the transportation infrastructure on genetic divergence of the roe deer population examined, but no impact on genetic diversity was detected. Based on the achieved genetic findings, recommendations for management implementation were made.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.