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

Rickettsia symbionts spread via mixed mode transmission, increasing female fecundity and sex ratio shift by host hormone modulating

Heritable symbionts are common among animals in nature, but the molecular mechanisms underpinning symbiont invasions of host populations have been elusive. In this study, we demonstrate the spread of Rickettsia in an invasive agricultural pest, the whitefly Bemisia tabaci Mediterranean (MED), across northeastern China from 2018 to 2023. Here, we show that the beneficial symbiont Rickettsia spreads by manipulating host hormone signals. Our analyses suggest that Rickettsia have been horizontally acquired by B. tabaci MED from another invasive whitefly B. tabaci Middle East-Asia Minor 1 during periods of coexistence. Rickettsia is transmitted maternally and horizontally from female B. tabaci MED individuals. Rickettsia infection enhances fecundity and results in female bias among whiteflies. Our findings reveal that Rickettsia infection stimulates juvenile hormone (JH) synthesis, in turn enhancing fecundity, copulation events, and the female ratio of the offspring. Consequently, Rickettsia infection results in increased whitefly fecundity and female bias by modulating the JH pathway. More female progeny facilitates the transmission of Rickettsia. This study illustrates that the spread of Rickettsia among invasive whiteflies in northeastern China is propelled by host hormone regulation. Such symbiont invasions lead to rapid physiological and molecular evolution in the host, influencing the biology and ecology of an invasive species.

Role of aminopeptidase N-like in the acquisition of begomoviruses by Bemisia tabaci, the whitefly vector

Sri Lankan cassava mosaic virus (SLCMV) is a prominent causative agent of cassava mosaic disease in Asia and relies on the whitefly Bemisia tabaci cryptic complex for its transmission. However, the molecular mechanisms involved in SLCMV transmission by B. tabaci have yet to be understood. In this study, we identified an aminopeptidase N-like protein (BtAPN) in B. tabaci Asia II 1, an efficient vector of SLCMV, which is involved in the SLCMV transmission process. Through the use of glutathione S-transferase pull-down assay and LC-MS/MS analysis, we demonstrated the interaction between BtAPN and the coat protein (CP) of SLCMV. This interaction was further confirmed in vitro, and we observed an induction of BtAPN gene expression following SLCMV infection. By interfering with the function of BtAPN, the quantities of SLCMV were significantly reduced in various parts of B. tabaci Asia II 1, including the whole body, midgut, hemolymph, and primary salivary gland. Furthermore, we discovered that BtAPN is conserved in B. tabaci Middle East-Asia Minor 1 (MEAM1) and interacts with the CP of tomato yellow leaf curl virus (TYLCV), a begomovirus known to cause severe damage to tomato production. Blocking BtAPN with antibody led to a significant reduction in the quantities of TYLCV in whitefly whole body and organs/tissues. These results demonstrate that BtAPN plays a generic role in interacting with the CP of begomoviruses and positively regulates their acquisition by the whitefly.

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

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

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

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

A shift of vector specificity acquired by a begomovirus through natural homologous recombination

Recombination is common in plant viruses such as geminiviruses, but the ecological and pathogenic consequences have been explored only in a few cases. Here, we found that a new begomovirus, tomato yellow leaf curl Shuangbai virus (TYLCSbV), probably originated from the recombination of Ageratum yellow vein China virus (AYVCNV) and tobacco curl shoot virus (TbCSV). Agrobacterium-mediated inoculation showed that TYLCSbV and AYVCNV have similar levels of infectivity on tomato and tobacco plants. However, the two viruses exhibit contrasting specificities for vector transmission, that is, TYLCSbV was efficiently transmitted by the whitefly Bemisia tabaci Mediterranean (MED) rather than by the whitefly B. tabaci Middle East-Asia Minor 1 (MEAM1), whereas AYVCNV was more efficiently transmitted by MEAM1. We also showed that the transmission efficiencies of TYLCSbV and AYVCNV are positively correlated with the accumulation of the viruses in whitefly whole bodies and organs/tissues. The key coat protein amino acids that determine their accumulation are between positions 147 and 256. Moreover, field surveys suggest that MED has displaced MEAM1 in some regions where TYLCSbV was collected. Viral competition assays indicated that TYLCSbV outcompeted AYVCNV when transmitted by MED, while the outcome was the opposite when transmitted by MEAM1. Our findings suggest that recombination has resulted in a shift of vector specificity that could provide TYLCSbV with a potential selective transmission advantage, and the population shift of whitefly cryptic species could have influenced virus evolution towards an extended trajectory of transmission.

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.

