Genome of the African cassava whitefly Bemisia tabaci and distribution and genetic diversity of cassava-colonizing whiteflies in Africa

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.

Diversity and functional characteristics of culturable bacterial endosymbionts from cassava whitefly biotype Asia II-5, Bemisia tabaci

Whitefly Bemisia tabaci, a carrier of cassava mosaic disease (CMD), poses a significant threat to cassava crops. Investigating culturable bacteria and their impact on whiteflies is crucial due to their vital role in whitefly fitness and survival. The whitefly biotype associated with cassava and transmitting CMD in India has been identified as Asia II 5 through partial mitochondrial cytochrome oxidase I gene sequencing. In this study, bacteria associated with adult B. tabaci feeding on cassava were extracted using seven different media. Nutrient Agar (NA), Soyabean Casein Digest Medium (SCDM), Luria Bertani agar (LBA), and Reasoner's 2A agar (R2A) media resulted in 19, 6, 4, and 4 isolates, respectively, producing a total of 33 distinct bacterial isolates. Species identification through 16SrRNA gene sequencing revealed that all isolates belonged to the Bacillota and Pseudomonadota phyla, encompassing 11 genera: Bacillus, Cytobacillus, Exiguobacterium, Terribacillus, Brevibacillus, Enterococcus, Staphylococcus, Brucella, Novosphingobium, Lysobacter, and Pseudomonas. All bacterial isolates were tested for chitinase, protease, siderophore activity, and antibiotic sensitivity. Nine isolates exhibited chitinase activity, 28 showed protease activity, and 23 displayed siderophore activity. Most isolates were sensitive to antibiotics such as Vancomycin, Streptomycin, Erythromycin, Kanamycin, Doxycycline, Tetracycline, and Ciprofloxacin, while they demonstrated resistance to Bacitracin and Colistin. Understanding the culturable bacteria associated with cassava whitefly and their functional significance could contribute to developing effective cassava whitefly and CMD control in agriculture.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03949-0.

In silico prediction of candidate gene targets for the management of African cassava whitefly (Bemisia tabaci, SSA1-SG1), a key vector of viruses causing cassava brown streak disease

Whiteflies (Bemisia tabaci sensu lato) have a wide host range and are globally important agricultural pests. In Sub-Saharan Africa, they vector viruses that cause two ongoing disease epidemics: cassava brown streak disease and cassava mosaic virus disease. These two diseases threaten food security for more than 800 million people in Sub-Saharan Africa. Efforts are ongoing to identify target genes for the development of novel management options against the whitefly populations that vector these devastating viral diseases affecting cassava production in Sub-Saharan Africa. This study aimed to identify genes that mediate osmoregulation and symbiosis functions within cassava whitefly gut and bacteriocytes and evaluate their potential as key gene targets for novel whitefly control strategies. The gene expression profiles of dissected guts, bacteriocytes and whole bodies were compared by RNAseq analysis to identify genes with significantly enriched expression in the gut and bacteriocytes. Phylogenetic analyses identified three candidate osmoregulation gene targets: two α-glucosidases, SUC 1 and SUC 2 with predicted function in sugar transformations that reduce osmotic pressure in the gut; and a water-specific aquaporin (AQP1) mediating water cycling from the distal to the proximal end of the gut. Expression of the genes in the gut was enriched 23.67-, 26.54- and 22.30-fold, respectively. Genome-wide metabolic reconstruction coupled with constraint-based modeling revealed four genes (argH, lysA, BCAT & dapB) within the bacteriocytes as potential targets for the management of cassava whiteflies. These genes were selected based on their role and essentiality within the different essential amino acid biosynthesis pathways. A demonstration of candidate osmoregulation and symbiosis gene targets in other species of the Bemisia tabaci species complex that are orthologs of the empirically validated osmoregulation genes highlights the latter as promising gene targets for the control of cassava whitefly pests by in planta RNA interference.

Gene identification and RNAi-silencing of p62/SQSTM1 in the vector Rhodnius prolixus reveals a high degree of sequence conservation but no apparent deficiency-related phenotypes in vitellogenic females

