Open Access
How to translate text using browser tools
1 December 2014 DNA Barcoding and Elucidation of Cryptic Diversity in Thrips (Thysanoptera)
K. B. Rebijith, R. Asokan, V. Krishna, H. H. Ranjitha, N. K. Krishna Kumar, V. V. Ramamurthy
Author Affiliations +

Accurate and timely identification of invasive insect pests underpins most biological endeavors ranging from biodiversity estimation to insect pest management. In this regard, identification of thrips, an invasive insect pest is important and challenging due to their complex life cycles, parthenogenetic mode of reproduction, sex and color morphs. In the recent years, DNA barcoding employing 5′ region of the mitochondrial Cytochrome Oxidase I (CO-I) gene has become a popular tool for species identification. In this study, we employed CO-I gene sequences for discriminating 151 species of thrips for the first time. Analyses of the intraspecific and intrageneric distances of the CO-I sequences ranged from 0.0 to 7.91% and 8.65% to 31.15% respectively. This study has revealed the existence of cryptic species in Thrips hawaiiensis (Morgan) (Thysanoptera: Thripidae) and Scirtothrips perseae Nakahara (Thysanoptera: Thripidae) for the first time, along with previously reported cryptic species such as Thrips palmi Karny (Thysanoptera: Thripidae), T. tabaci Lindeman, Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae), Scirtothrips dorsalis Hood. We are proposing, the feasibility of hosting an independent integrated taxonomy library for thrips and indicate that it can serve as an effective system for species identification, this approach could potentially play a key role in formulating effective insect pest management strategies.

Thrips (Thysanoptera) include major sap sucking insect pests, which limit crop productivity and nutritional security, by direct feeding or indirectly by transmitting plant pathogenic viruses (Rebijith et al. 2011). Among the known 6,000 species of Thrips (Thripidae) (Mound & Morris 2007) nearly 1% are pests of agricultural crops and the following 14 species are reported as vectors of Tospovirus: Thrips tabaci Lindeman, T. palmi Karny, T. setosus Moulton, Scirtothrips dorsalis Hood, Frankliniella occidentalis (Pergande), F. schultzei (Trybom), F. bispinosa (Morgan), F. cephalica (Crawford), F. fusca (Hinds), F. gemina Bagnall, F. intonsa (Trybom), F. zucchini Nakahara Monteiro, Ceratothripoides claratris (Shumsher) and Dictyothrips betae Uzel (Ullman et al. 1997; Mound 2005; Jones 2005; Pappu et al. 2009; Ciuffo et al. 2010; Hassani-Mehraban et al. 2010). Most of them are highly polyphagous, have overlapping host ranges, thus being collected together which makes their identification difficult, even though morphological identification keys and web-based identification systems are available. Difficulty in identification of thrips not only exists in the developmental stage, but also between polyphagous thrips species, e.g. Thrips flavus Schrank found to be morphologically very similar to Thrips palmi Karny (Glover et al. 2010). Morphological examination of Thrips to species level is restricted to adult specimens, as there are no adequate keys for identification of egg, larvae, or pupae (Kadirvel et al. 2013). Thrips are notorious for eliciting taxonomic problems mainly because of minute size, polymorphism (Murai et al. 2001), lack of solid morphological characters (Kadirvel et al. 2013), co-existence of different species on the same host plant, high intraspecific variations observed in thrips populations (Mound 2011) and need for taxonomic expertise (Brunner et al. 2004).

In order to implement an integrated pest management (IPM) strategy, a simple, accurate, general and easily applicable method is required to facilitate the identification of Thrips spp. Subsequent studies on thrips have demonstrated the utility of CO-I sequences being effective in identification of thrips species (Timm 2008; Glover et al. 2010) and also in unraveling cryptic species as in the case of F. occidentalis and Pseudophilothrips gandolfoi (Rugman-Jones et al. 2010; Mound et al. 2010).

Considering all the above mentioned factors, it is even more important to properly identify invasive quarantine pest species introduced at the ports-of-entry for early detection and risk analysis, where speed and accuracy are paramount (Glover et al. 2010). In this connection, Hebert et al. (2003a,b) proposed the concept of DNA barcoding, a powerful tool to identify all metazoan species employing a short standardized 658 bp fragment of the 5? end of the mitochondrial cytochrome oxidase-I (CO-I) gene. DNA barcoding can be employed as the most effective approach for molecular identification of species independent of life stages, sex and polymorphism (Asokan et al. 2011), discriminating cryptic species (Glover et al. 2010; Rebijith et al. 2013), biotypes (Shufran et al. 2000), haplotypes (Zhang et al. 2011) and host and geographic associated genetic differences (Brunner et al. 2004; Rebijith et al. 2011). DNA barcoding can also play an important role in insect pest management program, where both selection and timing of the management practices can be affected by polymorphism and host adaptation/ suitability (Rebijith et al. 2013; Brunner et al. 2004).

In this study, molecular data have been generated and acquired (from NCBI and BOLD) to investigate the use of CO-I DNA barcoding in exploring diversity of thrips and is the first attempt to provide some understanding of the relationship between different Thrips spp. Thus, the purpose of the present study was to discriminate 996 sequences (CO-I) representing 151 thrips species globally and to record the presence of cryptic species and host or geographic associated genetic forms among thrips, if any.

