Gene flow and population differentiation

Emotional Profiles of Anxiety, Depression, and Stress: Differences in School Anxiety

IntroductionThis study aims to confirm the existence of profiles according to the combinations of anxiety, depression, and stress, and looks to examine the differences between profiles according to the mean scores obtained in school anxiety.MethodsA total of 1,234 Spanish students at the secondary education level with an age range of 13-16 years old (M = 14.52; SD = 1.24) participated in the study by completing the abbreviated version of the Depression, Anxiety, and Stress Scale (DASS-21) and the School Anxiety Inventory.ResultsThe results showed positive, statistically significant, and moderate-sized correlations between all the variables analyzed. The Latent Profile Analysis identified four distinct profiles of depression, anxiety, and stress: Low DAS, Moderate DAS, High DAS, and Very High DAS. The results of the MANOVA showed statistically significant differences between these profiles regarding the school anxiety dimensions, with the profiles Very High DAS and Low DAS being the ones that reported, respectively, the highest and lowest levels in all the school anxiety components. Post hoc analyses revealed significant differences for the large part of profile comparisons, with there being large and moderate differences observed in the majority of cases (d = .30 and 1.66).ConclusionsThe results show the importance of considering social anxiety as a construct that is strongly associated with emotional problems such as depression, anxiety, and stress when developing effective actions to detect them and intervene with adolescents.

Extensive Admixture Among Karst-Obligate Salamanders Reveals Evidence of Recent Divergence and Gene Exchange Through Aquifers

Karst ecosystems often contain extraordinary biodiversity, but the complex underground aquifers of karst regions present challenges for assessing and conserving stygobiont diversity and investigating their evolutionary history. We examined the karst-obligate salamanders of the Eurycea neotenes species complex in the Edwards Plateau region of central Texas using population genomics data to address questions about population connectivity and the potential for gene exchange within the underlying aquifer system. The E. neotenes species complex has historically been divided into three nominal species, but their status, and spatial extent of species ranges, have remained uncertain. We discovered evidence of extensive admixture among species within the complex and with adjacent lineages. We observed relatively low levels of differentiation among all sampling localities which supports the hypothesis of recent divergence. Nominal taxonomy, aquifer region and geography each accounted for a modest amount of the overall population genomic variation; however, these predictors were largely confounded and difficult to disentangle. Importantly, current taxonomy of the three nominal species does not reflect the admixture apparent in clustering analyses. Inference of migration events revealed a complex pattern of gene exchange, suggesting that Eurycea salamanders have a dynamic history of dispersal through the aquifer system. These results highlight the need for greater understanding of how stygobiont populations are connected via dispersal and gene exchange through karst aquifers. These results also highlight the applicability of population genomics data as a powerful lever for investigating connectivity among populations in systems where direct detection of dispersal paths is difficult, as in underground, aquatic systems.

Ethnoprimatology reveals new extended distribution of critically endangered banded langur Presbytis femoralis (Martin, 1838) in Pahang, Malaysia: Insights from indigenous traditional knowledge and molecular analysis

The banded langur (Presbytis femoralis) is a critically endangered primate, which within Malaysia has not been known to extend significantly outside the state of Johor. Traditionally, distribution studies on this highly threatened primate have relied on conventional methods such as DNA identification, live counting, and camera trapping. However, ethnoprimatology offers an alternative approach to data collection, involving the active participation of indigenous and local communities possessing valuable knowledge and experience with local primate species. This study employed an integrated approach incorporating ethnoprimatology by utilizing pooled local expert opinion, local surveys, interviews, and fecal DNA analysis, resulting in a novel distribution range for the banded langur. The combination of expert opinions revealed this species' most optimistic distribution scenario across Johor and Pahang, inhabiting various ecosystems, including lowland forests, peat swamps, and human-modified landscapes. Further interviews and surveys conducted within the Orang Asli community in Tasik Chini and Tasek Bera have provided additional support for the revised distribution, documenting occurrences of banded langur utilization in indigenous practices, such as food consumption, cultural beliefs, medicinal applications, and craftsmanship. Phylogenetic analysis demonstrated genetic differentiation between populations in Johor and Pahang, with the populations in the southern part of Peninsular Malaysia likely serving as ancestral sources for other populations. Consequently, this study not only elucidated the updated distribution of banded langur through DNA records and direct observations but also established the efficacy of ethnoprimatology as a precursory tool for uncovering the present distribution patterns of other primate species in Malaysia.

Population structure and natural selection across a flower color polymorphism in the desert plant Encelia farinosa

PREMISE

Clines-or the geographic sorting of phenotypes across continual space-provide an opportunity to understand the interaction of dispersal, selection, and history in structuring polymorphisms.

METHODS

In this study, we combine field-sampling, genetics, climatic analyses, and machine learning to understand a flower color polymorphism in the wide-ranging desert annual Encelia farinosa.

RESULTS

We find evidence for replicated transitions in disk floret color from brown to yellow across spatial scales, with the most prominent cline stretching ~100 km from southwestern United States into México. Because population structure across the cline is minimal, selection is more likely than drift to have an important role in determining cline width.

CONCLUSIONS

Given that the cline aligns with a climatic transition but there is no evidence for pollinator preference for flower color, we hypothesize that floret color likely varies as a function of climatic conditions.

Variation in the Local Grey Mullet Populations (Mugil cephalus) on the Western Pacific Fringe

Background: Understanding population genetic structures is crucial for planning and implementing conservation programmes to preserve species' adaptive and evolutionary potential and thus ensure their long-term persistence. The grey mullet (Mugil cephalus) is a globally distributed coastal fish. Its populations in waters surrounding Taiwan on the western Pacific fringe are divided into at least two stocks (migratory and residential), but questions remain regarding their genetic divergence and gene flow. Methods and Results: To cast more light on this, allozyme variations at 21 presumptive gene loci of 1217 adult grey mullets from 15 localities in Japan, Taiwan and mainland China, and four gene loci from 1470 juveniles from three localities in Taiwan were used to investigate patterns of genetic variation. The mean expected heterozygosity (He) was 0.128-ranging from 0.031 (Matsu) to 0.442 (Kaoping)-and the mean observed heterozygosity (Ho) was 0.086-ranging from 0.017 (Kaohsiung) to 0.215 (Kaoping). Both AMOVA and the high overall mean FST of 0.252 indicated enormous genetic differentiation among populations and the positive mean value of FIS was 0.328, indicating a deficiency of heterozygotes. PCoA indicated that the samples of M. cephalus could be split into three groups and STRUCTURE analysis showed that all individuals were grouped into three genetic clusters. The results of mutation-drift equilibrium tests did not suggest that the populations experienced any recent genetic bottleneck. The results from all localities in the present investigation showed significant change in the GPI-A genotype frequencies with latitudes-e.g., increases in GPI-A*135/135 homozygote frequencies and GPI-A*100/100 frequencies were highly correlated with latitudinal cline. All migratory populations with the GPI-A genotype were almost exclusively the GPI-A*100/100 homozygote. During the life history of M. cephalus, the GPI-A*100/135 heterozygote frequency significantly decreases with age. Conclusions: Based on these data, we suggest that each GPI-A genotype represents trait combinations of higher fitness in some portions of the environment. Furthermore, the genotypic frequencies change in accordance with life stages, suggesting that selection occurs throughout the life span.

