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.

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.

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.

Low levels of genetic differentiation with isolation by geography and environment in populations of Drosophila melanogaster from across China.

The fruit fly, Drosophila melanogaster, is a model species in evolutionary studies. However, population processes of this species in East Asia are poorly studied. Here we examined the population genetic structure of D. melanogaster across China. There were 14 mitochondrial haplotypes with 10 unique ones out of 23 known from around the globe. Pairwise FST values estimated from 15 novel microsatellites ranged from 0 to 0.11, with geographically isolated populations showing the highest level of genetic uniqueness. STRUCTURE analysis identified high levels of admixture at both the individual and population levels. Mantel tests indicated a strong association between genetic distance and geographical distance as well as environmental distance. Full redundancy analysis (RDA) showed that independent effects of environmental conditions and geography accounted for 62.10% and 31.58% of the total explained genetic variance, respectively. When geographic variables were constrained in a partial RDA analysis, the environmental variables bio2 (mean diurnal air temperature range), bio13 (precipitation of the wettest month), and bio15 (precipitation seasonality) were correlated with genetic distance. Our study suggests that demographic history, geographical isolation, and environmental factors have together shaped the population genetic structure of D. melanogaster after its introduction into China.

Field-evolved resistance and cross-resistance of the two-spotted spider mite, Tetranychus urticae, to bifenazate, cyenopyrafen and SYP-9625.

The acaricide bifenazate acts as complex III inhibitor whereas cyenopyrafen and SYP-9625 act as complex II inhibitors. All these acaricides are commonly used to control two-spotted spider mite (TSSM), Tetranychus urticae Koch. We examined field-evolved and laboratory-selected resistance of TSSM to these three acaricides and determined cross-resistance among them. Six field populations of TSSM showed low levels of resistance to bifenazate with resistance ratios ranging from 2.20 to 10.65 compared to a susceptible strain. SYP-9625, structurally similar to cyenopyrafen, showed slightly higher activity to TSSMs but significant cross-resistance in both field populations and a laboratory-selected strain by SYP-9625. However, low levels of resistance to these two chemicals were found in field populations even when used for short time periods. Cross-resistance was not found between bifenazate and Complex II inhibitors, cyenopyrafen and SYP-9625, in both field populations and the laboratory-selected strain. Field-evolved resistance of TSSM to the tested acaricides is still low and should be delayed by the implementation of resistance management practices. Cross-resistance between cyenopyrafen and SYP-9625 is obvious, so they should not be used together in resistance management strategies based on mode of action rotation.

Toxicity of seven insecticides to different developmental stages of the whitefly Bemisia tabaci MED (Hemiptera: Aleyrodidae) in multiple field populations of China.

Chemical control is important in the management of the tobacco whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae). Susceptibility of B. tabaci to insecticides may vary among different developmental stages and geographical populations. In this study, we examined toxicity of seven commonly-used insecticides to B. tabaci MED in four field populations from China. Avermectin has high level of toxicity to all stages of B. tabaci MED in all four populations. Cyantraniliprole and sulfoxaflor have high toxicity to adults. Spirotetramat, cyantraniliprole and flonicamid have high toxicity to nymphs but not adults. Acetamiprid, cyantraniliprole and sulfoxaflor have high toxicity to eggs. However, the relative toxicity of B. tabaci MED to these chemicals varied across different populations, with little consistency in population differences across developmental stages. Our findings together with some instances where LC95 values were higher than field recommended dosages indicate field-evolved resistance to insecticides (such as thiamethoxam and sulfoxaflor) and stage-specific mechanisms that will influence effective control of B. tabaci MED by insecticides.

Preference and performance of the two-spotted spider mite Tetranychus urticae (Acari: Tetranychidae) on strawberry cultivars.

