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

An alternative Semi-Autotropic Hydroponics (SAH) substrate for cassava rapid propagation: A first study case

The expansion of Semi-Autotrophic Hydroponics technology to address the issue of multiplying and disseminating virus-free planting materials for vegetatively propagated crops is challenged by the utilization of imported substrate, namely, KlasmannTS3. In this study, we evaluated the growth parameters and cutting production of cassava genotypes during three subsequent plantlet production cycles using three single substrates, namely, KlasmannTS3 (K), vermiculite (V), and local peat (P), and three blended substrates. The blended substrates were a combination of 25% K and 75% P (K25P75), a combination of V and P at respective rates of 25% and 75% (V25P75), and respective rates of 10% and 90% (V10P90). All cuttings obtained in one plantlet production cycle were transplanted into the next. The multiplication rate of cutting from cycle 1 to 2 (R1) and cycle 2 to 3 (R2) was calculated as the ratios of the number of cuttings per the number of plantlets in each cycle. K and K25P75 led to similar R1 and R2, except with the genotype IBA961089A, where K25P75 led to a higher R1. Local peat and V solely showed similar cutting multiplication rates, and were lower than V25P75 and V10P90. Substrates with a higher cutting production also led to a higher plantlet height, leaf, and internode number. V and its combinations with local peat led to the densest plantlet root system. The performance of the substrates contrasted among the genotypes, but IBA961089A mostly outperformed the two other genotypes. We concluded that up to 75% of K and, to a lesser extent 75% of V, can be substituted by P without compromising cutting production. V and P should be combined instead of being used separately.

Genetic diversity and population structure of Uganda cassava germplasm

Cassava (Manihot esculenta Crantz) holds significant economic importance globally. Evaluating a diverse range of germplasm based on molecular characteristics not only enhances its preservation but also supports its utilization in breeding programs. In this study, we assessed genetic diversity and population structure among 155 cassava genotypes from Uganda using 5247 single nucleotide polymorphism (SNP) markers. Genotyping by sequencing (GBS) was employed for SNP discovery and to evaluate genetic diversity and population structure using the ADMIXTURE software. The cassava accessions comprised two populations: 49 accessions from Ugandan lines and 106 accessions resulting from crosses between South American and Ugandan lines. The average call rate of 96% was utilized to assess marker polymorphism. Polymorphic information content values of the markers ranged from 0.1 to 0.5 with an average of 0.4 which was moderately high. The principal component analysis (PCA) showed that the first two components captured ~ 24.2% of the genetic variation. The average genetic diversity was 0.3. The analysis of molecular variance (AMOVA) indicated that 66.02% and 33.98% of the total genetic variation occurred within accessions and between sub-populations, respectively. Five sub-populations were identified based on ADMIXTURE structure analysis (K = 5). Neighbor-joining tree and hierarchical clustering tree revealed the presence of three different groups which were primarily based on the source of the genotypes. The results suggested that there was considerable genetic variation among the cassava genotypes which is useful in cassava improvement and conservation efforts.

Fur microbiome as a putative source of symbiotic bacteria in sucking lice

Symbiosis between insects and bacteria has been established countless times. While it is well known that the symbionts originated from a variety of different bacterial taxa, it is usually difficult to determine their environmental source and a route of their acquisition by the host. In this study, we address this question using a model of Neisseriaceae symbionts in rodent lice. These bacteria established their symbiosis independently with different louse taxa (Polyplax, Hoplopleura, Neohaematopinus), most likely from the same environmental source. We first applied amplicon analysis to screen for candidate source bacterium in the louse environment. Since lice are permanent ectoparasites, often specific to the particular host, we screened various microbiomes associated with three rodent species (Microtus arvalis, Clethrionomys glareolus, and Apodemus flavicollis). The analyzed samples included fur, skin, spleen, and other ectoparasites sampled from these rodents. The fur microbiome data revealed a Neisseriaceae bacterium, closely related to the known louse symbionts. The draft genomes of the environmental Neisseriaceae, assembled from all three rodent hosts, converged to a remarkably small size of approximately 1.4 Mbp, being even smaller than the genomes of the related symbionts. Our results suggest that the rodent fur microbiome can serve as a source for independent establishment of bacterial symbiosis in associated louse species. We further propose a hypothetical scenario of the genome evolution during the transition of a free-living bacterium to the member of the rodent fur-associated microbiome and subsequently to the facultative and obligate louse symbionts.

