Spatio-temporal patterns of genetic change amongst populations of cassava Bemisia tabaci whiteflies driving virus pandemics in East and Central Africa

Functional role of geminivirus encoded proteins in the host: Past and present

During plant-pathogen interaction, plant exhibits a strong defense system utilizing diverse groups of proteins to suppress the infection and subsequent establishment of the pathogen. However, in response, pathogens trigger an anti-silencing mechanism to overcome the host defense machinery. Among plant viruses, geminiviruses are the second largest virus family with a worldwide distribution and continue to be production constraints to food, feed, and fiber crops. These viruses are spread by a diverse group of insects, predominantly by whiteflies, and are characterized by a single-stranded DNA (ssDNA) genome coding for four to eight proteins that facilitate viral infection. The most effective means to managing these viruses is through an integrated disease management strategy that includes virus-resistant cultivars, vector management, and cultural practices. Dynamic changes in this virus family enable the species to manipulate their genome organization to respond to external changes in the environment. Therefore, the evolutionary nature of geminiviruses leads to new and novel approaches for developing virus-resistant cultivars and it is essential to study molecular ecology and evolution of geminiviruses. This review summarizes the multifunctionality of each geminivirus-encoded protein. These protein-based interactions trigger the abrupt changes in the host methyl cycle and signaling pathways that turn over protein normal production and impair the plant antiviral defense system. Studying these geminivirus interactions localized at cytoplasm-nucleus could reveal a more clear picture of host-pathogen relation. Data collected from this antagonistic relationship among geminivirus, vector, and its host, will provide extensive knowledge on their virulence mode and diversity with climate change.

Factor-Dependent Internal Ribosome Entry Site and -1 Programmed Frameshifting Signal in the Bemisia-Associated Dicistrovirus 2

The dicistrovirus intergenic (IGR) IRES uses the most streamlined translation initiation mechanism: the IRES recruits ribosomes directly without using protein factors and initiates translation from a non-AUG codon. Several subtypes of dicistroviruses IRES have been identified; typically, the IRESs adopt two -to three overlapping pseudoknots with key stem-loop and unpaired regions that interact with specific domains of the ribosomal 40S and 60S subunits to direct translation. We previously predicted an atypical IGR IRES structure and a potential -1 programmed frameshift (-1 FS) signal within the genome of the whitefly Bemisia-associated dicistrovirus 2 (BaDV-2). Here, using bicistronic reporters, we demonstrate that the predicted BaDV-2 -1 FS signal can drive -1 frameshifting in vitro via a slippery sequence and a downstream stem-loop structure that would direct the translation of the viral RNA-dependent RNA polymerase. Moreover, the predicted BaDV-2 IGR can support IRES translation in vitro but does so through a mechanism that is not typical of known factorless dicistrovirus IGR IRES mechanisms. Using deletion and mutational analyses, the BaDV-2 IGR IRES is mapped within a 140-nucleotide element and initiates translation from an AUG codon. Moreover, the IRES does not bind directly to purified ribosomes and is sensitive to eIF2 and eIF4A inhibitors NSC1198983 and hippuristanol, respectively, indicating an IRES-mediated factor-dependent mechanism. Biophysical characterization suggests the BaDV-2 IGR IRES contains several stem-loops; however, mutational analysis suggests a model whereby the IRES is unstructured or adopts distinct conformations for translation initiation. In summary, we have provided evidence of the first -1 FS frameshifting signal and a novel factor-dependent IRES mechanism in this dicistrovirus family, thus highlighting the diversity of viral RNA-structure strategies to direct viral protein synthesis.

Resistance to Cassava Whitefly (Bemisia tabaci) among Eastern and Southern African Elite Cassava Genotypes

Cassava whitefly, Bemisia tabaci, directly damages cassava leaves by feeding on phloem, causing chlorosis and abscission, leading to a yield loss of up to 50%. The pest also causes indirect damage through sooty mold formation. Most Ugandan cassava varieties resist cassava mosaic disease (CMD) and tolerate cassava brown streak disease (CBSD), but little is known about their response to whitefly infestation. The main objective of this study was to identify cassava genotypes with putative resistance to whitefly in Uganda. This was conducted on 24 improved cassava varieties in three agro-ecological zones during the second rains of 2016. Monthly data were taken for adult and nymph counts, whitefly and sooty mold damage, and CMD and CBSD severities from 2 to 9 months after planting (MAPs). The results show that the whitefly population is highly significantly (p < 0.000) amongst varieties across the three agro-ecological zones. Mkumba consistently supported the low adult numbers and nymphs. The findings demonstrate the potential of the improved cassava varieties as sources of whitefly resistance for sustainable management.

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.

Microbiome diversity and composition in Bemisia tabaci SSA1-SG1 whitefly are influenced by their host's life stage

