Signatures of local adaptation along environmental gradients in a range-expanding damselfly (Ischnura elegans)

Reproductive character displacement: insights from genital morphometrics in damselfly hybrid zones

Reproductive Character Displacement (RCD) refers to the phenomenon of greater differences in reproductive characters between two species when they occur in sympatry compared to when they occur in allopatry to prevent maladaptive hybridization. We explored whether reinforcement of a mechanical barrier involved in the first contact point between male and female genital traits during copulation in the cross between Ischnura graellsii males and Ischnura elegans females has led to RCD, and whether it supports the lock-and-key hypothesis of genital evolution. We employed geometric morphometrics to analyze the shape and size of male and female genital traits, controlling for environmental and geographic factors. Consistent with an increase in mechanical isolation via reinforcement, we detected larger divergence in genital traits between the species in sympatry than in allopatry, and also stronger signal in females than in males. In the Northwest (NW) hybrid zone, we detected RCD in I. graellsii males and I. elegans females, while in the Northcentral (NC) hybrid zone we detected RCD only in I. elegans females and I. elegans males. The detection of RCD in both sexes of I. elegans was consistent with the lock-and-key hypothesis of genital evolution via female choice for conspecific males in this species. Our study highlights the importance of using geometric morphometrics to deal with the complexity of female reproductive structures while controlling for environmental and geographic factors to investigate RCD. This study contributes valuable insights into the dynamics of reproductive isolation mechanisms and genital coevolution.

Transposable elements in Drosophila montana from harsh cold environments

BACKGROUND

Substantial discoveries during the past century have revealed that transposable elements (TEs) can play a crucial role in genome evolution by affecting gene expression and inducing genetic rearrangements, among other molecular and structural effects. Yet, our knowledge on the role of TEs in adaptation to extreme climates is still at its infancy. The availability of long-read sequencing has opened up the possibility to identify and study potential functional effects of TEs with higher precision. In this work, we used Drosophila montana as a model for cold-adapted organisms to study the association between TEs and adaptation to harsh climates.

RESULTS

Using the PacBio long-read sequencing technique, we de novo identified and manually curated TE sequences in five Drosophila montana genomes from eco-geographically distinct populations. We identified 489 new TE consensus sequences which represented 92% of the total TE consensus in D. montana. Overall, 11-13% of the D. montana genome is occupied by TEs, which as expected are non-randomly distributed across the genome. We identified five potentially active TE families, most of them from the retrotransposon class of TEs. Additionally, we found TEs present in the five analyzed genomes that were located nearby previously identified cold tolerant genes. Some of these TEs contain promoter elements and transcription binding sites. Finally, we detected TEs nearby fixed and polymorphic inversion breakpoints.

CONCLUSIONS

Our research revealed a significant number of newly identified TE consensus sequences in the genome of D. montana, suggesting that non-model species should be studied to get a comprehensive view of the TE repertoire in Drosophila species and beyond. Genome annotations with the new D. montana library allowed us to identify TEs located nearby cold tolerant genes, and present at high population frequencies, that contain regulatory regions and are thus good candidates to play a role in D. montana cold stress response. Finally, our annotations also allow us to identify for the first time TEs present in the breakpoints of three D. montana inversions.

Local adaptation of Dromiciops marsupials (Microbiotheriidae) from southern South America: Implications for species management facing climate change

The two species of the microbiotheriid marsupial genus Dromiciops (Dromiciops bozinovici: "Panchos's monito del monte" and Dromiciops gliroides: "monito del monte") exhibit a marked latitudinal genetic differentiation. Nevertheless, it is unclear whether this differentiation results from neutral processes or can be explained, to some extent, by local adaptation to different environmental conditions. Here, we used an SNP panel gathered by Rad-seq and searched for footprints of local adaptation (putative loci under selection) by exploring genetic associations with environmental variables in the two species of Dromiciops in Chilean and Argentinean populations. We applied three methods for detecting outlier SNPs and two genotype-environment associations approaches to quantify associations between allelic frequencies and environmental variables. Both species display strong genetic structure. D. bozinovici exhibited three distinct genetic groups, marking the first report of such structuring in this species using SNPs. In contrast, D. gliroides displayed four genetic clusters, consistent with previous studies. Both species exhibited an association of their genetic structure with environmental variables. D. bozinovici exhibited significant associations of allelic frequencies with elevation, precipitation during the warmest periods, and seasonality in the thermal regime. For D. gliroides, genetic variation appeared to be associated with more variables than D. bozinovici, including precipitation and temperature-related variables, isothermality, and elevation. All the outlier SNPs were mapped to the D. gliroides reference genome to explore if they fell within functionally known genes. These results represent a necessary first step toward identifying the genome regions that harbor genes associated with climate adaptations in Dromiciops. Notably, we identified genes involved in various functions, including carbohydrate synthesis (ALG8), muscle and neuronal regulation (MEF2D), and stress responses (PTGES3). Ultimately, this study contributes valuable insights that can inform targeted conservation strategies aimed at preserving the genetic diversity of Dromiciops in the face of environmental challenges.

Assessing genome-wide adaptations associated with range expansion in the pink rice borer, Sesamia inferens

Understanding the genetic basis of adaptive evolution following habitat expansion can have important implications for pest management. The pink rice borer (PRB), Sesamia inferens (Walker), is a destructive pest of rice that was historically restricted to regions south of 34° N latitude in China. However, with changes in global climate and farming practices, the distribution of this moth has progressively expanded, encompassing most regions in North China. Here, 3 highly differentiated subpopulations were discovered using high-quality single-nucleotide polymorphism and structural variant datasets across China, corresponding to northern, southern China regions, and the Yunnan-Guizhou Plateau, with significant patterns of isolation by geographic and environmental distances. Our estimates of evolutionary history indicate asymmetric migration with varying population sizes across the 3 subpopulations. Selective sweep analyses estimated strong selection at insect cuticle glycine-rich cuticular protein genes which are associated with enhanced desiccation adaptability in the northern group, and at the histone-lysine-N-methyltransferase gene associated with range expansion and local adaptation in the Shandong population. Our findings have significant implications for the development of effective strategies to control this pest.

Selection despite low genetic diversity and high gene flow in a rapid island invasion of the bumblebee, Bombus terrestris