Occurrence and Distribution of Major Cassava Pests and Diseases in Cultivated Cassava Varieties in Western Kenya

Cassava is an important food crop in western Kenya, yet its production is challenged by pests and diseases that require routine monitoring to guide development and deployment of control strategies. Field surveys were conducted in 2022 and 2023 to determine the prevalence, incidence and severity of cassava mosaic disease (CMD) and cassava brown streak disease (CBSD), whitefly numbers and incidence of cassava green mite (CGM) in six counties of western Kenya. Details of the encountered cassava varieties were carefully recorded to determine the adoption of improved varieties. A total of 29 varieties were recorded, out of which 13 were improved, although the improved varieties were predominant in 60% of fields and the most widely grown variety was MM96/4271. The CMD incidence was higher in 2022 (26.4%) compared to 2023 (10.1%), although the proportion of CMD attributable to whitefly infection was greater (50.6%) in 2023 than in 2022 (18.0%). The CBSD incidence in 2022 was 6.4%, while in 2023 it was 4.1%. The CMD incidence was significantly lower (5.9%) for the improved varieties than it was for the local varieties (35.9%), although the CBSD incidence did not differ significantly between the improved (2.3%) and local varieties (9.7%). Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) were both detected. Most infections were single CBSV infections (82.9%), followed by single UCBSV (34.3%) and coinfection with both viruses (16.7%). Whiteflies were more abundant in 2023, in which 28% of the fields had super-abundant populations of >100/plant, compared to 5% in 2022. KASP SNP genotyping designated 92.8% of the specimens as SSA-ECA for 2022, while it was 94.4% for 2023. The cassava green mite incidence was 65.4% in 2022 compared to 79.9% in 2023. This study demonstrates that cassava viruses, whiteflies and cassava green mites continue to be important constraints to cassava production in western Kenya, although the widespread cultivation of improved varieties is reducing the impact of cassava viruses. The more widespread application of high-quality seed delivery mechanisms could further enhance the management of these pests/diseases, coupled with wider application of IPM measures for whiteflies and mites.

Genetic Diversity of Whiteflies Colonizing Crops and Their Associated Endosymbionts in Three Agroecological Zones of Cameroon

Bemisia tabaci (Gennadius) is as a major pest of vegetable crops in Cameroon. These sap-sucking insects are the main vector of many viruses infecting plants, and several cryptic species have developed resistance against insecticides. Nevertheless, there is very little information about whitefly species on vegetable crops and the endosymbionts that infect them in Cameroon. Here, we investigated the genetic diversity of whiteflies and their frequency of infection by endosymbionts in Cameroon. Ninety-two whitefly samples were collected and characterized using mitochondrial cytochrome oxidase I (mtCOI) markers and Kompetitive Allele Specific PCR (KASP). The analysis of mtCOI sequences of whiteflies indicated the presence of six cryptic species (mitotypes) of Bemisia tabaci, and two distinct clades of Bemisia afer and Trialeurodes vaporariorum. Bemisia tabaci mitotypes identified included: MED on tomato, pepper, okra, and melon; and SSA1-SG1, SSA1-SG2, SSA1-SG5, SSA3, and SSA4 on cassava. The MED mitotype predominated in all regions on the solanaceous crops, suggesting that MED is probably the main phytovirus vector in Cameroonian vegetable cropping systems. The more diverse cassava-colonizing B. tabaci were split into three haplogroups (SNP-based grouping) including SSA-WA, SSA4, and SSA-ECA using KASP genotyping. This is the first time that SSA-ECA has been reported in Cameroon. This haplogroup is predominant in regions currently affected by the severe cassava mosaic virus disease (CMD) and cassava brown streak virus disease (CBSD) pandemics. Three endosymbionts including Arsenophonus, Rickettsia, and Wolbachia were present in female whiteflies tested in this study with varying frequency. Arsenophonus, which has been shown to influence the adaptability of whiteflies, was more frequent in the MED mitotype (75%). Cardinium and Hamiltonella were absent in all whitefly samples. These findings add to the knowledge on the diversity of whiteflies and their associated endosymbionts, which, when combined, influence virus epidemics and responses to whitefly control measures, especially insecticides.

Identification and Characterization of Hibiscus mutabilis Varieties Resistant to Bemisia tabaci and Their Resistance Mechanisms

