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

Genomic variation in the American pika: signatures of geographic isolation and implications for conservation

Background

Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada.

Results

Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θ_W = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θ_W = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θ_W = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide.

Conclusions

Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

Double-digest RAD-sequencing: do pre- and post-sequencing protocol parameters impact biological results?

Next-generation sequencing technologies have opened a new era of research in population genetics. Following these new sequencing opportunities, the use of restriction enzyme-based genotyping techniques, such as restriction site-associated DNA sequencing (RAD-seq) or double-digest RAD-sequencing (ddRAD-seq), has dramatically increased in the last decade. From DNA sampling to SNP calling, the laboratory and bioinformatic parameters of enzyme-based techniques have been investigated in the literature. However, the impact of those parameters on downstream analyses and biological results remains less documented. In this study, we investigated the effects of sevral pre- and post-sequencing settings on ddRAD-seq results for two biological systems: a complex of butterfly species (Coenonympha sp.) and several populations of common beech (Fagus sylvatica). Our results suggest that pre-sequencing parameters (i.e., DNA quantity, number of PCR cycles during library preparation) have a significant impact on the number of recovered reads and SNPs, on the number of unique alleles and on individual heterozygosity. In the same way, we found that post-sequencing settings (i.e., clustering and minimum coverage thresholds) influenced loci reconstruction (e.g., number of loci, mean coverage) and SNP calling (e.g., number of SNPs; heterozygosity) but had only a marginal impact on downstream analyses (e.g., measure of genetic differentiation, estimation of individual admixture, and demographic inferences). In addition, replication analyses confirmed the reproducibility of the ddRAD-seq procedure. Overall, this study assesses the degree of sensitivity of ddRAD-seq data to pre- and post-sequencing protocols, and illustrates its robustness when studying population genetics.

Harnessing landscape genomics to identify future climate resilient genotypes in a desert annual

Local adaptation features critically in shaping species responses to changing environments, complicating efforts to revegetate degraded areas. Rapid climate change poses an additional challenge that could reduce fitness of even locally sourced seeds in restoration. Predictive restoration strategies that apply seeds with favourable adaptations to future climate may promote long-term resilience. Landscape genomics is increasingly used to assess spatial patterns in local adaption and may represent a cost-efficient approach for identifying future-adapted genotypes. To demonstrate such an approach, we genotyped 760 plants from 64 Mojave Desert populations of the desert annual Plantago ovata. Genome scans on 5,960 SNPs identified 184 potentially adaptive loci related to climate and satellite vegetation metrics. Causal modelling indicated that variation in potentially adaptive loci was not confounded by isolation by distance or isolation by habitat resistance. A generalized dissimilarity model (GDM) attributed spatial turnover in potentially adaptive loci to temperature, precipitation and NDVI amplitude, a measure of vegetation green-up potential. By integrating a species distribution model (SDM), we find evidence that summer maximum temperature may both constrain the range of P. ovata and drive adaptive divergence in populations exposed to higher temperatures. Within the species' current range, warm-adapted genotypes are predicted to experience a fivefold expansion in climate niche by midcentury and could harbour key adaptations to cope with future climate. We recommend eight seed transfer zones and project each zone into its relative position in future climate. Prioritizing seed collection efforts on genotypes with expanding future habitat represents a promising strategy for restoration practitioners to address rapidly changing climates.

Evolutionary responses of marine organisms to urbanized seascapes

Many of the world's major cities are located in coastal zones, resulting in urban and industrial impacts on adjacent marine ecosystems. These pressures, which include pollutants, sewage, runoff and debris, temperature increases, hardened shorelines/structures, and light and acoustic pollution, have resulted in new evolutionary landscapes for coastal marine organisms. Marine environmental changes influenced by urbanization may create new selective regimes or may influence neutral evolution via impacts on gene flow or partitioning of genetic diversity across seascapes. While some urban selective pressures, such as hardened surfaces, are similar to those experienced by terrestrial species, others, such as oxidative stress, are specific to aquatic environments. Moreover, spatial and temporal scales of evolutionary responses may differ in the ocean due to the spatial extent of selective pressures and greater capacity for dispersal/gene flow. Here, we present a conceptual framework and synthesis of current research on evolutionary responses of marine organisms to urban pressures. We review urban impacts on genetic diversity and gene flow and examine evidence that marine species are adapting, or are predicted to adapt, to urbanization over rapid evolutionary time frames. Our findings indicate that in the majority of studies, urban stressors are correlated with reduced genetic diversity. Genetic structure is often increased in urbanized settings, but artificial structures can also act as stepping stones for some hard-surface specialists, promoting range expansion. Most evidence for rapid adaptation to urban stressors comes from studies of heritable tolerance to pollutants in a relatively small number of species; however, the majority of marine ecotoxicology studies do not test directly for heritability. Finally, we highlight current gaps in our understanding of evolutionary processes in marine urban environments and present a framework for future research to address these gaps.