The two-spotted spider mite (TSSM), Tetranychus urticae Koch (Acari: Tetranychidae), is one of the most serious pests of strawberry worldwide. Understanding the preference of TSSM for particular cultivars of strawberry and performance on them helps identify host-plant resistance to this pest mite. In this study, we tested preference, developmental duration, fecundity and population levels of TSSM on 14 strawberry cultivars. TSSM showed strong preference for the Chinese cultivars of Yanxiang, Baixuegongzhu, and Jingtaoxiang. Development of TSSM on the cultivars varied from 32.32 to 36.82 days; it was longest on the cultivars Hongxiutianxiang and Baixuegongzhu, and shortest on Yanxiang, Jingzangxiang, and Darselect as well as on a wild variety (Wuye). TSSM had high fecundity on the cultivars Yanxiang, Taoxun, Hongxiutianxiang, Jingzangxiang, Albion and Baixuegongzhu as well as on Wuye, whereas egg production was lowest on Sweet Charlie, Portola, Akihime, and Benihoppe. After 28 days of plant infestation with 10 pairs of adults, the cultivars Yanxiang, Taoxun, Jingzangxiang, Jingtaoxiang, and Baixuegongzhu had the highest number of mites (> 1000 per plant), whereas mite numbers on Albion and Camarosa were low. The population size of TSSM was correlated with fecundity, but no correlation was found between other preference/performance measures. Our study suggests that a rapid increase of population size of TSSM on cultivars of strawberry is related to high fecundity, and also that there are substantial differences in preference and performance across cultivars.

Efficacy of carbon dioxide treatments for the control of the two-spotted spider mite, Tetranychus urticae, and treatment impact on plant seedlings.

To develop a new control method for the two-spotted spider mite (TSSM), Tetranychus urticae, we investigated the effect of controlled atmospheres of carbon dioxide (CO2) on TSSM mortality under different concentrations and treatment periods, and evaluated the impact of treatments on seedlings of five host plants of TSSM. Egg hatching rate of TSSM was reduced to 37.7, 5.4 or 0% after 24 h treatment involving concentrations of 16.7, 33.3 or 50%, respectively. Mobile stages (nymphs and adult) of TSSM were completely controlled after 24 h treatment at concentrations higher than 33.3%. After 4 h at concentrations of 33.3 or 50%, 1st-day survival rate for all mobile stages was 45.3 or 36.0%, respectively, whereas after 8 or 16 h treatments, all values were decreased to zero. Seedlings of four major host plants of TSSM (cucumber, eggplant, rape, green peppers) were damaged to varying degrees after 24 h at the three concentrations, but strawberry, another host plant, was not damaged. Cucumber suffered the most serious damage, resulting in wilting and death. In conclusion, controlled atmospheres of CO2 can kill TSSM, particularly at high concentrations and with long treatment times. It can be used to control TSSM on strawberry, but should be used cautiously on other host plants.

Nuclear and mitochondrial genomes of the plum fruit moth Grapholita funebrana.

The plum fruit moth Grapholita funebrana (Tortricidae, Lepidoptera) is an important pest of many wild and cultivated stone fruits and other plants in the family Rosaceae. Here, we assembled its nuclear and mitochondrial genomes using Illumina, Nanopore, and Hi-C sequencing technologies. The nuclear genome size is 570.9 Mb, with a repeat rate of 51.28%, and a BUCSO completeness of 97.7%. The karyotype for males is 2n = 56. We identified 17,979 protein-coding genes, 5,643 tRNAs, and 94 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomes provide resources to understand the genetics, ecology, and genome evolution of the tortricid moths.

Chromosome-level genome assembly of the western flower thrips Frankliniella occidentalis.

The western flower thrips Frankliniella occidentalis (Thysanoptera: Thripidae) is a global invasive species that causes increasing damage by direct feeding on crops and transmission of plant viruses. Here, we assemble a previously published scaffold-level genome into a chromosomal level using Hi-C sequencing technology. The assembled genome has a size of 302.58 Mb, with a contig N50 of 1533 bp, scaffold N50 of 19.071 Mb, and BUSCO completeness of 97.8%. All contigs are anchored on 15 chromosomes. A total of 16,312 protein-coding genes are annotated in the genome with a BUSCO completeness of 95.2%. The genome contains 492 non-coding RNA, and 0.41% of interspersed repeats. In conclusion, this high-quality genome provides a convenient and high-quality resource for understanding the ecology, genetics, and evolution of thrips.

Nuclear and mitochondrial genomes of Polylopha cassiicola: the first assembly in Chlidanotinae (Tortricidae).

Tortricidae is one of the largest families in Lepidoptera, including subfamilies of Tortricinae, Olethreutinae, and Chlidanotinae. Here, we assembled the gap-free genome for the subfamily Chlidanotinae using Illumina, Nanopore, and Hi-C sequencing from Polylopha cassiicola, a pest of camphor trees in southern China. The nuclear genome is 302.03 Mb in size, with 36.82% of repeats and 98.4% of BUCSO completeness. The karyotype is 2n = 44 for males. We identified 15412 protein-coding genes, 1052 tRNAs, and 67 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and one rRNA. These high-quality genomes provide valuable information for studying phylogeny, karyotypic evolution, and adaptive evolution of tortricid moths.

