The complete mitochondrial genome of Melon thrips, Thrips palmi (Thripinae): Comparative analysis

The melon thrips, Thrips palmi is a serious pest and vector for plant viruses on a wide range of economically important crops. DNA barcoding evidenced the presence of cryptic diversity in T. palmi and that warrants exhaustive molecular studies. Our present study is on decoding the first complete mitochondrial genome of T. palmi (15,333 bp) through next-generation sequencing (NGS). The T. palmi mt genome contains 37 genes, including 13 Protein coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA (tRNAs), and two control regions (CRs). The majority strand of T. palmi revealed 78.29% A+T content, and 21.72% G+C content with positive AT skew (0.09) and negative GC skew (-0.06). The ATN initiation codons were observed in 12 PCGs except for cox1 which have unique start codon (TTG). The relative synonymous codon usage (RSCU) analysis revealed Phe, Leu, Ile, Tyr, Asn, Lys and Met were the most frequently used amino acids in all PCGs. The codon (CGG) which is assigned to Arginine in most insects but absent in T. palmi. The Ka/Ks ratio ranges from 0.078 in cox1 to 0.913 in atp8. We observed the typical cloverleaf secondary structure in most of the tRNA genes with a few exceptions; absence of DHU stem and loop in trnV and trnS, absence of DHU loop in trnE, lack of T-arm and loop in trnN. The T. palmi gene order (GO) was compared with ancestral GO and observed an extensive gene arrangement in PCGs, tRNAs and rRNAs. The cox2 gene was separated from the gene block ‘cox2-trnL2’ in T. palmi as compared with the other thrips mt genomes, including ancestor GO. Further, the nad1, trnQ, trnC, trnL1, trnV, trnF, rrnS, and rrnL were inversely transpositioned in T. palmi GO. The gene blocks ‘trnQ-trnS2-trnD’ and ‘trnN-trnE-trnS1-trnL1’ seems to be genus specific. The T. palmi mt genome contained 24 intergenic spacer regions and 12 overlapping regions. The 62 bp of CR2 shows the similarity with CR1 indicating a possible duplication. The occurrence of multiple CRs in thrips mt genomes seems to be a derived trait which needs further investigation. Although, the study depicted extensive gene rearrangements in T. palmi mt genome, but the negative GC skew reflects only strand asymmetry. Both the ML and BI phylogenetic trees revealed the close relationships of Thrips with Scirtothrips as compared to Frankliniella. Thus, more mt genomes of the diverse thrips species are required to understand the in-depth phylogenetic and evolutionary relationships.


Introduction
The members of insect order Thysanoptera (commonly known as thrips) are usually tiny, fringe winged and are classified into nine families within two suborders [1]. The family Thripidae is the most diverse family and further divided into four subfamilies (Dendrothripinae, Panchaetothripinae, Sericothripinae and Thripinae). Thrips are widely distributed throughout the world and known by 6154 species, of which 739 species are reported from India [2]. Thrips are one of the major sucking pests and sole transmitters of Tospoviruses (Family Bunyaviridae) on a wide number of agricultural and horticultural crops [3]. Fifteen species are reported as vectors of tospoviruses, of which six species are known from India (Ceratothripoides claratris, Frankliniella occidentalis, Frankliniella schultzei, Scirtothrips dorsalis, Thrips palmi and Thrips tabaci) [4][5]. Due to their economic importance, the accurate identification of these species is the basic necessity; for the study of their disease transmission efficiency in crops, and implementation of the effective management strategies. The morphological identification is time consuming and challenging because of their cryptic behaviour and overlapping geographical distributions. Therefore, modern molecular tools such as DNA barcoding have been applied for thrips identification and phylogenetic studies [6][7][8].
Thrips palmi is commonly known as melon thrips, one of the major pest on agricultural crops and has been reported as vector for four tospoviruses; Calla lily chlorotic spot virus (CCSV), Groundnut bud necrosis virus (GBNV), Melon yellow spot virus (MYSV), Watermelon silver mottle virus (WSMV) [4]. The species is widely distributed and highly polyphagous and is often confused with the Thrips flavus and Thrips alatus [9]. Recent DNA barcoding studies on thrips, revealed four molecular operations taxonomic units (MOTUs) in T. palmi (T. palmi Ia1, T. palmi IIa1, T. palmi Ib1, and T. palmi Ib2) representing multiple cryptic species [8]. Considering these hitches, more and in-depth molecular data on this species is required to understand the cryptic speciation and evolutionary affiliations.
In the present study, we sequenced and characterized the complete mt genomes of the melon thrips, T. palmi under the family Thripidae using the next generation sequencing (NGS) technology and comparative analysis with other thrips mt genomes. The comparisons were based on gene arrangements, nucleotide composition, codon usage, evolutionary rates, and strand asymmetry etc. Further, to infer the phylogenetic relationships, 13 PCGs of T. palmi and other six thrips species mt genomes were analyzed using maximum likelihood (ML) and Bayesian inference (BI). This study will provide a better understanding of comparative mitochondrial genomics T. palmi with thysanopterans and other insects.