Climate change driven hyposalinity as a selective agent in the littoral mesoherbivore Idotea balthica

Climate change will include a decrease in seawater salinity in the Baltic Sea. We quantified the effects of the projected future desalination on survival of the early life stage of the littoral herbivore Idotea balthica. We collected egg-bearing Idotea from three range-margin Baltic Sea populations, we exposed half of each brood to either current (6‰) or future salinity (3.5‰). We genotyped a subsample of each brood to analyse patterns of allelic change and to identify genomic regions targeted by selection. The survival was overall reduced by hyposalinity and broods varied in response to hyposalinity implying genetic variation in tolerance, with a stronger decrease in genetic diversity in future salinity. Finally, we identified proteins with crucial roles in basic cellular functions. This study indicates that projected future northern Baltic Sea hyposalinity will not just hamper I. balthica survival, but its selective pressure may also affect genetic diversity and cell physiology.

Adaptive Evolution in Cities: Progress and Misconceptions

Current narratives suggest that urban adaptation - the adaptive evolution of organisms to cities - is pervasive across taxa and cities. However, in reviewing hundreds of studies, we find only six comprehensive examples of species adaptively evolving to urbanization. We discuss the utility and shortcomings of methods for studying urban adaptation. We then review diverse systems offering preliminary evidence for urban adaptation and outline a research program for advancing its study. Urban environments constitute diverse, interacting selective agents that test the limits of adaptation. Understanding urban adaptation therefore offers unique opportunities for addressing fundamental questions in evolutionary biology and for better conserving biodiversity in cities. However, capitalizing on these opportunities requires appropriate research methods and dissemination of accurate narratives.

Landscape genomics of Quercus lobata reveals genes involved in local climate adaptation at multiple spatial scales

Understanding how the environment shapes genetic variation provides critical insight about the evolution of local adaptation in natural populations. At multiple spatial scales and multiple geographic contexts within a single species, such information could address a number of fundamental questions about the scale of local adaptation and whether or not the same loci are involved at different spatial scales or geographic contexts. We used landscape genomic approaches from three local elevational transects and rangewide sampling to (a) identify genetic variation underlying local adaptation to environmental gradients in the California endemic oak, Quercus lobata; (b) examine whether putatively adaptive SNPs show signatures of selection at multiple spatial scales; and (c) map putatively adaptive variation to assess the scale and pattern of local adaptation. Of over 10 k single-nucleotide polymorphisms (SNPs) generated with genotyping-by-sequencing, we found signatures of natural selection by climate or local environment at over 600 SNPs (536 loci), some at multiple spatial scales across multiple analyses. Candidate SNPs identified with gene-environment tests (LFMM) at the rangewide scale also showed elevated associations with climate variables compared to the background at both rangewide and elevational transect scales with gradient forest analysis. Some loci overlap with those detected in other oak species, raising the question of whether the same loci might be involved in local climate adaptation in different congeneric species that inhabit different geographic contexts. Mapping landscape patterns of adaptive versus background genetic variation identified regions of marked local adaptation and suggests nonlinear association of candidate SNPs and environmental variables. Taken together, our results offer robust evidence for novel candidate genes for local climate adaptation at multiple spatial scales.

Whole genome resequencing reveals genomic regions associated with thermal adaptation in redband trout

Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole-genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri, a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency-based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%-75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests cerk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.

A comparative genomics multitool for scientific discovery and conservation

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity.

A whole-genome alignment of 240 phylogenetically diverse species of eutherian mammal—including 131 previously uncharacterized species—from the Zoonomia Project provides data that support biological discovery, medical research and conservation.