Autophagy and the ubiquitin-proteasome system (UPS) are important cellular mechanisms that coordinate protein degradation essential for proteostasis. P62/SQSTM1 is a receptor cargo protein able to deliver ubiquitinated targets to the proteasome proteolytic complex and/or to the autophagosome. In the insect vector of Chagas disease, Rhodnius prolixus, previous works have shown that the knockdown of different autophagy-related genes (ATGs) and ubiquitin-conjugating enzymes resulted in abnormal oogenesis phenotypes and embryo lethality. Here, we investigate the role of the autophagy/UPS adaptor protein p62 during the oogenesis and reproduction of this vector. We found that R. prolixus presents one isoform of p62 encoded by a non-annotated gene. The predicted protein presents the domain architecture anticipated for p62: PB1 (N-term), ZZ-finger, and UBA (C-term) domains, and phylogenetic analysis showed that this pattern is highly conserved within insects. Using parental RNAi, we found that although p62 is expressed in the ovary, midgut, and fat body of adult females, systemic silencing of this gene did not result in any apparent phenotypes under in-house conditions. The insects' overall levels of blood meal digestion, lifespan, yolk protein production, oviposition, and embryo viability were not altered when compared to controls. Because it is known that autophagy and UPS can undergo compensatory mechanisms, we asked whether the silencing of p62 was triggering adaptative changes in the expression of genes of the autophagy, UPS, and the unfolded protein response (UPR) and found that only ATG1 was slightly up regulated in the ovaries of silenced females. In addition, experiments to further investigate the role of p62 in insects previously silenced for the E1-conjugating enzyme (a condition known to trigger the upregulation of p62), also did not result in any apparent phenotypes in vitellogenic females.

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.

Knockdown of carnitine palmitoyltransferase I (CPT1) reduces fat body lipid mobilization and resistance to starvation in the insect vector Rhodnius prolixus

The energy stored in fatty acids is essential for several critical activities of insects, such as embryogenesis, oviposition, and flight. Rhodnius prolixus is an obligatory hematophagous hemipteran and vector of Chagas disease, and it feeds infrequently on very large blood meals. As digestion slowly occurs, lipids are synthesized and accumulate in the fat body, mainly as triacylglycerol, in lipid droplets. Between feeding bouts, proper mobilization and oxidation of stored lipids are crucial for survival, and released fatty acids are oxidized by mitochondrial β-oxidation. Carnitine palmitoyl transferase I (CPT1) is the enzyme that catalyzes the first reaction of the carnitine shuttle, where the activated fatty acid, acyl-CoA, is converted to acyl-carnitine to be transported into the mitochondria. Here, we investigated the role of CPT1 in lipid metabolism and in resistance to starvation in Rhodnius prolixus. The expression of the CPT1 gene (RhoprCpt1) was determined in the organs of adult females on the fourth day after a blood meal, and the flight muscle showed higher expression levels than the ovary, fat body, and anterior and posterior midgut. RhoprCpt1 expression in the fat body dramatically decreased after feeding, and started to increase again 10 days later, but no changes were observed in the flight muscle. β-oxidation rates were determined in flight muscle and fat body homogenates with the use of ³H-palmitate, and in unfed females, they were higher in the flight muscle. In the fat body, lipid oxidation activity did not show any variation before or at different days after feeding, and was not affected by the presence of etomoxir or malonyl-CoA. We used RNAi and generated RhoprCPT1-deficient insects, which surprisingly did not show a decrease in measured ³H-palmitate oxidation rates. However, the RNAi-knockdown females presented increased amounts of triacylglycerol and larger lipid droplets in the fat body, but not in the flight muscle. When subjected to starvation, these insects had a shorter lifespan. These results indicated that the inhibition of RhoprCpt1 expression compromised lipid mobilization and affected resistance to starvation.

Silencing of the 20S proteasomal subunit-α6 triggers full oogenesis arrest and increased mRNA levels of the selective autophagy adaptor protein p62/SQSTM1 in the ovary of the vector Rhodnius prolixus

The high reproductive rates of insects contribute significantly to their ability to act as vectors of a variety of vector-borne diseases. Therefore, it is strategically critical to find molecular targets with biotechnological potential through the functional study of genes essential for insect reproduction. The ubiquitin-proteasome system is a vital degradative pathway that contributes to the maintenance of regular eukaryotic cell proteostasis. This mechanism involves the action of enzymes to covalently link ubiquitin to proteins that are meant to be delivered to the 26S proteasome and broken down. The 26S proteasome is a large protease complex (including the 20S and 19S subcomplexes) that binds, deubiquitylates, unfolds, and degrades its substrates. Here, we used bioinformatics to identify the genes that encode the seven α and β subunits of the 20S proteasome in the genome of R. prolixus and learned that those transcripts are accumulated into mature oocytes. To access proteasome function during oogenesis, we conducted RNAi functional tests employing one of the 20S proteasome subunits (Prosα6) as a tool to suppress 20S proteasomal activity. We found that Prosα6 silencing resulted in no changes in TAG buildup in the fat body and unaffected availability of yolk proteins in the hemolymph of vitellogenic females. Despite this, the silencing of Prosα6 culminated in the impairment of oocyte maturation at the early stages of oogenesis. Overall, we discovered that proteasome activity is especially important for the signals that initiate oogenesis in R. prolixus and discuss in what manner further investigations on the regulation of proteasome assembly and activity might contribute to the unraveling of oogenesis molecular mechanisms and oocyte maturation in this vector.

