Genetic structure of Bemisia tabaci Med populations from home-range countries, inferred by nuclear and cytoplasmic markers: impact on the distribution of the insecticide resistance genes

Evaluation of the evolutionary genetics and population structure of Culex pipiens pallens in Shandong province, China based on knockdown resistance (kdr) mutations and the mtDNA-COI gene

BACKGROUND

Mosquitoes are important vectors for a range of diseases, contributing to high rates of morbidity and mortality in the human population. Culex pipiens pallens is dominant species of Culex mosquito in northern China and a major vector for both West Nile virus and Bancroftian filariasis. Insecticide application were largely applied to control the mosquito-mediated spread of these diseases, contributing to increasing rates of resistance in the mosquito population. The voltage-gated sodium channel (Vgsc) gene is the target site of pyrethroids, and mutations in this gene cause knockdown resistance (kdr). While these kdr mutations are known to be critical to pyrethroid resistance, their evolutionary origins remain poorly understood. Clarifying the origins of these mutations is potential to guide further vector control and disease prevention efforts. Accordingly, the present study was designed to study the evolutionary genetics of kdr mutations and their association with the population structure of Cx. p. pallens in Shandong province, China.

METHODS

Adult Culex females were collected from Shandong province and subjected to morphological identification under a dissection microscope. Genomic DNA were extracted from the collected mosquitoes, the Vgsc gene were amplified via PCR and sequenced to assess kdr allele frequencies, intron polymorphisms, and kdr codon evolution. In addition, population genetic diversity and related population characteristics were assessed by amplifying and sequencing the mitochondrial cytochrome C oxidase I (COI) gene.

RESULTS

Totally, 263 Cx. p. pallens specimens were used for DNA barcoding and sequencing analyses to assess kdr allele frequencies in nine Culex populations. The kdr codon L1014 in the Vgsc gene identified two non-synonymous mutations (L1014F and L1014S) in the analyzed population. These mutations were present in the eastern hilly area and west plain region of Shandong Province. However, only L1014F mutation was detected in the southern mountainous area and Dongying city of Shandong Province, where the mutation frequency was low. Compared to other cities, population in Qingdao revealed significant genetic differentiation. Spatial kdr mutation patterns are likely attributable to some combination of prolonged insecticide-mediated selection coupled with the genetic isolation of these mosquito populations.

CONCLUSIONS

These data suggest that multiple kdr alleles associated with insecticide resistance are present within the Cx. p. pallens populations of Shandong Province, China. The geographical distributions of kdr mutations in this province are likely that the result of prolonged and extensive insecticide application in agricultural contexts together with frequent mosquito population migrations. In contrast, the low-frequency kdr mutation detected in central Shandong Province populations may originate from the limited selection pressure in this area and the relative genetic isolation. Overall, the study compares the genetic patterns revealed by a functional gene with a neutral marker and demonstrates the combined impact of demographic and selection factors on population structure.

Predicting the insecticide-driven mutations in a crop pest insect: Evidence for multiple polymorphisms of acetylcholinesterase gene with potential relevance for resistance to chemicals

The silverleaf whitefly Bemisia tabaci (Gennadius, 1889) (Homoptera: Aleyrodidae) is a serious invasive herbivorous insect pest worldwide. The excessive use of pesticides has progressively selected B. tabaci specimens, reducing the effectiveness of the treatments, and ultimately ending in the selection of pesticide-resistant strains. The management of this crop pest has thus become challenging owing to the level of resistance to all major classes of recommended insecticides. Here, we used in silico techniques for detecting sequence polymorphisms in ace1 gene from naturally occurring B. tabaci variants, and monitor the presence and frequency of the detected putative mutations from 30 populations of the silverleaf whitefly from Egypt and Pakistan. We found several point mutations in ace1-type acetylcholinesterase (ace1) in the studied B. tabaci variants naturally occurring in the field. By comparing ace1 sequence data from an organophosphate-susceptible and an organophosphate-resistant strains of B. tabaci to ace1 sequence data retrieved from GenBank for that species and to nucleotide polymorphisms from other arthropods, we identified novel mutations that could potentially influence insecticide resistance. Homology modeling and molecular docking analyses were performed to determine if the mutation-induced changes in form 1 acetylcholinesterase (AChE1) structure could confer resistance to carbamate and organophosphate insecticides. Mutations had small effects on binding energy (ΔG_b) interactions between mutant AChE1 and insecticides; they altered the conformation of the peripheral anionic site of AChE1, and modified the enzyme surface, and these changes have potential effects on the target-site sensitivity. Altogether, the results from this study provide information on genic variants of B. tabaci ace1 for future monitoring insecticide resistance development and report a potential case of environmentally driven gene variations.

Next-generation molecular diagnostics (TaqMan qPCR and ddPCR) for monitoring insecticide resistance in Bemisia tabaci

BACKGROUND

Insecticide resistance has developed in several populations of the whitefly Bemisia tabaci worldwide and threatens to compromise the efficacy of chemical control. The molecular mechanisms underpinning resistance have been characterized and markers associated with the trait have been identified, allowing the development of diagnostics for individual insects.

RESULTS