Materials and Methods

Taxon Sampling and Morphological Identification

All the specimens were collected and stored in 95% ethanol during 2008–2012, and kept at -20 °C until processed. In each case, adults were examined using various characters described by Mound et al. 1996; Moritz et al. 2004. Prior to molecular work, all the thrips species were morphologically identified. In total, 151 thrips specimens representing 8 species were used for CO-I sequencing in this study (Table 1) and other sequences were retrieved from BOLD and NCBIGenBank (see Supplementary Table 1). A summary of the current scientific classification of each thysanopteran species is given in Supplementary Table 2. These supplementary tables are available online in Florida Entomologist 97(4) (2014) at

DNA Isolation and Polymerase Chain Reaction (PCR)

A single specimen of each species of Thrips was digested overnight at 60 °C in lysis solution (10 mM Tris-HCl- pH-7.60, 20 mM NaCl, 100 mM Na2 EDTA- pH-8.0, 1% Sarkosyl and 0.1 mg/mL proteinase K). DNA was extracted from the supernatant employing a non-destructive method (Mound and Morris, 2007a); while at the same time voucher specimens were mounted on glass slides and deposited with the National Pusa Collection (NPC), Indian Agricultural Research Institute (IARI), New Delhi, India. Standard protocols were employed for Polymerase Chain Reaction (PCR), cloning and sequencing (Hajibabaei et al. 2006).

Table 1.

Analysed samples of thrips species from the current study with descriptions of the sampling locations, Genbank accession numbers, names of thrips species and host plants, dates and voucher specimen details.










Table 2.

Intraspecific K2P genetic divergence among thrips species.


PCR was performed in a thermal cycler (ABIApplied Biosystems, Veriti, USA) using the following cycling parameters; an initial denaturation step at 94 °C for 5 min followed by 35 cycles at 94 °C for 30 s, an annealing step at 48 °C for 45 s, an extension step at 72 °C for 45 s and a final extension step at 72 °C for 20 min using the CO-I specific primers: LCO-1490 ; 5′-GGT CAA CAA ATC ATA AAG ATA TTG G-3′ and HCO-2198; 5′- TAA ACT TCA GGG TGA CCA AAA AAT CA-3′ (Hebert et al. 2003a; Hebert et al. 2003b). The total reaction volume of 25 µL contained 10 µM of each primer, 2.5 µL of 10X buffer, 2.5 mM MgCl2, 0.25 mM of each dNTP and 1.0 U of Taq DNA polymerase (Fermentas Life Sciences, United Kingdom). The amplified products were resolved on 1.2% agarose gel, stained with ethidium bromide (10µg/mL) and visualized in a gel documentation system (UVP).

Sequencing and Sequence analyses

Each PCR product was purified using Gel extraction Kit (NucleospinR Extract II, Macherey Nagel, Germany), cloned by ligation into PTZ57R/T vector (Fermentas Life Sciences, UK) and used to transform competent Escherichia coli (DH5á) cells. Blue- White colony, selection was carried out and plasmids were isolated using GenJET plasmid MiniPrep kit (Fermentas Life Sciences, United Kingdom), according to the manufacturer's protocol from the overnight culture of positive clones cultured in Luria Broth. Sequencing was performed in an automated sequencer (ABI prismR 3730 XL DNA Analyzer; Applied Biosystems, USA) using M13 universal primers, in both forward and reverse directions. A homology search was done using NCBI-BLAST ( and sequence alignment was performed using BioEdit version (Hall et al. 1999). In this study, we analyzed all the available thrips CO-I sequences from NCBI-GenBank for understanding and determining the intraspecific variations in each species (Meyer et al. 2005). All the sequences generated from this study are available at BOLD ( and in NCBI-GenBank. The complete details such as Thrips species, host plants, date of collection and voucher specimen number are given in Supplementary Table 1 along with the NCBI-GenBank extracted sequences.

CO-I sequences were aligned using the Clustal W program integrated in BioEdit.7.0 (Hall et al. 1999). The sequences were further analyzed employing MEGA.5.0 (Kumar et al. 1993) to obtain conspecific and congeneric distances, whilst Neighbor-Joining trees were constructed employing the Kimura-2-parameter (K2P) distance model (Kimura et al. 1980; Saitou et al. 1987). Node support were assessed with 1000 bootstrap pseudoreplicates.


The dataset consisted of CO-I sequences from 996 individuals representing 151 species of Thrips (Thysanoptera: Insecta), of which 151 sequences were produced in this study (Supplementary Table. 1, Supplementary Table. 2, Table. 1 respectively). Sequence analysis revealed that, 459 characters were variable among which 406 was parsimony informative. Evidence of nuclear copies was not found in any of the sequences subjected for analyses, which was supported by the absence of stop codons and the base composition was similar with no indels (Rebijith et al. 2012). Majority of the nucleotide substitutions occurred in the wobble position (third position) of the triplet codon (48.9%, 218 sites). Nucleotide frequencies were 29.46% (A), 39.38% (T), 13.87% (G) and 17.29% (C) and base composition found to be biased towards Adenine and Thymine, which together constituted 69.2% as is typical for other invertebrate genes (Wang et al. 2011; Rebijith et al. 2013).