ATP synthase affects lipid metabolism in the kissing bug Rhodnius prolixus beyond its role in energy metabolism

ATP synthase plays an essential role in mitochondrial metabolism, being responsible for the production of ATP in oxidative phosphorylation. However, recent results have shown that it may also be present in the cell membrane, involved in lipophorin binding to its receptors. Here, we used a functional genetics approach to investigate the roles of ATP synthase in lipid metabolism in the kissing bug Rhodnius prolixus. The genome of R. prolixus encodes five nucleotide-binding domain genes of the ATP synthase alpha and beta family, including the alpha and beta subunits of ATP synthase (RpATPSynA and RpATPSynB), and the catalytic and non-catalytic subunits of the vacuolar ATPase (RpVha68 and RpVha55). These genes were expressed in all analyzed organs, being their expression highest in the ovaries, fat body and flight muscle. Feeding did not regulate the expression of ATP synthases in the posterior midgut or fat body. Furthermore, ATP synthase is present in the fat body's mitochondrial and membrane fractions. RpATPSynB knockdown by RNAi impaired ovarian development and reduced egg-laying by approximately 85%. Furthermore, the lack of RpATPSynB increased the amount of triacylglycerol in the fat body due to increased de novo fatty acid synthesis and reduced transfer of lipids to lipophorin. RpATPSynA knockdown had similar effects, with altered ovarian development, reduced oviposition, and triacylglycerol accumulation in the fat body. However, ATP synthases knockdown had only a slight effect on the amount of ATP in the fat body. These results support the hypothesis that ATP synthase has a direct role in lipid metabolism and lipophorin physiology, which are not directly due to changes in energy metabolism.

Sidestepping Darwin: horizontal gene transfer from plants to insects

Horizontal transfer of genetic material (HT) is the passage of DNA between organisms by means other than reproduction. Increasing numbers of HT are reported in insects, with bacteria, fungi, plants, and insects acting as the main sources of these transfers. Here, we provide a detailed account of plant-to-insect HT events. At least 14 insect species belonging to 6 orders are known to have received plant genetic material through HT. One of them, the whitefly Bemisia tabaci (Middle East Asia Minor 1), concentrates most of these transfers, with no less than 28 HT events yielding 55 plant-derived genes in this species. Several plant-to-insect HT events reported so far involve gene families known to play a role in plant-parasite interactions. We highlight methodological approaches that may further help characterize these transfers. We argue that plant-to-insect HT is likely more frequent than currently appreciated and that in-depth studies of these transfers will shed new light on plant-insect interactions.

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.

Survey of cassava brown streak disease and association of factors influencing its epidemics in smallholder cassava cropping systems of coastal Kenya

Cassava productivity is threatened by viral diseases which have become the main phytosanitary problems in cassava farmers. Cassava brown streak disease (CBSD) is a devastating viral disease caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) which are transmitted by whitefly vectors and mainly disseminated through the use of infected stem cuttings as planting materials. The aim of the current study was to: (1) assess farmers' knowledge, perceptions on spread, causes and current management practices of CBSD; (2) determine the factors associated with farmers' satisfaction with cassava planting material; and (3) determine the distribution, incidence, and severity of CBSD and association of factors influencing the disease epidemics in smallholder cassava cropping systems in coastal Kenya. Information was collected using semi-structured questionnaire administered to 250 smallholder farmers through face-to-face interviews coupled with field visits to assess the incidence, severity and distribution of CBSD. Symptomatic and asymptomatic cassava leaf samples were collected for reverse transcription-polymerase chain reaction (RT-PCR) analysis of the causal viruses of CBSD. The results revealed that majority of the farmers (96.6%) could recognize CBSD symptoms on the roots, and only 11.5% could recognize the foliar symptoms of the disease. The cause of the disease was unknown to the farmers, with no effective management methods available to them. Majority of farmers (82.5%) recycled own cassava cuttings from previous season's crop as planting material followed by exchanging/borrowing from neighbors (67.5%). The field incidence of CBSD was highest in Kilifi (27.9%) followed by Kwale (24.7%) and Taita Taveta (10.8%), with severities ranging from 2 to 3 in the three Counties. RT-PCR analysis indicated that 91% of the symptomatic samples tested positive for either of the two viruses occurring either singly or as dual infection. Approximately 3.2% of the asymptomatic samples tested positive for only CBSV. Findings from this study demonstrates the need for awareness creation of farmers on the causes, spread and management practices to control CBSD and the importance of strengthening certified cassava seed systems to reduce the impact of the disease. The study provides base-line information imperative for development of management strategies of CBSD.

