Standardized molecular diagnostic tool for the identification of cryptic species within the Bemisia tabaci complex

Bemisia tabaci MEAM1 still remains the dominant species in open field crops in Brazil

Abstract Among Bemisia tabaci species, the invasive MEAM1 and MED species are key agricultural pests for many crops. In Brazil, most part of B. tabaci population outbreaks were associated with MEAM1, which, since 1990s quickly spread across the entire country. Later in 2014, the MED was identified in Brazil, initially more restricted to greenhouses, but suddenly reaching new areas in the South and Southeast open regions. Thus, our objective was to investigate the geographical distribution of MEAM1 and MED on open field crops in Brazil. MEAM1 is still the predominant species on open field crops such as soybean, cotton, and tomato. The sequencing of a cytochrome c oxidase subunit I (COI) gene fragment revealed a single haplotype of MEAM1, suggesting the establishment of a single MEAM1 strain in the country. The haplotypes found for MEAM1 and MED are genetically related to the globally dispersed strains, Jap1 and Mch1, respectively. Continuous monitoring of B. tabaci species is crucial because landscape alterations, climatic changes, and pest management methods may shift the B. tabaci species distribution and dominance in Brazilian crop areas.

Expression profile of CYP402C1 and its role in resistance to imidacloprid in the whitefly, Bemisia tabaci

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is a cosmopolitan insect pest causing serious damage to crop production. Cytochromes P450 (CYPs) of B. tabaci are widely known to be involved in the metabolic resistance to a variety of insecticides, continuously increasing the difficulty in controlling this pest. In this study, four P450 genes (CYP6CM1, CYP6CX1, CYP6CX3, and CYP402C1) in B. tabaci exhibited correlations with the resistance to imidacloprid. We have focused on trying to understand the function and metabolism capacity of CYP402C1. The expression profiles of CYP402C1 were examined by reverse transcription quantitative real-time PCR and fluorescence in situ hybridizations. Its role in resistance to imidacloprid was investigated by RNA interference, transgenic Drosophila melanogaster, and heterologous expression. The results showed that CYP402C1 was highly expressed in the active feeding stages of B. tabaci, such as nymphs and female adults. CYP402C1 was mainly expressed in midguts of nymphs and adults, especially in the filter chamber. Knockdown of CYP402C1 significantly decreased the resistance of B. tabaci to imidacloprid by 3.96-fold (50% lethal concentration: 186.46 versus 47.08 mg/L). Overexpression of CYP402C1 in a transgenic D. melanogaster line (Gal4 > UAS-CYP402C1) significantly increased the resistance to imidacloprid from 12.68- to 14.92-fold (129.01 and 151.80 mg/L versus 1925.14 mg/L). The heterologous expression of CYP402C1 showed a metabolism ability of imidacloprid (imidacloprid decreased by 12.51% within 2 h). This study provides new insights for CYP402C1 function in B. tabaci and will help develop new strategies in B. tabaci control and its insecticide resistance.

Screening for insecticide resistance in Australian field populations of Bemisia tabaci (Hemiptera: Aleyrodidae) using bioassays and DNA sequencing

BACKGROUND

Species within the Bemisia tabaci cryptic species complex can cause significant crop damage. We used high-throughput amplicon sequencing to identify the species composition and resistance allele genotypes in field populations from cotton fields in Australia. For selected populations, the resistance phenotype was determined in bioassays and compared with sequencing data.

RESULTS

A metabarcoding approach was used to analyse the species composition in 144 field populations collected between 2013 and 2021. Two mixed AUS I and MEAM1 populations were detected, whereas the remaining 142 populations consisted of MEAM1 only. High-throughput sequencing of organophosphate and pyrethroid resistance gene amplicons showed that the organophosphate resistance allele F331W was fixed (> 99%) in all MEAM1 populations, whereas the pyrethroid resistance allele L925I in the voltage-gated sodium channel gene was detected at varying frequencies [1.0%-7.0% (43 populations); 27.7% and 42.1% (two populations); 95%-97.5% (three populations)]. Neither organophosphate nor pyrethroid resistance alleles were detected in the AUS I populations. Pyrethroid bioassays of 85 MEAM1 field-derived populations detected no resistance in 51 populations, whereas 32 populations showed low frequency resistance, and 2 populations were highly resistant.

