Species-complex diversification and host-plant associations in Bemisia tabaci: A plant-defence, detoxification perspective revealed by RNA-Seq analyses

A comprehensive framework for the delimitation of species within the Bemisia tabaci cryptic complex, a global pest-species group

Identifying cryptic species poses a substantial challenge to both biologists and naturalists due to morphological similarities. Bemisia tabaci is a cryptic species complex containing more than 44 putative species; several of which are currently among the world's most destructive crop pests. Interpreting and delimiting the evolution of this species complex has proved problematic. To develop a comprehensive framework for species delimitation and identification, we evaluated the performance of distinct data sources both individually and in combination among numerous samples of the B. tabaci species complex acquired worldwide. Distinct datasets include full mitogenomes, single-copy nuclear genes, restriction site-associated DNA sequencing, geographic range, host speciation, and reproductive compatibility datasets. Phylogenetically, our well-supported topologies generated from three dense molecular markers highlighted the evolutionary divergence of species of the B. tabaci complex and suggested that the nuclear markers serve as a more accurate representation of B. tabaci species diversity. Reproductive compatibility datasets facilitated the identification of at least 17 different cryptic species within our samples. Native geographic range information provides a complementary assessment of species recognition, while the host range datasets provide low rate of delimiting resolution. We further summarized different data performances in species classification when compared with reproductive compatibility, indicating that combination of mtCOI divergence, nuclear markers, geographic range provide a complementary assessment of species recognition. Finally, we represent a model for understanding and untangling the cryptic species complexes based on the evidence from this study and previously published articles.

Functional response of Franklinothrips vespiformis (Thysanoptera: Aeolothripidae) to eggs and nymphs of Bemisia tabaci (Hemiptera: Aleyrodidae)

The Middle East Asia Minor 1 biotype of Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) is a greenhouse and field crop pest of global significance. The objective of this study was to assess the potential of the generalist predatory thrips, Franklinothrips vespiformis Crawford (Thysanoptera: Aeolothripidae), as a biological control agent for B. tabaci. This was achieved by determining the functional responses of F. vespiformis larvae and adults to the egg and nymphal stages of B. tabaci under laboratory conditions. Analyses consisted of 10 replicates of each predator and prey stage combination on bean leaf discs for a 24-h period. Following logistic regression analyses to determine the functional response type exhibited, response parameters were estimated with nonlinear least squares regression using Roger's equation. Results showed that F. vespiformis larvae and adults exhibited a Type II functional response when feeding on immature B. tabaci. The handling times (Th) of F. vespiformis larvae and adults were magnitudes higher for B. tabaci nymphs than they were for eggs, which were in part driven by the higher attack rates (a) observed on eggs. The maximum attack rate (T/Th) for B. tabaci eggs and nymphs exhibited by first-stage larvae, second-stage larvae, and adult F. vespiformis increased with increasing predator age. Results from this study suggest that F. vespiformis larvae and particularly adults are promising biological control agents for B. tabaci and are efficient predators at both low and high prey densities.

Whitefly-Plant Interactions: An Integrated Molecular Perspective

The rapid advances in available transcriptomic and genomic data and our understanding of the physiology and biochemistry of whitefly-plant interactions have allowed us to gain new and significant insights into the biology of whiteflies and their successful adaptation to host plants. In this review, we provide a comprehensive overview of the mechanisms that whiteflies have evolved to overcome the challenges of feeding on phloem sap. We also highlight the evolution and functions of gene families involved in host perception, evaluation, and manipulation; primary metabolism; and metabolite detoxification. We discuss the emerging themes in plant immunity to whiteflies, focusing on whitefly effectors and their sites of action in plant defense-signaling pathways. We conclude with a discussion of advances in the genetic manipulation of whiteflies and the potential that they hold for exploring the interactions between whiteflies and their host plants, as well as the development of novel strategies for the genetic control of whiteflies.

A nuclear receptor HR96-related gene underlies large trans-driven differences in detoxification gene expression in a generalist herbivore

The role, magnitude, and molecular nature of trans-driven expression variation underlying the upregulation of detoxification genes in pesticide resistant arthropod populations has remained enigmatic. In this study, we performed expression quantitative trait locus (eQTL) mapping (n = 458) between a pesticide resistant and a susceptible strain of the generalist herbivore and crop pest Tetranychus urticae. We found that a single trans eQTL hotspot controlled large differences in the expression of a subset of genes in different detoxification gene families, as well as other genes associated with host plant use. As established by additional genetic approaches including RNAi gene knockdown, a duplicated gene with a nuclear hormone receptor HR96-related ligand-binding domain was identified as causal for the expression differences between strains. The presence of a large family of HR96-related genes in T. urticae may enable modular control of detoxification and host plant use genes, facilitating this species' known and rapid evolution to diverse pesticides and host plants.

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

Begomoviruses: what is the secret(s) of their success?