What Can Be Learned by a Synoptic Review of Plant Disease Epidemics and Outbreaks Published in 2021?

A synoptic review of plant disease epidemics and outbreaks was made using two complementary approaches. The first approach involved reviewing scientific literature published in 2021, in which quantitative data related to new plant disease epidemics or outbreaks were obtained via surveys or similar methodologies. The second approach involved retrieving new records added in 2021 to the CABI Distribution Database, which contains over a million global geographic records of organisms from over 50,000 species. The literature review retrieved 186 articles, describing studies in 62 categories (pathogen species/species complexes) across more than 40 host species on six continents. Pathogen species with more than five articles were Bursaphelenchus xylophilus, 'Candidatus Liberibacter asiaticus', cassava mosaic viruses, citrus tristeza virus, Erwinia amylovora, Fusarium spp. complexes, F. oxysporum f. sp. cubense, Magnaporthe oryzae, maize lethal necrosis co-infecting viruses, Meloidogyne spp. complexes, Pseudomonas syringae pvs., Puccinia striiformis f. sp. tritici, Xylella fastidiosa, and Zymoseptoria tritici. Automated searches of the CABI Distribution Database identified 617 distribution records new in 2021 of 283 plant pathogens. A further manual review of these records confirmed 15 pathogens reported in new locations: apple hammerhead viroid, apple rubbery wood viruses, Aphelenchoides besseyi, Biscogniauxia mediterranea, 'Ca. Liberibacter asiaticus', citrus tristeza virus, Colletotrichum siamense, cucurbit chlorotic yellows virus, Erwinia rhapontici, Erysiphe corylacearum, F. oxysporum f. sp. cubense Tropical race 4, Globodera rostochiensis, Nothophoma quercina, potato spindle tuber viroid, and tomato brown rugose fruit virus. Of these, four pathogens had at least 25% of all records reported in 2021. We assessed two of these pathogens-tomato brown rugose fruit virus and cucurbit chlorotic yellows virus-to be actively emerging in/spreading to new locations. Although three important pathogens-'Ca. Liberibacter asiaticus', citrus tristeza virus, and F. oxysporum f. sp. cubense-were represented in the results of both our literature review and our interrogation of the CABI Distribution Database, in general, our dual approaches revealed distinct sets of plant disease outbreaks and new records, with little overlap. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

The Dynamic Distribution of Wolbachia and Rickettsia in AsiaII1 Bemisia tabaci

Wolbachia and Rickettsia are bacterial endosymbionts that can induce a number of reproductive abnormalities in their arthropod hosts. We screened and established the co-infection of Wolbachia and Rickettsia in Bemisia tabaci and compared the spatial and temporal distribution of Wolbachia and Rickettsia in eggs (3-120 h after spawning), nymphs, and adults of B. tabaci by qPCR quantification and fluorescent in situ hybridization (FISH). The results show that the titer of Wolbachia and Rickettsia in the 3-120 h old eggs showed a "w" patterned fluctuation, while the titers of Wolbachia and Rickettsia had a "descending-ascending descending-ascending" change process. The titers of Rickettsia and Wolbachia nymphal and the adult life stages of Asia II1 B. tabaci generally increased with the development of whiteflies. However, the location of Wolbachia and Rickettsia in the egg changed from egg stalk to egg base, and then from egg base to egg posterior, and finally back to the middle of the egg. These results will provide basic information on the quantity and localization of Wolbachia and Rickettsia within different life stages of B. tabaci. These findings help to understand the dynamics of the vertical transmission of symbiotic bacteria.

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.

Transcription dynamics of heat-shock proteins (Hsps) and endosymbiont titres in response to thermal stress in whitefly, Bemisia tabaci (Asia-I)