CONCLUSIONS

We demonstrate that high-throughput sequencing and bioassays are complementary approaches. The detection of target site mutations and the phenotypic provides a comprehensive analysis of the low-level resistance to pyrethroids that is present in Australian cotton farms. By contrast, a limited survey of whitefly populations from horticulture found evidence of high-level resistance against pyrethroids. Furthermore, we found that the F331W allele (linked to organophosphate resistance) is ubiquitous in Australian MEAM1. © 2022 Commonwealth of Australia. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda

Bean leaf beetle (BLB) (Ootheca mutabilis) has emerged as an important bean pest in Uganda, leading to devastating crop losses. There is limited information on the population genetic structure of BLB despite its importance. In this study, novel microsatellite DNA markers and the partial mitochondrial cytochrome oxidase subunit I (mtCOI) gene sequences were used to analyze the spatial population genetic structure, genetic differentiation and haplotype diversity of 86 O. mutabilis samples from 16 (districts) populations. We identified 19,356 simple sequence repeats (SSRs) (mono, di-, tri-, tetra-, penta-, and hexa-nucleotides) of which 81 di, tri and tetra-nucleotides were selected for primer synthesis. Five highly polymorphic SSR markers (4-21 alleles, heterozygosity 0.59-0.84, polymorphic information content (PIC) 50.13-83.14%) were used for this study. Analyses of the 16 O. mutabilis populations with these five novel SSRs found nearly all the genetic variation occurring within populations and there was no evidence of genetic differentiation detected for both types of markers. Also, there was no evidence of isolation by distance between geographical and genetic distances for SSR data and mtCOI data except in one agro-ecological zone for mtCOI data. Bayesian clustering identified a signature of admixture that suggests genetic contributions from two hypothetical ancestral genetic lineages for both types of markers, and the minimum-spanning haplotype network showed low differentiation in minor haplotypes from the most common haplotype with the most common haplotype occurring in all the 16 districts. A lack of genetic differentiation indicates unrestricted migrations between populations. This information will contribute to the design of BLB control strategies.

A high-throughput amplicon sequencing approach for population-wide species diversity and composition survey

Management of agricultural pests requires an understanding of pest species diversity, their interactions with beneficial insects and spatial-temporal patterns of pest abundance. Invasive and agriculturally important insect pests can build up very high populations, especially in cropping landscapes. Traditionally, sampling effort for species identification involves small sample sizes and is labour intensive. Here, we describe a multi-primer high throughput sequencing (HTS) metabarcoding method and associated analytical workflow for a rapid, intensive, high-volume survey of pest species compositions. We demonstrate our method using the taxonomically challenging Bemisia pest cryptic species complex as examples. The whiteflies Bemisia including the 'tabaci' species are agriculturally important capable of vectoring diverse plant viruses that cause diseases and crop losses. Our multi-primer metabarcoding HTS amplicon approach simultaneously process high volumes of whitefly individuals, with efficiency to detect rare (i.e., 1%) test-species, while our improved whitefly primers for metabarcoding also detected beneficial hymenopteran parasitoid species from whitefly nymphs. Field-testing our redesigned Bemisia metabarcoding primer sets across the Tanzania, Uganda and Malawi cassava cultivation landscapes, we identified the sub-Saharan Africa 1 Bemisia putative species as the dominant pest species, with other cryptic Bemisia species being detected at various abundances. We also provide evidence that Bemisia species compositions can be affected by host crops and sampling techniques that target either nymphs or adults. Our multi-primer HTS metabarcoding method incorporated two over-lapping amplicons of 472bp and 518bp that spanned the entire 657bp 3' barcoding region for Bemisia, and is particularly suitable to molecular diagnostic surveys of this highly cryptic insect pest species complex that also typically exhibited high population densities in heavy crop infestation episodes. Our approach can be adopted to understand species biodiversity across landscapes, with broad implications for improving trans-boundary biosecurity preparedness, thus contributing to molecular ecological knowledge and the development of control strategies for high-density, cryptic, pest-species complexes.