Mitonuclear and phenotypic discordance in an Atlantic Forest frog hybrid zone

Discordance between mitochondrial and nuclear DNA is common among animals and can be the result of a number of evolutionary processes, including incomplete lineage sorting and introgression. Particularly relevant in contact zones, mitonuclear discordance is expected because the mitochondrial genome is haploid and primarily uniparentally inherited, whereas nuclear loci are evolving at slower rates. In addition, when closely related taxa come together in hybrid zones, the distribution of diagnostic phenotypic characters and their concordance with the mitochondrial or nuclear lineages can also inform on historical and ongoing dynamics within hybrid zones. Overall, genetic and phenotypic discordances provide evidence for evolutionary divergence and processes that maintain boundaries among sister species or lineages. In this study, we characterized patterns of genetic and phenotypic variation in a contact zone between Cycloramphus dubius and Cycloramphus boraceiensis, two sister species of frogs endemic to the Atlantic Coastal Forest of Brazil. We examined genomic-scale nuclear diversification across 19 populations, encompassing the two parental forms and a contact zone between them. We compared the distribution of genomic DNA variability with that of a mitochondrial locus (16S) and two morphological traits (dorsal tubercles and body size). Our results reveal multiple divergent lineages with ongoing admixture. We detected discordance in patterns of introgression across the three data types. Cycloramphus dubius males are significantly larger than C. boraceiensis males, and we posit that competition among males in the hybrid zone, coupled with mate choice by females, may be one mechanism leading to patterns of introgression observed between the species.

Clinal variation in natural populations of Drosophila melanogaster: An old debate about natural selection and neutral processes

Distinguishing between environmental adaptations and neutral processes poses a challenge in population genetics and evolutionary studies, particularly when phenomena can be explained by both processes. Clines are genotypic or phenotypic characters correlated with environmental variables, because of that correlation, they are used as examples of spatially varying selection. At the same time, many genotypic clines can be explained by demographic history, like isolation by distance or secondary contact zones. Clines have been extensively studied in Drosophila melanogaster, especially in North America and Australia, where they are attributed to both differential selection and various demographic processes. This review explores existing literature supporting this conclusion and suggests new approaches to better understand the influence of these processes on clines. These innovative approaches aim to shed light on the longstanding debate regarding the importance of natural selection versus neutral processes in maintaining variation in natural populations.

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

Ecological niche contributes to the persistence of the western × glaucous-winged gull hybrid zone

Hybrid zones occur in nature when populations with limited reproductive barriers overlap in space. Many hybrid zones persist over time, and different models have been proposed to explain how selection can maintain hybrid zone stability. More empirical studies are needed to elucidate the role of ecological adaptation in maintaining stable hybrid zones. Here, we investigated the role of exogenous factors in maintaining a hybrid zone between western gulls (Larus occidentalis) and glaucous-winged gulls (L. glaucescens). We used ecological niche models (ENMs) and niche similarity tests to quantify and examine the ecological niches of western gulls, glaucous-winged gulls, and their hybrids. We found evidence of niche divergence between all three groups. Our results support the bounded superiority model, providing further evidence that exogenous selection favoring hybrids may be an important factor in maintaining this stable hybrid zone.

The Variability of Phytophthora infestans Isolates Collected from Estonian Islands in the Baltic Sea

Knowledge of a pathogen's genetic variability and population structure is of great importance to effective disease management. In this study, 193 isolates of Phytophthora infestans collected from three Estonian islands were characterized over 3 years using simple sequence repeat (SSR) marker data complemented by information on their mating type and resistance to metalaxyl. In combination with SSR marker data from samples in the neighboring Pskov region of Northwest Russia, the impact of regional and landscape structure on the level of genetic exchange was also examined. Among the 111 P. infestans isolates from Estonian islands, 49 alleles were detected among 12 SSR loci, and 59 SSR multilocus genotypes were found, of which 64% were unique. The genetic variation was higher among years than that among islands, as revealed by the analysis of molecular variance. The frequency of metalaxyl-resistant isolates increased from 9% in 2012 to 30% in 2014, and metalaxyl resistance was most frequent among A1 isolates. The test for isolation by distance among the studied regions was not significant, and coupled with the absence of genetic differentiation, the result revealed gene flow and the absence of local adaptation. The data are consistent with a sexual population in which diversity is driven by an annual germination of soilborne oospores. The absence of shared genotypes over the years has important implications when it comes to the management of diseases. Such population diversity can make it difficult to predict the nature of the outbreak in the coming year as the genetic makeup is different for each year.

Selection leads to remarkable variability in the outcomes of hybridisation across replicate hybrid zones

Hybrid zones have been viewed as an opportunity to see speciation in action. When hybrid zones are replicated, it is assumed that if the same genetic incompatibilities are maintaining reproductive isolation across all instances of secondary contact, those incompatibilities should be identifiable by consistent patterns in the genome. In contrast, changes in allele frequencies due to genetic drift should be idiosyncratic for each hybrid zone. To test this assumption, we simulated 20 replicates of each of 12 hybrid zone scenarios with varied genetic incompatibilities, rates of migration, selection and different starting population size ratios of parental species. We found remarkable variability in the outcomes of hybridisation in replicate hybrid zones, particularly with Bateson-Dobzhansky-Muller incompatibilities and strong selection. We found substantial differences among replicates in the overall genomic composition of individuals, including admixture proportions, inter-specific ancestry complement and number of ancestry junctions. Additionally, we found substantial variation in genomic clines among replicates at focal loci, regardless of locus-specific selection. We conclude that processes other than selection are responsible for some consistent outcomes of hybridisation, whereas selection on incompatibilities can lead to genomically widespread and highly variable outcomes. We highlight the challenge of mapping between pattern and process in hybrid zones and call attention to how selection against incompatibilities will commonly lead to variable outcomes. We hope that this study informs future research on replicate hybrid zones and encourages further development of statistical techniques, theoretical models and exploration of additional axes of variation to understand reproductive isolation.

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Temporal Variation in Selection Influences Microgeographic Local Adaptation