The sweet potato whitefly, Bemisia tabaci (Gennadius), is one of the several species complexes of whitefly that are currently significant agricultural pests. Bemisia tabaci infests more than 600 plant species and thrives under a wide range of temperature conditions. In addition to the direct damage caused by sucking plant sap, it vectors several plant viruses. Heat-shock proteins play a pivotal role in enabling the insect to extend its geographical location, survival, and reproduction under different stress conditions. B. tabaci harbours several endosymbionts under the genera Portiera, Rickettsia, Hamiltonella, Wolbachia, Arsenophonus, Cardinium, and Fritschea that directly or indirectly affect its fitness. By accelerating cuticle biosynthesis and sclerotisation, symbiotic microbes can reduce or enhance tolerance to extreme temperatures and detoxify heavy metals. Thus, symbionts or microbial communities can expand or constrain the abiotic niche space of their host and affect its ability to adapt to changing conditions. The present study delineates the effect of thermal stress on the expression of heat-shock genes and endosymbionts in B. tabaci. Studies of the expression level of heat-shock proteins with the help of quantitative real-time polymerase chain reaction (qRT-PCR) showed that heat- and cold-shock treatment fuels the increased expression of heat-shock proteins (Hsp40 and Hsp70). However, Hsp90 was not induced by a heat- and cold-shock treatment. A significant decrease in the relative titre of secondary endosymbionts, such as Rickettsia, Arsenophonus, and Wolbachia, were recorded in B. tabaci upon heat treatment. However, the titre of the primary symbiont, C. Portiera, was relatively unaffected by both cold and heat treatments. These results are indicative of the fact that Hsp genes and endosymbionts in B. tabaci are modulated in response to thermal stress, and this might be responsible for the adaptation of whitefly under changing climatic scenario.

Evidence of the association between the Q2 mitochondrial group of Bemisia tabaci MED species (Hemiptera: Aleyrodidae) and low competitive displacement capability

The whitefly, Bemisia tabaci (Gennadius), is one of the most serious agricultural pests worldwide. Bemisia tabaci is a cryptic species complex of more than 40 species among which the invasive MEAM1 and MED species are the most widespread and economically important. Both MEAM1 and MED present intraspecific genetic variability and some haplotypes are reported to be more invasive than others. MED can be further deconstructed into different genetic groups, including MED-Q1 and MED-Q2. However, distinct biological phenotypes discerning the different MED mitochondrial haplotypes are yet to be characterized. Competitive displacement and life-history trials were carried out between MED-Q2 and MEAM1 populations collected in Florida, USA. In addition, a phylogenetic analysis was carried out including populations from previous whitefly competitive displacement studies for identification and comparison of the MED mitochondrial groups. In contrast to other studies with MED-Q1, the MED-Q2 population from Florida is less likely to displace MEAM1 on pepper. In addition, both pepper and watermelon were a more favorable host to MEAM1 compared to MED-Q2 according to the life history trials.

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.

Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene Superfamily in Bemisia tabaci and Expression Pattern Analysis under Heat Shock

Simple Summary

Bemisia tabaci MED is an invasive pest that had caused considerable economic damage in the past decades. Its successful colonization is closely related to heat-shock proteins (HSPs), which are related to heat resistance. In this study, 33 BtaHsps were identified based on the sequenced genome of B. tabaci MED belonging to six HSP families, among which 22 BtaHsps were newly identified. Analysis of the secondary structure and evolutionary relationship showed that they were all closely related. In addition, BtaHsp90A3 of the HSP90 family was screened by analyzing the expression level changes of these genes under 42 °C heat shock and RNAi was performed on the BtaHsp90A3. The results showed that the silencing of BtaHsp90A3 is closely related to the heat resistance of B. tabaci MED. Taken together, this study conducted an in-depth identification of BtaHsps that clarifies their evolutionary relationships and their response to thermal stress in B. tabaci MED.

Abstract

The thermal tolerance of Bemisia tabaci MED, an invasive whitefly species with worldwide distribution, plays an important role in its ecological adaptation during the invasion process. Heat-shock proteins (HSPs) are closely related to heat resistance. In this study, 33 Hsps (BtaHsps) were identified based on sequenced genome of B. tabaci MED belonging to six HSP families, among which 22 Hsps were newly identified. The secondary structures of a further 22 BtaHsps were also predicted. The results of RT-qPCR showed that heat shock could affect the expression of 14 of the 22 Hsps newly identified in this study. Among them, the expression level of six Hsps increased under 42 °C treatment. As the unstudied gene, BtaHsp90A3 had the highest increase rate. Therefore, BtaHsp90A3 was chosen for the RNAi test, and silencing BtaHsp90A3 by RNAi decreased the survival rate of adult B. tabaci at 42 °C. The results indicated that only a few Hsps were involved in the thermal tolerance of host whitefly although many Hsps would response under heat stress. This study conducted a more in-depth and comprehensive identification that demonstrates the evolutionary relationship of BtaHsps and illustrates the response of BtaHsps under the influence of thermal stress in B. tabaci MED.