Hibiscus mutabilis, the city flower of Chengdu, is culturally significant and has nutritional and medicinal benefits. However, frequent infestations of Bemisia tabaci have caused economic losses. This study aimed to identify insect-resistant H. mutabilis varieties. Over two years, varieties like Jinqiusong, Zuiyun, and Zuifurong showed moderate to high resistance based on reproductive indices. Assessments of antixenosis and developmental impacts revealed that adult B. tabaci exhibited low selectivity toward these resistant varieties, indicating a strong repellent effect. Gas chromatography-mass spectrometry analysis identified volatile organic compounds, such as alcohols, alkanes, and terpenes. Notably, 2-ethylhexanol and 6-methylheptanol exhibited repellent properties. Using nontargeted metabolomics, this study compared the metabolite profiles of the insect-resistant variety Jinqiusong (JQS), moderately resistant Bairihuacai (BRHC), and highly susceptible Chongbanbai (CBB) post B. tabaci infestation. Fifteen key metabolites were linked to resistance, emphasizing the phenylpropanoid biosynthesis pathway as crucial in defense. These findings offer a theoretical foundation for breeding insect-resistant H. mutabilis varieties and developing eco-friendly strategies against B. tabaci infestations.

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.

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.

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.

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.

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.

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.

Transcriptional analysis of Bemisia tabaci MEAM1 cryptic species under the selection pressure of neonicotinoids imidacloprid, acetamiprid and thiamethoxam

BACKGROUND

Neonicotinoids are widely applied in the control of the destructive agricultural pest Bemisia tabaci, and resistance against these chemicals has become a common, severe problem in the control of whiteflies. To investigate the molecular mechanism underlying resistance against nenonicotinoids in whiteflies, RNA-seq technology was applied, and the variation in the transcriptomic profiles of susceptible whiteflies and whiteflies selected by imidacloprid, acetamiprid and thiamethoxam treatment was characterized.

RESULTS

A total of 90.86 GB of clean sequence data were obtained from the 4 transcriptomes. Among the 16,069 assembled genes, 584, 110 and 147 genes were upregulated in the imidacloprid-selected strain (IMI), acetamiprid-selected strain (ACE), and thiamethoxam (THI)-selected strain, respectively, relative to the susceptible strain. Detoxification-related genes including P450s, cuticle protein genes, GSTs, UGTs and molecular chaperone HSP70s were overexpressed in the selected resistant strains, especially in the IMI strain. Five genes were downregulated in all three selected resistant strains, including 2 UDP-glucuronosyltransferase 2B18-like genes (LOC 109030370 and LOC 109032577).

CONCLUSIONS

Ten generations of selection with the three neonicotinoids induced different resistance levels and gene expression profiles, mainly involving cuticle protein and P450 genes, in the three selected resistant whitefly strains. The results provide a reference for research on resistance and cross-resistance against neonicotinoids in B. tabaci.

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.

Influence of Olfaction in Host-Selection Behavior of the Cassava Whitefly Bemisia tabaci

Cassava is a vital food-security crop in Sub-Saharan Africa. Cassava crops are, however, severely affected by viral diseases transmitted by members of the whitefly species complex Bemisia tabaci. We have here investigated the role of olfaction in host selection behavior of the cassava whitefly B. tabaci SSA-ESA biotype. Surprisingly, we find that the whiteflies appear to make little use of olfaction to find their favored host. The cassava whitely shows a highly reduced olfactory system, both at the morphological and molecular level. Whitefly antennae possess only 15 sensilla with possible olfactory function, and from the genome we identified just a handful of candidate chemoreceptors, including nine tuning odorant receptors, which would afford the whitefly with one of the smallest olfactomes identified from any insect to date. Behavioral experiments with host and non-host plants, as well as with identified specific volatiles from these sources, suggest that the few input channels present are primarily tuned toward the identification of unwanted features, rather than favored ones, a strategy quite unlike most other insects. The demonstrated repellence effect of specific volatile chemicals produced by certain plants unflavored by whiteflies suggests that intercropping with these plants could be a viable strategy to reduce whitefly infestations in cassava fields.

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.

Optimizing Efficient RNAi-Mediated Control of Hemipteran Pests (Psyllids, Leafhoppers, Whitefly): Modified Pyrimidines in dsRNA Triggers

