Editorial Type:
Article Category: Research Article
 | 
Online Publication Date: 05 Jul 2021

Sequence Analysis of Mitochondrial Cytochrome Oxidase 1 from Bemisia tabaci (Hemiptera: Aleyrodidae) Populations in Iran

and
Page Range: 366 – 373
DOI: 10.18474/JES21-01
Save
Download PDF

Abstract

Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) has been phylogenetically divided into several putative species based on nucleotide sequence of mitochondrial cytochrome oxidase 1 (mtCOI). To confirm the putative species among B. tabaci populations in Iran, insects were collected from different regions and plant hosts in Khuzestan province of southwestern Iran, and their mtCOI sequences were determined. DNA polymorphism and phylogenetic analysis of the 10 mtCOI sequences obtained showed four haplotypes among the specimens from Khuzestan province that were identified as the putative species MEAM1-subcladeB. Moreover, specimens previously collected from Iran were found to be members of the MEAM1-subcladeB2. Twelve variable sites were detected throughout the Iran-originated mtCOI sequences contributing to their position on the phylogenetic tree. This is the first study reporting the putative species MEAM1-subcladeB from southwestern Iran.

Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) is a devastating phloem-feeding pest of many crops (Basu et al. 2019) that transmits several plant viruses (Gilbertson et al. 2015, Pinheiro-Lima et al. 2020). Based on the mitochondrial cytochrome oxidase 1 (mtCO1) sequence, 39–42 cryptic (putative) species have been described within the species (Kanakala and Ghanim 2019). In Iran, it has been showed that the putative species of Africa/Middle East/Asia Minor group and the Middle East/Asia Minor 1 subgroup is the dominant species (Shahbazi et al. 2014). However, there is no information on the putative species within Khuzestan province of southwestern Iran. The objective of this study was to determine the putative species among B. tabaci populations from the region. Moreover, the mtCOI sequences of Iranian specimens were subjected to DNA polymorphism and phylogenetic analysis.

Materials and Methods

Bemisia tabaci adults were collected from tobacco, tomato, cucumber, and cowpea fields in six counties of Khuzestan province during 2019–2020 (Table 1). Specimens were preserved individually in microtubes containing 70% alcohol and stored at -20°C. At each sampling, 20–30 whiteflies representing 10 distinct populations were collected. Morphologic characters of pupae from each population were used to confirm species identity (Hodges and Evans 2005).

Table 1 Bemisia tabaci specimens collected from different regions in Khuzestan province, southwestern Iran.
Table 1

Total DNA was extracted from the specimens using the Animal DNA Isolation Kit (Denazist, Iran) and stored at -20°C. Polymerase chain reaction (PCR) amplification was performed according to Masood et al. (2017). PCR products were then sequenced at BIONEER Corp. (Daejeon, South Korea).

Sequences were primarily checked using Chromas (ver. 2.6.6) (Technelysium, South Brisbane, Australia) software, and a complete sequence was obtained by assembling the reads of 3′ and 5′ termini using SeqMan Pro software (DNASTAR Lasergene, ver. 8). Nuclear mitochondrial sequences (NUMTs) and PCR-derived pseudogenes were excluded according to Kunz et al. (2019). BLASTn was performed, and the confirmed sequences were then submitted to GenBank under the accession numbers presented in Table 2. DNA polymorphism and haplotype analysis were performed using DnaSP software (ver. 6). Haplotype diversity and sampling variance were conducted according to Nei (1987). Nucleotide diversity (π), mean nucleotide differences per site, sampling variance, and standard deviation of the sequences were determined. Theta (h) per site from π, S, and g and mean the number of nucleotide differences (k) were also calculated. Diallelic model was applied to analyze insertion/deletion (InDel) polymorphism among the sequences, and possible InDel haplotypes were investigated.

Table 2 Mitochondrial cytochrome oxidase 1 sequences of B. tabaci specimens from different provinces in Iran.
Table 2