On species delimitation, hybridization and population structure of cassava whitefly in Africa

The Bemisia cassava whitefly complex includes species that cause severe crop damage through vectoring cassava viruses in eastern Africa. Currently, this whitefly complex is divided into species and subgroups (SG) based on very limited molecular markers that do not allow clear definition of species and population structure. Based on 14,358 genome-wide SNPs from 62 Bemisia cassava whitefly individuals belonging to sub-Saharan African species (SSA1, SSA2 and SSA4), and using a well-curated mtCOI gene database, we show clear incongruities in previous taxonomic approaches underpinned by effects from pseudogenes. We show that the SSA4 species is nested within SSA2, and that populations of the SSA1 species comprise well-defined south-eastern (Madagascar, Tanzania) and north-western (Nigeria, Democratic Republic of Congo, Burundi) putative sub-species. Signatures of allopatric incipient speciation, and the presence of a 'hybrid zone' separating the two putative sub-species were also detected. These findings provide insights into the evolution and molecular ecology of a highly cryptic hemipteran insect complex in African, and allow the systematic use of genomic data to be incorporated in the development of management strategies for this cassava pest.

Spatial distribution and colonization pattern of Bemisia tabaci in tropical tomato crops

BACKGROUND

In precision integrated pest management, management tactics are implemented only where and when needed, by identifying the sites where the pest population has reached economic thresholds. Tomato, Solanum lycopersicum (Linn.), is a vegetable cultivated worldwide, but its production is reduced by insect pests such as the whitefly, Bemisia tabaci (Genn.). To improve management, there is a need to understand B. tabaci spatial dynamics in tomato fields, which will elucidate colonization patterns and may improve management of this pest. Thus, this study was conducted to assess the spatial autocorrelation, distribution, and colonization patterns of B. tabaci in 19 commercial tomato fields through the growing season.

RESULTS

A total of 69 isotropic variograms were fit for B. tabaci. The insect distribution was aggregated with a strong level of spatial dependence. Ranges of spatial dependence varied from 0.53 to 19.05 m and 0.5 to 20 m for adults and nymphs, respectively. Overall, densities of adults and nymphs were higher and reached the economic threshold mainly at the field edges.

CONCLUSION

Our results suggest a colonization pattern for B. tabaci starting at the edges and spreading inwards in to the tomato fields. This study can improve B. tabaci management in tomato fields, especially scouting and decision-making to treat fields. Scouting for this pest should be directed to the field edges, with sample points at least 20 m apart from each other for independent insect counts. © 2020 Society of Chemical Industry.

Evaluating Acylsugars-Mediated Resistance in Tomato against Bemisia tabaci and Transmission of Tomato Yellow Leaf Curl Virus

