Genomic Shifts, Phenotypic Clines, and Fitness Costs Associated With Cold Tolerance in the Asian Tiger Mosquito.

Climatic variation is a key driver of genetic differentiation and phenotypic traits evolution, and local adaptation to temperature is expected in widespread species. We investigated phenotypic and genomic changes in the native range of the Asian tiger mosquito, Aedes albopictus. We first refine the phylogeographic structure based on genome-wide regions (1,901 double-digest restriction-site associated DNA single nucleotide polymophisms [ddRAD SNPs]) from 41 populations. We then explore the patterns of cold adaptation using phenotypic traits measured in common garden (wing size and cold tolerance) and genotype-temperature associations at targeted candidate regions (51,706 exon-capture SNPs) from nine populations. We confirm the existence of three evolutionary lineages including clades A (Malaysia, Thailand, Cambodia, and Laos), B (China and Okinawa), and C (South Korea and Japan). We identified temperature-associated differentiation in 15 out of 221 candidate regions but none in ddRAD regions, supporting the role of directional selection in detected genes. These include genes involved in lipid metabolism and a circadian clock gene. Most outlier SNPs are differently fixed between clades A and C, whereas clade B has an intermediate pattern. Females are larger at higher latitude yet produce no more eggs, which might favor the storage of energetic reserves in colder climate. Nondiapausing eggs from temperate populations survive better to cold exposure than those from tropical populations, suggesting they are protected from freezing damages but this cold tolerance has a fitness cost in terms of egg viability. Altogether, our results provide strong evidence for the thermal adaptation of A. albopictus across its wide temperature range.

The molecular signatures of compatible and incompatible pollination in Arabidopsis.

BACKGROUND: Fertilization in flowering plants depends on the early contact and acceptance of pollen grains by the receptive papilla cells of the stigma. Deciphering the specific transcriptomic response of both pollen and stigmatic cells during their interaction constitutes an important challenge to better our understanding of this cell recognition event. RESULTS: Here we describe a transcriptomic analysis based on single nucleotide polymorphisms (SNPs) present in two Arabidopsis thaliana accessions, one used as female and the other as male. This strategy allowed us to distinguish 80% of transcripts according to their parental origins. We also developed a tool which predicts male/female specific expression for genes without SNP. We report an unanticipated transcriptional activity triggered in stigma upon incompatible pollination and show that following compatible interaction, components of the pattern-triggered immunity (PTI) pathway are induced on the female side. CONCLUSIONS: Our work unveils the molecular signatures of compatible and incompatible pollinations both at the male and female side. We provide invaluable resource and tools to identify potential new molecular players involved in pollen-stigma interaction.

Live-cell imaging of early events following pollen perception in self-incompatible Arabidopsis thaliana.

Early events occurring at the surface of the female organ are critical for plant reproduction, especially in species with a dry stigma. After landing on the stigmatic papilla cells, the pollen hydrates and germinates a tube, which penetrates the cell wall and grows towards the ovules to convey the male gametes to the embryo sac. In self-incompatible species within the Brassicaceae, these processes are blocked when the stigma encounters an incompatible pollen. Based on the generation of self-incompatible Arabidopsis lines and by setting up a live imaging system, we showed that control of pollen hydration has a central role in pollen selectivity. The faster the pollen pumps water from the papilla during an initial period of 10 min, the faster it germinates. Furthermore, we found that the self-incompatibility barriers act to block the proper hydration of incompatible pollen and, when hydration is promoted by high humidity, an additional control prevents pollen tube penetration into the stigmatic wall. In papilla cells, actin bundles focalize at the contact site with the compatible pollen but not with the incompatible pollen, raising the possibility that stigmatic cells react to the mechanical pressure applied by the invading growing tube.

Combined Proteomic and Metabolomic Profiling of the Arabidopsis thaliana vps29 Mutant Reveals Pleiotropic Functions of the Retromer in Seed Development.

