Whole-genome single nucleotide polymorphism and mating compatibility studies reveal the presence of distinct species in sub-Saharan Africa Bemisia tabaci whiteflies

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.

Interaction of Whitefly Effector G4 with Tomato Proteins Impacts Whitefly Performance

The phloem-feeding insect Bemisia tabaci is an important pest, responsible for the transmission of several crop-threatening virus species. While feeding, the insect secretes a cocktail of effectors to modulate plant defense responses. Here, we present a set of proteins identified in an artificial diet on which B. tabaci was salivating. We subsequently studied whether these candidate effectors can play a role in plant immune suppression. Effector G4 was the most robust suppressor of an induced- reactive oxygen species (ROS) response in Nicotiana benthamiana. In addition, G4 was able to suppress ROS production in Solanum lycopersicum (tomato) and Capsicum annuum (pepper). G4 localized predominantly in the endoplasmic reticulum in N. benthamiana leaves and colocalized with two identified target proteins in tomato: REF-like stress related protein 1 (RSP1) and meloidogyne-induced giant cell protein DB141 (MIPDB141). Silencing of MIPDB141 in tomato reduced whitefly fecundity up to 40%, demonstrating that the protein is involved in susceptibility to B. tabaci. Together, our data demonstrate that effector G4 impairs tomato immunity to whiteflies by interfering with ROS production and via an interaction with tomato susceptibility protein MIPDB141. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Developing a predictive model for an emerging epidemic on cassava in sub-Saharan Africa

The agricultural productivity of smallholder farmers in sub-Saharan Africa (SSA) is severely constrained by pests and pathogens, impacting economic stability and food security. An epidemic of cassava brown streak disease, causing significant yield loss, is spreading rapidly from Uganda into surrounding countries. Based on sparse surveillance data, the epidemic front is reported to be as far west as central DRC, the world's highest per capita consumer, and as far south as Zambia. Future spread threatens production in West Africa including Nigeria, the world's largest producer of cassava. Using innovative methods we develop, parameterise and validate a landscape-scale, stochastic epidemic model capturing the spread of the disease throughout Uganda. The model incorporates real-world management interventions and can be readily extended to make predictions for all 32 major cassava producing countries of SSA, with relevant data, and lays the foundations for a tool capable of informing policy decisions at a national and regional scale.

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.

Genetic diversity, distribution, and structure of Bemisia tabaci whitefly species in potential invasion and hybridization regions of East Africa

Outbreaks of whitefly, Bemisia tabaci species in East and Central Africa, have become increasingly prevalent during the previous 25 years and are responsible for driving the spread of plant-virus diseases, such as cassava mosaic disease and cassava brown steak disease. Epidemics of these diseases have expanded their ranges over the same period, spreading from Uganda into other sub-Saharan African countries. It was hypothesised that a highly abundant 'invader' population of B. tabaci was responsible for spreading these diseases from Uganda to neighbouring countries and potentially hybridising with the resident cassava B. tabaci populations. Here, we test this hypothesis by investigating the molecular identities of the highly abundant cassava B. tabaci populations from their supposed origin in Uganda, to the northern, central, eastern and coastal regions of Tanzania. Partial mitochondrial cytochrome oxidase I (mtCOI) barcoding sequences and nuclear microsatellite markers were used to analyse the population genetic diversity and structure of 2734 B. tabaci collected from both countries and in different agroecological zones. The results revealed that: (i) the putative SSA1 species is structured according to countries, so differ between them. (ii) Restricted gene flow occurred between SSA1-SG3 and both other SSA1 subgroups (SG1 and SG2), even in sympatry, demonstrating strong barriers to hybridization between those genotypes. (iii) Not only B. tabaci SSA1-(SG1 and SG2) was found in highly abundant (outbreak) numbers, but B. tabaci SSA1-SG3 and the Indian Ocean (IO) species were also recorded in high numbers in several sites in Tanzania. (iv) The SSA1-(SG1 and SG2) species was distributed in both countries, but in Tanzania, the B. tabaci IO and SSA1-SG3 species predominated. These data confirm that multiple, local Tanzanian B. tabaci species produce highly abundant populations, independent of the spread of the putative invasive B. tabaci SSA1-(SG1 and SG2) populations.

The Bemisia tabaci cryptic (sibling) species group - imperative for a taxonomic reassessment