Epidemiological Analysis of Cassava Mosaic and Brown Streak Diseases, and Bemisia tabaci in the Comoros Islands

A comprehensive assessment of cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) was carried out in Comoros where cassava yield (5.7 t/ha) is significantly below the African average (8.6 t/ha) largely due to virus diseases. Observations from 66 sites across the Comoros Islands of Mwali, Ngazidja, and Ndzwani revealed that 83.3% of cassava fields had foliar symptoms of CBSD compared with 95.5% for CMD. Molecular diagnostics confirmed the presence of both cassava brown streak ipomoviruses (CBSIs) and cassava mosaic begomoviruses (CMBs). Although real-time RT-PCR only detected the presence of one CBSI species (Cassava brown streak virus, CBSV) the second species (Ugandan cassava brown streak virus, UCBSV) was identified using next-generation high-throughput sequencing. Both PCR and HTS detected the presence of East African cassava mosaic virus (EACMV). African cassava mosaic virus was not detected in any of the samples. Four whitefly species were identified from a sample of 131 specimens: Bemisia tabaci, B. afer, Aleurodicus dispersus, and Paraleyrodes bondari. Cassava B. tabaci comprised two mitotypes: SSA1-SG2 (89%) and SSA1-SG3 (11%). KASP SNP genotyping categorized 82% of cassava B. tabaci as haplogroup SSA-ESA. This knowledge will provide an important base for developing and deploying effective management strategies for cassava viruses and their vectors.

Diversity and Distribution of Whiteflies Colonizing Cassava in Eastern Democratic Republic of Congo

The present study characterizes Bemisia tabaci and Bemisia afer from cassava in eastern Democratic Republic of Congo (DRC). The Mitochondrial COI sequencing revealed the occurrence of six cassava B. tabaci mitotypes, which were designated into four haplogroups (SSA-ECA, SSA-CA, SSA2, and SSA-ESA) using KASP SNP genotyping. SSA-ECA (72%) was the most prevalent and occurred in the northern part of the surveyed area, in the Ituri and Nord/Sud-Kivu provinces, whilst SSA-CA (21%) was present in the south, primarily in Haut-Katanga. SSA-ECA was predominant in the areas of north-eastern DRC most severely affected by cassava brown streak disease and was also reported in the new outbreak area in Pweto territory, Haut-Katanga, in the south. Bemisia afer comprised two major clusters with 85.5% of samples in cluster one, while the rest were in cluster two, which has no reference sequence in GenBank. This study provides important information on the genetic diversity of B. tabaci and B. afer in eastern DRC. This knowledge will be used as a basis for further studies to understand and to identify the role of whitefly haplogroups, their population densities and consequences for virus epidemics and spread as well as leading to improved vector and virus management strategies.

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.

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

In sub-Saharan Africa cassava growing areas, two members of the Bemisia tabaci species complex termed sub-Saharan Africa 1 (SSA1) and SSA2 have been reported as the prevalent whiteflies associated with the spread of viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) pandemics. At the peak of CMD pandemic in the late 1990s, SSA2 was the prevalent whitefly, although its numbers have diminished over the last two decades with the resurgence of SSA1 whiteflies. Three SSA1 subgroups (SG1 to SG3) are the predominant whiteflies in East Africa and vary in distribution and biological properties. Mating compatibility between SSA1 subgroups and SSA2 whiteflies was reported as the possible driver for the resurgence of SSA1 whiteflies. In this study, a combination of both phylogenomic methods and reciprocal crossing experiments were applied to determine species status of SSA1 subgroups and SSA2 whitefly populations. Phylogenomic analyses conducted with 26 548 205 bp whole genome single nucleotide polymorphisms (SNPs) and the full mitogenomes clustered SSA1 subgroups together and separate from SSA2 species. Mating incompatibility between SSA1 subgroups and SSA2 further demonstrated their distinctiveness from each other. Phylogenomic analyses conducted with SNPs and mitogenomes also revealed different genetic relationships among SSA1 subgroups. The former clustered SSA1-SG1 and SSA1-SG2 together but separate from SSA1-SG3, while the latter clustered SSA1-SG2 and SSA1-SG3 together but separate from SSA1-SG1. Mating compatibility was observed between SSA1-SG1 and SSA1-SG2, while incompatibility occurred between SSA1-SG1 and SSA1-SG3, and SSA1-SG2 and SSA1-SG3. Mating results among SSA1 subgroups were coherent with phylogenomics results based on SNPs but not the full mitogenomes. Furthermore, this study revealed that the secondary endosymbiont-Wolbachia-did not mediate reproductive success in the crossing assays carried out. Overall, using genome wide SNPs together with reciprocal crossings assays, this study established accurate genetic relationships among cassava-colonizing populations, illustrating that SSA1 and SSA2 are distinct species while at least two species occur within SSA1 species.

