The complete mitochondrial genome of Aeolothrips indicus Bhatti, 1964 (Thysanoptera: Thripidae)

Abstract Here, we have generated the complete mitochondrial sequence of Aeolothrips indicus Bhatti, 1964. So far, this is the first largest mitogenome with 17,042 bp length in order Thysanoptera. It includes 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes along with three non-coding regions. AT composition of A. indicus is 72.5% (37.7% A and 34.8% T) and GC 27.5% (15.6% C and 11.9% G). The constructed phylogeny revealed the monophyly of family Aeolothripidae in the order Thysanoptera. The data would provide further insight into the evolution and phylogeny of the order Thysanoptera.

The family Aeolothripidae of order Thysanoptera comprises of less than 4% of thrips diversity worldwide, including 218 extant species in 24 genera (ThripsWiki 2021). Most of species in this family are phytophagous in nature, while a few are monophagous like Indothrips bhushani Bhatti is restricted to neem flowers (Azadirachta indica) for a very short period of time in India (Tyagi et al. 2008). They are also known to be obligate or facultative predator on minute arthropods (Mound and Marullo 1998). The genus Aeolothrips was first described by Haliday in 1836 with type Aeolothrips albicincta Haliday. This genus is known by 113 species across the world (ThripsWiki 2021), of which seven species reported from India (Tyagi and Kumar 2016). Aeolothrips indicus was described by Bhatti in 1964(Bhatti 1964, and also known to be endemic to India. This is the second mitogenome for the family Aeolothripidae and first largest genome in order Thysanoptera. Specimen of A. indicus was collected from the flowers and leaves of Mangifera indica in March 2020 from Nainital (29.28N, 79.27E), Uttarakhand, India. The nondestructive DNA extraction from a single specimen was done with the DNeasy Blood and Tissue Kit (QIAGEN, Hilden, Germany). Voucher specimen with registration no. 11139/H17 was deposited in the National Zoological Collections (NZC) at the Centre for DNA Taxonomy, Molecular Systematics Division, Zoological Survey of India, Kolkata, India (collection in-charge, Dr. Vikas Kumar, kumar.vikas@zsi.gov.in). The sequencing was done in Illumina platform (NovaSeq 6000) with 2 Â 150 base pair chemistry and the raw reads were assembled by GetOrganelle software (version 1.7.4) (Jin et al. 2020) using COI sequence of this species as a seed. The annotation of the gene boundaries was done by MITOS Web Server (Bernt et al. 2013), ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder), BLASTn, and BLASTp. Further, these boundaries were confirmed with the alignment of previously available thrips mitogenomes.
A. indicus (17,042 bp) mitogenome includes 37 genes with 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and three putative control regions. This is the second mitogenome in the order Thysanoptera which was reported with three control regions (Tyagi et al. 2020). Most of the genes were located on the majority strand except nad5, nad4, nad4L, trnH, and trnP. The AT content of the genome was 72.5% (37.7% of A and 34.8% of T) and GC 27.5% (11.9% of G and 15.6% of C). atp8, cox2, nad2, nad4L were used as ATA start codon; cox1, cox3, nad1, nad4, nad5, nad6 by ATT; and atp6, cytb, nad3 by ATG. TAA stop codon was used by all the PCGs with few exceptions, like TAG stop codon was used by atp8, nad2, nad4L, nad5. tRNAs were ranging from 74 bp (trnR) to 61 bp (trnS1) with typical cloverleaf secondary structure. The length of rrnL and rrnS was 1325 bp and 804 bp, respectively. Five overlapping regions (1-10 bp with a total of 18 bp) and 23 intergenic spacer regions (1-78 bp with a total of 330 bp) were detected in the mitogenome.
The TranslatorX tool (Abascal et al. 2010) was used for the PCGs alignment and SequenceMatrix v1.8 (Vaidya et al. 2011) for PCGs concatenation. The GTR þ IþG was detected as a best fit model in PartitionFinder version 2.1.1 with BIC criterion (Lanfear et al. 2017). Mr.Bayes ver. 3.2 (Ronquist et al. 2012) was used for the construction of the Bayesian inference (BI) and IQ tree web server (http://iqtree.cibiv.univie.ac.at/) for maximum-likelihood (ML) phylogenetic trees. Alloeorhynchus bakeri (Hemiptera) sequence (GenBank accession HM235722) was used as an out group for the phylogeny. The similar topology was revealed by both phylogenetic methods and their support values were superimposed on a BI tree in Figure 1. Both the phylogeny reflected the monophyly of the family Aeolothripidae with high support.

Disclosure statement
No potential conflict of interest was reported by the author(s).

Funding
This work was financially supported by Zoological Survey of India (ZSI) in-house project, 'National Faunal Genome Resources (NFGR)' and SERB funded project 'Delimiting Species Boundaries in Pest and Vector thrips (Thysanoptera: Thripidae) from India'. This work was also supported by Science and Engineering Research Board.

Data availability statement
The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at under the accession no.