The retromer is a multiprotein complex conserved from yeast to humans, which is involved in intracellular protein trafficking and protein recycling. Selection of cargo proteins transported by the retromer depends on the core retromer subunit composed of the three vacuolar protein sorting (VPS) proteins, namely VPS26, VPS29, and VPS35. To gain a better knowledge of the importance of the plant retromer in protein sorting, we carried out a comparative proteomic and metabolomic analysis of Arabidopsis thaliana seeds from the wild-type and the null-retromer mutant vps29. Here, we report that the retromer mutant displays major alterations in the maturation of seed storage proteins and synthesis of lipid reserves, which are accompanied by severely impaired seed vigor and longevity. We also show that the lack of retromer components is counterbalanced by an increase in proteins involved in intracellular trafficking, notably members of the Ras-related proteins in brain (RAB) family proteins. Our study suggests that loss of the retromer stimulates energy metabolism, affects many metabolic pathways, including that of cell wall biogenesis, and triggers an osmotic stress response, underlining the importance of retromer function in seed biology.

Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing.

The capacity of mosquitoes to resist insecticides threatens the control of diseases such as dengue and malaria. Until alternative control tools are implemented, characterizing resistance mechanisms is crucial for managing resistance in natural populations. Insecticide biodegradation by detoxification enzymes is a common resistance mechanism; however, the genomic changes underlying this mechanism have rarely been identified, precluding individual resistance genotyping. In particular, the role of copy number variations (CNVs) and polymorphisms of detoxification enzymes have never been investigated at the genome level, although they can represent robust markers of metabolic resistance. In this context, we combined target enrichment with high-throughput sequencing for conducting the first comprehensive screening of gene amplifications and polymorphisms associated with insecticide resistance in mosquitoes. More than 760 candidate genes were captured and deep sequenced in several populations of the dengue mosquito Ae. aegypti displaying distinct genetic backgrounds and contrasted resistance levels to the insecticide deltamethrin. CNV analysis identified 41 gene amplifications associated with resistance, most affecting cytochrome P450s overtranscribed in resistant populations. Polymorphism analysis detected more than 30,000 variants and strong selection footprints in specific genomic regions. Combining Bayesian and allele frequency filtering approaches identified 55 nonsynonymous variants strongly associated with resistance. Both CNVs and polymorphisms were conserved within regions but differed across continents, confirming that genomic changes underlying metabolic resistance to insecticides are not universal. By identifying novel DNA markers of insecticide resistance, this study opens the way for tracking down metabolic changes developed by mosquitoes to resist insecticides within and among populations.

Peroxisome extensions deliver the Arabidopsis SDP1 lipase to oil bodies.

Lipid droplets/oil bodies (OBs) are lipid-storage organelles that play a crucial role as an energy resource in a variety of eukaryotic cells. Lipid stores are mobilized in the case of food deprivation or high energy demands--for example, during certain developmental processes in animals and plants. OB degradation is achieved by lipases that hydrolyze triacylglycerols (TAGs) into free fatty acids and glycerol. In the model plant Arabidopsis thaliana, Sugar-Dependent 1 (SDP1) was identified as the major TAG lipase involved in lipid reserve mobilization during seedling establishment. Although the enzymatic activity of SDP1 is associated with the membrane of OBs, its targeting to the OB surface remains uncharacterized. Here we demonstrate that the core retromer, a complex involved in protein trafficking, participates in OB biogenesis, lipid store degradation, and SDP1 localization to OBs. We also report an as-yet-undescribed mechanism for lipase transport in eukaryotic cells, with SDP1 being first localized to the peroxisome membrane at early stages of seedling growth and then possibly moving to the OB surface through peroxisome tubulations. Finally, we show that the timely transfer of SDP1 to the OB membrane requires a functional core retromer. In addition to revealing previously unidentified functions of the retromer complex in plant cells, our work provides unanticipated evidence for the role of peroxisome dynamics in interorganelle communication and protein transport.

Cellular auxin homeostasis under high temperature is regulated through a sorting NEXIN1-dependent endosomal trafficking pathway.