Begomoviruses constitute an extremely successful group of emerging plant viruses transmitted by whiteflies of the Bemisia tabaci complex. Hosts include important vegetable, root, and fiber crops grown in the tropics and subtropics. Factors contributing to the ever-increasing diversity and success of begomoviruses include their predisposition to recombine their genomes, interaction with DNA satellites recruited throughout their evolution, presence of wild plants as a virus reservoir and a source of speciation, and extreme polyphagia and continuous movement of the insect vectors to temperate regions. These features as well as some controversial issues (replication in the insect vector, putative seed transmission, transmission by insects other than B. tabaci, and expansion of the host range to monocotyledonous plants) will be analyzed in this review.

Metabolization and sequestration of plant specialized metabolites in insect herbivores: Current and emerging approaches

Herbivorous insects encounter diverse plant specialized metabolites (PSMs) in their diet, that have deterrent, anti-nutritional, or toxic properties. Understanding how they cope with PSMs is crucial to understand their biology, population dynamics, and evolution. This review summarizes current and emerging cutting-edge methods that can be used to characterize the metabolic fate of PSMs, from ingestion to excretion or sequestration. It further emphasizes a workflow that enables not only to study PSM metabolism at different scales, but also to tackle and validate the genetic and biochemical mechanisms involved in PSM resistance by herbivores. This review thus aims at facilitating research on PSM-mediated plant-herbivore interactions.

Evolved transcriptional responses and their trade-offs after long-term adaptation of Bemisia tabaci to a marginally-suitable host

Although generalist insect herbivores can migrate and rapidly adapt to a broad range of host plants, they can face significant difficulties when accidentally migrating to novel and marginally-suitable hosts. What happens, both in performance and gene expression regulation, if these marginally-suitable hosts must be used for multiple generations before migration to a suitable host can take place, largely remains unknown. In this study, we established multigenerational colonies of the whitefly Bemisia tabaci, a generalist phloem-feeding species, adapted to a marginally-suitable host (habanero pepper) or an optimal host (cotton). We used reciprocal host tests to estimate the differences in performance of the populations on both hosts under optimal (30 oC) and mild-stressful (24 oC) temperature conditions, and documented the associated transcriptomic changes. The habanero pepper-adapted population greatly improved its performance on habanero pepper but did not reach its performance level on cotton, the original host. It also showed reduced performance on cotton, relative to the non-adapted population, and an antagonistic effect of the lower-temperature stressor. The transcriptomic data revealed that most of the expression changes, associated with long-term adaptation to habanero pepper, can be categorized as "evolved" with no initial plastic response. Three molecular functions dominated: enhanced formation of cuticle structural constituents, enhanced activity of oxidation-reduction processes involved in neutralization of phytotoxins and reduced production of proteins from the cathepsin B family. Taken together, these findings indicate that generalist insects can adapt to novel host plants by modifying the expression of a relatively small set of specific molecular functions.

Transcriptional responses are oriented towards different components of the rearing environment in two Drosophila sibling species

Background

The chance to compare patterns of differential gene expression in related ecologically distinct species can be particularly fruitful to investigate the genetics of adaptation and phenotypic plasticity. In this regard, a powerful technique such as RNA-Seq applied to ecologically amenable taxa allows to address issues that are not possible in classic model species. Here, we study gene expression profiles and larval performance of the cactophilic siblings Drosophila buzzatii and D. koepferae reared in media that approximate natural conditions and evaluate both chemical and nutritional components of the diet. These closely related species are complementary in terms of host-plant use since the primary host of one is the secondary of the other. D. koepferae is mainly a columnar cactus dweller while D. buzzatii prefers Opuntia hosts.

Results

Our comparative study shows that D. buzzatii and D. koepferae have different transcriptional strategies to face the challenges posed by their natural resources. The former has greater transcriptional plasticity, and its response is mainly modulated by alkaloids of its secondary host, while the latter has a more canalized genetic response, and its transcriptional plasticity is associated with the cactus species.

Conclusions

Our study unveils a complex pleiotropic genetic landscape in both species, with functional links that relate detox responses and redox mechanisms with developmental and neurobiological processes. These results contribute to deepen our understanding of the role of host plant shifts and natural stress driving ecological specialization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08745-9.

Characterization of transposable elements within the Bemisia tabaci species complex

Background

Whiteflies are agricultural pests that cause negative impacts globally to crop yields resulting at times in severe economic losses and food insecurity. The Bemisia tabaci whitefly species complex is the most damaging in terms of its broad crop host range and its ability to serve as vector for over 400 plant viruses. Genomes of whiteflies belonging to this species complex have provided valuable genomic data; however, transposable elements (TEs) within these genomes remain unexplored. This study provides the first accurate characterization of TE content within the B. tabaci species complex.

Results

This study identified that an average of 40.61% of the genomes of three whitefly species (MEAM1, MEDQ, and SSA-ECA) consists of TEs. The majority of the TEs identified were DNA transposons (22.85% average) while SINEs (0.14% average) were the least represented. This study also compared the TE content of the three whitefly genomes with three other hemipteran genomes and found significantly more DNA transposons and less LINEs in the whitefly genomes. A total of 63 TE superfamilies were identified to be present across the three whitefly species (39 DNA transposons, six LTR, 16 LINE, and two SINE). The sequences of the identified TEs were clustered which generated 5766 TE clusters. A total of 2707 clusters were identified as uniquely found within the whitefly genomes while none of the generated clusters were from both whitefly and non-whitefly TE sequences.

