Rapid and strong population genetic differentiation and genomic signatures of climatic adaptation in an invasive mealybug

Genomic investigations of successful invasions: the picture emerging from recent studies

Invasion biology aims to identify traits and mechanisms that contribute to successful invasions, while also providing general insights into the mechanisms underlying population expansion and adaptation to rapid climate and habitat changes. Certain phenotypic attributes have been linked to successful invasions, and the role of genetics has been critical in understanding adaptation of invasive species. Nevertheless, a comprehensive summary evaluating the most common evolutionary mechanisms associated with successful invasions across species and environments is still lacking. Here we present a systematic review of studies since 2015 that have applied genomic tools to investigate mechanisms of successful invasions across different organisms. We examine demographic patterns such as changes in genomic diversity at the population level, the presence of genetic bottlenecks and gene flow in the invasive range. We review mechanisms of adaptation such as selection from standing genetic variation and de novo mutations, hybridisation and introgression, all of which can have an impact on invasion success. This comprehensive review of recent articles on the genomic diversity of invasive species led to the creation of a searchable database to provide researchers with an accessible resource. Analysis of this database allowed quantitative assessment of demographic and adaptive mechanisms acting in invasive species. A predominant role of admixture in increasing levels of genetic diversity enabling molecular adaptation in novel habitats is the most important finding of our study. The "genetic paradox" of invasive species was not validated in genomic data across species and ecosystems. Even though the presence of genetic drift and bottlenecks is commonly reported upon invasion, a large reduction in genomic diversity is rarely observed. Any decrease in genetic diversity is often relatively mild and almost always restored via gene flow between different invasive populations. The fact that loci under selection are frequently detected suggests that adaptation to novel habitats on a molecular level is not hindered. The above findings are confirmed herein for the first time in a semi-quantitative manner by molecular data. We also point to gaps and potential improvements in the design of studies of mechanisms driving rapid molecular adaptation in invasive populations. These include the scarcity of comprehensive studies that include sampling from multiple native and invasive populations, identification of invasion sources, longitudinal population sampling, and the integration of fitness measures into genomic analyses. We also note that the potential of whole genome studies is often not exploited fully in predicting invasive potential. Comparative genomic studies identifying genome features promoting invasions are underrepresented despite their potential for use as a tool in invasive species control.

Limited population structure but signals of recent selection in introduced African Fig Fly (Zaprionus indianus) in North America

Invasive species have devastating consequences for human health, food security, and the environment. Many invasive species adapt to new ecological niches following invasion, but little is known about the early steps of adaptation. Here we examine population genomics of a recently introduced drosophilid in North America, the African Fig Fly, Zaprionus indianus. This species is likely intolerant of subfreezing temperatures and recolonizes temperate environments yearly. We generated a new chromosome-level genome assembly for Z. indianus. Using resequencing of over 200 North American individuals collected over four years in temperate Virginia, plus a single collection from subtropical Florida, we tested for signatures of recolonization, population structure, and adaptation within invasive populations. We show founding populations are sometimes small and contain close genetic relatives, yet temporal population structure and differentiation of populations is mostly absent across recurrent recolonization events. Although we find limited signals of genome-wide spatial or temporal population structure, we identify haplotypes on the X chromosome that are repeatedly differentiated between Virginia and Florida populations. These haplotypes show signatures of natural selection and are not found in African populations. We also find evidence for several large structural polymorphisms segregating within North America populations and show X chromosome evolution in invasive populations is strikingly different from the autosomes. These results show that despite limited population structure, populations may rapidly evolve genetic differences early in an invasion. Further uncovering how these genomic regions influence invasive potential and success in new environments will advance our understanding of how organisms evolve in changing environments.

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Climatic niche and range shifts of grey squirrels (Sciurus carolinensis Gmelin) in Europe: An invasive pest displacing native squirrels

BACKGROUND

As an invasive pest from North America, grey squirrels (GSs; Sciurus carolinensis Gmelin) are displacing native squirrels in Europe. However, the climatic niche and range dynamics of GSs in Europe remain largely unknown. Through niche and range dynamic models, we investigated climatic niche and range shifts between introduced GSs in Europe and native GSs in North America.

RESULTS

GSs in North America can survive in more variable climatic conditions and have much wider climatic niche breadth than do GSs in Europe. Based on climate, the potential range of GSs in Europe included primarily Britain, Ireland, and Italy, whereas the potential range of GSs in North America included vast regions of western and southern Europe. If GSs in Europe could occupy the same climatic niche space and potential range as GSs in North America, they would occupy an area ca. 2.45 times the size of their current range. The unfilling ranges of GSs in Europe relative to those of GSs in North America were primarily in France, Italy, Spain, Croatia, and Portugal.