Within the Bemisia tabaci group of cryptic whitefly species, many are damaging agricultural pests and plant-virus vectors, conferring upon this group the status of one of the world's top 100 most invasive and destructive species, affecting farmers' income and threatening their livelihoods. Studies on the microbiome of whitefly life stages are scarce, although their composition and diversity greatly influence whitefly fitness and development. We used high-throughput sequencing to understand microbiome diversity in different developmental stages of the B. tabaci sub-Saharan Africa 1 (SSA1-SG1) species of the whitefly from Uganda. Endosymbionts (Portiera, Arsenophonus, Wolbachia, and Hemipteriphilus were detected but excluded from further statistical analysis as they were not influenced by life stage using Permutational Multivariate Analysis of Variance Using Distance Matrices (ADONIS, p = 0.925 and Bray, p = 0.903). Our results showed significant differences in the meta microbiome composition in different life stages of SSA1-SG1. The diversity was significantly higher in eggs (Shannon, p = 0.024; Simpson, p = 0.047) than that in nymphs and pupae, while the number of microbial species observed by the amplicon sequence variant (ASV) was not significant (n(ASV), p = 0.094). At the phylum and genus levels, the dominant constituents in the microbiome changed significantly during various developmental stages, with Halomonas being present in eggs, whereas Bacillus and Caldalkalibacillus were consistently found across all life stages. These findings provide the first description of differing meta microbiome diversity in the life stage of whiteflies, suggesting their putative role in whitefly development.

Integrative transcriptomics reveals association of abscisic acid and lignin pathways with cassava whitefly resistance

BACKGROUND

Whiteflies are a global threat to crop yields, including the African subsistence crop cassava (Manihot esculenta). Outbreaks of superabundant whitefly populations throughout Eastern and Central Africa in recent years have dramatically increased the pressures of whitefly feeding and virus transmission on cassava. Whitefly-transmitted viral diseases threaten the food security of hundreds of millions of African farmers, highlighting the need for developing and deploying whitefly-resistant cassava. However, plant resistance to whiteflies remains largely poorly characterized at the genetic and molecular levels. Knowledge of cassava-defense programs also remains incomplete, limiting characterization of whitefly-resistance mechanisms. To better understand the genetic basis of whitefly resistance in cassava, we define the defense hormone- and Aleurotrachelus socialis (whitefly)-responsive transcriptome of whitefly-susceptible (COL2246) and whitefly-resistant (ECU72) cassava using RNA-seq. For broader comparison, hormone-responsive transcriptomes of Arabidopsis thaliana were also generated.

RESULTS

Whitefly infestation, salicylic acid (SA), jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) transcriptome responses of ECU72 and COL2246 were defined and analyzed. Strikingly, SA responses were largely reciprocal between the two cassava genotypes and we suggest candidate regulators. While susceptibility was associated with SA in COL2246, resistance to whitefly in ECU72 was associated with ABA, with SA-ABA antagonism observed. This was evidenced by expression of genes within the SA and ABA pathways and hormone levels during A. socialis infestation. Gene-enrichment analyses of whitefly- and hormone-responsive genes suggest the importance of fast-acting cell wall defenses (e.g., elicitor recognition, lignin biosynthesis) during early infestation stages in whitefly-resistant ECU72. A surge of ineffective immune and SA responses characterized the whitefly-susceptible COL2246's response to late-stage nymphs. Lastly, in comparison with the model plant Arabidopsis, cassava's hormone-responsive genes showed striking divergence in expression.

CONCLUSIONS

This study provides the first characterization of cassava's global transcriptome responses to whitefly infestation and defense hormone treatment. Our analyses of ECU72 and COL2246 uncovered possible whitefly resistance/susceptibility mechanisms in cassava. Comparative analysis of cassava and Arabidopsis demonstrated that defense programs in Arabidopsis may not always mirror those in crop species. More broadly, our hormone-responsive transcriptomes will also provide a baseline for the cassava community to better understand global responses to other yield-limiting pests/pathogens.

Cutting Dipping Application of Flupyradifurone against Cassava Whiteflies Bemisia tabaci and Impact on Its Parasitism in Cassava

The cassava whitefly Bemisia tabaci causes damage in cassava through the feeding and vectoring of plant viruses that cause cassava mosaic and cassava brown streak diseases. This study sought to explore the efficacy of cutting dipping in flupyradifurone for whitefly control and the impact of the mode of application on whitefly parasitism under farmer field conditions. The insecticide treatment significantly reduced adult whiteflies by 41%, nymphs by 64%, and cassava mosaic disease (CMD) incidence by 16% and increased root yield by 49%. The whitefly parasitism rate by Encarsia spp. parasitoids was 27.3 and 21.1%, while Eretmocerus spp. had 26.7 and 18.0% in control and flupyradifurone, respectively, and these differences were not significant. Electropenetrography recordings of whitefly feeding behaviour on flupyradifurone-treated plants showed significantly reduced probing activity and a delay in reaching the phloem as compared to the control. The findings from this study demonstrated that cassava cutting dipping in flupyradifurone significantly reduces whitefly numbers and cassava mosaic disease incidence, thus contributing to a significant root yield increase in cassava. Flupyradifurone applied through cutting dips does not significantly impact parasitism rates in cassava fields. Routine monitoring of parasitoids and predators in insecticide-treated versus control fields should be emphasized to determine the impact of pesticides on these beneficial non-target organisms.

Identifying potential insecticide resistance markers through genomic-level comparison of Bemisia tabaci (Gennadius) lines

The invasive whitefly (Bemisia tabaci) MED is one of the most economically damaging plant pests. The extensive use of insecticide over decades has led to that the invasive B. tabaci MED has developed resistance to a wide range of insecticide classes, but little is known about the genetic background associated with resistance. To this end, we conducted a comparative genome-wide analysis of single-base nucleotide polymorphisms between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line collected in 1976. First, low-coverage genome sequencings were conducted on DNA isolated from individual whiteflies. The sequencing results were evaluated using an available B. tabaci MED genome as a reference. Significant genetic differences were discovered between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line based on the principal component analyses. Top GO categories and KEGG pathways that might be involved in insecticide resistance development were identified, and several of them have not been previously associated with resistance. Additionally, we identified several genetic loci with novel variations including Cytochrome P450 monooxygenases (P450s), UDP-glucuronosyltransferases (UGTs), Glutathione S-transferases (GSTs), esterase, carboxyl-esterases (COE), ABC transporters, fatty acyl-CoA reductase, voltage-gated sodium channels, GABA receptor, and cuticle proteins (CPs) that were previously reported to have close associations with pesticide resistance in well-studied insect groups that provide an essential resource for the design of insecticide resistance-linked loci arrays insecticide. Our results was obtained solely on resequencing genome data sets, more pesticide bio-assays combined with omics datasets should be further used to verify the markers identified here.