TaqMan and Droplet Digital PCR (ddPCR) assays were developed and validated, in individual and pooled whitefly samples, respectively, for the following target-site mutations: the acetylcholinesterase (ace1) F331W mutation conferring organophosphate-resistance; the voltage-gated sodium channel (vgsc) mutations L925I and T929V conferring pyrethroid-resistance; and the acetyl-CoA carboxylase (acc) A2083V mutation conferring ketoenol-resistance. The ddPCR's limit of detection (LoD) was <0.2% (i.e. detection of one heterozygote whitefly in a pool of 249 wild-type individuals). The assays were applied in 11 B. tabaci field populations from four locations in Crete, Greece. The F331W mutation was detected to be fixed or close to fixation in eight of 11 B. tabaci populations, and at lower frequency in the remaining ones. The pyrethroid-resistance mutations were detected at very high frequencies. The A2083V spiromesifen resistance mutation was detected in eight of 11 populations (frequencies = 6.16-89.56%). Spiromesifen phenotypic resistance monitoring showed that the populations tested had variable levels of resistance, ranging from full susceptibility to high resistance. A strong spiromesifen-resistance phenotype-genotype (A2083V) correlation (r_s  = -0.839, P = 0.002) was observed confirming the ddPCR diagnostic value.

CONCLUSION

The ddPCR diagnostics developed in this study are a valuable tool to support evidence-based rational use of insecticides and resistance management strategies. © 2022 Society of Chemical Industry.

Distribution of invasive versus native whitefly species and their pyrethroid knock-down resistance allele in a context of interspecific hybridization

The invasion success of a species in an agrosystem is greatly influenced by environmental factors such as the use of insecticides, by the intrinsic evolutionary capabilities of the species, and also by interactions with resident species. On the island of La Réunion, the successive invasions of MEAM1 and MED whitefly species over the last 20 years have not only led an increased use of insecticides, but have also challenged the resident IO species. To trace the evolution of the 3 species, and the distribution of the kdr mutation (resistance to pyrethroid) in the para-type voltage-gated sodium channel, we genotyped 41 populations (using neutral nuclear markers) and look at the prevalence of the kdr allele. MEAM1 was predominantly present in agrosystems showing quasi fixation of the resistant kdr allele whereas IO was mainly in natural environments and did not have any resistant allele. Hybridization between the two former species was detected in low frequency but has not led to introgression of resistant alleles in the resident species so far. MED showed a limited distribution in agrosystems but all individuals displayed a resistant allele. These highly contrasting patterns of distribution and resistant mutations between invasive and resident whitefly species are further discussed.

Population Structure and Genetic Diversity of the Pepper Weevil (Coleoptera: Curculionidae) Using the COI Barcoding Region

The pepper weevil Anthonomus eugenii Cano (Coleoptera: Curculionidae) is a pest of economic importance for Capsicum species pepper in North America that attacks the reproductive structures of the plant. The insect is distributed across Mexico, the United States, and the Caribbean, and is occasionally found during the pepper growing season in southern Ontario, Canada. Continuous spread of the insect to new areas is partially the result of global pepper trade. Here, we describe the genetic diversity of the pepper weevil using the mitochondrial COI barcoding region across most of its geographic range. In this study, 44 (H1-H44) highly similar haplotypes were identified, the greatest number of haplotypes and haplotype diversity were observed among specimens from its native Mexico, followed by specimens from the United States. Unlike Mexico, a low haplotype diversity was found among specimens from Canada, the Dominican Republic, Italy, and the Netherlands. Out of these 44 haplotypes, 29 are reported for the first time. Haplotype diversity in the Canadian population suggests either multiple and continuous introductions of the pepper weevil into this area or a single introduction of genetically diverse individuals. We discuss the importance of such population genetic data in tailoring pepper weevil management programs, using Canada as an example.

The Costs and Benefits of Two Secondary Symbionts in a Whitefly Host Shape Their Differential Prevalence in the Field

Insects commonly harbor maternally inherited intracellular symbionts in nature, and the microbial partners often exert influence on host reproduction and fitness to promote their prevalence. Here, we investigated composition of symbionts and their biological effects in the invasive Bemisia tabaci MED species of a whitefly complex. Our field surveys revealed that populations of the MED whitefly, in addition to the primary symbiont Portiera, mainly contain two secondary symbionts Hamiltonella, which is nearly fixed in the host populations, and Cardinium with infection frequencies ranging from 0 to 86%. We isolated and established Cardinium-positive and Cardinium-free whitefly lines with a similar nuclear genetic background from a field population, and compared performance of the two whitefly lines. The infection of Cardinium incurred significant fitness costs on the MED whitefly, including reduction of fecundity and egg viability as well as delay in development. We then selectively removed Hamiltonella from the Cardinium-free whitefly line and compared performance of two whitefly lines, one harboring both Portiera and Hamiltonella and the other harboring only Portiera. While depletion of Hamiltonella had little or only marginal effects on the fecundity, developmental rate, and offspring survival, the Hamiltonella-free whitefly line produced very few female offspring, often reducing the progeny female ratio from about 50% to less than 1%. Our findings indicate that the varying costs and benefits of the association between these two symbionts and the MED whitefly may play an important role in shaping their differential prevalence in the field.

Insecticide resistance trait may contribute to genetic cluster change in Bemisia tabaci MED (Hemiptera: Aleyrodidae) as a potential driving force

BACKGROUND

Previously, we reported that the majority of the Bemisia tabaci Mediterranean (MED) populations converged from two dominant genetic clusters (cluster 1 and 2) to one (cluster 2) during 1 year in greenhouse tomatoes in Korea. To find possible mechanisms for this phenomenon, we investigated the concurrent changes in resistance traits of the two clusters for three insecticide classes (organophosphate, pyrethroid, and neonicotinoid).

RESULTS

