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

Multiplex Cas9-based excision of CLCuV betasatellite and DNA-A revealed reduction of viral load with asymptomatic cotton plants

MAIN CONCLUSION

Multiplexed Cas9-based genome editing of cotton resulted in reduction of viral load with asymptomatic cotton plants. In depth imaging of proteomic dynamics of resulting CLCuV betasatellite and DNA-A protein was also performed. The notorious  cotton leaf curl virus (CLCuV), which is transmitted by the sap-sucking insect whitefly, continuously damages cotton crops. Although the application of various toxins and RNAi has shown some promise, sustained control has not been achieved. Consequently, CRISPR_Cas9 was applied by designing multiplex targets against DNA-A (AC2 and AC3) and betasatellite (βC1) of CLCuV using CRISPR direct and ligating into the destination vector of the plant using gateway ligation method. The successful ligation of targets into the destination vector was confirmed by the amplification of 1049 bp using a primer created from the promoter and target, while restriction digestion using the AflII and Asc1 enzymes determined how compact the plasmid developed and the nucleotide specificity of the plasmid was achieved through Sanger sequencing. PCR confirmed the successful introduction of plasmid into CKC-1 cotton variety. Through Sanger sequencing and correlation with the mRNA expression of DNA-A and betasatellite in genome-edited cotton plants subjected to agroinfiltration of CLCuV infectious clone, the effectiveness of knockout was established. The genome-edited cotton plants demonstrated edited efficacy of 72% for AC2 and AC3 and 90% for the (βC1) through amplicon sequencing, Molecular dynamics (MD) simulations were used to further validate the results. Higher RMSD values for the edited βC1 and AC3 proteins indicated functional loss caused by denaturation. Thus, CRISPR_Cas9 constructs can be rationally designed using high-throughput MD simulation technique. The confidence in using this technology to control plant virus and its vector was determined by the knockout efficiency and the virus inoculation assay.

Sentinel plot surveillance of cotton leaf curl disease in Pakistan- a case study at the cultivated cotton-wild host plant interface

A sentinel plot case study was carried out to identify and map the distribution of begomovirus-betasatellite complexes in sentinel plots and commercial cotton fields over a four-year period using molecular and high-throughput DNA 'discovery' sequencing approaches. Samples were collected from 15 study sites in the two major cotton-producing areas of Pakistan. Whitefly- and leafhopper-transmitted geminiviruses were detected in previously unreported host plant species and locations. The most prevalent begomovirus was cotton leaf curl Kokhran virus-Burewala (CLCuKoV-Bu). Unexpectedly, a recently recognized recombinant, cotton leaf curl Multan virus-Rajasthan (CLCuMuV-Ra) was prevalent in five of 15 sites. cotton leaf curl Alabad virus (CLCuAlV) and cotton leaf curl Kokhran virus-Kokhran, 'core' members of CLCuD-begomoviruses that co-occurred with CLCuMuV in the 'Multan' epidemic were detected in one of 15 sentinel plots. Also identified were chickpea chlorotic dwarf virus and 'non-core' CLCuD-begomoviruses, okra enation leaf curl virus, squash leaf curl virus, and tomato leaf curl New Delhi virus. Cotton leaf curl Multan betasatellite (CLCuMuB) was the most prevalent CLCuD-betasatellite, and less commonly, two 'non-core' betasatellites. Recombination analysis revealed previously uncharacterized recombinants among helper virus-betasatellite complexes consisting of CLCuKoV, CLCuMuV, CLCuAlV and CLCuMuB. Population analyses provided early evidence for CLCuMuV-Ra expansion and displacement of CLCuKoV-Bu in India and Pakistan from 2012-2017. Identification of 'core' and non-core CLCuD-species/strains in cotton and other potential reservoirs, and presence of the now predominant CLCuMuV-Ra strain are indicative of ongoing diversification. Investigating the phylodynamics of geminivirus emergence in cotton-vegetable cropping systems offers an opportunity to understand the driving forces underlying disease outbreaks and reconcile viral evolution with epidemiological relationships that also capture pathogen population shifts.

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.

Native and Non-Native Bemisia tabaci NAFME Haplotypes Can Be Implicated in Dispersal of Endemic and Introduced Begomoviruses in Oman