Comparative analysis of salicylic acid levels and gene expression in resistant, tolerant, and susceptible cassava varieties following whitefly-mediated SLCMV infection

Sri Lankan cassava mosaic virus (SLCMV), the primary pathogen responsible for cassava mosaic disease in cassava plantations, is transmitted via infected cutting stems and the whitefly vector, Bemisia tabaci. To obtain better insights into the defense mechanism of cassava against SLCMV, whiteflies were used to induce SLCMV infection for activating the salicylic acid (SA) signaling pathway, which triggers the innate immune system. The study aimed to investigate the specific interactions between viruliferous whiteflies and SA accumulation in resistant (C33), tolerant (Kasetsart 50; KU50), and susceptible (Rayong 11) cassava cultivars by infecting with SLCMV. Leaf samples were collected at various time points, from 1 to 7 days after inoculation (dai). The SA levels were quantified by gas chromatography-mass spectrometry and validated by quantitative reverse transcription polymerase chain reaction. The SA levels increased in KU50 and C33 plants at 2 and 3 dai, respectively, but remained undetected in Rayong11 plants. The expression of PR-9e, PR-7f5, SPS1, SYP121, Hsf8, and HSP90 increased in infected C33 plants at 4 dai, whereas that of KU50 plants decreased immediately at 2 dai, and that of Rayong11 plants increased at 1 dai but gradually decreased thereafter. These findings strongly indicate that SA plays a crucial role in regulating antiviral defense mechanisms, especially in SLCMV-resistant plants. Altogether, the findings provide valuable insights into the mechanisms underlying the activation of SA-mediated anti-SLCMV defense pathways, and the resistance, tolerance, and susceptibility of cassava, which can aid future breeding programs aimed at enhancing SLCMV resistance.

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.

Cassava begomovirus species diversity changes during plant vegetative cycles

Cassava is a root crop important for global food security and the third biggest source of calories on the African continent. Cassava production is threatened by Cassava mosaic disease (CMD), which is caused by a complex of single-stranded DNA viruses (family: Geminiviridae, genus: Begomovirus) that are transmitted by the sweet potato whitefly (Bemisia tabaci). Understanding the dynamics of different cassava mosaic begomovirus (CMB) species through time is important for contextualizing disease trends. Cassava plants with CMD symptoms were sampled in Lake Victoria and coastal regions of Kenya before transfer to a greenhouse setting and regular propagation. The field-collected and greenhouse samples were sequenced using Illumina short-read sequencing and analyzed on the Galaxy platform. In the field-collected samples, African cassava mosaic virus (ACMV), East African cassava mosaic virus (EACMV), East African cassava mosaic Kenya virus (EACMKV), and East African cassava mosaic virus-Uganda variant (EACMV-Ug) were detected in samples from the Lake Victoria region, while EACMV and East African mosaic Zanzibar virus (EACMZV) were found in the coastal region. Many of the field-collected samples had mixed infections of EACMV and another begomovirus. After 3 years of regrowth in the greenhouse, only EACMV-like viruses were detected in all samples. The results suggest that in these samples, EACMV becomes the dominant virus through vegetative propagation in a greenhouse. This differed from whitefly transmission results. Cassava plants were inoculated with ACMV and another EACMV-like virus, East African cassava mosaic Cameroon virus (EACMCV). Only ACMV was transmitted by whiteflies from these plants to recipient plants, as indicated by sequencing reads and copy number data. These results suggest that whitefly transmission and vegetative transmission lead to different outcomes for ACMV and EACMV-like viruses.

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.

Manipulation of Whitefly Behavior by Plant Viruses

Whiteflies of the Bemisia tabaci complex transmit hundreds of plant viruses belonging to the genera Begomovirus and Crinivirus, among others. Tripartite interactions of whitefly-virus-plant frequently occur during virus infection and transmission. Specifically, virus transmission-related behavior of whitefly, such as preference and feeding, may be altered by viruses and thus exert significant impacts on the outcome of virus spread and epidemics. Here, we provide an overview on the current understanding of the manipulation of whitefly behavior by plant viruses. Plant viruses can significantly modulate whitefly preference and feeding behavior, either directly or in a plant-mediated manner. In general, non-viruliferous whiteflies tend to prefer virus-infected plants, and viruliferous whiteflies are more likely to prefer uninfected plants. In most cases, virus infection of plants and/or whitefly seems to exhibit positive or no effects on whitefly feeding on plants. The significance and evolution of these patterns are then discussed. Finally, we suggest several future directions of research, such as the exploration of temporal dynamics and the dissection of underlying mechanisms of virus-induced changes in whitefly behavior.

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.

Evaluation of the Efficacy of Flupyradifurone against Bemisia tabaci on Cassava in Tanzania

