Characterization of the complete mitochondrial genome of Cerogria popularis Borchmann, 1936 (Coleoptera: Tenebrionidae)

Abstract In this study, we sequenced the complete mitochondrial genome (mitogenome) of Cerogria popularis Borchmann, 1936 based on the Illumina platform. The circular DNA molecule is 16,175 bp in size, including the 37 typical animal mitochondrial genes and a non-coding control region. All 37 genes are arranged in the same order as the previously reported most mitogenomes of Tenebrionidae. All PCGs initiate with standard start codon ATN (ATA/T/G/C) except cox1 with the special start codon AAT. Most PCGs terminate with TAA/G, whereas cox1, atp6, nad5, and nad4 end with its incomplete form T-. All the 22 tRNAs have the typical clover-leaf structure except for trnS1. The rrnL and rrnS genes are 1,250 and 737 bp in length, with an AT content of 82.6 and 84.5%, respectively. The phylogenetic tree supports the monophyly of the included tenebrionid subfamilies represented by more than one species. Furthermore, the sister relationship between Lagriinae and other tenebrionid subfamilies is recovered.

The subfamily Lagriinae belongs to the family Tenebrionidae of Coleoptera, and to date there are more than 2,200 species (subspecies) described worldwide (Zhou and Chen 2014). Its range and taxonomic status have been debated throughout the history of research (Zhu 2003). This study firstly sequenced the complete mitochondrial genome (mitogenome) of Cerogria popularis Borchmann, 1936, which will be useful for exploring the phylogenetic status of Lagriinae.
The mitogenome of C. popularis (GenBank accession no. MZ699994) is 16,175 bp in size, which comprises 13 proteincoding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding control region (length: 1,692 bp). All 37 genes are arranged in the same order as the previously sequenced most mitogenomes of Tenebrionidae (Hong et al. 2020). The overall nucleotide composition of the newly sequenced mitogenome is 41.9% T, 10.5% C, 38.3% A, and 9.3% G, with an AT content of 80.2%. All PCGs initiate with standard start codon ATN (ATA/T/G/C) except cox1 with the special start codon AAT. Most PCGs terminate with TAA/G, whereas cox1, atp6, nad5, and nad4 end with its incomplete form T-. All tRNAs exhibit a typical cloverleaf secondary structure except for trnS1, in which the dihydrouridine (DHU) arm is replaced by a simple loop. The rrnL and rrnS genes are 1,250 and 737 bp in length, with an AT content of 82.6 and 84.5%, respectively. The control region is located between rrnS and trnI, which is 1,692 bp in length.
In this study, based on the nucleotide data of 13 PCGs from 17 Tenebrionidae species and two outgroup taxa from Cerambycidae, we reconstructed the ML phylogenetic tree (Figure 1). The phylogenetic analyses strongly support the monophyly of the included tenebrionid subfamilies represented by more than one species (Alleculinae, Stenochiinae, Tenebrioninae, and Pimeliinae) (BS 95), similar to the recent study (Hong et al. 2020), but our taxon sample is much larger. In Tenebrionidae, the relationships among included subfamilies are inferred as (Lagriinae þ (Pimeliinae þ (Tenebrioninae þ (Stenochiinae þ (Alleculinae þ Diaperinae))))). Furthermore, the sister relationship between Lagriinae and other tenebrionid subfamilies is recovered (BS ¼ 100).

Ethical approval
Experiments were performed following the recommendations of the Ethics Committee of Guizhou Academy of Forestry. These policies were enacted according to the Chinese Association for the Laboratory Animal Sciences and the Institutional Animal Care and Use Committee (IACUC) protocols.

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