The advantages from exogenously applied RNAi biopesticides have yet to be realized in through commercialization due to inconsistent activity of the dsRNA trigger, and the activity level of RNAi suppression. This has prompted research on improving delivery methods for applying exogenous dsRNA into plants and insects for the management of pests and pathogens. Another aspect to improve RNAi activity is the incorporation of modified 2′-F pyrimidine nucleotides into the dsRNA trigger. Modified dsRNA incorporating 32–55% of the 2′-F- nucleotides produced improved RNAi activity that increased insect mortality by 12–35% greater than non-modified dsRNA triggers of the same sequence. These results were repeatable across multiple Hemiptera: the Asian citrus psyllid (Diaphorina citri, Liviidae); whitefly (Bemisia tabaci, Aleyroididae); and the glassy-winged sharpshooter (Homalodisca vitripennis, Cicadellidae). Studies using siRNA with modified 2′-F- pyrimidines in mammalian cells show they improved resistance to degradation from nucleases, plus result in greater RNAi activity, due to increase concentrations and improved binding affinity to the mRNA target. Successful RNAi biopesticides of the future will be able to increase RNAi repeatability in the field, by incorporating modifications of the dsRNA, such as 2′-F- pyrimidines, that will improve delivery after applied to fruit trees or crop plants, with increased activity after ingestion by insects. Costs of RNA modification have decreased significantly over the past few years such that biopesticides can now compete on pricing with commercial chemical products.

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.

High-Density Genetic Linkage Map Construction and White Rot Resistance Quantitative Trait Loci Mapping for Genus Vitis Based on Restriction Site-Associated DNA Sequencing

Grape white rot (Coniothyrium diplodiella) is a major fungal disease affecting grape yield and quality. Quantitative trait locus (QTL) analysis is an important method for studying important horticultural traits of grapevine. This study was conducted to construct a high-density map and conduct QTL mapping for grapevine white rot resistance. A mapping population with 177 genotypes was developed from interspecific hybridization of a white rot-resistant cultivar (Vitis vinifera × V. labrusca 'Zhuosexiang') and white rot-susceptible cultivar (V. vinifera 'Victoria'). Single-nucleotide polymorphism (SNP) markers were developed by restriction site-associated DNA sequencing. The female, male, and integrated maps contained 2,501, 4,110, and 6,249 SNP markers with average genetic distances of adjacent markers of 1.25, 0.77, and 0.50 cM, respectively. QTL mapping was conducted based on white rot resistance identification of 177 individuals in July and August of 2017 and 2018. Notably, one stable QTL related to white rot resistance was detected and located on linkage group LG14. The phenotypic variance ranged from 12.93 to 13.43%. An SNP marker (chr14_3929380), which cosegregated with white rot resistance, was discovered and shows potential for use in marker-assisted selection to generate new grapevine cultivars with resistance to white rot.

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

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

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

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

SUPPLEMENTARY INFORMATION

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

Host plants and Wolbachia shape the population genetics of sympatric herbivore populations

Changing climate and land-use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host-plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome-wide SNP data revealed two major genetic clusters-psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non-crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co-occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host-plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non-crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

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.

Exploring the Use of High-Resolution Melting Analysis and Helicase-Dependent Amplification for Discrimination of Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species and Trialeurodes vaporariorum

The sweetpotato whitefly, Bemisia tabaci (Gennadius) (Hemiptera; Aleyrodidae), and greenhouse whitefly, Trialeurodes vaporariorum (Westwood) (Hemiptera: Aleyrodidae), are highly problematic plant pests and virus vectors with worldwide distributions. Identification of whitefly species is typically accomplished by observation of distinct morphological characters; however, because of morphological inconsistency and indistinguishability, the discrimination of B. tabaci species variants is dependent on molecular techniques based on genetic differences. New assays were designed for the detection of B. tabaci A, B, and Q mitotype groups, and T. vaporariorum. Specific primer sets were designed for amplification of the mitochondrial cytochrome c oxidase I gene of the four targets to perform in end-point PCR, real-time PCR coupled to high-resolution melting analysis (HRM), and the isothermal helicase-dependent amplification (HDA). Primer specificities were validated using end-point PCR, then tested in HRM and HDA. Bemisia tabaci A, B, and Q mitotypes, and T. vaporariorum-targeted primer sets discriminately amplified specimens of different populations within their target whitefly group. These tests provide three novel discrimination assays for the high-consequence, exotic B. tabaci B and Q groups, along with the native B. tabaci A group and T. vaporariorum.

