Bemisia tabaci (Hemiptera: Aleyrodidae) is a haplo-diploid species of a sap-feeding insect belonging to the group of insects commonly known as whiteflies. From earlier analyses of mitochondrial and ribosomal markers, the Asia-Pacific region was construed to have three major indigenous races and a large collection of genotypes with no clear association with any race. However, recent research using microsatellites has demonstrated that B. tabaci from the Asia-Pacific region may be split into six genetic populations with little or no gene flow between them. These genetic populations bared only superficial similarity to the mitochondrial and ribosomal-defined races. Further, four of the six can be further split into two subpopulations, which again showed little evidence of gene flow between them. While the patterns reflected a strong geographic structure, physical barriers alone could not explain all the observed structure. Differential host plant utilization explained some of the substructure, but could not explain the overall structure. The roles of mating interference and Wolbachia in developing the genetic structure were considered. The lack of gene flow between genetic populations and some subpopulations further suggested that the barriers were either sufficiently impermeable to immigration or that reproductive isolation and competitive interactions were sufficiently strong to prevent gene flow. If the latter was the case, there might be as many as 10 morphologically indistinguishable species indigenous to the Asia-Pacific region.
Key words: Bemisia tabaci, Asia-Pacific region, whiteflies, mitochondrial and ribosomal markers
Bemisia tabaci (Hemiptera: Aleyrodidae) is a haplo-diploid species of a sap-feeding insect belonging to the group of insects commonly known as whiteflies. For field crops, it is a pest in warm to hot climates between 30° north and south of the equator, where it can complete a generation in two to three weeks. It prefers short-lived herbaceous hosts including numerous dicotyledonous agricultural and horticultural species. In these systems, crops usually mature in less than six months. A typical cropping regime involves successive overlapping planting. After harvest the crop is destroyed, necessitating dispersal to younger crops or weeds. During the life of a crop, infestation involves an initial colonization phase, followed by local, within crop dispersal and population increase, and then a final emigration phase as the crop finishes (for review, see De Barro 1995).
Bemisia tabaci first drew major attention in the mid-1980s when severe crop losses were observed in Israel and the southern USA (for review, see De Barro 1995). Subsequently, major outbreaks and invasions have occurred on all continents except Antarctica. The B. tabaci responsible for the Israel and USA outbreaks demonstrated high levels of insecticide resistance, transmitted previously unknown begomoviruses, and induced a physiological change in squash known as "squash silverleafing" and was found to exhibit a unique electromorph pattern for general esterases and was designated as the B-biotype. The indigenous, American B. tabaci exhibited a different esterase pattern and was named the A-biotype (Costa & Brown 1991; Perring et al. 1992, Burban et al. 1992, Costa et al. 1993). Since then, genetic variation within this species has been subject to numerous studies using allozymes especially esterases (Coats et al. 1994, Brown et al. 1995a, Gunning et al. 1995), RAPD PCR (Gawel & Bartlett 1993; De Barro & Driver 1997), AFLP (Cervera et al. 2000), mitochondrial CO1 (Frohlich et al. 1999, Brown 2000, Viscarret et al. 2003, Berry et al. 2004), and ribosomal ITS1 (De Barro et al. 2000a; Abdullahi et al. 2003; Wu et al. 2003). All have shown considerable variation, with genotypes within particular geographic regions generally showing greater relatedness (Frohlich et al. 1999, De Barro et al. 2000a, Berry et al. 2004, De Barro et al. 2005a). As well as genetic variation, there is considerable phenotypic variability with regard to the relative ability of different genotypes to transmit begomoviruses (Bedford et al. 1994, Brown & Bird 1995), the rate of development (Wang & Tsai 1996), ability to utilize different hosts (Burban et al. 1992, Brown & Bird 1995, Bedford et al. 1994), and the ability to induce physiological changes in some hosts (Costa & Brown 1991, Cohen et al. 1992, Bedford et al. 1994). This has led to the subdivision of the species into a series of biotypes, of which at least 20 have so far been described (biotypes A through to T, see Perring 2001 for review). In reality, many of these biotypes are based primarily on genetic information with little or no biological data to underpin them and certainly no comprehensive comparison across the different biotypes so as to determine the biotic factors that define them. Further, not only have the biological bounds for each biotype not been determined, the genetic bounds have similarly been ignored; as a consequence it is not known where one biotype begins and another ends.
The processes that operate to create the observed genetic structure have not been addressed to a large extent until very recently. Brown et al. (1995a) suggested that host plant associations might play an important role. This was based on Bird (1957) which noted in Puerto Rico the existence of the narrowly specific "Jatropha" and the more polyphagous "Sida" races of B. tabaci while Burban et al. (1992) and later, Abdullahi et al. (2003) and Berry et al. (2004) observed in Ivory Coast and elsewhere in sub-Saharan Africa a race considered largely confined to cassava. These led to the development of the host-races concept (Bird 1957, Brown et al. 1995a), but apart from these examples, there are no other well-documented examples of such narrow host specificity. De Barro et al. (2005a) reanalyzed the available genetic data drawn from Frohlich et al. (1999), De Barro et al. (2000), Abdullahi et al. (2003), Viscarret et al. (2003), and Berry et al. (2004), and concluded that while evidence for widespread host-based genetic structuring was scant, there was considerable evidence for allopatric divergence resulting from geographic separation.
Another question is whether these races consist of a genetically homogenous mix or whether there is substructure that to some extent reflects the disturbed and highly heterogenous agricultural landscapes that these insects occupy. Our current understanding of genetic relatedness in B. tabaci has, until very recently, been based on the mitochondrial CO1 and ribosomal ITS1 gene regions. Based on mitochondrial CO1, the genetically distinct races identified in De Barro et al. (2005a) diverge from each other by 15-22% (Brown et al. 1995a, Frohlich et al. 1999). Using a "ball park" estimate of mutation rates of 2-5% per 1 million years for mitochondrial DNA, suggests a common origin well before the evolution of either agricultural (~10,000 y bp) or horticultural (~5,000 y bp) systems (Wells 2002). As such, it is unknown whether these gene regions provide sufficient resolution to identify substructure that may be useful from a management perspective. An understanding of how B. tabaci is structured may also offer insights into the relationship between population structure on a regional scale and management in agricultural production systems. Through this, the relative importance of various components of the system in terms of factors such as sources of re-infestation and the flow of insecticide resistance, could be assessed.
De Barro et al. (2005a) in their reanalysis of the phylogenetic relationships between genotypes of B. tabaci concluded that the association of races with particular geographic regions, the lack of evidence for interbreeding between representatives of the different races to produce fertile female offspring, and the lack of evidence for generalized host-based segregation suggested that allopatric rather than sympatric divergence was the process that best explained the observed genetic structure. On the basis of their analyses they concluded that there were six major races, Asia, Bali, Australia, sub-Saharan Africa, Mediterranean/Asia Minor/Africa, and New World as well as a large collection of genotypes from the Asian region with no clear association with any of the races.
To try and resolve the uncertainty of association within the cluster of Asian genotypes where the genetic relatedness remained unresolved, microsatellites (Fig. 1-A) (De Barro et al. 2003, De Barro 2005 in press) have recently been developed and applied. The resulting structures for B. tabaci from the Asia-Pacific region not surprisingly varied from that obtained using both ribosomal ITS1 and CO1 (De Barro et al. 2000a, De Barro et al. 2005a). This is to be expected given that not only did microsatellites uncover layers of structure that had previously remained hidden, they also enabled the relationships of the large number of "unresolved" genotypes identified in De Barro et al. (2005a) to be determined. The relationships suggest that instead of three races and an unresolved cluster of Asian genotypes indigenous to the Asia-Pacific Region, there are six genetic populations with minimal, if any, gene flow between them. In particular, the Asia race is subdivided across four genetic populations; the Australia race splits into two genetic populations, one Australian and the other Indonesian; while Bali has genotypes from central Java added to it. The relationship between the six genetic populations has a marked geographical association; the key exceptions, Nauru and FSM, are incursions resulting from the recent movement of live plant material (PJ De Barro unpublished data).
The relationship between samples within genetic populations (Fig. 1-B) has a similarly strong geographical association in several of the groups. Here, Group 1 segregated into a mainland Asian subpopulation and an Indonesian subpopulation; Group 2 split into eastern and western Asian subpopulations; and Group 4 split into a central Java subpopulation and an east Java/Bali subpopulation. Group 5, the Australian genetic population, has a less marked degree of separation between the northern and southern Australian subpopulations as there is a high level of migration and gene flow between the two. The remaining genetic populations either showed substructure not related to geographic separation as in the case of Group 3 or showed no evidence of substructure at all as in Group 6.
For each of the genetic populations, there is little evidence of gene flow between populations. Only two individuals from Group 1 and five from Group 2 had alleles that were likely to have been derived from a different genetic population, suggesting there was very little genetic exchange between genetic populations. This was similarly reflected by the low level of immigration (seven individuals all found in Bangladesh4) seen in Group 2 and its absence in all other genetic populations. Therefore, with the exception of Bangladesh4, individuals from any given sample location were represented by a single genetic population rather than by a mix. The question though is why mixed infestations are uncommon and how this strong geographical structure is maintained. Physical barriers that prevent movement are one possible answer. These would explain why Australia is a unique group and why mainland Asia is, for the most part, distinct from the Indonesian populations. However, it becomes increasingly difficult to invoke an argument for physical barriers to explain the complex structure that appears in mainland Asia or in Java where four different genetic populations occur in apparently unmixed states. An alternative explanation lies in the competitive interactions between genetic populations which effectively drive one to extinction. Reitz & Trumble (2002) in their review indicate that individuals belonging to one race have the capacity to suppress those from another through mating interference. Suppression of one genetic group by another has been described previously and has been suggested as the cause for one race of B. tabaci overwhelming another (Perring et al. 1994, De Barro & Hart 2000, Pascual & Callejas 2004, De Barro et al. 2005b, Zang et al. 2005). In the interaction between the B and A biotypes, Perring et al. (1994) found that male B biotypes blocked courtship by the A biotype. Similarly, mating interference was also observed by De Barro & Hart (2000) and Pascual & Callejas (2004) in interactions between the AN and B, and B and Q biotypes, respectively. Further, De Barro et al. (2005b) determined that this interaction had a density-dependent relationship where at densities greater than 20:1 (AN biotype: B biotype) the B biotype was unable to establish unless it was provided with a host plant mix which contained one host that only it was able to utilize. Therefore, if one assumes that the host range of each of the various genetic populations in Asia is fairly equivalent, and there is no data to suggest otherwise, then it is possible that density-dependent competitive interactions will operate to exclude an invader belonging to a different genetic population. What is not clear from this study is at what scale competitive interactions operate and at what length of time is taken for competitive interactions to result in local extinction.
The grouping of Indonesian and mainland Asian genotypes into Group 1 is a little unexpected given the overall genetic separation observed between mainland Asian and Indonesian samples. Collections spanning the period 1997_2003 suggest that the appearance of these genotypes in Indonesia is recent as none were detected prior to 1999 despite extensive surveys (PJ De Barro, unpublished data). What is particularly interesting about this observation is that the detection of a mainland Asian B. tabaci in Indonesia coincided with the first detections of a begomovirus disease of chilli known locally as "Bulai Amerika" (Sulandari et al. 2001, Rahayu 2004). While circumstantial, the simultaneous first detections of the two suggest a strong possibility of a recent incursion from mainland Asia into Indonesia. Further, the observation that the whitefly reproduces on a host that does not appear to be utilized by other Indonesian B. tabaci suggests that this may be a fulfillment of the prerequisite for establishment predicted in De Barro et al. (2005b).
While geographical distance between populations was in most cases the key element behind the observed patterns, there was evidence, be it limited, that host plant choice was an additional factor influencing genetic structure. In the Hemiptera, host-based genetic structuring of populations has been observed in several species of aphids including Sitobion avenae (Sunnucks et al. 1997), Acyrthosiphon pisum (Simon et al. 2003) and Aphis gossypii (Vanlerberghe-Masutti & Chavigy 1998) as well as in the diaspidid scale Aspidiotus cryptomeriae (Miyanoshita & Tatsuki 2001) and the membracid leaf hopper, Enchenopa binotata (Rodriguez et al. 2004).
In this study, host plant plays no apparent role in defining the limits between the six genetic populations identified. However, there are two instances where the host plant can be considered to be contributing to subpopulation genetic structure. In the case of Group 3, a genetic population consisting of individuals from Vietnam had two subpopulations that occurred within a few meters of each other, one associated with poinsettia, Euphorbia pulcherima, and the other with vegetables. Bemisia tabaci is known for its capacity to disperse up to distances of at least 7 km (Byrne 1999) so the distances between the poinsettia and the vegetables, a few meters, were well within the whitefly's capacity to disperse. While the capacity of individuals from the two subpopulations to utilize the two hosts would need to be tested experimentally, the observation does support a role for bidirectional, host-based isolation. The only other association with a host plant is seen in Group 1. Here, the Indonesian subpopulation is associated with billygoat weed, Ageratum conyzoides, irrespective of the distance between collection sites. However, whether geographic separation between mainland Asia and Indonesia or the host plant itself is driving the structure cannot be determined.
The lack of evidence for the host plant playing a major role in population genetic structure may help explain why the Bird (1957) "Jatropha" race and the Burban et al. (1992) and Abdullahi et al. (2003) cassavarace are not supported by more, rigorous examples. There are numerous claims for monophagy, but these appear to be based on the failure to make to the transition from "not known to be a host" to "known not to be a host." The example of the sub-Saharan Africa cassava race of Burban et al. (1992) and Abdullahi et al. (2003) is perhaps a case in point. The association of B. tabaci with cassava can only have arisen in the past few hundred years as cassava is exotic to this region having been introduced in the late 16th century from the New World. The whitefly itself is not derived from the New World as it is genetically distinct from the New World race (Frohlich et al. 1999, De Barro et al. 2000, De Barro et al. 2005a). In the first place, the levels of divergence for CO1 of up to 16.5% identified between the cassava and non-cassava populations suggest (Berry et al. 2004) the two diverged (based on 5% per million years) around three million years ago. Further, Berry et al. (2004) showed considerable evidence for allopatric genetic structure with the cassava-feeding B. tabaci from West Africa differing from those from southern Africa by at least 2%, and an Ivory Coast individual differing from the West African clade by approximately 6%. Finally, Berry et al. (2004) identified five distinct cassava-feeding genetic groups of sub-Saharan B. tabaci all separated by steps of between 2% and 4% divergence, differences that suggest separations of several 100,000 years. For a host plant that has been introduced only 400 years ago to be driving this structure would require simultaneous shifts from indigenous hosts to cassava in five genetically distinct groups in at least three separate geographic areas. Perhaps a more plausible explanation is that the majority of sub-Saharan B. tabaci have a capacity to utilize cassava as well as a range of other hosts that have not as yet been discovered and that factors other than host plant explain the observed structure. In conclusion, the evidence points overwhelmingly to B. tabaci having the ability to utilize multiple hosts, with many of these being common to most B. tabaci genotypes.
The association of genetic populations with particular geographic regions suggests isolation is the key mechanism driving the range of genetic diversity observed in B. tabaci. The explicit testing of the capacity of members of the various ITS1 races or microsatellite groupings from this study to interbreed has not been undertaken. However, studies involving B. tabaci from races not included in this study suggest a very limited capacity to interbreed (Liu et al. 1992, Costa et al. 1993, Bedford et al. 1994, Byrne et al. 1995, Ronda et al. 1999, De Barro & Hart 2000, Moya et al. 2001, Maruthi et al. 2001, Pascual & Callejas 2004, Liansheng et al. 2005, Zhang et al. 2005). This, coupled with the lack of evidence for gene flow and the geographic separation of race, point toward allopatric divergence as the process that best explains the observed relationships. In contrast, evidence for sympatric divergence was weak and there was certainly no compelling basis for assuming segregation based on differential host utilization, which was an important factor for determining population structure.
There has been a long-running debate as to whether B. tabaci is a single species or a complex of closely related taxa. Frohlich et al. (1999) concluded that B. tabaci was "species complex" although it was not clear from the study as to how this conclusion was reached. Since then, the term species complex has filtered into the Bemisia literature. The first comprehensive assessment of the species complex idea was undertaken by De Barro et al. (2005a) who concluded that there was insufficient data to support the step of raising the races identified using ITS1 and CO1 to species status as in many cases, the molecular distances among the races were trivial even though in some cases the ability to interbreed had been lost. This study however demonstrates that while the genetic differences using ITS1 may be considered small, the evidence for gene flow between genetic populations is scant. While only genotypes considered indigenous to the Asia-Pacific Region were included in this study, the microsatellite data indicated that they may be segregated into six reproductively isolated genetic populations. On this basis, one can argue that the three Asia-Pacific races, Asia, Australia and Bali, plus the unresolved cluster of Asian genotypes identified using ITS1 represent six distinct species (De Barro et al. 2005). The lack of gene flow between subpopulations further suggests that the barriers were either sufficiently impermeable to immigration or that reproductive isolation and competitive interactions were sufficiently strong to prevent gene flow between subpopulations. If the latter is the case, it suggests that there may be as many as 10 morphologically indistinguishable species indigenous to the Asia-Pacific Region.
Bemisia tabaci is well known for its capacity to develop resistance against insecticides (De Barro 1995). As to why this species seems so adept in this regard is unknown, however, the marked deficiency of heterozygotes observed in (De Barro 2005 in press) may provide some insight. The lack of heterozygotes cannot be explained by the differing numbers of chromosomes between males (haploid) and females (diploid) as only females were used in the study. Heterozygote deficits can be caused by null alleles (Callen et al. 1993); however, as indicated earlier, these are not thought to be common or pervasive. Wahlund effects are another source of deficiencies. These occur when populations are subdivided unequally with regard to allele frequencies so that random matings only occur within a portion of the population. The deficiency becomes apparent when the various subpopulations are sampled as a single unit. While we can not discount this possibility entirely, our sampling regime, the relative mobility of the adults, and the capacity of females to mate repeatedly would make this unlikely. Heterozygote deficiencies can also be attributable to infections with sex-ratio-distorting bacteria such as Wolbachia. Wolbachia have been found in B. tabaci belonging to the sub-Saharan African and New World races (Zchori-Fein & Brown 2002) as well as from some Old World populations in Asia and the Mediterranean (Nirgianaki et al. 2003). A screen of the collections that made up this study (data not shown) found that the collections belonging to each group other than those from Group 2 were positive for Wolbachia. All were infected with the same strain of Wolbachia belonging to the Con/Rug subgroup (Nirgianaki et al. 2003) and therefore, it was unlikely that Wolbachia was responsible for the structure observed in this study. A more likely explanation relates to the ecology of the pest which, in turn, may provide insight into the capacity of the species to rapidly develop insecticide resistance. Bemisia tabaci is a colonizer of short-lived hosts and if relatively few genotypes colonize any given crop, the chances of inbreeding are high and would lead to the observed deficiencies. A high level of inbreeding in combination with a strongly r selected species in which males are haploid would create an optimal environment for the rapid selection of resistant genotypes.
This broad genetic survey of B. tabaci suggests that the genetic populations are genetically heterogenous with respect to distance and, to a much lesser extent, to the host plant. The lack of gene flow between genetic populations and subpopulations supports the view that these populations are reproductively isolated and when combined with data on mating incompatibility suggests that B. tabaci is a species complex containing numerous sibling species. The data suggest that an understanding of how different B. tabaci genotypes are distributed across the landscape could be a very useful tool in the interpretation of the pest's ecology at farm and regional scales.
Figure 1 Genetic Groups (a, Groups 1_6, Left to Right) and Genetic Sub-Groups (B) of Bemisia Tabaci Identified Using Microsatellites. the Its1 Groupings Appear at the Top of the Figure.
Download the PDF. of this document, 191,308 bytes (187 KB).