Genome-wide analyses of the Bemisia tabaci species complex reveal contrasting patterns of admixture and complex demographic histories

Functional response of Franklinothrips vespiformis (Thysanoptera: Aeolothripidae) to eggs and nymphs of Bemisia tabaci (Hemiptera: Aleyrodidae)

The Middle East Asia Minor 1 biotype of Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) is a greenhouse and field crop pest of global significance. The objective of this study was to assess the potential of the generalist predatory thrips, Franklinothrips vespiformis Crawford (Thysanoptera: Aeolothripidae), as a biological control agent for B. tabaci. This was achieved by determining the functional responses of F. vespiformis larvae and adults to the egg and nymphal stages of B. tabaci under laboratory conditions. Analyses consisted of 10 replicates of each predator and prey stage combination on bean leaf discs for a 24-h period. Following logistic regression analyses to determine the functional response type exhibited, response parameters were estimated with nonlinear least squares regression using Roger's equation. Results showed that F. vespiformis larvae and adults exhibited a Type II functional response when feeding on immature B. tabaci. The handling times (Th) of F. vespiformis larvae and adults were magnitudes higher for B. tabaci nymphs than they were for eggs, which were in part driven by the higher attack rates (a) observed on eggs. The maximum attack rate (T/Th) for B. tabaci eggs and nymphs exhibited by first-stage larvae, second-stage larvae, and adult F. vespiformis increased with increasing predator age. Results from this study suggest that F. vespiformis larvae and particularly adults are promising biological control agents for B. tabaci and are efficient predators at both low and high prey densities.

A temperature-dependent phenology model for Bemisia tabaci MEAM1 (Hemiptera: Aleyrodidae)

The sweetpotato whitefly, Bemisia tabaci (Gennadius) Middle East-Asia Minor 1 (MEAM1), is widespread across tropical and subtropical regions, affecting hundreds of cultivated and wild plant species. Because the species transmits a variety of viruses, the whitefly has become one of the most economically significant insect pests in the world. Determining a pest's population growth potential as a function of temperature is critical for understanding a species population dynamics, predicting the potential range of the species and its associated diseases, and designing adaptive pest management strategies. The life history of B. tabaci MEAM1 was studied in life-table experiments at 7 constant temperatures ranging from 12 to 35 °C. Nonlinear equations were fitted to development, mortality, and reproduction data and compiled into an overall phenology rate-summation model using Insect Life Cycle Modeling (ILCYM) software, to simulate life-table parameters based on temperature. Life tables of B. tabaci MEAM1 observed at naturally variable temperature in La Molina, Lima, during different seasons, covering the entire temperature range of the species' predicted performance curve, were used to validate the model. Simulations predicted population growth within temperature between 13.9 and 33.4 °C, revealing a maximum finite rate of population increase (λ = 1.163), with a generation time of 33.3 days at 26.4 °C. Predicted species performance agreed well when compared against observed life tables and published data. The process-based physiological model presented here for B. tabaci MEAM1 should prove useful to predict the potential spatial distribution of the species based on temperature and to adjust pest control measures taking different population growth potentials due to prevailing temperature regimes into account.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Complex multiple introductions drive fall armyworm invasions into Asia and Australia

The fall armyworm (FAW) Spodoptera frugiperda is thought to have undergone a rapid 'west-to-east' spread since 2016 when it was first identified in western Africa. Between 2018 and 2020, it was recorded from South Asia (SA), Southeast Asia (SEA), East Asia (EA), and Pacific/Australia (PA). Population genomic analyses enabled the understanding of pathways, population sources, and gene flow in this notorious agricultural pest species. Using neutral single nucleotide polymorphic (SNP) DNA markers, we detected genome introgression that suggested most populations in this study were overwhelmingly C- and R-strain hybrids (n = 252/262). SNP and mitochondrial DNA markers identified multiple introductions that were most parsimoniously explained by anthropogenic-assisted spread, i.e., associated with international trade of live/fresh plants and plant products, and involved 'bridgehead populations' in countries to enable successful pest establishment in neighbouring countries. Distinct population genomic signatures between Myanmar and China do not support the 'African origin spread' nor the 'Myanmar source population to China' hypotheses. Significant genetic differentiation between populations from different Australian states supported multiple pathways involving distinct SEA populations. Our study identified Asia as a biosecurity hotspot and a FAW genetic melting pot, and demonstrated the use of genome analysis to disentangle preventable human-assisted pest introductions from unpreventable natural pest spread.