Fifty-six sequences of mtCOI, including 4 sequences obtained here and 52 sequences from GenBank, were used to perform a phylogenetic analysis. The sequences were aligned by CLC Main Workbench software (ver. 7.6.2), and the maximum likelihood algorithm using the general time reversible model test with 100 bootstrap replicates was used to construct nucleotide-based phylogenetic trees. A corresponding sequence from Trialeurodes abutilonea Haldeman was used as an outgroup.

Results

The examined pupal specimens from the populations were identified as B. tabaci. The partial mtCOI gene (780 bp) was amplified using the specific primers, and four haplotypes (H) were identified among the sequences (Table 3). The sequences from Khuzestan province exhibited the least values of polymorphism parameters including the number of polymorphic sites (S), haplotype diversity (HD), Pi, Theta (per site) from Eta (Thetaeta), Theta (per site) from S (Theta-W), and K. In contrast, the highest values of S, HD, Pi, Thetaeta, Theta-W, and K were observed in the population from Kerman, Hormozgan, Khorasan, and Yazd provinces. The sequences from Fars province were placed in the second rank relative to the values. No InDel polymorphic event was observed among the sequences.

Table 3 DNA polymorphism parameters of B. tabaci populations examined in this study using DNASP software (ver. 6).
Table 3

All sequences obtained herein were highly identical (99.88%) to a sample from Pakistan (LN835385). The second highest identity (99.84%) was found with an isolate of MEAM1-subcladeB from the UAE (DQ133382). Additionally, a 99.69% identity was observed against three Iraqi isolates (HM070413, KX679575, and KX679577) and one isolate from the Kingdom of Saudi Arabia (GU086344). Two Syrian isolates (KP342512 and AB473559) exhibited high levels of identity (99.74 and 99.61%, respectively) to the sequences determined here. Isolates from Yemen and Kuwait (GU086343 and GU086346, respectively) showed a 98.96% identity to Khuzestan-originated sequences. A putative species MEAM1-subcladeB2 from Pakistan (GU977267) differed by 21 nucleotides from the sequences obtained herein. The lowest identity (95.06%) was observed against an isolate of the putative species MEAM2 from Uganda (KX570778). The sequences obtained here were identified as putative species MEAM1-subcladeB.

The phylogenetic tree showed the position of the putative species MEAM1-subcladeB (MW161484, MW161487, MW161488, and MW161490) from the present study (Fig. 1). These sequences, together with the isolate from the UAE (DQ133382), formed a separate cluster on the tree. The Iranian haplotypes reported by Shahbazi et al. (2014) formed two separate clusters on the tree. One cluster containing four isolates (JN559744, JN559748, JN559749, and JN559750) were clustered together with an Iraqi MEAM1-subcladeB2 species (KX679577). The other cluster consisted of two isolates (JN559751 and JN559753), which were placed close to a MEAM1-subcladeB2 species from Syria (AB473559 and KP342512). Another Iranian haplotype reported by Shoorcheh et al. (2008) (EU547771) was clustered with a MEAM1-subcladeB2 species from Pakistan (GU977267). Similarly, one unpublished Iranian specimen (MT038431) was placed close to the MEAM1-subcladeB2 species from Pakistan (GU977267) (Fig. 1).

Fig. 1Fig. 1Fig. 1
Fig. 1 Circular cladogram created using multiple alignments of 56 mitochondrial cytochrome oxidase 1 sequences from Bemisia tabaci specimens all around the world. The phylogenetic tree was constructed using maximum likelihood and bootstrap (100 replicates) by CLC Main Work-bench 7.6.2. The four Iranian B. tabaci haplotypes B obtained in this study are marked with an asterisk. The MEAM1 and MEAM2 sequences were obtained from GenBank database. The corresponding sequence of Trialeurodes abutilonea were used as outgroups (Og). All other sequences are available at CSIRO (De Barro and Boykin 2013).

Citation: Journal of Entomological Science 56, 3; 10.18474/JES21-01

Multiple alignments of the sequences showed 12 nucleotide positions in which genetic variation had occurred. The sequences from Fars province had seven nucleotide sites (positions: 139, 193, 242, 468, 493, 624, and 643) with genetic variation. The sequences from Khuzestan province demonstrated three variable nucleotide sites at positions 58, 548, and 606. One sequence from Khorasan province (EU547770) showed the highest number (four) of nucleotide positions (58, 369, 548, 606, and 692) with genetic variation, whereas another sequence from the same province (MT038431) exhibited no variation. Additionally, a specimen from Yazd province (EU547771) showed only one nucleotide position (692) in which genetic variation was observed. The remaining specimens showed no variation throughout the sequence.