Neighbor Joining Analysis

The CO-I dataset resulted in a single NJ tree representing 151 thrips species, which formed distinct haplotype clusters (Fig. 1). The intraspecific and intrageneric sequence divergence ranged from 0.0 to 7.91% and 8.65% to 31.15%, respectively (Table 2 and Table 3, respectively). This discrete barcoding gap between intra and inter specific distances (Hebert et al. 2004) allowed us to clearly distinguish all the thrips species employed in this study. Besides being direct pests, many of the thrips species are known to vector Tospoviruses (Mound et al. 1996). In this regard, it is important to analyze the diversity within thrips species such as, Thrips palmi, T. tabaci, Frankliniella occidentalis, Scirtothrips dorsalis, F. schultzei, Thrips hawaeiinsis, S. perseae and Chirothrips manicatus, etc. For the first time, we were able to record the existence of possible cryptic species within 2 Thrips spp. viz. T. hawaiiensis and S. perseae along with previously reported cryptic species viz., T. palmi (Rebijith et al. 2011; Glover et al. 2010), T. tabaci (Brunner et al. 2004), F. occidentalis (Rugman- Jones et al. 2010), S. dorsalis (Rebijith et al. 2011; Kadirvel et al. 2013) based on the ‘CO-I- 10X barcoding Rule' proposed by Hebert et al. 2004. All of them were supported by the calculated intra and inter specific genetic divergence (Table. 4) for different sub clusters as shown in Fig. 2 (AF).

Thrips palmi

Thrips palmi was represented by 145 specimens forming 2 lineages. One lineage (n = 102) was consisted of the Indian T. palmi population collected from various host plants, and the second lineage (n = 43) represented the world population with T. palmi representatives from China, Dominican Republic, Thailand, Japan and United Kingdom. There was much less variation (0.4% to 0.7%) within these lineages, however 8.2% difference (hereafter, D) was observed between them (Fig. 2A, Table 4), which is indicative of geographically isolated cryptic species.

Thrips tabaci

A total of 146 of T. tabaci specimens were analyzed and found to form 2 clades. The first clade comprised a T. tabaci population representing various geographical locations across the globe from different host plants. However, the second clade represented a population collected only from the tobacco plant, Nicotiana tabacum. As an indicator of cryptic species, there was less variation evident within these clades (0.3% to 0.9%), whereas 9.0% D observed between them (Fig. 2B, Table 4).

Scirtothrips dorsalis

Scirtothrips dorsalis was represented by 39 specimens forming 4 lineages, among which no perfect lineages - neither with geographical locations nor host associated genetic differences - were formed. However, very little variation (0 to 0.96%) was observed within these lineages, and 1.84% to 10.07% D observed between them, which was indicative of the cryptic speciation in S. dorsalis (Fig. 2C, Table 4).

Frankliniella occidentalis

Frankliniella occidentalis was represented by 224 specimens, which formed 5 lineages of which none was > 500 bp in length. However, the genetic variation within and between the major 2 clades were 0 to 0.10% and 1.11 to 6.2%, respectively (Fig. 2D, Table 4).

Thrips hawaiiensis

Thrips hawaiiensis was represented by 22 specimens forming 5 lineages (Fig. 2E) with genetic variation 0 to 0.5% within these lineages and 3.4% to 9.5% between them (Table 4).

Scirtothrips perseae

Scirtothrips perseae was represented by 47 specimens forming 3 lineages, of which the first, second and third were represented by 6, 4 and 37 specimens respectively, with intra and inter subcluster values ranging from 0.07% to 0.91% and 5.1 to 11.6%, respectively (Fig. 2F, Table 4).


Rapid or timely and accurate identification of invasive pests, such as thrips, is important and chal lenging worldwide, as these pests can cause crop damage either by direct feeding or by transmitting plant pathogens (German et al. 1992; Hebert et al. 2004). In this regard, classical taxonomy has its own strength, however DNA barcoding employing CO-I has the added advantage of not being limited by polymorphism, sex, and life stage of the target species (Rebijith et al. 2013). Our study is the first attempt to examine a large number of thrips species on a global scale, which include large number of genera represented by multiple species, where we expected that lower intra specific distances may cause problems in delimiting species. However, barcoding employing the CO-I gene allowed us to accurately discriminate all 151 species of thrips. Thus, CO-I DNA barcodes proved to be an invaluable tool for delimiting thrips species, an approach complementing classical taxonomy in the context of effective plant quarantine and biological control initiatives (Rugman-Jones et al. 2006).

Fig. 1.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing clusters of thrips species for CO-I sequences. Distinct clades for 151 species of thrips can be seen in the figure, in which both T. palmi and T. tabaci consist of two groups and 4 other species viz.T. hawaiiensis, S. dorsalis, S. perseae and F. occidentalis are each marked with a dark circle as an indicator of possible cryptic species. Dendrograms for these 5 species are shown in Figs. 2A, 2B, 2C, 2D and 2F, respectively. Bootstrap value > 65% only showed in the figure.


Table 3.

Intrageneric K2P genetic divergence within different genera of thysanoptera.


Genetic Divergence