Rapture-ready darters: Choice of reference genome and genotyping method (whole-genome or sequence capture) influence population genomic inference in Etheostoma

Researchers studying nonmodel organisms have an increasing number of methods available for generating genomic data. However, the applicability of different methods across species, as well as the effect of reference genome choice on population genomic inference, remain difficult to predict in many cases. We evaluated the impact of data type (whole-genome vs. reduced representation) and reference genome choice on data quality and on population genomic and phylogenomic inference across several species of darters (subfamily Etheostomatinae), a highly diverse radiation of freshwater fish. We generated a high-quality reference genome and developed a hybrid RADseq/sequence capture (Rapture) protocol for the Arkansas darter (Etheostoma cragini). Rapture data from 1,900 individuals spanning four darter species showed recovery of most loci across darter species at high depth and consistent estimates of heterozygosity regardless of reference genome choice. Loci with baits spanning both sides of the restriction enzyme cut site performed especially well across species. For low-coverage whole-genome data, choice of reference genome affected read depth and inferred heterozygosity. For similar amounts of sequence data, Rapture performed better at identifying fine-scale genetic structure compared to whole-genome sequencing. Rapture loci also recovered an accurate phylogeny for the study species and demonstrated high phylogenetic informativeness across the evolutionary history of the genus Etheostoma. Low cost and high cross-species effectiveness regardless of reference genome suggest that Rapture and similar sequence capture methods may be worthwhile choices for studies of diverse species radiations.

Genomic Prediction of (Mal)Adaptation Across Current and Future Climatic Landscapes

Signals of local adaptation have been found in many plants and animals, highlighting the heterogeneity in the distribution of adaptive genetic variation throughout species ranges. In the coming decades, global climate change is expected to induce shifts in the selective pressures that shape this adaptive variation. These changes in selective pressures will likely result in varying degrees of local climate maladaptation and spatial reshuffling of the underlying distributions of adaptive alleles. There is a growing interest in using population genomic data to help predict future disruptions to locally adaptive gene-environment associations. One motivation behind such work is to better understand how the effects of changing climate on populations’ short-term fitness could vary spatially across species ranges. Here we review the current use of genomic data to predict the disruption of local adaptation across current and future climates. After assessing goals and motivationsunderlying the approach, we review the main steps and associated statistical methods currently in use and explore our current understanding of the limits and future potential of using genomics to predict climate change (mal)adaptation.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction

Background

Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes.

Results

We generated 96.7 × 10⁹ 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 10⁹ paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females.

Conclusions

The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

Structure and Sequence of the Sex Determining Locus in Two Wild Populations of Nile Tilapia

In domesticated strains of the Nile tilapia, phenotypic sex has been linked to genetic variants on linkage groups 1, 20 and 23. This diversity of sex-loci might reflect a naturally polymorphic sex determination system in Nile tilapia, or it might be an artefact arising from the process of domestication. Here, we searched for sex-determiners in wild populations from Kpandu, Lake Volta (Ghana-West Africa), and from Lake Koka (Ethiopia-East Africa) that have not been subjected to any genetic manipulation. We analysed lab-reared families using double-digest Restriction Associated DNA sequencing (ddRAD) and analysed wild-caught males and females with pooled whole-genome sequencing (WGS). Strong sex-linked signals were found on LG23 in both populations, and sex-linked signals with LG3 were observed in Kpandu samples. WGS uncovered blocks of high sequence coverage, suggesting the presence of B chromosomes. We confirmed the existence of a tandem amh duplication in LG23 in both populations and determined its breakpoints between the oaz1 and dot1l genes. We found two common deletions of ~5 kb in males and confirmed the presence of both amhY and amh∆Y genes. Males from Lake Koka lack both the previously reported 234 bp deletion and the 5 bp frameshift-insertion that creates a premature stop codon in amh∆Y.

EXFI: Exon and splice graph prediction without a reference genome

For population genetic studies in nonmodel organisms, it is important to use every single source of genomic information. This paper presents EXFI, a Python pipeline that predicts the splice graph and exon sequences using an assembled transcriptome and raw whole-genome sequencing reads. The main algorithm uses Bloom filters to remove reads that are not part of the transcriptome, to predict the intron-exon boundaries, to then proceed to call exons from the assembly, and to generate the underlying splice graph. The results are returned in GFA1 format, which encodes both the predicted exon sequences and how they are connected to form transcripts. EXFI is written in Python, tested on Linux platforms, and the source code is available under the MIT License at https://github.com/jlanga/exfi.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