A novel butenolide insecticide-flupyradifurone (Sivanto SL 200)-was evaluated for efficacy against cassava-colonizing Bemisia tabaci whitefly under laboratory, screenhouse and field conditions. LC50 values from leaf disc spray assays were comparable for both flupyradifurone (12.7 g a.i/100 L) and imidacloprid (12.6 g a.i/100 L). Both insecticides caused high levels of adult whitefly mortality in leaf disc and leaf dip assays when compared to untreated controls. In screenhouse-based trials, longer soaking (60 min) with flupyradifurone or imidacloprid was more effective than shorter soaking durations (15 or 30 min). In field spraying experiments, flupyradifurone significantly reduced whiteflies, and both insecticides demonstrated powerful knockdown effects on whitefly adult abundances over a period up to 24 h. Single cutting dip application of flupyradifurone reduced whitefly adult abundance by 2 to 6 times, and nymphs by 2 to 13 times. Lower whitefly abundances resulting from insecticide application reduced the incidence of CMD or CBSD. In addition, in field experiments, whiteflies were fewer during the long rainy season (Masika) and on cassava variety Mkuranga1. The findings from this study demonstrate that cutting dips with flupyradifurone could be incorporated as a management tactic against cassava whiteflies. This would ideally be combined in an IPM strategy with other cassava virus and virus vector management tactics including host-plant resistance, phytosanitation and the use of clean seed.

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.

Modelling cassava production and pest management under biotic and abiotic constraints

We summarise modelling studies of the most economically important cassava diseases and arthropods, highlighting research gaps where modelling can contribute to the better management of these in the areas of surveillance, control, and host-pest dynamics understanding the effects of climate change and future challenges in modelling. For over 30 years, experimental and theoretical studies have sought to better understand the epidemiology of cassava diseases and arthropods that affect production and lead to considerable yield loss, to detect and control them more effectively. In this review, we consider the contribution of modelling studies to that understanding. We summarise studies of the most economically important cassava pests, including cassava mosaic disease, cassava brown streak disease, the cassava mealybug, and the cassava green mite. We focus on conceptual models of system dynamics rather than statistical methods. Through our analysis we identified areas where modelling has contributed and areas where modelling can improve and further contribute. Firstly, we identify research challenges in the modelling developed for the surveillance, detection and control of cassava pests, and propose approaches to overcome these. We then look at the contributions that modelling has accomplished in the understanding of the interaction and dynamics of cassava and its' pests, highlighting success stories and areas where improvement is needed. Thirdly, we look at the possibility that novel modelling applications can achieve to provide insights into the impacts and uncertainties of climate change. Finally, we identify research gaps, challenges, and opportunities where modelling can develop and contribute for the management of cassava pests, highlighting the recent advances in understanding molecular mechanisms of plant defence.

The role of pathogen-mediated insect superabundance in the East African emergence of a plant virus

One of the major crops for food security is cassava. Superabundant Bemisia tabaci whitefly, comprising unusually high landscape populations of the insect, have been implicated in cassava virus emergence. Studies have been unable to select from several hypotheses, however, as to the dynamic drivers of superabundant whitefly associated with the emergence in East Africa of severe cassava mosaic disease. One possibility is that pathogenic modification of infected plants can itself increase the growth of insect vector colonies on infected plants.Through the modelling of population processes at the landscape scale we introduce a framework for analysing patterns in the association of disease and insect waves.Our analyses demonstrate the role of pathogen-mediated insect superabundance in a plant disease invasion. Synthesis. An elevated abundance of insects at the landscape scale is frequently implicated in invasions of the plant pathogens that they carry. We advance ecological understanding of plant disease invasions by showing how landscape data can be used to investigate the causes of insect vector superabundance.

Early detection of plant virus infection using multispectral imaging and spatial-spectral machine learning

Cassava brown streak disease (CBSD) is an emerging viral disease that can greatly reduce cassava productivity, while causing only mild aerial symptoms that develop late in infection. Early detection of CBSD enables better crop management and intervention. Current techniques require laboratory equipment and are labour intensive and often inaccurate. We have developed a handheld active multispectral imaging (A-MSI) device combined with machine learning for early detection of CBSD in real-time. The principal benefits of A-MSI over passive MSI and conventional camera systems are improved spectral signal-to-noise ratio and temporal repeatability. Information fusion techniques further combine spectral and spatial information to reliably identify features that distinguish healthy cassava from plants with CBSD as early as 28 days post inoculation on a susceptible and a tolerant cultivar. Application of the device has the potential to increase farmers' access to healthy planting materials and reduce losses due to CBSD in Africa. It can also be adapted for sensing other biotic and abiotic stresses in real-world situations where plants are exposed to multiple pest, pathogen and environmental stresses.

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

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

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.

Phylogeographic and SNPs Analyses of Bemisia tabaci B Mitotype Populations Reveal Only Two of Eight Haplotypes Are Invasive

Simple Summary

The whitefly Bemisia tabaci taxon consists of an undefined number of morphologically identical genetic variants of which only a few, including the B, harbor invasive haplotypes. These haplotypes have potential to upsurge and become important pests and plant virus vectors in irrigated agroecosystems worldwide. In the 1980s, unprecedented outbreaks associated with the B variant were reported worldwide, however, the precise origin(s) of the invasive haplotypes has not been determined. In this study, available B. tabaci mitochondrial gene sequences were examined for patterns of conserved single nucleotide changes (SNPs). The whitefly sequence records represented North Africa-Middle Eastern habitats, the proposed B variant center of origin, and distant locales recently invaded by haplotype(s) of the B variant. Unexpectedly, the analysis revealed eight SNPs groups (haplotypes) demonstrating that the genetic architecture of the B mitoype is more complex than previously recognized. Also, the distribution patterns of the eight B haplotypes were tightly linked to well-defined eco-geographic regions, suggesting the different groups have diversified by geographic isolation. Contrary to claims that collectively, the B variant is invasive, only two of the eight haplotypic groups have established in geographical locations outside of their zone of endemism.

Abstract