Discussion

Here, we report for the first time the mtCOI sequences of B. tabaci populations from Khuzestan province, southwestern Iran. We demonstrated that these sequences belong to the putative species MEAM1-subcladeB. Additionally, the sequence analysis of Iranian specimens has led to the determination of their subclade (B and B2), which had not been previously reported.

The Africa/Middle East–Asia Minor genetic group was found to be the main putative species in Iran that also has been reported from neighboring countries including Iraq (Kareem et al. 2020) and Pakistan (Masood et al. 2017). The MEAM1 group that contains two subclades (B and B2) has been observed in tropical and subtropical regions globally. The species have been considered invasive and cosmopolitan (De Barro et al. 2011). It appears that MEAM1-subcladeB and MEAM1-subcladeB2 are the dominant putative species in Iran. The MEAM1-subcladeB presented here was identical to the previous sequence (EU547770) reported from Khorasan province (Shoorcheh et al. 2008). It had shown a genetic variation among the sequences determined by Shoorcheh et al. (2008) resulting in formation of a separate branch on the phylogenetic tree (Fig. 1). Accordingly, this sequence was the only Iranian isolate with a close phylogenetic relationship to the sequences obtained here (Fig. 1). Although Fars and Khuzestan specimens showed no diversity in their putative species (MEAM1-subcladeB2 and -subcladeB, respectively), no location-dependent distribution of the putative species was observed. Accordingly, a single putative species (MEAM1-subcladeB2) has been reported from different provinces (Yazd and Khorasan) (Shoorcheh et al. 2008). This phenomenon also has been shown in other studies (De Barro et al. 2011). More specimens are required to better conceptualize the distribution of the putative species in the country.

DNA polymorphism analysis showed the highest rate of polymorphism among the sequences reported by Shoorcheh et al. (2008) (Table 2). This might be because of a relatively wider sampling span covering four provinces. In contrast, the sequences from a single province showed a relatively low diversity (Table 2). Furthermore, no host-dependent genetic diversity was found among the specimens collected from Khuzestan as B. tabaci individuals on different hosts shared identical nucleotide sequences. These results were consistent with the results of previous studies (Kareem et al. 2020, Masood et al. 2017, Shahbazi et al. 2014, Shoorcheh et al. 2008).

In Iran, control of B. tabaci is highly dependent on chemical practices; however, resistance has developed to at least six groups of pesticides (Basij et al. 2017). The alternative control practice, biocontrol, also has been challenged by the symbiont-mediated immune responses of B. tabaci individuals to their parasitoid (Mahadav et al. 2008). Considering the presence of subcladeB2 individuals, and possibly its symbionts, this might lead to the development of resistance to chemical compounds and biocontrol agents. Further experiments for detection of the symbionts are required to estimate the risk of pesticide/parasitoid-resistance, virus transmission, and insect fitness.

Acknowledgments

The authors thank the Vice-Consular Office for Research and Technology in Agricultural Sciences and Natural Resources University of Khuzestan. This study was supported by Agricultural Sciences and Natural Resources University of Khuzestan Grant 971.51.