Plasticity of cold and heat stress tolerance induced by hardening and acclimation in the melon thrips.

Extreme temperatures threaten species under climate change and can limit range expansions. Many species cope with changing environments through plastic changes. This study tested phenotypic changes in heat and cold tolerance under hardening and acclimation in the melon thrips, Thrips palmi Karny (Thysanoptera: Thripidae), an agricultural pest of many vegetables. We first measured the critical thermal maximum (CTmax) of the species by the knockdown time under static temperatures and found support for an injury accumulation model of heat stress. The inferred knockdown time at 39 °C was 82.22 min. Rapid heat hardening for 1 h at 35 °C slightly increased CTmax by 1.04 min but decreased it following exposure to 31 °C by 3.46 min and 39 °C by 6.78 min. Heat acclimation for 2 and 4 days significantly increased CTmax at 35 °C by 1.83, and 6.83 min, respectively. Rapid cold hardening at 0 °C and 4 °C for 2 h, and cold acclimation at 10 °C for 3 days also significantly increased cold tolerance by 6.09, 5.82, and 2.00 min, respectively, while cold hardening at 8 °C for 2 h and acclimation at 4 °C and 10 °C for 5 days did not change cold stress tolerance. Mortality at 4 °C for 3 and 5 days reached 24.07 % and 43.22 % respectively. Our study showed plasticity for heat and cold stress tolerance in T. palmi, but the thermal and temporal space for heat stress induction is narrower than for cold stress induction.

A chromosome-level genome for the flower thrips Frankliniella intonsa.

The flower thrips Frankliniella intonsa (Thysanoptera: Thripidae) is a common insect found in flowers of many plants. Sometimes, F. intonsa causes damage to crops through direct feeding and transmission of plant viruses. Here, we assembled a chromosomal level genome of F. intonsa using the Illumina, Oxford Nanopore (ONT), and Hi-C technologies. The assembled genome had a size of 209.09 Mb, with a contig N50 of 997 bp, scaffold N50 of 13.415 Mb, and BUSCO completeness of 92.5%. The assembled contigs were anchored on 15 chromosomes. A set of 14,109 protein-coding genes were annotated in the genome with a BUSCO completeness of 95.0%. The genome contained 491 non-coding RNA and 0.57% of interspersed repeats. This high-quality genome provides a valuable resource for understanding the ecology, genetics, and evolution of F. intonsa, as well as for controlling thrips pests.

Chromosome-level genome assembly of the Japanese sawyer beetle Monochamus alternatus.

The Japanese sawyer beetle Monochamus alternatus (Coleoptera: Cerambycidae) is a pest in pine forests and acts as a vector for the pine wood nematode Bursaphelenchus xylophilus, which causes the pine wilt disease. We assembled a high-quality genome of M. alternatus at the chromosomal level using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled genome is 767.12 Mb, with a scaffold N50 of 82.0 Mb. All contigs were assembled into ten pseudo-chromosomes. The genome contains 63.95% repeat sequences. We identify 16, 284 protein-coding genes in the genome, of which 11,244 were functionally annotated. The high-quality genome of M. alternatus provides an invaluable resource for the biological, ecological, and genetic study of this beetle and opens new avenues for understanding the transmission of pine wood nematode by insect vectors.

A comprehensive assessment of insecticide resistance mutations in source and immigrant populations of the diamondback moth Plutella xylostella (L.).