AbstractEcological heterogeneity can lead to local adaptation when populations exhibit fitness trade-offs among habitats. However, the degree to which local adaptation is affected by the spatial and temporal scale of environmental variation is poorly understood. A multiyear reciprocal transplant experiment was performed with populations of the annual plant Leptosiphon parviflorus living on adjacent serpentine and nonserpentine soil. Local adaptation over this small geographic scale was observed, but there were differences in the temporal variability of selection across habitats. On serpentine soil, the local population had a consistently large survival advantage, presumably as a result of the temporal stability in selection imposed by soil cation content. In contrast, a fecundity advantage was observed for the sandstone population on its native soil type but only in the two study years with the highest rainfall. A manipulative greenhouse experiment demonstrated that the fitness advantage of the sandstone population in its native soil type depends critically on water availability. The temporal variability in local adaptation driven by variation in precipitation suggests that continued drought conditions have the potential to erode local adaptation in these populations. These results show how different selective factors can influence spatial and temporal patterns of variation in fitness trade-offs.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Using Stable Isotopes to Determine Natal Origin and Feeding Habits of the Army Cutworm Moth, Euxoa auxiliaris (Lepidoptera: Noctuidae)

The army cutworm, Euxoa auxiliaris (Grote), is a migratory noctuid that is both an agricultural pest and an important late-season food source for grizzly bears, Ursus arctos horribilis (Linnaeus, Carnivora: Ursidae), within the Greater Yellowstone Ecosystem. Beyond the confirmation of the moths' seasonal, elevational migration in the mid-1900s, little else has been documented about their migratory patterns. To address this missing ecological component, we examined (1) migratory routes during their spring and fall migratory periods throughout their natal range, the Great Plains, and (2) natal origin at two of their summering ranges using stable hydrogen (δ2H) analyses of wings from samples collected within the areas of interest. Stable carbon (δ13C) and stable nitrogen (δ15N) analyses of wings were used to evaluate larval feeding habits of the migrants and agricultural intensity of natal origin sites, respectively. Results suggest that, rather than migrating exclusively east to west, army cutworm moths are also migrating north to south during their spring migration. Moths did not exhibit natal origin site fidelity when returning to the Great Plains. Migrants collected from the Absaroka Range had the highest probability of natal origin in Alberta, British Columbia, Saskatchewan, the most southern region of the Northwest Territories, and second highest probability of origin in Montana, Wyoming, and Idaho. Migrants collected in the Lewis Range had the highest probability of origin in the same provinces of Canada. Results suggest that migrants of the Absaroka Range fed exclusively on C3 plants as larvae and rarely fed in heavily fertilized agroecosystems.

Integrative characterization of genetic and phenotypic differentiation in an ant species complex with strong hierarchical population structure and low dispersal abilities

Low dispersal, occurrence of asexual reproduction and geographic discontinuity increase genetic differentiation between populations, which ultimately can lead to speciation. In this work, we used a multidisciplinary framework to characterize the genetic and phenotypic differentiation between and within two cryptic ant species with restricted dispersal, Cataglyphis cursor and C. piliscapa and used behavioral experiments to test for reproductive isolation. Their distribution is segregated by the Rhône River and they have been traditionally distinguished only by hair numbers, although a statistical assessment is still lacking. We found strong genetic (microsatellites, nuclear and mitochondrial sequences), morphological (number of hairs, tibia length, male genitalia) and chemical (cuticular hydrocarbons) differentiation not only between species but also among localities within species. However, inter-specific differentiation was slightly higher than intra-specific differentiation for most markers. Overall, this pattern could either reflect reproductive isolation or could result from a longer period of geographic isolation between species than among localities within species without necessarily involving reproductive isolation. Interestingly, our behavioral experiments showed an absence of mating between species associated to a higher aggressiveness of workers towards heterospecific males. This suggests that sexual selection may, at least partially, fuel reproductive isolation. We also showed that cuticular hydrocarbons, mtDNA sequences and number of hairs provide reliable criteria allowing species discrimination. Overall, this species complex offers a case study to further investigate varying stages of a speciation continuum by estimating reproductive isolation between pairs of localities varying by their level of genetic differentiation.

Analysis of genetic diversity and population structure of Babesia gibsoni

Babesia gibsoni is a tick-borne apicomplexan protozoan causing canine babesiosis. This parasite has diploid sexual reproduction in ticks, during which genetic exchanges can occur leading to increased genetic diversity, which is an important factor in adapting to environmental changes. Exploring the genetic variation of B. gibsoni population can provide a foundation for understanding the patterns of disease transmission and developing babesiosis control strategies. Partial 18S rRNA fragment sequences were obtained from 11 B. gibsoni isolates collected from different regions in China and 117 publicly available sequences were from 12 geographical areas including China. The genetic variation, demographic expansion and population structure were examined. A total of 34 haplotypes were identified among B. gibsoni populations. Analysis of molecular variance, pairwise Fst and structure analysis showed that high genetic variation within populations, low genetic differentiation and obvious mixture haplotype were apparent in a single continent, but higher genetic differentiation was detected across different continents. Neutrality tests implied that B. gibsoni populations had experienced population extension. These findings will contribute to understand the genetics and evolution of B. gibsoni and will be useful for formulating effective management strategies to prevent and control this parasite.

Selection on early survival does not explain germination rate clines in Mimulus cardinalis

PREMISE

Many traits covary with environmental gradients to form phenotypic clines. While local adaptation to the environment can generate phenotypic clines, other nonadaptive processes may also. If local adaptation causes phenotypic clines, then the direction of genotypic selection on traits should shift from one end of the cline to the other. Traditionally, genotypic selection on non-Gaussian traits like germination rate have been hampered because it is challenging to measure their genetic variance.

