The genomics of adaptation

Genomic and Phenotypic Adaptations of Rattus tanezumi to Cold Limit Its Further Northward Expansion and Range Overlap with R. norvegicus

Global climate change has led to shifts in the distribution ranges of many terrestrial species, promoting their migration from lower altitudes or latitudes to higher ones. Meanwhile, successful invaders have developed genetic adaptations enabling the colonization of new environments. Over the past 40 years, Rattus tanezumi (RT) has expanded into northern China (Northwest and North China) from its southern origins. We studied the cold adaptation of RT and its potential for northward expansion by comparing it with sympatric Rattus norvegicus (RN), which is well adapted to cold regions. Through population genomic analysis, we revealed that the invading RT rats have split into three distinct populations: the North, Northwest, and Tibetan populations. The first two populations exhibited high genetic diversity, while the latter population showed remarkably low genetic diversity. These rats have developed various genetic adaptations to cold, arid, hypoxic, and high-UV conditions. Cold acclimation tests revealed divergent thermoregulation between RT and RN. Specifically, RT exhibited higher brown adipose tissue activity and metabolic rates than did RN. Transcriptome analysis highlighted changes in genes regulating triglyceride catabolic processes in RT, including Apoa1 and Apoa4, which were upregulated, under selection and associated with local adaptation. In contrast, RN showed changes in carbohydrate metabolism genes. Despite the cold adaptation of RT, we observed genotypic and phenotypic constraints that may limit its ability to cope with severe low temperatures farther north. Consequently, it is less likely that RT rats will invade and overlap with RN rats in farther northern regions.

Whole-genome resequencing reveals genetic diversity and adaptive evolution in Chinese honeybee (Apis cerana cerana) in Guizhou, China

Introduction: Guizhou Province, characterized by complex and diverse geographic and climatic environments, has rich genetic resources for the Chinese honeybee (Apis cerana cerana) and is one of the main bee-producing areas in China. However, research on the genetic diversity of Chinese honeybee in the Guizhou region is very limited, despite implications for conservation of biodiversity. Methods: In this study, we analyzed the genetic diversity, differentiation, and selection signals based on 116 Chinese honeybees from 12 regions in Guizhou Province using whole-genome sequencing. Results: We identified 1,400,430 high-quality SNPs across all samples. A population structure analysis revealed two independent genetic subgroups of Chinese honeybees in Guizhou, a Yunnan-Guizhou Plateau population in western Guizhou and a hilly-mountainous population in eastern Guizhou. The average nucleotide diversity (Pi) ranged from 0.00138 to 0.00161 and average expected heterozygosity (He) ranged from 0.2592 to 0.2604. The average genetic differentiation index (F ST) for Chinese honeybees in pairwise comparisons of 12 regions ranged from 0.0094 to 0.0293. There was clear genetic differentiation between the western plateau and the eastern hilly mountainous areas of Guizhou; however, F ST values between the eastern and western populations ranged from 0.0170 to 0.0293, indicating a low degree of differentiation. A genome-wide scan revealed a number of genes under selection in the Yunnan-Guizhou Plateau environment. These genes were related to growth and development, reproduction, and cold resistance, and several candidate genes involved in environmental adaptation were identified, including CTR, MAPK, MAST, HSF, and MKKK. Discussion: The results of the present study provide important theoretical bases for the conservation, evaluation, development, and utilization of genetic resources for Chinese honeybees in the Guizhou region and for further investigations of environmental adaptation and underlying mechanisms in the species.

First genome-wide data from Italian European beech (Fagus sylvatica L.): Strong and ancient differentiation between Alps and Apennines

The European beech (Fagus sylvatica L.) is one of the most widespread forest trees in Europe whose distribution and intraspecific diversity has been largely shaped by repeated glacial cycles. Previous studies, mainly based on palaeobotanical evidence and a limited set of chloroplast and nuclear genetic markers, highlighted a complex phylogeographic scenario, with southern and western Europe characterized by a rather heterogeneous genetic structure, as a result of recolonization from different glacial refugia. Despite its ecological and economic importance, the genome of this broad-leaved tree has only recently been assembled, and its intra-species genomic diversity is still largely unexplored. Here, we performed whole-genome resequencing of nine Italian beech individuals sampled from two stands located in the Alpine and Apennine mountain ranges. We investigated patterns of genetic diversity at chloroplast, mitochondrial and nuclear genomes and we used chloroplast genomes to reconstruct a temporally-resolved phylogeny. Results allowed us to test European beech differentiation on a whole-genome level and to accurately date their divergence time. Our results showed comparable, relatively high levels of genomic diversity in the two populations and highlighted a clear differentiation at chloroplast, mitochondrial and nuclear genomes. The molecular clock analysis indicated an ancient split between the Alpine and Apennine populations, occurred between the Günz and the Riss glaciations (approximately 660 kyrs ago), suggesting a long history of separation for the two gene pools. This information has important conservation implications in the context of adaptation to ongoing climate changes.