Irrigated agriculture and global trade expansion have facilitated diversification and spread of begomoviruses (Geminiviridae), transmitted by the Bemisia tabaci (Gennadius) cryptic species. Oman is situated on major crossroads between Africa and South Asia, where endemic/native and introduced/exotic begomoviruses occur in agroecosystems. The B. tabaci 'B mitotype' belongs to the North Africa-Middle East (NAFME) cryptic species, comprising at least eight endemic haplotypes, of which haplotypes 6 and/or 8 are recognized invasives. Prevalence and associations among native and exotic begomoviruses and NAFME haplotypes in Oman were investigated. Nine begomoviral species were identified from B. tabaci infesting crop or wild plant species, with 67% and 33% representing native and exotic species, respectively. Haplotypes 2, 3, and 5 represented 31%, 3%, and 66% of the B. tabaci population, respectively. Logistic regression and correspondence analyses predicted 'strong'- and 'close' virus-vector associations involving haplotypes 5 and 2 and the exotic chili leaf curl virus (ChiLCV) and endemic tomato yellow leaf curl virus-OM, respectively. Patterns favor a hypothesis of relaxed virus-vector specificity between an endemic haplotype and the introduced ChiLCV, whereas the endemic co-evolved TYLCV-OM and haplotype 2 virus-vector relationship was reinforced. Thus, in Oman, at least one native haplotype can facilitate the spread of endemic and introduced begomoviruses.

Genetic diversity of cryptic species of Bemisia tabaci in Asia

Bemisia tabaci is a species complex consisting of various genetically different cryptic species worldwide. To understand the genetic characteristics and geographic distribution of cryptic species of B. tabaci in Asia, we conducted an extensive collection of B. tabaci samples in ten Asian countries (Bangladesh, Indonesia, Japan, Korea, Myanmar, Nepal, Philippines, Singapore, Taiwan, and Vietnam) from 2013 to 2020 and determined 56 different partial sequences of mitochondrial cytochrome oxidase subunit I (COI) DNA. In addition, information on 129 COI sequences of B. tabaci identified from 16 Asian countries was downloaded from the GenBank database. Among the total 185 COI sequences of B. tabaci, the sequence variation reached to 19.68%. In addition, there were 31 cryptic species updated from 16 countries in Asia, that is, Asia I, Asia I India, Asia II (1-13), Asia III, Asia IV, Asia V, China 1-6, MEAM (1, 2, K), MED, Australia/Indonesia, Japan (1 and 2). Further, MED cryptic species consisted of 2 clades, Q1 and Q2. This study provides updated information to understand the genetic variation and geographic diversity of B. tabaci in Asia.

Hsp90 Gene Is Required for Mi-1-Mediated Resistance of Tomato to the Whitefly Bemisia tabaci

The Mi-1 gene of tomato (Solanum lycopersicum) confers resistance against some nematodes and insects, but the resistance mechanisms differ depending on the harmful organism, as a hypersensitive reaction (HR) occurs only in the case of nematodes. The gene Rme1 is required for Mi-1-mediated resistance to nematodes, aphids, and whiteflies, and several additional proteins also play a role in this resistance. Among them, the involvement of the chaperone HSP90 has been demonstrated in Mi-1-mediated resistance for aphids and nematodes, but not for whiteflies. In this work, we studied the implication of the Hsp90 gene in the Mi-1 resistance against the whitefly Bemisia tabaci by means of Tobacco rattle virus (TRV)-based virus-induced gene silencing (VIGS). The silencing of the Hsp90 gene in tomato Motelle plants carrying the Mi-1 gene resulted in a decrease in resistance to whiteflies, as oviposition values were significantly higher than those on non-silenced plants. This decrease in resistance was equivalent to that caused by the silencing of the Mi-1 gene itself. Infiltration with the control TRV vector did not alter Mi-1 mediated resistance to B. tabaci. Similar to the Mi-1 gene, silencing of Hsp90-1 occurs partially, as silenced plants showed a significant but not complete suppression of gene expression. Thus, our results demonstrate the requirement of Hsp90 in the Mi-1-mediated resistance to B. tabaci and reinforce the hypothesis of a common model for this resistance to nematodes and insects.

Predicting the insecticide-driven mutations in a crop pest insect: Evidence for multiple polymorphisms of acetylcholinesterase gene with potential relevance for resistance to chemicals

