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

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

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

Differential Transmission of Old and New World Begomoviruses by Middle East-Asia Minor 1 (MEAM1) and Mediterranean (MED) Cryptic Species of Bemisia tabaci.

Middle East-Asia Minor 1 (MEAM1) and Mediterranean (MED) are two of the most invasive members of the sweetpotato whitefly, Bemisia tabaci, cryptic species complexes and are efficient vectors of begomoviruses. Bemisia tabaci MEAM1 is the predominant vector of begomoviruses in open-field vegetable crops in the southeastern United States. However, recently B. tabaci MED also has been detected in the landscape outside of greenhouses in Florida and Georgia. This study compared the transmission efficiency of one Old-World (OW) and two New-World (NW) begomoviruses prevalent in the southeastern United States, viz., tomato yellow leaf curl virus (TYLCV), cucurbit leaf crumple virus (CuLCrV), and sida golden mosaic virus (SiGMV) between B. tabaci MEAM1 and B. tabaci MED. Bemisia tabaci MEAM1 efficiently transmitted TYLCV, CuLCrV, or SiGMV, whereas B. tabaci MED only transmitted TYLCV. Percent acquisition and retention of OW TYLCV following a 72 h acquisition access period was significantly higher for B. tabaci MED than B. tabaci MEAM1. In contrast, B. tabaci MEAM1 acquired and retained significantly more NW bipartite begomoviruses, CuLCrV or SiGMV, than B. tabaci MED. Quantitative analysis (qPCR) of virus DNA in whitefly internal tissues revealed reduced accumulation of CuLCrV or SiGMV in B. tabaci MED than in B. tabaci MEAM1. Fluorescent in situ hybridization (FISH) showed localization of CuLCrV or SiGMV in the midgut of B. tabaci MED and B. tabaci MEAM1. However, localization of CuLCrV or SiGMV was only observed in the primary salivary glands of B. tabaci MEAM1 and not B. tabaci MED. TYLCV localization was observed in all internal tissues of B. tabaci MEAM1 and B. tabaci MED. Overall, results demonstrate that both B. tabaci MEAM1 and B. tabaci MED are efficient vectors of OW TYLCV. However, for the NW begomoviruses, CuLCrV and SiGMV, B. tabaci MEAM1 seems to a better vector.

A Review on Transcriptional Responses of Interactions between Insect Vectors and Plant Viruses.

This review provides a synopsis of transcriptional responses pertaining to interactions between plant viruses and the insect vectors that transmit them in diverse modes. In the process, it attempts to catalog differential gene expression pertinent to virus-vector interactions in vectors such as virus reception, virus cell entry, virus tissue tropism, virus multiplication, and vector immune responses. Whiteflies, leafhoppers, planthoppers, and thrips are the main insect groups reviewed, along with aphids and leaf beetles. Much of the focus on gene expression pertinent to vector-virus interactions has centered around whole-body RNA extraction, whereas data on virus-induced tissue-specific gene expression in vectors is limited. This review compares transcriptional responses in different insect groups following the acquisition of non-persistent, semi-persistent, and persistent (non-propagative and propagative) plant viruses and identifies parallels and divergences in gene expression patterns. Understanding virus-induced changes in vectors at a transcriptional level can aid in the identification of candidate genes for targeting with RNAi and/or CRISPR editing in insect vectors for management approaches.

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.

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.

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.

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.

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.

Specific and Spillover Effects on Vectors Following Infection of Two RNA Viruses in Pepper Plants.