Unveiling biogeographic patterns in the worldwide distributed Ceratitis capitata (medfly) using populations genomics and microbiome composition

Invasive species are among the most important, growing threats to food security and agricultural systems. The Mediterranean medfly, Ceratitis capitata, is one of the most damaging representatives of a group of rapidly expanding species in the Tephritidae family, due to their wide host range and high invasiveness potential. Here, we used restriction site-associated DNA sequencing (RADseq) to investigate the population genomic structure and phylogeographic history of medflies collected from six sampling sites, including Africa (South Africa), the Mediterranean (Spain, Greece), Latin America (Guatemala, Brazil) and Australia. A total of 1,907 single nucleotide polymorphisms (SNPs) were used to identify two genetic clusters separating native and introduced ranges, consistent with previous findings. In the introduced range, all individuals were assigned to one genetic cluster except for those in Brazil, which showed introgression of an additional genetic cluster that also appeared in South Africa, and which could not be previously identified using microsatellite markers. Moreover, we assessed the microbial composition variations in medfly populations from selected sampling sites using amplicon sequencing of the 16S ribosomal RNA (V4 region). Microbiome composition and structure were highly similar across geographic regions and host plants, and only the Brazilian specimens showed increased diversity levels and a unique composition of its microbiome compared to other sampling sites. The unique SNP patterns and microbiome features in the Brazilian specimens could point to a direct migration route from Africa with subsequent adaptation of the microbiota to the specific conditions present in Brazil. These findings significantly improve our understanding of the evolutionary history of the global medfly invasions and their adaptation to newly colonised environments.

The Adipokinetic Peptides of Hemiptera: Structure, Function, and Evolutionary Trends

The Hemiptera comprise the most species-rich order of the hemimetabolous insects. Members of a number of superfamilies, most notably especially the more basal ones such as white flies, psyllids and aphids, belong to the most destructive agricultural insects known worldwide. At the other end of the phylogenetic tree are hemipterans that are notorious medical pests (e.g. kissing bugs). Most of the hemipteran species are good flyers, and lipid oxidation plays a pivotal role to power the contraction of flight muscles and, in aquatic water bugs, also deliver the ATP for the extensive swimming action of the leg muscles. Mobilization of stored lipids (mostly triacylglycerols in the fat body) to circulating diacylglycerols in the hemolymph is regulated by a set of small neuropeptides, the adipokinetic hormones (AKHs). We searched the literature and publicly available databases of transcriptomes and genomes to present here AKH sequences from 191 hemipteran species. Only few of these peptides were sequenced via Edman degradation or mass spectrometry, and even fewer were characterized with molecular biology methods; thus, the majority of the AKHs we have identified by bioinformatics are merely predicted sequences at this stage. Nonetheless, a total of 42 AKH primary sequences are assigned to Hemiptera. About 50% of these structures occur also in other insect orders, while the remaining 50% are currently unique for Hemiptera. We find 9 novel AKHs not shown to be synthesized before in any insect. Most of the hemipteran AKHs are octapeptides (28) but there is an impressive number of decapeptides (12) compared to other speciose orders such as Diptera and Lepidoptera. We attempt to construct a hypothetical molecular peptide evolution of hemipteran AKHs and find quite a bit of overlapping with current phylogenetic ideas of the Hemiptera. Lastly, we discuss the possibility to use the sequence of the aphid AKH as lead peptide for the research into a peptide mimetic fulfilling criteria of a green insecticide.

Distribution of invasive versus native whitefly species and their pyrethroid knock-down resistance allele in a context of interspecific hybridization