Invasive species are predicted to adjust their morphological, physiological and life-history traits to adapt to their non-native environments. Although a loss of genetic variation during invasion may restrict local adaptation, introduced species often thrive in novel environments. Despite being founded by just a few individuals, Bombus terrestris (Hymenoptera: Apidae) has in less than 30 years successfully spread across the island of Tasmania (Australia), becoming abundant and competitive with native pollinators. We use RADseq to investigate what neutral and adaptive genetic processes associated with environmental and morphological variation allow B. terrestris to thrive as an invasive species in Tasmania. Given the widespread abundance of B. terrestris, we expected little genetic structure across Tasmania and weak signatures of environmental and morphological selection. We found high gene flow with low genetic diversity, although with significant isolation-by-distance and spatial variation in effective migration rates. Restricted migration was evident across the mid-central region of Tasmania, corresponding to higher elevations, pastural land, low wind speeds and low precipitation seasonality. Tajima's D indicated a recent population expansion extending from the south to the north of the island. Selection signatures were found for loci in relation to precipitation, wind speed and wing loading. Candidate loci were annotated to genes with functions related to cuticle water retention and insect flight muscle stability. Understanding how a genetically impoverished invasive bumblebee has rapidly adapted to a novel island environment provides further understanding about the evolutionary processes that determine successful insect invasions, and the potential for invasive hymenopteran pollinators to spread globally.

Genomic differentiation and niche divergence in the Hetaerina americana (Odonata) cryptic species complex

The evolution of reproductive barriers, that is, the speciation process, implies the limitation of gene flow between populations. Different patterns of genomic differentiation throughout the speciation continuum may provide insights into the causal evolutionary forces of species divergence. In this study, we analysed a cryptic species complex of the genus Hetaerina (Odonata). This complex includes H. americana and H. calverti; however, in H. americana two highly differentiated genetic groups have been previously detected, which, we hypothesize, may correspond to different species with low morphological variation. We obtained single nucleotide polymorphism (SNP) data for 90 individuals belonging to the different taxa in the complex and carried out differentiation tests to identify genetic isolation. The results from STRUCTURE and discriminant analysis of principal components (DAPC), based on almost 5000 SNPs, confirmed the presence of three highly differentiated taxa. Also, we found FST values above 0.5 in pairwise comparisons, which indicates a considerable degree of genetic isolation among the suggested species. We also found low climatic niche overlap among all taxa, suggesting that each group occurs at specific conditions of temperature, precipitation and elevation. We propose that H. americana comprises two cryptic species, which may be reproductively isolated by ecological barriers related to niche divergence, since the morphological variation is minimal and, therefore, mechanical barriers are probably less effective compared to other related species such as H. calverti.

The genomic response to urbanization in the damselfly Ischnura elegans

The complex and rapid environmental changes brought about by urbanization pose significant challenges to organisms. The multifaceted effects of urbanization often make it difficult to define and pinpoint the very nature of adaptive urban phenotypes. In such situations, scanning genomes for regions differentiated between urban and non-urban populations may be an attractive approach. Here, we investigated the genomic signatures of adaptation to urbanization in the damselfly Ischnura elegans sampled from 31 rural and urban localities in three geographic regions: southern and northern Poland, and southern Sweden. Genome-wide variation was assessed using more than 370,000 single nucleotide polymorphisms (SNPs) genotyped by ddRADseq. Associations between SNPs and the level of urbanization were tested using two genetic environment association methods: Latent Factors Mixed Models and BayPass. While we found numerous candidate SNPs and a highly significant overlap between candidates identified by the two methods within the geographic regions, there was a distinctive lack of repeatability between the geographic regions both at the level of individual SNPs and of genomic regions. However, we found "synapse organization" at the top of the functional categories enriched among the genes located in the proximity of the candidate urbanization SNPs. Interestingly, the overall significance of "synapse organization" was built up by the accretion of different genes associated with candidate SNPs in different geographic regions. This finding is consistent with the highly polygenic nature of adaptation, where the response may be achieved through a subtle adjustment of allele frequencies in different genes that contribute to adaptive phenotypes. Taken together, our results point to a polygenic adaptive response in the nervous system, specifically implicating genes involved in synapse organization, which mirrors the findings from several genomic and behavioral studies of adaptation to urbanization in other taxa.

Statistical methods to identify mechanisms in studies of eco-evolutionary dynamics

While the reciprocal effects of ecological and evolutionary dynamics are increasingly recognized as an important driver for biodiversity, detection of such eco-evolutionary feedbacks, their underlying mechanisms, and their consequences remains challenging. Eco-evolutionary dynamics occur at different spatial and temporal scales and can leave signatures at different levels of organization (e.g., gene, protein, trait, community) that are often difficult to detect. Recent advances in statistical methods combined with alternative hypothesis testing provides a promising approach to identify potential eco-evolutionary drivers for observed data even in non-model systems that are not amenable to experimental manipulation. We discuss recent advances in eco-evolutionary modeling and statistical methods and discuss challenges for fitting mechanistic models to eco-evolutionary data.

Global patterns of genomic and phenotypic variation in the invasive harlequin ladybird

BACKGROUND

The harlequin ladybird Harmonia axyridis (Coleoptera: Coccinellidae), native to Asia, has been introduced to other major continents where it has caused serious negative impacts on local biodiversity. Though notable advances to understand its invasion success have been made during the past decade, especially with then newer molecular tools, the conclusions reached remain to be confirmed with more advanced genomic analyses and especially using more samples from larger geographical regions across the native range. Furthermore, although H. axyridis is one of the best studied invasive insect species with respect to life history traits (often comparing invasive and native populations), the traits responsible for its colonization success in non-native areas warrant more research.

RESULTS

Our analyses of genome-wide nuclear population structure indicated that an eastern Chinese population could be the source of all non-native populations and revealed several putatively adaptive candidate genomic loci involved in body color variation, visual perception, and hemolymph synthesis. Our estimates of evolutionary history indicate (1) asymmetric migration with varying population sizes across its native and non-native range, (2) a recent admixture between eastern Chinese and American populations in Europe, (3) signatures of a large progressive, historical bottleneck in the common ancestors of both populations and smaller effective sizes of the non-native population, and (4) the southwest origin and subsequent dispersal routes within its native range in China. In addition, we found that while two mitochondrial haplotypes-Hap1 and Hap2 were dominant in the native range, Hap1 was the only dominant haplotype in the non-native range. Our laboratory observations in both China and USA found statistical yet slight differences between Hap1 and Hap2 in some of life history traits.