Mixed infection of plant viruses is ubiquitous in nature and can affect virus-plant-vector interactions differently than single virus infection. While several studies have examined virus-virus interactions involving mixed virus infection, relatively few have examined effects of mixed virus infection on vector preference and fitness, especially when multiple vectors are involved. This study explored how single and mixed viral infection of a non-persistently transmitted cucumber mosaic virus (CMV) and propagative and persistently-transmitted tomato spotted wilt orthotospovirus (TSWV) in pepper, Capsicum annum L., influenced the preference and fitness of their vectors, the green peach aphid, Myzus persicae (Sulzer), and the tobacco thrips, Frankliniella fusca (Hinds), respectively. In general, mixed infected plants exhibited severe symptoms compared with individually infected plants. An antagonistic interaction between the two viruses was observed when CMV titer was reduced following mixed infection with TSWV in comparison with the single infection. TSWV titer did not differ between single and mixed infection. Myzus persicae settling preference and median developmental were not significantly different between CMV and/or TSWV-infected and non-infected plants. Moreover, M. persicae fecundity did not differ between CMV-infected and non-infected pepper plants. However, M. persicae fecundity was substantially greater on TSWV-infected plants than non-infected plants. Myzus persicae fecundity on mixed-infected plants was significantly lower than on singly-infected and non-infected plants. Frankliniella fusca fecundity was higher on CMV and/or TSWV-infected pepper plants than non-infected pepper plants. Furthermore, F. fusca-induced feeding damage was higher on TSWV-infected than on CMV-infected, mixed-infected, or non-infected pepper plants. Overall, our results indicate that the effects of mixed virus infection on vectors were not different from those observed following single virus infection. Virus-induced host phenotype-modulated effects were realized on both specific and non-specific vectors, suggesting crosstalk involving all vectors and viruses in this pathosystem. The driving forces of these interactions need to be further examined. The effects of interactions between two viruses and two vectors towards epidemics of one or both viruses also need to be examined.

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

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

African ancestry of New World, Bemisia tabaci-whitefly species.

Bemisia tabaci whitefly species are some of the world's most devastating agricultural pests and plant-virus disease vectors. Elucidation of the phylogenetic relationships in the group is the basis for understanding their evolution, biogeography, gene-functions and development of novel control technologies. We report here the discovery of five new Sub-Saharan Africa (SSA) B. tabaci putative species, using the partial mitochondrial cytochrome oxidase 1 gene: SSA9, SSA10, SSA11, SSA12 and SSA13. Two of them, SSA10 and SSA11 clustered with the New World species and shared 84.8‒86.5% sequence identities. SSA10 and SSA11 provide new evidence for a close evolutionary link between the Old and New World species. Re-analysis of the evolutionary history of B. tabaci species group indicates that the new African species (SSA10 and SSA11) diverged from the New World clade c. 25 million years ago. The new putative species enable us to: (i) re-evaluate current models of B. tabaci evolution, (ii) recognise increased diversity within this cryptic species group and (iii) re-estimate divergence dates in evolutionary time.

Genetic diversity and geographic distribution of Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) genotypes associated with cassava in East Africa.

The genetic variability of whitefly (Bemisia tabaci) species, the vectors of cassava mosaic begomoviruses (CMBs) in cassava growing areas of Kenya, Tanzania, and Uganda, was investigated through comparison of partial sequences of the mitochondria cytochrome oxidase I (mtCOI) DNA in 2010/11. Two distinct species were obtained including sub-Saharan Africa 1 (SSA1), comprising of two sub-clades (I and II), and a South West Indian Ocean Islands (SWIO) species. Among the SSA1, sub-clade I sequences shared a similarity of 97.8-99.7% with the published Uganda 1 genotypes, and diverged by 0.3-2.2%. A pairwise comparison of SSA1 sub-clade II sequences revealed a similarity of 97.2-99.5% with reference southern Africa genotypes, and diverged by 0.5-2.8%. The SSA1 sub-clade I whiteflies were widely distributed in East Africa (EA). In comparison, the SSA1 sub-clade II whiteflies were detected for the first time in the EA region, and occurred predominantly in the coast regions of Kenya, southern and coast Tanzania. They occurred in low abundance in the Lake Victoria Basin of Tanzania and were widespread in all four regions in Uganda. The SWIO species had a sequence similarity of 97.2-97.7% with the published Reunion sequence and diverged by 2.3-2.8%. The SWIO whiteflies occurred in coast Kenya only. The sub-Saharan Africa 2 whitefly species (Ug2) that was associated with the severe CMD pandemic in Uganda was not detected in our study.