CONCLUSION

Our observations implied that GSs in Europe have significant invasion potential, and that range projections based on their occurrence records in Europe may underestimate their invasion risk. Given that small niche shifts between GSs in Europe and in North America could lead to large range shifts, niche shifts could be a sensitive indicator in invasion risk assessment. The identified unfilling ranges of the GS in Europe should be prioritized in combating GS invasions in the future. © 2023 Society of Chemical Industry.

Chronosequence of invasion reveals minimal losses of population genomic diversity, niche expansion, and trait divergence in the polyploid, leafy spurge

Rapid evolution may play an important role in the range expansion of invasive species and modify forecasts of invasion, which are the backbone of land management strategies. However, losses of genetic variation associated with colonization bottlenecks may constrain trait and niche divergence at leading range edges, thereby impacting management decisions that anticipate future range expansion. The spatial and temporal scales over which adaptation contributes to invasion dynamics remain unresolved. We leveraged detailed records of the ~130-year invasion history of the invasive polyploid plant, leafy spurge (Euphorbia virgata), across ~500 km in Minnesota, U.S.A. We examined the consequences of range expansion for population genomic diversity, niche breadth, and the evolution of germination behavior. Using genotyping-by-sequencing, we found some population structure in the range core, where introduction occurred, but panmixia among all other populations. Range expansion was accompanied by only modest losses in sequence diversity, with small, isolated populations at the leading edge harboring similar levels of diversity to those in the range core. The climatic niche expanded during most of the range expansion, and the niche of the range core was largely non-overlapping with the invasion front. Ecological niche models indicated that mean temperature of the warmest quarter was the strongest determinant of habitat suitability and that populations at the leading edge had the lowest habitat suitability. Guided by these findings, we tested for rapid evolution in germination behavior over the time course of range expansion using a common garden experiment and temperature manipulations. Germination behavior diverged from the early to late phases of the invasion, with populations from later phases having higher dormancy at lower temperatures. Our results suggest that trait evolution may have contributed to niche expansion during invasion and that distribution models, which inform future management planning, may underestimate invasion potential without accounting for evolution.

Niche and Range Shifts of the Fall Webworm (Hyphantria cunea Dury) in Europe Imply Its Huge Invasion Potential in the Future

The fall webworm (Hyphantria cunea Dury) has a strong impact on agricultural systems in Europe. However, its invasive potential, which was inherited from its native niche in North America, remains unknown. Here, we investigated the climatic niche and range shifts of the fall webworm in Europe and compared them with those in native North America, then assessed the worms' invasive potential in Europe. Compared with the fall webworm in Europe, those in North America survived in more diverse climatic conditions, which was closely associated with their broader niche and larger potential ranges in Europe. If the fall webworm in Europe could exploit the native niche inherited from those in North America to adapt to climatic conditions in Europe, their potential ranges in Europe could be 5.5-fold those based on the niche as introduced in Europe. The potentially unfilled ranges of the fall webworm in Europe were mainly detected in vast regions of Europe, excluding Norway, Sweden, Finland, North Russia, Hungary, Croatia, Romania, and Ukraine, suggesting that, without strict control, these vast regions might be preferably invaded by the fall webworm in Europe in the future. Therefore, strict control against its invasion is needed. Given that small niche shifts in this invasive insect could result in large range shifts, the niche shifts represent a more sensitive indicator of invasion risk than range shifts.

Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects

Hypogeococcus pungens is a species complex native to southern South America that is composed of at least five putative species, each one specialized in the use of different host plants. Two of these undescribed species were registered as invasive in Central and North America: Hyp-C is a cactophagous mealybug that became an important pest that threatens endemic cactus species in Puerto Rico, and Hyp-AP feeds on Amaranthaceae and Portulacaceae hosts, but does not produce severe damage to the host plants. We quantified genomic variation and investigated the demographic history of both invasive species by means of coalescent-based simulations using high throughput sequencing data. We also evaluated the incidence of host plant infestation produced by both species and used an ecological niche modeling approach to assess potential distribution under current and future climatic scenarios. Our genetic survey evinced the footprints of strong effective population size reduction and signals of genetic differentiation among populations within each species. Incidence of plant attacks varied between species and among populations within species, with some host plant species preferred over others. Ecological niche modeling suggested that under future climatic scenarios both species would expand their distribution ranges in Puerto Rico. These results provide valuable information for the design of efficient management and control strategies of the Puerto Rican cactus pest and shed light on the evolutionary pathways of biological invasions.