Unraveling the Basis of Neonicotinoid Resistance in Whitefly Species Complex: Role of Endosymbiotic Bacteria and Insecticide Resistance Genes

Bemisia tabaci (whitefly) is one of the most detrimental agricultural insect pests and vectors of many plant viruses distributed worldwide. Knowledge of the distribution patterns and insecticide resistance of this cryptic species is crucial for its management. In this study, genetic variation of mitochondrial cytochrome oxidase subunit 1 (MtCoI) gene of B. tabaci was analyzed followed by a study of the infection profile of various endosymbionts in 26 whitefly populations collected from West Bengal, India. Phylogenetic analysis revealed Asia I as the major cryptic species (65.38%), followed by Asia II 5, China 3, and Asia II 7, which were diversified into 20 different haplotypes. In addition to the primary endosymbiont (C. poriera), each of the four whitefly species showed a variable population of three secondary endosymbionts, majorly Arsenophonus with the highest infection rate (73.07%), followed by Wolbachia and Rickettsia. Further phylogenetic analyses revealed the presence of two subgroups of Arsenophonus, viz., A1 and A2, and one each in Wolbachia (W1) and Rickettsia (R3). Resistance to thiamethoxam, imidacloprid, and acetamiprid insecticides was analyzed for a clear picture of pesticide resistance status. The highest susceptibility was noted toward thiamethoxam (LC₅₀ = 5.36 mg/L), followed by imidacloprid and acetamiprid. The whitefly population from Purulia and Hooghly districts bearing Asia II 7 and Asia II 5 cryptic species, respectively, shows maximum resistance. The differences in mean relative titer of four symbiotic bacteria among field populations varied considerably; however, a significant positive linear correlation was observed between the resistance level and relative titer of Arsenophonus and Wolbachia in the case of imidacloprid and thiamethoxam, while only Wolbachia was found in case of acetamiprid. Expression analysis demonstrated differential upregulation of insecticide resistance genes with Purulia and Hooghly populations showing maximally upregulated P450 genes. Moreover, thiamethoxam and imidacloprid resistance ratio (RR) showed a significant correlation with CYP6CM1, CYP6DZ7, and CYP4C64 genes, while acetamiprid RR correlated with CYP6CX1, CYP6DW2, CYP6DZ7, and CYP4C64 genes. Taken together, these findings suggested that P450 mono-oxygenase and symbiotic bacteria together affected whitefly resistance to neonicotinoids. Hence, a symbiont-oriented management programme could be a better alternative to control or delay resistance development in whitefly and can be used for pesticide clean-up in an agricultural field.

Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda

Bean leaf beetle (BLB) (Ootheca mutabilis) has emerged as an important bean pest in Uganda, leading to devastating crop losses. There is limited information on the population genetic structure of BLB despite its importance. In this study, novel microsatellite DNA markers and the partial mitochondrial cytochrome oxidase subunit I (mtCOI) gene sequences were used to analyze the spatial population genetic structure, genetic differentiation and haplotype diversity of 86 O. mutabilis samples from 16 (districts) populations. We identified 19,356 simple sequence repeats (SSRs) (mono, di-, tri-, tetra-, penta-, and hexa-nucleotides) of which 81 di, tri and tetra-nucleotides were selected for primer synthesis. Five highly polymorphic SSR markers (4-21 alleles, heterozygosity 0.59-0.84, polymorphic information content (PIC) 50.13-83.14%) were used for this study. Analyses of the 16 O. mutabilis populations with these five novel SSRs found nearly all the genetic variation occurring within populations and there was no evidence of genetic differentiation detected for both types of markers. Also, there was no evidence of isolation by distance between geographical and genetic distances for SSR data and mtCOI data except in one agro-ecological zone for mtCOI data. Bayesian clustering identified a signature of admixture that suggests genetic contributions from two hypothetical ancestral genetic lineages for both types of markers, and the minimum-spanning haplotype network showed low differentiation in minor haplotypes from the most common haplotype with the most common haplotype occurring in all the 16 districts. A lack of genetic differentiation indicates unrestricted migrations between populations. This information will contribute to the design of BLB control strategies.

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.

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

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