The taxonomy of the Bemisia tabaci cryptic species group remains a challenge due to the lack of morphological differentiation and porous species boundaries among its members. Additionally, it is unclear whether B. tabaci consists of several species in evolutionary stasis with limited morphological change or is the result of a recent adaptive radiation characterized by great ecological diversity but little morphological divergence. Here, a historical overview of the development of the nomenclature used to classify B. tabaci is provided covering changes after synonymizing several species in 1957 until recent insights gained from whole-genome sequencing data. The article discusses the limitations of using a 3.5% mtCOI threshold and argues that a 1% nuclear divergence cutoff better reflects ecological and biogeographic species boundaries. Finally, a plan of action is outlined for naming B. tabaci species using a Latin binomial system in accordance with the International Comission on Zoological Nomenclature (ICZN) regulations.

Distribution of Bemisia tabaci in different agro-ecological regions in Uganda and the threat of vector-borne pandemics into new cassava growing areas

Previous studies in sub-Saharan Africa have showed the spread of cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) pandemics into different cassava growing regions by high Bemisia tabaci populations. Studies did indicate that there were stark differences in some whitefly species, yet they have not looked extensively across agroecologies. Members of B. tabaci species complex termed sub-Saharan Africa 1 (SSA1) and SSA2 have been linked to the spread of CMD and CBSD viruses. During the period of a severe CMD pandemic in the 1990s, SSA2 was the most predominant until the resurgence of SSA1, particularly SSA1-subgroup1 (SSA1-SG1) from the early 2000s to date. Cassava being a drought resilient crop has become an important food security crop and has been introduced into new areas and regions. Considering the role B. tabaci in the spread of cassava virus pandemics into neighboring regions, we investigated the genetic diversity and distribution of B. tabaci in nine different agro-ecological regions of Uganda in 2017. Adult whiteflies were collected from cassava and 33 other host plants from cassava-growing areas, those with limited cassava and areas with no cassava, where it is being introduced as a food security crop. The partial sequences of the mitochondrial cytochrome oxidase 1 (mtCO1) gene (657 bp) were used to determine the phylogenetic relationships between the sampled B. tabaci. Cassava B. tabaci SSA1 (-SG1, -SG2, -Hoslundia (previously called SSA1-SG1/2), -SG3), SSA2 and SSA3; non-cassava B. tabaci SSA6, SSA10, SSA11, SSA12, SSA13, MED-ASL, MED-Q1, MEAM1, Indian Ocean; and other Bemisia species, Bemisia afer and Bemisia Uganda1 were identified in the study. SSA3, one of the key B. tabaci species that occurs on cassava in West Africa, was identified for the first time in Uganda. The SSA1-SG1 was widely distributed, predominated on cassava and was found on 17 other host-plants. The ability of SSA1-SG1 to exist in environments with limited or no cassava growing poses the risk of continued spread of virus pandemics. Therefore, measures must be put in place to prevent the introduction of diseased materials into new areas, since the vectors exist.

Genetic diversity of cryptic species of Bemisia tabaci in Asia

Bemisia tabaci is a species complex consisting of various genetically different cryptic species worldwide. To understand the genetic characteristics and geographic distribution of cryptic species of B. tabaci in Asia, we conducted an extensive collection of B. tabaci samples in ten Asian countries (Bangladesh, Indonesia, Japan, Korea, Myanmar, Nepal, Philippines, Singapore, Taiwan, and Vietnam) from 2013 to 2020 and determined 56 different partial sequences of mitochondrial cytochrome oxidase subunit I (COI) DNA. In addition, information on 129 COI sequences of B. tabaci identified from 16 Asian countries was downloaded from the GenBank database. Among the total 185 COI sequences of B. tabaci, the sequence variation reached to 19.68%. In addition, there were 31 cryptic species updated from 16 countries in Asia, that is, Asia I, Asia I India, Asia II (1-13), Asia III, Asia IV, Asia V, China 1-6, MEAM (1, 2, K), MED, Australia/Indonesia, Japan (1 and 2). Further, MED cryptic species consisted of 2 clades, Q1 and Q2. This study provides updated information to understand the genetic variation and geographic diversity of B. tabaci in Asia.

Spatio-temporal changes in endosymbiont diversity and composition in the African cassava whitefly, Bemisia tabaci SSA1