Whitefly endosymbionts: IPM opportunity or tilting at windmills?

Whiteflies are sap-sucking insects responsible for high economic losses. They colonize hundreds of plant species and cause direct feeding damage and indirect damage through transmission of devastating viruses. Modern agriculture has seen a history of invasive whitefly species and populations that expand to novel regions, bringing along fierce viruses. Control efforts are hindered by fast virus transmission, insecticide-resistant populations, and a wide host range which permits large natural reservoirs for whiteflies. Augmentative biocontrol by parasitoids while effective in suppressing high population densities in greenhouses falls short when it comes to preventing virus transmission and is ineffective in the open field. A potential source of much needed novel control strategies lays within a diverse community of whitefly endosymbionts. The idea to exploit endosymbionts for whitefly control is as old as identification of these bacteria, yet it still has not come to fruition. We review where our knowledge stands on the aspects of whitefly endosymbiont evolution, biology, metabolism, multitrophic interactions, and population dynamics. We show how these insights are bringing us closer to the goal of better integrated pest management strategies. Combining most up to date understanding of whitefly-endosymbiont interactions and recent technological advances, we discuss possibilities of disrupting and manipulating whitefly endosymbionts, as well as using them for pest control.

Biotechnological interventions for the sustainable management of a global pest, whitefly (Bemisia tabaci)

Whiteflies (Bemisia tabaci) are polyphagous invasive hemipteran insects that cause serious losses of important crops by directly feeding on phloem sap and transmitting pathogenic viruses. These insects have emerged as a major threat to global agriculture and food security. Chemically synthesized insecticides are currently the only option to control whiteflies, but the ability of whiteflies to evolve resistance against insecticides has made the management of these insects very difficult. Natural host-plant resistance against whiteflies identified in some crop plants has not been exploited to a great extent. Genetic engineering approaches, such as transgenics and RNA interference (RNAi), are potentially useful for the control of whiteflies. Transgenic plants harboring insecticidal toxins/lectins developed via nuclear or chloroplast transformation are a promising vehicle for whitefly control. Double-stranded RNAs (dsRNAs) of several insect genes, delivered either through microinjection into the insect body cavity or orally via an artificial diet and transiently or stably expressed in transgenic plants, have controlled whiteflies in model plants and in some crops at the laboratory level, but not at the field level. In this review, we highlight the merits and demerits of each delivery method along with strategies for sustained delivery of dsRNAs via fungal entomopathogen/endosymbiont or nontransgenic RNAi approaches, foliar sprays, root absorption or nanocarriers as well as the factors affecting efficient RNAi and their biosafety issues. Genome sequencing and transcriptome studies of whitefly species are facilitating the selection of appropriate genes for RNAi and gene-editing technology for the efficient and resilient management of whiteflies and their transmitted viruses.

Distribution and Molecular Diversity of Whitefly Species Colonizing Cassava in Kenya

Simple Summary

The whitefly, Bemisia tabaci (Gennadium, Hemiptera) is a crop pest and plant-virus vector known to transmit more than 300 plant viruses. Among other plant viruses, the vector transmits viruses that cause the two major devastating viral diseases of cassava in sub-Saharan Africa namely cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). In order to achieve effective implementation of whitefly management programs including prevention of spread of the species, accurate species identification is vital. Morphological identification approaches toward whitefly species limits the capacity to accurately identify new species, especially the presence of cryptic species such as the numerous B. tabaci genetic groups. Using the mitochondrial DNA cytochrome oxidase 1 (mtCO1) sequences, four distinct whitefly species namely Bemisia tabaci, Bemisia afer, Aleurodicus dispersus and Paraleyrodes bondari were identified in samples collected from major cassava growing regions in Kenya. The study presents the first report of P. bondari (Bondar’s nesting whitefly) on cassava in Kenya. We found three B. tabaci genetic groups of SSA1, SSA2 and Indian Ocean (IO) putative species colonizing cassava in Kenya. The information is useful to inform better management strategies of the whitefly vectors to reduce the impact of cassava viral diseases, which continue to be a threat to food security in major cassava growing regions.

Abstract