Since the resistance mutation frequencies in regional samples were either high (i.e. the voltage-sensitive sodium channel L925I/T929V mutations and the F392 acetylcholinesterase 1 mutation) or zero (the nicotinic acetylcholine receptor R81T mutation), no meaningful correlation between the resistance allele frequency and genetic cluster was deduced. However, the actual resistance levels to all three insecticide classes were significantly higher in cluster 2 than in cluster 1, suggesting that cluster 2 has a higher resistance potential. Furthermore, thiamethoxam treatment to the mixed population of clusters 1 and 2 over three generations exhibited a strong tendency of population change from cluster 1 to cluster 2.

CONCLUSION

Our results demonstrated that the insecticide resistance trait is one of the driving forces for rapid genetic cluster change in B. tabaci MED populations. © 2021 Society of Chemical Industry.

Distribution and Genetic Variability of Bemisia tabaci Cryptic Species (Hemiptera: Aleyrodidae) in Italy

Bemisia tabaci is a key pest of horticultural, fibre and ornamental crops worldwide, primarily as a vector of plant viruses. In Italy, B. tabaci has established since the 1980s-1990s in southern regions as well as in Sicily and Sardinia. Recent reports of infestations in some areas of central Italy prompted a new survey to assess the whitefly distribution in the country as well as to update the species and haplotype composition of the populations present in southern Italy and in the main islands. The survey confirmed that B. tabaci is nowadays established in central Italy even at more northern latitudes than those noticed before. Most of the specimens collected throughout the country belonged to the Mediterranean (MED) species. The MEDQ1 and Q2 haplogroups were prevailing in open-field and greenhouse cultivations, respectively, except in Sardinia where only Q1 specimens were found on a wide range of crops and weeds. Population genetics analyses showed that several MEDQ1 haplotypes currently occur in Italy and their distribution is unrelated to evident temporal and geographic trends, except for a new genetic variant which seems to have originated in Sardinia. The MED species is known to better adapt to insecticide treatments and high temperatures, and its northward spread in Italy may have been favoured by the intensive agricultural practices and steady increase in both winter and summer temperatures occurring in the last few decades. The extensive presence of B. tabaci in Italy proves that a strict surveillance for possible new outbreaks of whitefly-transmitted viruses should be addressed to a range of sites that are expanding northwards.

What Is the Spatial Extent of a Bemisia tabaci Population?

Effective pest management depends on basic knowledge about insect dispersal patterns and gene flow in agroecosystems. The globally invasive sweet potato whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is considered a weak flier whose life history nonetheless predisposes it to frequent dispersal, but the scale over which populations exchange migrants, and should therefore be managed, is uncertain. In this review, we synthesize the emergent literature on B. tabaci population genetics to address the question: What spatial scales define B. tabaci populations? We find that within-species genetic differentiation among sites is often low, and evidence of population structuring by host plant or geography is rare. Heterozygote deficits prevail among populations, indicating that migrants from divergent populations are frequently sampled together. Overall, these results suggest that there is high ongoing gene flow over large spatial extents. However, genetic homogeneity typical of recently invading populations could obscure power to detect real isolation among populations. Genome-wide data collected systematically across space and time could distinguish signatures of invasion history from those of ongoing gene flow. Characterizing the spatial extent of B. tabaci populations could reveal whether insecticide rotations can be tailored to specific commodities or if coordination across linked commodities and regions is justified.

Origin of resistance to pyrethroids in the redlegged earth mite (Halotydeus destructor) in Australia: repeated local evolution and migration

BACKGROUND

Halotydeus destructor is a major pest of crops and pastures across southern parts of Australia. This invasive mite has been chemically controlled for over 50 years, but resistance to synthetic pyrethroids and organophosphates is developing. Understanding processes behind the emerging resistance is important for effective management efforts. We undertook a ddRAD pool-sequencing approach to analyse genome-wide single nucleotide polymorphism variation in H. destructor population samples at two scales: local resistance across a set of fields, and regional resistance across their Australian range, along with toxicology bioassays to screen for pyrethroid resistance.

RESULTS

Spatial patterns of genomic variation and resistance at a local scale indicated that genetic similarity among samples was more closely correlated with distance along roads and fence-lines than with straight-line geographic distance. This pattern was particularly strong in resistant samples, which were also more related than susceptible samples, suggesting local spread of resistance within an area after it emerged. By contrast, regional data suggest resistance has emerged repeatedly within parts of Australia. Our de novo annotation of the H. destructor draft genome sequence and Bayesian analysis identified several candidate loci strongly associated with population-level resistance to pyrethroids, located in genomic regions that code for transmembrane transport and signalling proteins that have previously been linked to insecticide resistance in other arthropods.

CONCLUSION

Our findings highlight multiple independent evolutionary events leading to resistance in H. destructor, and demonstrate the utility and cost-effectiveness of a cross-population, genome-wide association study to reveal processes underlying adaptive evolution in a non-model invasive species. © 2019 Society of Chemical Industry.

What has changed in the outbreaking populations of the severe crop pest whitefly species in cassava in two decades?

High populations of African cassava whitefly (Bemisia tabaci) have been associated with epidemics of two viral diseases in Eastern Africa. We investigated population dynamics and genetic patterns by comparing whiteflies collected on cassava in 1997, during the first whitefly upsurges in Uganda, with collections made in 2017 from the same locations. Nuclear markers and mtCOI barcoding sequences were used on 662 samples. The composition of the SSA1 population changed significantly over the 20-year period with the SSA1-SG2 percentage increasing from 0.9 to 48.6%. SSA1-SG1 and SSA1-SG2 clearly interbreed, confirming that they are a single biological species called SSA1. The whitefly species composition changed: in 1997, SSA1, SSA2 and B. afer were present; in 2017, no SSA2 was found. These data and those of other publications do not support the 'invader' hypothesis. Our evidence shows that no new species or new population were found in 20 years, instead, the distribution of already present genetic clusters composing SSA1 species have changed over time and that this may be in response to several factors including the introduction of new cassava varieties or climate changes. The practical implications are that cassava genotypes possessing both whitefly and disease resistances are needed urgently.