The sweetpotato whitefly, Bemisia tabaci, is a major pest of cultivated tomato. Whitefly feeding-related injuries and transmission of viruses including tomato yellow leaf curl virus (TYLCV) cause serious losses. Management strategy includes planting resistant cultivars/hybrids. However, TYLCV resistance is incomplete and whiteflies on TYLCV-resistant cultivars/hybrids are managed by insecticides. Acylsugars'-mediated resistance against whiteflies has been introgressed from wild solanums into cultivated tomato. This study evaluated acylsugar-producing tomato lines with quantitative trait loci (QTL) containing introgressions from Solanum pennellii LA716, known to alter acylsugars' levels or chemistry. Evaluated acylsugar-producing lines were the benchmark line CU071026, QTL6/CU071026-a CU071026 sister line with QTL6, and three other CU071026 sister lines with varying QTLs-FA2/CU71026, FA7/CU071026, and FA2/FA7/CU071026. Non-acylsugar tomato hybrid Florida 47 (FL47) was also evaluated. Acylsugars' amounts in FA7/CU071026 and FA2/FA7/CU071026 were 1.4 to 2.2 times greater than in other acylsugar-producing lines. Short chain fatty acid, i-C5, was dominant in all acylsugar-producing lines. Long chain fatty acids, n-C10 and n-C12, were more abundant in FA7/CU071026 and FA2/FA7/CU071026 than in other acylsugar-producing lines. Whiteflies preferentially settled on non-acylsugar hybrid FL47 leaves over three out of five acylsugar-producing lines, and whiteflies settled 5 to 85 times more on abaxial than adaxial leaf surface of FL47 than on acylsugar-producing lines. Whiteflies' survival was 1.5 to 1.9 times lower on acylsugar-producing lines than in FL47. Nevertheless, whiteflies' developmental time was up to 12.5% shorter on acylsugar-producing lines than on FL47. TYLCV infection following whitefly-mediated transmission to acylsugar-producing lines was 1.4 to 2.8 times lower than FL47, and TYLCV acquisition by whiteflies from acylsugar-producing lines was up to 77% lower than from FL47. However, TYLCV accumulation in acylsugar-producing lines following infection and TYLCV loads in whiteflies upon acquisition from acylsugar-producing lines were not different from FL47. Combining TYLCV resistance with acylsugars'-mediated whitefly resistance in cultivated tomato could substantially benefit whiteflies and TYLCV management.

Secondary endosymbiont diversity of Bemisia tabaci and its parasitoids

Cotton whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is one of the most important insect pests worldwide. It is known as a species complex consisting of at least 40 cryptic species. Although there are substantial data regarding species composition, parasitoids and endosymbionts of B. tabaci, data on relationship between the pest, parasitoids and endosymbionts are very restricted. Therefore, in this study, secondary endosymbionts in populations of B. tabaci and their parasitoids collected from Turkey and the USA were determined by PCR-based DNA analysis. Whitefly populations in Turkey represented both Mediterranean (MED) and Middle East-Asia Minor1 (MEAM1) genotypes from single or mixed populations of both genotypes. Arsenophonus, Rickettsia and Wolbachia were found in MED, while Hamiltonella and Rickettsia in MEAM1. Whitefly populations collected from Arizona were all MEAM1 and dually infected with Hamiltonella and Rickettsia. The aphelinid parasitoids Encarsia lutea and Eretmocerus mundus predominated in all Turkish populations. While almost all En. lutea populations were infected with Wolbachia, no endosymbionts were detected in any Er. mundus. Parasitoid species and the pattern of secondary endosymbiont infection in Arizona populations were different with Rickettsia detected only from Encarsia sophia while both Rickettsia and Wolbachia were found in Eretmocerus species. As a result, four secondary endosymbionts, namely, Rickettsia, Hamiltonella, Arsenophonus and Wolbachia, were detected from B.tabaci and its parasitoids. Among them only Wolbachia and Rickettsia were found in both the pest and parasitoids. It is conclude that further studies should be pursued to determine effect of these endosymbionts on biology of the parasitoids and success in biological control of B. tabaci.

Climate model for seasonal variation in Bemisia tabaci using CLIMEX in tomato crops

The whitefly, Bemisia tabaci, is considered one of the most important pests for tomato Solanum lycopersicum. The population density of this pest varies throughout the year in response to seasonal variation. Studies of seasonality are important to understand the ecological dynamics and insect population in crops and help to identify which seasons have the best climatic conditions for the growth and development of this insect species. In this research, we used CLIMEX to estimate the seasonal abundance of a species in relation to climate over time and species geographical distribution. Therefore, this research is designed to infer the mechanisms affecting population processes, rather than simply provide an empirical description of field observations based on matching patterns of meteorological data. In this research, we identified monthly suitability for Bemisia tabaci, with the climate models, for 12 commercial tomato crop locations through CLIMEX (version 4.0). We observed that B. tabaci displays seasonality with increased abundance in tomato crops during March, April, May, June, October and November (first year) and during March, April, May, September and October (second year) in all monitored areas. During this period, our model demonstrated a strong agreement between B. tabaci density and CLIMEX weekly growth index (GIw), which indicates significant reliability of our model results. Our results may be useful to design sampling and control strategies, in periods and locations when there is high suitability for B. tabaci.