The direction of genital asymmetry is expressed stochastically in internally fertilizing anablepid fishes

Animal genitalia vary considerably across taxa, with divergence in many morphological traits, including striking departures from symmetry. Different mechanisms have been proposed to explain this diversity, mostly assuming that at least some of the phenotypic variation is heritable. However, heritability of the direction of genital asymmetry has been rarely determined. Anablepidae are internally fertilizing fish where the anal fin of males has been modified into an intromittent organ that transfers sperm into the gonopore of females. Males of anablepid fishes exhibit asymmetric genitalia, and both left- and right-sided individuals are commonly found at similar proportions within populations (i.e. antisymmetry). Although this polymorphism was described over a century ago, there have been no attempts to determine if genital asymmetry has a genetic basis and whether the different morphs are accumulating genetic differences, as might be expected since in some species females have also asymmetric gonopores and thereby can only be fertilized by compatible asymmetric males. We address this issue by combining breeding experiments with genome-wide data (ddRAD markers) in representative species of the two anablepid genera with asymmetric genitalia: Anableps and Jenynsia. Breeding experiments showed that all offspring were asymmetric, but their morphotype (i.e. right- or left-sided) was independent of parental morphotype, implying that the direction of asymmetry does not have a strong genetic component. Consistent with this conclusion, association analyses based on approximately 25 000 SNPs did not identify markers significantly associated with the direction of genital asymmetry and there was no evidence of population structure between left- and right-sided individuals. These results suggest that the direction of genital asymmetry in anablepid fishes might be stochastic, a commonly observed pattern in species with antisymmetry in morphological traits.

Combining genotype, phenotype, and environmental data to delineate site-adjusted provenance strategies for ecological restoration

Despite the importance of climate-adjusted provenancing to mitigate the effects of environmental change, climatic considerations alone are insufficient when restoring highly degraded sites. Here we propose a comprehensive landscape genomic approach to assist the restoration of moderately disturbed and highly degraded sites. To illustrate it we employ genomic data sets comprising thousands of single nucleotide polymorphisms from two plant species suitable for the restoration of iron-rich Amazonian Savannas. We first use a subset of neutral loci to assess genetic structure and determine the genetic neighbourhood size. We then identify genotype-phenotype-environment associations, map adaptive genetic variation, and predict adaptive genotypes for restoration sites. Whereas local provenances were found optimal to restore a moderately disturbed site, a mixture of genotypes seemed the most promising strategy to recover a highly degraded mining site. We discuss how our results can help define site-adjusted provenancing strategies, and argue that our methods can be more broadly applied to assist other restoration initiatives.

Disease swamps molecular signatures of genetic-environmental associations to abiotic factors in Tasmanian devil (Sarcophilus harrisii) populations

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors (i.e., disease). The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species' extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment (i.e., climatic variables, elevation, vegetation cover) and/or DFTD. We employ genetic-environment association analyses using 6886 SNPs from 3287 individuals sampled pre- and post-disease arrival across the devil's geographic range. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 365 unique loci linked to 71 genes, suggesting local adaptation to abiotic environment. The majority of candidate loci detected pre-DFTD are not detected post-DFTD arrival. Several post-DFTD candidate loci are associated with disease prevalence and were in linkage disequilibrium with genes involved in tumor suppression and immune response. Loss of apparent signal of abiotic local adaptation post-disease suggests swamping by strong selection resulting from the rapid onset of DFTD.

Genome-wide association analyses based on whole-genome sequencing of Protosalanx hyalocranius provide insights into sex determination of Salangid fishes