Unraveling the hierarchical genetic structure of tea green leafhopper, Matsumurasca onukii, in East Asia based on SSRs and SNPs

Matsumurasca onukii (Matsuda, R. (1952). Oyo-Kontyu Tokyo, 8(1): 19-21), one of the dominant pests in major tea production areas in Asia, currently is known to occur in Japan, Vietnam, and China, and severely threatens tea production, quality, and international trade. To elucidate the population genetic structure of this species, 1633 single nucleotide polymorphisms (SNPs) and 18 microsatellite markers (SSRs) were used to genotype samples from 27 sites representing 18 geographical populations distributed throughout the known range of the species in East Asia. Analyses of both SNPs and SSRs showed that M. onukii populations in Yunnan exhibit high-genetic differentiation and structure compared with the other populations. The Kagoshima (JJ) and Shizuoka (JS) populations from Japan were separated from populations from China by SNPs, but clustered with populations from Jinhua (JH), Yingde (YD), Guilin (GL), Fuzhou (FZ), Hainan (HQ), Leshan (CT), Chongqing (CY), and Zunyi (ZY) tea plantations in China and the Vietnamese Vinh Phuc (VN) population based on the SSR data. In contrast, CT, CY, ZY, and Shaanxi (SX) populations clustered together based on SNPs, but were separated by SSRs. Both marker datasets identified significant geographic differentiation among the 18 populations. Various environmental and anthropogenic factors, including geographical barriers to migration, human transport of hosts (Camellia sinesis [L.] O. Kuntze) and adaptation of M. onukii to various local climatic zones possibly account for the rapid spread of this pest in Asia. The results demonstrate that SNPs from high-throughput genotyping data can be used to reveal subtle genetic substructure at broad scales in r-strategist insects.

Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops

Weevils, classified in the family Curculionidae (true weevils), constitute a group of phytophagous insects of which many species are considered significant pests of crops. Within this family, the red palm weevil (RPW), Rhynchophorus ferrugineus, has an integral role in destroying crops and has invaded all countries of the Middle East and many in North Africa, Southern Europe, Southeast Asia, Oceania, and the Caribbean Islands. Simple sequence repeats (SSRs), also termed microsatellites, have become the DNA marker technology most applied to study population structure, evolution, and genetic diversity. Although these markers have been widely examined in many mammalian and plant species, and draft genome assemblies are available for many species of true weevils, very little is yet known about SSRs in weevil genomes. Here we carried out a comparative analysis examining and comparing the relative abundance, relative density, and GC content of SSRs in previously sequenced draft genomes of nine true weevils, with an emphasis on R. ferrugineus. We also used Illumina paired-end sequencing to generate draft sequence for adult female RPW and characterized it in terms of perfect SSRs with 1–6 bp nucleotide motifs. Among weevil genomes, mono- to trinucleotide SSRs were the most frequent, and mono-, di-, and hexanucleotide SSRs exhibited the highest GC content. In these draft genomes, SSR number and genome size were significantly correlated. This work will aid our understanding of the genome architecture and evolution of Curculionidae weevils and facilitate exploring SSR molecular marker development in these species.

Biology, ecology, and management of cotton mealybug Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae)

The cotton mealybug Phenacoccus solenopsis Tinsley (Hemiptera: Pseudococcidae), is a highly invasive and harmful pest. It causes considerable loss of cotton crops in China, India and Pakistan. Little is known about its bionomics since it was first recorded in Pakistan and India in 2005. Rapid spread of this pest worldwide has accelerated research on its biology, ecology and management. The P. solenopsis has a short life cycle, and optimal temperatures lead to an increase in the number of generations per year, which is a serious threat to cotton crop production. Cotton mealybug is native to the USA, although it has now spread to >43 countries. Insecticidal control is the primary and dominant practice for this pest, and its resistance to commonly used insecticides is increasing. Biocontrol agents have strong potential for the management of nymphal instar stages. We read >250 articles related to our review title and finally reviewed recent advances in the understanding of P. solenopsis biology, ecology and control approaches, aiming to highlight integrated and biological management practices of this pest. © 2021 Society of Chemical Industry.

Local climate adaptation and gene flow in the native range of two co-occurring fruit moths with contrasting invasiveness