The Bemisia tabaci cryptic species contains 39 known mitotypes of which the B and Q are best recognized for having established outside their extant endemic range. In the 1980s, previously uncharacterized haplotype(s) of the B mitotype rapidly established in tropical and subtropical locales distant from their presumed center of origin, leading to displacement of several native mitotypes and extreme damage to crops and other vegetation particularly in irrigated agroecosystems. To trace the natural and evolutionary history of the invasive B haplotypes, a phylo-biogeographic study was undertaken. Patterns of single nucleotide polymorphisms (SNPs) and signatures potentially indicative of geographic isolation were investigated using a globally representative mitochondrial cytochrome oxidase I gene (mtCOI) sequence database. Eight haplotype groups within the North Africa-Middle East (NAFME) region were differentiated, NAFME 1–8. The NAFME 1–3 haplotypes were members of the same population that is associated with warm desert climate niches of the Arabian Peninsula and east coastal Africa-Ethiopia. The NAFME 4 and 5 haplotypes are endemic to warm and cold semi-arid niches delimited by the Irano-Turanian floristic region, itself harboring extensive biodiversity. Haplotypes 6 and 7 co-occurred in the Middle East along eastern Mediterranean Sea landmasses, while NAFME 8 was found to be endemic to Cyprus, Turkey, and desert micro-niches throughout Egypt and Israel. Contrary to claims that collectively, the B mitotype is invasive, NAFME 6 and 8 are the only haplotypes to have established in geographical locations outside of their zone of endemism.

Expansion of the Cassava Brown Streak Disease Epidemic in Eastern Democratic Republic of Congo

Cassava plays a key role in ensuring food security and generating income for smallholder farmers throughout Central Africa, particularly in the Democratic Republic of Congo (DRC). This status is threatened, however, by cassava brown streak disease (CBSD), which has expanded its incidence and range in eastern DRC. The study described here comprises the first extensive assessment of temporal change in the occurrence of CBSD and its causal viruses in DRC, based on surveys conducted during 2016 and 2018. Cassava fields were inspected in Ituri, Nord-Kivu, Sud-Kivu, Tanganyika, and Haut-Katanga provinces within eastern DRC to record foliar incidence and severity of CBSD. Leaf samples were collected for virus detection and species-level identification. New occurrences of CBSD, confirmed by virus diagnostic tests, were recorded in two provinces (Haut-Katanga and Sud-Kivu) and nine previously unaffected territories, covering an area of >62,000 km², and at up to 900 km from locations of previously published reports of CBSD in DRC. Overall, average CBSD incidence within fields was 13.2% in 2016 and 16.1% in 2018. In the new spread zone of Haut-Katanga, incidence increased from 1.7 to 15.9%. CBSD is now present in provinces covering 321,000 km², which is approximately 14% of the total area of DRC. This represents a major expansion of the CBSD epidemic, which was only recorded from one province (Nord-Kivu) in 2012. Both Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus were detected in Ituri, Nord-Kivu, and Sud-Kivu, but only CBSV was detected in Haut-Katanga. Overall, these results confirm the increasing threat that CBSD poses to cassava production in DRC and describe an important expansion in the African pandemic of CBSD.

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.

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.

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.

Mitogenomic analysis of diversity of key whitefly pests in Kenya and its implication to their sustainable management

Whiteflies (Hemiptera: Aleyrodidae) are devastating agricultural pests of economic importance vectoring pathogenic plant viruses. Knowledge on their diversity and distribution in Kenya is scanty, limiting development of effective sustainable management strategies. The present study is aimed at identifying whitefly pest species present in Kenya across different agroecological zones and establish predictive models for the most abundant species in Africa. Whiteflies were sampled in Kenya from key crops known to be severely infested and identified using 16S rRNA markers and complete mitochondrial genomes. Four whitefly species were identified: Aleyrodes proletella, Aleurodicus dispersus, Bemisia afer and Trialeurodes vaporariorum, the latter being the most dominant species across all the agroecology. The assembly of complete mitogenomes and comparative analysis of all 13 protein coding genes confirmed the identities of the four species. Furthermore, prediction spatial models indicated high climatic suitability of T. vaporariorum in Africa, Europe, Central America, parts of Southern America, parts of Australia, New Zealand and Asia. Consequently, our findings provide information to guide biosecurity agencies on protocols to be adopted for precise identification of pest whitefly species in Kenya to serve as an early warning tool against T. vaporariorum invasion into unaffected areas and guide appropriate decision-making on their management.

Kawuki R. Cassava Brown Streak Disease Response and Association With Agronomic Traits in Elite Nigerian Cassava Cultivars

Cassava mosaic geminiviruses (CMGs) and cassava brown streak viruses (CBSVs) cause the highest yield losses in cassava production in Africa. In particular, cassava brown streak disease (CBSD) is and continues to be a significant constraint to optimal cassava production in Eastern and Southern Africa. While CBSD has not been reported in West Africa, its recent rapid spread and damage to cassava productivity in Eastern, and Southern Africa is alarming. The aim of this study was to evaluate Nigerian cassava genotypes in order to determine their responses to CBSD, in the event that it invades Nigeria, the world's largest cassava producer. The study gathered information on whether useful CBSD resistance alleles are present in the elite Nigerian cassava accessions. A total of 1,980 full-sib cassava seedlings from 106 families were assessed in the field at the seedling stage for a year. A subset of 569 clones were selected and assessed for another year at the clonal stage in Namulonge, central Uganda, a known hotspot for CBSD screening. Results indicated that foliar and root incidences and severities varied significantly (p ≤ 0.01, p ≤ 0.001) except for CBSD foliar incidence at 6 months (CBSD_6i ). Highest and lowest plot-based heritability estimates for CBSD were registered for CBSD root severity (CBSD _rs ) (0.71) and CBSD_6i (0.5). Positive and highly significant correlations were noted between CBSD root incidence (CBSD _ri ) and CBSD _rs (r = 0.90^***). Significant positive correlations were also noted between CBSD foliar severity at 3 months (CBSD_3s ) and CBSD foliar incidence at 6 months (CBSD_6i ) (r = 0.77^***), CBSD_3s and CBSD _rs (r = 0.35^***). Fresh root weight (Fresh _RW ) negatively correlated with CBSD _ri and CBSD _rs , respectively (r = -0.21^*** and r = -0.22^***). Similarly, CBSD_3s correlated negatively with cassava mosaic disease severity at 3 (CMD_3s ) and 6 months (CMD_6s ), respectively (r = -0.25^*** and r = -0.21^***). Fifteen clones were selected using a non-weighted summation selection index for further screening. In conclusion, results revealed that the elite Nigerian accessions exhibited significant susceptibility to CBSD within 2 years of evaluation period. It is expected that this information will aid future breeding decisions for the improvement of CBSD resistance among the Nigerian cassava varieties.

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.