Population genetic structure of Bemisia tabaci MED (Hemiptera: Aleyrodidae) in Korea

The sweet potato whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is a major agricultural pest that causes economic damages worldwide. In particular, B. tabaci MED (Mediterranean) has resulted in serious economic losses in tomato production of Korea. In this study, 1,145 B. tabaci MED females from 35 tomato greenhouses in different geographic regions were collected from 2016 to 2018 (17 populations in 2016, 13 in 2017, and five in 2018) and analyzed to investigate their population genetic structures using eight microsatellite markers. The average number of alleles per population (NA) ranged from 2.000 to 5.875, the expected heterozygosity (HE) ranged from 0.218 to 0.600, the observed heterozygosity (HO) ranged from 0.061 to 0.580, and the fixation index inbreeding coefficient (FIS) ranged from -0.391 to 0.872 over the three years of the study. Some significant correlation (p < 0.05) was present between genetic differentiations (FST) and geographical distance, and a comparatively high proportion of variation was found among the B. tabaci MED populations. The B. tabaci MED populations were divided into two well-differentiated genetic clusters within different geographic regions. Interestingly, its genetic structures converged into one genetic cluster during just one year. The reasons for this genetic change were speculated to arise from different fitness, insecticide resistance, and insect movement by human activities.

Comparative transcriptome analysis reveals genetic diversity in the endosymbiont Hamiltonella between native and exotic populations of Bemisia tabaci from Brazil

The whitefly, Bemisia tabaci, is a species complex of more than 40 cryptic species and a major agricultural pest. It causes extensive damage to plants mainly by transmitting plant viruses. There is still a lack of genomic data available for the different whitefly species found in Brazil and their bacterial endosymbionts. Understanding the genetic and transcriptomic composition of these insect pests, the viruses they transmit and the microbiota is crucial to sustainable solutions for farmers to control whiteflies. Illumina RNA-Seq was used to obtain the transcriptome of individual whiteflies from 10 different populations from Brazil including Middle East-Asia Minor 1 (MEAM1), Mediterranean (MED) and New World 2 (NW2). Raw reads were assembled using CLC Genomics Workbench and subsequently mapped to reference genomes. We obtained whitefly complete mitochondrial genomes and draft genomes from the facultative bacterial endosymbiont Hamiltonella for further phylogenetic analyses. In addition, nucleotide sequences of the GroEL chaperonin gene from Hamiltonella from different populations were obtained and analysed. There was concordance in the species clustering using the whitefly complete mitogenome and the mtCOI gene tree. On the other hand, the phylogenetic analysis using the 12 ORF's of Hamiltonella clustered the native species NW2 apart from the exotics MEAM1 and MED. In addition, the amino acid analysis of GroEL chaperonin revealed a deletion only in Hamiltonella infecting NW2 among whiteflies populations analysed which was further confirmed by PCR and Sanger sequencing. The genomic data obtained in this study will aid understanding the functions that Hamiltonella may have in whitefly biology and serve as a reference for further studies regarding whiteflies in Brazil.

An integrative approach to discovering cryptic species within the Bemisia tabaci whitefly species complex

Bemisia tabaci is a cryptic whitefly-species complex that includes some of the most damaging pests and plant-virus vectors of a diverse range of food and fibre crops worldwide. We combine experimental evidence of: (i) differences in reproductive compatibility, (ii) hybrid verification using a specific nuclear DNA marker and hybrid fertility confirmation and (iii) high-throughput sequencing-derived mitogenomes, to show that the "Mediterranean" (MED) B. tabaci comprises at least two distinct biological species; the globally invasive MED from the Mediterranean Basin and the "African silver-leafing" (ASL) from sub-Saharan Africa, which has no associated invasion records. We demonstrate that, contrary to its common name, the "ASL" does not induce squash silver-leafing symptoms and show that species delimitation based on the widely applied 3.5% partial mtCOI gene sequence divergence threshold produces discordant results, depending on the mtCOI region selected. Of the 292 published mtCOI sequences from MED/ASL groups, 158 (54%) are low quality and/or potential pseudogenes. We demonstrate fundamental deficiencies in delimiting cryptic B. tabaci species, based solely on partial sequences of a mitochondrial barcoding gene. We advocate an integrative approach to reveal the true species richness within cryptic species complexes, which is integral to the deployment of effective pest and disease management strategies.

Diversity and Phylogenetic Analyses of Bacterial Symbionts in Three Whitefly Species from Southeast Europe