CONCLUSIONS

Our study on H. axyridis provides new insights into its invasion processes into other major continents from its native Asian range, reconstructs a geographic range evolution across its native region China, and tentatively suggests that its invasiveness may differ between mitochondrial haplotypes.

The genomic and epigenetic footprint of local adaptation to variable climates in kiwifruit

A full understanding of adaptive genetic variation at the genomic level will help address questions of how organisms adapt to diverse climates. Actinidia eriantha is a shade-tolerant species, widely distributed in the southern tropical region of China, occurring in spatially heterogeneous environments. In the present study we combined population genomic, epigenomic, and environmental association analyses to infer population genetic structure and positive selection across a climatic gradient, and to assess genomic offset to climatic change for A. eriantha. The population structure is strongly shaped by geography and influenced by restricted gene flow resulting from isolation by distance due to habitat fragmentation. In total, we identified 102 outlier loci and annotated 455 candidate genes associated with the genomic basis of climate adaptation, which were enriched in functional categories related to development processes and stress response; both temperature and precipitation are important factors driving adaptive variation. In addition to single-nucleotide polymorphisms (SNPs), a total of 27 single-methylation variants (SMVs) had significant correlation with at least one of four climatic variables and 16 SMVs were located in or adjacent to genes, several of which were predicted to be involved in plant response to abiotic or biotic stress. Gradient forest analysis indicated that the central/east populations were predicted to be at higher risk of future population maladaptation under climate change. Our results demonstrate that local climate factors impose strong selection pressures and lead to local adaptation. Such information adds to our understanding of adaptive mechanisms to variable climates revealed by both population genome and epigenome analysis.

Rapid introgression of non-native alleles following hybridization between a native Anolis lizard species and a cryptic invader across an urban landscape

Invasive species can impact native populations through competition, predation, habitat alteration, and disease transmission, but also genetically through hybridization. Potential outcomes of hybridization span the continuum from extinction to hybrid speciation and can be further complicated by anthropogenic habitat disturbance. Hybridization between the native green anole lizard (Anolis carolinensis) and a morphologically similar invader (A. porcatus) in south Florida provides an ideal opportunity to study interspecific admixture across a heterogeneous landscape. We used reduced-representation sequencing to describe introgression in this hybrid system and to test for a relationship between urbanization and non-native ancestry. Our findings indicate that hybridization between green anole lineages was probably a limited, historic event, producing a hybrid population characterized by a diverse continuum of ancestry proportions. Genomic cline analyses revealed rapid introgression and disproportionate representation of non-native alleles at many loci and no evidence for reproductive isolation between parental species. Three loci were associated with urban habitat characteristics; urbanization and non-native ancestry were positively correlated, although this relationship did not remain significant when accounting for spatial nonindependence. Ultimately, our study demonstrates the persistence of non-native genetic material even in the absence of ongoing immigration, indicating that selection favouring non-native alleles can override the demographic limitation of low propagule pressure. We also note that not all outcomes of admixture between native and non-native species should be considered intrinsically negative. Hybridization with ecologically robust invaders can lead to adaptive introgression, which may facilitate the long-term survival of native populations otherwise unable to adapt to anthropogenically mediated global change.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

ChIP-seq profiling of H3K4me3 and H3K27me3 in an invasive insect, Bactrocera dorsalis

Introduction: While it has been suggested that histone modifications can facilitate animal responses to rapidly changing environments, few studies have profiled whole-genome histone modification patterns in invasive species, leaving the regulatory landscape of histone modifications in invasive species unclear. Methods: Here, we screen genome-wide patterns of two important histone modifications, trimethylated Histone H3 Lysine 4 (H3K4me3) and trimethylated Histone H3 Lysine 27 (H3K27me3), in adult thorax muscles of a notorious invasive pest, the Oriental fruit fly Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), using Chromatin Immunoprecipitation with high-throughput sequencing (ChIP-seq). Results: We identified promoters featured by the occupancy of H3K4me3, H3K27me3 or bivalent histone modifications that were respectively annotated with unique genes key to muscle development and structure maintenance. In addition, we found H3K27me3 occupied the entire body of genes, where the average enrichment was almost constant. Transcriptomic analysis indicated that H3K4me3 is associated with active gene transcription, and H3K27me3 is mostly associated with transcriptional repression. Importantly, we identified genes and putative motifs modified by distinct histone modification patterns that may possibly regulate flight activity. Discussion: These findings provide the first evidence of histone modification signature in B. dorsalis, and will be useful for future studies of epigenetic signature in other invasive insect species.

Using landscape genomics to assess local adaptation and genomic vulnerability of a perennial herb Tetrastigma hemsleyanum (Vitaceae) in subtropical China

Understanding adaptive genetic variation of plant populations and their vulnerabilities to climate change are critical to preserve biodiversity and subsequent management interventions. To this end, landscape genomics may represent a cost-efficient approach for investigating molecular signatures underlying local adaptation. Tetrastigma hemsleyanum is, in its native habitat, a widespread perennial herb of warm-temperate evergreen forest in subtropical China. Its ecological and medicinal values constitute a significant revenue for local human populations and ecosystem. Using 30,252 single nucleotide polymorphisms (SNPs) derived from reduced-representation genome sequencing in 156 samples from 24 sites, we conducted a landscape genomics study of the T. hemsleyanum to elucidate its genomic variation across multiple climate gradients and genomic vulnerability to future climate change. Multivariate methods identified that climatic variation explained more genomic variation than that of geographical distance, which implied that local adaptation to heterogeneous environment might represent an important source of genomic variation. Among these climate variables, winter precipitation was the strongest predictor of the contemporary genetic structure. F _ST outlier tests and environment association analysis totally identified 275 candidate adaptive SNPs along the genetic and environmental gradients. SNP annotations of these putatively adaptive loci uncovered gene functions associated with modulating flowering time and regulating plant response to abiotic stresses, which have implications for breeding and other special agricultural aims on the basis of these selection signatures. Critically, modelling revealed that the high genomic vulnerability of our focal species via a mismatch between current and future genotype-environment relationships located in central-northern region of the T. hemsleyanum's range, where populations require proactive management efforts such as assistant adaptation to cope with ongoing climate change. Taken together, our results provide robust evidence of local climate adaption for T. hemsleyanum and further deepen our understanding of adaptation basis of herbs in subtropical China.