Global Cropland Connectivity: A Risk Factor for Invasion and Saturation by Emerging Pathogens and Pests

The geographic pattern of cropland is an important risk factor for invasion and saturation by crop-specific pathogens and arthropods. Understanding cropland networks supports smart pest sampling and mitigation strategies. We evaluate global networks of cropland connectivity for key vegetatively propagated crops (banana and plantain, cassava, potato, sweet potato, and yam) important for food security in the tropics. For each crop, potential movement between geographic location pairs was evaluated using a gravity model, with associated uncertainty quantification. The highly linked hub and bridge locations in cropland connectivity risk maps are likely priorities for surveillance and management, and for tracing intraregion movement of pathogens and pests. Important locations are identified beyond those locations that simply have high crop density. Cropland connectivity risk maps provide a new risk component for integration with other factors-such as climatic suitability, genetic resistance, and global trade routes-to inform pest risk assessment and mitigation.

High-density genetic linkage map construction and cane cold hardiness QTL mapping for Vitis based on restriction site-associated DNA sequencing

Background

Cold hardiness is an important agronomic trait and can significantly affect grape production and quality. Until now, there are no reports focusing on cold hardiness quantitative trait loci (QTL) mapping. In this study, grapevine interspecific hybridisation was carried out with the maternal parent ‘Cabernet sauvignon’ and paternal parent ‘Zuoyouhong’. A total of 181 hybrid offspring and their parents were used as samples for restriction-site associated DNA sequencing (RAD). Grapevine cane phloem and xylem cold hardiness of the experimental material was detected using the low-temperature exotherm method in 2016, 2017 and 2018. QTL mapping was then conducted based on the integrated map.

Results

We constructed a high-density genetic linkage map with 16,076, 11,643, and 25,917 single-nucleotide polymorphism (SNP) markers anchored in the maternal, paternal, and integrated maps, respectively. The average genetic distances of adjacent markers in the maps were 0.65 cM, 0.77 cM, and 0.41 cM, respectively. Colinearity analysis was conducted by comparison with the grape reference genome and showed good performance. Six QTLs were identified based on the phenotypic data of 3 years and they were mapped on linkage group (LG) 2, LG3, and LG15. Based on QTL results, candidate genes which may be involved in grapevine cold hardiness were selected.

Conclusions

High-density linkage maps can facilitate grapevine fine QTL mapping, genome comparison, and sequence assembly. The cold hardiness QTL mapping and candidate gene discovery performed in this study provide an important reference for molecular-assisted selection in grapevine cold hardiness breeding.

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.

Two Distinct Genotypes of Spissistilus festinus (Say, 1830) (Hemiptera, Membracidae) in the United States Revealed by Phylogenetic and Morphological Analyses

Spissistilus festinus (Say, 1830) (Hemiptera: Membracidae) is a frequent pest of leguminous crops in the Southern United States, and a vector of grapevine red blotch virus. There is currently no information on the genetic diversity of S. festinus. In this study, populations of S. festinus were collected in 2015-2017 from various crops and geographic locations in the United States, and fragments of the mitochondrial cytochrome C oxidase 1 (mt-COI) gene and the nuclear internal transcribed spacer 2 (ITS2) region were characterized by polymerase chain reaction and sequencing. Maximum-likelihood and Bayesian analyses of the mt-COI and ITS2 sequences yielded similar phylogenetic tree topologies, revealing two distinct genetic S. festinus lineages with all of the specimens from California comprising one phylogenetic clade, alongside a single GenBank entry from Arizona, and all specimens from the Southeastern United States comprising a statistically-supported distinct clade, regardless of host and year of collection. The mt-COI gene fragment showed up to 10.8% genetic distance between the two phylogenetic clades. These results suggest the existence of two genotypes within S. festinus in the United States. The only distinct morphological trait between the two genotypes was a less elevated pronotum in the representative specimens from California, compared to the representative specimens from the Southeastern United States. Since this phenotypic feature is inconspicuous, a diagnostic polymerase chain reaction targeting a variable region of the mt-COI fragment was developed to reliably distinguish between the specimens of the two genotypes of S. festinus and to facilitate their specific identification.

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.