This study is the first to characterize TEs found within different B. tabaci species and has created a standardized annotation workflow that could be used to analyze future whitefly genomes.

Conclusion

This study is the first to characterize the landscape of TEs within the B. tabaci whitefly species complex. The characterization of these elements within the three whitefly genomes shows that TEs occupy significant portions of B. tabaci genomes, with DNA transposons representing the vast majority. This study also identified TE superfamilies and clusters of TE sequences of potential interest, providing essential information, and a framework for future TE studies within this species complex.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00270-6.

Dominance of Asia II 1 species of Bemisia tabaci in Pakistan and beyond

Globally, Whitefly (Bemisia tabaci) is one of the most important insect pests of crops that causes huge economical losses. The current study was designed to exclusively screen the B. tabaci species in the cotton field of Pakistan during 2017-2020 and have to conduct comparative analysis of B. tabaci species in Asia where Asia II 1 has been reported. A total of 5142 B. tabaci sequences of mitochondrial cytochrome oxidase 1 (mtCO1) from Asian countries were analyzed to determine the species and their distribution in the region. Our analysis over time and space showed that Asia II 1 has gradually dominated over Asia 1 in Punjab Province and over both Asia 1 and MEAM1 in Sindh Province. Asia has been divided into three regions i.e., South Asia (2524 sequences), Southeast Asia (757 sequences) and East Asia (1569 sequences) and dominance of different species of B. tabaci has been determined by calculating the relative percentage of each species. Interestingly, Asia II 1 has been found dominant in the neighboring region (northern zone) of India and also being dominant in its central zone. The dominance of Asia II 1 in Pakistan and northern India explains whitefly epidemic being reported in recent years.

Cytotype Affects the Capability of the Whitefly Bemisia tabaci MED Species To Feed and Oviposit on an Unfavorable Host Plant

The acquisition of nutritional obligate primary endosymbionts (P-symbionts) allowed phloemo-phageous insects to feed on plant sap and thus colonize novel ecological niches. P-symbionts often coexist with facultative secondary endosymbionts (S-symbionts), which may also influence their hosts' niche utilization ability. The whitefly Bemisia tabaci is a highly diversified species complex harboring, in addition to the P-symbiont "Candidatus Portiera aleyrodidarum," seven S-symbionts whose roles remain poorly understood. Here, we compare the phenotypic and metabolic responses of three B. tabaci lines differing in their S-symbiont community, reared on three different host plants, hibiscus, tobacco, or lantana, and address whether and how S-symbionts influence insect capacity to feed and produce offspring on those plants. We first show that hibiscus, tobacco, and lantana differ in their free amino acid composition. Insects' performance, as well as free amino acid profile and symbiotic load, were shown to be plant dependent, suggesting a critical role for the plant nutritional properties. Insect fecundity was significantly lower on lantana, indicating that it is the least favorable plant. Remarkably, insects reared on this plant show a specific amino acid profile and a higher symbiont density compared to the two other plants. In addition, this plant was the only one for which fecundity differences were observed between lines. Using genetically homogeneous hybrids, we demonstrate that cytotype (mitochondria and symbionts), and not genotype, is a major determinant of females' fecundity and amino acid profile on lantana. As cytotypes differ in their S-symbiont community, we propose that these symbionts may mediate their hosts' suitable plant range. IMPORTANCE Microbial symbionts are universal in eukaryotes, and it is now recognized that symbiotic associations represent major evolutionary driving forces. However, the extent to which symbionts contribute to their hosts' ecological adaptation and subsequent diversification is far from being fully elucidated. The whitefly Bemisia tabaci is a sap feeder associated with multiple coinfecting intracellular facultative symbionts. Here, we show that plant species simultaneously affect whiteflies' performance, amino acid profile, and symbiotic density, which could be partially explained by differences in plant nutritional properties. We also demonstrate that, on lantana, the least favorable plant used in our study, whiteflies' performance is determined by their cytotype. We propose that the host plant utilization in B. tabaci is influenced by its facultative symbiont community composition, possibly through its impact on the host dietary requirements. Altogether, our data provide new insights into the impact of intracellular microorganisms on their animal hosts' ecological niche range and diversification.

Performance of two species of whiteflies is unaffected by glucosinolate profile in Brassica plants

BACKGROUND

While plant glucosinolates are known to impart resistance to many insects, their role in the interactions between plants and many phloem-feeding insects such as whiteflies are poorly understood. The whitefly Bemisia tabaci complex comprises many cryptic species that differ in the ability to utilize Brassica plants. However, whether Brassica plants-specific traits such as glucosinolates determine differences of whiteflies in colonizing Brassica plants remains in question.

RESULTS

We first observed performance of two whitefly species MEAM1 and Asia II 3, which differ obviously in their ability to colonize Brassica plants, on four cultivars of three Brassica species that vary in glucosinolate profile. We found that the life history characteristics of each of the two whitefly species seems to be only marginally affected by cultivar. We next used wild-type Arabidopsis plants and mutants defective in glucosinolate biosynthesis or hydrolysis to explore the effects of glucosinolates on the whitefly. We found that fecundity and development of immature stages of neither of the two whitefly species differ significantly between wild-type and mutants.