Mosaic translation hypothesis: chimeric polypeptides produced via multiple ribosomal frameshifting as a basis for adaptability

How many different proteins can be produced from a single spliced transcript? Genome annotation projects overlook the coding potential of reading frames other than that of the reference open reading frames (refORFs). Recently, alternative open reading frames (altORFs) and their translational products, alternative proteins, have been shown to carry out important functions in various organisms. AltORFs overlapping refORFs or other altORFs in a different reading frame may be involved in one fundamental mechanism so far overlooked. A few years ago, it was proposed that altORFs may act as building blocks for chimeric (mosaic) polypeptides, which are produced via multiple ribosomal frameshifting events from a single mature transcript. We adopt terminology from that earlier discussion and call this mechanism mosaic translation. This way of extracting and combining genetic information may significantly increase proteome diversity. Thus, we hypothesize that this mechanism may have contributed to the flexibility and adaptability of organisms to a variety of environmental conditions. Specialized ribosomes acting as sensors probably played a central role in this process. Importantly, mosaic translation may be the main source of protein diversity in genomes that lack alternative splicing. The idea of mosaic translation is a testable hypothesis, although its direct demonstration is challenging. Should mosaic translation occur, we would currently highly underestimate the complexity of translation mechanisms and thus the proteome.

Belgica antarctica (Diptera: Chironomidae): A natural model organism for extreme environments

Belgica antarctica (Diptera: Chironomidae), a brachypterous midge endemic to the maritime Antarctic, was first described in 1900. Over more than a century of study, a vast amount of information has been compiled on the species (3 750 000 Google search results as of January 10, 2021), encompassing its ecology and biology, life cycle and reproduction, polytene chromosomes, physiology, biochemistry and, increasingly, omics. In 2014, B. antarctica's genome was sequenced, further boosting research. Certain developmental stages can be cultured successfully in the laboratory. Taken together, this wealth of information allows the species to be viewed as a natural model organism for studies of adaptation and function in extreme environments.

Chromosome-level de novo genome assembly of Telopea speciosissima (New South Wales waratah) using long-reads, linked-reads and Hi-C