Ethics statement
There is no need of specific permission for the collection and usage of the insects in this study because these insects are common pests of crops. Neither endangered nor protected species were involved in the study.

Sample collection, and DNA extraction
The specimens of T. palmi were collected from the Odisha state of India from Eggplant (Solanum melongena). The specimens were morphologically identified by the second author (K. T) with the available taxonomic keys [9], and preserved in absolute ethyl alcohol at −80˚C in Centre for DNA Taxonomy, Molecular Systematics Division, Zoological Survey of India, Kolkata. The DNeasy DNA Extraction kit (Qiagen, Valencia, CA) was used to extract the genomic DNA and the concentration was measured on a quantified by Qubit fluorometer (Thermo Fisher Scientific, MA, USA) using a dsDNA high-sensitivity kit with the standard protocol.

Mitochondrial genome sequencing and assembly
The complete mt genome sequencing, assembly and annotation were carried out at the Genotypic Technology Pvt. Ltd. Bangalore, India (http://www.genotypic.co.in/). Whole genome sequencing (WGS) library was prepared with Illumina-compatible NEXTflex Rapid DNA sequencing kit (BIOO Scientific, Austin, Texas, U.S.A.). The DNA was sheared using a Covaris S2 sonicator (Covaris, Woburn, Massachusetts, USA) to generate approximate fragment size distribution of 200 bp to 600 bp. The fragment size distribution was checked on Agilent TapeStation and subsequently purified using Highprep magnetic beads (Magbio). Purified fragments were end-repaired, adenylated and ligated to Illumina multiplex barcode adaptors. Adapter-ligated DNA was purified using Highprep beads and the resulted fragments were amplified for eight cycles of PCR using Illumina-compatible primers. The final PCR product was purified, followed by library quality control check. Illumina-compatible sequencing library was quantified by Qubit fluorometer (Thermo Fisher Scientific, MA, USA) and its fragment size distribution were analyzed by Agilent 2200 Tapestation and Agilent Bioanalyzer. Libraries with adapter contamination were pooled and cleaned up with Highprep magnetic beads (Magbio) and then sequenced on Nextseq 500 150X2 chemistry. Raw Sequences were trimmed and filtered using NGS Toolkit. The adapter contamination and low-quality reads with base N's or more than 70% of the bases with a quality score <20 also had been removed. The acquired high quality~18 million reads were screened out using the Burrows-Wheeler Alignment (BWA) tool [20]. Out of 18 million reads, 0.10% (~1.8 million) of the reads got aligned, then assembled with SPAdes 3.9.0 [21], using default parameters considering T. imaginis mitochondrial genome (AF335993) as reference contig. The aligned reads were considered for the denovo mitochondrial genome of T. palmi.

Genome annotation, visualization, and comparative analysis
The assembled mt genome was annotated by using MITOS web-server (http://mitos.bioinf. uni-leipzig.de/index.py) to estimate the position of PCGs, tRNAs, rRNAs and their secondary structures. The boundaries of PCGs and rRNAs was confirmed manually by nucleotidenucleotide BLAST (BLASTn), protein-protein BLAST (BLASTp), and open reading frame finder (ORF finder) in the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/orffinder/). The PCGs were translated into putative proteins using the invertebrate mitochondrial DNA genetic code. The ClustalX program was used to assign the initiation and termination codons in comparison with other thrips reference sequences [22]. MEGAX was used for alignment of the homologous sequences of T. palmi with other thrips species [23]. The complete annotated T. palmi mt genome was prepared using the Sequin submission tool (http://www.ncbi.nlm.nih.gov/Sequin/), for acquiring the accession number from GenBank database. The circular representation of T. palmi mt genome was drawn by CGView online server (http://stothard.afns.ualberta.ca/cgview_ server/) with default parameters [24]. The assembled T. palmi mt genome was compared with the other thrips mt genomes to calculate the nucleotide composition, Relative Synonymous Codon Usage (RSCU), AT-GC skew, non-synonymous (Ka) and synonymous (Ks) substitutions, gene rearrangements, secondary structure of tRNAs and CRs etc. The nucleotide composition and RSCU were determined using MEGAX [23]. The following formula was used to calculate the skew: AT skew = (A−T) / (A+T) and GC skew = (G−C) / (G+C) [25]. The sequence substitution saturation analysis of PCGs was calculated by DAMBE5 software [26]. The overlapping and intergenic spacer regions of thrips mt genome genes were compared in terms of length and locations. Further the homology of CRs in T. palmi and other Thrips mt genomes were determined through sequence alignment using Clustal Omega [22]. The nucleotide sequences of each PCG were aligned based on the amino acid sequences in TranslatorX [27] with MAFFT algorithm using GBlocks parameters. The ratios of non-synonymous substitutions (Ka) and synonymous (Ks) substitutions were estimated in DnaSP6.0 [28]. The secondary structures of tRNAs were anticipated by the MITOS web server (http://mitos.bioinf.uni-leipzig.de) and confirmed by the tRNAscan-SE (http:// lowelab.ucsc.edu/tRNAscan-SE/) [29] and ARWEN 1.2 [30]. The RNAstructure version 6.0.1 was used to predict the secondary structure of CRs [31]. The gene arrangements, pairwise comparisons of gene orders, and evolutionary pathways of T. palmi mt genome were evaluated by CREx web tool (Common Interval Rearrangement Explorer) [32]. This webbased program uses transpositions, reverse transpositions, reversals and tandem-duplication-random-loss (TDRL) events to find the parsimonious rearrangement events in a phylogenetic hypothesis. CREx analysis based on common intervals was implemented to determine the gene arrangement scenarios between the T. palmi gene order (GO) with the ancestral A. bakeri (Hemiptera) and other thrips species.

Phylogenetic analysis
Six complete mt genomes of five thrips species were retrieved from GenBank on April 1 st , 2018 for phylogenetic inference (S1 Table). The A. bakeri mt genome was used in the dataset as an out-group [33]. The PCGs were individually aligned with the TranslatorX online platform using the MAFFT algorithm [27] with the GBlocks parameters and default settings. The dataset of all PCGs was concatenated using SequenceMatrix v1.7.845 to form 9696 bp dataset [34]. The optimal substitution model was estimated by jModelTest [35]. The dataset was analyzed using maximum likelihood (ML) implemented in RaxML, and Bayesian inference (BI) implemented in MrBayes 3.2 [36]. For likelihood analysis, the bootstrap analysis of 1,000 replicates was performed in RaxML [37] with GTR+G+I as a best fit model. For BI analysis, two simultaneous runs of 12 million generations were conducted for the dataset using GTR+G+I model and trees were sampled in every 1000 generations, with the first 25% discarded as burnin. The BI analysis was stopped after reaching the stationary phase and average standard deviation of split frequencies below 0.01. The phylogenetic tree was visualized and edited using Fig-Tree v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/) [38].

Genome structure and nucleotide composition
The complete mt genome of T. palmi (GenBank accession no. MH253898) was 15,333 base pairs (bp) in length. This is the second largest mt genome in size among all previously generated mt genomes of insect order Thysanoptera (S1 Table). The T. palmi mt genome was characterized by 37 genes, including 13 PCGs, large and small ribosomal RNA genes (rrnL and rrnS), 22 transfer RNA genes (tRNAs) and two control regions (CRs) (Fig 1). Among all genes, 31 genes were detected on the majority (light) strand and six genes on the minority (heavy) strand ( Table 1). The nucleotide composition of T. palmi mt genome was 78.29% A+T content (42.71% A + 35.58% T) and 21.72% G+C content (10.14% G + 11.58% C) ( Table 2). The A+T content was higher in control regions (85.81%) followed by tRNAs (80.26%), rRNAs (80.21%) and PCGs (77.43%). The sequence similarity searches of T. palmi mt genome showed the PCGs are shown as purple arrows, rRNA genes as pink arrows, tRNA genes as peach color arrows and CR regions as gray rectangles. The GC content is plotted using a black sliding window, as the deviation from the average GC content of the entire sequence. The GC-skew is plotted using a colored sliding window (green and orchid color), as the deviation from the average GC skew of the entire sequence. The figure was drawn using the CGView online server (http://stothard.afns.ualberta.ca/cgview_server/) with default parameters. The species photograph was taken by the second author (KT) using Leica Microscope DM1000 with Leica software application suite (LAS EZ) and edited manually in Adobe Photoshop CS 8.0. https://doi.org/10.1371/journal.pone.0199404.g001 Complete mitochondrial genome of T. palmi highest similarity (81%) with T. imaginis followed by F. occidentalis (74%), F. intonsa (73%), S. dorsalis SA1 (72%), S. dorsalis EA1 (71%) and A. obscurus (72%) in GenBank.

Protein coding genes and Relative Synonymous Codon Usage
The T. palmi mt genome was represented by 13 PCGs (atp6, atp8, cox1, cox2, cox3, cytb, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6). The ATN initiation codons (six with ATA, four with ATT and two with ATG) was observed in 12 PCGs except cox1 (TTG). The TTG start codon for cox1 is unique in T. palmi as ATN start codon is observed for cox1 in all other thrips mt genomes (S2 Table). The PCGs were terminated with TAA stop codon except nad4L (TAG) and atp8 with an incomplete stop codon. The incomplete termination codons are common in other insect mt genomes which are assumed to be recovered by the post transcriptional polyadenylation [39][40][41]. The average A+T content of T. palmi PCGs was 77.43% with the  10.52% Guanine (G), and 12.52% Cytosine (C) (S1 Fig, S3 Table). Further, the comparative analysis of PCGs in all seven mt genomes of six thrips species showed that the adenine (A) and thiamine (T) composition was higher than guanine (G) and cytosine (C) ( Table 2). The RSCU data analysis of 3,676 codons in 13 PCGs of T. palmi revealed that Phenylalanine (Phe), Leucine (Leu), Isoleucine (Ile), Tyrosine (Tyr), Asparagine (Asn), Lysine (Lys), and Methionine (Met) were the most frequently used amino acids (S2 Fig,  S4 Table). Comparative RSCU analysis of seven thrips mt genomes showed that Phe, Leu, Ile, Tyr, Asn, Lys and Met were most frequent amino acids with TTT (Phe), TTA (Leu), ATT (Ile), TAT (Tyr), AAT (Asn) and AAA (Lys), ATA (Met) were the most frequently used codons. In contrast, almost all the frequently used codons were ended with A/T, which may lead to the A and T bias in thrips species mt genomes. The codon CGG (Arg) and GCG (Ala) were absent in T. palmi and S. dorsalis SA1 respectively, while these codons were present in other thrips mt genomes. Both of the missing codons were preferred to end with 'G' in the third codon position ( Table).

Non-synonymous and synonymous substitutions
The non-synonymous and synonymous substitutions (Ka/Ks) ratio is an indicator for investigating the selective pressure and evolutionary relations of the homogenous or heterogeneous species [11]. It was reported that, (i) the Ka/Ks>1 for positive selection, (ii) Ka/Ks = 1 for neutral mutation, and (iii) Ka/Ks<1 for negative selection [42][43][44]. The Ka/Ks ratio ranges from 0.078±0.02 in cox1 to 0.913±0.40 in atp8 gene and the resulted following order: cox1<cox3<cytb<cox2<atp6<nad1<nad5<nad4L<nad4<nad3<nad6<nad2<atp8. This result indicated that the 13 PCGs of all thrips mt genomes including T. palmi were evolving under purifying selection (Fig 2A). Comparative analysis of the Ka/Ks ratio among 13 PCGs of thrips species showed that atp8 (1.7) and nad6 (1.3) genes are evolving under positive/relaxed selection with reference to F. intonsa and S. dorsalis SA1 respectively (Fig 2B). Further, the Ka/ Ks was <1 for the remaining PCGs of T. palmi with reference to other thrips species, suggested that the mutations were replaced by synonymous substitutions. The lowest Ka/Ks ratio was observed in cox1 gene representing less variations in amino acids and hence had been evidenced as a potential molecular marker for thrips species identification and phylogenetic analysis [6][7][8].

Ribosomal RNA and transfer RNA
The T. palmi mt genome comprises two rRNA genes as observed in other insect mt genomes. The large ribosomal gene (rrnL) was 1111 bp long, and located between trnV and cox1. Further, the small ribosomal gene (rrnS) was 728 bp long, and located between trnF and atp8 gene. The A+T content of two rRNA was 80.04% in T. palmi which is the highest in comparison with other thrips mt genomes. Both AT skew (0.18) and GC skew (0.08) of T. palmi were positive, that is also similar to other previously sequenced thrips mt genomes. The locations of rrnL and rrnS were upstream of cox1 and atp8 gene, the arrangements seem to be conserved in insect order Thysanoptera. The T. palmi mt genome contained 22 tRNAs (ranging from 57 to 68 bp in length) with a total length of 1,393 bp. Nineteen tRNA genes were coded by the majority strand and three (trnY, trnP and trnH) by the minority strand. The A+T content of tRNAs was 85.81% with positive AT skew (0.05) and GC skew (0.3) ( Table 2). Most of the tRNA showed the typical cloverleaf secondary structure; absence of DHU stem and loop was observed in trnV and trnS; absence of DHU loop in trnE; lack of TCC loop in trnN; variation in TCC arm and loop in other tRNAs (S5 Fig). The absence of DHU stem and loop in trnV was consistent in all the thrips mt genomes. The mismatched base pairs (G-U wobble pairs) were observed in seven tRNAs of T. palmi mt genome. These wobble mismatches were observed in trnL2 and trnL1 (in DHU arm), trnA, trnS1, and trnL2 (in acceptor arm), trnS2 and trnT (in TCC arm).

Overlapping and intergenic spacer regions
The T. palmi mt genome had 24 intergenic spacer regions with a total of 520 bp, varying from 1 to 99 bp in length. There are 14 major intergenic spacers of >10 bp in length were observed  trnP and CR2, cytb and trnI, trnA and trnM, trnD and trnS2. Further, the comparative studies showed that the longest intergenic spacer 150 bp were between trnH and nad4 in S. dorsalis EA1. The T. palmi mt genome contained 11 overlapping regions with a total length of 78 bp. The comparative analysis showed that, the highest (15 overlapping regions of 86 bp) were observed in S. dorsalis SA1. The smallest overlapping region (1 bp) was observed in five positions in T. palmi: nad3 and cox1, trnK and trnG, trnF and trnA, atp8 and rrnS, trnQ and  atp6. The largest overlapping region (50 bp) was observed between nad2 and trnW in T. palmi. However, the largest overlapping region (66 bp) was observed between rrnL and trnS2 in F. occidentalis in comparative studies (S6 Table).

Control regions
The CRs is usually characterized by five conserved elements, such as (i) a polyT stretch at the 5 0 end; (ii) a [TA(A)]n-like stretch (iii) a stem and loop structure (iv) a TATA motif and a G (A)nT motif flanking the stem and loop structure (v) a G+A-rich sequence downstream of the stem and loop structure [45][46][47]. Both CRs of T. palmi were detected with four conserved elements except G+A-rich sequence (Fig 3). The control region (CR) in mt genomes is most crucial region, which regulates transcription and replication [48]. Duplicate CRs has been documented in many insect species, including thrips [15][16][17][18]. The study resulted in the duplication of CRs in T. palmi mt genome; 489 bp CR1 lied between trnR and nad5, 138 bp CR2 between the trnL1 and trnP. Further, we observed 62 bp sequence similarity of CR2 with CR1 indicating a possible duplication. The total length of CR1 and CR2 was 627 bp, which was higher than A. obscurus (145 bp) and F. occidentalis (595bp). T. palmi had two CRs, similar to T. imaginis, and S. dorsalis EA1. The presence of multiple CRs in thrips species that are serious pests and vectors for tospoviruses, seems to be a genomic apomorphy over other common Complete mitochondrial genome of T. palmi thrips species and the ancestral species. However, more thrips mt genomes are needed to be studied to confirm, whether the duplication of CR in thrips species is linked with the association of tospoviruses.

Gene arrangement
The gene arrangements can be categorized by following aspects, (i) transpositions, (ii) inversions, (iii) inverse transpositions, and (iv) the tandem duplication random loss operation (TDRL) [49][50][51]. The gene arrangement of T. palmi was evaluated by comparing the common intervals with A.bakeri GO, which was assumed to be the putative ancestor of hexapods [19,33]. The CREx identified seven operations in the evolution of gene arrangement in T. palmi, including three inversions and four TDRLs with two sets of alternative scenarios (S6 Fig). The CREx analysis detected inversions of trnF, trnC and 'nad1-rrnS' gene block, in both scenarios. All the genes in T. palmi mt genome were encoded by the majority strand except three PCGs (nad5, nad4, nad4L) and three tRNAs (trnH, trnY and trnP). T. palmi GO showed intense gene rearrangements of 11 PCGs, 22tRNAs, and two rRNAs as compared with the ancestral GO. Most of the genes were transpositioned, while eight genes (nad1, trnL1, trnF, trnQ, trnC, trnV, rrnS, and rrnL) were inversely transpositioned as compared with the ancestral GO. Further, T. palmi GO was also compared other thrips mt genomes GO and resulted that nad2 was conserved between trnY and trnW in all thrips species. The cox2 was separated in T. palmi GO from 'cox2-trnL2' gene block, which was conserved in other thrips species including ancestral GO (Fig 4). The comparison of T. palmi GO with other thrips mt genomes further revealed trnD and trnR were translocated and separated from 'trnD-cox3' gene block in T. plami and T. imaginis. The gene blocks 'cox3-nad2' and 'atp6-nad5' have been identified as the most frequently rearranged position in thrips mt genomes, while 'nad5-trnH-nad4-nad4L' seems to be an ancestral gene block. The tRNAs gene blocks 'trnQ-trnS2-trnD' and 'trnN-trnE-trnS1-trnL1' might be specific to genus Thrips and Frankliniella respectively.

Strand asymmetry
The strand asymmetry is a remarkable feature of mitochondrial genomes denoted by the AT skew and GC skew on the majority strand (encoded with more genes) [52][53][54]. Positive AT skew (A>T) and negative GC skew (C>G) is usually observed in most of the insect's mitochondrial genome with few exceptions; where strand asymmetry is reversed, indicating A<T and C<G on the majority strand [46][47]50]. A recent study on the mechanism of strand asymmetry in insects confirmed that reversal of strand asymmetry was due to the inversion of replication origin in CRs [54]. We calculated and compared the AT and GC skews of T. palmi with other thrips mt genomes. The T. palmi mt genome showed positive AT skew (0.09) and negative GC skew (-0.07) found to be similar with most of the thrips mt genomes. The AT skew of complete mt genomes in other thrips species, ranged from -0.02 (A. obscurus) to 0.15 (T. imaginis), while the GC skew varies from -0.11 (T. imaginis) to 0.01 (F. occidentalis). The occurrence of Guanine (Gs) is more than Cytosine (Cs) with positive GC skew was observed in A. obscurus revealed the reversal of strand asymmetry (Table 2). We have also observed and compared the skew value of all PCGs in thrips mt genomes including T. palmi. The present study depicted that duplication of CRs and structural elements might regulate the extensive gene rearrangements in thrips mt genomes.

Phylogenetic analysis
The Maximum likelihood (ML) and Bayesian Inference (BI) phylogenetic trees were constructed by using 13 PCGs. The A. bakeri mt genome was used as the out-group in the phylogenetic analysis. The phylogenetic trees generated using both the methods resulted similar topologies (Fig 5). The tree clustering revealed that species under genus Frankliniella and Thrips, were clustered under the respective genus clade. The phylogenetic analysis showed that genus Thrips is more closely related to genus Scirtothrips as compared to genus Frankliniella. The genus Thrips and Frankliniella forms two large genus groups 'Thrips genus groups' and 'Frankliniella genus group' respectively, which were supposed to be closely related based on Complete mitochondrial genome of T. palmi the assumed homology for paired ctenidia on abdominal segments V-VIII [55]. However, Mound 2002 suggested that, these two genus-groups are not closely related based on the chaetotaxy of the abdomen and head [55]. Further, the close relationship between Scirtothrips and Thrips cannot be supported with their morphology. Till date, the comprehensive mt genomes data of thrips species is in its early stages and thus warrant generation of more mt genomes data of diverse thrips species from different hierarchical level to understand the phylogenetic and evolutionary relationships among them. Complete mitochondrial genome of T. palmi S1