METHODS

Here we used quantitative genetics and reciprocal transplants to test whether a previously discovered cline in germination rate showed additional signatures of adaptation in the scarlet monkeyflower (Mimulus cardinalis). We measured genotypic and population level covariation between germination rate and early survival, a component of fitness. We developed a novel discrete log-normal model to estimate genetic variance in germination rate.

RESULTS

Contrary to our adaptive hypothesis, we found no evidence that genetic variation in germination rate contributed to variation in early survival. Across populations, southern populations in both gardens germinated earlier and survived more.

CONCLUSIONS

Southern populations have higher early survival but it is not caused by faster germination. This pattern is consistent with nonadaptive forces driving the phenotypic cline in germination rate, but future work will need to assess whether there is selection at other life stages. This statistical framework should help expand quantitative genetic analyses for other waiting-time traits.

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

Phenology and plasticity can prevent adaptive clines in thermal tolerance across temperate mountains: The importance of the elevation-time axis

Critical thermal limits (CTmax and CTmin) decrease with elevation, with greater change in CTmin, and the risk to suffer heat and cold stress increasing at the gradient ends. A central prediction is that populations will adapt to the prevailing climatic conditions. Yet, reliable support for such expectation is scant because of the complexity of integrating phenotypic, molecular divergence and organism exposure. We examined intraspecific variation of CTmax and CTmin, neutral variation for 11 microsatellite loci, and micro- and macro-temperatures in larvae from 11 populations of the Galician common frog (Rana parvipalmata) across an elevational gradient, to assess (1) the existence of local adaptation through a PST-FST comparison, (2) the acclimation scope in both thermal limits, and (3) the vulnerability to suffer acute heat and cold thermal stress, measured at both macro- and microclimatic scales. Our study revealed significant microgeographic variation in CTmax and CTmin, and unexpected elevation gradients in pond temperatures. However, variation in CTmax and CTmin could not be attributed to selection because critical thermal limits were not correlated to elevation or temperatures. Differences in breeding phenology among populations resulted in exposure to higher and more variable temperatures at mid and high elevations. Accordingly, mid- and high-elevation populations had higher CTmax and CTmin plasticities than lowland populations, but not more extreme CTmax and CTmin. Thus, our results support the prediction that plasticity and phenological shifts may hinder local adaptation, promoting thermal niche conservatism. This may simply be a consequence of a coupled variation of reproductive timing with elevation (the "elevation-time axis" for temperature variation). Mid and high mountain populations of R. parvipalmata are more vulnerable to heat and cool impacts than lowland populations during the aquatic phase. All of this contradicts some of the existing predictions on adaptive thermal clines and vulnerability to climate change in elevational gradients.

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26-taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.

Empirical and philosophical problems with the subspecies rank

Species-level taxonomy derives from empirical sources (data and techniques) that assess the existence of spatiotemporal evolutionary lineages via various species "concepts." These concepts determine if observed lineages are independent given a particular methodology and ontology, which relates the metaphysical species concept to what "kind" of thing a species is in reality. Often, species concepts fail to link epistemology back to ontology. This lack of coherence is in part responsible for the persistence of the subspecies rank, which in modern usage often functions as a placeholder between the evolutionary events of divergence or collapse of incipient species. Thus, prospective events like lineages merging or diverging require information from unknowable future information. This is also conditioned on evidence that the lineage already has a detectably distinct evolutionary history. Ranking these lineages as subspecies can seem attractive given that many lineages do not exhibit intrinsic reproductive isolation. We argue that using subspecies is indefensible on philosophical and empirical grounds. Ontologically, the rank of subspecies is either identical to that of species or undefined in the context of evolutionary lineages representing spatiotemporally defined individuals. Some species concepts more inclined to consider subspecies, like the Biological Species Concept, are disconnected from evolutionary ontology and do not consider genealogy. Even if ontology is ignored, methods addressing reproductive isolation are often indirect and fail to capture the range of scenarios linking gene flow to species identity over space and time. The use of subspecies and reliance on reproductive isolation as a basis for an operational species concept can also conflict with ethical issues governing the protection of species. We provide a way forward for recognizing and naming species that links theoretical and operational species concepts regardless of the magnitude of reproductive isolation.

A novel analytical framework to quantify co-gradient and countergradient variation

Spatial covariance between genotypic and environmental influences on phenotypes (Cov_GE ) can result in the nonrandom distribution of genotypes across environmental gradients and is a potentially important factor driving local adaptation. However, a framework to quantify the magnitude and significance of Cov_GE has been lacking. We develop a novel quantitative/analytical approach to estimate and test the significance of Cov_GE from reciprocal transplant or common garden experiments, which we validate using simulated data. We demonstrate how power to detect Cov_GE changes over a range of experimental designs. We confirm an inverse relationship between gene-by-environment interactions (GxE) and Cov_GE , as predicted by first principles, but show how phenotypes can be influenced by both. The metric provides a way to measure how phenotypic plasticity covaries with genetic differentiation and highlights the importance of understanding the dual influences of Cov_GE and GxE on phenotypes in studies of local adaptation and species' responses to environmental change.

Multivariate selection mediated by aridity predicts divergence of drought-resistant traits along natural aridity gradients of an invasive weed

Geographical variation in the environment underpins selection for local adaptation and evolutionary divergence among populations. Because many environmental conditions vary across species' ranges, identifying the specific environmental variables underlying local adaptation is profoundly challenging. We tested whether natural selection mediated by aridity predicts clinal divergence among invasive populations of capeweed (Arctotheca calendula) that established and spread across southern Australia during the last two centuries. Using common garden experiments with two environmental treatments (wet and dry) that mimic aridity conditions across capeweed's invasive range, we estimated clinal divergence and effects of aridity on fitness and multivariate phenotypic selection in populations sampled along aridity gradients in Australia. We show that: (1) capeweed populations have relatively high fitness in aridity environments similar to their sampling locations; (2) the magnitude and direction of selection strongly differs between wet and dry treatments, with drought stress increasing the strength of selection; and (3) differences in directional selection between wet and dry treatments predict patterns of clinal divergence across the aridity gradient, particularly for traits affecting biomass, flowering phenology and putative antioxidant expression. Our results suggest that aridity-mediated selection contributes to trait diversification among invasive capeweed populations, possibly facilitating the expansion of capeweed across southern Australia.