Tracing the Invasion and Expansion Characteristics of the Flatid Planthopper, Metcalfa pruinosa (Hemiptera: Flatidae), in Korea Using Mitochondrial DNA Sequences

The flatid planthopper, Metcalfa pruinosa (Hemiptera: Flatidae), which is an invasive species, is widespread in Korea. We sequenced a fragment of the COI from 536 individuals collected mainly in Korea and the European countries and combined these sequence data with the public data, totaling 830 individuals worldwide. The identification of one shared haplotype only between Korea and the USA, the presence of this haplotype only in the North-West region of Korea, and the highest haplotype diversity in this region suggested that the North-West region is another point of entry in addition to the South-East region, which is the presumed sole point of entry to Korea. Furthermore, it suggested that North-West entry involves the M. pruinosa originating from the USA. In an effort to find further variable regions in the mitochondrial genome, one region provided substantially increased variability compared to that of the fragment of COI. F ST estimation, PCoA, and BAPS analysis, using the concatenated sequences of COI and the newly detected variable region to infer the expansion pattern in Korea, indicates that the main highway, running obliquely between the North-West and South-East regions, appears to be responsible for the current population genetic structure of M. pruinosa in Korea, facilitating gene flow through this highway traffic.

Improving climate suitability for Bemisia tabaci in East Africa is correlated with increased prevalence of whiteflies and cassava diseases

Projected climate changes are thought to promote emerging infectious diseases, though to date, evidence linking climate changes and such diseases in plants has not been available. Cassava is perhaps the most important crop in Africa for smallholder farmers. Since the late 1990’s there have been reports from East and Central Africa of pandemics of begomoviruses in cassava linked to high abundances of whitefly species within the Bemisia tabaci complex. We used CLIMEX, a process-oriented climatic niche model, to explore if this pandemic was linked to recent historical climatic changes. The climatic niche model was corroborated with independent observed field abundance of B. tabaci in Uganda over a 13-year time-series, and with the probability of occurrence of B. tabaci over 2 years across the African study area. Throughout a 39-year climate time-series spanning the period during which the pandemics emerged, the modelled climatic conditions for B. tabaci improved significantly in the areas where the pandemics had been reported and were constant or decreased elsewhere. This is the first reported case where observed historical climate changes have been attributed to the increase in abundance of an insect pest, contributing to a crop disease pandemic.

Disease Pandemics and Major Epidemics Arising from New Encounters between Indigenous Viruses and Introduced Crops

Virus disease pandemics and epidemics that occur in the world’s staple food crops pose a major threat to global food security, especially in developing countries with tropical or subtropical climates. Moreover, this threat is escalating rapidly due to increasing difficulties in controlling virus diseases as climate change accelerates and the need to feed the burgeoning global population escalates. One of the main causes of these pandemics and epidemics is the introduction to a new continent of food crops domesticated elsewhere, and their subsequent invasion by damaging virus diseases they never encountered before. This review focusses on providing historical and up-to-date information about pandemics and major epidemics initiated by spillover of indigenous viruses from infected alternative hosts into introduced crops. This spillover requires new encounters at the managed and natural vegetation interface. The principal virus disease pandemic examples described are two (cassava mosaic, cassava brown streak) that threaten food security in sub-Saharan Africa (SSA), and one (tomato yellow leaf curl) doing so globally. A further example describes a virus disease pandemic threatening a major plantation crop producing a vital food export for West Africa (cacao swollen shoot). Also described are two examples of major virus disease epidemics that threaten SSA’s food security (rice yellow mottle, groundnut rosette). In addition, brief accounts are provided of two major maize virus disease epidemics (maize streak in SSA, maize rough dwarf in Mediterranean and Middle Eastern regions), a major rice disease epidemic (rice hoja blanca in the Americas), and damaging tomato tospovirus and begomovirus disease epidemics of tomato that impair food security in different world regions. For each pandemic or major epidemic, the factors involved in driving its initial emergence, and its subsequent increase in importance and geographical distribution, are explained. Finally, clarification is provided over what needs to be done globally to achieve effective management of severe virus disease pandemics and epidemics initiated by spillover events.

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 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.

The potential geographical distribution and phenology of Bemisia tabaci Middle East/Asia Minor 1, considering irrigation and glasshouse production