DNA barcoding has become an important tool in species identification, and has improved our ability to understand diversity among populations through to higher level taxa (Puckridge et al. 2013). Researchers employed a 2% genetic divergence cutoff value for species identification. However, many exceptional cases are reported with lower interspecific genetic divergences, yet most of the cases were still able to be correctly dissected out amongst most of the species (Pfannenstiel et al. 2008). By employing this 2% cutoff criterion for species delimitation, we aptly identified all thrips species in this study.

The mean K2P genetic distance values found for conspecific and congeneric comparison were 3.5% and 17.3% respectively, which were found to be on par with previous studies (Glover et al. 2010). Additionally, we analyzed a large number of closely related species, e.g., species of Genus Thrips, our observed mean congeneric divergence value was slightly smaller than in previous studies (Glover et al. 2010). This could be either due to the possible recent species radiation of Thrips similar to that of fresh water fish fauna (Albert et al. 2011; Bermingham et al. 1998; Montoya- Burgos et al. 2003), or the possible evolutionary rate variation of CO-I among different species employed in this study.

Advantages of DNA Barcoding in Taxonomy

Classical taxonomy employing morphological characters cannot be used for all species of thrips in all life stages because of insufficient phenotypic variation (Brunner et al. 2004). In addition, the presence of unusual morphological forms of species on different hosts, minute size, co-existence of different species on the same host, complex life cycles, color morphs, and parthenogenetic mode of reproduction, etc., add to the difficulties of precise identification. Furthermore, morphological examination of thrips to species requires adult specimens as there are no reliable keys for identification of immature stages, but even these are often difficult for a non-expert to use. In this regard, DNA barcoding can be an added advantage and an effective tool for molecular species identification (Brunner et al. 2004; Glover et al. 2010; Lee et al. 2010), elucidation of biotypes or cryptic species, and host associated genetic differences (Shufran et al. 2000; Brunner et al. 2004) and species discovery (Foottit et al. 2010) in insects. Additionally, DNA barcoding could also be used in the identification of unknown thrips species that co-exist in a cropping system, since one particular Tospovirus disease may sometimes be vectored by more than one species of thrips (Wijkamp et al. 1995; Kadirvel et al. 2013; Amin et al. 1981).

Table 4.

List of species with higher intraspecific k2p divergences.


Fig. 2A.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing clustering of T. palmi species for CO-I sequences. Two distinct groups can be seen Group-I (Thrips palmi- RWD, represents the T. palmi populations from various countries viz. Japan, China, Dominican Republic, Thailand and U.K.) and Group-II (Thrips palmi- IND, which clearly associated with Indian subcontinent), with 100% bootstrap support.


In the recent past, DNA barcoding has become an effective tool in revealing cryptic and potentially new species, which has increased our knowledge of biodiversity (Rebijith et al. 2013; Puckridge et al. 2013). In this study, 7 species showed conspecific genetic divergence values ≥ 2%, and were subdivided into further sub clusters (Fig. 3 a– h). Hebert et al. (2004) proposed the ‘CO-I DNA 10X barcoding rule', whereby identification of cryptic species is possible if the 2 lineages diverge by 10 times or more the average intra sub cluster variability within these lineages. In this study, such a sub cluster analysis has revealed tight clusters with significant larger mean values (2.57% to 18.02%), and smaller mean values within subclusters (0.00 to 3.28%) (Table 5). Such deep conspecific divergence has been reported previously in DNA barcoding and most of the entities are considered as cryptic species (Handfield & Handfield et al. 2006; Smith et al. 2006; Gomez et al. 2007; Pfenninger et al. 2007; April et al. 2011). The possible reasons for such higher conspecific genetic divergences are (i) representation of a novel species (ii) misidentification of species, and (iii) host/geographic preference of the organism (Brunner et al. 2004). However, according to Shao et al. (2003), the 10X estimation may be too low for Thrips, since they are known to have greater mitochondrial gene rearrangements and molecular evolution.

On the other hand, 2 Thrips species, viz. T. palmi and T. tabaci, require special attention because of their geographic and host associated genetic differences. Thips palmi formed 2 lineages; of which one lineage is unique to India (IND) and second lineage represents the rest of the world (ROW) (Rebijith et al. 2011). Geographic isolation and genetic drift can contribute strongly to intraspecific phylogeographic structure (Avise et al. 1987). Such allopatric lineages reinforce the fact that both these populations (IND and ROW) have independent evolutionary histories, which could be explained by restricted gene flow due to many physical and chemical barriers (April et al. 2011). Whereas in the case of T. tabaci, the 2 biotypes viz. ‘tabaci' and ‘communis', can be clearly distinguished with a character on the abdomen of the second larval stage (Zawirska et al. 1976). Brunner et al. (2004) proved the existence of 3 host associated genetic groups, viz. L1, L2 (Leak plant- ‘communis type' sensu Zawirska) and T (Tobacco plant-‘tabaci type') and proposed that when host fidelity is perfect, reproductive isolation is complete. In recent studies Jacobson et al. (2013) grouped T. tabaci into two entities based on their reproductive mode. However, our results showed 2 lineages, one with tobacco plants and the second is with various host plants such as onion, leak, etc., on par with previous studies (Zawirska et al. 1976; Brunner et al. 2004). In a nutshell, we conclude that all these cases represent cryptic species, but both T. palmi and T. tabaci can be subdivided and are strong candidates for novel species.