Phylogenomic analysis of the bowfin (Amia calva) reveals unrecognized species diversity in a living fossil lineage

The Bowfin (Amia calva), as currently recognized, represents the sole living member of the family Amiidae, which dates back to approximately 150 Ma. Prior to 1896, 13 species of extant Bowfins had been described, but these were all placed into a single species with no rationale or analysis given. This situation has persisted until the present day, with little attention given to re-evaluation of those previously described nominal forms. Here, we present a phylogenomic analysis based on over 21,000 single nucleotide polymorphisms (SNPs) from 94 individuals that unambiguously demonstrates the presence of at least two independent evolutionary lineages within extant Amia populations that merit species-level standing, as well as the possibility of two more. These findings not only expand the recognizable species diversity in an iconic, ancient lineage, but also demonstrate the utility of such methods in addressing previously intractable questions of molecular systematics and phylogeography in slowly evolving groups of ancient fishes.

Parallel adaptation prompted core-periphery divergence of Ammopiptanthus mongolicus

Range expansion requires peripheral populations to shift adaptive optima to breach range boundaries. Opportunities for range expansion can be assessed by investigating the associations of core-periphery environmental and genetic differences. This study investigates differences in the core-periphery adaptation of Ammopiptanthus mongolicus, a broad-leaved evergreen shrub species in a relatively homogeneous temperate Asian desert environment, to explore the environmental factors that limit the expansion of desert plants. Temperate deserts are characterized by severe drought, a large diurnal temperature range, and seasonality. Long-standing adaptation to the harsh desert environment may confine the genetic diversity of A. mongolicus, despite its distribution over a wide range of longitude, latitude, and altitude. Since range edges defined by climate niches may have different genetic responses to environmental extremes, we compared genome-wide polymorphisms between nine environmental core populations and ten fragmented peripheral populations to determine the "adaptive peripheral" populations. At least four adaptive peripheral populations had similar genetic-environmental association patterns. High elevations, summer drought, and winter cold were the three main determinants of converging these four adaptive peripheral populations. Elevation mainly caused similar local climates among different geographic regions. Altitudinal adaptation resulting from integrated environmental-genetic responses was a breakthrough in breaching niche boundaries. These peripheral populations are also located in relatively humid and warmer environments. Relaxation of the drought and cold constraints facilitated the genetic divergence of these peripheral populations from the core population's adaptive legacy. We conclude that pleiotropic selection synchronized adaptative divergence to cold and drought vs. warm and humid environments between the core and peripheral populations. Such parallel adaptation of peripheral populations relies on selection under a background of abundant new variants derived from the core population's standing genetic variation, i.e., integration of genetic surfing and local adaptation.

The importance of eco-evolutionary dynamics for predicting and managing insect range shifts

Evolutionary change impacts the rate at which insect pests, pollinators, or disease vectors expand or contract their geographic ranges. Although evolutionary changes, and their ecological feedbacks, strongly affect these risks and associated ecological and economic consequences, they are often underappreciated in management efforts. Greater rigor and scope in study design, coupled with innovative technologies and approaches, facilitates our understanding of the causes and consequences of eco-evolutionary dynamics in insect range shifts. Future efforts need to ensure that forecasts allow for demographic and evolutionary change and that management strategies will maximize (or minimize) the adaptive potential of range-shifting insects, with benefits for biodiversity and ecosystem services.

Restricted X chromosome introgression and support for Haldane's rule in hybridizing damselflies

Contemporary hybrid zones act as natural laboratories for the investigation of species boundaries and may shed light on the little understood roles of sex chromosomes in species divergence. Sex chromosomes are considered to function as a hotspot of genetic divergence between species; indicated by less genomic introgression compared to autosomes during hybridization. Moreover, they are thought to contribute to Haldane's rule, which states that hybrids of the heterogametic sex are more likely to be inviable or sterile. To test these hypotheses, we used contemporary hybrid zones of Ischnura elegans, a damselfly species that has been expanding its range into the northern and western regions of Spain, leading to chronic hybridization with its sister species Ischnura graellsii. We analysed genome-wide SNPs in the Spanish I. elegans and I. graellsii hybrid zone and found (i) that the X chromosome shows less genomic introgression compared to autosomes, and (ii) that males are underrepresented among admixed individuals, as predicted by Haldane's rule. This is the first study in Odonata that suggests a role of the X chromosome in reproductive isolation. Moreover, our data add to the few studies on species with X0 sex determination system and contradict the hypothesis that the absence of a Y chromosome causes exceptions to Haldane's rule.

Temperature-dependent genetics of thermotolerance between yeast species

Many traits of industrial and basic biological interest arose long ago, and manifest now as fixed differences between a focal species and its reproductively isolated relatives. In these systems, extant individuals can hold clues to the mechanisms by which phenotypes evolved in their ancestors. We harnessed yeast thermotolerance as a test case for such molecular-genetic inferences. In viability experiments, we showed that extant Saccharomyces cerevisiae survived at temperatures where cultures of its sister species S. paradoxus died out. Then, focusing on loci that contribute to this difference, we found that the genetic mechanisms of high-temperature growth changed with temperature. We also uncovered an enrichment of low-frequency variants at thermotolerance loci in S. cerevisiae population sequences, suggestive of a history of non-neutral selective forces acting at these genes. We interpret these results in light of models of the evolutionary mechanisms by which the thermotolerance trait arose in the S. cerevisiae lineage. Together, our results and interpretation underscore the power of genetic approaches to explore how an ancient trait came to be.

Genomic associations with poxvirus across divergent island populations in Berthelot's pipit