Bemisia tabaci (Gennadius), Trialeurodes vaporariorum (Westwood), and Siphoninus phillyreae (Haliday) are whitefly species that harm agricultural crops in many regions of the world. These insects live in close association with bacterial symbionts that affect host fitness and adaptation to the environment. In the current study, we surveyed the infection of whitefly populations in Southeast Europe by various bacterial symbionts and performed phylogenetic analyses on the different symbionts detected. Arsenophonus and Hamiltonella were the most prevalent symbionts in all three whitefly species. Rickettsia was found to infect mainly B. tabaci, while Wolbachia mainly infected both B. tabaci and S. phillyreae. Furthermore, Cardinium was rarely found in the investigated whitefly populations, while Fritschea was never found in any of the whitefly species tested. Phylogenetic analyses revealed a diversity of several symbionts (e.g., Hamiltonella, Arsenophonus, Rickettsia), which appeared in several clades. Reproductively isolated B. tabaci and T. vaporariorum shared the same (or highly similar) Hamiltonella and Arsenophonus, while these symbionts were distinctive in S. phillyreae. Interestingly, Arsenophonus from S. phillyreae did not cluster with any of the reported sequences, which could indicate the presence of Arsenophonus, not previously associated with whiteflies. In this study, symbionts (Wolbachia, Rickettsia, and Cardinium) known to infect a wide range of insects each clustered in the same clades independently of the whitefly species. These results indicate horizontal transmission of bacterial symbionts between reproductively isolated whitefly species, a mechanism that can establish new infections that did not previously exist in whiteflies.

A landscape genetic analysis of important agricultural pest species in Tunisia: The whitefly Bemisia tabaci

Combining landscape ecology and genetics provides an excellent framework to appreciate pest population dynamics and dispersal. The genetic architectures of many species are always shaped by environmental constraints. Because little is known about the ecological and genetic traits of Tunisian whitefly populations, the main objective of this work is to highlight patterns of biodiversity, genetic structure and migration routes of this pest. We used nuclear microsatellite loci to analyze B. tabaci populations collected from various agricultural areas across the country and we determine their biotype status. Molecular data were subsequently interpreted in an ecological context supplied from a species distribution model to infer habitat suitability and hereafter the potential connection paths between sampling localities. An analysis of landscape resistance to B. tabaci genetic flow was thus applied to take into account habitat suitability, genetic relatedness and functional connectivity of habitats within a varied landscape matrix. We shed light on the occurrence of three geographically delineated genetic groups with high levels of genetic differentiation within each of them. Potential migration corridors of this pest were then established providing significant advances toward the understanding of genetic features and the dynamic dispersal of this pest. This study supports the hypothesis of a long-distance dispersal of B. tabaci followed by infrequent long-term isolations. The Inference of population sources and colonization routes is critical for the design and implementation of accurate management strategies against this pest.

Population genetic structure of the biological control agent Macrolophus pygmaeus in Mediterranean agroecosystems

Biological control of agricultural pests relies on knowledge of agroecosystem functionality, particularly when affected by the use of mass-produced biological agents. Incorporating pre- and/or post-release information such as genetic diversity and structure on these agents using molecular-based approaches could advance our knowledge of how they perform in agroecosystems. We evaluated the population genetics of Macrolophus pygmaeus, the most widely used predatory mirid against many arthropod pests of greenhouse crops in the Mediterranean region, using the mitochondrial Cytb sequence and microsatellite data, and population genetics and phylogeny approaches. We investigated commercially mass-produced insects (i.e., commercial insects either mass-reared in the laboratory for many generations, or purchased by farmers and released in the greenhouses) and "wild" insects (i.e., that occur naturally outside or are collected in nature for release in the greenhouses). The mirids were mainly collected in agroecosystems in which solanaceous plants are grown in northern Spain, southern France and Greece. Both molecular markers and approaches distinguished 2 genetically differentiated populations. The less genetically diverse population, hereafter named the "commercial" strain included all individuals from laboratory mass-rearings and most releases of commercially bred individuals. The most genetically diverse population mainly comprised individuals originating from noncultivated environments, or from releases of "wild" individuals. Rare examples of hybridization between M. pygmaeus from the 2 populations were observed and asymmetric gene flow was revealed. These findings provide new insights into what happens to M. pygmaeus released in the agroecosystems we studied, and show that it is possible to monitor some commercial strains.

Analysis of Species, Subgroups, and Endosymbionts of Bemisia tabaci (Hemiptera: Aleyrodidae) From Southwestern Cotton Fields in Turkey

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is one of the most important insect pests worldwide including Turkey. Although there are substantial data regarding species composition of Turkish B. tabaci populations, the situation is still not clear and further investigations are needed. Therefore, in this study, species and subgroups of B. tabaci collected from cotton fields in southwestern part of Turkey (Antalya, Aydın, Denizli, and Muğla) were determined using microsatellite analysis, AluI-based mtCOI polymerase chain reaction-random length polymorphism, and sequencing. Secondary endosymbionts were also determined using diagnostic species-specific PCR. Middle East Asia Minor 1 (MEAM1), Mediterranean (MED) Q1, and MED Q2 were the species and subgroups found in this study. The MED species (85.3%) were found to be more dominant than MEAM1. Species status of B. tabaci varied depending on the location. Although all samples collected from Aydın were found to be Q1, three species and subgroups were found in Muğla. Secondary endosymbionts varied according to species and subgroups. Arsenophonus was found only from Q2, while Hamiltonella was detected in MEAM1 and Q1. In addition, high Rickettsia and low Wolbachia infections were detected in MEAM1 and Q1 populations, respectively. In conclusion, for the first time, we report the presence and symbiotic communities of Q1 from Turkey. We also found that the symbiont complement of the Q1 is more congruent with Q1 from Greece than other regions of the world, which may have some interesting implications for movement of this invasive subgroup.

New putative cryptic species detection and genetic network analysis of Bemisia tabaci (Hempitera: Aleyrodidae) in China based on mitochondrial COI sequences

The whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is a cryptic species complex and widely distributed throughout tropical and subtropical regions. To understand the B. tabaci cryptic species diversity in China more comprehensively, in the year 2014 and 2016, a large-scale sampling was conducted from the famous biodiversity hotspot of China, Yunnan province. Mitochondrial cytochrome oxidase I gene sequences were used to identify new putative cryptic species. Phylogenetic analyses were performed using Bayesian methods to evaluate the position of new cryptic species in the context of the B. tabaci diversity in Asia. Two new cryptic species, China 5 and Asia V were identified. In total, 19 B. tabaci cryptic species are present in China, two invasive (MED and MEAM1) and 17 indigenous. A new sibling species of B. tabaci was first defined and reported. Based on the mtCOI sequences and haplotype network analyses, the genetic diversity of MED was far higher than MEAM1. We confirmed the exotic MED was originated from the western Mediterranean regions and first invaded into Yunnan, China. The genetic structures of other four indigenous species (Asia I, Asia II 1, Asia II 6, and China 1) with relatively wide distribution ranges in China were also discussed.

Pyrethroid resistance is associated with a kdr-type mutation (L1014F) in the potato tuber moth Tecia solanivora

BACKGROUND

The Guatemalan potato tuber moth, Tecia solanivora, has been the most important pest species in Hispanico-American potato fields since its first record on potatoes in 1956 in Guatemala. This insect pest has been spreading to other parts of the world, including the Canary Islands in Europe. Tuber moth control relies heavily on the use of insecticides, including pyrethroids. Here, we assessed the likelihood of control failures and performed concentration-response bioassays in five Colombian strains of T. solanivora to evaluate their susceptibilities to the pyrethroid permethrin.

RESULTS

Evidence of control failures was observed in four strains tested, which exhibited moderate resistance levels (i.e. ranging from 5.4- to 24.4-fold). However, no spatial dependence was observed between the permethrin LC₅₀ values and the geographic distances among the tuber moth strains. In order to evaluate whether permethrin resistance was mediated by potential mutations in the para-type sodium channels of T. solanivora, the IIS4-IIS6 region of the para gene was PCR amplified and sequenced from the five strains tested. As demonstrated across a range of different arthropod species that exhibited knockdown resistance (kdr), we observed a single point substitution (L1014F) at high frequencies in the para gene of all four resistant strains.

CONCLUSION

This is the first identification of a target-site-alteration-based resistance in the Guatemalan potato tuber moth T. solanivora, which is widespread and exhibits high frequencies among geographically distant strains, indicating that pyrethroids are probably becoming ineffective for the control of this pest species. © 2016 Society of Chemical Industry.

Global Population Structure of a Worldwide Pest and Virus Vector: Genetic Diversity and Population History of the Bemisia tabaci Sibling Species Group

The whitefly Bemisia tabaci sibling species (sibsp.) group comprises morphologically indiscernible lineages of well-known exemplars referred to as biotypes. It is distributed throughout tropical and subtropical latitudes and includes the contemporary invasive haplotypes, termed B and Q. Several well-studied B. tabaci biotypes exhibit ecological and biological diversity, however, most members are poorly studied or completely uncharacterized. Genetic studies have revealed substantial diversity within the group based on a fragment of the mitochondrial cytochrome oxidase I (mtCOI) sequence (haplotypes), with other tested markers being less useful for deep phylogenetic comparisons. The view of global relationships within the B. tabaci sibsp. group is largely derived from this single marker, making assessment of gene flow and genetic structure difficult at the population level. Here, the population structure was explored for B. tabaci in a global context using nuclear data from variable microsatellite markers. Worldwide collections were examined representing most of the available diversity, including known monophagous, polyphagous, invasive, and indigenous haplotypes. Well-characterized biotypes and other related geographic lineages discovered represented highly differentiated genetic clusters with little or no evidence of gene flow. The invasive B and Q biotypes exhibited moderate to high levels of genetic diversity, suggesting that they stemmed from large founding populations that have maintained ancestral variation, despite homogenizing effects, possibly due to human-mediated among-population gene flow. Results of the microsatellite analyses are in general agreement with published mtCOI phylogenies; however, notable conflicts exist between the nuclear and mitochondrial relationships, highlighting the need for a multifaceted approach to delineate the evolutionary history of the group. This study supports the hypothesis that the extant B. tabaci sibsp. group contains ancient genetic entities and highlights the vast cryptic diversity throughout the genome in the group.

Landscape genetic structure and evolutionary genetics of insecticide resistance gene mutations in Anopheles sinensis

Background

Anopheles sinensis is one of the most abundant vectors of malaria and other diseases in Asia. Vector control through the use of insecticides is the front line control method of vector-borne diseases. Pyrethroids are the most commonly used insecticides due to their low toxicity to vertebrates and low repellency. However, the extensive use of insecticides has imposed strong selection pressure on mosquito populations for resistance. High levels of resistance to pyrethroid insecticides and various mutations and haplotypes in the para sodium channel gene that confers knockdown resistance (kdr) have been detected in An. sinensis. Despite the importance of kdr mutations in pyrethroid resistance, the evolutionary origin of the kdr mutations is unknown. This study aims to examine the evolutionary genetics of kdr mutations in relation to spatial population genetic structure of An. sinensis.

Methods

Adults or larvae of Anopheles sinensis were collected from various geographic locations in China. DNA was extracted from individual mosquitoes. PCR amplification and DNA sequencing of the para-type sodium channel gene were conducted to analyse kdr allele frequency distribution, kdr codon upstream and downstream intron polymorphism, population genetic diversity and kdr codon evolution. The mitochondrial cytochrome c oxidase COI and COII genes were amplified and sequenced to examine population variations, genetic differentiation, spatial population structure, population expansion and gene flow patterns.

Results