The whitefly, Bemisia tabaci (Gennadium, Hemiptera) has been reported to transmit viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) in many parts of sub-Saharan Africa (SSA). Currently, there is limited information on the distribution, species and haplotype composition of the whitefly populations colonizing cassava in Kenya. A study was conducted in the major cassava growing regions of Kenya to address this gap. Analyses of mitochondrial DNA cytochrome oxidase 1 (mtCO1) sequences revealed the presence of four distinct whitefly species: Bemisia tabaci, Bemisia afer, Aleurodicus dispersus and Paraleyrodes bondari in Kenya. The B. tabaci haplotypes were further resolved into SSA1, SSA2 and Indian Ocean (IO) putative species. The SSA1 population had three haplogroups of SSA1-SG1, SSA-SG2 and SSA1-SG3. Application of KASP genotyping grouped the Bemisia tabaci into two haplogroups namely sub-Saharan Africa East and Southern Africa (SSA-ESA) and sub-Saharan Africa East and Central Africa (SSA-ECA). The study presents the first report of P. bondari (Bondar’s nesting whitefly) on cassava in Kenya. Bemisia tabaci was widely distributed in all the major cassava growing regions in Kenya. The increased detection of different whitefly species on cassava and genetically diverse B. tabaci mitotypes indicates a significant influence on the dynamics of cassava virus epidemics in the field. The study highlights the need for continuous monitoring of invasive whitefly species population on cassava for timely application of management practices to reduce the impact of cassava viral diseases and prevent potential yield losses.

Why has permanent control of cassava brown streak disease in Sub-Saharan Africa remained a dream since the 1930s?

Effective control of ipomoviruses that cause cassava brown streak disease (CBSD) in Africa has remained problematic despite eight remarkable decades (1930-2021) of research efforts. Molecular mechanisms underlying resistance breakdown in genetically improved cassava are still unknown. The vast genetic diversity of cassava brown streak viruses, which is crucial for the improvement of routine reverse transcription polymerase chain reaction (RT-qPCR) assays in CBSD-endemic regions of Africa, is controversial and underrepresented. From a molecular epidemiology viewpoint, this review discusses the reasons for why permanent control of CBSD is difficult in the modern era, even with the presence of diverse in silico and omics tools, recombinant DNA, and high throughput next-generation sequencing technologies. Following an extensive nucleotide data search in the National Centre for Biotechnology Information (NCBI) database and a literature review in PubMed and Scopus, we report that genomic data of 87.62% (474/541) strains of cassava brown streak virus are missing due to poor sequencing capacity in Africa. The evolution dynamics of viral virulence and pathogenicity has not yet been fully explored from the available 67 (12.38%) genomic sequences, owing to poor bioinformatics capacity. Tanzania and Zambia have the highest and lowest disease inoculum pressure, correspondingly. Knowledge gaps in molecular biology and the overall molecular pathogenesis of CBSD viruses impede effective disease control in Africa. Recommendations for possible solutions to the research questions, controversies, and hypotheses raised in this study serve as a roadmap for the invention of more effective CBSD control methods.

Two Deoxythymidine Triphosphate Synthesis-Related Genes Regulate Obligate Symbiont Density and Reproduction in the Whitefly Bemisia tabaci MED

Deoxythymidine triphosphate (dTTP) is essential for DNA synthesis and cellular growth in all organisms. Here, genetic capacity analysis of the pyrimidine pathway in insects and their symbionts revealed that dTTP is a kind of metabolic input in several host insect/obligate symbiont symbiosis systems, including Bemisia tabaci MED/Candidatus Portiera aleyrodidarum (hereafter Portiera). As such, the roles of dTTP on both sides of the symbiosis system were investigated in B. tabaci MED/Portiera. Dietary RNA interference (RNAi) showed that suppressing dTTP production significantly reduced the density of Portiera, significantly repressed the expression levels of horizontally transferred essential amino acid (EAA) synthesis-related genes, and significantly decreased the reproduction of B. tabaci MED adults as well as the hatchability of their offspring. Our results revealed the regulatory role of dTTP in B. tabaci MED/Portiera and showed that dTTP synthesis-related genes could be potential targets for controlling B. tabaci as well as other sucking pests.

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.

What Is the Spatial Extent of a Bemisia tabaci Population?

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

Genetic and comparative mapping of Lupinus luteus L. highlight syntenic regions with major orthologous genes controlling anthracnose resistance and flowering time

Anthracnose susceptibility and ill-adapted flowering time severely affect Lupinus luteus yield, which has high seed protein content, is excellent for sustainable agriculture, but requires genetic improvement to fulfil its potential. This study aimed to (1) develop a genetic map; (2) define collinearity and regions of synteny with Lupinus angustifolius; and (3) map QTLs/candidate genes for anthracnose resistant and flowering time. A few linkage groups/genomic regions tended to be associated with segregation distortion, but did not affect the map. The developed map showed collinearity, and syntenic regions with L. angustifolius. Major QTLs were mapped in syntenic regions. Alleles from the wild parent and cultivar, explained 75% of the phenotypic variance for anthracnose resistance and 83% for early flowering, respectively. Marker sequences flanking the QTLs showed high homology with the Lanr1 gene and Flowering-locus-T of L. angustifolius. This suggests orthologous genes for both traits in the L. luteus genome. The findings are remarkable, revealing the potential to combine early flowering/anthracnose resistant in fulfilling yield capacity in L. luteus, and can be a major strategy in the genetic improvement and usage of this species for sustainable protein production. Allele sequences and PCR-marker tagging of these genes are being applied in marker assisted selection.