The Bemisia tabaci species complex is one of the most important pests of open field and protected cropping globally. Within this complex, one species (Middle East Asia Minor 1, B. tabaci MEAM1, formerly biotype B) has been especially problematic, invading widely and spreading a large variety of plant pathogens, and developing broad spectrum pesticide resistance. Here, we fit a CLIMEX model to the distribution records of B. tabaci MEAM1, using experimental observations to calibrate its temperature responses. In fitting the model, we consider the effects of irrigation and glasshouses in extending its potential range. The validated niche model estimates its potential distribution as being considerably broader than its present known distribution, especially in the Americas, Africa and Asia. The potential distribution of the fitted model encompasses the known distribution of B. tabaci sensu lato, highlighting the magnitude of the threat posed globally by this invasive pest species complex and the viruses it vectors to open field and protected agriculture.

What is pathogen-mediated insect superabundance?

When increasing abundance of insect vectors is manifest across multiple fields of a crop at the landscape scale, the phenomenon is sometimes referred to as insect superabundance. The phenomenon may reflect environmental factors (i.e. environmentally mediated insect superabundance, EMiS), including climatic change. A number of pathogens, however, are also known to modify the quality of infected plants as a resource for their insect vectors. In this paper, we term increasing vector abundance when associated with pathogen modification of plants as pathogen-mediated insect superabundance (henceforth PMiS). We investigate PMiS using a new epidemiological framework. We formalize a definition of PMiS and indicate the epidemiological mechanism by which it is most likely to arise. This study is motivated by the occurrence of a particularly destructive cassava virus epidemic that has been associated with superabundant whitefly populations in sub-Saharan Africa. Our results have implications for how PMiS can be distinguished from EMiS in field data. Above all, they represent a timely foundation for further investigations into the association between insect superabundance and plant pathogens.

Genotype by environment cultivar evaluation for cassava brown streak disease resistance in Tanzania

Cassava brown steak disease (CBSD), caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), is the most important biotic constraint to cassava production in East and Central Africa. Concerted efforts are required to prevent further spread into West Africa as well as to reduce losses in areas already affected. The study reported here was part of a five-country (Kenya, Malawi, Mozambique, Tanzania and Uganda) programme that aimed to identify superior cassava cultivars resistant to CBSD and to disseminate them widely in the region. Seventeen tissue-cultured and virus-tested cultivars were evaluated in Tanzania across nine sites with diverse CBSD inoculum conditions. Experiments were planted using an alpha-lattice design and assessments were made of surrounding inoculum pressure, CBSD foliar and root incidence and root yield at harvest. There were large differences in CBSD infection between sites, with greatest spread recorded from the north-western Lake (Victoria) zone. Differences were driven by Bemisia tabaci whitefly vector abundance and CBSD inoculum pressure. Both CBSV and UCBSV were almost equally represented in cassava fields surrounding experimental plots, although CBSV predominated in the north-west whilst UCBSV was more frequent in coastal and southern sites. However, the incidence of CBSV was much greater than that of UCBSV in initially virus-free experimental plots, suggesting that CBSV is more virulent. Cultivars could be categorised into three groups based on the degree of CBSD symptom expression in shoots and roots. The seven cultivars (F10_30R2, Eyope, Mkumba, Mkuranga1, Narocass1, Nase3 and Orera) in the most resistant category each had shoot and root incidences of less than 20%. Fresh root yield differed between sites and cultivars, but there was no genotype by environment interaction for this trait, probably attributable to the large fertility and soil moisture differences between sites. Susceptible cultivars and the local check performed well in the absence of CBSD pressure, highlighting the importance of exploiting quality and yield traits of local landraces in breeding programmes. Overall, our results emphasized the importance of applying a balanced strategy for CBSD management. This should use both improved and local germplasm resources to generate high yielding cultivars for specific end-user traits, and combine the deployment of improved cultivars with phytosanitary control measures including the use of healthy planting material and planting during periods of reduced CBSD infection.

Screening for Resistance in Farmer-Preferred Cassava Cultivars from Ghana to a Mixed Infection of CBSV and UCBSV

Cassava brown streak disease (CBSD) caused by the Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) is a threat to cassava production in Africa. The potential spread of CBSD into West Africa is a cause for concern, therefore screening for resistance in farmer-preferred genotypes is crucial for effective control and management. We multiplied a selection of eleven cassava cultivars grown by farmers in Ghana to test their response to a mixed infection of CBSV (TAZ-DES-01) and UCBSV (TAZ-DES-02) isolates using a stringent top-cleft graft inoculation method. Virus titers were quantified in the inoculated scions and cuttings propagated from the inoculated scions to assess virus accumulation and recovery. All cultivars were susceptible to the mixed infection although their response and symptom development varied. In the propagated infected scions, CBSV accumulated at higher titers in leaves of eight of the eleven cultivars. Visual scoring of storage roots from six-month-old virus-inoculated plants revealed the absence of CBSD-associated necrosis symptoms and detectable titers of CBSVs in the cultivar, IFAD. Although all eleven cultivars supported the replication of CBSV and UCBSV in their leaves, the absence of virus replication and CBSD-associated symptoms in the roots of some cultivars could be used as criteria to rapidly advance durable CBSD tolerance using breeding and genetic engineering approaches.

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.

Genome-wide analyses of cassava Pathogenesis-related (PR) gene families reveal core transcriptome responses to whitefly infestation, salicylic acid and jasmonic acid

Background

Whiteflies are a threat to cassava (Manihot esculenta), an important staple food in many tropical/subtropical regions. Understanding the molecular mechanisms regulating cassava’s responses against this pest is crucial for developing control strategies. Pathogenesis-related (PR) protein families are an integral part of plant immunity. With the availability of whole genome sequences, the annotation and expression programs of the full complement of PR genes in an organism can now be achieved. An understanding of the responses of the entire complement of PR genes during biotic stress and to the defense hormones, salicylic acid (SA) and jasmonic acid (JA), is lacking. Here, we analyze the responses of cassava PR genes to whiteflies, SA, JA, and other biotic aggressors.