Identification of sex determination system and sex-determining genes have important implications in conservation, ecology and evolution. However, much remains to be discovered about the evolution of different sexual determination systems in teleost fishes, of which the mechanisms of sex determination are remarkably variable. In the present study, the whole genomes of 20 males and 20 females of a Salangid fish, Protosalanx hyalocranius, were sequenced and genome wide association analyses were conducted to uncover its sex determination system and putative sex-determining genes. A total of 150 SNPs were significantly associated with sex, which showed high differentiation between sexes (F_ST ranged from 0.245 to 0.556). Of the 150 sex-associated SNPs, 76 SNPs displayed sex specificity with even coverage of depth and were female heterogametic, which suggested a ZZ/ZW sex determination system. Interestingly, one scaffold containing sex-specific SNPs displayed synteny to the sex chromosome of medaka. Annotations of sex-associated loci suggested that both transcriptional regulators (e.g., FOX genes) and secreted hormones and their receptors might be involved in the sex determination/differentiation of P. hyalocranius. More strikingly, we found a nonsense mutation in one copy of GALNT homology gene of all females, which suggested that "Z dosage" effect might play a vital role in the processes of sex determination/differentiation. These sex-specific loci could be a valuable resource for further research on sex determination of Salangid fishes and the results could contribute to the understanding of sex determination mechanisms and the evolution of sex chromosome in teleost fishes.

Simulation with RADinitio improves RADseq experimental design and sheds light on sources of missing data

Restriction-site associated DNA sequencing (RADseq) has become a powerful and versatile tool in modern population genomics, enabling large-scale evolutionary and genomic analyses in otherwise inaccessible biological systems. With its widespread use, different variants on the protocol have been developed to suit specific experimental needs. Researchers face the challenge of choosing the optimal molecular and sequencing protocols for their reduced representation experimental design, an often-complicated process. Strategic errors can lead to biased data generation that has reduced power to answer biological questions. Here, we present RADinitio, simulation software for the selection and optimization of RADseq experiments via the generation of sequencing data that behave similarly to empirical sources. RADinitio provides an evolutionary simulation of populations, implementation of various RADseq protocols with customizable parameters, and thorough assessment of missing data. We test the efficacy of the software using different RAD protocols across several organisms, highlighting the importance of protocol selection on the magnitude and quality of data acquired. Additionally, we test the effects of RAD library preparation and sequencing on allelic dropout, observing that library preparation and sequencing often contributes more to missing alleles than population-level variation.

Demography and adaptation promoting evolutionary transitions in a mammalian genus that diversified during the Pleistocene

Species that evolved in temperate regions during the Pleistocene experienced periods of extreme climatic transitions. Consequent population fragmentation and dynamics had the potential to generate small, isolated populations where the influence of genetic drift would be expected to be strong. We use comparative genomics to assess the evolutionary influence of historical demographics and natural selection through a series of transitions associated with the formation of the genus Capreolus, speciation within this genus during the Quaternary and during divergence among European roe deer (C. capreolus) populations. Our analyses were facilitated by the generation of a new high-coverage reference genome for the Siberian roe deer (C. pygargus). We find progressive reductions in effective population size (N_e ), despite very large census sizes in modern C. capreolus populations and show that low N_e has impacted the C. capreolus genome, reducing diversity and increasing linkage disequilibrium. Even so, we find evidence for natural selection shared among C. capreolus populations, including a historically documented founder population that has been through a severe bottleneck. During each phylogenetic transition there is evidence for selection (from dN/dS and nucleotide diversity tests), including at loci associated with diapause (delayed embryonic development), a phenotype restricted to this genus among the even-toed ungulates. Together these data allow us to assess expectations for the origin and diversification of a mammalian genus during a period of extreme environmental change.

Within-Generation Polygenic Selection Shapes Fitness-Related Traits across Environments in Juvenile Sea Bream

Understanding the genetic underpinnings of fitness trade-offs across spatially variable environments remains a major challenge in evolutionary biology. In Mediterranean gilthead sea bream, first-year juveniles use various marine and brackish lagoon nursery habitats characterized by a trade-off between food availability and environmental disturbance. Phenotypic differences among juveniles foraging in different habitats rapidly appear after larval settlement, but the relative role of local selection and plasticity in phenotypic variation remains unclear. Here, we combine phenotypic and genetic data to address this question. We first report correlations of opposite signs between growth and condition depending on juvenile habitat type. Then, we use single nucleotide polymorphism (SNP) data obtained by Restriction Associated DNA (RAD) sequencing to search for allele frequency changes caused by a single generation of spatially varying selection between habitats. We found evidence for moderate selection operating at multiple loci showing subtle allele frequency shifts between groups of marine and brackish juveniles. We identified subsets of candidate outlier SNPs that, in interaction with habitat type, additively explain up to 3.8% of the variance in juvenile growth and 8.7% in juvenile condition; these SNPs also explained significant fraction of growth rate in an independent larval sample. Our results indicate that selective mortality across environments during early-life stages involves complex trade-offs between alternative growth strategies.