High-temperature-mediated adaptation in plant architecture is linked to the increased synthesis of the phytohormone auxin, which alters cellular auxin homeostasis. The auxin gradient, modulated by cellular auxin homeostasis, plays an important role in regulating the developmental fate of plant organs. Although the signaling mechanism that integrates auxin and high temperature is relatively well understood, the cellular auxin homeostasis mechanism under high temperature is largely unknown. Using the Arabidopsis thaliana root as a model, we demonstrate that under high temperature, roots counterbalance the elevated level of intracellular auxin by promoting shootward auxin efflux in a PIN-FORMED2 (PIN2)-dependent manner. Further analyses revealed that high temperature selectively promotes the retrieval of PIN2 from late endosomes and sorts them to the plasma membrane through an endosomal trafficking pathway dependent on SORTING NEXIN1. Our results demonstrate that recycling endosomal pathway plays an important role in facilitating plants adaptation to increased temperature.

Mechanisms governing the endosomal membrane recruitment of the core retromer in Arabidopsis.

The retromer complex localizes to endosomal membranes and is involved in protein trafficking. In mammals, it is composed of a dimer of sorting nexins and of the core retromer consisting of vacuolar protein sorting (VPS)26, VPS29, and VPS35. Although homologs of these proteins have been identified in plants, how the plant retromer functions remains elusive. To better understand the role of VPS components in the assembly and function of the core retromer, we characterize here Arabidopsis vps26-null mutants. We show that impaired VPS26 function has a dramatic effect on VPS35 levels and causes severe phenotypic defects similar to those observed in vps29-null mutants. This implies that functions of plant VPS26, VPS29, and VPS35 are tightly linked. Then, by combining live-cell imaging with immunochemical and genetic approaches, we report that VPS35 alone is able to bind to endosomal membranes and plays an essential role in VPS26 and VPS29 membrane recruitment. We also show that the Arabidopsis Rab7 homolog RABG3f participates in the recruitment of the core retromer to the endosomal membrane by interacting with VPS35. Altogether our data provide original information on the molecular interactions that mediate assembly of the core retromer in plants.

Pre-selecting resistance against individual Bti Cry toxins facilitates the development of resistance to the Bti toxins cocktail.

The bioinsecticide Bacillus thuringiensis subsp. israelensis is a larvicide used worldwide for mosquito control, which contains three Cry toxins and one Cyt toxin. We investigated for the first time in Aedes aegypti (1) the evolution of resistance and cross-resistance of strains selected with each Cry toxin, and (2) the effect of pre-selection with Cry toxin on the evolution of resistance to a mix of Bti toxins. Cross resistance was higher between Cry4Ba and Cry11Aa than between Cry4Aa and either Cry4Ba or Cry11Aa, suggesting both common and specific mechanisms of resistance. Pre-selecting resistance to each Cry toxins facilitated the development of resistance to the full Bti toxins cocktail.

Analyses of sorting nexins reveal distinct retromer-subcomplex functions in development and protein sorting in Arabidopsis thaliana.

Sorting nexins (SNXs) are conserved eukaryotic proteins that associate with three types of vacuolar protein sorting (VPS) proteins to form the retromer complex. How SNXs act in this complex and whether they might work independently of the retromer remains elusive. Here, we show by genetic and cell imaging approaches that the Arabidopsis thaliana SNX1 protein recruits SNX2 at the endosomal membrane, a process required for SNX1-SNX2 dimer activity. We report that, in contrast with the mammalian retromer, SNXs are dispensable for membrane binding and function of the retromer complex. We also show that VPS retromer components can work with or independently of SNXs in the trafficking of seed storage proteins, which reveals distinct functions for subcomplexes of the plant retromer. Finally, we provide compelling evidence that the combined loss of function of SNXs and VPS29 leads to embryo or seedling lethality, underlining the essential role of these proteins in development.

AtSNX1 defines an endosome for auxin-carrier trafficking in Arabidopsis.

Polarized cellular distribution of the phytohormone auxin and its carriers is essential for normal plant growth and development. Polar auxin transport is maintained by a network of auxin influx (AUX) and efflux (PIN) carriers. Both auxin transport and PIN protein cycling between the plasma membrane and endosomes require the activity of the endosomal GNOM; however, intracellular routes taken by these carriers remain largely unknown. Here we show that Arabidopsis thaliana SORTING NEXIN 1 (AtSNX1) is involved in the auxin pathway and that PIN2, but not PIN1 or AUX1, is transported through AtSNX1-containing endosomes. We demonstrate that the snx1-null mutant exhibits multiple auxin-related defects and that loss of function of AtSNX1 severely enhances the phenotype of a weak gnom mutant. In root cells, we further show that AtSNX1 localizes to an endosomal compartment distinct from GNOM-containing endosomes, and that PIN2 accumulates in this compartment after treatment with the phosphatidylinositol-3-OH kinase inhibitor wortmannin or after a gravity stimulus. Our data reveal the existence of a novel endosomal compartment involved in PIN2 endocytic sorting and plant development.