Invasive species pose increasing threats to global biodiversity and ecosystems. While previous studies have characterized successful invaders based on ecological traits, characteristics related to evolutionary processes have rarely been investigated. Here we compared gene flow and local adaptation using demographic analyses and outlier tests in two co-occurring moth pests across their common native range of China, one of which (the peach fruit moth, Carposina sasakii) has maintained its native distribution, while the other (the oriental fruit moth, Grapholita molesta) has expanded its range globally during the past century. We found that both species showed a pattern of genetic differentiation and an evolutionary history consistent with a common southwestern origin and northward expansion in their native range. However, for the noninvasive species, genetic differentiation was closely aligned with the environment, and there was a relatively low level of gene flow, whereas in the invasive species, genetic differentiation was associated with geography. Genome scans indicated stronger patterns of climate-associated loci in the noninvasive species. While strong local adaptation and reduced gene flow across its native range may have decreased the invasiveness of C. sasakii, this requires further validation with additional comparisons of invasive and noninvasive species across their native range.

Genomic signatures of admixture and selection are shared among populations of Zaprionus indianus across the western hemisphere

Introduced species have become an increasingly common component of biological communities around the world. A central goal in invasion biology is therefore to identify the demographic and evolutionary factors that underlie successful introductions. Here we use whole genome sequences, collected from populations in the native and introduced range of the African fig fly, Zaprionus indianus, to quantify genetic relationships among them, identify potential sources of the introductions, and test for selection at different spatial scales. We find that geographically widespread populations in the western hemisphere are genetically more similar to each other than to lineages sampled across Africa, and that these populations share a mixture of alleles derived from differentiated African lineages. Using patterns of allele‐sharing and demographic modelling we show that Z. indinaus have undergone a single expansion across the western hemisphere with admixture between African lineages predating this expansion. We also find support for selection that is shared across populations in the western hemisphere, and in some cases, with a subset of African populations. This suggests either that parallel selection has acted across a large part of Z. indianus's introduced range; or, more parsimoniously, that Z. indianus has experienced selection early on during (or prior‐to) its expansion into the western hemisphere. We suggest that the range expansion of Z. indianus has been facilitated by admixture and selection, and that management of this invasion could focus on minimizing future admixture by controlling the movement of individuals within this region rather than between the western and eastern hemisphere.

Ecological phylogeography and coalescent models suggest a linear population expansion of Anastrepha fraterculus (Diptera: Tephritidae) in southern South America

This work is a first approach to an integrated view of the genetics, ecology and dispersion patterns of Anastrepha fraterculus in southern South America. We studied the association of genetic variation with geographical patterns and environmental variables to provide insight into the crucial factors that drive the structure and dynamics of fly populations. Data from a 417 bp mitochondrial COII gene fragment from seven Argentinian populations and one South Brazilian population (from five ecoregions grouped in three biomes) were used to identify population clusters using a model-based Bayesian phylogeographical and ecological clustering approach. The sequences were also analysed under a coalescent model to evaluate historical demographic changes. We identified 19 different haplotypes and two clusters differing in all the environmental covariables. The assumption of neutral evolution and constant population size was rejected, and the population growth parameters suggested a linear population expansion starting 2500 years before present. The most likely ancestral location is Posadas, from where A. fraterculus would have expanded southwards and westwards in Argentina. This result is consistent with Holocene changes and anthropic factors related to the expansion of the Tupí–Guaraní culture, 3000–1500 years before present.

Comparison of microsatellite distribution patterns in twenty-nine beetle genomes

Simple sequence repeats (SSRs) represent an important source of genetic variation that provides a basis for adaptation to different environments in organisms. In this study, we examined the distribution patterns of SSRs in twenty-nine beetle genomes and carried out Gene Ontology (GO) analysis of CDSs embedded with perfect SSRs (P-SSRs). The results demonstrated that imperfect SSRs (I-SSRs) represented the most abundant SSR category in beetle genomes and in different genomic regions (CDS, exon, and intron regions). The numbers of P-SSRs, I-SSRs, compound SSRs, and variable number tandem repeats were positively correlated with beetle genome size, whereas neither the frequency nor the density of the SSRs was correlated with genome size. Moreover, our results demonstrated that common genomic features of P-SSRs within the same suborder or family of Coleoptera were rare. Mono-, di-, tri-, or tetranucleotide SSRs were the most abundant P-SSR categories in beetle genomes. The preferred predominant repeat motif among the mononucleotide P-SSRs was (A)n, but the most frequent repeat motifs for other length classes varied differentially among these genomes. Furthermore, the P-SSR type with the highest GC content differed in the beetle genomes and in different genomic regions. CV (coefficient of variability) analysis demonstrated that the repeat copy numbers of P-SSRs presented relatively higher variation in introns than in CDSs and exons. The GO terms of CDSs containing P-SSRs for molecular functions were mainly enriched in "binding" and "transcription". Our findings will be useful for studying the functional roles of microsatellite heterogeneity in beetle adaptation.