Fig. 2B.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing clustering of T. tabaci species for CO-I sequences. Two distinct groups can be seen Group-I (Thrips tabaci- ‘tabaci type'- represents the T. tabaci populations from host plant tobacco, Nicotiana tabacum and Group-II (Thrips tabaci- ‘communis type'- which clearly associated with various host plants such as onion, leek, etc.) with 100% bootstrap support.


DNA Barcoding in Biosecurity

One of the major threats to crop production and productivity is the spread of invasive species such as thrips, aphids, whiteflies and planthoppers, etc., within and outside the country. Many of the monitoring programs are challenged by both quantity and quality of materials at the port of entry as well as the taxonomic breadth of intercepted organisms. Furthermore, morphological examination of majority of agricultural invasive pests such as thrips to species level is usually restricted to (i) adult specimens as there are no adequate keys for eggs, larvae or pupae (ii) polymorphism and (iii) sexual forms. At this juncture, CO-I based molecular identification has an added advantage of not being limited by any of the above factors, and DNA barcoding has shown success in discriminating arthropods of quarantine importance such as Tephritidae, Lymantriidae (Armstrong et al. 2005), Spodoptera spp. (Nagoshi et al. 2011), Liriomyza spp. (Scheffer et al. 2006) and Thysanoptera (Qiao et al. 2012). Thus, DNA barcoding can play a vital role in the international bio-security as an additional diagnostic protocol for quarantine pests at the port of entry (FAO 2006). DNA barcoding technique has been employed in New Zealand since 2005 for the identification of Tephritidae and Lymantriidae (Armstrong et al. 2005). Thrips species used in the current study were differentiated on the basis of DNA barcodes, which proved to be an invaluable tool for identification of Thysanoptera invasive insect pests, an approach complementing classical taxonomy. Having all these advantages, DNA barcoding was included as an effective diagnostic tool for invasive insect pest in the International Plant Protection Convention (IPPC) (FAO 2006).

Fig. 2C.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing clustering of S. dorsalis species for CO-I sequences. Four distinct groups can be seen, among which group-III corresponds to S. dorsalis collected on various host plants from China. Sequences generated from this study formed group I along with other sequences from Japan and Thailand.


Implications of DNA Barcoding for Insect Pest Management

Thrips continue to pose dual problems by either direct feeding or by transmitting plant pathogenic tospoviruses in both field and green house conditions (Mound et al. 1996; Rebijith et al. 2011). Control of thrips with insecticides is a difficult task due to their parthenogenetic mode of reproduction (Rebijith et al. 2013), mode of life cycles (Brunner et al. 2004), etc. Eggs and pupae are protected either in leaf tissue or in soil/leaf litter and the larvae & adults are protected within buds and flowers (Glover et al. 2010). Yet, farmers employ insecticides as a primary control measure for the management of thrips, which could ultimately lead to the insecticide resistance.

Fig. 2D.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing five distinct clusters for F. occidentalis CO-I sequences. Number in bracket indicates, the individuals formed that cluster.


Fig. 2E.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing five clusters for T. hawaiiensis CO-I sequences.


Plant disease management requires accurate identification of species to understand the biology, population structure and ecology of the species (Rebijith et al. 2013). In this connection, T. palmi and T. tabaci, which infest watermelon to which they transmit Watermelon Bud Necrosis Virus, WBNV, and onion to which they transmitting Iris Yellow Spot Virus, IYSV, respectively, demand quick control measures using pesticides in order to limit the spread of these most potent viruses. On the other hand, Thrips flavus and Thrips nigropilosus (morphologically very similar to T. palmi and T. tabaci), can be managed effectively by employing biological control agents known as entomopathogenic nematodes (EPNs), such as Thripenema fuscum n. sp. (Asokan and Rebijith, unpublished data).

Reducing unnecessary pesticide usage through proper species identification can save growers money, and reduce chances of development insecticide resistance, as well as be more environmentally protective.


In this study, DNA barcoding proved to be an effective tool that can be employed for species identification, elucidation of cryptic species, biotypes and also in the discovery of new species. We trust that our work will serve as a rapid, precise, independent identification approach for the discrimination of thrips species of different life stages and color morphs, both for the species presently studied, and in the future, for other pest species of agricultural, horticultural and forestry interest and importance. This will in turn help in elucidation of the epidemiology of tospoviruses, their management and serve as a potentially valuable tool in quarantine at ports-of-entry. Moreover, as our study has revealed, the existence of cryptic thrips species in Thrips hawaiiensis and Scirtothrips perseae and shows that further studies on the evolution of these particular species (and doubtless others too) are required before we can certain that we are dealing with sensu stricto taxa rather than sensu lato taxa.

Fig. 2F.

Neighbor-Joining tree with bootstrap support (2,000 replicates) showing three clusters for S. perseae CO-I sequences.