Going global - genomic insights into insect invasions

The spread of invasive insect pests is becoming an increasing problem for agriculture globally. We discuss a number of invasive insects, already of major economic significance that have recently expanded their range to become truly global threats. These include the noctuid moths Helicoverpa and Spodoptera, whose caterpillars have long been among the worst pests in their native Old and New World habitats, respectively, and the whitefly Bemisia, a major vector of plant virus diseases. Importantly, genomic resources for these species have recently become available, allowing research to move beyond the restrictions imposed by earlier approaches limited to a single or few mitochondrial and nuclear markers, to employ genome-wide genotyping and resequencing protocols. These studies have shown hybridisation within the various species complexes, identified regions under selection in agricultural environments, and enable monitoring of genes important as biosecurity risks through introgression into established populations free of the genes. In all cases studied, global trade has emerged as the probable cause of insect spread, making it ever more important that biosecurity protocols and agencies work with researchers to make the most effective use of emerging genomic resources and tools.

Incidence of Bemisia tabaci Gennadius (Aleyrodidae, Aleyrodinae) in four Atlantic Forest fragments

ABSTRACT: Bemisia tabaci has been known in Brazil since 1923, in cotton, soybean, tomato and bean crops, and it leads to economic losses due to the transmission of different viruses. Studies focused on analyzing the incidence of this pest in forest environments remain scarce in the literature. Thus, the aims of the current study were to investigate the incidence of B. tabaci in four Atlantic forest fragments in Rio de Janeiro state, Brazil, as well as to record its hosts in these environments. An area of approximately one hectare was delimited for each Atlantic forest fragment based on images provided by the Brazilian Institute of Geography and Statistics (IBGE). Plants presenting whiteflies were tagged throughout monthly inspections conducted in each forest fragment in order to be checked during the following visits. Collection criteria were adopted to assure sampling standardization in each fragment. The mean number of whiteflies was analyzed through the Kruskal-Wallis non-parametric test, which was followed by Dunn’s post-hoc test, both at 5% significance level. The total of 216 collections were conducted and 1,126 B. tabaci fourth-instar nymphs were collected in four Atlantic forest fragments over 24 months. The following hosts were recorded for the first time: Caesalpinia pluviosa Benth, Tradescantia zebrina Linn., Impatiens walleriana Linn., Cupania sp. and Talisia esculenta Klotzsch & O. Berg. The dry period (winter) enabled the increase of B. tabaci populations in the herein investigated forest fragments.

African cassava whitefly, Bemisia tabaci, cassava colonization preferences and control implications

Cassava is a staple food for people across sub-Saharan Africa. Over the last 20 years, there has been an increased frequency of outbreaks and crop damage in this region caused by the cassava-adapted Bemisia tabaci putative species. Little is known about when and why B. tabaci adults move and colonize new cassava crops, especially in farming systems that contain a mixture of cultivar types and plant ages. Here, we assessed experimentally whether the age and variety of cassava affected the density of B. tabaci. We also tested whether the age and variety of the source cassava field affected the variety preference of B. tabaci when they colonized new cassava plants. We placed uninfested potted "sentinel" plants of three cassava varieties (Nam 130, Nase 14, and Njule Red) in source fields containing one of two varieties (Nam 130 or Nase 14) and one of three age classes (young, medium, or old). After two weeks, the numbers of nymphs on the sentinel plants were used as a measure of colonization. Molecular identification revealed that the B. tabaci species was sub-Saharan Africa 1 (SSA1). We found a positive correlation between the density of nymphs on sentinel plants and the density of adults in the source field. The density of nymphs on the sentinels was not significantly related to the age of the source field. Bemisia tabaci adults did not preferentially colonize the sentinel plant of the same variety as the source field. There was a significant interactive effect, however, between the source and sentinel variety that may indicate variability in colonization. We conclude that managing cassava source fields to reduce B. tabaci abundance will be more effective than manipulating nearby varieties. We also suggest that planting a "whitefly sink" variety is unlikely to reduce B. tabaci SSA1 populations unless fields are managed to reduce B. tabaci densities using other integrative approaches.