Sap-sucking insects, including whiteflies, are amongst the most devastating and widely distributed organisms on the planet. They are often highly invasive and endosymbiont communities within these insects help them adapt to new or changing environments. Bemisia tabaci (Gennadius; Hemiptera: Aleyrodidae) whitefly species are vectors of more than 500 known plant-viruses and harbour highly diverse endosymbionts communities. To date, however, whitefly–endosymbiont interactions, community structure and their spatio-temporal changes are still poorly understood. In this study, we investigated the spatio-temporal changes in the composition and diversity of bacterial endosymbionts in the agricultural crop pest whitefly species, Bemisia tabaci sub-Saharan Africa 1-subgroup 1 and 2 (SSA1-SG1 and SSA1-SG2). 16S rRNA amplicon sequencing analysis was carried out to characterise endosymbiont compositionsin field-collected SSA1 (SSA1-SG1 and SSA1-SG2) populations infesting cassava in Uganda in 1997 and 2017. We detected Portiera, Arsenophonus, Wolbachia, Hamiltonella and Hemipteriphilus, with Arsenophonus and Wolbachia infections being predominant. Hemipteriphilus and Hamiltonella frequencies were very low and were detected in seven and two samples, respectively. Bacterial diversity based on three independent parameters including Simpson index, number of haplotypes and Bray–Curtis dissimilarity matrix was significantly higher in 1997 than in 2017. This period also coincided with the advent of super-abundant cassava-whitefly populations on cassava crops in Uganda. We discuss how endosymbionts may influence the biology and behaviour of whiteflies leading to population explosions.

Cytochrome P450 Gene, CYP6CX3, Is Involved in the Resistance to Cyantraniliprole in Bemisia tabaci

Bemisia tabaci is an important agricultural sucking pest, and it develops serious resistance to various insecticides. Although cytochrome P450 was involved in the resistance to cyantraniliprole, limited studies have been conducted on B. tabaci. In the present study, piperonyl butoxide significantly increased the toxicity of cyantraniliprole. P450 activities in two resistant populations were 1.97- and 2.17-fold higher than that in the susceptible population. Among 79 P450 genes, CYP6CX3 expressions in two resistant populations were 3.08- and 3.67-fold higher than that in the susceptible population. When CYP6CX3 was knocked down, the toxicity of cyantraniliprole increased significantly. The LC₅₀ value of cyantraniliprole to the Drosophila melanogaster line overexpressing B. tabaci CYP6CX3 increased 7.34-fold. The content of cyantraniliprole was decreased by 25.74 ± 4.27% after mixing with CYP6CX3 and CPR for 2 h. These results suggested that the overexpression of CYP6CX3 was likely involved in the resistance to cyantraniliprole in B. tabaci.

Differential Transcriptional Responses in Two Old World Bemisia tabaci Cryptic Species Post Acquisition of Old and New World Begomoviruses

Begomoviruses are transmitted by several cryptic species of the sweetpotato whitefly, Bemisia tabaci (Gennadius), in a persistent and circulative manner. Upon virus acquisition and circulative translocation within the whitefly, a multitude of molecular interactions occur. This study investigated the differentially expressed transcript profiles associated with the acquisition of the Old World monopartite begomovirus, tomato yellow leaf curl virus (TYLCV), and two New World bipartite begomoviruses, sida golden mosaic virus (SiGMV) and cucurbit leaf crumple virus (CuLCrV), in two invasive B. tabaci cryptic species, Middle East-Asia Minor 1 (MEAM1) and Mediterranean (MED). A total of 881 and 559 genes were differentially expressed in viruliferous MEAM1 and MED whiteflies, respectively, compared with their non-viruliferous counterparts, of which 146 genes were common between the two cryptic species. For both cryptic species, the number of differentially expressed genes (DEGs) associated with TYLCV and SiGMV acquisition were higher compared with DEGs associated with CuLCrV acquisition. Pathway analysis indicated that the acquisition of begomoviruses induced differential changes in pathways associated with metabolism and organismal systems. Contrasting expression patterns of major genes associated with virus infection and immune systems were observed. These genes were generally overexpressed and underexpressed in B. tabaci MEAM1 and MED adults, respectively. Further, no specific expression pattern was observed among genes associated with fitness (egg production, spermatogenesis, and aging) in viruliferous whiteflies. The weighted gene correlation network analysis of viruliferous B. tabaci MEAM1 and MED adults identified different hub genes potentially implicated in the vector competence and circulative tropism of viruses. Taken together, the results indicate that both vector cryptic species and the acquired virus species could differentially affect gene expression.

Morphology-based identification of Bemisia tabaci cryptic species puparia via embedded group-contrast convolution neural network analysis