The invasion success of a species in an agrosystem is greatly influenced by environmental factors such as the use of insecticides, by the intrinsic evolutionary capabilities of the species, and also by interactions with resident species. On the island of La Réunion, the successive invasions of MEAM1 and MED whitefly species over the last 20 years have not only led an increased use of insecticides, but have also challenged the resident IO species. To trace the evolution of the 3 species, and the distribution of the kdr mutation (resistance to pyrethroid) in the para-type voltage-gated sodium channel, we genotyped 41 populations (using neutral nuclear markers) and look at the prevalence of the kdr allele. MEAM1 was predominantly present in agrosystems showing quasi fixation of the resistant kdr allele whereas IO was mainly in natural environments and did not have any resistant allele. Hybridization between the two former species was detected in low frequency but has not led to introgression of resistant alleles in the resident species so far. MED showed a limited distribution in agrosystems but all individuals displayed a resistant allele. These highly contrasting patterns of distribution and resistant mutations between invasive and resident whitefly species are further discussed.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest’s global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

This population genomics study identifies the complex multiple introduction history of Spodoptera frugiperda (fall armyworm) from the Americas, into Africa and Asia. This provides new insight into the ‘east-to-west’ directionality of gene flow, and suggests ample genomic exchange at the nuclear level.

Dominance of Asia II 1 species of Bemisia tabaci in Pakistan and beyond

Globally, Whitefly (Bemisia tabaci) is one of the most important insect pests of crops that causes huge economical losses. The current study was designed to exclusively screen the B. tabaci species in the cotton field of Pakistan during 2017-2020 and have to conduct comparative analysis of B. tabaci species in Asia where Asia II 1 has been reported. A total of 5142 B. tabaci sequences of mitochondrial cytochrome oxidase 1 (mtCO1) from Asian countries were analyzed to determine the species and their distribution in the region. Our analysis over time and space showed that Asia II 1 has gradually dominated over Asia 1 in Punjab Province and over both Asia 1 and MEAM1 in Sindh Province. Asia has been divided into three regions i.e., South Asia (2524 sequences), Southeast Asia (757 sequences) and East Asia (1569 sequences) and dominance of different species of B. tabaci has been determined by calculating the relative percentage of each species. Interestingly, Asia II 1 has been found dominant in the neighboring region (northern zone) of India and also being dominant in its central zone. The dominance of Asia II 1 in Pakistan and northern India explains whitefly epidemic being reported in recent years.

A taxonomic outline of the Poecilimon affinis complex (Orthoptera) using the geometric morphometric approach

The genus Poecilimon contains 145 species, widely distributed in the Palaearctic, among which the Poecilimon ornatus group has the greatest diversity in the Balkans. Despite several revisions of the genus, the systematics of the species group, and in particular, of the taxa associated with the species Poecilimon affinis, is still unsolved. Due to morphological similarity, P. affinis with its subspecies, P. nonveilleri and P. pseudornatus form the Poecilimon affinis complex. The aim of this study is to test the hypotheses of an outlined species complex, namely the P. affinis complex, within the P. ornatus group using morphological data. Geometric analysis was conducted to explore variation in the structure of the male tegmen, ovipositor, male cercus, and male pronotum. The number of teeth and stridulatory file measurements provided additional information on morphological variation within the complex. A phylogenetic tree based on the cytochrome c oxidase subunit I gene (COI) was used for comparison with the morphological data. Canonical variate analysis showed that male tegmen and male cercus are good morphostructures to distinguish the taxa belonging to the P. affinis complex from other species in the P. ornatus group. This may confirm our assumption for the designation of the P. affinis complex. The results of the principal component analysis of stridulatory file measurements, molecular data, and CVA of the ovipositor suggest adding two additional species to the complex: P. ornatus and P. hoelzeli.

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 sequencing to detect mutations associated with resistance to insecticides and Bt proteins in Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a major pest native to the Americas that has recently invaded the Old World. Point mutations in the target-site proteins acetylcholinesterase-1 (ace-1), voltage-gated sodium channel (VGSC) and ryanodine receptor (RyR) have been identified in S. frugiperda as major resistance mechanisms to organophosphate, pyrethroid and diamide insecticides respectively. Mutations in the adenosine triphosphate-binding cassette transporter C2 gene (ABCC2) have also been identified to confer resistance to Cry1F protein. In this study, we applied a whole-genome sequencing (WGS) approach to identify point mutations in the target-site genes in 150 FAW individuals collected from China, Malawi, Uganda and Brazil. This approach revealed three amino acid substitutions (A201S, G227A and F290V) of S. frugiperda ace-1, which are known to be associated with organophosphate resistance. The Brazilian population had all three ace-1 point mutations and the 227A allele (mean frequency = 0.54) was the most common. Populations from China, Malawi and Uganda harbored two of the three ace-1 point mutations (A201S and F290V) with the 290V allele (0.47-0.58) as the dominant allele. Point mutations in VGSC (T929I, L932F and L1014F) and RyR (I4790M and G4946E) were not detected in any of the 150 individuals. A novel 12-bp insertion mutation in exon 15 of the ABCC2 gene was identified in some of the Brazilian individuals but absent in the invasive populations. Our results not only demonstrate robustness of the WGS-based genomic approach for detection of resistance mutations, but also provide insights for improvement of resistance management tactics in S. frugiperda.