Sexual dimorphism in phenotypic plasticity and persistence under environmental change: An extension of theory and meta-analysis of current data

Populations must adapt to environmental changes to remain viable. Both evolution and phenotypic plasticity contribute to adaptation, with plasticity possibly being more important for coping with rapid change. Adaptation is complex in species with separate sexes, as the sexes can differ in the strength or direction of natural selection, the genetic basis of trait variation, and phenotypic plasticity. Many species show sex differences in plasticity, yet how these differences influence extinction susceptibility remains unclear. We first extend theoretical models of population persistence in changing environments and show that persistence is affected by sexual dimorphism for phenotypic plasticity, trait genetic architecture, and sex-specific selection. Our models predict that female-biased adaptive plasticity-particularly in traits with modest-to-low cross-sex genetic correlations-typically promotes persistence, though we also identify conditions where sexually monomorphic or male-biased plasticity promotes persistence. We then perform a meta-analysis of sex-specific plasticity under manipulated thermal conditions. Although examples of sexually dimorphic plasticity are widely observed, systematic sex differences are rare. An exception-cold resistance-is systematically female-biased and represents a trait wherein sexually dimorphic plasticity might elevate population viability in changing environments. We discuss our results in light of debates about the roles of evolution and plasticity in extinction susceptibility.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

The evolution of thermal performance in native and invasive populations of Mimulus guttatus

The rise of globalization has spread organisms beyond their natural range, allowing further opportunity for species to adapt to novel environments and potentially become invaders. Yet, the role of thermal niche evolution in promoting the success of invasive species remains poorly understood. Here, we use thermal performance curves (TPCs) to test hypotheses about thermal adaptation during the invasion process. First, we tested the hypothesis that if species largely conserve their thermal niche in the introduced range, invasive populations may not evolve distinct TPCs relative to native populations, against the alternative hypothesis that thermal niche and therefore TPC evolution has occurred in the invasive range. Second, we tested the hypothesis that clines of TPC parameters are shallower or absent in the invasive range, against the alternative hypothesis that with sufficient time, standing genetic variation, and temperature‐mediated selection, invasive populations would re‐establish clines found in the native range in response to temperature gradients. To test these hypotheses, we built TPCs for 18 native (United States) and 13 invasive (United Kingdom) populations of the yellow monkeyflower, Mimulus guttatus. We grew clones of multiple genotypes per population at six temperature regimes in growth chambers. We found that invasive populations have not evolved different thermal optima or performance breadths, providing evidence for evolutionary stasis of thermal performance between the native and invasive ranges after over 200 years post introduction. Thermal optimum increased with mean annual temperature in the native range, indicating some adaptive differentiation among native populations that was absent in the invasive range. Further, native and invasive populations did not exhibit adaptive clines in thermal performance breadth with latitude or temperature seasonality. These findings suggest that TPCs remained unaltered post invasion, and that invasion may proceed via broad thermal tolerance and establishment in already climatically suitable areas rather than rapid evolution upon introduction.

Species delimitation and mitonuclear discordance within a species complex of biting midges

The inability to distinguish between species can be a serious problem in groups responsible for pathogen transmission. Culicoides biting midges transmit many pathogenic agents infecting wildlife and livestock. In North America, the C. variipennis species complex contains three currently recognized species, only one of which is a known vector, but limited species-specific characters have hindered vector surveillance. Here, genomic data were used to investigate population structure and genetic differentiation within this species complex. Single nucleotide polymorphism data were generated for 206 individuals originating from 17 locations throughout the United States and Canada. Clustering analyses suggest the occurrence of two additional cryptic species within this complex. All five species were significantly differentiated in both sympatry and allopatry. Evidence of hybridization was detected in three different species pairings indicating incomplete reproductive isolation. Additionally, COI sequences were used to identify the hybrid parentage of these individuals, which illuminated discordance between the divergence of the mitochondrial and nuclear datasets.

Genetic traces of hybrid zone movement across a fragmented habitat

Theoretical and empirical studies suggest that the structure and position of hybrid zones can change over time. Evidence for moving hybrid zones has been directly inferred by repeated sampling over time, or indirectly through the detection of genetic footprints left by the receding species and the resulting asymmetric patterns of introgression across markers. We here investigate a hybrid zone formed by two subspecies of the Iberian golden-striped salamander, Chioglossa lusitanica, using a panel of 35 nuclear loci (31 SNPs and 4 allozymes) and one mitochondrial locus in a transect in central Portugal. We found concordant and coincident clines for most of the nuclear loci (n=22, 63%), defining a narrow hybrid zone of ca. 6 km wide, with the centre positioned ca. 15 km south of the Mondego river. Asymmetric introgression was observed at another 14 loci. Their clines are displaced towards the north, with positions located either close to the Mondego river (n=6), or further northwards (n=8). We interpret these profiles as genetic traces of the southward displacement of C. lusitanica lusitanica by C. l. longipes over the wider Mondego river valley. We noted the absence of significant linkage disequilibrium and we inferred low levels of effective selection per locus against hybrids, suggesting that introgression in the area of species replacement occurred under a neutral diffusion process. A species distribution model suggests that the C. lusitanica hybrid zone coincides with a narrow corridor of fragmented habitat. From the position of the displaced clines, we infer that patches of locally suitable habitat trapped some genetic variants that became disassociated from the southward moving hybrid zone. This study highlights the influence of habitat availability on hybrid zone movement.

The Role of Directed Dispersal in Driving Genetic and Morphological Structure in Invasive Smallmouth Bass

Dispersal is an essential life-history trait crucial to species persistence and diversification. This is particularly important in spatiotemporal fluctuating environments such as freshwater habitats, where species movement is confined to the dendritic network and wetted boundaries. To persist in such fluctuating environments, a species can modify, adaptively and plastically, its phenotypic variation to better match the environment or escape via directed dispersal to a more suitable habitat (i.e., matching habitat choice). We use the invasive smallmouth bass, Micropterus dolomieu, sampled at 10 km intervals, to assess the effect of directed dispersal on the fine scale genetic and phenotypic variation in populations of M. dolomieu along a river course. Gene flow was used as a proxy for dispersal. By comparing population genetic structure, morphological variation (of linear traits and geometric landmarks), and environmental heterogeneity, we discovered a clear correlation between environmental variation and morphological traits. Although isolation by distance seemed to have shaped the overall genetic pattern detected among the populations, the strong genetic structuring observed within the Ratel tributary appeared to be non-random. These results provide novel insights into the potential mechanisms promoting the spread and establishment of invasive species and the possible influence multiple introductions may have on fine scale genetic structuring.