Emerging Frontiers in the Study of Molecular Evolution

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

Early diversification and permeable species boundaries in the Mediterranean firs

BACKGROUND AND AIMS

Inferring the evolutionary relationships of species and their boundaries is critical in order to understand patterns of diversification and their historical drivers. Despite Abies (Pinaceae) being the second most diverse group of conifers, the evolutionary history of Circum-Mediterranean firs (CMFs) remains under debate.

METHODS

We used restriction site-associated DNA sequencing (RAD-seq) on all proposed CMF taxa to investigate their phylogenetic relationships and taxonomic status.

KEY RESULTS

Based on thousands of genome-wide single nucleotide polymorphisms (SNPs), we present here the first formal test of species delimitation, and the first fully resolved, complete species tree for CMFs. We discovered that all previously recognized taxa in the Mediterranean should be treated as independent species, with the exception of Abies tazaotana and Abies marocana. An unexpectedly early pulse of speciation in the Oligocene-Miocene boundary is here documented for the group, pre-dating previous hypotheses by millions of years, revealing a complex evolutionary history encompassing both ancient and recent gene flow between distant lineages.

CONCLUSIONS

Our phylogenomic results contribute to shed light on conifers' diversification. Our efforts to resolve the CMF phylogenetic relationships help refine their taxonomy and our knowledge of their evolution.

Genomic divergence reveals unique populations among Indian Yaks

The present study focused upon identification of genome-wide SNPs through the reduced representation approach and to study the genomic divergence of the Indian yak populations. A total of 80 samples belonging to Arunachali yak (N = 20), Himachali yak (N = 20), Ladakhi yak (N = 20) and Sikkimi yak (N = 20) of India were used in the study. The results of the study revealed a total of 579575 high quality SNPs along with 50319 INDELs in the Indian yaks. The observed heterozygosity was found to be high in Himachali yak, followed by Arunachali yak, Ladakhi yak and Sikkimi yaks. The Sikkimi yaks was found to be genetically distant, followed by Ladakhi yaks which was observed to have some few individuals from Arunachali and Himachali yaks. Arunachali and Himachali yaks are found to get clustered together and are genetically similar. The study provides evidence about the genomic diversity in the Indian yak populations and information generated in the present study may help to formulate a suitable breeding plan for endangered Indian yaks. Moreover, the unique yak populations identified in the study will further help to focus attention for future characterization and prioritization of the animals for conservation purposes through the ddRAD approach.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

Stacking up RADSeq assembly programs: From complete hit to completely abysmal

Decreasing sequencing costs have driven a rapid expansion of novel genotyping methods. One of these methods is the exploitation of restriction enzyme cut sites to generate genome-wide but reduced representation sequencing libraries (RRLs), alternatively termed genotyping by sequencing or restriction-site associated DNA sequencing. Without a reference genome, the resulting short sequence reads must be assembled de novo. There are many possible assembly programs, most not explicitly developed for RRL data, and we know little of their effectiveness. In this issue of Molecular Ecology Resources, LaCava et al. (2020) systematically evaluate six commonly used programs and two commonly varied parameters for complete and accurate assembly of RRLs, using simulated double digests of Homo sapiens and Arabidopsis thaliana genomes with varied mutation rates and types. The authors find substantial variation in performance across assembly programs. The most consistently high-performing assembler is infrequently used in their literature survey (CD-HIT; Li and Godzik, 2006), while several others fail to produce complete, accurate assemblies under many conditions. LaCava et al. additionally recommend best practices in parameter choice and evaluation of future assembly programs-advice that molecular ecologists working to assemble sequences of all kinds should take to heart.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

Strong Selection Against Early Generation Hybrids in Joshua Tree Hybrid Zone Not Explained by Pollinators Alone