South West and North Central Nigeria: Assessment of cassava mosaic disease and field status of African cassava mosaic virus and East African cassava mosaic virus

Cassava mosaic disease (CMD), caused by cassava mosaic begomoviruses (CMBs), is a major threat to cassava production in Nigeria. The predominant CMBs in Nigeria are African cassava mosaic virus (ACMV), East African cassava mosaic virus (EACMV) and East African cassava mosaic Cameroon virus (EACMCV), which are transmitted through infected stem cuttings and whitefly vectors. This study was conducted in 2015 and 2017 to assess the epidemiology of CMD and the current distribution of CMBs in cassava farms in South West (SW) and North Central (NC) Nigeria. A survey of cassava farms was undertaken, and samples representative of disease symptoms were collected and assessed using molecular techniques. A total of 184 and 328 cassava farms were sampled in 2015 and 2017, respectively. CMD incidence for both regions surveyed was 43.80 and 12.25% in 2015 and 2017, respectively. Fields in SW recorded a higher incidence rate in 2015 (SW: 45.11%, NC: 42.47%), while the reverse occurred in 2017 (SW: 10.90%, NC: 14.01%). Overall, the CMD incidence in Benue State (NC) was significantly higher than other locations surveyed in both years. CMD symptom severity and mean whitefly population were higher in SW Nigeria in the two survey years. ACMV was widespread across both zones, occurring in 79.1% (453/613) and 54.8% (386/704) of cassava leaf samples analysed in 2015 and 2017, respectively. EACMV was detected in only 6.0% (37/613) and 4.7% (33/704) of all cassava leaf samples analysed in 2015 and 2017, respectively. Overall, a higher proportion of infected samples were found in NC in both 2015 (NC: 85.2%, SW: 75.4%) and 2017 (NC: 73.6%, SW: 45.2%). Detection using strain-specific primers revealed that 97% of EACMV positive samples were indeed infected by the EACMCV strain of the virus. As previously reported, samples with mixed infections showed a higher symptom severity than samples with single ACMV or EACMV infections. This study provides an update to the distribution of CMBs in SW and NC Nigeria and will be useful in development of monitoring and management strategies for the disease in both regions.

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.

Emergence of Maruca vitrata as a Major Pest of Food Legumes and Evolution of Management Practices in Asia and Africa

Legume pod borer, Maruca vitrata, has emerged as a major pest on food legumes in Asia and Africa. It is an oligophagous pest, feeding on over 70 species in Fabaceae. We examine the species complex in Asia, Africa, Oceania, and the Americas, with an emphasis on molecular taxonomy. Studies on pheromone production and perception suggest the existence of pheromone polymorphism, especially in Asia and Africa. No Maruca-resistant varieties are available in the major food legumes including cowpea, pigeonpea, mungbean, and yard-long bean. Legume growers use chemical pesticides indiscriminately, leading to the development of pesticide resistance. However recent developments in habitat management, classical biocontrol with more efficient parasitoids, biopesticides, and judicious use of insecticides pave the way for sustainable management of M. vitrata, which can reduce the pesticide misuse. Active engagement of the private sector and policy makers can increase the adoption of integrated pest management approaches in food legumes.

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.

Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control