Minimum sample sizes for invasion genomics: Empirical investigation in an invasive whitefly

Analysis of population genetics provides insights into the evolutionary processes, among which the sample size choice is per se a crucial issue in the analysis. Genome-wide high-throughput techniques based on RADseq have been increasingly used in studies on the population genomics of invasive species. However, there is little information available regarding optimal sample sizes for analyzing population genomics of invasive species. In this study, we first use type IIB endonucleases restriction site-associated DNA (2b-RAD) to mine thousands of single nucleotide polymorphisms (SNPs) for native and introduced populations in Q1 clade (SPB and 17JN) and Q2 clade (ISQ and UAS0601) of the whitefly, Bemisia tabaci (Gennadius) MED (also known as B. tabaci biotype Q). Then, we used resampling techniques to create simulated populations with a random subset of individuals and 3,000 SNPs to determine how many individuals should be sampled for accurate estimates of intra- and interpopulation genetic diversity. We calculated the intrapopulation genetic diversity parameters (unbiased expected heterozygosity, observed heterozygosity, and the number of effect alleles) and pairwise genetic differentiation F ST; finally, an ad hoc statistic, ΔK, was used to determine the optimal value. Our results showed that a sample size greater than four individuals (n ≥ 4) has little impact on estimates of genetic diversity within whitefly populations; moreover, precise estimate of F ST can be easily achieved at a very small simple size (n = 3 or 4). Our results will provide in-depth understanding of the optimization of sampling schemes in population genomics of invasive species.

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.

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.

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

The whitefy Bemisia tabaci, a species complex consisting of many morphologically indistinguishable species divided into distinct clades, is one of the most globally important agricultural pests and plant virus vectors. Cassava-colonizing B. tabaci transmits viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). Half of all cassava plants in Africa are affected by these viral diseases, resulting in annual production losses of more than US$ 1 billion. Here we report the draft genome of the cassava whitefly B. tabaci Sub-Saharan Africa - East and Central Africa (SSA-ECA), the super-abundant population that has been associated with the rapid spread of viruses causing the pandemics of CMD and CBSD. The SSA-ECA genome assembled from Illumina short reads has a total size of 513.7 Mb and a scaffold N50 length of 497 kb, and contains 15,084 predicted protein-coding genes. Phylogenetic analysis suggests that SSA-ECA diverged from MEAM1 around 5.26 million years ago. A comprehensive genetic analysis of cassava-colonizing B. tabaci in Africa was also conducted, in which a total of 243 whitefly specimens were collected from 18 countries representing all major cassava-growing regions in the continent and genotyped using NextRAD sequencing. Population genomic analyses confirmed the existence of six major populations linked by gene flow and inferred the distribution patterns of these populations across the African continent. The genome of SSA-ECA and the genetic findings provide valuable resources and guidance to facilitate whitefly research and the development of strategies to control cassava viral diseases spread by whiteflies.

Going global - genomic insights into insect invasions

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

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.

Whitefly-transmitted viruses threatening cassava production in Africa

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

Investigation of Chinese Wolfberry (Lycium spp.) Germplasm by Restriction Site-Associated DNA Sequencing (RAD-seq)

Chinese wolfberry (Lycium spp.) is an important edible and medicinal plant, with a long cultivation history. The genetic relationships among wild Lycium species and landraces have been unclear for a number of reasons, which has hindered the breeding of modern Chinese wolfberry cultivars. In this study, we collected 19 accessions of Chinese wolfberry germplasm, and constructed the genetic relationship based on RAD-seq markers. We obtained 30.32 Gb of clean data, with the average value of each sample being 1.596 Gb. The average mapping rate was 85.7%, and the average coverage depth was 6.76 X. The phylogeny results distinguished all accessions clearly. All the studied landraces shared their most recent common ancestor with L. barbarum, which indicated that L. barbarum may be involved in cultivation of these landraces. The relationship of some landraces, namely the ‘Ningqi’ series, ‘Qingqi-1’ and ‘Mengqi-1,’ has been supported by the phylogeny results, while the triploid wolfberry was shown to be based on a hybrid between ‘Ningqi-1’ and a tetraploid wolfberry. This study uncovered the genetic background of Chinese wolfberry, and developed the foundation for species classification, accession identification and protection, and the production of hybrid cultivars of wolfberry.