Coevolution frequently plays an important role in diversification, but the role of obligate pollination mutualisms in the maintenance of hybrid zones has rarely been investigated. Like most members of the genus Yucca, the two species of Joshua tree (Yucca brevifolia and Yucca jaegeriana) are involved in a tightly coevolved mutualism with yucca moths. There is strong evidence of a history of coevolution between Joshua trees and their moth pollinators. We use a geographic clines approach in the Joshua tree hybrid zone to ask if selection by the moths may currently contribute to maintaining separation between these species. We compare genomic, phenotypic, and pollinator frequency clines to test whether pollinators maintain the hybrid zone or follow it as passive participants. The results reveal dramatic overlapping genomic and pollinator clines, consistent with a narrow hybrid zone maintained by strong selection. Wider phenotypic clines and a chloroplast genomic cline displaced opposite the expected direction suggest that pollinators are not the main source of selection maintaining the hybrid zone. Rather, it seems that high levels of reproductive isolation, likely acting through multiple barriers and involving many parts of the genome, keep the hybrid zone narrow.

Genomic assessment of local adaptation in dwarf birch to inform assisted gene flow

When populations of a rare species are small, isolated and declining under climate change, some populations may become locally maladapted. Detecting this maladaptation may allow effective rapid conservation interventions, even if based on incomplete knowledge. Population maladaptation may be estimated by finding genome-environment associations (GEA) between allele frequencies and environmental variables across a local species range, and identifying populations whose allele frequencies do not fit with these trends. We can then design assisted gene flow strategies for maladapted populations, to adjust their allele frequencies, entailing lower levels of intervention than with undirected conservation action. Here, we investigate this strategy in Scottish populations of the montane plant dwarf birch (Betula nana). In genome-wide restriction site-associated single nucleotide polymorphism (SNP) data, we found 267 significant associations between SNP loci and environmental variables. We ranked populations by maladaptation estimated using allele frequency deviation from the general trends at these loci; this gave a different prioritization for conservation action than the Shapely Index, which seeks to preserve rare neutral variation. Populations estimated to be maladapted in their allele frequencies at loci associated with annual mean temperature were found to have reduced catkin production. Using an environmental niche modelling (ENM) approach, we found annual mean temperature (35%), and mean diurnal range (15%), to be important predictors of the dwarf birch distribution. Intriguingly, there was a significant correlation between the number of loci associated with each environmental variable in the GEA and the importance of that variable in the ENM. Together, these results suggest that the same environmental variables determine both adaptive genetic variation and species range in Scottish dwarf birch. We suggest an assisted gene flow strategy that aims to maximize the local adaptation of dwarf birch populations under climate change by matching allele frequencies to current and future environments.

Evolutionary genetics in insect phenotypic radiations: the value of a comparative genomic approach

Expanding genome sequencing and transgenic technologies are enabling the discovery of genes driving phenotypic diversity across insect taxa. Limitations in downstream functional genetic approaches, however, have been an obstacle for developing non-model systems for evolutionary genetics. Phenotypically diverse radiations, such as those exhibiting convergence and divergence as a result of mimicry, are ideal for evolutionary genetics as they can lead to insights using comparative genomic approaches alone. The varied and repeated instances of phenotypes in highly polymorphic systems allow assessment of whether similar loci are repeatedly targeted by selection and can inform how alleles sort across lineages. Comparative genomics of these taxa can be used to decipher components of gene regulatory networks, dissect regulatory regions, and validate genes.

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Accuracy of de novo assembly of DNA sequences from double-digest libraries varies substantially among software

Advances in DNA sequencing have made it feasible to gather genomic data for non-model organisms and large sets of individuals, often using methods for sequencing subsets of the genome. Several of these methods sequence DNA associated with endonuclease restriction sites (various RAD and GBS methods). For use in taxa without a reference genome, these methods rely on de novo assembly of fragments in the sequencing library. Many of the software options available for this application were originally developed for other assembly types and we do not know their accuracy for reduced representation libraries. To address this important knowledge gap, we simulated data from the Arabidopsis thaliana and Homo sapiens genomes and compared de novo assemblies by six software programs that are commonly used or promising for this purpose (ABySS, CD-HIT, Stacks, Stacks2, Velvet and VSEARCH). We simulated different mutation rates and types of mutations, and then applied the six assemblers to the simulated data sets, varying assembly parameters. We found substantial variation in software performance across simulations and parameter settings. ABySS failed to recover any true genome fragments, and Velvet and VSEARCH performed poorly for most simulations. Stacks and Stacks2 produced accurate assemblies of simulations containing SNPs, but the addition of insertion and deletion mutations decreased their performance. CD-HIT was the only assembler that consistently recovered a high proportion of true genome fragments. Here, we demonstrate the substantial difference in the accuracy of assemblies from different software programs and the importance of comparing assemblies that result from different parameter settings.