Cryptic taxa have often been observed in the form of host‐associated species that diverged as the result of adaptation to alternate host plants. Untangling cryptic diversity in species complexes that encompass invasive species is a mandatory task for pest management. Moreover, investigating the evolutionary history of a species complex may help to understand the drivers of their diversification. The mealybug Hypogeococcus pungens was believed to be a polyphagous species from South America and has been reported as a pest devastating native cacti in Puerto Rico, also threatening cactus diversity in the Caribbean and North America. There is neither certainty about the identity of the pest nor the source population from South America. Recent studies pointed to substantial genetic differentiation among local populations, suggesting that H. pungens is a species complex. In this study, we used a combination of genome‐wide SNPs and mtDNA variation to investigate species diversity within H. pungens sensu lato to establish host plant ranges of each one of the putative members of the complex, to evaluate whether the pattern of host plant association drove diversification in the species complex, and to determine the source population of the Puerto Rican cactus pest. Our results suggested that H. pungens comprises at least five different species, each one strongly associated with specific host plants. We also established that the Puerto Rican cactus pest derives from southeastern Brazilian mealybugs. This is an important achievement because it will help to design reliable strategies for biological control using natural enemies of the pest from its native range.

Arthropod Invasions Versus Soybean Production in Brazil: A Review

Soybean production in Brazil has been markedly affected by invasions of non-native arthropod species that feed on the crop, severely impacting biodiversity, food security, health, and economic development. Data on soybean production losses and increase in insecticide usage over the last two decades have not been explored in association with past invasion events, and the dynamics underlying the recent blitz of invasive species into Brazil remain largely unclear. We provide a review of arthropod invasions in the Brazilian soybean agroecosystem since 1990, indicating that the introductions of Bemisia tabaci (Gennadius) MEAM1 (Hemiptera: Aleyrodidae), Tetranychus urticae (Koch) (Acari: Tetranychidae), and Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) are likely correlated with periods of increase in insecticide usage for soybean production. Using these three cases as examples, we review factors that could lead to increased likelihood of future invasions by particular pests, outlining four possible criteria to evaluate potential invasiveness of non-native arthropods: likelihood of entry, likelihood of establishment, biological features of the species, and availability of control measures. Spodoptera litura (F.) (Lepidoptera: Noctuidae) and Aphis glycines (Matsumura) (Hemiptera: Sternorrhynca) are examples of highly damaging soybean pests, related to one or more of these factors, that could be introduced into Brazil over the next years and which could lead to problematic scenarios. Melanagromyza sojae (Zehnter) (Diptera: Agromyzidae) also meets these criteria and has successfully invaded and colonized Brazilian soybean fields in recent years. Our review identifies current issues within soybean pest management in Brazil and highlights the need to adopt management measures to offset future costs and minimize lost revenue.

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.

Understanding Admixture: Haplodiploidy to the Rescue

Hybridization has broad evolutionary consequences, from fueling or counteracting speciation to facilitating adaptation to novel environments. Hybridization and subsequent introgression appear widespread along the tree of life. However, our understanding of how distinct evolutionary forces shape admixed genomes and the fate of introgressed genetic variants remains scarce. Most admixture research in animals has focused on diploid organisms. We propose that haplodiploid organisms can help resolve open questions about the genomic consequences of hybridization in natural populations. The ploidy difference between haploid males and diploid females, the availability of genome-wide male haplotypes, and ongoing cases of admixture make haplodiploid organisms promising models to improve our knowledge with regards to the evolution of hybrid genomes.

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.

Whole genome sequencing of Asia II 1 species of whitefly reveals that genes involved in virus transmission and insecticide resistance have genetic variances between Asia II 1 and MEAM1 species

Background