The silverleaf whitefly Bemisia tabaci (Gennadius, 1889) (Homoptera: Aleyrodidae) is a serious invasive herbivorous insect pest worldwide. The excessive use of pesticides has progressively selected B. tabaci specimens, reducing the effectiveness of the treatments, and ultimately ending in the selection of pesticide-resistant strains. The management of this crop pest has thus become challenging owing to the level of resistance to all major classes of recommended insecticides. Here, we used in silico techniques for detecting sequence polymorphisms in ace1 gene from naturally occurring B. tabaci variants, and monitor the presence and frequency of the detected putative mutations from 30 populations of the silverleaf whitefly from Egypt and Pakistan. We found several point mutations in ace1-type acetylcholinesterase (ace1) in the studied B. tabaci variants naturally occurring in the field. By comparing ace1 sequence data from an organophosphate-susceptible and an organophosphate-resistant strains of B. tabaci to ace1 sequence data retrieved from GenBank for that species and to nucleotide polymorphisms from other arthropods, we identified novel mutations that could potentially influence insecticide resistance. Homology modeling and molecular docking analyses were performed to determine if the mutation-induced changes in form 1 acetylcholinesterase (AChE1) structure could confer resistance to carbamate and organophosphate insecticides. Mutations had small effects on binding energy (ΔG_b) interactions between mutant AChE1 and insecticides; they altered the conformation of the peripheral anionic site of AChE1, and modified the enzyme surface, and these changes have potential effects on the target-site sensitivity. Altogether, the results from this study provide information on genic variants of B. tabaci ace1 for future monitoring insecticide resistance development and report a potential case of environmentally driven gene variations.

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.

Characterization of the Asian Citrus Psyllid-‘Candidatus Liberibacter Asiaticus’ Pathosystem in Saudi Arabia Reveals Two Predominant CLas Lineages and One Asian Citrus Psyllid Vector Haplotype

In Saudi Arabia (SA), the citrus greening disease is caused by ‘Candidatus Liberibacter asiaticus’ (CLas) transmitted by the Asian citrus psyllid (ACP) Diaphorina citri. The origin and route(s) of the ACP-CLas pathosystem invasion in SA have not been studied. Adult ACP were collected from citrus trees in SA and differentiated by analysis of the mitochondrial cytochrome oxidase I (mtCOI) and nuclear copper transporting protein (atox1) genes. A phylogenetic analysis of the Wolbachia spp. surface protein (wsp) gene was used to identify the ACP-associated Wolbachia spp. A phylogenetic analysis of the atox1 and mtCOI gene sequences revealed one predominant ACP haplotype most closely related to the Indian subcontinent founder populations. The detection and identification of CLas in citrus trees were carried out by polymerase chain reaction (PCR) amplification and sequencing of the 16S rDNA gene. The CLas-integrated prophage genomes were sequenced, annotated, and used to differentiate CLas populations. The ML and ASTRAL trees reconstructed with prophages type 1 and 2 genome sequences, separately and concatenated, resolved two major lineages, CLas-1 and -2. The CLas-1 clade, reported here for the first time, consisted of isolates from SA isolates and Pakistan. The CLas-2 sequences formed two groups, CLas-2-1 and -2-2, previously the ‘Asiatic’ and ‘Floridian’ strains, respectively. Members of CLas-2-1 originated from Southeast Asia, the USA, and other worldwide locations, while CLas-2-2 was identified only in Florida. This study provides the first snapshot into the status of the ACP-CLas pathosystem in SA. In addition, the results provide new insights into the pathosystem coevolution and global invasion histories of two ACP-CLas lineages with a predicted center of origin in South and Southeast Asia, respectively.

Earlier than expected introductions of the Bemisia tabaci B mitotype in Brazil reveal an unprecedented, rapid invasion history

During 1991, in Brazil, the presence of the exotic Bemisia tabaci B mitotype was reported in São Paulo state. However, the duration from the time of initial introduction to population upsurges is not known. To investigate whether the 1991 B mitotype outbreaks in Brazil originated in São Paulo or from migrating populations from neighboring introduction sites, country-wide field samples of B. tabaci archived from 1989-2005 collections were subjected to analysis of mitochondrial cytochrome oxidase I (mtCOI) and nuclear RNA-binding protein 15 (RP-15) sequences. The results of mtCOI sequence analysis identified all B. tabaci as the NAFME 8 haplotype of the B mitotype. Phylogenetic analyses of RP-15 sequences revealed that the B mitotype was likely a hybrid between a B type parent related to a haplotype Ethiopian endemism (NAFME 1-3), and an unidentified parent from the North Africa-Middle East (NAF-ME) region. Results provide the first evidence that this widely invasive B mitotype has evolved from a previously undocumented hybridization event. Samples from Rio de Janeiro (1989) and Ceará state (1990), respectively, are the earliest known B mitotype records in Brazil. A simulated migration for the 1989 introduction predicted a dispersal rate of 200-500 km/year, indicating that the population was unlikely to have reached Ceará by 1990. Results implicated two independent introductions of the B mitotype in Brazil in 1989 and 1990, that together were predicted to have contributed to the complete invasion of Brazil in only 30 generations.