Understanding the mechanisms and genes that enable animal populations to adapt to pathogens is important from an evolutionary, health and conservation perspective. Berthelot's pipit (Anthus berthelotii) experiences extensive and consistent spatial heterogeneity in avian pox infection pressure across its range of island populations, thus providing an excellent system with which to examine how pathogen-mediated selection drives spatial variation in immunogenetic diversity. Here we test for evidence of genetic variation associated with avian pox at both an individual and population-level. At the individual level, we find no evidence that variation in MHC class I and TLR4 (both known to be important in recognising viral infection) was associated with pox infection within two separate populations. However, using genotype-environment association (Bayenv) in conjunction with genome-wide (ddRAD-seq) data, we detected strong associations between population-level avian pox prevalence and allele frequencies of single nucleotide polymorphisms (SNPs) at a number of sites across the genome. These sites were located within genes involved in cellular stress signalling and immune responses, many of which have previously been associated with responses to viral infection in humans and other animals. Consequently, our analyses indicates that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlights the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.

Genetic divergence and range expansion in a western North Pacific coral

Coral poleward range expansions have recently been observed in response to warming oceans. Range expansion can lead to reduced genetic diversity and increased frequency of deleterious mutations that were rare in core populations, potentially limiting the ability for adaptation and persistence in novel environments. Successful expansions that overcome these founder effects and colonize new habitat have been attributed to multiple introductions from different sources, hybridization with native populations, or rapid adaptive evolution. Here, we investigate population genomic patterns of the reef-building coral Acropora hyacinthus along a latitudinal cline that includes a well-established range expansion front in Japan using 2b-RAD sequencing. A total of 184 coral samples were collected across seven sites spanning from ~24°N to near its northern range front at ~33°N. We uncover the presence of three cryptic lineages of A. hyacinthus, which occupy discrete reefs within this region. Only one lineage is present along the expansion front and we find evidence for its historical occupation of marginal habitats. Within this lineage we also find evidence of bottleneck pressures associated with expansion events including higher clonality, increased linkage disequilibrium, and lower genetic diversity in range edge populations compared to core populations. Asymmetric migration between populations was also detected with lower migration from edge sites. Lastly, we describe genomic signatures of local adaptation potentially attributed to lower winter temperatures experienced at the more recently expanded northern populations. Together these data illuminate the genomic consequences of range expansion in a coral and highlight how adaptation to discrete environments along expansion fronts may facilitate further range expansion in this temperate coral lineage.

Phenotypic plasticity is aligned with phenological adaptation on both micro- and macroevolutionary timescales

In seasonally variable environments, phenotypic plasticity in phenology may be critical for adaptation to fluctuating environmental conditions. Using an 18-generation longitudinal dataset from natural damselfly populations, we show that phenology has strongly advanced. Individual fitness data suggest this is likely an adaptive response towards a temperature-dependent optimum. A laboratory experiment revealed that developmental plasticity qualitatively matches the temperature dependence of selection, partially explaining observed advance in phenology. Expanding our analysis to the macroevolutionary level, we use a database of over 1-million occurrence records and spatiotemporally matched temperature data from 49 Swedish Odonate species to infer macroevolutionary dynamics of phenology. Phenological plasticity was more closely aligned with adaptation for species that have recently colonised northern latitudes, but with higher phenological mismatch at lower latitudes. Our results show that phenological plasticity plays a key role in microevolutionary dynamics within a single species, and such plasticity may have facilitated post-Pleistocene range expansion in this insect clade.

Population Genomics Analysis with RAD, Reprised: Stacks 2

Restriction enzymes have been one of the primary tools in the population genetics toolkit for 50 years, being coupled with each new generation of technology to provide a more detailed view into the genetics of natural populations. Restriction site-Associated DNA protocols, which joined enzymes with short-read sequencing technology, have democratized the field of population genomics, providing a means to assay the underlying alleles in scores of populations. More than 10 years on, the technique has been widely applied across the tree of life and served as the basis for many different analysis techniques. Here, we provide a detailed protocol to conduct a RAD analysis from experimental design to de novo analysis-including parameter optimization-as well as reference-based analysis, all in Stacks version 2, which is designed to work with paired-end reads to assemble RAD loci up to 1000 nucleotides in length. The protocol focuses on major points of friction in the molecular approaches and downstream analysis, with special attention given to validating experimental analyses. Finally, the protocol provides several points of departure for further analysis.

Large-scale genome-wide study reveals climate adaptive variability in a cosmopolitan pest

Understanding the genetic basis of climatic adaptation is essential for predicting species' responses to climate change. However, intraspecific variation of these responses arising from local adaptation remains ambiguous for most species. Here, we analyze genomic data from diamondback moth (Plutella xylostella) collected from 75 sites spanning six continents to reveal that climate-associated adaptive variation exhibits a roughly latitudinal pattern. By developing an eco-genetic index that combines genetic variation and physiological responses, we predict that most P. xylostella populations have high tolerance to projected future climates. Using genome editing, a key gene, PxCad, emerged from our analysis as functionally temperature responsive. Our results demonstrate that P. xylostella is largely capable of tolerating future climates in most of the world and will remain a global pest beyond 2050. This work improves our understanding of adaptive variation along environmental gradients, and advances pest forecasting by highlighting the genetic basis for local climate adaptation.

Wolbachia-driven selective sweep in a range expanding insect species

Background

Evolutionary processes can cause strong spatial genetic signatures, such as local loss of genetic diversity, or conflicting histories from mitochondrial versus nuclear markers. Investigating these genetic patterns is important, as they may reveal obscured processes and players. The maternally inherited bacterium Wolbachia is among the most widespread symbionts in insects. Wolbachia typically spreads within host species by conferring direct fitness benefits, and/or by manipulating its host reproduction to favour infected over uninfected females. Under sufficient selective advantage, the mitochondrial haplotype associated with the favoured maternally-inherited symbiotic strains will spread (i.e. hitchhike), resulting in low mitochondrial genetic variation across the host species range.

Method

The common bluetail damselfly (Ischnura elegans: van der Linden, 1820) has recently emerged as a model organism for genetics and genomic signatures of range expansion during climate change. Although there is accumulating data on the consequences of such expansion on the genetics of I. elegans, no study has screened for Wolbachia in the damselfly genus Ischnura. Here, we present the biogeographic variation in Wolbachia prevalence and penetrance across Europe and Japan (including samples from 17 populations), and from close relatives in the Mediterranean area (i.e. I. genei: Rambur, 1842; and I. saharensis: Aguesse, 1958).

Results