BACKGROUND: The diamondback moth (DBM) Plutella xylostella has developed resistance to almost all insecticides used to control it. Populations of DBM in temperate regions mainly migrate from annual breeding areas. However, the distribution pattern of insecticide resistance of DBM within the context of long-distance migration remains unclear. RESULTS: In this study, we examined the frequency of 14 resistance mutations for 52 populations of DBM collected in 2010, 2011, 2017 and 2018 across China using a high-throughput KASP genotyping method. Mutations L1041F and T929I conferring pyrethroid resistance, and mutations G4946E and E1338D conferring chlorantraniliprole resistance were near fixation in most populations, whereas resistant alleles of F1020S, M918I, A309V and F1845Y were uncommon or absent in most populations. Resistance allele frequencies were relatively stable among different years, although the frequency of two mutations decreased. Principal component analysis based on resistant allele frequencies separated a southern population as an outlier, whereas the immigrants clustered with other populations, congruent with the migration pattern of northern immigrants coming from the Sichuan area of southwestern China. Most resistant mutations deviated from Hardy-Weinberg equilibrium due to a lower than expected frequency of heterozygotes. The deviation index of heterozygosity for resistant alleles was significantly higher than the index obtained from single nucleotide polymorphisms across the genome. These findings suggest heterogeneous selection pressures on resistant mutations. CONCLUSION: Our results provide a picture of resistant mutation patterns in DBM shaped by insecticide usage and migration of this pest, and highlight the widespread distribution of resistance alleles in DBM. © 2022 Society of Chemical Industry.

A high-throughput KASP assay provides insights into the evolution of multiple resistant mutations in populations of the two-spotted spider mite Tetranychus urticae across China.

BACKGROUND: The two-spotted spider mite (TSSM), Tetranychus urticae (Acari: Tetranychidae), is a cosmopolitan phytophagous pest in agriculture and horticulture. It has developed resistance to many acaricides by target-site mutations. Understanding the status and evolution of resistant mutations in the field is essential for resistance management. Here, we applied a high-throughput Kompetitive allele-specific polymerase chain reaction (KASP) method for detecting six mutations conferring resistance to four acaricides of the TSSM. We genotyped 3274 female adults of TSSM from 43 populations collected across China in 2017, 2020, and 2021. RESULTS: The KASP genotyping of 24 testing individuals showed 99% agreement with Sanger sequencing results. KASP assays showed that most populations had a high frequency of mutations conferring avermectin (G314D and G326E) and pyridaben (H92R) resistance. The frequency of mutation conferring bifenazate (A269V and G126S) and etoxazole (I1017F) resistance was relatively low. Multiple mutations were common in the TSSM, with 70.2% and 24.6% of individuals having 2-6 and 7-10 of 10 possible resistant alleles, respectively. No loci were linked in most populations among the six mutations, indicating the development of multiple resistance is mainly by independent selection. However, G314D and I1017F on the nuclear genome deviated from Hardy-Weinberg equilibrium in most populations, indicating significant selective pressure on TSSM populations by acaricides or fitness cost of the mutations in the absence of acaricide selection. CONCLUSION: Our study revealed that the high frequency of TSSMs evolved multiple resistant mutations in population and individual levels by independent selection across China, alarming for managing multiple-acaricides resistance. © 2023 Society of Chemical Industry.

KASP genotyping and semi-quantitation of G275E mutation in the α6 subunit of Thrips palmi nAChR gene conferring spinetoram resistance.

BACKGROUND: Pesticide resistance is a long-standing and growing problem in the chemical control of invertebrate pests. Molecular diagnostic methods can facilitate pesticide resistance management by accurately and efficiently detecting resistant mutations and their frequency. In this study, the kompetitive allele specific PCR (KASP) approach, a technology for high-throughput single nucleotide polymorphism (SNP) genotyping, is validated as a useful method for characterizing genotypes at a pesticide-resistance locus for the first time. We focus on the spinetoram resistance mutation of G275E in the nicotinic acetylcholine receptor alpha 6 (nAChR α6) subunit gene of Thrips palmi. RESULTS: Of the 341 individuals of Thrips palmi tested, 98.24% were successfully genotyped, with 100% concordance with Sanger sequencing results. We then quantitatively mixed genomic DNA of known genotypes to establish 21 DNA mixtures with a resistant allele frequency ranging from 0 to 100% at steps of 5%. The linear discriminant analysis (LDA) showed that 75.8% of original grouped cases were correctly classified; six groups had no overlap in membership (resistant allele frequency: 0%, 5%, 10-75%, 80-85%, 90-95%, and 100%). When we chose 11 pooled samples with 10% steps for LDA, 84.4% of original grouped cases were correctly classified; seven groups had no overlap in membership (0%, 10%, 20-30%, 40-70%, 80%, 90%, 100%). The results indicated that KASP applied to pooled samples may provide a semi-quantitative estimate of resistance. CONCLUSIONS: Our study points to the suitability of KASP for high-throughput genotyping of genotypes affecting pesticide resistance and semi-quantitative assessments of resistance allele frequencies in populations. © 2023 Society of Chemical Industry.

