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

Identifying a potentially invasive population in the native range of a species: The enlightenment from the phylogeography of the yellow spotted stink bug, Erthesina fullo (Hemiptera: Pentatomidae)

The yellow spotted stink bug (YSSB), Erthesina fullo (Thunberg, 1783) is an important Asian pest that has recently successfully invaded Europe and an excellent material for research on the initial stage of biological invasion. Here, we reported the native evolutionary history, recent invasion history, and potential invasion threats of YSSB for the first time based on population genetic methods [using double digest restriction-site associated DNA (ddRAD) data and mitochondrial COI and CYTB] and ecological niche modelling. The results showed that four lineages (east, west, southwest, and Hainan Island) were established in the native range with a strong east-west differentiation phylogeographical structure, and the violent climate fluctuation might cause population divergence during the Middle and Upper Pleistocene. In addition, land bridges and monsoon promote dispersal and directional genetic exchanging between island populations and neighboring continental populations. The east lineage (EA) was identified as the source of invasion in Albania. EA had the widest geographical distribution among all other lineages, with a star-like haplotype network with the main haplotype as the core. It also had a rapid population expansion history, indicating that the source lineage might have stronger diffusion ability and adaptability. Our findings provided a significant biological basis for fine tracking of invasive source at the lineage or population level and promote early invasion warning of potential invasive species on a much subtler lineage level.

Testing for adaptive changes linked to range expansion following a single introduction of the fall webworm

Adaptive evolution following colonization can affect the impact of invasive species. The fall webworm (FWW) invaded China 40 years ago through a single introduction event involving a severe bottleneck and subsequently diverged into two genetic groups. The well-recorded invasion history of FWW, coupled with a clear pattern of genetic divergence, provides an opportunity to investigate whether there is any sign of adaptive evolution following the invasion. Based on genome-wide SNPs, we identified genetically separated western and eastern groups of FWW and correlated spatial variation in SNPs with geographical and climatic factors. Geographical factors explained a similar proportion of the genetic variation across all populations compared with climatic factors. However, when the two population groups were analysed separately, environmental factors explained more variation than geographical factors. SNP outliers in populations of the western group had relatively stronger response to precipitation than temperature-related variables. Functional annotation of SNP outliers identified genes associated with insect cuticle protein potentially related to desiccation adaptation in the western group and genes associated with lipase biosynthesis potentially related to temperature adaptation in the eastern group. Our study suggests that invasive species may maintain the evolutionary potential to adapt to heterogeneous environments despite a single invasion event. The molecular data suggest that quantitative trait comparisons across environments would be worthwhile.

Molecular characterization and functional analysis of the ecdysone receptor isoform (EcR) from the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae)

20-Hydroxyecdysone (20E) plays a vital role in a series of biological processes, via the nuclear receptors, EcR/USP by activating the ecdysone regulatory cascade. To clarify the role of EcR during the development of Grapholita molesta, the complementary DNA of ecdysone receptor isoform B1 (GmEcR-B1) was obtained from the transcriptome of G. molesta and verified by PCR. Alignment analysis revealed that the deduced protein sequence of GmEcR-B1 was highly homologous to EcR proteins identified in other lepidopteran species, especially the EcR-B1 isoform in Spodoptera litura. Quantitative real-time PCR showed that GmEcRs was expressed at all test developmental stages, and the expression level of GmEcRs was relatively higher during the period of the 3rd day of fifth instar larvae to 2nd of pupa than those in other stages. Moreover, the messenger RNA of GmEcRs was much more strongly expressed in the Malpighian tubule and epidermis than those in other tissues, which suggests that this gene may function in a tissue-specific manner during larval development. Silencing of GmEcRs could significantly downregulate the transcriptional level of ecdysone-inducible genes and result in increased mortality during metamorphosis and prolonged prepupal duration. Taken together, the present results indicate that GmEcRs may directly or indirectly affect the development of G. molesta.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

Genomes of the cosmopolitan fruit pest Bactrocera dorsalis (Diptera: Tephritidae) reveal its global invasion history and thermal adaptation

INTRODUCTION

The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural pests worldwide, with highly debated species delimitation, origin, and global spread routes.

OBJECTIVES