CONCLUSION

The data suggest that glucosinolates may have little effect on the oviposition by adults and the survival and development of immature stages of MEAM1 and Asia II 3 whiteflies. The marked differences in colonizing Brassica crops between the two whitefly species are likely due to plant traits other than glucosinolates. © 2021 Society of Chemical Industry.

Comparative genome and transcriptome analyses reveal innate differences in response to host plants by two color forms of the two-spotted spider mite Tetranychus urticae

Background

The two-spotted spider mite, Tetranychus urticae, is a major agricultural pest with a cosmopolitan distribution, and its polyphagous habits provide a model for investigating herbivore-plant interactions. There are two body color forms of T. urticae with a different host preference. Comparative genomics and transcriptomics are used here to investigate differences in responses of the forms to host plants at the molecular level. Biological responses of the two forms sourced from multiple populations are also presented.

Results

We carried out principal component analysis of transcription changes in three red and three green T. urticae populations feeding on their original host (common bean), and three hosts to which they were transferred: cotton, cucumber and eggplant. There were differences among the forms in gene expression regardless of their host plant. In addition, different changes in gene expression were evident in the two forms when responding to the same host transfer. We further compared biological performance among populations of the two forms after feeding on each of the four hosts. Fecundity of 2-day-old adult females showed a consistent difference between the forms after feeding on bean. We produced a 90.1-Mb genome of the red form of T. urticae with scaffold N50 of 12.78 Mb. Transcriptional profiles of genes associated with saliva, digestion and detoxification showed form-dependent responses to the same host and these genes also showed host-specific expression effects.

Conclusions

Our research revealed that forms of T. urticae differ in host-determined transcription responses and that there is form-dependent plasticity in the transcriptomic responses. These differences may facilitate the extreme polyphagy shown by spider mites, although fitness differences on hosts are also influenced by population differences unrelated to color form.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07894-7.

Assessing the diversity of whiteflies infesting cassava in Brazil

Background

The necessity of a competent vector for transmission is a primary ecological factor driving the host range expansion of plant arthropod-borne viruses, with vectors playing an essential role in disease emergence. Cassava begomoviruses severely constrain cassava production in Africa. Curiously, begomoviruses have never been reported in cassava in South America, the center of origin for this crop. It has been hypothesized that the absence of a competent vector in cassava is the reason why begomoviruses have not emerged in South America.

Methods

We performed a country-wide whitefly diversity study in cassava in Brazil. Adults and/or nymphs of whiteflies were collected from sixty-six cassava fields in the main agroecological zones of the country. A total of 1,385 individuals were genotyped based on mitochondrial cytochrome oxidase I sequences.

Results

A high species richness was observed, with five previously described species and two putative new ones. The prevalent species were Tetraleurodes acaciae and Bemisia tuberculata, representing over 75% of the analyzed individuals. Although we detected, for the first time, the presence of Bemisia tabaci Middle East-Asia Minor 1 (BtMEAM1) colonizing cassava in Brazil, it was not prevalent. The species composition varied across regions, with fields in the Northeast region showing a higher diversity. These results expand our knowledge of whitefly diversity in cassava and support the hypothesis that begomovirus epidemics have not occurred in cassava in Brazil due to the absence of competent vector populations. However, they indicate an ongoing adaptation process of BtMEAM1 to cassava, increasing the likelihood of begomovirus emergence in this crop.

Spotlight on the Roles of Whitefly Effectors in Insect-Plant Interactions

The Bemisia tabaci species complex (whitefly) causes enormous agricultural losses. These phloem-feeding insects induce feeding damage and transmit a wide range of dangerous plant viruses. Whiteflies colonize a broad range of plant species that appear to be poorly defended against these insects. Substantial research has begun to unravel how phloem feeders modulate plant processes, such as defense pathways, and the central roles of effector proteins, which are deposited into the plant along with the saliva during feeding. Here, we review the current literature on whitefly effectors in light of what is known about the effectors of phloem-feeding insects in general. Further analysis of these effectors may improve our understanding of how these insects establish compatible interactions with plants, whereas the subsequent identification of plant defense processes could lead to improved crop resistance to insects. We focus on the core concepts that define the effectors of phloem-feeding insects, such as the criteria used to identify candidate effectors in sequence-mining pipelines and screens used to analyze the potential roles of these effectors and their targets in planta. We discuss aspects of whitefly effector research that require further exploration, including where effectors localize when injected into plant tissues, whether the effectors target plant processes beyond defense pathways, and the properties of effectors in other insect excretions such as honeydew. Finally, we provide an overview of open issues and how they might be addressed.

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.

Activation and detoxification of cassava cyanogenic glucosides by the whitefly Bemisia tabaci