Whitefly-transmitted viruses threatening cassava production in Africa

Emerging plant viruses are one of the greatest problems facing crop production worldwide, and have severe consequences in the developing world where subsistence farming is a major source of food production, and knowledge and resources for management are limited. In Africa, evolution of two viral disease complexes, cassava mosaic begomoviruses (CMBs) (Geminiviridae) and cassava brown streak viruses (CBSVs) (Potyviridae), have resulted in severe pandemics that continue to spread and threaten cassava production. Identification of genetically diverse and rapidly evolving CMBs and CBSVs, extensive genetic variation in the vector, Bemisia tabaci (Hemiptera: Aleyrodidae), and numerous secondary endosymbiont profiles that influence vector phenotypes suggest that complex local and regional vector-virus-plant-environment interactions may be driving the evolution and epidemiology of these viruses.

The Trouble with MEAM2: Implications of Pseudogenes on Species Delimitation in the Globally Invasive Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species Complex

Molecular species identification using suboptimal PCR primers can over-estimate species diversity due to coamplification of nuclear mitochondrial (NUMT) DNA/pseudogenes. For the agriculturally important whitefly Bemisia tabaci cryptic pest species complex, species identification depends primarily on characterization of the mitochondrial DNA cytochrome oxidase I (mtDNA COI) gene. The lack of robust PCR primers for the mtDNA COI gene can undermine correct species identification which in turn compromises management strategies. This problem is identified in the B. tabaci Africa/Middle East/Asia Minor clade which comprises the globally invasive Mediterranean (MED) and Middle East Asia Minor I (MEAM1) species, Middle East Asia Minor 2 (MEAM2), and the Indian Ocean (IO) species. Initially identified from the Indian Ocean island of Réunion, MEAM2 has since been reported from Japan, Peru, Turkey and Iraq. We identified MEAM2 individuals from a Peruvian population via Sanger sequencing of the mtDNA COI gene. In attempting to characterize the MEAM2 mitogenome, we instead characterized mitogenomes of MEAM1. We also report on the mitogenomes of MED, AUS, and IO thereby increasing genomic resources for members of this complex. Gene synteny (i.e., same gene composition and orientation) was observed with published B. tabaci cryptic species mitogenomes. Pseudogene fragments matching MEAM2 partial mtDNA COI gene exhibited low frequency single nucleotide polymorphisms that matched low copy number DNA fragments (<3%) of MEAM1 genomes, whereas presence of internal stop codons, loss of expected stop codons and poor primer annealing sites, all suggested MEAM2 as a pseudogene artifact and so not a real species.

Genome-wide analyses of the Bemisia tabaci species complex reveal contrasting patterns of admixture and complex demographic histories

Once considered a single species, the whitefly, Bemisia tabaci, is a complex of numerous morphologically indistinguishable species. Within the last three decades, two of its members (MED and MEAM1) have become some of the world's most damaging agricultural pests invading countries across Europe, Africa, Asia and the Americas and affecting a vast range of agriculturally important food and fiber crops through both feeding-related damage and the transmission of numerous plant viruses. For some time now, researchers have relied on a single mitochondrial gene and/or a handful of nuclear markers to study this species complex. Here, we move beyond this by using 38,041 genome-wide Single Nucleotide Polymorphisms, and show that the two invasive members of the complex are closely related species with signatures of introgression with a third species (IO). Gene flow patterns were traced between contemporary invasive populations within MED and MEAM1 species and these were best explained by recent international trade. These findings have profound implications for delineating the B. tabaci species status and will impact quarantine measures and future management strategies of this global pest.