Our study intended to (i) resolve the taxonomic uncertainties between B. dorsalis and B. carambolae, (ii) reveal the population structure and global invasion routes of B. dorsalis across Asia, Africa, and Oceania, and (iii) identify genomic regions that are responsible for the thermal adaptation of B. dorsalis.

METHODS

Based on a high-quality chromosome-level reference genome assembly, we explored the population relationship using a genome-scale single nucleotide polymorphism dataset generated from the resequencing data of 487 B. dorsalis genomes and 25 B. carambolae genomes. Genome-wide association studies and silencing using RNA interference were used to identify and verify the candidate genes associated with extreme thermal stress.

RESULTS

We showed that B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide: (i) from Northern India to Northern Southeast Asia, then to Southern Southeast Asia; (ii) from Northern India to Northern Southeast Asian, then to China and Hawaii; and (iii) from Southern India toward the African mainland, then to Madagascar, which is mainly facilitated by human activities including trade and immigration. Twenty-seven genes were identified by a genome-wide association study to be associated with 11 temperature bioclimatic variables. The Cyp6a9 gene may enhance the thermal adaptation of B. dorsalis and thus boost its invasion, which tended to be upregulated at a hardening temperature of 38 °C. Functional verification using RNA interference silencing against Cyp6a9, led to the specific decrease in Cyp6a9 expression, reducing the survival rate of dsRNA-feeding larvae exposed to extreme thermal stress of 45 °C after heat hardening treatments in B. dorsalis.

CONCLUSION

This study provides insights into the evolutionary history and genetic basis of temperature adaptation in B. dorsalis.

The impact of habitat loss on molecular signatures of coevolution between an iconic butterfly (Alcon blue) and its host plant (Marsh gentian)

Habitat loss is threatening natural communities worldwide. Small and isolated populations suffer from inbreeding and genetic drift, which jeopardize their long-term survival and adaptive capacities. However, the consequences of habitat loss for reciprocal coevolutionary interactions remain poorly studied. In this study, we investigated the effects of decreasing habitat patch size and connectivity associated with habitat loss on molecular signatures of coevolution in the Alcon blue butterfly (Phengaris alcon) and its most limited host, the marsh gentian (Gentiana pneumonanthe). Because reciprocal coevolution is characterized by negative frequency-dependent selection as a particular type of balancing selection, we investigated how signatures of balancing selection vary along a gradient of patch size and connectivity, using single nucleotide polymorphisms (SNPs). We found that signatures of coevolution were unaffected by patch characteristics in the host plants. On the other hand, more pronounced signatures of coevolution were observed in both spatially isolated and in large Alcon populations, together with pronounced spatial variation in SNPs that are putatively involved in coevolution. These findings suggest that habitat loss can facilitate coevolution in large butterfly populations through limiting swamping of locally beneficial alleles by maladaptive ones. We also found that allelic richness (Ar) of the coevolutionary SNPs is decoupled from neutral Ar in the butterfly, indicating that habitat loss has different effects on coevolutionary as compared with neutral processes. We conclude that this specialized coevolutionary system requires particular conservation interventions aiming at generating a spatial mosaic of both connected and of isolated habitat to maintain coevolutionary dynamics.

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.

Genetic analyses reveal regional structure and demographic expansion of the predominant tea pest Empoasca onukii (Hemiptera: Cicadellidae) in China

BACKGROUND

The tea green leafhopper, Empoasca onukii Matsuda, is the most destructive insect pest of tea plantations in East Asia. Despite its economic importance and previous studies on this species, it remains unclear as to how this small-sized pest can have such wide range.

RESULTS

By sequencing three mitochondrial genes and 17 microsatellite loci, we revealed the regional structure and demographic expansion of 59 E. onukii populations in China. Bayesian analysis of population genetic structure (BAPS) on microsatellites identified four genetic groups with spatial discontinuities, while analysis on mitochondrial genes inferred five nested and differentiated clusters. Both the Mantel test and the generalized linear model indicated a significant pattern of isolation by geographic distance in E. onukii populations. Based on the approximate Bayesian computation approach, E. onukii was found to have originated from southwestern China and expanded northward and eastward. While MIGRATE-N and Bayesian stochastic search variable selection (BSSVS) procedure in BEAST confirmed the possible eastward and northward dispersal from Yunnan, they also detected more gene flow from the derived populations in central and southeastern China.