Characterization of transposable elements within the Bemisia tabaci species complex

Background

Whiteflies are agricultural pests that cause negative impacts globally to crop yields resulting at times in severe economic losses and food insecurity. The Bemisia tabaci whitefly species complex is the most damaging in terms of its broad crop host range and its ability to serve as vector for over 400 plant viruses. Genomes of whiteflies belonging to this species complex have provided valuable genomic data; however, transposable elements (TEs) within these genomes remain unexplored. This study provides the first accurate characterization of TE content within the B. tabaci species complex.

Results

This study identified that an average of 40.61% of the genomes of three whitefly species (MEAM1, MEDQ, and SSA-ECA) consists of TEs. The majority of the TEs identified were DNA transposons (22.85% average) while SINEs (0.14% average) were the least represented. This study also compared the TE content of the three whitefly genomes with three other hemipteran genomes and found significantly more DNA transposons and less LINEs in the whitefly genomes. A total of 63 TE superfamilies were identified to be present across the three whitefly species (39 DNA transposons, six LTR, 16 LINE, and two SINE). The sequences of the identified TEs were clustered which generated 5766 TE clusters. A total of 2707 clusters were identified as uniquely found within the whitefly genomes while none of the generated clusters were from both whitefly and non-whitefly TE sequences.

This study is the first to characterize TEs found within different B. tabaci species and has created a standardized annotation workflow that could be used to analyze future whitefly genomes.

Conclusion

This study is the first to characterize the landscape of TEs within the B. tabaci whitefly species complex. The characterization of these elements within the three whitefly genomes shows that TEs occupy significant portions of B. tabaci genomes, with DNA transposons representing the vast majority. This study also identified TE superfamilies and clusters of TE sequences of potential interest, providing essential information, and a framework for future TE studies within this species complex.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00270-6.

A novel microRNA regulates cooperation between symbiont and a laterally acquired gene in the regulation of pantothenate biosynthesis within Bemisia tabaci whiteflies

Horizontally transferred genes (HTGs) play a key role in animal symbiosis, and some horizontally transferred genes or proteins are highly expressed in specialized host cells (bacteriocytes). However, it is not clear how HTGs are regulated, but miRNAs are prime candidates given their previously demonstrated roles in symbiosis and impacts on the expression of host genes. A horizontally acquired PanBC that is highly expressed in whitefly bacteriocytes can cooperate with an obligate symbiont Portiera for pantothenate production, facilitating whitefly performance and Portiera titer. Here, we found that a whitefly miRNA novel-m0780-5p was up-regulated and its target panBC was down-regulated in Portiera-eliminated whiteflies. This miRNA was located to the cytoplasmic region of whitefly bacteriocytes. A novel-m0780-5p agomir injection reduced the expression of PanBC in whitefly bacteriocytes, while a novel-m0780-5p antagomir injection enhanced PanBC expression. Agomir injection also reduced the pantothenate level, Portiera titer and whitefly performance. Supplementation with pantothenate restored Portiera titer and the fitness of agomir-injected whiteflies. Thus, we demonstrated that a whitefly miRNA regulates panBC-mediated host-symbiont collaboration required for pantothenate synthesis, benefiting the whitefly-Portiera symbiosis. Both panBC and novel-m0780-5p are present in the genomes of six B. tabaci species. The expression of a novel miRNA in multiple B. tabaci species suggests that the miRNA evolved after panBC acquisition, and allowed this gene to be more tightly regulated. Our discovery provides the first account of a HTG being regulated by a miRNA from the host genome, and suggests key roles for interactions between miRNAs and HTGs in the functioning of symbiosis.

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.

Epidemiological assessment of cassava mosaic disease in Burkina Faso