Whiteflies (Bemisia tabaci) are phloem sap-sucking pests that because of their broad host range and ability to transmit viruses damage crop plants worldwide. B. tabaci are now known to be a complex of cryptic species that differ from each other in many characteristics such as mode of interaction with viruses, invasiveness, and resistance to insecticides. Asia II 1 is an indigenous species found on the Indian sub-continent and south-east Asia while the species named as Middle East Asia Minor 1 (MEAM1), likely originated from the Middle-East and has spread worldwide in recent decades. The purpose of this study is to find genomic differences between these two species.

Results

Sequencing of the nuclear genome of Asia II 1 with Illumina HiSeq and MiSeq generated 198.90 million reads that covers 88% of the reference genome. The sequence comparison with MEAM1 identified 2,327,972 SNPs and 202,479 INDELs. In Total, 1294 genes were detected with high impact variants. The functional analysis revealed that some of the genes are involved in virus transmission including 4 genes in Tomato yellow leaf curl virus (TYLCV) transmission, 96 in Tomato crinivirus (ToCV) transmission, and 14 genes in insecticide resistance.

Conclusions

These genetic differences between Asia II 1 and MEAM1 may underlie the major biological differences between the two species such as virus transmission, insecticide resistance, and range of host plants. The present study provides new genomic data and information resources for Asia II 1 that will not only contribute to the species delimitation of whitefly, but also help in conceiving future research studies to develop more targeted management strategies against whitefly.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5877-9) contains supplementary material, which is available to authorized users.

Mitochondrial DNA genomes of five major Helicoverpa pest species from the Old and New Worlds (Lepidoptera: Noctuidae)

Five species of noctuid moths, Helicoverpa armigera, H. punctigera, H. assulta, H. zea, and H. gelotopoeon, are major agricultural pests inhabiting various and often overlapping global distributions. Visual identification of these species requires a great deal of expertise and misidentification can have repercussions for pest management and agricultural biosecurity. Here, we report on the complete mitochondrial genomes of H. assulta assulta and H. assulta afra, H. gelotopoeon, H. punctigera, H. zea, and H. armigera armigera and H. armigera conferta' assembled from high-throughput sequencing data. This study significantly increases the mitogenome resources for these five agricultural pests with sequences assembled from across different continents, including an H. armigera individual collected from an invasive population in Brazil. We infer the phylogenetic relationships of these five Helicoverpa species based on the 13 mitochondrial DNA protein-coding genes (PCG's) and show that two publicly available mitogenomes of H. assulta (KP015198 and KR149448) have been misidentified or incorrectly assembled. We further consolidate existing PCR-RFLP methods to cover all five Helicoverpa pest species, providing an updated method that will contribute to species differentiation and to future monitoring efforts of Helicoverpa pest species across different continents. We discuss the value of Helicoverpa mitogenomes to assist with species identification in view of the context of the rapid spread of H. armigera in the New World. With this work, we provide the molecular resources necessary for future studies of the evolutionary history and ecology of these species.

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.

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.

An integrative approach to discovering cryptic species within the Bemisia tabaci whitefly species complex

Bemisia tabaci is a cryptic whitefly-species complex that includes some of the most damaging pests and plant-virus vectors of a diverse range of food and fibre crops worldwide. We combine experimental evidence of: (i) differences in reproductive compatibility, (ii) hybrid verification using a specific nuclear DNA marker and hybrid fertility confirmation and (iii) high-throughput sequencing-derived mitogenomes, to show that the "Mediterranean" (MED) B. tabaci comprises at least two distinct biological species; the globally invasive MED from the Mediterranean Basin and the "African silver-leafing" (ASL) from sub-Saharan Africa, which has no associated invasion records. We demonstrate that, contrary to its common name, the "ASL" does not induce squash silver-leafing symptoms and show that species delimitation based on the widely applied 3.5% partial mtCOI gene sequence divergence threshold produces discordant results, depending on the mtCOI region selected. Of the 292 published mtCOI sequences from MED/ASL groups, 158 (54%) are low quality and/or potential pseudogenes. We demonstrate fundamental deficiencies in delimiting cryptic B. tabaci species, based solely on partial sequences of a mitochondrial barcoding gene. We advocate an integrative approach to reveal the true species richness within cryptic species complexes, which is integral to the deployment of effective pest and disease management strategies.