CONCLUSION

Our results suggest that the current distribution and structure of E. onukii is complicatedly influenced by human activities of cultivation, wide dissemination of tea in ancient China as well as recent transportation of tea seedlings for establishing new tea plantations. Insights into genetic differentiation and demographic expansion patterns from this study provide an important basis for the development of area-wide management of the E. onukii populations. © 2022 Society of Chemical Industry.

The Role of Genetic Factors in the Differential Invasion Success of Two Spartina Species in China

Biological invasions have become one of the greatest threats to global biodiversity and ecosystem conservation. Most previous studies have revealed how successful invasive species adapt to new environments and climate change through phenotypic and genetic evolution. Some researchers suggested that understanding unsuccessful or less successful biological invasions might be important for understanding the relationships between invasion adaptability and climate factors. We compared the sexual reproduction ability, genetic diversity, and gene × environment interaction in two intentionally introduced alien species in China (Spartina anglica and Spartina alterniflora) based on restriction site-associated DNA (RAD) sequencing. After more than 50 years, the distribution of S. alterniflora has rapidly expanded, while S. anglica has experienced extreme dieback. A total of 212,939 single nucleotide polymorphisms (SNPs) for the two Spartina species were used for analysis. The multilocus genotype (MLG) analysis revealed that clonal reproduction was the prevalent mode of reproduction in both species, indicating that a change in the mode of reproduction was not the key factor enabling successful invasion by Spartina. All genetic diversity indicators (He, Ho, π) in S. alterniflora populations were at least two times higher than those in S. anglica populations, respectively (p < 0.001). Furthermore, the population genetic structure and stronger patterns of climate-associated loci provided support for rapid adaptive evolution in the populations of S. alterniflora in China. Altogether, our results highlight the importance of genetic diversity and local adaptation, which were driven by multiple source populations, in increasing the invasiveness of S. alterniflora.

Population genomic signatures of the oriental fruit moth related to the Pleistocene climates

The Quaternary climatic oscillations are expected to have had strong impacts on the evolution of species. Although legacies of the Quaternary climates on population processes have been widely identified in diverse groups of species, adaptive genetic changes shaped during the Quaternary have been harder to decipher. Here, we assembled a chromosome-level genome of the oriental fruit moth and compared genomic variation among refugial and colonized populations of this species that diverged in the Pleistocene. High genomic diversity was maintained in refugial populations. Demographic analysis showed that the effective population size of refugial populations declined during the penultimate glacial maximum (PGM) but remained stable during the last glacial maximum (LGM), indicating a strong impact of the PGM rather than the LGM on this pest species. Genome scans identified one chromosomal inversion and a mutation of the circadian gene Clk on the neo-Z chromosome potentially related to the endemicity of a refugial population. In the colonized populations, genes in pathways of energy metabolism and wing development showed signatures of selection. These different genomic signatures of refugial and colonized populations point to multiple impacts of Quaternary climates on adaptation in an extant species.

Population Genetic Structure of the Invasive Spotted Alfalfa Aphid Therioaphis trifolii (Hemiptera: Aphididae) in China Inferred From Complete Mitochondrial Genomes

Biological invasions represent a natural rapid evolutionary process in which invasive species may present a major threat to biodiversity and ecosystem integrity. Analyzing the genetic structure and demographic history of invaded populations is critical for the effective management of invasive species. The spotted alfalfa aphid (SAA) Therioaphis trifolii is indigenous in the Mediterranean region of Europe and Africa and has invaded China, causing severe damages to the alfalfa industry. However, little is known about its genetic structure and invasion history. In this study, we obtained 167 complete mitochondrial genome sequences from 23 SAA populations across China based on high-throughput sequencing and performed population genetic and phylogenomic analyses. High haplotype diversity and low nucleotide diversity were found in SAA populations in China with distinct genetic structures, i.e., all populations diverged into three phylogenetic lineages. Demographic history analyses showed a recent expansion of the SAA population, consistent with the recent invasion history. Our study indicated that SAA may have invaded through multiple introduction events during commercial trades of alfalfa, although this needs further validation by nuclear markers.