Our data reveal (a) multiple Wolbachia-strains, (b) potential transfer of the symbiont through hybridization, (c) higher infection rates at higher latitudes, and (d) reduced mitochondrial diversity in the north-west populations, indicative of hitchhiking associated with the selective sweep of the most common strain. We found low mitochondrial haplotype diversity in the Wolbachia-infected north-western European populations (Sweden, Scotland, the Netherlands, Belgium, France and Italy) of I. elegans, and, conversely, higher mitochondrial diversity in populations with low penetrance of Wolbachia (Ukraine, Greece, Montenegro and Cyprus). The timing of the selective sweep associated with infected lineages was estimated between 20,000 and 44,000 years before present, which is consistent with the end of the last glacial period about 20,000 years.

Conclusions

Our findings provide an example of how endosymbiont infections can shape spatial variation in their host evolutionary genetics during postglacial expansion. These results also challenge population genetic studies that do not consider the prevalence of symbionts in many insects, which we show can impact geographic patterns of mitochondrial genetic diversity.

Environmental association modelling with loci under divergent selection predicts the distribution range of a lizard

During the historical building of a species range, individual colonizers have to confront different ecological challenges, and the capacity of the species to broaden its range may depend on the total amount of adaptive genetic variation supplied by evolution. We set out to increase our understanding of what defines a range and the role of underlying genetics by trying to predict an entire species' range from the geographical distribution of its genetic diversity under selection. We sampled five populations of the western Mediterranean lizard Psammodromus algirus that inhabit a noticeable environmental gradient of temperature and precipitation. We correlated the genotypes of 95 individuals (18-20 individuals per population) for 21 SNPs putatively under selection with environmental scores on a bioclimatic gradient, using 1 × 1 km²  grid cells as sampling units. By extrapolating the resulting model to all possible combinations of alleles, we inferred all the geographic cells that were theoretically suitable for a given amount of genetic variance under selection. The inferred distribution range overlapped to a large extent with the realized range of the species (77.46% of overlap), including an accurate prediction of internal gaps and range borders. Our results suggest an adaptability threshold determined by the amount of genetic variation available that would be required to warrant adaptation beyond a certain limit of environmental variation. These results support the idea that the expansion of a species' range can be ultimately linked to the arising of new variants under selection (either newly selected variants from standing genetic variation or innovative mutations under selection).

The genomic signal of local environmental adaptation in Aedes aegypti mosquitoes

Local adaptation is important when predicting arthropod-borne disease risk because of its impacts on vector population fitness and persistence. However, the extent that vector populations are adapted to the environment generally remains unknown. Despite low population structure and high gene flow in Aedes aegypti mosquitoes across Panama, excepting the province of Bocas del Toro, we identified 128 candidate SNPs, clustered within 17 genes, which show a strong genomic signal of local environmental adaptation. This putatively adaptive variation occurred across fine geographical scales with the composition and frequency of candidate adaptive loci differing between populations in wet tropical environments along the Caribbean coast and dry tropical conditions typical of the Pacific coast. Temperature and vegetation were important predictors of adaptive genomic variation in Ae. aegypti with several potential areas of local adaptation identified. Our study lays the foundations of future work to understand whether environmental adaptation in Ae. aegypti impacts the arboviral disease landscape and whether this could either aid or hinder efforts of population control.

Genome assembly, sex-biased gene expression and dosage compensation in the damselfly Ischnura elegans

The evolution of sex chromosomes, and patterns of sex-biased gene expression and dosage compensation, are poorly known among early winged insects such as odonates. We assembled and annotated the genome of Ischnura elegans (blue-tailed damselfly), which, like other odonates, has a male-hemigametic sex-determining system (X0 males, XX females). By identifying X-linked genes in I. elegans and their orthologs in other insect genomes, we found homologies between the X chromosome in odonates and chromosomes of other orders, including the X chromosome in Coleoptera. Next, we showed balanced expression of X-linked genes between sexes in adult I. elegans, i.e. evidence of dosage compensation. Finally, among the genes in the sex-determining pathway only fruitless was found to be X-linked, while only doublesex showed sex-biased expression. This study reveals partly conserved sex chromosome synteny and independent evolution of dosage compensation among insect orders separated by several hundred million years of evolutionary history.

Latitudinal clines in sexual selection, sexual size dimorphism and sex-specific genetic dispersal during a poleward range expansion

Range expansions can be shaped by sex differences in behaviours and other phenotypic traits affecting dispersal and reproduction. Here, we investigate sex differences in morphology, behaviour and genomic population differentiation along a climate-mediated range expansion in the common bluetail damselfly (Ischnura elegans) in northern Europe. We sampled 65 sites along a 583-km gradient spanning the I. elegans range in Sweden and quantified latitudinal gradients in site relative abundance, sex ratio and sex-specific shifts in body size and mating status (a measure of sexual selection). Using single nucleotide polymorphism (SNP) data for 426 individuals from 25 sites, we further investigated sex-specific landscape and climatic effects on neutral genetic connectivity and migration patterns. We found evidence for sex differences associated with the I. elegans range expansion, namely (a) increased male body size with latitude, but no latitudinal effect on female body size, resulting in reduced sexual dimorphism towards the range limit, (b) a steeper decline in male genetic similarity with increasing geographic distance than in females, (c) male-biased genetic migration propensity and (d) a latitudinal cline in migration distance (increasing migratory distances towards the range margin), which was stronger in males. Cooler mean annual temperatures towards the range limit were associated with increased resistance to gene flow in both sexes. Sex ratios became increasingly male biased towards the range limit, and there was evidence for a changed sexual selection regime shifting from favouring larger males in the south to favouring smaller males in the north. Our findings suggest sex-specific spatial phenotype sorting at the range limit, where larger males disperse more under higher landscape resistance associated with cooler climates. The combination of latitudinal gradients in sex-biased dispersal, increasing male body size and (reduced) sexual size dimorphism should have emergent consequences for sexual selection dynamics and the mating system at the expanding range front. Our study illustrates the importance of considering sex differences in the study of range expansions driven by ongoing climate change.

Evolution of cold tolerance and thermal plasticity in life history, behaviour and physiology during a poleward range expansion

