The tobacco whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae), also known as the cotton, sweet potato, or silverleaf whitefly, is a pest of several plants including ornamental plants and agricultural crops grown around the world. Whiteflies damage crops through feeding and transmitting various plant pathogenic viruses (Liu et al., 2014). B. tabaci is currently known as a complex composed of several cryptic species. Morphologically indistinguishable populations of B. tabaci can exhibit variable molecular features and biological behavior. The host plant ranging from efficiency in plant virus transmission, ability to cause phytotoxicity, fecundity and insecticide resistance are traits that have been used to distinguish the so-called biotypes of B. tabaci (Barbosa et al., 2014).
The protein and DNA polymorphisms of B. tabaci can be combined with studies of biological characteristics to distinguish their biotypes by using experimental or technological approaches like electrophoresis of allozymes, analysis of randomly amplified polymorphic DNA (RAPDs) and nucleic acid sequence comparisons of nuclear or mitochondrial DNA markers (Calvert et al., 2001).
Nowadays, the more reproducible and informative method available to determine the genetic affiliation of a B. tabaci individual is based on molecular data, using markers such as Random Amplified Polymorphic DNA (RAPD) and Internal Transcribed Sequence (ITS) (Rabello et al., 2008), Sequence Characterized from Amplified Regions (SCAR) (Queiroz et al., 2016) and mitochondrial genes such as cytochrome oxidase 1 (De Barro, 2005; Marubayashi et al., 2012).
Molecular markers of B. tabaci have been recognized as a useful tool to understand the genetic dynamics of this species in the field and may correlate some molecular profiles with biological features of populations established in agricultural areas. However, since this tool has not yet been used on a large regional scale to evaluate local populations of B. tabaci on different crops, the objective of this study was to present the molecular profile and genetic variability of B. tabaci collected in several crops and localities in Brazil.
Obtaining DNA
Adults B. tabaci collected from crops in several locations in Brazil were identified according to morphological criteria (Gregory and Gregory, 2005) and maintained in 70% ethanol at -20°C. DNA for molecular analysis was obtained by macerating individual female adults in 60 μL extraction buffer (Tris-HCl 10 mM pH 8, EDTA 1 mM, Triton X-100 0.3%, proteinase K μg.mLâ»1). The mash was incubated for 30 min at 65°C followed by boiling for 10 min. The final homogenate was stored at -20°C until use.
Insect molecular identification
OPA-10 and OPA-13 oligonucleotides were used for molecular identification of whitefly biotypes and other whitefly species according to the methodology established by Lima et al. (2000). The default molecular weight profile of the fragments observed by RAPD was established for the B biotypes and the Brazilian native biotype (BR). Any sample whose fragment profile was different from those previously established was subjected to a new morphological identification to identify the biotype of Bemisia or other whitefly species.
Analysis of genetic variability of whitefly populations occurring in different crops in Brazil
For the purpose of analysis of genetic variability, five individuals in populations of B. tabaci were collected on crops in the US, but also on crops established in several locations in Brazil (Table 1), with the oligonucleotides OPA-02, OPA-03, OPA-05, OPA-10, OPA-11, OPA-13, OPA-15 and OPA-20.
A set of 17 samples (samples 74-93) (Table 1) from 10 Brazilian states and 10 different crops, and two samples from California, USA (cotton and tomato) were selected for analysis of genetic variability. Prior molecular identification had determined that 16 samples collected in Brazil and one sample from California belonged to biotype.
A sample identified as biotype A from California and a sample identified as the native biotype of B. tabaci occurring in Brazil were used as external samples to the group.
Comparison of the fragment profiles of the various whitefly biotypes was performed using five individuals of B. tabaci collected on various crops distributed in several Brazilian regions, and compared with biotypes coming from other regions of the world (Table 2) using the oligonucleotides OPA-04, OPA-10, OPA-11, OPA-13 and OPA-15.
RAPD reaction
Amplification reactions were carried out in 30 μL of a mixture containing 24.9 μL milliQ autoclaved water, 3.0 μL 10X buffer (60 mM Tris-HCl pH 8.8, 500 mM KCl and 20 mM MgCl2, Amersham, CA, USA), 1.2 μL of an oligonucleotide of random sequence (OPA-02, OPA-03, OPA-04, OPA-05, OPA-10, OPA-11, OPA-13, OPA-15 or OPA-20) (Operon Technologies, Inc.) at a concentration of 10 μM, 0.6 μL of 10 mM dNTP and 1 U of Taq DNA polymerase (Amersham, CA, USA) and 4 μL DNA (20 ng).
Genomic DNA amplification conditions
Amplifications were carried out in a thermal cycler (MJ Research PTC-100) having an initial denaturation step of 3 min at 94°C, programmed to 45 cycles of 1 min at 93°C denaturation, annealing for 1 min at 35°C and extension for 2 min at 72°C, and a final extension step of 5 min at 72 °C.
Visualization of DNA fragments by electrophoresis
The amplification products from the RAPD reactions were visualized on 1.5% agarose gel submerged in 1X TBE buffer (90 mM Tris-borate and 1 mM EDTA) for 3 h at 160 V, photographed and stored in the Eagle eye system (Stratagene).
In all gels, molecular weight markers (100 bp Ladder - Invitrogen) were used to determine molecular weight of the fragments amplified by the RAPD technique.
Statistical analysis
Five whiteflies were selected, at random, from the pool of individuals of each sample used in each experiment for molecular studies. The procedure described previously was repeated three times to confirm the data obtained.
The polymorphism among individuals in the population was determined from the fragments present in the electrophoretic gels. Then, a presence/absence matrix (1 or 0) was generated for each fragment. In case of doubt, the number 9 was used as standard. The similarity matrix was calculated using the Dice coefficient and UPGMA analysis produced a dendrogram which showed a grouping of individuals using the NTSYS program version 2.2 pc (Rohlf, 1993).
After that, the values obtained and tabulated in a spreadsheet were subjected to analysis of molecular variance (AMOVA) for statistical determination of the possible sources of variability found in populations by applying specific algorithms using Arlequin program version 2000 (Schneider et al., 2000).
Molecular identification
Insect samples from different crops originating from diverse locations had their molecular profiles determined using the OPA-13 oligonucleotide according to Lima et al. (2000). The fragment patterns produced by five individuals from each location were compared with the B. tabaci whitefly standard previously determined to be biotype B or the native biotype found in Brazil, referred as BR (Figure 1).
Using the oligonucleotide OPA-13 characteristic, RAPD fragment profiles for both biotype B and for the native biotype found in Brazil were identified and used as standard. For the biotype B, two fragments of 1000 and 1700 bp were found biotype B and for the native biotype (BR), 500, 1000 and 1500 bp fragments were found. The B. tabaci sample classification shown in Table 1 was made based on the standards for these biotypes. An RAPD fragment profile of 800 and 1400 bp was observed in two samples. After further investigation, it was observed that the molecular profile corresponded to Trialeurodes vaporariorum. This profile was identified in samples from Águas Mornas (SC) and Brasilia (DF).
This strategy allowed the analysis of 100 samples in a short period of time, noting that there was a predominance of individuals belonging to the B biotype of B. tabaci, found in 78% of samples, followed by the Brazilian native biotype (BR) of B. tabaci (4%). Other whiteflies were also identified using this methodology corresponding to 3% of the samples (1% Aleurodicus cocois and 2% T. vaporariorum). Some samples (11%), however, could not be identified using the established standards.
In addition, two samples previously identified as biotype Q of B. tabaci from Spain and Morocco had their electrophoretic profiles determined. The same was done with a sample from Nigeria containing the B. tabaci cassava biotype (Figure 2).
The Q biotype individuals from Spain and Morocco presented a characteristic RAPD fragment of 1350 bp. For the cassava biotype originated from Nigeria, there was a 950 bp fragment. Differentiation of the biotypes Q and cassava is based on the absence of the 1350 bp fragment in the latter. Regarding the established standards, it was observed that the whitefly individuals from the United States and previously identified as biotype B presented the same profile as the other whitefly individuals from Brazil, also identified as biotype B. This was confirmed by the presence of three DNA fragments in all samples (420, 850 and 1300 bp). Analyzing the profile of whitefly individuals previously identified as biotype A from the United States, a very similar pattern of DNA fragments was observed when compared to the native Brazilian biotype. The samples analyzed in this study show that, currently, the B biotype of B. tabaci predominates in all regions of Brazil. Previous identification results, described by Lima et al. (2002) showed that the Brazilian native biotype was present in several states such as Mato Grosso, Minas Gerais, Rio de Janeiro, Goiás, and the Federal District. However, this biotype currently appears to be restricted to certain regions and crops. These results demonstrate the high capacity for dispersion and invasion of the B biotype of B. tabaci in different regions of Brazil occurring in habitats previously occupied by the native Brazilian biotype. In addition, there apparently does not exist a host crop preference, indicating the extreme adaptability of this pest to different nutritional conditions.
The RAPD molecular profiles obtained in this study confirmed the descriptions previously set forth by Lima et al. (2002) for the native biotype, as well as permitting the monitoring of dispersion and occurrence of whitefly populations belonging to other species or biotypes.
The molecular strategy used can address the need to detect the presence of the B biotype and other biotypes in Brazilian agriculture from samples stored in alcohol relatively quickly and accurately in order to prepare plant protection reports.
Analysis of the genetic variability of whitefly population in Brazil
From the information generated by the RAPD marker profiles, a study of the genetic variability of populations from different crops in the Northeast, Midwest and Southeastern regions of Brazil, as well as from the United States, was carried out for comparison purposes. By submitting the DNA of individuals from ten whitefly populations to RAPD analysis with in decameric oligonucleotides, differences were observed
RAPD fragment profiles for the various samples (Table 2).
The ten RAPD oligonucleotides produced 103 loci and each sample produced 76.8±5.6 polymorphic fragments, on average. Regarding the number of polymorphic fragments, it was observed that sample 88 produced the highest number (82.7%) and the population represented by sample 86 produced the smallest number (63.5%). Each initiator generated 13±3.2 fragments per sample (Table 3).
The OPA-02 oligonucleotide generated an average of 15.4±1.8 polymorphic RAPD fragments for the study of genetic variability of whitefly populations. The OPA-05 primer produced the lowest average of polymorphic fragments. From the RAPD marker analysis, it was observed that the selected oligonucleotides showed potential for use in the generation of molecular profiles of various samples of B. tabaci (Table 4).
Once more, the molecular analysis confirmed the potential use of the OPA-10 and OPA-13 oligo-nucleotides to identify biotypes of B. tabaci, confirming the information by Lima et al. (2000) on the application of these oligonucleotides to discriminate the B biotype of B. tabaci. Noteworthy is also the potential of the oligonucleotides OPA-05, OPA-11 and OPA-15 in the identification of B. tabaci biotype B. Then, the binary data generated by RAPD oligonucleotides using B. tabaci samples 74-93 were used to generate a dendrogram and for analyses of molecular variance (AMOVA). The analyses showed differences among the analyzed samples. In the generation of groupings of whitefly samples by UPGMA the formation of specific groups for each location and crop was observed (Figure 3).
The dendrogram analysis showed that biotype A of B. tabaci from the United States differed from the other biotypes as it showed around 8% of genetic similarity.
For the remaining populations previously identified by the molecular method as biotype B, separation was observed among them based on location and crop of collection, so that it is not possible to establish a pattern among these parameters. Regarding B. tabaci native Brazilian biotype, it was observed that this did not form a separate grouping from the B biotype. However, the population supposedly defined as a native biotype, showed 25% similarity with the B biotype. This biotype was collected mainly on cotton, suggesting that the biotype presents greater preference for that crop. In addition, the initial molecular identification using the OPA-10 oligonucleotide showed that the native Brazil biotype had a similar molecular profile to biotype A from the United States. However, subjecting these samples to an analysis with more oligonucleotides showed a different grouping compared to biotype A.
From the information obtained in the dendrogram, an analysis of molecular variance (AMOVA) was carried out in order to better understand the dynamics of these whitefly populations. Analyzing the differences among biotypes A, B and the native Brazilian biotype, it was observed that 22.99% of the variation originated between the biotypes, 76.78% between samples within the biotypes, and 0.22% within samples from the same biotype. Following AMOVA analysis with all the biotypes grouped, it was observed that 99.68% of the variation was between the samples and 0.32% within the samples.
The analysis of only samples initially identified as biotype B of B. tabaci indicated a value of 99.75% between the samples and 0.25% variation within the samples. Comparing the A biotype and the native Brazilian biotype, it was found that 25.67% of the variation was between the biotypes, 74.33% between samples within the group, and 0% within the biotypes samples. The comparison between the B biotypes and the native biotype found in Brazil indicates a variation of 16.79% between groups, 83.04% between samples within the group, and 0.17% within the samples.
Regarding the B-biotype samples, the AMOVA values indicate differences between the samples of this whitefly biotype. In other words, despite the molecular identification revealing molecular markers for this specific biotype, it was observed that the samples belonging to this biotype have distinct gene pools. This was confirmed due to the high variation found among the samples.
This phenomenon has also been described by Silva et al. (2009) using RAPD markers. The authors collected whiteflies on okra, bean and pepper crops in São Luís (MA), and from the 96 RAPD molecular markers generated, a dendrogram was constructed which showed that populations on okra, beans and peppers were grouped according to the host crops. These data corroborate the findings of this study. However, in this work it was possible to analyze a larger number of locations and describe the phenomenon of grouping of whitefly populations as a function of crop that may be a phenomenon independent of location, suggesting a pattern of dispersion and allocation of this biotype in the field. This phenomenon can be explained partly by the reproductive strategy of the species, which has parthenogenetic arrhenotokous cycles, generating haploid males. Thus, the data obtained in the dendrogram indicate that each population established on a given crop is composed of clones due to the asexual reproduction of B. tabaci. This observation is supported by the low variability within the samples, observed to be less than 1%.
Helmi (2011) used random primers to assess the genetic variability of B. tabaci occurring in six hosts, that is, eggplant, tomato, cotton, watermelon, squash and okra. The results obtained showed that the whitefly may have different genotypes on adaptations to certain host plant species in Egypt, suggesting host specificity amongst under study B. tabaci populations based on RAPD molecular markers. The clustering pattern observed in the dendrogram showed that distinct genotypes exist among populations collected within the Egyptian region. After these information the author concludes these differences may be influencing some biological features like as virus vectoring capabilities and insecticide susceptibilities for example.
Furthermore, the availability of diet and the location of establishment of the sample may exert some in-fluence on the ecological dynamics of these individuals. Bernays (1999) presents B. tabaci behavior parameters for the selection of the host plant. Experiments with B. tabaci adult females showed that when a mixture of host plants is found in an area, the insect performance is decreased when compared to the same experiment using only a plant known as the best host. In mixtures of host plants, whiteflies spend less time on the leaf surface, moving most of the time and more often to new feeding places, spending less time at a given nutrition point. The author suggests that B. tabaci eventually perceives alternative feeding locations and that the observed performance must be related to the insect’s information processing mechanisms and decision-making. The data pointed by the author suggest that the observed behavior is the result of a characteristic loss of attention when the insect is facing a complex mixture of sensory signals. This can result in distraction and consequent reduction in decision making by the insect.
Finally, the phenomenon observed is due probably to the fact that the clonal populations established in a given area when at a low migration frequency and under selection pressure caused by chemical agents can contribute to the selection of a profile of molecular markers that are found specifically in a given B. tabaci population.
Using molecular analysis, it was observed that gene flow between two biotypes of the same population is smaller than between populations of identical biotypes (Moya et al., 2001). The present work confirmed this finding, since it was possible to identify specific molecular markers for tested biotypes independently of the sample origin. Thus, molecular techniques based on DNA proved useful to confirm the existence of biotypes and the degree of separation between them. For epidemiological and agro-economic purposes, it is essential to understand the intimate relationship developed between B. tabaci local populations and especially any differentiation among populations that is related to the host (Burban et al., 1992).
This information may serve as a basis for studies aiming to identify which factors influence the population distribution in a given area, but also, it may be useful in determining whitefly chemical control methods in relation to a specific host crop. It becomes very important to obtain information with respect to the sources of variation occurring within the B biotype of B. tabaci, since resistant populations of whitefly can be quickly selected after chemical insecticide applications.