Three non-synonymous mutations (L1014F, L1014C, and L1014S) were detected at the kdr codon L1014 of para-type sodium channel gene. A patchy distribution of kdr mutation allele frequencies from southern to central China was found. Near fixation of kdr mutation was detected in populations from central China, but no kdr mutations were found in populations from southwestern China. More than eight independent mutation events were detected in the three kdr alleles, and at least one of them evolved multiple times subsequent to their first divergence. Based on sequence analysis of the mitochondrial COI and COII genes, significant and large genetic differentiation was detected between populations from southwestern China and central China. The patchy distribution of kdr mutation frequencies is likely a consequence of geographic isolation in the mosquito populations and the long-term insecticide selection.

Conclusion

Our results indicate multiple origins of the kdr insecticide-resistant alleles in An. sinensis from southern and central China. Local selection related to intense and prolonged use of insecticide for agricultural purposes, as well as frequent migrations among populations are likely the explanations for the patchy distribution of kdr mutations in China. On the contrary, the lack of kdr mutations in Yunnan and Sichuan is likely a consequence of genetic isolation and absence of strong selection pressure. The present study compares the genetic patterns revealed by a functional gene with a neutral marker and demonstrates the combined impact of demographic and selection factors on population structure.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1513-6) contains supplementary material, which is available to authorized users.

Further insight into reproductive incompatibility between putative cryptic species of the Bemisia tabaci whitefly complex

The whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae), with its global distribution and extensive genetic diversity, is now known to be a complex of over 35 cryptic species. However, a satisfactory resolution of the systematics of this species complex is yet to be achieved. Here, we designed experiments to examine reproductive compatibility among species with different levels of mitochondrial cytochrome oxidase I (mtCOI) divergence. The data show that putative species with mtCOI divergence of >8% between them consistently exhibited complete reproductive isolation. However, two of the putative species, Asia II 9 and Asia II 3, with mtCOI divergence of 4.47% between them, exhibited near complete reproductive compatibility in one direction of their cross, and partial reproductive compatibility in the other direction. Together with some recent reports on this topic from the literature, our data indicates that, while divergence in the mtCOI sequences provides a valid molecular marker for species delimitation in most clades, more genetic markers and more sophisticated molecular phylogeny will be required to achieve adequate delimitation of all species in this whitefly complex. While many attempts have been made to examine the reproductive compatibility among genetic groups of the B. tabaci complex, our study represents the first effort to conduct crossing experiments with putative species that were chosen with considerations of their genetic divergence. In light of the new data, we discuss the best strategy and protocols to conduct further molecular phylogenetic analysis and crossing trials, in order to reveal the overall pattern of reproductive incompatibility among species of this whitefly complex.

Transcription analysis of neonicotinoid resistance in Mediterranean (MED) populations of B. tabaci reveal novel cytochrome P450s, but no nAChR mutations associated with the phenotype

Background

Bemisia tabaci is one of the most damaging agricultural pests world-wide. Although its control is based on insecticides, B. tabaci has developed resistance against almost all classes of insecticides, including neonicotinoids.

Results

We employed an RNA-seq approach to generate genome wide expression data and identify genes associated with neonicotinoid resistance in Mediterranean (MED) B. tabaci (Q1 biotype). Twelve libraries from insecticide resistant and susceptible whitefly populations were sequenced on an Illumina Next-generation sequencing platform, and genomic sequence information of approximately 73 Gbp was generated.

A reference transcriptome was built by de novo assembly and functionally annotated. A total of 146 P450s, 18 GSTs and 23 CCEs enzymes (unigenes) potentially involved in the detoxification of xenobiotics were identified, along with 78 contigs encoding putative target proteins of six different insecticide classes. Ten unigenes encoding nicotinic Acetylcholine Receptors (nAChR), the target of neoinicotinoids, were identified and phylogenetically classified. No nAChR polymorphism potentially related with the resistant phenotypes, was observed among the studied strains.

DE analysis revealed that among the 550 differentially (logFC > 1) over-transcribed unigenes, 52 detoxification enzymes were over expressed including unigenes with orthologues in P450s, GSTs, CCE and UDP-glucuronosyltransferases.

Eight P450 unigenes belonging to clades CYP2, CYP3 and CYP4 were highly up-regulated (logFC > 2) including CYP6CM1, a gene already known to confer imidacloprid resistance in B. tabaci. Using quantitative qPCRs, a larger screening of field MED B. tabaci from Crete with known neonicotinoid phenotype was performed to associate expression levels of P450s with resistance levels. Expression levels of five P450s, including CYP6CM1, were found associated with neonicotinoid resistance. However, a significant correlation was found only in CYP303 and CYP6CX3, with imidacloprid and acetamiprid respectively.

Conclusion

Our work has generated new toxicological data and genomic resources which will significantly enrich the available dataset and substantially facilitate the molecular studies in MED B. tabaci. No evidence of target site neonicotinoid resistance has been found. Eight P450 unigenes, including CYP6CM1, were found significantly over-expressed in resistant B. tabaci. This study suggests at least two novel P450s (CYP303 and CYP6CX3) as candidates for their functional characterization as detoxification mechanisms of neonicotinoid resistance in B. tabaci.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2161-5) contains supplementary material, which is available to authorized users.

Development of a lateral flow test to detect metabolic resistance in Bemisia tabaci mediated by CYP6CM1, a cytochrome P450 with broad spectrum catalytic efficiency