Many species that are moving polewards encounter novel thermal regimes to which they have to adapt. Therefore, rapid evolution of thermal tolerance and of thermal plasticity in fitness-related traits in edge populations can be crucial for the success and speed of range expansions. We tested for adaptation in cold tolerance and in life history, behavioural and physiological traits and their thermal plasticity during a poleward range expansion. We reconstructed the thermal performance curves of life history (survival, growth and development rates), behaviour (food intake) and cold tolerance (chill coma recovery time) in the aquatic larval stage of the damselfly Ischnura elegans that is currently showing a poleward range expansion in northern Europe. We studied larvae from three edge and three core populations using a common-garden experiment. Consistent with the colder annual temperatures, larvae at the expansion front evolved an improved cold tolerance. The edge populations showed no overall (across temperatures) evolution of a faster life history that would improve their range-shifting ability. Moreover, consistent with damselfly edge populations from colder latitudes, edge populations evolved at the highest rearing temperature (28°C) a faster development rate, likely to better exploit the rare periods with higher temperatures. This was associated with a higher food intake and a lower metabolic rate. In conclusion, our results suggest that the edge populations rapidly evolved adaptive changes in trait means and thermal plasticity to the novel thermal conditions at the edge front. Our results highlight the importance of considering besides trait plasticity and the evolution of trait means, also the evolution of trait plasticity to improve forecasts of responses to climate change.

Bioclimatic Modelling Identifies Suitable Habitat for the Establishment of the Invasive European Paper Wasp (Hymenoptera: Vespidae) across the Southern Hemisphere

Species distribution models (SDMs) are tools used by ecologists to help predict the spread of invasive species. Information provided by these models can help direct conservation and biosecurity efforts by highlighting areas likely to contain species of interest. In this study, two models were created to investigate the potential range expansion of Polistes dominula Christ (Hymenoptera: Vespidae) in the southern hemisphere. This palearctic species has spread to invade North and South America, South Africa, Australia, and more recently New Zealand. Using the BIOCLIM and MAXENT modelling methods, regions that were suitable for P. dominula were identified based on climate data across four regions in the southern hemisphere. In South America areas of central Chile, eastern Argentina, parts of Uruguay, and southern Brazil were identified as climatically suitable for the establishment of P. dominula. Similarly, southern parts of South Africa and Australia were identified by the model to be suitable as well as much of the North Island and east of the South Island of New Zealand. Based on outputs from both models, significant range expansion by P. dominula is possible across its more southern invaded ranges.

Signatures of selection in a recent invasion reveal adaptive divergence in a highly vagile invasive species

A detailed understanding of population genetics in invasive populations helps us to identify drivers of successful alien introductions. Here, we investigate putative signals of selection in Australian populations of invasive common starlings, Sturnus vulgaris, and seek to understand how these have been influenced by introduction history. We used reduced representation sequencing to determine population structure, and identify Single Nucleotide Polymorphisms (SNPs) that are putatively under selection. We found that since their introduction into Australia, starling populations have become genetically differentiated despite the potential for high levels of dispersal, and that starlings have responded to selective pressures imposed by a wide range of environmental conditions across their geographic range. Isolation by distance appears to have played a strong role in determining genetic substructure across the starling's Australian range. Analyses of candidate SNPs that are putatively under selection indicated that aridity, precipitation and temperature may be important factors driving adaptive variation across the starling's invasive range in Australia. However, we also noted that the historic introduction regime may leave footprints on sites flagged as being under adaptive selection, and encourage critical interpretation of selection analyses in non-native populations.

Association of genetic and climatic variability in giant sequoia, Sequoiadendron giganteum, reveals signatures of local adaptation along moisture-related gradients

Uncovering the genetic basis of local adaptation is a major goal of evolutionary biology and conservation science alike. In an era of climate change, an understanding of how environmental factors shape adaptive diversity is crucial to predicting species response and directing management. Here, we investigate patterns of genomic variation in giant sequoia, an iconic and ecologically important tree species, using 1,364 bi-allelic single nucleotide polymorphisms (SNPs). We use an F _ST outlier test and two genotype-environment association methods, latent factor mixed models (LFMMs) and redundancy analysis (RDA), to detect complex signatures of local adaptation. Results indicate 79 genomic regions of potential adaptive importance, with limited overlap between the detection methods. Of the 58 loci detected by LFMM, 51 showed strong correlations to a precipitation-driven composite variable and seven to a temperature-related variable. RDA revealed 24 outlier loci with association to climate variables, all of which showed strongest relationship to summer precipitation. Nine candidate loci were indicated by two methods. After correcting for geographic distance, RDA models using climate predictors accounted for 49% of the explained variance and showed significant correlations between SNPs and climatic factors. Here, we present evidence of local adaptation in giant sequoia along gradients of precipitation and provide a first step toward identifying genomic regions of adaptive significance. The results of this study will provide information to guide management strategies that seek to maximize adaptive potential in the face of climate change.

Contrasting patterns of local adaptation along climatic gradients between a sympatric parasitic and autotrophic tree species

Sympatric tree species are subject to similar climatic drivers, posing a question as to whether they display comparable adaptive responses. However, no study has explicitly examined local adaptation of co-occurring parasitic and autotrophic plant species to the abiotic environment. Here we test the hypotheses that a generalist parasitic tree would display a weaker signal of selection and that genomic variation would associate with fewer climatic variables (particularly precipitation) but have similar spatial patterns to a sympatric autotrophic tree species. To test these hypotheses, we collected samples from 17 sites across the range of two tree species, the hemiparasite Nuytsia floribunda (n = 264) and sympatric autotroph Melaleuca rhaphiophylla (n = 272). We obtained 5,531 high-quality genome-wide single nucleotide polymorphisms (SNPs) for M. rhaphiophylla and 6,727 SNPs for N. floribunda using DArTseq genome scan technology. Population differentiation and environmental association approaches were used to identify signals of selection. Generalized dissimilarly modelling was used to detect climatic and spatial patterns of local adaptation across climatic gradients. Overall, 322 SNPs were identified as putatively adaptive for the autotroph, while only 57 SNPs were identified for the parasitic species. We found genomic variation to associate with different sets of bioclimatic variables for each species, with precipitation relatively less important for the parasite. Spatial patterns of predicted adaptive variability were different and indicate that co-occurring species with disparate life history traits may not respond equally to selective pressures (i.e., temperature and precipitation). Together, these findings provide insight into local adaptation of sympatric parasitic and autotrophic tree species to abiotic environments.

Thermal evolution of life history and heat tolerance during range expansions toward warmer and cooler regions

