The Complete Mitochondrial Genome of endemic giant tarantula, Lyrognathus crotalus (Araneae: Theraphosidae) and comparative analysis

The complete mitochondrial genome of Lyrognathus crotalus is sequenced, annotated and compared with other spider mitogenomes. It is 13,865 bp long and featured by 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and a control region (CR). Most of the PCGs used ATN start codon except cox3, and nad4 with TTG. Comparative studies indicated the use of TTG, TTA, TTT, GTG, CTG, CTA as start codons by few PCGs. Most of the tRNAs were truncated and do not fold into the typical cloverleaf structure. Further, the motif (CATATA) was detected in CR of nine species including L. crotalus. The gene arrangement of L. crotalus compared with ancestral arthropod showed the transposition of five tRNAs and one tandem duplication random loss (TDRL) event. Five plesiomophic gene blocks (A-E) were identified, of which, four (A, B, D, E) retained in all taxa except family Salticidae. However, block C was retained in Mygalomorphae and two families of Araneomorphae (Hypochilidae and Pholcidae). Out of 146 derived gene boundaries in all taxa, 15 synapomorphic gene boundaries were identified. TreeREx analysis also revealed the transposition of trnI, which makes three derived boundaries and congruent with the result of the gene boundary mapping. Maximum likelihood and Bayesian inference showed similar topologies and congruent with morphology, and previously reported multi-gene phylogeny. However, the Gene-Order based phylogeny showed sister relationship of L. crotalus with two Araneomorphae family members (Hypochilidae and Pholcidae) and other Mygalomorphae species.