Cotton whitefly, Bemisia tabaci (Genn.) (Homoptera: Aleyrodidae) is a major sucking pest in many agricultural and horticultural cropping systems globally. The frequent use of insecticides of different mode of action classes resulted in populations resisting treatments used to keep numbers under economic damage thresholds. Recently it was shown that resistance to neonicotinoids such as imidacloprid is linked to the over-expression of CYP6CM1, a cytochrome P450 monooxygenase detoxifying imidacloprid and other neonicotinoid insecticides when recombinantly expressed in insect cells. However over-expression of CYP6CM1 is also known to confer cross-resistance to pymetrozine, an insecticide not belonging to the chemical class of neonicotinoids. In addition we were able to demonstrate by LC-MS/MS analysis the metabolisation of pyriproxyfen by recombinantly expressed CYP6CM1. Based on our results CYP6CM1 is one of the most versatile detoxification enzymes yet identified in a pest of agricultural importance, as it detoxifies a diverse range of chemical classes used to control whiteflies. Therefore we developed a field-diagnostic antibody-based lateral flow assay which detects CYP6CM1 protein at levels providing resistance to neonicotinoids and other insecticides. The ELISA based test kit can be used as a diagnostic tool to support resistance management strategies based on the alternation of different modes of action of insecticides.

Nuclear and cytoplasmic differentiation among Mediterranean populations of Bemisia tabaci: testing the biological relevance of cytotypes

BACKGROUND

The taxonomy of the species complex Bemisia tabaci is still an unresolved issue. Recently, phylogenetic analysis based on mtCOI identified 31 cryptic species. However, mitochondrial diversity is observed within these species, associated with distinct symbiotic bacterial communities forming associations, which here are called cytotypes. The authors investigated the biological significance of two cytotypes (Q1 and Q2) belonging to the Mediterranean species, which have only been found in allopatry in the Western Mediterranean to date. Sampling was done over a few years in Western Europe, and sympatric situations were found that allowed their reproductive compatibility to be tested in the field with the use of microsatellites.

RESULTS

The field survey indicated that, in spite of its recent introduction, Q2 is well established in France and Spain, where it coexists with Q1. Microsatellite data showed that, in allopatry, Q1 and Q2 are highly differentiated, while there is little or no genetic differentiation when they coexist in sympatry, suggesting a high rate of hybridisation. Crossing experiments in the lab confirmed their interfertility.

CONCLUSION

Q1 and Q2 hybridise, which confirms that they belong to the same species, in spite of the high degree of genetic differentiation at both the cytoplasmic and nuclear levels, and also suggests that their symbiotic bacteria do not prevent hybridisation.

Agroecosystems shape population genetic structure of the greenhouse whitefly in Northern and Southern Europe

BACKGROUND

To predict further invasions of pests it is important to understand what factors contribute to the genetic structure of their populations. Cosmopolitan pest species are ideal for studying how different agroecosystems affect population genetic structure within a species at different climatic extremes. We undertook the first population genetic study of the greenhouse whitefly (Trialeurodes vaporariorum), a cosmopolitan invasive herbivore, and examined the genetic structure of this species in Northern and Southern Europe. In Finland, cold temperatures limit whiteflies to greenhouses and prevent them from overwintering in nature, and in Greece, milder temperatures allow whiteflies to inhabit both fields and greenhouses year round, providing a greater potential for connectivity among populations. Using nine microsatellite markers, we genotyped 1274 T. vaporariorum females collected from 18 greenhouses in Finland and eight greenhouses as well as eight fields in Greece.

RESULTS

Populations from Finland were less diverse than those from Greece, suggesting that Greek populations are larger and subjected to fewer bottlenecks. Moreover, there was significant population genetic structure in both countries that was explained by different factors. Habitat (field vs. greenhouse) together with longitude explained genetic structure in Greece, whereas in Finland, genetic structure was explained by host plant species. Furthermore, there was no temporal genetic structure among populations in Finland, suggesting that year-round populations are able to persist in greenhouses.

CONCLUSIONS

Taken together our results show that greenhouse agroecosystems can limit gene flow among populations in both climate zones. Fragmented populations in greenhouses could allow for efficient pest management. However, pest persistence in both climate zones, coupled with increasing opportunities for naturalization in temperate latitudes due to climate change, highlight challenges for the management of cosmopolitan pests in Northern and Southern Europe.

Global distribution and origin of target site insecticide resistance mutations in Tetranychus urticae

The control of Tetranychus urticae, a worldwide agricultural pest, is largely dependent on pesticides. However, their efficacy is often compromised by the development of resistance. Recent molecular studies identified a number of target site resistance mutations, such as G119S, A201S, T280A, G328A, F331W in the acetylcholinesterase gene, L1024V, A1215D, F1538I in the voltage-gated sodium channel gene, G314D and G326E in glutamate-gated chloride channel genes, G126S, I136T, S141F, D161G, P262T in the cytochrome b and the I1017F in the chitin synthase 1 gene. We examined their distribution, by sequencing the relevant gene fragments in a large number of T. urticae collections from a wide geographic range. Our study revealed that most of the resistance mutations are spread worldwide, with remarkably variable frequencies. Furthermore, we analyzed the variability of the ace locus, which has been subjected to longer periods of selection pressure historically, to investigate the evolutionary origin of ace resistant alleles and determine whether they resulted from single or multiple mutation events. By sequencing a 1540 bp ace fragment, encompassing the resistance mutations and downstream introns in 139 T. urticae individuals from 27 countries, we identified 6 susceptible and 31 resistant alleles which have arisen from at least three independent mutation events. The frequency and distribution of these ace haplotypes varied geographically, suggesting an interplay between different mutational events, gene flow and local selection.