Our sincere thanks are due to Dr. Kaomud Tyagi, NBAII, Bangalore for morphological identification of all the thrips species used in the current study. Financial support for field and laboratory work was provided by the Indian Council for Agricultural Research (ICAR), New Delhi through ‘Out Reach Program on Management of Sucking Pests on Horticultural Crops' under 11th plan. The authors are grateful to Chaitanya, B. N and Roopa, H. K for the technical assistance and thoughtful comments. This work is a part of the Ph. D thesis of the senior author.

References Cited


J. S. Albert , and R. E. Reis 2011. Historical Biogeography of Neotropical Freshwater Fishes. Berkeley, Los Angeles, London: Univ. California Press. Google Scholar


P. W. Amin , D. V. R. Reddy , and A. M. Ghanekar 1981. Transmission of tomato spotted wilt virus, causal agent of bud necrosis of peanut, by Scirtothrips dorsalis and Frankliniella schultzei. Plant Dis. 65: 663–665. Google Scholar


J. April , R. L. Mayden , R. Hanner , and L. Bernatchez 2011. Genetic calibration of species diversity among North America's freshwater fishes. Proc. Natl. Acad. Sci. USA. 108(26): 10602–10607. Google Scholar


K. F. Armstrong , and S. L. Ball 2005. DNA barcodes for biosecurity: Invasive species identification. Phil. Trans. R. Soc. London Series, B. 360: 1813–1823. Google Scholar


R. Asoka n , K. B. Rebijith , S. K. Singh , A. S. Sidhu , S. Siddharthan , P. K. Karanth , R. Ellango , and V. V. Ramamurthy 2011. Molecular identification and phylogeny of Bactrocera species (Diptera: Tephritidae). Florida Entomol. 94(4): 1026–1035. Google Scholar


J. C. Avise , J. Arnold , M. R. Ball , E. Bermingham , T. Lamb , J. E. Neigel et al. 1987. Intraspecific phylogeography: The mitochondrial bridge between population genetics and systematics. Annu. Rev. Ecol. Syst. 18: 489–522. Google Scholar


E. Bermingham , and E. P. Martin 1998. Comparative mtDNA phylogeography of Neotropical freshwater fishes: Testing shared history to infer the evolutionary landscape of lower Central America. Mol. Ecol. 7: 499–517. Google Scholar


P. C. Brunner , E. K. Chatzivassilious , N. I. Katis , and J. E. Frey 2004. Host-associated genetic differentiation in Thrips tabaci (Insecta: Thysanoptera), as determined from mtDNA sequence data. Heredity 93: 364–370. Google Scholar


M. Ciuffo , G. C. Mautino , L. Bosco , M. Turina , and L. Tavella 2010. Identification of Dictyothrips betae as the vector of Polygonum ring spot virus. Ann. Appl. Biol. 157: 299–307. Google Scholar


FAO. 2006. Global forest resources assessment 2005: Progress towards sustainable forest management. Forestry Paper 147. Available from: Accessed X-2007). U. N. Food and Agric. Org. Rome. Google Scholar


R. G. Foottit , H. E. L. Maw , K. S. Pike , and R. H. Miller 2010. The identity of Pentalonia nigronervosa Coquerel and P. caladii van der Goot (Hemiptera: Aphididae) based on molecular and morphometric analysis. Zootaxa 2358: 25–38. Google Scholar


T. L. German , D. E. Ullman , and J. W. Moyer 2006. Tospoviruses: Diagnosis, molecular biology, phylogeny and vector relationships. Annu. Rev. Phytopathol. 30: 315–348. Google Scholar


R. H. Glover , D. W. Collins , K. Walsh , and N. Boonham 2010. Assessment of loci for DNA barcoding in the genus Thrips (Thysanoptera: Thripidae). Mol. Ecol. Res. 10: 51–59. Google Scholar


A. Gomez , P. J. Wright , D. H. Lunt , J. M. Cancino , G. R. Carvalho , and R. N. Hughes 2007. Mating trails validate the use of DNA barcoding to reveal cryptic speciation of a marine bryozoans taxon. Proc. Biol. Sci. 274: 199–207. Google Scholar


M. Hajibabaei , D. H. Janzen , J. M. Burns , W. Hallwachs , and P. D. N. Hebert 2006. DNA barcodes distinguish species of tropical Lepidoptera. Proc. Natl. Acad. Sci. USA 103(4): 968–971. Google Scholar


T. A. Hall 1999. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids. Symp. Ser. 41: 95–98. Google Scholar


D. Handfield , and L. Handfield 2006. A new species of Plusia (Lepidoptera : Noctuidae) from North America. Canadian Entomol. 138: 853–859. Google Scholar


A. Hassani-Mehraban , M. Botermans , J. T. Verhoeven , E. Meekes , J. Saaijer , D. Peters , R. Goldbach , and R. Kormelink 2010. A distinct tospovirus causing necrotic streak on Alstroemeria sp. in Colombia. Arch. Virol. 155: 423–428. Google Scholar


P. D. N. Hebert , A. Cywinska , S. L. Ball , and J. R. Dewaard 2003a. Biological identifications through DNA barcodes. Proc. R. Soc. Series B. 270: 313–321. Google Scholar


P. D. N. Hebert , S. Ratnasignham , and J. R. Dewaar d 2003b. Barcoding animal life: Cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Series B 270 (Suppl. 1): S96–S99. Google Scholar