Receptor kinase signalling in plants and animals: distinct molecular systems with mechanistic similarities.

Plant genomes encode large numbers of receptor kinases that are structurally related to the tyrosine and serine/threonine families of receptor kinase found in animals. Here, we describe recent advances in the characterisation of several of these plant receptor kinases at the molecular level, including the identification of receptor complexes, small polypeptide ligands and cytosolic proteins involved in signal transduction and receptor downregulation. Phylogenetic analysis indicates that plant receptor kinases have evolved independently of the receptor kinase families found in animals. This hypothesis is supported by functional studies that have revealed differences between receptor kinase signalling in plants and animals, particularly concerning their interactions with cytosolic proteins. Despite these dissimilarities, however, plant and animal receptor kinases share many common features, such as their single membrane-pass structure, their inclusion in membrane-associated complexes, the involvement of dimerisation and trans autophosphorylation in receptor activation, and the existence of inhibitors and phosphatases that downregulate receptor activity. These points of convergence may represent features that are essential for a functional receptor-kinase signalling system.

KATANIN and cortical microtubule organization have a pivotal role in early pollen tube guidance.

Following pollen deposition on the receptive surface of the stigma, pollen germinates a tube that carries male gametes toward the ovule where fertilization occurs. As soon as it emerges from the pollen grain, the pollen tube has to be properly guided through the pistil tissues so as to reach the ovule and ensure double fertilization. Chemical attractants, nutrients as well as receptor kinase-dependent signaling pathways have been implicated in this guidance. Recently, we showed in Arabidopsis that the microtubule severing enzyme KATANIN, by acting both on cortical microtubule (CMT) dynamics and cellulose microfibril (CMF) deposition, conferred particular mechanical properties to the papilla cell wall that act as active guidance factors. Here we confirm the importance of KATANIN and CMT orientation in pollen tube directionality by examining another katanin mutant.

Molecular bases of P450-mediated resistance to the neonicotinoid insecticide imidacloprid in the mosquito Ae. aegypti.

Resistance to chemical insecticides including pyrethroids, the main insecticide class used against mosquitoes, has re-kindled interest in the use of neonicotinoids. In this context, the present study aimed to characterize the molecular basis of neonicotinoid resistance in the mosquito Aedes aegypti. Resistance mechanisms were studied by combining transcriptomic and genomic data obtained from a laboratory strain selected at the larval stage after 30 generations of exposure to imidacloprid (Imida-R line). After thirty generations of selection, larvae of the Imida-R line showed an 8-fold increased resistance to imidacloprid and a significant cross-tolerance to the pyrethroids permethrin and deltamethrin. Cross-resistance to pyrethroids was only observed in adults when larvae were previously exposed to imidacloprid suggesting a low but inducible expression of resistance alleles at the adult stage. Resistance of the Imida-R line was associated with a slower larval development time in females. Multiple detoxification enzymes were over-transcribed in larvae in association with resistance including the P450s CYP6BB2, CYP9M9 and CYP6M11 previously associated with pyrethroid resistance. Some of them together with their redox partner NADPH P450 reductase were also affected by non-synonymous mutations associated with resistance. Combining genomic and transcriptomic data allowed identifying promoter variations associated with the up-regulation of CYP6BB2 in the resistant line. Overall, these data confirm the key role of P450s in neonicotinoid resistance in Ae. aegypti and their potential to confer cross-resistance to pyrethroids, raising concerns about the use of neonicotinoids for resistance management in this mosquito species.

Impact of selection regime and introgression on deltamethrin resistance in the arbovirus vector Aedes aegypti - a comparative study between contrasted situations in New Caledonia and French Guiana.