Telopea speciosissima, the New South Wales waratah, is an Australian endemic woody shrub in the family Proteaceae. Waratahs have great potential as a model clade to better understand processes of speciation, introgression and adaptation, and are significant from a horticultural perspective. Here, we report the first chromosome-level genome for T. speciosissima. Combining Oxford Nanopore long-reads, 10x Genomics Chromium linked-reads and Hi-C data, the assembly spans 823 Mb (scaffold N50 of 69.0 Mb) with 97.8% of Embryophyta BUSCOs "Complete". We present a new method in Diploidocus (https://github.com/slimsuite/diploidocus) for classifying, curating and QC-filtering scaffolds, which combines read depths, k-mer frequencies and BUSCO predictions. We also present a new tool, DepthSizer (https://github.com/slimsuite/depthsizer), for genome size estimation from the read depth of single-copy orthologues and estimate the genome size to be approximately 900 Mb. The largest 11 scaffolds contained 94.1% of the assembly, conforming to the expected number of chromosomes (2n = 22). Genome annotation predicted 40,158 protein-coding genes, 351 rRNAs and 728 tRNAs. We investigated CYCLOIDEA (CYC) genes, which have a role in determination of floral symmetry, and confirm the presence of two copies in the genome. Read depth analysis of 180 "Duplicated" BUSCO genes using a new tool, DepthKopy (https://github.com/slimsuite/depthkopy), suggests almost all are real duplications, increasing confidence in the annotation and highlighting a possible need to revise the BUSCO set for this lineage. The chromosome-level T. speciosissima reference genome (Tspe_v1) provides an important new genomic resource of Proteaceae to support the conservation of flora in Australia and further afield.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

Genome sequencing and population genomics modeling provide insights into the local adaptation of weeping forsythia

Understanding the genetic basis underlying the local adaptation of nonmodel species is a fundamental goal in evolutionary biology. In this study, we explored the genetic mechanisms of the local adaptation of Forsythia suspensa using genome sequence and population genomics data obtained from specific-locus amplified fragment sequencing. We assembled a high-quality reference genome of weeping forsythia (Scaffold N50 = 7.3 Mb) using ultralong Nanopore reads. Then, genome-wide comparative analysis was performed for 15 natural populations of weeping forsythia across its current distribution range. Our results revealed that candidate genes associated with local adaptation are functionally correlated with solar radiation, temperature and water variables across heterogeneous environmental scenarios. In particular, solar radiation during the period of fruit development and seed drying after ripening, cold, and drought significantly contributed to the adaptive differentiation of F. suspensa. Natural selection exerted by environmental factors contributed substantially to the population genetic structure of F. suspensa. Our results supported the hypothesis that adaptive differentiation should be highly pronounced in the genes involved in signal crosstalk between different environmental variables. Our population genomics study of F. suspensa provides insights into the fundamental genetic mechanisms of the local adaptation of plant species to climatic gradients.

Genomes reveal drastic and recurrent phenotypic divergence in firetip skipper butterflies (Hesperiidae: Pyrrhopyginae)

Biologists marvel at the powers of adaptive convergence, when distantly related animals look alike. While mimetic wing patterns of butterflies have fooled predators for millennia, entomologists inferred that mimics were distant relatives despite similar appearance. However, the obverse question has not been frequently asked. Who are the close relatives of mimetic butterflies and what are their features? As opposed to close convergence, divergence from a non-mimetic relative would also be extreme. When closely related animals look unalike, it is challenging to pair them. Genomic analysis promises to elucidate evolutionary relationships and shed light on molecular mechanisms of divergence. We chose the firetip skipper butterfly as a model due to its phenotypic diversity and abundance of mimicry. We sequenced and analysed whole genomes of nearly 120 representative species. Genomes partitioned this subfamily Pyrrhopyginae into five tribes (1 new), 23 genera and, additionally, 22 subgenera (10 new). The largest tribe Pyrrhopygini is divided into four subtribes (three new). Surprisingly, we found five cases where a uniquely patterned butterfly was formerly placed in a genus of its own and separately from its close relatives. In several cases, extreme and rapid phenotypic divergence involved not only wing patterns but also the structure of the male genitalia. The visually striking wing pattern difference between close relatives frequently involves disappearance or suffusion of spots and colour exchange between orange and blue. These differences (in particular, a transition between unspotted black and striped wings) happen recurrently on a short evolutionary time scale, and are therefore probably achieved by a small number of mutations.

Tigers of the World: Genomics and Conservation

Of all the big cats, or perhaps of all the endangered wildlife, the tiger may be both the most charismatic and most well-recognized flagship species in the world. The rapidly changing field of molecular genetics, particularly advances in genome sequencing technologies, has provided new tools to reconstruct what characterizes a tiger. Here we review how applications of molecular genomic tools have been used to depict the tiger's ancestral roots, phylogenetic hierarchy, demographic history, morphological diversity, and genetic patterns of diversification on both temporal and geographical scales. Tiger conservation, stabilization, and management are important areas that benefit from use of these genome resources for developing survival strategies for this charismatic megafauna both in situ and ex situ.

Mitochondria, sex and variation in routine metabolic rate

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (ṀO_{2 routine} ), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with ṀO_{2 routine} , and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with ṀO_{2 routine} involved different SNPs if females and males were analysed separately, and there was a sex-dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex-specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in ṀO_{2 routine} . Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait-associated SNPs, suggesting that epistatic effects on ṀO_{2 routine} are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex-specific.

The genomic basis of adaptation to high-altitude habitats in the eastern honey bee (Apis cerana)

The eastern honey bee (Apis cerana) is of central importance for agriculture in Asia. It has adapted to a wide variety of environmental conditions across its native range in southern and eastern Asia, which includes high-altitude regions. eastern honey bees inhabiting mountains differ morphologically from neighbouring lowland populations and may also exhibit differences in physiology and behaviour. We compared the genomes of 60 eastern honey bees collected from high and low altitudes in Yunnan and Gansu provinces, China, to infer their evolutionary history and to identify candidate genes that may underlie adaptation to high altitude. Using a combination of F_ST -based statistics, long-range haplotype tests and population branch statistics, we identified several regions of the genome that appear to have been under positive selection. These candidate regions were strongly enriched for coding sequences and had high haplotype homozygosity and increased divergence specifically in highland bee populations, suggesting they have been subjected to recent selection in high-altitude habitats. Candidate loci in these genomic regions included genes related to reproduction and feeding behaviour in honey bees. Functional investigation of these candidate loci is necessary to fully understand the mechanisms of adaptation to high-altitude habitats in the eastern honey bee.

Tolerance traits related to climate change resilience are independent and polygenic

The resilience of organisms to climate change through adaptive evolution is dependent on the extent of genetically based variation in key phenotypic traits and the nature of genetic associations between them. For aquatic animals, upper thermal tolerance and hypoxia tolerance are likely to be a important determinants of sensitivity to climate change. To determine the genetic basis of these traits and to detect associations between them, we compared naturally occurring populations of two subspecies of Atlantic killifish, Fundulus heteroclitus, that differ in both thermal and hypoxia tolerance. Multilocus association mapping demonstrated that 47 and 35 single nucleotide polymorphisms (SNPs) explained 43.4% and 51.9% of variation in thermal and hypoxia tolerance, respectively, suggesting that genetic mechanisms underlie a substantial proportion of variation in each trait. However, no explanatory SNPs were shared between traits, and upper thermal tolerance varied approximately linearly with latitude, whereas hypoxia tolerance exhibited a steep phenotypic break across the contact zone between the subspecies. These results suggest that upper thermal tolerance and hypoxia tolerance are neither phenotypically correlated nor genetically associated, and thus that rates of adaptive change in these traits can be independently fine-tuned by natural selection. This modularity of important traits can underpin the evolvability of organisms to complex future environmental change.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback

Heterogeneous genetic divergence can accumulate across the genome when populations adapt to different habitats while still exchanging alleles. How long does diversification take and how much of the genome is affected? When divergence occurs in parallel from standing genetic variation, how often are the same haplotypes involved? We explore these questions using restriction site-associated DNA sequencing genotyping data and show that broad-scale genomic repatterning, fueled by copious standing variation, can emerge in just dozens of generations in replicate natural populations of threespine stickleback fish (Gasterosteus aculeatus). After the catastrophic 1964 Alaskan earthquake, marine stickleback colonized newly created ponds on seismically uplifted islands. We find that freshwater fish in these young ponds differ from their marine ancestors across the same genomic segments previously shown to have diverged in much older lake populations. Outside of these core divergent regions the genome shows no population structure across the ocean-freshwater divide, consistent with strong local selection acting in alternative environments on stickleback populations still connected by significant gene flow. Reinforcing this inference, a majority of divergent haplotypes that are at high frequency in ponds are detectable in the sea, even across great geographic distances. Building upon previous population genomics work in this model species, our data suggest that a long history of divergent selection and gene flow among stickleback populations in oceanic and freshwater habitats has maintained polymorphisms of alternatively adapted DNA sequences that facilitate parallel evolution.

Estimating Time to the Common Ancestor for a Beneficial Allele

The haplotypes of a beneficial allele carry information about its history that can shed light on its age and the putative cause for its increase in frequency. Specifically, the signature of an allele's age is contained in the pattern of variation that mutation and recombination impose on its haplotypic background. We provide a method to exploit this pattern and infer the time to the common ancestor of a positively selected allele following a rapid increase in frequency. We do so using a hidden Markov model which leverages the length distribution of the shared ancestral haplotype, the accumulation of derived mutations on the ancestral background, and the surrounding background haplotype diversity. Using simulations, we demonstrate how the inclusion of information from both mutation and recombination events increases accuracy relative to approaches that only consider a single type of event. We also show the behavior of the estimator in cases where data do not conform to model assumptions, and provide some diagnostics for assessing and improving inference. Using the method, we analyze population-specific patterns in the 1000 Genomes Project data to estimate the timing of adaptation for several variants which show evidence of recent selection and functional relevance to diet, skin pigmentation, and morphology in humans.

Revisiting the phylogeography, demography and taxonomy of the frog genus Ptychadena in the Ethiopian highlands with the use of genome-wide SNP data

Understanding the diversification of biological lineages is central to evolutionary studies. To properly study the process of speciation, it is necessary to link micro-evolutionary studies with macro-evolutionary mechanisms. Micro-evolutionary studies require proper sampling across a taxon's range to adequately infer genetic diversity. Here we use the grass frogs of the genus Ptychadena from the Ethiopian highlands as a model to study the process of lineage diversification in this unique biodiversity hotspot. We used thousands of genome-wide SNPs obtained from double digest restriction site associated DNA sequencing (ddRAD-seq) in populations of the Ptychadena neumanni species complex from the Ethiopian highlands in order to infer their phylogenetic relationships and genetic structure, as well as to study their demographic history. Our genome-wide phylogenetic study supports the existence of approximately 13 lineages clustered into 3 species groups. Our phylogenetic and phylogeographic reconstructions suggest that those endemic lineages diversified in allopatry, and subsequently specialized to different habitats and elevations. Demographic analyses point to a continuous decrease in the population size across the majority of lineages and populations during the Pleistocene, which is consistent with a continuous period of aridification that East Africa experienced since the Pliocene. We discuss the taxonomic implications of our analyses and, in particular, we warn against the recent practice to solely use Bayesian species delimitation methods when proposing taxonomic changes.

Investigating the genomic basis of discrete phenotypes using a Pool-Seq-only approach: New insights into the genetics underlying colour variation in diverse taxa

While large-scale genomic approaches are increasingly revealing the genetic basis of polymorphic phenotypes such as colour morphs, such approaches are almost exclusively conducted in species with high-quality genomes and annotations. Here, we use Pool-Seq data for both genome assembly and SNP frequency estimation, followed by scanning for F_ST outliers to identify divergent genomic regions. Using paired-end, short-read sequencing data from two groups of individuals expressing divergent phenotypes, we generate a de novo rough-draft genome, identify SNPs and calculate genomewide F_ST differences between phenotypic groups. As genomes generated by Pool-Seq data are highly fragmented, we also present an approach for super-scaffolding contigs using existing protein-coding data sets. Using this approach, we reanalysed genomic data from two recent studies of birds and butterflies investigating colour pattern variation and replicated their core findings, demonstrating the accuracy and power of a Pool-Seq-only approach. Additionally, we discovered new regions of high divergence and new annotations that together suggest novel parallels between birds and butterflies in the origins of their colour pattern variation.

flyDIVaS: A Comparative Genomics Resource for Drosophila Divergence and Selection

With arguably the best finished and expertly annotated genome assembly, Drosophila melanogaster is a formidable genetics model to study all aspects of biology. Nearly a decade ago, the 12 Drosophila genomes project expanded D. melanogaster’s breadth as a comparative model through the community-development of an unprecedented genus- and genome-wide comparative resource. However, since its inception, these datasets for evolutionary inference and biological discovery have become increasingly outdated, outmoded, and inaccessible. Here, we provide an updated and upgradable comparative genomics resource of Drosophila divergence and selection, flyDIVaS, based on the latest genomic assemblies, curated FlyBase annotations, and recent OrthoDB orthology calls. flyDIVaS is an online database containing D. melanogaster-centric orthologous gene sets, CDS and protein alignments, divergence statistics (% gaps, d_N, d_S, d_N/d_S), and codon-based tests of positive Darwinian selection. Out of 13,920 protein-coding D. melanogaster genes, ∼80% have one aligned ortholog in the closely related species, D. simulans, and ∼50% have 1–1 12-way alignments in the original 12 sequenced species that span over 80 million yr of divergence. Genes and their orthologs can be chosen from four different taxonomic datasets differing in phylogenetic depth and coverage density, and visualized via interactive alignments and phylogenetic trees. Users can also batch download entire comparative datasets. A functional survey finds conserved mitotic and neural genes, highly diverged immune and reproduction-related genes, more conspicuous signals of divergence across tissue-specific genes, and an enrichment of positive selection among highly diverged genes. flyDIVaS will be regularly updated and can be freely accessed at www.flydivas.info. We encourage researchers to regularly use this resource as a tool for biological inference and discovery, and in their classrooms to help train the next generation of biologists to creatively use such genomic big data resources in an integrative manner.

Genomic Response to Selection for Predatory Behavior in a Mammalian Model of Adaptive Radiation

If genetic architectures of various quantitative traits are similar, as studies on model organisms suggest, comparable selection pressures should produce similar molecular patterns for various traits. To test this prediction, we used a laboratory model of vertebrate adaptive radiation to investigate the genetic basis of the response to selection for predatory behavior and compare it with evolution of aerobic capacity reported in an earlier work. After 13 generations of selection, the proportion of bank voles (Myodes [=Clethrionomys] glareolus) showing predatory behavior was five times higher in selected lines than in controls. We analyzed the hippocampus and liver transcriptomes and found repeatable changes in allele frequencies and gene expression. Genes with the largest differences between predatory and control lines are associated with hunger, aggression, biological rhythms, and functioning of the nervous system. Evolution of predatory behavior could be meaningfully compared with evolution of high aerobic capacity, because the experiments and analyses were performed in the same methodological framework. The number of genes that changed expression was much smaller in predatory lines, and allele frequencies changed repeatably in predatory but not in aerobic lines. This suggests that more variants of smaller effects underlie variation in aerobic performance, whereas fewer variants of larger effects underlie variation in predatory behavior. Our results thus contradict the view that comparable selection pressures for different quantitative traits produce similar molecular patterns. Therefore, to gain knowledge about molecular-level response to selection for complex traits, we need to investigate not only multiple replicate populations but also multiple quantitative traits.

Evolution of stickleback in 50 years on earthquake-uplifted islands

How rapidly can animal populations in the wild evolve when faced with sudden environmental shifts? Uplift during the 1964 Great Alaska Earthquake abruptly created freshwater ponds on multiple islands in Prince William Sound and the Gulf of Alaska. In the short time since the earthquake, the phenotypes of resident freshwater threespine stickleback fish on at least three of these islands have changed dramatically from their oceanic ancestors. To test the hypothesis that these freshwater populations were derived from oceanic ancestors only 50 y ago, we generated over 130,000 single-nucleotide polymorphism genotypes from more than 1,000 individuals using restriction site-associated DNA sequencing (RAD-seq). Population genomic analyses of these data support the hypothesis of recent and repeated, independent colonization of freshwater habitats by oceanic ancestors. We find evidence of recurrent gene flow between oceanic and freshwater ecotypes where they co-occur. Our data implicate natural selection in phenotypic diversification and support the hypothesis that the metapopulation organization of this species helps maintain a large pool of genetic variation that can be redeployed rapidly when oceanic stickleback colonize freshwater environments. We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.

Oviposition but Not Sex Allocation Is Associated with Transcriptomic Changes in Females of the Parasitoid Wasp Nasonia vitripennis

Linking the evolution of the phenotype to the underlying genotype is a key aim of evolutionary genetics and is crucial to our understanding of how natural selection shapes a trait. Here, we consider the genetic basis of sex allocation behavior in the parasitoid wasp Nasonia vitripennis using a transcriptomics approach. Females allocate offspring sex in line with the local mate competition (LMC) theory. Female-biased sex ratios are produced when one or a few females lay eggs on a patch. As the number of females contributing offspring to a patch increases, less female-biased sex ratios are favored. We contrasted the transcriptomic responses of females as they oviposit under conditions known to influence sex allocation: foundress number (a social cue) and the state of the host (parasitized or not). We found that when females encounter other females on a patch or assess host quality with their ovipositors, the resulting changes in sex allocation is not associated with significant changes in whole-body gene expression. We also found that the gene expression changes produced by females as they facultatively allocate sex in response to a host cue and a social cue are very closely correlated. We expanded the list of candidate genes associated with oviposition behavior in Nasonia, some of which may be involved in fundamental processes underlying the ability to facultatively allocate sex, including sperm storage and utilization.

Signatures of natural selection between life cycle stages separated by metamorphosis in European eel

BACKGROUND

Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels.

RESULTS

We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle.

CONCLUSIONS

The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Using genomics to characterize evolutionary potential for conservation of wild populations

Genomics promises exciting advances towards the important conservation goal of maximizing evolutionary potential, notwithstanding associated challenges. Here, we explore some of the complexity of adaptation genetics and discuss the strengths and limitations of genomics as a tool for characterizing evolutionary potential in the context of conservation management. Many traits are polygenic and can be strongly influenced by minor differences in regulatory networks and by epigenetic variation not visible in DNA sequence. Much of this critical complexity is difficult to detect using methods commonly used to identify adaptive variation, and this needs appropriate consideration when planning genomic screens, and when basing management decisions on genomic data. When the genomic basis of adaptation and future threats are well understood, it may be appropriate to focus management on particular adaptive traits. For more typical conservations scenarios, we argue that screening genome-wide variation should be a sensible approach that may provide a generalized measure of evolutionary potential that accounts for the contributions of small-effect loci and cryptic variation and is robust to uncertainty about future change and required adaptive response(s). The best conservation outcomes should be achieved when genomic estimates of evolutionary potential are used within an adaptive management framework.

Fast evolution from precast bricks: genomics of young freshwater populations of threespine stickleback Gasterosteus aculeatus

Adaptation is driven by natural selection; however, many adaptations are caused by weak selection acting over large timescales, complicating its study. Therefore, it is rarely possible to study selection comprehensively in natural environments. The threespine stickleback (Gasterosteus aculeatus) is a well-studied model organism with a short generation time, small genome size, and many genetic and genomic tools available. Within this originally marine species, populations have recurrently adapted to freshwater all over its range. This evolution involved extensive parallelism: pre-existing alleles that adapt sticklebacks to freshwater habitats, but are also present at low frequencies in marine populations, have been recruited repeatedly. While a number of genomic regions responsible for this adaptation have been identified, the details of selection remain poorly understood. Using whole-genome resequencing, we compare pooled genomic samples from marine and freshwater populations of the White Sea basin, and identify 19 short genomic regions that are highly divergent between them, including three known inversions. 17 of these regions overlap protein-coding genes, including a number of genes with predicted functions that are relevant for adaptation to the freshwater environment. We then analyze four additional independently derived young freshwater populations of known ages, two natural and two artificially established, and use the observed shifts of allelic frequencies to estimate the strength of positive selection. Adaptation turns out to be quite rapid, indicating strong selection acting simultaneously at multiple regions of the genome, with selection coefficients of up to 0.27. High divergence between marine and freshwater genotypes, lack of reduction in polymorphism in regions responsible for adaptation, and high frequencies of freshwater alleles observed even in young freshwater populations are all consistent with rapid assembly of G. aculeatus freshwater genotypes from pre-existing genomic regions of adaptive variation, with strong selection that favors this assembly acting simultaneously at multiple loci.

Adaptation to novel environments is a keystone of evolution. There is only a handful of natural and experimental systems in which the process of adaptation has been studied in detail, and each studied system brings its own surprises with regard to the number of loci involved, dynamics of adaptation, extent of interactions between loci and of parallelism between different adapting populations. The threespine stickleback is an excellent model organism for evolutionary studies. Marine-derived freshwater populations of this species have consistently acquired a specific set of morphological, physiological and behavioral traits allowing them to reside in freshwater for their whole lifespan. Previous studies identified several genomic regions responsible for this adaptation. Here, using whole-genome sequencing, we compare the allele frequencies at such regions in four derived freshwater populations of known ages: two natural, and two artificially established in 1978. Knowledge of population ages allows us to infer the strength of selection that acted at these loci. Adaptation of threespine stickleback to freshwater is typically fast, and is driven by strong selection favoring pre-existing alleles that are likely present in the ancestral marine population at low frequencies; however, some of the adaptation may also be due to young population-specific alleles.

Population genomics of marine fishes: next-generation prospects and challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

Dynamics of copy number variation in host races of the pea aphid

Copy number variation (CNV) makes a major contribution to overall genetic variation and is suspected to play an important role in adaptation. However, aside from a few model species, the extent of CNV in natural populations has seldom been investigated. Here, we report on CNV in the pea aphid Acyrthosiphon pisum, a powerful system for studying the genetic architecture of host-plant adaptation and speciation thanks to multiple host races forming a continuum of genetic divergence. Recent studies have highlighted the potential importance of chemosensory genes, including the gustatory and olfactory receptor gene families (Gr and Or, respectively), in the process of host race formation. We used targeted resequencing to achieve a very high depth of coverage, and thereby revealed the extent of CNV of 434 genes, including 150 chemosensory genes, in 104 individuals distributed across eight host races of the pea aphid. We found that CNV was widespread in our global sample, with a significantly higher occurrence in multigene families, especially in Ors. We also observed a decrease in the gene probability of being completely duplicated or deleted (CDD) with increase in coding sequence length. Genes with CDD variants were usually more polymorphic for copy number, especially in the P450 gene family where toxin resistance may be related to gene dosage. We found that Gr were overrepresented among genes discriminating host races, as were CDD genes and pseudogenes. Our observations shed new light on CNV dynamics and are consistent with CNV playing a role in both local adaptation and speciation.

Genome-wide single-generation signatures of local selection in the panmictic European eel

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Sympatric incipient speciation of spiny mice Acomys at "Evolution Canyon," Israel

Does the paucity of empirical evidence of sympatric speciation in nature reflect reality, despite theoretical support? Or is it due to inappropriate searches in nature with overly restrictive assumptions and an incorrect null hypothesis? Spiny mice, Acomys, described here at Evolution Canyon (EC) incipiently and sympatrically speciate owing to microclimatic interslope divergence. The opposite slopes at EC vary dramatically, physically and biotically, representing the dry and hot south-facing slope savannoid-African continent ["African" slope (AS)], abutting with the north-facing slope forested south-European continent ["European" slope (ES)]. African-originated spiny mice, of the Acomys cahirinus complex, colonized Israel 30,000 y ago based on fossils. Genotypically, we showed significantly higher genetic diversity of mtDNA and amplified fragment length polymorphism of Acomys on the AS compared with the ES. This is also true regionally across Israel. In complete mtDNA, 25% of the haplotypes at EC were slope-biased. Phenotypically, the opposite slope's populations also showed adaptive morphology, physiology, and behavior divergence paralleling regional populations across Israel. Preliminary tests indicate slope-specific mate choices. Colonization of Acomys at the EC first occurred on the AS and then moved to the ES. Strong slope-specific natural selection (both positive and negative) overrules low interslope gene flow. Both habitat slope selection and mate choices suggest ongoing incipient sympatric speciation. We conclude that Acomys at the EC is ecologically and genetically adaptively, incipiently, sympatrically speciating on the ES owing to adaptive microclimatic natural selection.

Differential strengths of positive selection revealed by hitchhiking effects at small physical scales in Drosophila melanogaster

The long time scale of adaptive evolution makes it difficult to directly observe the spread of most beneficial mutations through natural populations. Therefore, inferring attributes of beneficial mutations by studying the genomic signals left by directional selection is an important component of population genetics research. One kind of signal is a trough in nearby neutral genetic variation due to selective fixation of initially rare alleles, a phenomenon known as "genetic hitchhiking." Accumulated evidence suggests that a considerable fraction of substitutions in the Drosophila genome results from positive selection, most of which are expected to have small selection coefficients and influence the population genetics of sites in the immediate vicinity. Using Drosophila melanogaster population genomic data, we found that the heterogeneity in synonymous polymorphism surrounding different categories of coding fixations is readily observable even within 25 bp of focal substitutions, which we interpret as the result of small-scale hitchhiking effects. The strength of natural selection on different sites appears to be quite heterogeneous. Particularly, neighboring fixations that changed amino acid polarities in a way that maintained the overall polarities of a protein were under stronger selection than other categories of fixations. Interestingly, we found that substitutions in slow-evolving genes are associated with stronger hitchhiking effects. This is consistent with the idea that adaptive evolution may involve few substitutions with large effects or many substitutions with small effects. Because our approach only weakly depends on the numbers of recent nonsynonymous substitutions, it can provide a complimentary view to the adaptive evolution inferred by other divergence-based evolutionary genetic methods.

Evolution in an acidifying ocean

Ocean acidification poses a global threat to biodiversity, yet species might have the capacity to adapt through evolutionary change. Here we summarize tools available to determine species' capacity for evolutionary adaptation to future ocean change and review the progress made to date with respect to ocean acidification. We focus on two key approaches: measuring standing genetic variation within populations and experimental evolution. We highlight benefits and challenges of each approach and recommend future research directions for understanding the modulating role of evolution in a changing ocean.

Accuracy of allele frequency estimation using pooled RNA-Seq

For nonmodel organisms, genome-wide information that describes functionally relevant variation may be obtained by RNA-Seq following de novo transcriptome assembly. While sequencing has become relatively inexpensive, the preparation of a large number of sequencing libraries remains prohibitively expensive for population genetic analyses of nonmodel species. Pooling samples may be then an attractive alternative. To test whether pooled RNA-Seq accurately predicts true allele frequencies, we analysed the liver transcriptomes of 10 bank voles. Each sample was sequenced both as an individually barcoded library and as a part of a pool. Equal amounts of total RNA from each vole were pooled prior to mRNA selection and library construction. Reads were mapped onto the de novo assembled reference transcriptome. High-quality genotypes for individual voles, determined for 23,682 SNPs, provided information on 'true' allele frequencies; allele frequencies estimated from the pool were then compared with these values. 'True' frequencies and those estimated from the pool were highly correlated. Mean relative estimation error was 21% and did not depend on expression level. However, we also observed a minor effect of interindividual variation in gene expression and allele-specific gene expression influencing allele frequency estimation accuracy. Moreover, we observed strong negative relationship between minor allele frequency and relative estimation error. Our results indicate that pooled RNA-Seq exhibits accuracy comparable with pooled genome resequencing, but variation in expression level between individuals should be assessed and accounted for. This should help in taking account the difference in accuracy between conservatively expressed transcripts and these which are variable in expression level.

Molecular genetics and genomics generate new insights into invertebrate pest invasions

Invertebrate pest invasions and outbreaks are associated with high social, economic, and ecological costs, and their significance will intensify with an increasing pressure on agricultural productivity as a result of human population growth and climate change. New molecular genetic and genomic techniques are available and accessible, but have been grossly underutilized in studies of invertebrate pest invasions, despite that they are useful tools for applied pest management and for understanding fundamental features of pest invasions including pest population demographics and adaptation of pests to novel and/or changing environments. Here, we review current applications of molecular genetics and genomics in the study of invertebrate pest invasions and outbreaks, and we highlight shortcomings from the current body of research. We then discuss recent conceptual and methodological advances in the areas of molecular genetics/genomics and data analysis, and we highlight how these advances will further our understanding of the demographic, ecological, and evolutionary features of invertebrate pest invasions. We are now well equipped to use molecular data to understand invertebrate dispersal and adaptation, and this knowledge has valuable applications in agriculture at a time when these are critically required.