Surveys were conducted in 2016 and 2017 across the main cassava-growing regions of Burkina Faso to assess the status of cassava mosaic disease (CMD) and to determine the virus strains causing the disease, using field observation and phylogenetic analysis. CMD incidence varied between regions and across years but was lowest in Hauts-Bassins (6.0%, 2016 and 5.4%, 2017) and highest in Centre-Sud (18.5%, 2016) and in Boucle du Mouhoun (51.7%, 2017). The lowest CMD severity was found in Est region (2.0) for both years and the highest in Sud-Ouest region (3.3, 2016) and Centre-Sud region (2.8, 2017). The CMD infection was primarily associated with contaminated cuttings in all regions except in Hauts-Bassins, where whitefly-borne infection was higher than cuttings-borne infection in 2016. PCR screening of 687 samples coupled with sequence analysis revealed the presence of African cassava mosaic-like (ACMV-like) viruses and East African cassava mosaic-like (EACMV-like) viruses as single infections at 79.5% and 1.1%, respectively. Co-infections of ACMV-like and EACMV-like viruses were detected in 19.4% of the tested samples. In addition, 86.7% of the samples positive for EACMV-like virus were found to be positive for East African cassava mosaic Cameroon virus (EACMCMV). Phylogenetic analysis revealed the segregation of cassava mosaic geminiviruses (CMGs) from Burkina Faso into three clades specific to ACMV, African cassava mosaic Burkina Faso virus (ACMBFV), and EACMCMV, confirming the presence of these viruses. The results of this study show that EACMCMV occurrence may be more prevalent in Burkina Faso than previously thought.

Lysine provisioning by horizontally acquired genes promotes mutual dependence between whitefly and two intracellular symbionts

Horizontal gene transfer is widespread in insects bearing intracellular symbionts. Horizontally transferred genes (HTGs) are presumably involved in amino acid synthesis in sternorrhynchan insects. However, their role in insect-symbiont interactions remains largely unknown. We found symbionts Portiera, Hamiltonella and Rickettsia possess most genes involved in lysine synthesis in the whitefly Bemisia tabaci MEAM1 although their genomes are reduced. Hamiltonella maintains a nearly complete lysine synthesis pathway. In contrast, Portiera and Rickettsia require the complementation of whitefly HTGs for lysine synthesis and have lysE, encoding a lysine exporter. Furthermore, each horizontally transferred lysine gene of ten B. tabaci cryptic species shares an evolutionary origin. We demonstrated that Hamiltonella did not alter the titers of Portiera and Rickettsia or lysine gene expression of Portiera, Rickettsia and whiteflies. Hamiltonella also did not impact on lysine levels or protein localization in bacteriocytes harboring Portiera and ovaries infected with Rickettsia. Complementation with whitefly lysine synthesis HTGs rescued E. coli lysine gene knockout mutants. Silencing whitefly lysA in whiteflies harboring Hamiltonella reduced lysine levels, adult fecundity and titers of Portiera and Rickettsia without influencing the expression of Hamiltonella lysA. Furthermore, silencing whitefly lysA in whiteflies lacking Hamiltonella reduced lysine levels, adult fecundity and titers of Portiera and Rickettsia in ovarioles. Therefore, we, for the first time, demonstrated an essential amino acid lysine synthesized through HTGs is important for whitefly reproduction and fitness of both obligate and facultative symbionts, and it illustrates the mutual dependence between whitefly and its two symbionts. Collectively, this study reveals that acquisition of horizontally transferred lysine genes contributes to coadaptation and coevolution between B. tabaci and its symbionts.

Bemisia tabaci Vesicle-Associated Membrane Protein 2 Interacts with Begomoviruses and Plays a Role in Virus Acquisition

Begomoviruses cause substantial losses to agricultural production, especially in tropical and subtropical regions, and are exclusively transmitted by members of the whitefly Bemisia tabaci species complex. However, the molecular mechanisms underlying the transmission of begomoviruses by their whitefly vector are not clear. In this study, we found that B. tabaci vesicle-associated membrane protein 2 (BtVAMP2) interacts with the coat protein (CP) of tomato yellow leaf curl virus (TYLCV), an emergent begomovirus that seriously impacts tomato production globally. After infection with TYLCV, the transcription of BtVAMP2 was increased. When the BtVAMP2 protein was blocked by feeding with a specific BtVAMP2 antibody, the quantity of TYLCV in B. tabaci whole body was significantly reduced. BtVAMP2 was found to be conserved among the B. tabaci species complex and also interacts with the CP of Sri Lankan cassava mosaic virus (SLCMV). When feeding with BtVAMP2 antibody, the acquisition quantity of SLCMV in whitefly whole body was also decreased significantly. Overall, our results demonstrate that BtVAMP2 interacts with the CP of begomoviruses and promotes their acquisition by whitefly.

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.