Two-component plant defenses such as cyanogenic glucosides are produced by many plant species, but phloem-feeding herbivores have long been thought not to activate these defenses due to their mode of feeding, which causes only minimal tissue damage. Here, however, we report that cyanogenic glycoside defenses from cassava (Manihot esculenta), a major staple crop in Africa, are activated during feeding by a pest insect, the whitefly Bemisia tabaci, and the resulting hydrogen cyanide is detoxified by conversion to beta-cyanoalanine. Additionally, B. tabaci was found to utilize two metabolic mechanisms to detoxify cyanogenic glucosides by conversion to non-activatable derivatives. First, the cyanogenic glycoside linamarin was glucosylated 1–4 times in succession in a reaction catalyzed by two B. tabaci glycoside hydrolase family 13 enzymes in vitro utilizing sucrose as a co-substrate. Second, both linamarin and the glucosylated linamarin derivatives were phosphorylated. Both phosphorylation and glucosidation of linamarin render this plant pro-toxin inert to the activating plant enzyme linamarase, and thus these metabolic transformations can be considered pre-emptive detoxification strategies to avoid cyanogenesis.

Whitefly Network Analysis Reveals Gene Modules Involved in Host Plant Selection, Development and Evolution

Whiteflies are Hemipterans that typically feed on the undersides of plant leaves. They cause severe damage by direct feeding as well as transmitting plant viruses to a wide range of plants. However, it remains largely unknown which genes play a key role in development and host selection. In this study, weighted gene co-expression network analysis was applied to construct gene co-expression networks in whitefly. Nineteen gene co-expression modules were detected from 15560 expressed genes of whitefly. Combined with the transcriptome data of salivary glands and midgut, we identified three gene co-expression modules related to host plant selection. These three modules contain genes related to host-plant recognition, such as detoxification genes, chemosensory genes and some salivary gland-associated genes. Results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses elucidated the following pathways involved in these modules: lysosome, metabolic and detoxification pathways. The modules related to the development contain two co-expression modules; moreover, the genes were annotated to the development of chitin-based cuticle. This analysis provides a basis for future functional analysis of genes involved in host-plant recognition.

Virus-induced phytohormone dynamics and their effects on plant-insect interactions

Attacks on plants by both viruses and their vectors is common in nature. Yet the dynamics of the plant-virus-vector tripartite system, in particular the effects of viral infection on plant-insect interactions, have only begun to emerge in the last decade. Viruses can modulate the interactions between insect vectors and plants via the jasmonate, salicylic acid and ethylene phytohormone pathways, resulting in changes in fitness and viral transmission capacity of their insect vectors. Virus infection of plants may also modulate other phytohormones, such as auxin, gibberellins, cytokinins, brassinosteroids and abscisic acid, with yet undefined consequences on plant-insect interactions. Moreover, virus infection in plants may incur changes to other plant traits, such as nutrition and secondary metabolites, that potentially contribute to virus-associated, phytohormone-mediated manipulation of plant-insect interactions. In this article, we review the research progress, discuss issues related to the complexity and variability of the viral modulation of plant interactions with insect vectors, and suggest future directions of research in this field.

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.

Within-Season Changes in Land-Use Impact Pest Abundance in Smallholder African Cassava Production Systems

Cassava (Manihot esculenta Crantz), an important commercial and food security crop in East and Central Africa, continues to be adversely affected by the whitefly Bemisia tabaci. In Uganda, changes in smallholder farming landscapes due to crop rotations can impact pest populations but how these changes affect pest outbreak risk is unknown. We investigated how seasonal changes in land-use have affected B. tabaci population dynamics and its parasitoids. We used a large-scale field experiment to standardize the focal field in terms of cassava age and cultivar, then measured how Bemisia populations responded to surrounding land-use change. Bemisia tabaci Sub-Saharan Africa 1 (SSA1) was identified using molecular diagnostics as the most prevalent species and the same species was also found on surrounding soybean, groundnut, and sesame crops. We found that an increase in the area of cassava in the 3-7-month age range in the landscape resulted in an increase in the abundance of the B. tabaci SSA1 on cassava. There was a negative relationship between the extent of non-crop vegetation in the landscape and parasitism of nymphs suggesting that these parasitoids do not rely on resources in the non-crop patches. The highest abundance of B. tabaci SSA1 nymphs in cassava fields occurred at times when landscapes had large areas of weeds, low to moderate areas of maize, and low areas of banana. Our results can guide the development of land-use strategies that smallholder farmers can employ to manage these pests.

The molecular mechanisms that determine different degrees of polyphagy in the Bemisia tabaci species complex