Effects of fungicides on fitness and Buchnera endosymbiont density in Aphis gossypii.

BACKGROUND: Several agricultural fungicides are known to affect insect pests directly and these effects may be transgenerational and mediated through impacts on endosymbionts, providing opportunities for pest control. The cotton aphid Aphis gossypii is a polyphagous pest that can cause large crop yield losses. Here, we tested the effects of three fungicides, pyraclostrobin, trifloxystrobin and chlorothalonil, on the fitness and Buchnera endosymbiont of A. gossypii. RESULTS: The formulations of trifloxystrobin and pyraclostrobin, and the active ingredient of pyraclostrobin produced dose-dependent mortality in A. gossypii, whereas there was no dose-dependent mortality for chlorothalonil. The formulations of trifloxystrobin and pyraclostrobin significantly reduced the lifespan and fecundity of A. gossypii, and increased the density of Buchnera in the parental generation but not the (unexposed) F1 . When the active ingredient of pyraclostrobin was tested, the lifespan of the F0 generation was also reduced, but the density of Buchnera was not, indicating that non-insecticidal chemicals in the fungicide formulation may affect the density of the endosymbiont of A. gossypii. There was no transgenerational effect of the active ingredient of pyraclostrobin on the lifespan and Buchnera of (unexposed) F1 . CONCLUSIONS: Our results suggest that formulations of two strobilurin fungicides have immediate impacts on the fitness of A. gossypii, and chemicals in the formulation impact the density of the primary Buchnera endosymbiont. Our study highlights the potential effects of non-insecticidal chemicals of fungicides on aphid pests and their primary endosymbionts but direct connections between fitness and Buchnera densities remain unclear. © 2023 Society of Chemical Industry.

Chromosome-level genome assembly of an agricultural pest Zeugodacus tau (Diptera: Tephritidae).

The fruit fly Zeugodacus tau (Diptera: Tephritidae) is a major pest of melons and other cucurbits in Southeast Asia. In this study, we used Illumina, Nanopore, and Hi-C sequencing technologies to assemble a reference genome of Z. tau at the chromosomal level. The assembled genome was 421.79 Mb and consisted of six chromosomes (one X-chromosome + five autosomes). The contig N50 was 4.23 Mb. We identified 20,922 protein-coding genes, of which 17,251 (82.45%) were functionally annotated. Additionally, we found 247 rRNAs, 435 tRNAs, 67 small nuclear RNAs, and 829 small RNAs in the genome. Repetitive elements accounted for 55.30 Mb (13.15%) of the genome. This high-quality genome assembly is valuable for evolutionary and genetic studies of Z. tau and its relative species.

Interspecific and intraspecific variation in susceptibility of two co-occurring pest thrips, Frankliniella occidentalis and Thrips palmi, to nine insecticides.

BACKGROUND: Field control of pest thrips mainly relies on insecticides, but the toxicity of insecticides can vary among thrips species and populations. In this study, we examined the susceptibility of multiple field populations of two thrips pests, Frankliniella occidentalis, and Thrips palmi, that often co-occur on vegetables, to nine insecticides belonging to seven subgroups. RESULTS: The highest level of variation in susceptibility among F. occidentalis populations was for spinetoram (73.92 fold difference between most resistant and most susceptible population), followed by three neonicotinoids (8.06-15.99 fold), while among T. palmi populations, it was also for spinetoram (257.19 fold), followed by emamectin benzoate, sulfoxaflor, and acetamiprid (23.64-45.50 fold). These findings suggest evolved resistance to these insecticides in some populations of the two thrips. One population of F. occidentalis had a particularly high level of resistance overall, being the most resistant for five of the nine insecticides tested. Likewise, a population of T. palmi had high resistance to all nine insecticides, again suggesting the evolution of resistance to multiple chemicals. For F. occidentalis, the LC95 values of most populations were higher than the field-recommended dosage for all insecticides except chlorfenapyr and emamectin benzoate. For several T. palmi populations, the LC95 values also tended to be higher than recommended dosages, except in the case of emamectin benzoate and spinetoram. CONCLUSIONS: Our study found interspecific and intraspecific variations in the susceptibility of two thrips to nine insecticides and multiple resistance in some populations, highlighting the need for ongoing monitoring and resistance management. © 2023 Society of Chemical Industry.