Genetic Diversity of Mitochondrial DNA of Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) Associated with Cassava and the Occurrence of Cassava Mosaic Disease in Zambia

Bemisia tabaci is an important vector of cassava brown streak viruses and cassava mosaic begomoviruses, the causal agents of cassava brown streak disease and cassava mosaic disease (CMD), respectively. A study was carried out to determine the genetic variability of B. tabaci associated with cassava and the occurrence of CMD in Zambia in 2013 and 2015. Phylogenetic analysis showed the presence of only the sub-Saharan Africa 1 (SSA1) genetic group in Zambia. The SSA1 population had three population subgroups (SGs): SSA1-SG1, SSA1-SG2 and SSA1-SG3. All three SSA1 population subgroups occurred in Western Province. However, only SSA1-SG3 occurred in Eastern Province, while only SSA1-SG1 occurred in North Western and Luapula Provinces. Adult B. tabaci were most abundant in Western Province in 2013 (11.1/plant) and 2015 (10.8/plant), and least abundant (0.2/plant) in Northern Province in both 2013 and 2015. CMD was prevalent in all seven provinces surveyed, with the highest incidence recorded in Lusaka Province in both 2013 (78%) and 2015 (83.6%), and the lowest in Northern Province in both 2013 (26.6%) and 2015 (29.3%). Although SSA1-SG1 occurred at greater abundances than the other subgroups, there was no direct association demonstrated between whitefly subgroup and incidence of CMD. Establishing which B. tabaci genetic groups and populations are associated with CMD and their distribution in the country is a key factor in guiding the development of CMD control strategies for cassava-dependent households.

Debugging: Strategies and Considerations for Efficient RNAi-Mediated Control of the Whitefly Bemisia tabaci

The whitefly Bemisia tabaci is a globally important pest that is difficult to control through insecticides, transgenic crops, and natural enemies. Post-transcriptional gene silencing through RNA interference (RNAi) has shown potential as a pest management strategy against B. tabaci. While genomic data and other resources are available to create highly effective customizable pest management strategies with RNAi, current applications do not capitalize on species-specific biology. This lack of specificity has the potential to have substantial ecological impacts. Here, we discuss both short- and long-term considerations for sustainable RNAi pest management strategies for B. tabaci, focusing on the need for species specificity incorporating both life history and population genetic considerations. We provide a conceptual framework for selecting sublethal target genes based on their involvement in physiological pathways, which has the greatest potential to ameliorate unintended negative consequences. We suggest that these considerations allow an integrated pest management approach, with fewer negative ecological impacts and reduced likelihood of the evolution of resistant populations.

Whitefly genomes contain ribotoxin coding genes acquired from plants

Ribosome inactivating proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine residue in the conserved sarcin/ricin loop of 28S rRNA. These enzymes are widely distributed among plants and bacteria. Previously, we have described for the first time RIP genes in mosquitoes belonging to the Culicidae family. We showed that these genes are derived from a single event of horizontal gene transfer (HGT) from a prokaryotic donor. Mosquito RIP genes are evolving under purifying selection, strongly suggesting that these toxins have acquired a functional role. In this work, we show the existence of two RIP encoding genes in the genome of the whitefly Bemisia tabaci, a hemiptera species belonging to the Aleyrodidae family distantly related to mosquitoes. Contamination artifacts were ruled out analyzing three independent B. tabaci genome databases. In contrast to mosquito RIPs, whitefly genes harbor introns and according to transcriptomic evidence are transcribed and spliced. Phylogeny and the taxonomic distribution strongly support that whitefly RIP genes are derived from an independent HGT event from a plant source. These results, along with our previous description of RIPs in Diptera, suggest that the acquired genes are functional in these insects and confer some fitness advantage.

KASP Genotyping as a Molecular Tool for Diagnosis of Cassava-Colonizing Bemisia tabaci