References Cited

  • Basij, M., Talebi K., Ghadamyari M., Hosseininaveh V. and SalamiS.A. 2017. Status of resistance of Bemisia tabaci (Hemiptera: Aleyrodidae) to neonicotinoids in Iran and detoxification by cytochrome P450-dependent monooxygenases.Neotrop. Entomol.46: 115124.
  • Basu, A.N. 2019. Bemisia tabaci (Gennadius): Crop Pest and the Principal Whitefly Vector of Plant Viruses,
    1st ed
    .
    CRC Press
    ,
    Boca Raton, FL
    .
  • De Barro, P.J. and BoykinL.M. 2013. Global Bemisia dataset release version 31 Dec 2012. v1.
    CSIRO Data Collection
    ,
    Australia
    .
  • De Barro, P.J., Liu S.S., Boykin L.M. and DinsdaleA.B. 2011. Bemisia tabaci: a statement of species status.Annu. Rev. Entomol.56: 119.
  • Gilbertson, R.L., Batuman O., Webster C.G. and AdkinsS. 2015. Role of the insect supervectors Bemisia tabaci and Frankliniella occidentalis in the emergence and global spread of plant viruses.Ann. Rev. Virol.2: 6793.
  • Hodges, G.S. and EvansG.A. 2005. An identification guide to the whiteflies (Hemiptera: Aleyrodidae) of the southeastern United States.Florida Entomol.88: 518534.
  • Kanakala, S. and GhanimM. 2019. Global genetic diversity and geographical distribution of Bemisia tabaci and its bacterial endosymbionts.PloS One14: e0213946.
  • Kareem, A.A., Logan S.A., Port G. and WolffK. 2020. Bemisia tabaci in Iraq: Population structure, endosymbiont diversity and putative species.J. Appl. Entomol.144: 297307.
  • Kunz, D., Tay W.T., Elfekih S., Gordon K.H.J. and De BarroP.J. 2019. Take out the rubbish–removing NUMTs and pseudogenes from the Bemisia tabaci cryptic species mtCOI database. BioRxivdoi.org/10.1101/724765.
  • Mahadav, A., Gerling D., Gottlieb Y., Czosnek H. and GhanimM. 2008. Parasitization by the wasp Eretmocerus mundus induces transcription of genes related to immune response and symbiotic bacteria proliferation in the whitefly Bemisia tabaci.BMC Genomics9: 342.
  • Masood, M., Amin I., Hassan I., Mansoor S., Brown J.K. and BriddonR.W. 2017. Diversity and distribution of cryptic species of the Bemisia tabaci (Hemiptera: Aleyrodidae) complex in Pakistan.J. Econ. Entomol.110: 22952300.
  • Mugerwa, H., Seal S., Wang H.L., Patel M.V., Kabaalu R., Omongo C.A., Alicai T., Tairo F., Ndunguru J., Sseruwagi P. and ColvinJ. 2018. African ancestry of New World, Bemisia tabaci-whitefly species.Sci. Rep.8: 111.
  • Nei, M. 1987. Molecular Evolutionary Genetics.
    Columbia Univ. Press
    ,
    New York
    .
  • Pinheiro-Lima, B., Pereira-Carvalho R.C., Alves-Freitas D.M., Kitajima E.W., Vidal A.H., Lacorte C., Godinho M.T., Fontenele R.S., Faria J.C., Abreu E.F. and VarsaniA. 2020. Transmission of the Bean-Associated Cytorhabdovirus by the Whitefly Bemisia tabaci MEAM1.Viruses.12: 1028.
  • Shahbazi, M., Behjatnia S.A.A., Alichi M., Roumi V. and IzadpanahK. 2014. Biotypes of Bemisia tabaci from Fars Province.Iran. J. Plant Prot. Sci.45: 201211.
  • Shoorcheh, H.R., Kazemi B., Manzari S., Brown J.K. and SarafraziA. 2008. Genetic variation and mtCOI phylogeny for Bemisia tabaci (Hemiptera, Aleyrodidae) indicate that the ‘B'biotype predominates in Iran.J. Pest Sci.81: 199.
Fig. 1
Fig. 1

Circular cladogram created using multiple alignments of 56 mitochondrial cytochrome oxidase 1 sequences from Bemisia tabaci specimens all around the world. The phylogenetic tree was constructed using maximum likelihood and bootstrap (100 replicates) by CLC Main Work-bench 7.6.2. The four Iranian B. tabaci haplotypes B obtained in this study are marked with an asterisk. The MEAM1 and MEAM2 sequences were obtained from GenBank database. The corresponding sequence of Trialeurodes abutilonea were used as outgroups (Og). All other sequences are available at CSIRO (De Barro and Boykin 2013).


Contributor Notes

Corresponding author (email: zandi@asnrukh.ac.ir).
Received: 01 Jan 2021
Accepted: 31 Jan 2021
  • Download PDF