BACKGROUND: Pyrethroid insecticides such as deltamethrin have been massively used against Aedes aegypti leading to the spread of resistance alleles worldwide. In an insecticide resistance management context, we evaluated the temporal dynamics of deltamethrin resistance using two distinct populations carrying resistant alleles at different frequencies. Three different scenarios were followed: a continuous selection, a full release of selection, or a repeated introgression with susceptible individuals. The responses of each population to these selection regimes were measured across five generations by bioassays and by monitoring the frequency of knockdown resistance (kdr) mutations and the transcription levels and copy number variations of key detoxification enzymes. RESULTS: Knockdown resistance mutations, overexpression and copy number variations of detoxification enzymes as a mechanism of metabolic resistance to deltamethrin was found and maintained under selection across generations. On comparison, the release of insecticide pressure for five generations did not affect resistance levels and resistance marker frequencies. However, introgressing susceptible alleles drastically reduced deltamethrin resistance in only three generations. CONCLUSION: The present study confirmed that strategies consisting to stop deltamethrin spraying are likely to fail when the frequencies of resistant alleles are too high and the fitness cost associated to resistance is low. In dead-end situations like in French Guiana where alternative insecticides are not available, alternative control strategies may provide a high benefit for vector control, particularly if they favor the introgression of susceptible alleles in natural populations. © 2021 Society of Chemical Industry.

A genomic amplification affecting a carboxylesterase gene cluster confers organophosphate resistance in the mosquito Aedes aegypti: From genomic characterization to high-throughput field detection.

By altering gene expression and creating paralogs, genomic amplifications represent a key component of short-term adaptive processes. In insects, the use of insecticides can select gene amplifications causing an increased expression of detoxification enzymes, supporting the usefulness of these DNA markers for monitoring the dynamics of resistance alleles in the field. In this context, the present study aims to characterize a genomic amplification event associated with resistance to organophosphate insecticides in the mosquito Aedes aegypti and to develop a molecular assay to monitor the associated resistance alleles in the field. An experimental evolution experiment using a composite population from Laos supported the association between the over-transcription of multiple contiguous carboxylesterase genes on chromosome 2 and resistance to multiple organophosphate insecticides. Combining whole genome sequencing and qPCR on specific genes confirmed the presence of a ~100-Kb amplification spanning at least five carboxylesterase genes at this locus with the co-existence of multiple structural duplication haplotypes. Field data confirmed their circulation in South-East Asia and revealed high copy number polymorphism among and within populations suggesting a trade-off between this resistance mechanism and associated fitness costs. A dual-color multiplex TaqMan assay allowing the rapid detection and copy number quantification of this amplification event in Ae. aegypti was developed and validated on field populations. The routine use of this novel assay will improve the tracking of resistance alleles in this major arbovirus vector.

Variability of the self-incompatibility reaction in Brassica oleracea L. with S 15 haplotype.

Self-incompatibility (SI) is thought to have played a key role in the evolution of species as it promotes their outcrossing through the recognition and rejection of self-pollen grains. In most species, SI is under the control of a complex, multiallelic S-locus. The recognition system is associated with quantitative variations of the strength of the SI reaction; the origin of these variations is still not elucidated. To define the genetic regulations involved, we studied the variability of the SI response in homozygous S 15 S 15 plants in cauliflower. These plants were obtained from a self-progeny of a self-compatible (SC) plant heterozygous for S 15, which was generated after five selfing generations from one strongly self-incompatible initial plant. We found a continuous phenotypic variation for SI response in the offspring plants homozygous for the S 15 haplotype, from the strict SI reaction to self-compatibility, with a great proportion of the plants being partially self-compatible (PSC). Decrease in SI levels was also observed during the life of the flower. The number of pollen tubes passing through the stigma barrier was higher when counted 3 or 5 days after pollination than one day after pollination. Analysis of the expression of the two key genes regulating self-pollen recognition in cauliflower, the S-locus receptor kinase (SRK) and S-locus cysteine-rich (SCR/SP11) genes, revealed that self-compatibility or PSC was associated with decreased SRK or SCR/SP11 expression. Our work shows the particularly high level of phenotypic plasticity of the SI response associated with certain S-haplotypes in cauliflower.

KATANIN-dependent mechanical properties of the stigmatic cell wall mediate the pollen tube path in Arabidopsis.