P. D. N. Hebert , E. H. Penton , J. M. Burns , D. H. Janzen , and W. Hallwachs 2004. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc. Natl. Acad. Sci. USA 101: 14812–14817. Google Scholar


D. R. Jones 2005. Plant viruses transmitted by thrips. European J. Plant Pathol. 113: 119–157. Google Scholar


P. Kadirvel , R. Srinivasan , Y. C. Hsu , F. C. Su , and R. De La Pena 2013. Application of cytochrome oxidase I sequences for phylogenetic analysis and identification of thrips species occurring on vegetable crops. J. Econ. Entomol. 106(1): 408–418. Google Scholar


M. Kimura 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16: 11–120. Google Scholar


S. Kumar , K. Tamura , and M. Nei 1993. MEGA (Molecular Evolutionary Genetics Analysis). Pennsylvania State Univ., University Park, PA. Google Scholar


W. Lee , H. Kim , J. Lim , H. Choi , Y. Kim , J. Ji , R. G. Foottit , and S. Lee 2010. Barcoding aphids (Hemiptera: Aphididae) of the Korean Peninsula: Updating the global data set. Mol. Ecol. Res. 11: 1–6. Google Scholar


C. P. Meyer , and G. Paulay 2005. DNA Barcoding: error rates based on comprehensive sampling. PLoS Biology 3(12), e422. Google Scholar


J. I. Montoya-Burgos 2003. Historical biogeography of the catfish genus Hypostomus (Siluriformes: Loricariidae), with implications on the diversification of Neotropical ichthyofauna. Mol. Ecol. 12: 1855–1867. Google Scholar


G. Moritz. , L. A. Mound , D. C. Morris , and A. Goldarazena 2004. Pest thrips of the world: Visual and molecular identification of pest thrips (CDROM), CBIT, Univ. Brisbane, ISBN-1-86499-781-8.  Google Scholar


L. A. Mound 1996. The Thysanoptera vector species of Tospoviruses. Acta Hort. 431: 298–309. Google Scholar


L. A. Mound 2005. Thysanoptera: Diversity and interactions. Annu. Rev. Entomol. 50: 247–269. Google Scholar


L. A. Mound , and D. C. Morris 2007a. A new thrips pest of Myoporum cultivars in California, in a new genus of leaf-galling Australian Phlaeothripidae (Thysanoptera). Zootaxa 1495: 35–45. Google Scholar


L. A. Mound , and D. C. Morris 2007b. The insect order Thysanoptera: Classification versus systematics. Zootaxa 1668: 395–411. Google Scholar


L. A. Mound , G. S. Wheeler , and D. A. Williams 2010. Resolving cryptic species with morphology and DNA: Thrips as a potential biocontrol agent of Brazilian peppertree, with a new species and overview of Pseudophilothrips (Thysanoptera). Zootaxa 2432: 59–68. Google Scholar


L. A. Mound , and R. Zur Strassen 2001. The genus Scirtothrips (Thysanoptera: Thripidae) in Mexico: A critique of the review by Johansen and Mojica- Guzman. Florida Entomol. 40: 133–141. Google Scholar


T. Murai , and S. Toda 2001. Variations of Thrips tabaci in color and size, pp. 377–378 In Proc. 7th Intl. Symp. Thysanoptera. Reggio Calabria, Italy. Google Scholar


R. N. Nagoshi , J. Brambila , and R. L. Meagher 2011. Use of DNA barcodes to identify invasive armyworm Spodoptera species in Florida. J. Insect Sci. 11: 154 available online: Google Scholar


H. R. Pappu , R. A. C. Jones , and R. K. Jain 2009. Global status of tospovirus epidemics in diverse cropping systems: successes achieved and challenges ahead. Virus Res. 141: 219–236. Google Scholar


L. H. G. Pereira , R. Hanner , F. Foresti , and C. Oliveira 2013. Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna? BMC Genetics 14: 20. Google Scholar


R. S. Pfannenstiel 2008. Spider predators of Lepidopteran eggs in South Texas field crops. Biol. Control 46: 202–208. Google Scholar


M. Pfenninger , C. Nowak , C. Kley , D. Steinke , and B. Streit 2007. Utility of DNA taxonomy and barcoding for the inference of larval community structure in morphologically cryptic Chironomus (Diptera) species. Mol. Ecol. 16: 1957–1968. Google Scholar


M. Puckridge , N. Andreakis , S. A. Appleyard , and R. D. Ward 2013. Cryptic diversity in flathead fishes (Scorpaeniformes: Platycephalidae) across the Indo-West Pacific uncovered by DNA barcoding. Mol. Ecol. Res. 13: 32–42. Google Scholar


W. N. Qiao , F. H. Wan , A. B. Zhang , L. Min , and G. F. Zhang 2012. Application of DNA barcoding technology for species identification of common thrips (Insecta: Thysanoptera) in China. Acta Entomol. Sinica 55: 344–356. Google Scholar


K. B. Rebijith , R. Asokan , N. K. Krishna Kumar , V. Krishna and V. V. Ramamurthy 2011. Development of species-specific markers and Molecular differences in mtDNA of Thrips palmi Karny and Scirtothrips dorsalis Hood (Thripidae: Thysanoptera), vectors of tospoviruses (Bunyaviridae) in India. Entomol. News 122(3): 201–213. Google Scholar