Bemisia tabaci is a cryptic species complex that requires the use of molecular tools for identification. The most widely used approach for achieving this is the partial sequencing of the mitochondrial DNA cytochrome oxidase I gene (COI). A more reliable single nucleotide polymorphism (SNP)-based genotyping approach, using Nextera restriction-site-associated DNA (NextRAD) sequencing, has demonstrated the existence of six major haplogroups of B. tabaci on cassava in Africa. However, NextRAD sequencing is costly and time-consuming. We, therefore, developed a cheaper and more rapid diagnostic using the Kompetitive Allele-Specific PCR (KASP) approach. Seven sets of primers were designed to distinguish the six B. tabaci haplogroups based on the NextRAD data. Out of the 152 whitefly samples that were tested using these primer sets, 151 (99.3%) produced genotyping results consistent with NextRAD. The KASP assay was designed using NextRAD data on whiteflies from cassava in 18 countries across sub-Saharan Africa. This assay can, therefore, be routinely used to rapidly diagnose cassava B. tabaci by laboratories that are researching or monitoring this pest in Africa. This is the first study to develop an SNP-based assay to distinguish B. tabaci whiteflies on cassava in Africa, and the first application of the KASP technique for insect identification.

Cassava whitefly species in eastern Nigeria and the threat of vector-borne pandemics from East and Central Africa

Bemisia tabaci (sensu latu) is a group of >40 highly cryptic whitefly species that are of global agricultural importance, both as crop pests and plant-virus vectors. Two devastating cassava diseases in East and Central Africa are spread by abundant populations of one of these species termed Sub-Saharan Africa 1 (SSA1). There is a substantive risk that these whitefly-borne pandemics will continue to spread westwards and disrupt cassava production for millions of smallholder farmers in West Africa. We report here, therefore, the first comprehensive survey of cassava B. tabaci in eastern Nigeria, a West African region likely to be the first affected by the arrival of these whitefly-borne pandemics. We found one haplotype comprising 32 individuals with 100% identical mtCO1 sequence to the East African SSA1 populations (previously termed SSA1-SG1) and 19 mtCO1 haplotypes of Sub-Saharan Africa 3 (SSA3), the latter being the most prevalent and widely distributed B. tabaci species in eastern Nigeria. A more divergent SSA1 mtCO1 sequence (previously termed SSA1-SG5) was also identified in the region, as were mtCO1 sequences identifying the presence of the MED ASL B. tabaci species and Bemisia afer. Although B. tabaci SSA1 was found in eastern Nigeria, they were not present in the high abundances associated with the cassava mosaic (CMD) and cassava brown streak disease (CBSD) pandemics of East and Central Africa. Also, no severe CMD or any CBSD symptoms were found in the region.

Genetic Diversity of Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) Colonizing Sweet Potato and Cassava in South Sudan

Bemisia tabaci (Gennadius) is a polyphagous, highly destructive pest that is capable of vectoring viruses in most agricultural crops. Currently, information regarding the distribution and genetic diversity of B. tabaci in South Sudan is not available. The objectives of this study were to investigate the genetic variability of B. tabaci infesting sweet potato and cassava in South Sudan. Field surveys were conducted between August 2017 and July and August 2018 in 10 locations in Juba County, Central Equatoria State, South Sudan. The sequences of mitochondrial DNA cytochrome oxidase I (mtCOI) were used to determine the phylogenetic relationships between sampled B. tabaci. Six distinct genetic groups of B. tabaci were identified, including three non-cassava haplotypes (Mediterranean (MED), Indian Ocean (IO), and Uganda) and three cassava haplotypes (Sub-Saharan Africa 1 sub-group 1 (SSA1-SG1), SSA1-SG3, and SSA2). MED predominated on sweet potato and SSA2 on cassava in all of the sampled locations. The Uganda haplotype was also widespread, occurring in five of the sampled locations. This study provides important information on the diversity of B. tabaci species in South Sudan. A comprehensive assessment of the genetic diversity, geographical distribution, population dynamics, and host range of B. tabaci species in South Sudan is vital for its effective management.

Nuclear Orthologs Derived from Whole Genome Sequencing Indicate Cryptic Diversity in the Bemisia tabaci (Insecta: Aleyrodidae) Complex of Whiteflies

The Bemisia tabaci complex of whiteflies contains globally important pests thought to contain cryptic species corresponding to geographically structured phylogenetic clades. Although mostly morphologically indistinguishable, differences have been shown to exist among populations in behavior, plant virus vector capacity, ability to hybridize, and DNA sequence divergence. These differences allow for certain populations to become invasive and cause great economic damage in a monoculture setting. Although high mitochondrial DNA divergences have been reported between putative conspecifics of the B. tabaci species complex, there is limited data that exists across the whole genome for this group. Using data from 2184 orthologs obtained from whole genome sequencing (Illumina), a phylogenetic analysis using maximum likelihood and coalescent methodologies was completed on ten individuals of the B. tabaci complex. In addition, automatic barcode gap discovery methods were employed, and results suggest the existence of five species. Although the divergences of the mitochondrial cytochrome oxidase I gene are high among members of this complex, nuclear divergences are much lower in comparison. Single-copy orthologs from whole genome sequencing demonstrate divergent population structures among members of the B. tabaci complex and the sequences provide an important resource to aid in future genomic studies of the group.