Plastic and quantitative genetic divergence mirror environmental gradients among wild, fragmented populations of Impatiens capensis

PREMISE

Habitat fragmentation generates molecular genetic divergence among isolated populations, but few studies have assessed phenotypic divergence and fitness in populations where the genetic consequences of habitat fragmentation are known. Phenotypic divergence could reflect plasticity, local adaptation, and/or genetic drift.

METHODS

We examined patterns and potential drivers of phenotypic divergence among 12 populations of jewelweed (Impatiens capensis) that show strong molecular genetic signals of isolation and drift among fragmented habitats. We measured morphological and reproductive traits in both maternal plants within natural populations and their self-fertilized progeny grown together in a common garden. We also quantified environmental divergence between home sites and the common garden.

RESULTS

Populations with less molecular genetic variation expressed less maternal phenotypic variation. Progeny in the common garden converged in phenotypes relative to their wild mothers but retained among-population differences in morphology, survival, and reproduction. Among-population phenotypic variance was 3-10× greater in home sites than in the common garden for 6 of 7 morphological traits measured. Patterns of phenotypic divergence paralleled environmental gradients in ways suggestive of adaptation. Progeny resembled their mothers less as the environmental distance between their home site and the common garden increased.

CONCLUSIONS

Despite strong molecular signatures of isolation and drift, phenotypic differences among these Impatiens populations appear to reflect both adaptive quantitative genetic divergence and plasticity. Quantifying the extent of local adaptation and plasticity and how these covary with molecular and phenotypic variation help us predict when populations may lose their adaptive capacity.

Spatial population structure: patterns of adaptation in populations of the water hyacinth grasshopper Cornops aquaticum (Bruner 1906)

The water hyacinth grasshopper Cornops aquaticum (Bruner 1906) constitutes an appropriate model to assess phenotypic and karyotypic variability in the context of differentiation and adaptation of insect populations because it occurs over a wide latitudinal range. This study represents a general analysis of phenotype, karyotype and molecular variation in native populations of C. aquaticum in South America. This is also relevant because this insect is considered a promising biological control agent of water hyacinth, a native South American aquatic plant but a pest in South Africa. Along Paraná and Uruguay River Basins, body size correlated negatively with latitude, and positively so with temperature and rainfall in both sexes. To test whether the chromosomal and phenotypic patterns were adaptive, we compared them with neutral microsatellite loci variation in populations from the medium and lower course of the Paraná River. Firstly, the lack of pairwise association between karyotype and phenotype distance matrixes with that of neutral loci suggested non-neutrality. Secondly, phenotypic differentiation for all morphometric traits (PST) was significantly larger than molecular differentiation (FST), indicating a prevailing divergence selection effect on the observed phenotypic patterns. Finally, the phenotypic and genotypic spatial structures - inferred from Bayesian approaches - were discordant: neutral genetic structure clustered together most populations except for the two southernmost, downstream ones, whereas phenotypic spatial structure groups together all the deltaic populations and singles out the two northernmost ones. The results suggest directional selection leading to higher centric fusion frequencies in the downstream populations and favouring morphometric optimal differences in relation to the environment.

Mito-nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern "Rocky Mountain'' color form with ferruginous mid-abdominal segments (B. m. melanopygus) and a southern "Pacific'' form with black mid-abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color-controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito-nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.

Parallel morphological evolution and habitat-dependent sexual dimorphism in cave- vs. surface populations of the Asellus aquaticus (Crustacea: Isopoda: Asellidae) species complex

Studying parallel evolution (repeated, independent evolution of similar phenotypes in similar environments) is a powerful tool to understand environment-dependent selective forces. Surface-dwelling species that repeatedly and independently colonized caves provide unique models for such studies. The primarily surface-dwelling Asellus aquaticus species complex is a good candidate to carry out such research, because it colonized several caves in Europe. By comparing 17 functional morphological traits between six cave and nine surface populations of the A. aquaticus species complex, we investigated population divergence in morphology and sexual dimorphism. We found habitat-dependent population divergence in 10 out of 17 traits, likely reflecting habitat-driven changes in selection acting on sensory systems, feeding, grooming, and antipredator mechanisms. Sexual dimorphism was present in 15 traits, explained by sexual selection acting on male traits important in male-male agonistic behavior or mate guarding and fecundity selection acting on female traits affecting offspring number and nursing. In eight traits, the degree of sexual dimorphism was habitat dependent. We conclude that cave-related morphological changes are highly trait- and function-specific and that the strength of sexual/fecundity selection strongly differs between cave and surface habitats. The considerable population variation within habitat type warrants further studies to reveal cave-specific adaptations besides the parallel patterns.

Assortative mating and epistatic mating-trait architecture induce complex movement of the crow hybrid zone

Hybrid zones provide a window into the evolutionary processes governing species divergence. Yet, the contribution of mate choice to the temporal and spatial stability of hybrid zones remains poorly explored. Here, we investigate the effects of assortative mating on hybrid-zone dynamics by means of a mathematical model parameterized with phenotype and genotype data from the hybrid zone between all-black carrion and gray-coated hooded crows. In the best-fit model, narrow clines of the two mating-trait loci were maintained by a moderate degree of assortative mating inducing pre- and postzygotic isolation via positive frequency-dependent selection. Epistasis between the two loci induced hybrid-zone movement in favor of alleles conveying dark plumage followed by a shift in the opposite direction favoring gray-coated phenotypes ∼ 1 200 generations after secondary contact. Unlinked neutral loci diffused near-unimpeded across the zone. These results were generally robust to the choice of matching rule (self-referencing or parental imprinting) and effects of genetic drift. Overall, this study illustrates under which conditions assortative mating can maintain steep clines in mating-trait loci without generalizing to genome-wide reproductive isolation. It further emphasizes the importance of the genetic mating-trait architecture for spatio-temporal hybrid-zone dynamics.

On the Trail of Spatial Patterns of Genetic Variation