K. B. Rebijith , R. Asokan , N. K. Krishna Kumar , V. Krishna , B. N. Chaitanya , and V. V. Ramamurthy 2013. DNA barcoding and elucidation of cryptic aphid species (Hemiptera: Aphididae) in India. Bull. Entomol. Res. 103: 601–610. Google Scholar


K. B. Rebijith , R. Asokan , N. K. Krishna Kumar , K. K. Srikumar , V. V. Ramamurthy , and P. Shivarama Bhat 2012. DNA barcoding and development of species-specific markers for the identification of tea mosquito bugs (Miridae: Heteroptera) in India. Environ. Entomol. 41(5): 1239–1245. Google Scholar


R. J. Rowley , J.C. Kostelnick , D. Braaten , X. Li , and J. Meisel 2007. Risk of rising sea level to population and land area. Eos, Trans. American Geophys. Union: 88(9): 105–107. Google Scholar


P. F. Rugman-Jones , M. S. Hoddle , and R. Stouthamer 2010. Nuclear-mitochondrial barcoding exposes the global pest western flower thrips (Thysanoptera: Thripidae) as two sympatric cryptic species in its native California. J. Econ. Entomol. 103: 877–886. Google Scholar


P. F. Rugman-Jones , M. S. Hoddle , L. A. Mound , and R. Stouthamer 2006. Molecular identification key for pest species of Scirtothrips (Thysanoptera: Thripidae). J. Econ. Entomol. 99: 1813–1819. Google Scholar


N. Saitou , and N. Nei 1987. The Neighbor-Joining method: A new method for reconstructing phylogenetic trees. Mol. Bio. Evol. 4: 406–425. Google Scholar


S. J. Scheffer , M. L. Lewis , and R. C. Joshi 2006. DNA barcoding applied to invasive leafminers (Diptera: Agromyzidae) in the Philippines. Ann. Entomol. Soc. America 99: 204–210. Google Scholar


R. Shao , M. Dowton , A. Murrel , and S. C. Barker 2003. Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects. Mol. Biol. Evol. 20: 1612– 1619. Google Scholar


K. A. Shufran , J. D. Burd , J. A. Anstead , and G. Lushai 2000. Mitochondrial DNA sequence divergence among greenbug (Homoptera: Aphididae) biotypes: evidence for host-adapted races. Insect Mol. Biol. 9: 179–184. Google Scholar


M. A. Smith , N. E. Woodley , D. H. Janzen , W. Hallwachs , and P. D. N. Hebert 2006. DNA barcodes reveal cryptic host-specificity within the presumed polyphagous members of a genus of parasitoid flies (Diptera: Tachinidae). Proc. Natl. Acad. Sci. USA 103: 3657–3662. Google Scholar


A. E. Timm , M. Stiller , and J. E. Frey 2008. A molecular identification key for economically important thrips species (Thysanoptera: Thripidae) in southern Africa. African Entomol. 16: 68–75. Google Scholar


S. Toda , and S. Komazaki 2002. Identification of thrips species (Thysanoptera: Thripidae) on Japa nese fruit trees by polymerase chain reaction and restriction fragment length polymorphism of the ribosomal ITS2 region. Bull. Entomol. Res. 92: 359–363. Google Scholar


D. E. Ullman , J. L. Sherwood , and T. L. German 1997. Thrips as vectors of plant pathogens, pp. 539–565 In T. Lewis [ed.], Thrips as Crop Pests. CAB Intl. New York, NY. Google Scholar


J. F. Wang , L. Y. Jiang , and G. X. Qiao 2011. Use of a mitochondrial COI sequence to identify species of the subtribe Aphidina (Hemiptera, Aphididae). ZooKeys 122: 1–17. Google Scholar


I. R. Wijkamp , R. Almarza , Goldbach , and D. Peters 1995.Distinct levels of specificity in thrips transmission of tospoviruses. Phytopathology 85: 1069– 1074. Google Scholar


I. Zawirska 1976. Untersuchungen uber zwei biologische Typen von Thrips tabaci Lind. (Thysanoptera, Thripidae) in der VR Polen. Arch. Phytopathol. Pflanzenschutz 12: 411–422. Google Scholar


G. F. Zhang , X. Q. Meng , L. Min , W. N. Qiao , and F. H. Wan 2011. Rapid diagnosis of the invasive species, Frankliniella occidentalis (Pergande): a species-specific COI marker. J. Appl. Entomol. (doi:10.1111/j.1439-0418. 2011.01661.x). Google Scholar


[1] Supplementary material for this article in Florida Entomologist 97(4) (2014) is online at

K. B. Rebijith, R. Asokan, V. Krishna, H. H. Ranjitha, N. K. Krishna Kumar, and V. V. Ramamurthy "DNA Barcoding and Elucidation of Cryptic Diversity in Thrips (Thysanoptera)," Florida Entomologist 97(4), 1328-1347, (1 December 2014).
Published: 1 December 2014
código de barras del ADN
distancias intra e interespecíficas
DNA barcoding
intra and inter-specific distances
Back to Top