New insights from RADseq data on differentiation in the Hottentot golden mole species complex from South Africa

Golden moles (Family Chrysochloridae) are small subterranean mammals, endemic to sub-Saharan Africa, and many of the 21 species are listed as threatened on the IUCN Red List. Most species have highly restricted ranges; however two species, the Hottentot golden mole (Amblysomus hottentotus) and the Cape golden mole (Chrysochloris asiatica) have relatively wide ranges. We recently uncovered cryptic diversity within A. hottentotus, through a phylogeographic analysis of this taxon using two mitochondrial gene regions and a nuclear intron. To further investigate this cryptic diversity, we generated nuclear SNP data from across the genome of A. hottentotus, by means of double-digest restriction-site associated DNA sequencing (ddRADSeq), and mapped reads to the Cape golden mole genome. We conducted a phylogenetic analysis and investigated population differentiation. Our results support the distinctiveness of A. h. meesteri. Furthermore, we provide evidence from nuclear SNPs in support of our previous finding that Central coastal samples represent a unique cryptic lineage that is highly divergent from A. h. pondoliae farther south. Although mtDNA suggests that Umtata may represent a unique lineage sister to A. h. longiceps, mito-nuclear discordance from our RADseq data indicate that these samples may instead be closer to A. h. pondoliae, and therefore may not represent a distinct lineage. We stress the importance of recognizing that understudied populations, such as that of Umtata, may represent populations or ESUs under threat and in need of conservation attention. We present a high-quality filtered SNP dataset, comprising thousands of SNPs, which may serve as a useful resource for future golden mole studies. We have thus added to the growing body of research demonstrating the power and utility of RADseq to investigate population differentiation.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Replicated Landscape Genomics Identifies Evidence of Local Adaptation to Urbanization in Wood Frogs

Native species that persist in urban environments may benefit from local adaptation to novel selection factors. We used double-digest restriction-side associated DNA (RAD) sequencing to evaluate shifts in genome-wide genetic diversity and investigate the presence of parallel evolution associated with urban-specific selection factors in wood frogs (Lithobates sylvaticus). Our replicated paired study design involved 12 individuals from each of 4 rural and urban populations to improve our confidence that detected signals of selection are indeed associated with urbanization. Genetic diversity measures were less for urban populations; however, the effect size was small, suggesting little biological consequence. Using an FST outlier approach, we identified 37 of 8344 genotyped single nucleotide polymorphisms with consistent evidence of directional selection across replicates. A genome-wide association study analysis detected modest support for an association between environment type and 12 of the 37 FST outlier loci. Discriminant analysis of principal components using the 37 FST outlier loci produced correct reassignment for 87.5% of rural samples and 93.8% of urban samples. Eighteen of the 37 FST outlier loci mapped to the American bullfrog (Rana [Lithobates] catesbeiana) genome, although none were in coding regions. This evidence of parallel evolution to urban environments provides a powerful example of the ability of urban landscapes to direct evolutionary processes.

Genome-wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab-adapted and wave-adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome-wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome-wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait-associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab-adapted and wave-adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.

Ecological speciation in sympatric palms: 3. Genetic map reveals genomic islands underlying species divergence in Howea

Although it is now widely accepted that speciation can occur in the face of continuous gene flow, with little or no spatial separation, the mechanisms and genomic architectures that permit such divergence are still debated. Here, we examined speciation in the face of gene flow in the Howea palms of Lord Howe Island, Australia. We built a genetic map using a novel method applicable to long-lived tree species, combining it with double digest restriction site-associated DNA sequencing of multiple individuals. Based upon various metrics, we detected 46 highly differentiated regions throughout the genome, four of which contained genes with functions that are particularly relevant to the speciation scenario for Howea, specifically salt and drought tolerance.