The accurate determination of the spatial trends on the variability of a species’ gene pool is essential to elucidate the underlying demographic-evolutionary events, thus helping to unravel the microevolutionary history of the population under study. Herein we present a new software called GenoCline, mainly addressed to detect genetic clines from allele, haplotype, and genome-wide data. This program package allows identifying the geographic orientation of clinal genetic variation through a system of iterative rotation of a virtual coordinate axis. Besides, GenoCline can perform complementary analyses to explore the potential origin of the genetic clines observed, including spatial autocorrelation, isolation by distance, centroid method, multidimensional scaling and Sammon projection. Among the advantages of this software is the ease in data entry and potential interconnection with other programs. Genetic and geographic data can be entered in spreadsheet table formatting (.xls), whereas genome-wide data can be imported in Eigensoft format. Genetic frequencies can also be exported in a format compatible with other programs dealing with population genetic and evolutionary biology analyses. All illustrations of results are saved in.svg format so that there will be high quality and easily editable vectorial graphs available for the researcher. Being implemented in Java, GenoCline is highly portable, thus working in different operating systems.

Thermal tolerance and preference are both consistent with the clinal distribution of house fly proto-Y chromosomes

Selection pressures can vary within localized areas and across massive geographical scales. Temperature is one of the best studied ecologically variable abiotic factors that can affect selection pressures across multiple spatial scales. Organisms rely on physiological (thermal tolerance) and behavioral (thermal preference) mechanisms to thermoregulate in response to environmental temperature. In addition, spatial heterogeneity in temperatures can select for local adaptation in thermal tolerance, thermal preference, or both. However, the concordance between thermal tolerance and preference across genotypes and sexes within species and across populations is greatly understudied. The house fly, Musca domestica, is a well-suited system to examine how genotype and environment interact to affect thermal tolerance and preference. Across multiple continents, house fly males from higher latitudes tend to carry the male-determining gene on the Y chromosome, whereas those from lower latitudes usually have the male determiner on the third chromosome. We tested whether these two male-determining chromosomes differentially affect thermal tolerance and preference as predicted by their geographical distributions. We identify effects of genotype and developmental temperature on male thermal tolerance and preference that are concordant with the natural distributions of the chromosomes, suggesting that temperature variation across the species range contributes to the maintenance of the polymorphism. In contrast, female thermal preference is bimodal and largely independent of congener male genotypes. These sexually dimorphic thermal preferences suggest that temperature-dependent mating dynamics within populations could further affect the distribution of the two chromosomes. Together, the differences in thermal tolerance and preference across sexes and male genotypes suggest that different selection pressures may affect the frequencies of the male-determining chromosomes across different spatial scales.

Ecological basis and genetic architecture of crypsis polymorphism in the desert clicker grasshopper (Ligurotettix coquilletti)

Color polymorphic species can offer exceptional insight into the ecology and genetics of adaptation. Although the genetic architecture of animal coloration is diverse, many color polymorphisms are associated with large structural variants and maintained by biotic interactions. Grasshoppers are notably polymorphic in both color and karyotype, which makes them excellent models for understanding the ecological drivers and genetic underpinnings of color variation. Banded and uniform morphs of the desert clicker grasshopper (Ligurotettix coquilletti) are found across the western deserts of North America. To address the hypothesis that predation maintains local color polymorphism and shapes regional crypsis variation, we surveyed morph frequencies and tested for covariation with two predation environments. Morphs coexisted at intermediate frequencies at most sites, consistent with local balancing selection. Morph frequencies covaried with the appearance of desert substrate-an environment used only by females-suggesting that ground-foraging predators are major agents of selection on crypsis. We next addressed the hypothesized link between morph variation and genome structure. To do so, we designed an approach for detecting inversions and indels using only RADseq data. The banded morph was perfectly correlated with a large putative indel. Remarkably, indel dominance differed among populations, a rare example of dominance evolution in nature.

The evolutionary genomics of species' responses to climate change

Climate change is a threat to biodiversity. One way that this threat manifests is through pronounced shifts in the geographical range of species over time. To predict these shifts, researchers have primarily used species distribution models. However, these models are based on assumptions of niche conservatism and do not consider evolutionary processes, potentially limiting their accuracy and value. To incorporate evolution into the prediction of species' responses to climate change, researchers have turned to landscape genomic data and examined information about local genetic adaptation using climate models. Although this is an important advancement, this approach currently does not include other evolutionary processes-such as gene flow, population dispersal and genomic load-that are critical for predicting the fate of species across the landscape. Here, we briefly review the current practices for the use of species distribution models and for incorporating local adaptation. We next discuss the rationale and theory for considering additional processes, reviewing how they can be incorporated into studies of species' responses to climate change. We summarize with a conceptual framework of how manifold layers of information can be combined to predict the potential response of specific populations to climate change. We illustrate all of the topics using an exemplar dataset and provide the source code as potential tutorials. This Perspective is intended to be a step towards a more comprehensive integration of population genomics with climate change science.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

Cold adaptation drives population genomic divergence in the ecological specialist, Drosophila montana

Detecting signatures of ecological adaptation in comparative genomics is challenging, but analysing population samples with characterised geographic distributions, such as clinal variation, can help identify genes showing covariation with important ecological variation. Here, we analysed patterns of geographic variation in the cold-adapted species Drosophila montana across phenotypes, genotypes and environmental conditions and tested for signatures of cold adaptation in population genomic divergence. We first derived the climatic variables associated with the geographic distribution of 24 populations across two continents to trace the scale of environmental variation experienced by the species, and measured variation in the cold tolerance of the flies of six populations from different geographic contexts. We then performed pooled whole genome sequencing of these six populations, and used Bayesian methods to identify SNPs where genetic differentiation is associated with both climatic variables and the population phenotypic measurements, while controlling for effects of demography and population structure. The top candidate SNPs were enriched on the X and fourth chromosomes, and they also lay near genes implicated in other studies of cold tolerance and population divergence in this species and its close relatives. We conclude that ecological adaptation has contributed to the divergence of D. montana populations throughout the genome and in particular on the X and fourth chromosomes, which also showed highest interpopulation F_ST . This study demonstrates that ecological selection can drive genomic divergence at different scales, from candidate genes to chromosome-wide effects.

Spatial population genetics in heavily managed species: Separating patterns of historical translocation from contemporary gene flow in white-tailed deer