Successful fertilization in angiosperms depends on the proper trajectory of pollen tubes through the pistil tissues to reach the ovules. Pollen tubes first grow within the cell wall of the papilla cells, applying pressure to the cell. Mechanical forces are known to play a major role in plant cell shape by controlling the orientation of cortical microtubules (CMTs), which in turn mediate deposition of cellulose microfibrils (CMFs). Here, by combining imaging, genetic and chemical approaches, we show that isotropic reorientation of CMTs and CMFs in aged Col-0 and katanin1-5 (ktn1-5) papilla cells is accompanied by a tendency of pollen tubes to coil around the papillae. We show that this coiled phenotype is associated with specific mechanical properties of the cell walls that provide less resistance to pollen tube growth. Our results reveal an unexpected role for KTN1 in pollen tube guidance on the stigma by ensuring mechanical anisotropy of the papilla cell wall.

Flowering plants produce small particles known as pollen that ­ with the help of the wind, bees and other animals ­ carry male sex cells (sperm) to female sex cells (eggs) contained within flowers. When a grain of pollen lands on the female organ of a flower, called the pistil, it gives rise to a tube that grows through the pistil towards the egg cells at the base. The surface of the pistil is covered in a layer of long cells named papillae. Like most plant cells, the papillae are surrounded by a rigid structure known as the cell wall, which is mainly composed of strands known as microfibrils. The pollen tube exerts pressure on a papilla to allow it to grow through the cell wall towards the base of the pistil. Previous studies have shown that the pistil produces signals that guide pollen tubes to the eggs. However, it remains unclear how pollen tubes orient themselves on the surface of papillae to grow in the right direction through the pistil. Riglet et al. combined microscopy, genetic and chemical approaches to study how pollen tubes grow through the surface of the pistils of a small weed known as Arabidopsis thaliana. The experiments showed that an enzyme called KATANIN conferred mechanical properties to the cell walls of papillae that allowed pollen tubes to grow towards the egg cells, and also altered the orientation of the microfibrils in these cell walls. In A. thaliana plants that were genetically modified to lack KATANIN the pollen tubes coiled around the papillae and sometimes grew in the opposite direction to where the eggs were. KATANIN is known to cut structural filaments inside the cells of plants, animals and most other living things. By revealing an additional role for KATANIN in regulating the mechanical properties of the papilla cell wall, these findings indicate this enzyme may also regulate the mechanical properties of cells involved in other biological processes.

Combining genetic crosses and pool targeted DNA-seq for untangling genomic variations associated with resistance to multiple insecticides in the mosquito Aedes aegypti.

In addition to combating vector-borne diseases, studying the adaptation of mosquitoes to insecticides provides a remarkable example of evolution-in-action driving the selection of complex phenotypes. Actually, most resistant mosquito populations show multi-resistance phenotypes as a consequence of the variety of insecticides employed and of the complexity of selected resistance mechanisms. Such complexity makes the identification of alleles conferring resistance to specific insecticides challenging and prevents the development of molecular assays to track them in the field. Here we showed that combining simple genetic crosses with pool targeted DNA-seq can enhance the specificity of resistance allele's detection while maintaining experimental work and sequencing effort at reasonable levels. A multi-resistant population of the mosquito Aedes aegypti was exposed to three distinct insecticides (deltamethrin, bendiocarb and fenitrothion), and survivors to each insecticide were crossed with a susceptible strain to generate three distinct lines. F2 individuals from each line were then segregated based on their survival to two insecticide doses. Hundreds of genes covering all detoxifying enzymes and insecticide targets together with more than 7,000 intergenic regions equally spread over mosquito genome were sequenced from pools of F0 and F2 individuals unexposed or surviving insecticide. Differential coverage analysis identified 39 detoxification enzymes showing an increased gene copy number in association with resistance. Combining an allele frequency filtering approach with a Bayesian F ST-based genome scan identified multiple genomic regions showing strong selection signatures together with 50 nonsynonymous variations associated with resistance. This study provides a simple and cost-effective approach to improve the specificity of resistance allele's detection in multi-resistant populations while reducing false positives frequently arising when comparing populations showing divergent genetic backgrounds. The identification of novel DNA resistance markers opens new opportunities for improving the tracking of insecticide resistance in the field.

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.