Results

The cassava genome possesses 14 of the 17 plant PR families, with a total of 447 PR genes. A cassava PR gene nomenclature is proposed. Phylogenetic relatedness of cassava PR proteins to each other and to homologs in poplar, rice and Arabidopsis identified cassava-specific PR gene family expansions. The temporal programs of PR gene expression in response to the whitefly (Aleurotrachelus socialis) in four whitefly-susceptible cassava genotypes showed that 167 of the 447 PR genes were regulated after whitefly infestation. While the timing of PR gene expression varied, over 37% of whitefly-regulated PR genes were downregulated in all four genotypes. Notably, whitefly-responsive PR genes were largely coordinately regulated by SA and JA. The analysis of cassava PR gene expression in response to five other biotic stresses revealed a strong positive correlation between whitefly and Xanthomonas axonopodis and Cassava Brown Streak Virus responses and negative correlations between whitefly and Cassava Mosaic Virus responses. Finally, certain associations between PR genes in cassava expansions and response to biotic stresses were observed among PR families.

Conclusions

This study represents the first genome-wide characterization of PR genes in cassava. PR gene responses to six biotic stresses and to SA and JA are demonstrably different to other angiosperms. We propose that our approach could be applied in other species to fully understand PR gene regulation by pathogens, pests and the canonical defense hormones SA and JA.

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.

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

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

Genetic Diversity of Bemisia tabaci (Hemiptera: Aleyrodidae) Species Complex Across Malaysia

The tobacco whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is a cryptic species complex with members capable of inducing huge economic losses. Precise identification of members of this complex proves essential in managing existing populations and preventing new incursions. Despite records of serious outbreaks of this pest in Malaysia little is known about species status of B. tabaci in this region. To address this, a comprehensive sampling of B. tabaci from different host plants was conducted in 10 states of Malaysia from 2010 to 2012. Members of the complex were identified by sequencing partial mitochondrial cytochrome oxidase subunit I (mtCOI) gene and constructing a Bayesian phylogenetic tree. Seven putative species were identified including Asia I, Mediterranean (MED), China 1, China 2, Asia II 6, Asia II 7, and Asia II 10. The most important finding of the study is the identification of the invasive MED species from locations without previous records of this species. All putative species except Asia I and MED are recorded from Malaysia for the first time. This study provided the first introductory map of B. tabaci species composition in Malaysia and emphasizes the urgent need for further studies to assess the status of MED invasion in this country.

Impact of Host Plant Species and Whitefly Species on Feeding Behavior of Bemisia tabaci

Whiteflies of the Bemisia tabaci species complex are economically important pests of cassava. In Africa, they cause greatest damage through vectoring viruses responsible for cassava mosaic disease and cassava brown streak disease. Several cryptic species from the B. tabaci complex colonize cassava and neighboring crops, but the feeding interactions between the different crops and B. tabaci species are unknown. The electrical penetration graph (EPG) technique makes it possible to conduct detailed feeding studies of sap-sucking insects by creating an electric circuit through the insect and the plant. The apparatus measures the voltage fluctuations while the wired-up insect feeds and produces graphs that describe feeding behavior. We utilized EPG to explore the feeding behavior of cassava-colonizing whiteflies (SSA1-SG3) on cassava, sweet potato, tomato, and cotton; and sweet potato-colonizing whiteflies (MED and IO) on cassava and sweet potato. Results show that: (1) feeding of SSA1-SG3 is not restricted to cassava. The least preferred host for SSA1-SG3 was tomato, where probing was delayed by 99 min compared to 10 min on other hosts, furthermore mean duration of phloem ingestion events was 36 min compared to 260 min on cassava. (2) Feeding of MED on cassava appeared to be non-functional, as it was characterized by short total phloem ingestion periods (<1 h) and few, short ingestion events, in contrast to feeding on sweet potato which was characterized by long phloem ingestion periods (>5 h). (3) Wire diameter affects the feeding in a statistically and practically significant manner. Implications for whitefly control and studies of host whitefly resistance are discussed.

Unravelling the Genetic Diversity among Cassava Bemisia tabaci Whiteflies Using NextRAD Sequencing

Bemisia tabaci threatens production of cassava in Africa through vectoring viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). B. tabaci sampled from cassava in eight countries in Africa were genotyped using NextRAD sequencing, and their phylogeny and population genetics were investigated using the resultant single nucleotide polymorphism (SNP) markers. SNP marker data and short sequences of mitochondrial DNA cytochrome oxidase I (mtCOI) obtained from the same insect were compared. Eight genetically distinct groups were identified based on mtCOI, whereas phylogenetic analysis using SNPs identified six major groups, which were further confirmed by PCA and multidimensional analyses. STRUCTURE analysis identified four ancestral B. tabaci populations that have contributed alleles to the six SNP-based groups. Significant gene flows were detected between several of the six SNP-based groups. Evidence of gene flow was strongest for SNP-based groups occurring in central Africa. Comparison of the mtCOI and SNP identities of sampled insects provided a strong indication that hybrid populations are emerging in parts of Africa recently affected by the severe CMD pandemic. This study reveals that mtCOI is not an effective marker at distinguishing cassava-colonizing B. tabaci haplogroups, and that more robust SNP-based multilocus markers should be developed. Significant gene flows between populations could lead to the emergence of haplogroups that might alter the dynamics of cassava virus spread and disease severity in Africa. Continuous monitoring of genetic compositions of whitefly populations should be an essential component in efforts to combat cassava viruses in Africa.