Approximately 100 years ago, unregulated harvest nearly eliminated white-tailed deer (Odocoileus virginianus) from eastern North America, which subsequently served to catalyze wildlife management as a national priority. An extensive stock-replenishment effort soon followed, with deer broadly translocated among states as a means of re-establishment. However, an unintended consequence was that natural patterns of gene flow became obscured and pretranslocation signatures of population structure were replaced. We applied cutting-edge molecular and biogeographic tools to disentangle genetic signatures of historical management from those reflecting spatially heterogeneous dispersal by evaluating 35,099 single nucleotide polymorphisms (SNPs) derived via reduced-representation genomic sequencing from 1143 deer sampled statewide in Arkansas. We then employed Simpson's diversity index to summarize ancestry assignments and visualize spatial genetic transitions. Using sub-sampled transects across these transitions, we tested clinal patterns across loci against theoretical expectations of their response under scenarios of re-colonization and restricted dispersal. Two salient results emerged: (A) Genetic signatures from historic translocations are demonstrably apparent; and (B) Geographic filters (major rivers; urban centers; highways) now act as inflection points for the distribution of this contemporary ancestry. These results yielded a statewide assessment of contemporary population structure in deer as driven by historic translocations as well as ongoing processes. In addition, the analytical framework employed herein to effectively decipher extant/historic drivers of deer distribution in Arkansas is also applicable for other biodiversity elements with similarly complex demographic histories.

Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae)

BACKGROUND

Elucidating the effects of geography and selection on genetic variation is critical for understanding the relative importance of adaptation in driving differentiation and identifying the environmental factors underlying its occurrence. Adaptive genetic variation is common in tree species, especially widely distributed long-lived species. Pseudotaxus chienii can occupy diverse habitats with environmental heterogeneity and thus provides an ideal material for investigating the process of population adaptive evolution. Here, we characterize genetic and expression variation patterns and investigate adaptive genetic variation in P. chienii populations.

RESULTS

We generated population transcriptome data and identified 13,545 single nucleotide polymorphisms (SNPs) in 5037 unigenes across 108 individuals from 10 populations. We observed lower nucleotide diversity (π = 0.000701) among the 10 populations than observed in other gymnosperms. Significant negative correlations between expression diversity and nucleotide diversity in eight populations suggest that when the species adapts to the surrounding environment, gene expression and nucleotide diversity have a reciprocal relationship. Genetic structure analyses indicated that each distribution region contains a distinct genetic group, with high genetic differentiation among them due to geographical isolation and local adaptation. We used FST outlier, redundancy analysis, and latent factor mixed model methods to detect molecular signatures of local adaptation. We identified 244 associations between 164 outlier SNPs and 17 environmental variables. The mean temperature of the coldest quarter, soil Fe and Cu contents, precipitation of the driest month, and altitude were identified as the most important determinants of adaptive genetic variation. Most candidate unigenes with outlier signatures were related to abiotic and biotic stress responses, and the monoterpenoid biosynthesis and ubiquitin-mediated proteolysis KEGG pathways were significantly enriched in certain populations and deserve further attention in other long-lived trees.

CONCLUSIONS

Despite the strong population structure in P. chienii, genomic data revealed signatures of divergent selection associated with environmental variables. Our research provides SNPs, candidate unigenes, and biological pathways related to environmental variables to facilitate elucidation of the genetic variation in P. chienii in relation to environmental adaptation. Our study provides a promising tool for population genomic analyses and insights into the molecular basis of local adaptation.

Localized environmental heterogeneity drives the population differentiation of two endangered and endemic Opisthopappus Shih species

BACKGROUND

Climate heterogeneity not only indirectly shapes the genetic structures of plant populations, but also drives adaptive divergence by impacting demographic dynamics. The variable localized climates and topographic complexity of the Taihang Mountains make them a major natural boundary in Northern China that influences the divergence of organisms distributed across this region. Opisthopappus is an endemic genus of the Taihang Mountains that includes only two spatially partitioned species Opisthopappus longilobus and Opisthopappus taihangensis. For this study, the mechanisms behind the genetic variations in Opisthopappus populations were investigated.

RESULTS

Using SNP and InDel data coupled with geographic and climatic information, significant genetic differentiation was found to exist either between Opisthopappus populations or two species. All studied populations were divided into two genetic groups with the differentiation of haplotypes between the groups. At approximately 17.44 Ma of the early Miocene, O. taihangensis differentiated from O. longilobus under differing precipitation regimes due to the intensification of the Asian monsoon. Subsequently, intraspecific divergence might be induced by the dramatic climatic transformation from the mid- to late Miocene. During the Pleistocene period, the rapid uplift of the Taihang Mountains coupled with violent climatic oscillations would further promote the diversity of the two species. Following the development of the Taihang Mountains, its complex topography created geographical and ecological heterogeneity, which could lead to spatiotemporal isolation between the Opisthopappus populations. Thus the adaptive divergence might occur within these intraspecific populations in the localized heterogeneous environment of the Taihang Mountains.

CONCLUSIONS

The localized environmental events through the integration of small-scale spatial effects impacted the demographic history and differentiation mechanism of Opisthopappus species in the Taihang Mountains. The results provide useful information for us to understand the ecology and evolution of organisms in the mountainous environment from population and species perspective.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

A population-genomic approach for estimating selection on polygenic traits in heterogeneous environments

Strong selection can cause rapid evolutionary change, but temporal fluctuations in the form, direction and intensity of selection can limit net evolutionary change over longer time periods. Fluctuating selection could affect molecular diversity levels and the evolution of plasticity and ecological specialization. Nonetheless, this phenomenon remains understudied, in part because of analytical limitations and the general difficulty of detecting selection that does not occur in a consistent manner. Herein, I fill this analytical gap by presenting an approximate Bayesian computation (ABC) method to detect and quantify fluctuating selection on polygenic traits from population genomic time-series data. I propose a model for environment-dependent phenotypic selection. The evolutionary genetic consequences of selection are then modelled based on a genotype-phenotype map. Using simulations, I show that the proposed method generates accurate and precise estimates of selection when the generative model for the data is similar to the model assumed by the method. The performance of the method when applied to an evolve-and-resequence study of host adaptation in the cowpea seed beetle (Callosobruchus maculatus) was more idiosyncratic and depended on specific analytical choices. Despite some limitations, these results suggest the proposed method provides a powerful approach to connect the causes of (variable) selection to traits and genome-wide patterns of evolution. Documentation and open-source computer software (fsabc) implementing this method are available from github (https://github.com/zgompert/fsabc.git).

Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.