The whitefly Bemisia tabaci is a closely related group of >35 cryptic species that feed on the phloem sap of a broad range of host plants. Species in the complex differ in their host-range breadth, but the mechanisms involved remain poorly understood. We investigated, therefore, how six different B. tabaci species cope with the environmental unpredictability presented by a set of four common and novel host plants. Behavioral studies indicated large differences in performances on the four hosts and putative specialization of one of the species to cassava plants. Transcriptomic analyses revealed two main insights. First, a large set of genes involved in metabolism (>85%) showed differences in expression between the six species, and each species could be characterized by its own unique expression pattern of metabolic genes. However, within species, these genes were constitutively expressed, with a low level of environmental responsiveness (i.e., to host change). Second, within each species, sets of genes mainly associated with the super-pathways "environmental information processing" and "organismal systems" responded to the host switching events. These included genes encoding for proteins involved in sugar homeostasis, signal transduction, membrane transport, and immune, endocrine, sensory and digestive responses. Our findings suggested that the six B. tabaci species can be divided into four performance/transcriptomic "Types" and that polyphagy can be achieved in multiple ways. However, polyphagy level is determined by the specific identity of the metabolic genes/pathways that are enriched and overexpressed in each species (the species' individual metabolic "tool kit").

Short term transcriptional responses of P450s to phytochemicals in insects and mites

Cytochrome P450 monooxygenases (P450s) play a key role in the detoxification of phytochemicals in arthropod herbivores. We present here an overview of recent progress in understanding the breadth and specificity of gene expression plasticity of P450s in response to phytochemicals. We discuss experimental setups and new findings in mechanisms of P450 regulation. Whole genome transcriptomic analysis of arthropod herbivores, either after direct administration of phytochemicals or after host plant shifts, allowed to integrate various levels of chemical complexity and lead to the unbiased identification of responsive P450 genes. However, despite progress in identification of inducible P450s, the link between induction and metabolism is still largely unexplored, and to what extent the overall response is biologically functional should be further investigated. In the near future, such studies will be more straightforward as forward and reverse genetic tools become more readily available.

Genetic diversity of whitefly (Bemisia spp.) on crop and uncultivated plants in Uganda: implications for the control of this devastating pest species complex in Africa

Over the past three decades, highly increased whitefly (Bemisia tabaci) populations have been observed on the staple food crop cassava in eastern Africa and associated with ensuing viral disease pandemics and food insecurity. Increased whitefly numbers have also been observed in other key agricultural crops and weeds. Factors behind the population surges on different crops and their interrelationships are unclear, although in cassava they have been associated with specific populations within the Bemisia tabaci species complex known to infest cassava crops in Africa. This study carried out an in-depth survey to understand the distribution of B. tabaci populations infesting crops and uncultivated plant hosts in Uganda, a centre of origin for this pest complex. Whitefly samples were collected from 59 identified plant species and 25 unidentified weeds in a countrywide survey. Identities of 870 individual adult whiteflies were determined through mitochondrial cytochrome oxidase 1 sequences (651 bp) in the 3' barcode region used for B. tabaci systematics. Sixteen B. tabaci and five related whitefly putative species were identified based on > 4.0% nucleotide divergence, of which three are proposed as novel B. tabaci putative species and four as novel closely related whitefly species. The most prevalent whiteflies were classified as B. tabaci MED-ASL (30.5% of samples), sub-Saharan Africa 1 (SSA1, 22.7%) and Bemisia Uganda1 (12.1%). These species were also indicated to be the most polyphagous occurring on 33, 40 and 25 identified plant species, respectively. Multiple (≥ 3) whitefly species occurred on specific crops (bean, eggplant, pumpkin and tomato) and weeds (Sida acuta and Ocimum gratissimum). These plants may have increased potential to act as reservoirs for mixed infections of whitefly-vectored viruses. Management of whitefly pest populations in eastern Africa will require an integration of approaches that consider their degree of polyphagy and a climate that enables the continuous presence of crop and uncultivated plant hosts.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10340-021-01355-6.

Multiple invasions of a generalist herbivore-Secondary contact between two divergent lineages of Nezara viridula Linnaeus in Australia

The presence of distinct evolutionary lineages within herbivorous pest insect taxa requires close attention. Scientific understanding, biosecurity planning and practice, and pest management decision-making each suffer when such situations remain poorly understood. The pest bug Nezara viridula Linnaeus has been recorded from numerous host plants and has two globally distributed mitochondrial (mtDNA) lineages. These mtDNA lineages co-occur in few locations globally, and the consequences of their divergence and recent secondary contact have not been assessed. We present evidence that both mtDNA lineages of N. viridula are present in Australia and their haplotype groups have a mostly separate distribution from one another. The north-western population has only Asian mtDNA haplotypes, and the population with an eastern distribution is characterized mostly by European mtDNA haplotypes. Haplotypes of both lineages were detected together at only one site in the north of eastern Australia, and microsatellite data indicate that this secondary contact has resulted in mating across the lineages. Admixture and the movement of mtDNA haplotypes outside of this limited area of overlap has not, however, been extensive. Some degree of mating incompatibility or differences in the climatic requirements and tolerances of the two lineages, and perhaps a combination of these influences, might limit introgression and the movement of individuals, but this needs to be tested. This work provides the foundation for further ecological investigation of the lineages of N. viridula, particularly the consequences of admixture on the ecology of this widespread pest. We propose that for now, the Asian and European lineages of N. viridula would best be investigated as subspecies, so that "pure" and admixed populations of this bug can each be considered directly with respect to management and research priorities.

Transcriptome analysis and comparison reveal divergence between the Mediterranean and the greenhouse whiteflies

Both the Mediterranean (MED) species of the Bemisia tabaci whitefly complex and the greenhouse whitefly (Trialeurodes vaporariorum, TV) are important agricultural pests. The two species of whiteflies differ in many aspects such as morphology, geographical distribution, host plant range, plant virus transmission, and resistance to insecticides. However, the molecular basis underlying their differences remains largely unknown. In this study, we analyzed the genetic divergences between the transcriptomes of MED and TV. In total, 2,944 pairs of orthologous genes were identified. The average identity of amino acid sequences between the two species is 93.6%. The average nonsynonymous (Ka) and synonymous (Ks) substitution rates and the ratio of Ka/Ks of the orthologous genes are 0.0389, 2.23 and 0.0204, respectively. The low average Ka/Ks ratio indicates that orthologous genes tend to be under strong purified selection. The most divergent gene classes are related to the metabolisms of xenobiotics, cofactors, vitamins and amino acids, and this divergence may underlie the different biological characteristics between the two species of whiteflies. Genes of differential expression between the two species are enriched in carbohydrate metabolism and regulation of autophagy. These findings provide molecular clues to uncover the biological and molecular differences between the two species of whiteflies.

Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control

Cryptic taxa have often been observed in the form of host‐associated species that diverged as the result of adaptation to alternate host plants. Untangling cryptic diversity in species complexes that encompass invasive species is a mandatory task for pest management. Moreover, investigating the evolutionary history of a species complex may help to understand the drivers of their diversification. The mealybug Hypogeococcus pungens was believed to be a polyphagous species from South America and has been reported as a pest devastating native cacti in Puerto Rico, also threatening cactus diversity in the Caribbean and North America. There is neither certainty about the identity of the pest nor the source population from South America. Recent studies pointed to substantial genetic differentiation among local populations, suggesting that H. pungens is a species complex. In this study, we used a combination of genome‐wide SNPs and mtDNA variation to investigate species diversity within H. pungens sensu lato to establish host plant ranges of each one of the putative members of the complex, to evaluate whether the pattern of host plant association drove diversification in the species complex, and to determine the source population of the Puerto Rican cactus pest. Our results suggested that H. pungens comprises at least five different species, each one strongly associated with specific host plants. We also established that the Puerto Rican cactus pest derives from southeastern Brazilian mealybugs. This is an important achievement because it will help to design reliable strategies for biological control using natural enemies of the pest from its native range.

Nutritional Relationship between Bemisia tabaci and Its Primary Endosymbiont, Portiera aleyrodidarum, during Host Plant Acclimation

Simple Summary

Plant sap-sucking insects commonly have established mutualistic relationships with bacteria that live within their bodies and often provide nutrients that are lacking in the insect’s diet. The sweet potato whitefly (Bemisia tabaci) harbors one primary and up to seven secondary endosymbiotic bacteria. The primary endosymbiont of B. tabaci is already known to play a critical role in providing necessary nutrients for B. tabaci. Our objective was to study the relationship among B. tabaci, its primary endosymbiont, and the host plant through the effects of host plant shifting and acclimation, that is, physiological adjustments as an insect becomes accustomed to a new host plant over several generations. The results showed that host shifting from Chinese kale to cotton plants led to a decrease in the fecundity of B. tabaci in the first generation, which was restored after 10 generations of acclimation, and that its developmental time was also decreased by the tenth generation. Furthermore, essential amino acid biosynthesis genes of its primary endosymbiont were differentially regulated after B. tabaci had become acclimated to cotton plants. We speculate that the primary endosymbiont has a close nutritional relationship with B. tabaci during host plant acclimation.

Abstract

Plant sap-sucking insects commonly have established mutualistic relationships with endosymbiotic bacteria that can provide nutrients lacking in their diet. Bemisia tabaci harbors one primary endosymbiont, Portiera aleyrodidarum, and up to seven secondary endosymbionts, including Hamiltonella defensa and Rickettsia sp. Portiera aleyrodidarum is already known to play a critical role in providing necessary nutrients for B. tabaci. In the present study, the relationship among B. tabaci, its primary endosymbiont, and the host plant were examined through the effects of host plant shifting and acclimation. Bemisia tabaci was transferred from Chinese kale to four different host plants, and the effects on both its performance and the expression levels of nutrient-related genes of P. aleyrodidarum were analyzed. The results showed that host shifting from Chinese kale to cotton plants led to a decrease in the performance of B. tabaci in the first generation, which was restored after 10 generations of acclimation. Furthermore, the expression levels of essential amino acid biosynthesis genes of P. aleyrodidarum were found to be differentially regulated after B. tabaci had acclimated to the cotton plants. Host plant shifting and acclimation to cucumber, poinsettia, and tomato plants did not affect the fecundity of B. tabaci and the expression levels of most examined genes. We speculate that P. aleyrodidarum may help B. tabaci improve its performance and acclimate to new hosts and that P. aleyrodidarum has a close nutritional relationship with its host during host plant acclimation.

Insecticide resistance and fitness cost in Bemisia tabaci (Hemiptera: Aleyrodidae) invasive and resident species in La Réunion Island

BACKGROUND

Global and intensive use of insecticides has led to the emergence and rapid evolution of resistance in the major pest Bemisia tabaci (Gennadius). In La Réunion, an island of the South West Indian Ocean, three whitefly species coexist, two of which are predominant, the indigenous Indian Ocean (IO) and the invasive Middle East Asia Minor 1 (MEAM1) species. To assess the resistance level of both of these species to acetamiprid and pymetrozine, whitefly populations were sampled at 15 collection sites located all over the island in agroecosystems and natural areas, and tested using leaf-dip bioassays. We also investigated the potential cost of resistance to acetamiprid by measuring six fitness-related traits for MEAM1 populations that displayed different resistance levels.

RESULTS

IO was mainly found in natural areas and was susceptible to both acetamiprid and pymetrozine. MEAM1 populations displayed evidence of high resistance to pymetrozine, whereas resistance to acetamiprid was more variable. No fitness-related costs were associated with this resistance in MEAM1 populations.

CONCLUSION

This is the first assessment of the susceptibility to insecticides for B. tabaci IO species. For the time being, no resistance to the tested insecticides has evolved in this species despite (i) its presence in agroecosystems and their surroundings, and (ii) its close proximity to, and possible hybridization with, the MEAM1 species. In contrast, with continuous selection pressure of insecticide treatments and in the absence of fitness cost to resistance, the invasive exotic species MEAM1 will continue to threaten agriculture in La Réunion. © 2019 Society of Chemical Industry.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Inside out: microbiota dynamics during host-plant adaptation of whiteflies

While most insect herbivores are selective feeders, a small proportion of them feed on a wide range of plants. This polyphagous habit requires overcoming a remarkable array of defenses, which often necessitates an adaptation period. Efforts for understanding the mechanisms involved mostly focus on the insect’s phenotypic plasticity. Here, we hypothesized that the adaptation process might partially rely on transient associations with bacteria. To test this, we followed in a field-like experiment, the adaptation process of Bemisia tabaci, a generalist sap feeder, to pepper (a less-suitable host), after switching from watermelon (a suitable host). Amplicon sequencing of 16S rRNA transcripts from hundreds of dissected guts revealed the presence of active “core” and “transient” bacterial communities, dominated by the phyla Proteobacteria, Actinobacteria, and Firmicutes, and increasing differences between populations grown on watermelon and pepper. Insects grown on pepper for over two generations presented a significant increase in specific genera, mainly Mycobacterium, with a predicted enrichment in degradative pathways of xenobiotics and secondary metabolites. This result correlated with a significant increase in the insect’s survival on pepper. Taken together, our findings suggest that gut-associated bacteria can provide an additional flexible metabolic “tool-box” to generalist sap feeders for facilitating a quick host switching process.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Generalism in Nature…The Great Misnomer: Aphids and Wasp Parasitoids as Examples

In the present article we discuss why, in our view, the term ‘generalism’ to define the dietary breadth of a species is a misnomer and should be revised by entomologists/ecologists with the more exact title relating to the animal in question’s level of phagy—mono-, oligo, or polyphagy. We discard generalism as a concept because of the indisputable fact that all living organisms fill a unique ecological niche, and that entry and exit from such niches are the acknowledged routes and mechanisms driving ecological divergence and ultimately speciation. The term specialist is probably still useful and we support its continuing usage simply because all species and lower levels of evolutionary diverge are indeed specialists to a large degree. Using aphids and parasitoid wasps as examples, we provide evidence from the literature that even some apparently highly polyphagous agricultural aphid pest species and their wasp parasitoids are probably not as polyphagous as formerly assumed. We suggest that the shifting of plant hosts by herbivorous insects like aphids, whilst having positive benefits in reducing competition, and reducing antagonists by moving the target organism into ‘enemy free space’, produces trade-offs in survival, involving relaxed selection in the case of the manicured agro-ecosystem.

Species-complex diversification and host-plant associations in Bemisia tabaci: A plant-defence, detoxification perspective revealed by RNA-Seq analyses

Insect–plant associations and their role in diversification are mostly studied in specialists. Here, we aimed to identify macroevolution patterns in the relationships between generalists and their host plants that have the potential to promote diversification. We focused on the Bemisia tabaci species complex containing more than 35 cryptic species. Mechanisms for explaining this impressive diversification have focused so far on allopatric forces that assume a common, broad, host range. We conducted a literature survey which indicated that species in the complex differ in their host range, with only few showing a truly broad one. We then selected six species, representing different phylogenetic groups and documented host ranges. We tested whether differences in the species expression profiles of detoxification genes are shaped more by their phylogenetic relationships or by their ability to successfully utilize multiple hosts, including novel ones. Performance assays divided the six species into two groups of three, one showing higher performance on various hosts than the other (the lower performance group). The same grouping pattern appeared when the species were clustered according to their expression profiles. Only species placed in the lower performance group showed a tendency to lower the expression of multiple genes. Taken together, these findings bring evidence for the existence of a common detoxification “machinery,” shared between species that can perform well on multiple hosts. We raise the possibility that this “machinery” might have played a passive role in the diversification of the complex, by allowing successful migration to new/novel environments, leading, in some cases, to fragmentation and speciation.