The Bemisia tabaci species complex is a group of tropical–subtropical hemipterans, some species of which have achieved global distribution over the past 150 years. Several species are regarded currently as among the world’s most pernicious agricultural pests, causing a variety of damage types via direct feeding and plant-disease transmission. Long considered a single variable species, genetic, molecular and reproductive compatibility analyses have revealed that this “species” is actually a complex of between 24 and 48 morphologically cryptic species. However, determinations of which populations represent distinct species have been hampered by a failure to integrate genetic/molecular and morphological species–diagnoses. This, in turn, has limited the success of outbreak-control and eradication programs. Previous morphological investigations, based on traditional and geometric morphometric procedures, have had limited success in identifying genetic/molecular species from patterns of morphological variation in puparia. As an alternative, our investigation focused on exploring the use of a deep-learning convolution neural network (CNN) trained on puparial images and based on an embedded, group-contrast training protocol as a means of searching for consistent differences in puparial morphology. Fifteen molecular species were selected for analysis, all of which had been identified via DNA barcoding and confirmed using more extensive molecular characterizations and crossing experiments. Results demonstrate that all 15 species can be discriminated successfully based on differences in puparium morphology alone. This level of discrimination was achieved for laboratory populations reared on both hairy-leaved and glabrous-leaved host plants. Moreover, cross-tabulation tests confirmed the generality and stability of the CNN discriminant system trained on both ecophenotypic variants. The ability to identify B. tabaci species quickly and accurately from puparial images has the potential to address many long-standing problems in B. tabaci taxonomy and systematics as well as playing a vital role in ongoing pest-management efforts. [Aleyrodidae; entomology; Hemiptera; machine learning; morphometrics; pest control; systematics; taxonomy; whiteflies.]

Contribution of historical herbarium small RNAs to the reconstruction of a cassava mosaic geminivirus evolutionary history

Emerging viral diseases of plants are recognised as a growing threat to global food security. However, little is known about the evolutionary processes and ecological factors underlying the emergence and success of viruses that have caused past epidemics. With technological advances in the field of ancient genomics, it is now possible to sequence historical genomes to provide a better understanding of viral plant disease emergence and pathogen evolutionary history. In this context, herbarium specimens represent a valuable source of dated and preserved material. We report here the first historical genome of a crop pathogen DNA virus, a 90-year-old African cassava mosaic virus (ACMV), reconstructed from small RNA sequences bearing hallmarks of small interfering RNAs. Relative to tip-calibrated dating inferences using only modern data, those performed with the historical genome yielded both molecular evolution rate estimates that were significantly lower, and lineage divergence times that were significantly older. Crucially, divergence times estimated without the historical genome appeared in discordance with both historical disease reports and the existence of the historical genome itself. In conclusion, our study reports an updated time-frame for the history and evolution of ACMV and illustrates how the study of crop viral diseases could benefit from natural history collections.

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.

Is High Whitefly Abundance on Cassava in Sub-Saharan Africa Driven by Biological Traits of a Specific, Cryptic Bemisia tabaci Species?

In East Africa, the prevalent Bemisia tabaci whiteflies on the food security crop cassava are classified as sub-Saharan Africa (SSA) species. Economically damaging cassava whitefly populations were associated with the SSA2 species in the 1990s, but more recently, it has been to SSA1 species. To investigate whether biological traits (number of first instar nymphs, emerged adults, proportion of females in progeny and development time) of the cassava whitefly species are significant drivers of the observed field abundance, our study determined the development of SSA1 sub-group (SG) 1 (5 populations), SG2 (5 populations), SG3 (1 population) and SSA2 (1 population) on cassava and eggplant under laboratory conditions. SSA1-(SG1-SG2) and SSA2 populations' development traits were similar. Regardless of the host plant, SSA1-SG2 populations had the highest number of first instar nymphs (60.6 ± 3.4) and emerged adults (50.9 ± 3.6), followed by SSA1-SG1 (55.5 ± 3.2 and 44.6 ± 3.3), SSA2 (45.8 ± 5.7 and 32.6 ± 5.1) and the lowest were SSA1-SG3 (34.2 ± 6.1 and 32.0 ± 7.1) populations. SSA1-SG3 population had the shortest egg-adult emergence development time (26.7 days), followed by SSA1-SG1 (29.1 days), SSA1-SG2 (29.6 days) and SSA2 (32.2 days). Regardless of the whitefly population, development time was significantly shorter on eggplant (25.1 ± 0.9 days) than cassava (34.6 ± 1.0 days). These results support that SSA1-(SG1-SG2) and SSA2 B. tabaci can become highly abundant on cassava, with their species classification alone not correlating with observed abundance and prevalence.

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.