Species' range edges are expanding to both warmer and cooler regions. Yet, no studies directly compared the changes in range-limiting traits within the same species during both types of range expansions. To increase our mechanistic understanding of range expansions, it is crucial to disentangle the contributions of plastic and genetic changes in these traits. The aim of this study was to test for plastic and evolutionary changes in heat tolerance, life history, and behavior, and compare these during range expansions toward warmer and cooler regions. Using laboratory experiments we reconstructed the thermal performance curves (TPCurves) of larval life history (survival, growth, and development rates) and larval heat tolerance (CTmax) across two recent range expansions from the core populations in southern France toward a warmer (southeastern Spain) and a cooler (northwestern Spain) region in Europe by the damselfly Ischnura elegans. First-generation larvae from field-collected mothers were reared across a range of temperatures (16°-28°C) in incubators. The range expansion to the warmer region was associated with the evolution of a greater ability to cope with high temperatures (increased mean and thermal plasticity of CTmax), faster development, and, in part, a faster growth, indicating a higher time constraints caused by a shorter time frame available for larval development associated with a transition to a greater voltinism. Our results thereby support the emerging pattern that plasticity in heat tolerance alone is inadequate to adapt to new thermal regimes. The range expansion to the cooler region was associated with faster growth indicating countergradient variation without a change in CTmax. The evolution of a faster growth rate during both range expansions could be explained by a greater digestive efficiency rather than an increased food intake. Our results highlight that range expansions to warmer and cooler regions can result in similar evolutionary changes in the TPCurves for life history, and no opposite changes in heat tolerance.

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.

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus)

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site-associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.

The odonate phenotypic database, a new open data resource for comparative studies of an old insect order

We present The Odonate Phenotypic Database (OPD): an online data resource of dragonfly and damselfly phenotypes (Insecta: Odonata). Odonata is a relatively small insect order that currently consists of about 6400 species belonging to 32 families. The database consists of multiple morphological, life-history and behavioral traits, and biogeographical information collected from literature sources. We see taxon-specific phenotypic databases from Odonata and other organismal groups as becoming an increasing valuable resource in comparative studies. Our database has phenotypic records for 1011 of all 6400 known odonate species. The database is accessible at http://www.odonatephenotypicdatabase.org/, and a static version with an information file about the variables in the database is archived at Dryad.

Characterization of Dendrolimus houi Lajonquiere (Lepidoptera: Lasiocampidae) Transcriptome across All Life Stages

Dendrolimus houi Lajonquiere is a phytophagous caterpillar infesting many economically important coniferous tree species in China, causing serious economic and ecological environment losses. Based on previous research, it has one generation per year in South China and East China in contrast to two generations per year in Yunnan province in southwestern China. The species is potentially resilient to climatic extremes in these regions with the eggs and 1st instar larvae surviving in the winter (5 °C), older instar larvae and pupae surviving high temperatures in the summer (35 °C), suggesting some temperature stress tolerance during different developmental stages. However, little is known in this species at the genetic and genomic level. In this study, we used high throughput sequencing to obtain transcriptome data from different developmental stages (eggs, 1st-3rd instar larvae, 4th-5th instar larvae, 6th-7th instar larvae, pupae, male and female adults), which were collected from Fujian province. In total, we obtained approximately 90 Gb of data, from which 33,720 unigenes were assembled and 17,797 unigenes were annotated. We furtherly analyzed the differentially expressed genes (DGEs) across all stages, the largest number between the eggs and 1st instar larvae stage and gene expression varied significantly in different developmental stages. Furthermore, 4138 SSR genes and 114,977 SNP loci were screened from transcriptome data. This paper will be a foundation for further study towards improved integrated pest management strategies for this species.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

Predicting the success of an invader: Niche shift versus niche conservatism

Invasive species can encounter environments different from their source populations, which may trigger rapid adaptive changes after introduction (niche shift hypothesis). To test this hypothesis, we investigated whether postintroduction evolution is correlated with contrasting environmental conditions between the European invasive and source ranges in the Asian tiger mosquito Aedes albopictus. The comparison of environmental niches occupied in European and source population ranges revealed more than 96% overlap between invasive and source niches, supporting niche conservatism. However, we found evidence for postintroduction genetic evolution by reanalyzing a published ddRADseq genomic dataset from 90 European invasive populations using genotype-environment association (GEA) methods and generalized dissimilarity modeling (GDM). Three loci, among which a putative heat-shock protein, exhibited significant allelic turnover along the gradient of winter precipitation that could be associated with ongoing range expansion. Wing morphometric traits weakly correlated with environmental gradients within Europe, but wing size differed between invasive and source populations located in different climatic areas. Niche similarities between source and invasive ranges might have facilitated the establishment of populations. Nonetheless, we found evidence for environmental-induced adaptive changes after introduction. The ability to rapidly evolve observed in invasive populations (genetic shift) together with a large proportion of unfilled potential suitable areas (80%) pave the way to further spread of Ae. albopictus in Europe.

Sampling strategy optimization to increase statistical power in landscape genomics: A simulation-based approach

An increasing number of studies are using landscape genomics to investigate local adaptation in wild and domestic populations. Implementation of this approach requires the sampling phase to consider the complexity of environmental settings and the burden of logistical constraints. These important aspects are often underestimated in the literature dedicated to sampling strategies. In this study, we computed simulated genomic data sets to run against actual environmental data in order to trial landscape genomics experiments under distinct sampling strategies. These strategies differed by design approach (to enhance environmental and/or geographical representativeness at study sites), number of sampling locations and sample sizes. We then evaluated how these elements affected statistical performances (power and false discoveries) under two antithetical demographic scenarios. Our results highlight the importance of selecting an appropriate sample size, which should be modified based on the demographic characteristics of the studied population. For species with limited dispersal, sample sizes above 200 units are generally sufficient to detect most adaptive signals, while in random mating populations this threshold should be increased to 400 units. Furthermore, we describe a design approach that maximizes both environmental and geographical representativeness of sampling sites and show how it systematically outperforms random or regular sampling schemes. Finally, we show that although having more sampling locations (between 40 and 50 sites) increase statistical power and reduce false discovery rate, similar results can be achieved with a moderate number of sites (20 sites). Overall, this study provides valuable guidelines for optimizing sampling strategies for landscape genomics experiments.