
==== Front
PLoS One
PLoS One
plos
PLOS ONE
1932-6203
Public Library of Science San Francisco, CA USA

10.1371/journal.pone.0310559
PONE-D-24-09889
Research Article
Biology and Life Sciences
Biochemistry
Bioenergetics
Energy-Producing Organelles
Mitochondria
Biology and Life Sciences
Cell Biology
Cellular Structures and Organelles
Energy-Producing Organelles
Mitochondria
Biology and life sciences
Biochemistry
Nucleic acids
RNA
Non-coding RNA
Transfer RNA
Biology and Life Sciences
Zoology
Entomology
Insects
Beetles
Biology and Life Sciences
Organisms
Eukaryota
Animals
Invertebrates
Arthropoda
Insects
Beetles
Biology and Life Sciences
Zoology
Animals
Invertebrates
Arthropoda
Insects
Beetles
Biology and Life Sciences
Evolutionary Biology
Evolutionary Systematics
Phylogenetics
Phylogenetic Analysis
Biology and Life Sciences
Taxonomy
Evolutionary Systematics
Phylogenetics
Phylogenetic Analysis
Computer and Information Sciences
Data Management
Taxonomy
Evolutionary Systematics
Phylogenetics
Phylogenetic Analysis
Biology and Life Sciences
Genetics
Genomics
Animal Genomics
Invertebrate Genomics
Biology and Life Sciences
Zoology
Entomology
Insects
Biology and Life Sciences
Organisms
Eukaryota
Animals
Invertebrates
Arthropoda
Insects
Biology and Life Sciences
Zoology
Animals
Invertebrates
Arthropoda
Insects
Biology and Life Sciences
Biochemistry
Nucleotides
Anticodons
Biology and life sciences
Biochemistry
Nucleic acids
RNA
Non-coding RNA
Transfer RNA
Anticodons
Biology and Life Sciences
Evolutionary Biology
Evolutionary Systematics
Phylogenetics
Animal Phylogenetics
Biology and Life Sciences
Taxonomy
Evolutionary Systematics
Phylogenetics
Animal Phylogenetics
Computer and Information Sciences
Data Management
Taxonomy
Evolutionary Systematics
Phylogenetics
Animal Phylogenetics
Biology and Life Sciences
Zoology
Animal Phylogenetics
The first report of complete mitogenomes of two endangered species of genus Propomacrus (Coleoptera: Scarabaeidae: Euchirinae) and phylogenetic implications
Mitogenomes of two endangered species of genus propomacrus and phylogenetic implications
Yi Chuanhui Conceptualization Funding acquisition Project administration Supervision Writing – review & editing 1
Shu Xu Formal analysis Methodology Writing – original draft 2
Wang Lingmin Formal analysis Methodology Writing – original draft 2
Yin Jing Formal analysis Methodology Writing – original draft 2
Wang Youhui Data curation Formal analysis 2
Wang Yuchen Formal analysis 2
Zhang Honghui Formal analysis 2
He Qiuju Conceptualization Project administration Supervision Writing – review & editing 2 *
https://orcid.org/0000-0002-8752-7221
Zhao Min Conceptualization Project administration Supervision Writing – review & editing 3 4 *
1 Yunnan Institute of Biological Diversity, Southwest Forestry University, Kunming, Yunnan Province, P.R. China
2 Department of Forest Conservation, College of Biodiversity Conservation, Southwest Forestry University, Kunming, Yunnan Province, P.R. China
3 Division of Resource Insects, Institute of Highland Forest Science, Chinese Academy of Forestry, Kunming, Yunnan Province, P.R. China
4 Key Laboratory of Breeding and Utilization of Resource Insects, National Forestry and Grassland Administration, Kunming, Yunnan Province, P.R. China
Bhardwaj Pankaj Editor
Central University of Punjab, INDIA
Competing Interests: The authors have declared that no competing interests exist.

* E-mail: 13678797631@163.com (QH); mzhao@caf.ac.cn (MZ)
18 9 2024
2024
19 9 e031055911 3 2024
3 9 2024
© 2024 Yi et al
2024
Yi et al
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

To understand the mitochondrial genome structure of two endangered and long-armed scarab beetles, Propomacrus davidi and Propomacrus bimucronatus, their complete mitogenomes were sequenced for the first time in this study. The complete mitogenomes of P. davidi and P. bimucronatus were 18, 042 bp and 18, 104 bp in length, respectively. The gene orders of their mitogenomes were highly consistent with other Coleopteran species, and the typical ATN was used as the start codon in most protein coding genes. The incomplete stop codon T was used in cox1, cox2, and nad5, and TAN was used as a complete stop codon in most protein coding genes. All predicted tRNAs could form a typical cloverleaf secondary structure, except that trnS1 lacked the dihydrouridine arm. Based on the maximum likelihood and the Bayesian inference methods, phylogenetic trees of 50 species were reconstructed. The results showed that P. davidi, P. bimucronatus, Cheirotonus jansoni and Cheirotonus gestroi clustered in the same branch, and were the most closely related. The results supported that subfamily Euchirinae is a monophyletic group of Scarabaeidae, which was consistent with the morphological classification. These molecular data enriched the complete mitogenome database of Euchirinae, and improved our understanding of the phylogenetic relationship and evolutionary characteristics of these two endangered species.

Yunnan Agricultural Basic Research Joint Project 202301BD070001-004 Yi Chuanhui http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 31660631 Yi Chuanhui This research was funded by Yunnan Agricultural Basic Research Joint Project (grant number 202301BD070001-004) and the National Natural Science Foundation of China (grant number 31660631), awarded to Chuan-Hui Yi. The funding source did not have any influence on the study design, analyses, interpretation of data. Data AvailabilityAll sequences' files are available from the GenBank (https://www.ncbi.nlm.nih.gov/genbank/) database, and the assembled mitochondrial genome sequence and annotation of Propomacrus davidi and P. bimucronatus described here are accessible via GenBank IDs NC_070351 and NC_070352.
Data Availability

All sequences' files are available from the GenBank (https://www.ncbi.nlm.nih.gov/genbank/) database, and the assembled mitochondrial genome sequence and annotation of Propomacrus davidi and P. bimucronatus described here are accessible via GenBank IDs NC_070351 and NC_070352.
==== Body
pmcIntroduction

The genus Propomacrus Newman, 1837 belongs to the subfamily Euchirinae in the order Coleoptera. It only includes three species, Propomacrus bimucronatus (Pallas, 1781), P. davidi Deyrolle, 1874 and P. muramotoae Fujioka, 2007. Propomacrus cypriacus Alexis & Makris, 2002 has been a subspecies of P. bimucronatus by genetic evidence [1]. The species of Propomacrus are giant beetles, and mainly found in the alpine regions of Asia, the Near East and southeastern Europe [1, 2], with high habitat requirements, sensitivity to environmental changes, poor adaptability, weak reproductive capacity, and long larval stage [3–5]. In recent years, due to the intensification of human activity intervention and serious damage to forests, the distribution areas of Propomacrus have declined sharply and some species have become endangered. P. davidi is endemic to China and only distributed in Jiangxi and Fujian provinces. It includes two subspecies, Propomacrus davidi davidi Deyrolle, 1874 and Propomacrus davidi fujianensis Wu, 2008 [6–9]. The beetle listed as one of the Grade II Key Protected Wild Animals based on the endangered status in China (https://www.forestry.gov.cn/main/5461/20210205/122418860831352.html). P. bimucronatus shows an intermittent distribution in Cyprus, Greece, Iran, Israel, Syria, Turkey, and the southern and central Balkans [1, 3, 10]. It was listed as Near Threatened (NT) in the European Red List of Saproxylic Beetles [11]. The adults of P. davidi and P. bimucronatus feed on the sap of trees and rotten fruits. The larvae feed on rotten wood and humus soil, which can accelerate the decomposition of dead wood and litter in forests [1]. They play an important role in the material cycle of forest ecosystems and have considerable ecological value. Meanwhile, the two species are famous ornamental insects because of their peculiar appearance and good ornamental value [12, 13]. Recent studies on the two beetles P. davidi and P. bimucronatus has mainly focused on morphological classification, biology, and ecology [2, 3, 9, 10, 14], only a few studies involved molecular biology [2, 8], and only cox1 gene and 16S of the two species were sequenced [1].

The family Scarabaeidae, within the order Coleoptera, comprises over 30,000 described species, with over 21 subfamilies [15]. These subfamilies exhibit distinct divisions in their primary feeding ecologies, broadly categorized into saprophagous and phytophagous lineages [15, 16]. The former mainly feed on decaying matter, animal dung, etc., including the subfamilies Scarabaeinae and Aphodiinae, while the latter typically consume plant fruits, flowers, pollen, tree sap, wood (including deadwood), leaves, etc. This group includes the four main subfamilies Melolonthinae, Rutelinae, Cetoniinae, and Dynastinae, as well as some less numerous subfamilies such as Euchirinae.

Euchirinae, a subfamily of Scarabaeidae, is an under-studied group with uncertain taxonomic status [3, 10, 17, 18], which includes three genera and 15 species and seven subspecies worldwide, namely 10 species and one subspecies of Cheirotonus Hope, 1841, two species and one subspecies of Euchirus Burmeister and Schaum, 1840, and three species and two subspecies of Propomacrus Newman, 1837. The phylogenetic study of Euchirinae has been conducted poorly, and it has been mainly based on morphological characteristics, rarely using molecular characteristics. Although its monophyly has been confirmed, the classification at the subfamily level remains uncertain [3, 19, 20]. Some scholars believe that Euchirinae is a subfamily of the family Scarabaeidae [10, 21, 22], but the evidence is needs to be more substantial [18, 23]. The subfamily is included in the subfamily Melolonthinae as tribe Euchirini sometimes [3] or in a clade composed of several representatives of Melolonthinae, Rutelinae, Dynastinae, as well as Cetoniinae [6, 23], or an independent family Euchiridae [1, 20]. Therefore, it is very important to develop reliable markers to reconstruct the phylogenetic relationship of the subfamily Euchirinae.

As a complete organelle genome, the mitochondrial genome (mitogenome) has the characteristics of maternal inheritance, small molecular weight, fast evolution rate and low recombination level [24]. Thus, it is a good marker for studying the phylogenetic relationships between species and genera. Therefore, complete mitogenomes have been used as a common molecular marker in insect systematics studies [25, 26]. The mitogenome of insect is a closed circular and double-stranded DNA molecule of 14–20 kb in length [27, 28]. As of May 2024, the number of mitochondrial genomes of Scarabaeidae insects stored in GenBank has exceeded 100, primarily from seven subfamilies: Scarabaeinae, Melolonthinae, Rutelinae, Cetoniinae, Aphodiinae, Dynastinae, and Euchirinae. Mitochondrial gene characteristics of over 30 Scarabaeidae species have been described and used for phylogenetic analysis [6, 29]. Just as mitochondrial genomes of coleopteran insects tend to be conservative [30], with the exception of some species in the subfamilies Dynastinae and Scarabaeinae that show gene rearrangements, the mitochondrial genome structure, base composition, codon usage, and gene arrangement of other Scarabaeidae species exhibit high similarity to other coleopterans [31, 32]. However, only three species of Euchirinae, Euchirus longimanus, Cheirotonus gestroi, and C. jansoni, are recorded in GenBank [6, 20, 29]. More complete mitogenome data of the Euchirinae species are needed for phylogenetic analysis.

In this study, the complete mitochondrial genomes of two species of genus Propomacrus (Coleoptera: Scarabaeidae: Euchirinae) were sequenced and analyzed for the first time. Combining these data with other available mitogenomes from GenBank, the phylogenetic relationships of Scarabaeoidea were reconstructed, which provides new insights into the relationships among Euchirinae and the other major subfamilies or families. These findings will be useful for better conserving endangered groups in the future.

Materials and methods

Specimen materials

The adult specimens of P. davidi davidi and P. bimucronatus were obtained from beetle breeding enthusiast in Kunming, Yunnan Province, China. Their ancestors were roughly mapped to the Longhushan Geopark in Jiangxi, China and Izmir, Turkey, respectively. They were reared from one instar in October 2019 and emerged in August 2020. Adult samples were stored in absolute ethanol and kept in the herbarium of Southwest Forestry University.

Mitochondrial genome sequencing

The total DNA of the samples was extracted from the muscles of the adult’s thorax using a Magnetic Animal Tissue Genomic DNA Kit (DP341) (TIANGEN, Beijing, China). After the quality detection by 1.0% agarose gel electrophoresis, the DNA was sent to Baiai Gene Information Technology Co., Ltd. (Xi’an, China) for sequencing. Whole Genome Shotgun (WGS) was used to construct the library, and next generation sequencing was used to perform Paired-end (PE) sequencing on the Illumina sequencing platform.

Sequence splicing and annotation

High quality next-generation sequencing data were assembled using SPAdes v3.9.0 [33] and contig and scaffold sequences were constructed. After comparing the sequence with the NT library on NCBI BLAST (BLAST v2.2.31+), the mitochondrial sequence of the splicing results was picked out. Combined with the reference sequences C. jansoni (GenBank ID: NC_023246) [34] and C. gestroi (NC_046890) [20], the mitochondrial sequences of splicing results were integrated. Collinearity analysis was performed using mummer v3.1 [35] to determine the positional relationship between contigs and fill the gap between contigs. The results were corrected using Pilon v1.18 software [36] to obtain the final mitochondrial sequence. The complete mitogenome sequence was uploaded to the MITOS website (http://mitos2.bioinf.uni-leipzig.de/index.py) for functional annotation [37]. The secondary structure of tRNA was predicted by tRNA Scan-SE [38] and MITOS web [37], and it was manually corrected by the SVG online editor. The complete mitogenome map was drawn using OGDRAW v1.3.1 [39].

Phylogenetic analyses

Based on the mitogenome sequence of 48 species from four families of Scarabaeoidea, published in the NCBI database, two species from Staphylinidae and Byturidae were selected as outgroups (Table 1). The species selected from the superfamily Scarabaeoidea include 25 species from six subfamilies of the family Scarabaeidae, including Cetoniinae, Dynastinae, Rutelinae, Melolonthinae, Euchirinae, and Scarabaeinae, as well as two species from the subfamilies Hybosorinae and Ceratocanthinae of the family Hybosoridae, two species from the subfamily Passalinae of the family Passalidae, and 19 species from the subfamily Lucaninae of the family Lucanidae. The selected species were used to construct a phylogenetic tree with the two species and outgroups. PhyloSuite v1.2.2 aligned the complete mitochondrial genome sequence [40]. IQ-TREE v1.6.12 [41] and MrBayes v3.2.7a [42] were used to construct phylogenetic trees of 13 PCGs based on the maximum likelihood (ML) and Bayesian inference (BI) methods.

10.1371/journal.pone.0310559.t001 Table 1 Taxa used for phylogenetic analysis in this study.

Superfamily	Family	Subfamily	Species	GenBank ID	Length (bp)	Reference	
Scarabaeoidea	Scarabaeidae	Cetoniinae	Glycyphana fulvistemma	NC_063847	16,701	N/A	
			Gametis jucunda	NC_063846	17,980	N/A	
			Coenochilus striatus	NC_065313	15,480	N/A	
			Trichius succinctus	NC_063849	18,358	N/A	
			Dicronorhina derbyana	OK484300	16,609	[31]	
			Eudicella smithii	OK484302	16,712	[31]	
			Eudicella quadrimaculata	OK484299	16,690	[31]	
			Protaetia brevitarsis	NC_023453	20,319	[43]	
			Osmoderma opicum	NC_030778	15,341	[44]	
			Mecynorhina polyphemus	OK484305	16,422	[31]	
		Dynastinae	Eupatorus gracilicornis	NC_065036	18,391	N/A	
			Eupatorus hardwickei	ON817147	18,494	N/A	
			Eupatorus sukkiti	ON756028	18,445	N/A	
			Megasoma mars	OK484311	16,983	[31]	
			Trichogomphus mongol	OM161966	17,377	N/A	
		Rutelinae	Anomala russiventris	NC_065310	15,601	[45]	
			Popillia japonica	NC_038115	16,541	[6]	
		Melolonthinae	Sophrops subrugatus	NC_065314	16,409	[45]	
			Rhopaea magnicornis	NC_013252	17,522	[46]	
		Euchirinae	Cheirotonus gestroi	MN893347	16,899	[20]	
			Cheirotonus jansoni	NC_023246	17,249	[34]	
			Propomacrus davidi	NC_070351	18,042	This study	
			Propomacrus bimucronatus	NC_070352	18,104	This study	
		Scarabaeinae	Dichotomius schiffleri	NC_039689	14,802	[47]	
			Copris tripartitus	NC_045923	15,457	[32]	
	Hybosoridae	Hybosorinae	Hybosorus sp.	JX412807	12,639	N/A	
		Ceratocanthinae	Ceratocanthus sp.	JX412772	12,579	N/A	
	Passalidae	Passalinae	Leptaulax koreanus	NC_044848	18,730	[48]	
			Ophrygonius sp.	NC_060602	16,195	N/A	
	Lucanidae	Lucaninae	Lucanus cervus	MN_580549	20,109	[49]	
			Lucanus kanoi kanoi	MK238348	17,991	N/A	
			Lucanus chengyuani	MK878514	16,926	[50]	
			Lucanus fortunei	NC_044961	16,591	[51]	
			Lucanus maculifemoratus taiwanus	MK250907	21,174	N/A	
			Macrodorcas seguyi	NC_038212	17,955	[52]	
			Serrognathus platymelus	NC_044096	16,790	[53]	
			Dorcus hansi	NC_043928	18,130	N/A	
			Prosopocoilus astacoides	NC_050851	17,746	N/A	
			Prosopocoilus bulbosus	MK134566	16,642	N/A	
			Prosopocoilus laterotarsus	NC_065362	17,333	[54]	
			Prosopocoilus gracilis	NC_027580	16,736	[55]	
			Rhaetus westwoodii	MG159815	18,131	[56]	
			Kirchnerius guangxii	NC_048957	14,562	[53]	
			Nigidius miwai	NC_063663	18,462	N/A	
			Nigidius sinicus	NC_065111	18,394	N/A	
			Prismognathus prossi	NC_044962	15,984	[51]	
			Figulus binodulus	NC_045102	16,261	[57]	
			Odontolabis cuvera fallaciosa	MF908524	19,614	[58]	
Staphylinoidea	Staphylinidae	Staphylininae	Creophilus maxillosus	NC_061687	17,395	N/A	
Cucujoidea	Byturidae	-	Byturus ochraceus	NC_036267	15,143	N/A	
Note:“-” stands for not available; “N/A” means unpublished.

Results and discussion

Mitogenome composition

The complete mitogenome lengths of P. davidi and P. bimucronatus were 18,042 bp and 18,104 bp, respectively, which were longer than the 16, 899 bp of C. gestroi in the same subfamily [20]. Most metazoan mitogenomes are approximately 16 kb in size [59, 60]. Longer mitochondrial genes may be characteristic of Propomacrus species, which is closely related to the length of control region sequence. The control regions of P. davidi and P. bimucronatus are 765–1,202 bp longer than in C. gestroi and C. jansoni in same subfamily. The A + T content was 69.49% for P. davidi, and 67.91% for P. bimucronatus, with high AT-skew. The gene composition and sequence were consistent with those of ancestral insects, including 13 protein coding genes (PCGs), 22 tRNAs, two rRNAs, and a control region (also known as A + T-rich region or D-loop) [33, 61]. The J-strand generally contains 23 genes (nine PCGs, 14 tRNAs), and the N-strand includes 14 genes (four PCGs, eight tRNAs, and two rRNAs). There was an intergenic spacer or gene overlap between adjacent genes of the two mitogenomes, and no gene rearrangement was found. The number of intergenic spacers in the whole mitogenomes of P. davidi and P. bimucronatus was seven overlapping regions. The sum of their lengths was 34 bp and 36 bp, respectively. The intergenic spacers seem to mid-length in Coleoptera, for example, nearly twice that of the Tenebrionidae species [62] but shorter than Coenochilus striatus [45]. However, the longest intergenic spacers were both 19 bp and located between trnS2 and nad1. There were 15 regions without intervals and overlaps. The number of overlapping region was 15, the sum of length was 40 bp, and the longest overlapping region was 8 bp, which was located between trnC and trnW (Fig 1 and Table 2), consistent with the results of some other Coleoptera species with only slight differences in length, such as the region of 9 bp and 8 bp of Chrysochroa fulgidissima and Coomaniella copipes, respectively [63, 64]. The gene arrangement, nucleotide composition, and codon usage may be used as a piece of very helpful information to assist the classification and infer evolutionary relationships.

10.1371/journal.pone.0310559.g001 Fig 1 Circular maps of the mitogenome of Propomacrus davidi and P. bimucronatus.

The direction of gene transcription is shown as an arrow. Dark blue represents transfer RNA (tRNA), yellow represents nicotinamide adenine dinucleotide (nadh), purple represents cytochrome c oxidase (cox) and cytochrome b (cytb), green represents adenosine triphosphate ATP synthase, and red represents ribosomal RNA (rRNA).

10.1371/journal.pone.0310559.t002 Table 2 Organization of the mitogenomes of Propomacrus davidi and P. bimucronatus.

Gene	Direction	Position	Start codon	Stop codon	Anticodon	Intergenic	
trnI	J	1–68 / 1–67			GAU / GAU	-3 / -3	
trnQ	N	66–134 / 65–133			UUG / UUG	-1 / -1	
trnM	J	134–202 / 133–201			CAU / CAU	0 / 0	
nad2	J	203–1210 / 202–1209	ATT / ATT	TAA / TAA		5 / 5	
trnW	J	1216–1283 / 1215–1282			UCA / UCA	-8 / -8	
trnC	N	1276–1343 / 1275–1341			GCA / GCA	5 / 7	
trnY	N	1349–1435 / 1349–1412			GUA / GUA	-2 / -2	
cox1	J	1434–2967 / 1411–2944		T / T		0 / 0	
trnL2	J	2968–3032 / 2945–3008			UAA / UAA	0 / 0	
cox2	J	3033–3720 / 3009–3696	ATA / ATA	T / T		0 / 0	
trnK	J	3721–3791 / 3697–3767			CUU / CUU	-1 / -1	
trnD	J	3791–3854 / 3767–3830			GUC / GUC	0 / 0	
atp8	J	3855–4010 / 3831–3986	ATT / ATT	TAA / TAA		-7 / -7	
atp6	J	4004–4675 / 3980–4651	ATG / ATG	TAA / TAA		-1 / -1	
cox3	J	4675–5461 / 4651–5437	ATG / ATG	TAA / TAA		0 / 0	
trnG	J	5462–5525 / 5438–5502			UCC / UCC	0 / 0	
nad3	J	5526–5879 / 5503–5856	ATT / ATT	TAA / TAA		-2 / -2	
trnA	J	5878–5940 / 5855–5917			UGC / UGC	0 / 0	
trnR	J	5941–6005 / 5918–5982			UCG / UCG	0 / 0	
trnN	J	6006–6071 / 5983–6048			GUU / GUU	0 / 0	
trnS1	J	6072–6139 / 6049–6116			UCU / UCU	0 / 0	
trnE	J	6140–6202 / 6117–6179			UUC / UUC	-2 / -2	
trnF	N	6201–6265 / 6178–6242			GAA / GAA	0 / 0	
nad5	N	6266–7979 / 6243–7956	ATT / ATT	T / T		1 / 1	
trnH	N	7981–8047 / 7958–8022			GUG / GUG	-1 / -1	
nad4	N	8047–9384 / 8022–9359	ATG / ATG	TAA / TAA		-7 / -7	
nad4L	N	9378–9668 / 9353–9643	ATG / ATG	TAA / TAA		2 / 2	
trnT	J	9671–9734 / 9646–9709			UGU / UGU	-1 / -1	
trnP	N	9734–9800 / 9709–9774			UGG / UGG	1 / 1	
nad6	J	9802–10302 / 9776–10276	ATT / ATC	TAA / TAA		-1 / -1	
cytb	J	10302–11444 / 10276–11418	ATG / ATG	TAG / TAG		-2 / -2	
trnS2	J	11443–11507 / 11417–11482			UGA / UGA	19 / 19	
nad1	N	11527–12477 / 11502–12452	ATT / ATT	TAA / TAA		1 / 1	
trnL1	N	12479–12541 / 12454–12515			UAG / UAG	0 / 0	
rrnL	N	12542–13835 / 12516–13808				0 / 0	
trnV	N	13836–13904 / 13809–13877			UAC / UAC	-1 / -1	
rrnS	N	13904–14688 / 13877–14661				0 / 0	
A+T rich region		14689–18042 / 14662–18104					
Note: N and J indicate that the gene was located in the minor (N) and major (J) strand. The ‘/’ indicated that these from left to right were P. davidi and P. bimucronatus.

Mitochondrial genome PCGs

The total length of the 13 PCGs was 11,137 bp for both P. davidi and P. bimucronatus, and the A + T content was 68.38% and 66.74%, respectively, these two species showed a negative AT-skew (−0.146, −0.163, respectively) and a negative GC-skew (−0.049 and −0.019, respectively) in PCGs (Table 3). The longest and shortest PCGs were nad5 (1,714 bp) and atp8 (156 bp). The combined 13 PCGs of P. davidi or P. bimucronatus mitogenomes encoded a total of 3,702 codons (excluding stop codons); the most frequently used codon was UUA (Leu2), and the relative synonymous codon usage (RSCU) was 2.43 and 2.41, respectively (Fig 2). The most frequently used amino acids were Ile, Phe, Leu1 and Leu2 (Fig 3).

10.1371/journal.pone.0310559.g002 Fig 2 Relative synonymous codon usage (RSCU) in mitochondrial genomes of Propomacrus davidi and P. bimucronatus.

10.1371/journal.pone.0310559.g003 Fig 3 Codon distribution of Propomacrus davidi and P. bimucronatus.

The codon families are represented on the X-axis.

10.1371/journal.pone.0310559.t003 Table 3 Base composition and skewness of Propomacrus davidi and P. bimucronatus.

Species	Items	Size(bp)	A%	T%	G%	C%	A+T%	G+C%	AT skew	GC skew	
Propomacrus davidi	Genome	18042	39.06	30.43	9.94	20.57	69.49	30.51	0.124	-0.348	
PCGs	11137	29.19	39.19	15.03	16.58	68.38	31.61	-0.146	-0.049	
rRNA	2079	33.48	41.65	17.46	7.41	75.13	24.87	-0.109	0.404	
tRNAs	1478	37.69	36.27	14.55	11.50	73.96	26.05	0.019	0.117	
A+T region	3354	40.16	27.58	9.21	23.05	67.74	32.26	0.186	-0.429	
P. bimucronatus	Genome	18104	36.74	31.16	11.96	20.13	67.91	32.09	0.082	-0.255	
PCGs	11137	27.92	38.82	16.31	16.95	66.74	33.26	-0.163	-0.019	
rRNA	2078	33.64	40.66	17.61	8.08	74.30	25.69	-0.094	0.371	
tRNAs	1450	38.28	36.07	13.93	11.72	74.35	25.65	0.030	0.086	
A+T region	3443	35.52	29.63	13.53	21.32	65.15	34.85	0.090	-0.224	
Note: AT skew = (A − T)/ (A + T), GC skew = (G − C)/ (G + C).

The start codon of cox1 was not detected in the two species, which may be related to the use of the non-canonical start codons of this gene [61, 65]; nad6 genes were ATT and ATC, respectively; nad2, atp8, nad3, nad5, and nad1 were typical ATT; cox2 gene was ATA; atp6, cox3, nad4, nad4L, and cytb were typical ATG. The cox1, cox2, and nad5 used the incomplete stop codon T, the complete stop codon of cytb was TAG, and the complete stop codons for other genes were the typical TAA (Table 2).

The start codons of P. davidi and P. bimucronatus differ from those of other Coleoptera species. In Coleoptera, except for the cox1 gene, most protein-coding genes use a typical ATN as the start codon. However, incomplete stop codons (T or TA) are common in the mitochondrial genome of metazoans [66].

The mitochondrial tRNA and rRNA genes

The total length of 22 tRNA genes in the mitochondrial genomes of P. davidi and P. bimucronatus was 1,478 bp and 1,450 bp, respectively. Most tRNAs exhibit the typical cloverleaf structures, except trnS1 was missing the DHU arm and the anticodon was changed from GCU to UCU. Although most arthropods use a GCU anticodon in tRNA-Ser(AGN), almost all beetle mitogenomes published so far have the UCU anticodon for this tRNA [61]. From the known complete mitogenome of Coleoptera insects, Polyphage uses UCU as an anticodon, and the other suborders of Archostemata, Myxophaga, and Adiphaga are GCU [66]. This suggests that the anticodon may be a molecular synapomorphy of Coleoptera [61, 66].

In the folding process of the cloverleaf structure, there were 26 pairs of G-U mismatches in P. davidi (Fig 4), of which nine pairs were located in the amino acid acceptor arms of trnQ, trnC (3), trnY, trnA, trnF, trnH and trnV; seven pairs were located in the DHU arms of trnQ (2), trnC, trnK, trnG, trnH and trnP, and seven pairs were located in the TΨC arms of trnQ, trnC, trnS1, trnH, trnP (2), and trnV; three pairs of anticodon arms were located at trnI, trnD, and trnH. There were 24 pairs of G-U mismatches in P. bimucronatus (Fig 5), of which eight pairs were located on the amino acid acceptor arm of trnC (3), trnY, trnA, trnH, and trnV (2); seven pairs were located on the DHU arm of trnQ (2), trnC, trnK, trnG, trnH, and trnP; eight pairs were located on the TΨC arm of trnQ, trnY, trnN, trnS1, trnP (3), and trnV, and one pair was located on the anticodon arm of trnH. The unpaired base A appeared in the anticodon arm of trnS1 in the two species.

10.1371/journal.pone.0310559.g004 Fig 4 Secondary structure of tRNA in mitochondrial genes of Propomacrus davidi.

10.1371/journal.pone.0310559.g005 Fig 5 Secondary structure of tRNA in mitochondrial genes of the Propomacrus bimucronatus.

rrnL genes of the two beetles were located between the trnL1 and trnV, and rrnS genes were located between trnV and control region, as in other insects [31, 34, 61], and rrnL had greater variation than rrnS. The sequence lengths of the rrnL (Fig 6A) and rrnS (Fig 6B) of P. davidi were 1,294 bp and 785 bp, respectively, and the A + T contents were 76.43% and 72.99%, respectively. The sequence lengths of the rrnL (Fig 7A) and rrnS (Fig 7B) of P. bimucronatus were 1,293 bp and 785 bp, respectively, and the A + T contents were 76.1% and 71.34%. The secondary structures of rrnL in the two species were predicted to have 42 helical structures and five domains (I, II, IV, V, VI). The domain III deletion of rrnL is a typical feature of arthropods [31, 34, 67–69]. The secondary structures of rrnS have 26 helical structures and three domains.

10.1371/journal.pone.0310559.g006 Fig 6 a. Predicted secondary structure of the rrnL gene in Propomacrus davidi. The inferred Watson Crick bonds are represented by lines, while the GU bond is represented by dots. b. Predicted secondary structure of the rrnS gene in Propomacrus davidi. The inferred Watson Crick bonds are represented by lines, while the GU bond is represented by dots.

10.1371/journal.pone.0310559.g007 Fig 7 a. Predicted secondary structure of the rrnL gene in Propomacrus bimucronatus. The inferred Watson Crick bonds are represented by lines, while the GU bond is represented by dots. b. Predicted secondary structure of the rrnS gene in the Propomacrus bimucronatus. The inferred Watson Crick bonds are represented by lines, while the GU bond is represented by dots.

Control region

The mitogenome control region is a non-coding region of the sequence in its genome, also known as the A+T-rich region, which might be the key to characterizing typical CRs in mitogenomes [70, 71]. The control region is thought to regulate the replication and transcription of the mitogenome [72, 73]. The total lengths of the A+T-rich region of P. davidi and P. bimucronatus mitogenomes were 3,354 bp and 3,443 bp, respectively, both located between rrnS and trnI. However, the difference between these two species in CRs was not limited by these features and varied considerably. The T content in the control region is always higher than that in the coding region, which is a well-documented pattern in the insect mitogenomes [74]. The A + T contents were 67.74% for P. davidi and 65.15% for P. bimucronatus; these two species showed a positive AT-skew (0.186 and 0.090, respectively) and a negative GC-skew (−0.429 and −0.224, respectively). There were two pairs of repeat units in the 15,475–16,362 bp region in the mitogenome of P. davidi. The interval between the two pairs of repeating units was 80 bp, and the repeat sequence length for P. davidi and P. bimucronatus was 164 bp and 210 bp, respectively. The interval between the 164 bp repeat sequences was 53 bp, and the interval between the 210 bp repeat sequences was 7 bp. In the P. bimucronatus mitogenome, there was a pair of repeat units with a length of 446 bp at 15, 504–16,149 bp region, and the length of the overlapping region between the repeat sequences was 246 bp (Fig 8). The control region of the two species also has one poly-T with a length of 20 bp. There were many microsatellite structures in the control region. P. davidi had a total of 12 (TA)n structures, including one (TA)8, three (TA)5, one (TA)7, four (TA)3, two (TA) 4 and one (TA)10. P. bimucronatus has nine (TA)n structures, one (TA)8, one (TA)7, and seven (TA)3.

10.1371/journal.pone.0310559.g008 Fig 8 The control region of mitochondrial genome of Propomacrus davidi and P. bimucronatus.

Non-repetitive unit (red), overlapping regions (green), intergenic spacer (purple), repeat unit (grey).

Phylogenetic analysis

Based on the 13 PCG datasets of 50 species, ML and BI trees with almost the identical topology were generated (Fig 9). The difference between the two trees is only the position of Mecynorhina polyphemus and Dicronorhina derbyana, and the two species have only a slight exchange of positions. The phylogenetic relationships at the subfamily level of Scarabaeidae constructed by the study are: ((((Cetoniinae + (Dynastinae + Rutelinae)) + Melolonthinae) + Euchirinae) + Scarabaeinae. Among the 25 species of Scarabaeidae, the four species of Euchirinae (P. davidi, P. bimucronatus, C. jansoni, C. gestroi) clustered into the same branch in the phylogenetic tree and were the most closely related, showing Euchirinae monophyly and a nested position within the Scarabaeidae, supporting the morphological characteristics conclusion of Šípek et al. (2011). At the same time, the results continue to support the phylogenetic relationship among the subfamilies Cetoniinae, Dynastinae, Rutelinae, Melolonthinae, and Euchirinae [6, 16, 29, 31, 75, 76]. In addition, the unclear phylogenetic relationships at the genus level also indicate to us that the data are insufficient. We need more basic research to support the species boundaries and evolutionary relationships of this charismatic-looking beetle.

10.1371/journal.pone.0310559.g009 Fig 9 Phylogenetic tree of Scarabaeoidea produced from maximum likelihood (ML) (A) and Bayesian inference (BI) (B) based on 13 PCGs. Bayesian posterior probabilities are indicated on branches. The two species sequenced in this study are highlighted in red.

The superfamily Scarabaeoidea constructed by ML (Fig 9A) and BI (Fig 9B) methods are: ((Scarabaeidae + Hybosoridae) + Lucanidae) + Passalidae. Scarabaeidae and Hybosoridae are sister groups, and the formation of a clade between the two families is supported strongly. Although it is different from the phylogenetic relationship of the families of Scarabaeoidea based on studies from nuclear protein-coding genes [77, 78], which indicated that the Lucanidae did not form a sister-group relationship with the Scarabaeidae and Hybosoridae, the results support the monophyly of the four families and are consistent with results of morphological characteristics studies [17, 18].

Conclusion

In this study, the complete mitogenome sequences of P. davidi and P. bimucronatus were sequenced and analyzed; the total lengths of the two species were 18,042 bp and 18,104 bp, respectively, which are longer than those of most metazoans, and the genomes showed an obvious AT-skew. Each gene composition and sequence were consistent with those of ancestral insects, including 13 PCGs, 22 tRNAs, two rRNAs, a control region, and no gene rearrangements or deletions. The mitogenome of the two species uses UCU as an anticodon. The secondary structure models of rrnL and rrnS were similar to those of other insect genes, and rrnL had greater variation than rrnS. The domain III deletion of rrnL is a typical feature of arthropods. The T content in the control region is always higher than in the coding region, a well-documented pattern in the insect mitogenome. The total lengths of the A+T-rich region of P. davidi and P. bimucronatus mitogenomes were 3,354 bp and 3,443 bp, respectively, both located between rrnS and trnI.

The analysis based on the ML and BI method of 13 PCGs produced a well-resolved framework supporting the conclusion that Euchirinae is a monophyletic subfamily of Scarabaeidae. Mitogenome analysis is useful for comparing genetic relationships at classificatory levels. More research is needed to reconstruct comprehensive phylogenies to better resolve the phylogenetic relationships of Euchirinae.
==== Refs
References

1 Sfenthourakis S , Hadjiconstantis M , Makris C , Dimitriou AD . Revisiting the saproxylic beetle ‘Propomacrus cypriacus Alexis & Makris, 2002’ (Coleoptera: Euchiridae) using molecular, morphological and ecological data. J Nat Hist. 2017; 51 (17–18 ): 1021–1034. doi: 10.1080/00222933.2017.1319521
2 Lumaret JP , Tauzin P . Le genre Propomacrus Newman, 1837. Donnees biologiques et morphologie larvaire (Coleoptera, Scarabaeoidea, Euchiridae). Nouv Revue Ent (N.S.). 1992; 9 (2 ): 173–180.
3 Šípek P , Janšta P , Král D . Immature stages of Euchirinae (Coleoptera: Scarabaeoidea): genera Cheirotonus and Propomacrus with comments on their phylogeny based on larval and adult characters. Invertebr Syst. 2011; 25 (4 ): 282–302. doi: 10.1071/IS11028
4 Bai M , Li WZ , Yang XK . Near extinct species-Euchirinae. Sci World. 2007; (6 ): 54. (In Chinese)
5 Zhang WW . Wild China. Euchirinae: Insect ’gibbon’. Sci World. 2013; (6 ): 42–45. (In Chinese)
6 Yang WZ , Zhang Y , Feng SQ , Liu LJ , Li ZH . The first complete mitochondrial genome of the Japanese beetle Popillia japonica (Coleoptera: Scarabaeidae) and its phylogenetic implications for the superfamily Scarabaeoidea. Int J Biol Macromol. 2018; 118 (Part B ): 1406–1413. doi: 10.1016/j.ijbiomac.2018.06.131 29959994
7 Wu L , Wu ZQ . Propomacrus davidi fujianensis ssp. Nov., A new Euchirinae from Fujian (Coleoptera, Euchiridae). Acta Zootaxonomica Sinica. 2008; 33 (4 ): 827–828. (In Chinese)
8 Zhang F , Hon B , Wang YZ , Li YM , Chen ZJ . Sequencing and phylogenetic analysis of the complete mitochondrial genome of Scythropus yasumatsui (Coleoptera: Curculionidae). Acta Ecol Sinica. 2019; 62 (11 ): 1305–1314. doi: 10.16380/j.kcxb.2019.11.008 (In Chinese)
9 Zhang JD , Li TQ , Zhou FW , Wang JR , Yang JR , Xu ZH . Preliminary Study on the Morphological Characters and Breeding of Propomacrus davidi Deyrolle. J Yunnan Agr Univ (Nat Sci). 2019; 34 (1 ): 175–179. (In Chinese)
10 Young RM . Euchirinae (Coleoptera: Scarabaeidae) of the world: distribution and taxonomy. Coleopterists Bull. 1989; 43 (3 ): 205–236.
11 Alexander K , Nieto A . European Red List of saproxylic beetles. Luxembourg: Publications Office of the European Union. 2010. http://ec.europa.eu/environment/nature/conservation/species
12 Yi CH , Chen Y , He QJ , Zhou Y , Wang L , Chen XM . Morphometric characteristics Study on Cherirotonus gestroi Pouillaud. J Northwest For Univ. 2015; 30 (2 ): 154–157. (In Chinese)
13 Liu SL , Liu LJ , Wang D , Wu FS , Zhou SZ . A female description of Cheirotonus jansoni (Coleoptera, Euchiridae) from China. Hubei Agr Sci. 2022; 61 (1 ): 89–90. doi: 10.14088/j.cnki.issn0439-8114.2022.01.016 (In Chinese)
14 Joharchi O , Halliday B , Beyzavi G . A new species of the genus Promacrolaelaps (Acari: Laelapidae) associated with Propomacrus bimucronatus (Pallas) (Coleoptera: Scarabaeidae) in Iran. Zootaxa. 2013; 3641 (4 ): 379–383. doi: 10.11646/zootaxa.3641.4.4 26287092
15 Gunter NL , Weir TA , Slipinksi A , Bocak L , Cameron SL . If Dung Beetles (Scarabaeidae: Scarabaeinae) Arose in Association with Dinosaurs, Did They Also Suffer a Mass Co-Extinction at the K-Pg Boundary? PLoS ONE. 2016; 11 (5 ): e0153570. doi: 10.1371/journal.pone.0153570 27145126
16 Guo SB , Lin XY , Song N . Mitochondrial phylogenomics reveals deep relationships of scarab beetles (Coleoptera, Scarabaeidae). PLoS ONE. 2022; 17 (12 ): e0278820. doi: 10.1371/journal.pone.0278820 36512580
17 Grebennikov VV , Scholtz CH . The basal phylogeny of Scarabaeoidea (Insecta: Coleoptera) inferred from larval morphology. Invertebr Syst. 2004; 18 (3 ): 321–348. doi: 10.1071/IS03013
18 Smith AB , Hawks DC , Heraty JM . An Overview of the Classification and Evolution of the Major Scarab Beetle Clades (Coleoptera: Scarabaeoidea) Based on Preliminary Molecular Analyses. Coleopterists Bull. 2006; 60 (5 ): 35–46. doi: 10.1649/0010-065X(2006)60[35:AOOTCA]2.0.CO;2
19 Browne J , Scholtz CH . A phylogeny of the families of Scarabaeoidea (Coleoptera). Syst. Entomol. 1999; 24 (1 ): 51–84. doi: 10.1046/j.1365-3113.1999.00067.x
20 Yang C , Zhu EJ , He QJ , Yi CH , Hu SJ , Wang XB . Complete mitochondrial genome of the Endangered long-armed scarab Cheirotonus gestroi (Coleoptera: Euchiridae). Mitochondrial DNA B. 2020; 5 (1 ): 869–870. doi: 10.1080/23802359.2020.1715888 33366789
21 Scholz CH , Grebennikov VV . Scarabaeiformia Crowson, 1960. In Handbook of Zoology. A Natural History of the Phyla of the Animal Kingdom. Vol. 4 Arthropoda: Insecta, Part 38. Coleoptera, Beetles Vol. 1: Morphology and Systematics (Archostemata, Adephaga, Myxophaga, Polyhaga partim; Beutel R.G and Leschen R.A.B ., eds; Walter. de Gruyter. 2005; pp: 345–425. doi: 10.1515/9783110904550.345
22 Bezděk A . Eucheirinae In Catalogue of Palaearctic Coleoptera. Vol. 3 . Scarabaeoidea—Scirtoidea—Dascilloidea—Buprestoidea—Byrrhoidea. Stenstrup: Apollo Books; Löbl I . and Smetana A . eds. 2006; pp. 181.
23 Ahrens D. The phylogeny of Sericini and their position within the Scarabaeidae based on morphological characters (Coleoptera: Scarabaeidae). Syst Entomol. 2006; 31 (1 ): 113–144. doi: 10.1111/j.1365-3113.2005.00307.x
24 Xing K , Kang C , Zhao XJ , Zhao F . Sequencing the complete mitochondrial genome of Nysius ericae(Schilling) (Hemiptera:Lygaeidae) and phylogenetic analysis of Lygaeoidea. J Environ Entomol. 2023; 45 (1 ): 142–154. (In Chinese)
25 Simon C , Frati F , Beckenbach A , Crespi B , Liu H , Flook P . Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994; 87 (6 ): 651–701. doi: 10.1093/aesa/87.6.651
26 Hebert PD , Ratnasingham S , deWaard JR . Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc Biol Sci. 2003; 270 (1 ): 96–99. doi: 10.1098/rsbl.2003.0025 12952648
27 Boore JL . Animal mitochondrial genomes. Nucleic Acids Res. 1999; 27 (8 ): 1767–1780. doi: 10.1093/nar/27.8.1767 10101183
28 Kim KG , Hong MY , Kim MJ , Im HH , Kim MI , Bae CH , et al . Complete mitochondrial genome sequence of the yellow-spotted long-horned beetle Psacothea hilaris (Coleoptera: Cerambycidae) and phylogenetic analysis among coleopteran insects. Mol Cells. 2009; 27 (4 ): 429–441. doi: 10.1007/s10059-009-0064-5 19390824
29 Yu YY , Lee SH , Jin LJ , Bai M . Mitochondrial genome of Euchirus longimanus (Coleoptera: Scarabaeidae) and phylogenetic relationship of Euchirini. J Asia-Pac Entomol. 2023; 26 (4 ): 102156. doi: 10.1016/j.aspen.2023.102156
30 Cameron SL . Insect mitochondrial genomics: implications for evolution and phylogeny. Annu Rev Entomol. 2014; 59 : 95–117. doi: 10.1146/annurev-ento-011613-162007 24160435
31 Ayivi SPG , Tong Y , Storey KB , Yu DN , Zhang JY . The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae. Insects. 2021; 12 (11 ): 1025. doi: 10.3390/insects12111025 34821825
32 Jeong JS , Kim MJ , Kim I . The mitochondrial genome of the dung beetle, Copris tripartitus, with mitogenomic comparisons within Scarabaeidae (Coleoptera). Int J Biol Macromol. 2020; 144 (1 ): 874–891. doi: 10.1016/j.ijbiomac.2019.09.165 31726138
33 Bankevich A , Nurk S , Antipov D , Gurevich AA , Dvorkin M , Kulikov AS , et al . SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012; 19 (5 ): 455–477. doi: 10.1089/cmb.2012.0021 22506599
34 Shao LL , Huang DY , Sun XY , Hao JS , Cheng CH , Zhang W , et al . Complete mitochondrial genome sequence of Cheirotonus jansoni (Coleoptera: Scarabaeidae). Genet Mol Res. 2014; 13 (1 ): 1047–1058. doi: 10.4238/2014.February.20.6 24634126
35 Kurtz S , Phillippy A , Delcher AL , Smoot M , Shumway M , Antonescu C , et al . Versatile and open software for comparing large genomes. Genome Biol. 2004; 5 : R12. doi: 10.1186/gb-2004-5-2-r12 14759262
36 Walker BJ , Abeel T , Shea T , Priest M , Abouelliel A , Sakthikumar S , et al . Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE. 2014; 9 (11 ): e112963. doi: 10.1371/journal.pone.0112963 25409509
37 Bernt M , Donath A , Jühling F , Externbrink F , Florentz C , Fritzsch G , et al . MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013; 69 (2 ): 313–319. doi: 10.1016/j.ympev.2012.08.023 22982435
38 Lowe TM , Eddy SR . TRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997; 25 (5 ): 955–964. doi: 10.1093/nar/25.5.955 9023104
39 Stephan G , Pascal L , Ralph Bock . OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 2019; 47 (W1 ): 59–64. doi: 10.1093/nar/gkz238 30949694
40 Zhang D , Gao FL , Jakovlić I , Zou H , Zhang J , Li WX , et al . PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol Ecol Resour. 2020a; 20 (1 ): 348–355. doi: 10.1111/1755-0998.13096 31599058
41 Nguyen LT , Schmidt HA , von Haeseler A , Minh BQ . IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015; 32 (1 ): 268–274. doi: 10.1093/molbev/msu300 25371430
42 Ronquist F , Huelsenbeck JP . MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003; 19 (12 ): 1572–1574. doi: 10.1093/bioinformatics/btg180 12912839
43 Kim MJ , Im HH , Lee KY , Han YS , Kim I . Complete mitochondrial genome of the whiter-spotted flower chafer, Protaetia brevitarsis (Coleoptera: Scarabaeidae). Mitochondr DNA. 2014; 25 (3 ): 177–178. doi: 10.3109/19401736.2013.792064 23841598
44 Kim MJ , Jeong SY , Jeong J C , Kim SS , Kim I . Complete mitochondrial genome of the endangered flower chafer Osmoderma opicum (Coleoptera: Scarabaeidae). Mitochondrial DNA B. 2016; 1 (1 ): 148–149. doi: 10.1080/23802359.2016.1144104 33473440
45 Li YJ , Nie RE , Lu YY , Lee S , Zhao ZY , Wu L , et al . Seven new mitochondrial genomes of phytophagous scarab beetles (Coleoptera: Scarabaeidae) and phylogenetic implications. Zootaxa. 2022; 5138 (3 ): 324–338. doi: 10.11646/zootaxa.5138.3.6 36095833
46 Cameron SL , Sullivan J , Song HJ , Miller KB , Whiting MF . A mitochondrial genome phylogeny of the Neuropterida (lace-wings, alderflies and snakeflies) and their relationship to the other holometabolous insect orders. Zool Scr. 2009; 38 (6 ): 575–590. doi: 10.1111/j.1463-6409.2009.00392.x
47 Amorim IC , Melo AS , Cruz GAS , Wallau GL , Moura RC . Dichotomius (Luederwaldtinia) schiffleri (Coleoptera: Scarabaeidae) mitochondrial genome and phylogenetic relationships within the superfamily Scarabaeoidea. Mitochondrial DNA B. 2017; 2 (2 ), 887–888. doi: 10.1080/23802359.2017.1407695 33474024
48 Lee SG , Choi KS , Kim CJ , Jeon HB , Kim IK . Complete mitochondrial genome of Leptaulax koreanus (Coleoptera: Passalidae), a Korean endemic bess beetle. Mitochondrial DNA B. 2019; 4 (1 ): 105–106. doi: 10.1080/23802359.2018.1536492
49 Chen D , Liu J , Bartolozzi L , Wan X . The complete mitochondrial genome of stag beetle Lucanus cervus (Coleoptera: Lucanidae) and phylogenetic analysis. PeerJ. 2019; 7 : e8274. doi: 10.7717/peerj.8274 31875157
50 Wang LJ , Huang JP , Shiao SF , Ko HP , Sung CH . Characterisation of the complete mitochondrial genome of Lucanus chengyuani (Coleoptera: Lucanidae). Mitochondrial DNA B, 2019; 4 (2 ): 3460–3462. doi: 10.1080/23802359.2019.1674703 33366039
51 Liu J , Cao YY , Zhou SJ , Chen YJ , Wan X . Complete mitochondrial genome of Prismognathus prossi (Coleoptera: Lucanidae) with phylogenetic implications. Entomol Fennica. 2019; 30 (2 ): 90–96. doi: 10.33338/ef.82927
52 Chen YJ , Liu J , Cao YY , Zhou SJ , Wan X . Two new complete mitochondrial genomes of Dorcus stag beetles (Coleoptera, Lucanidae). Genes genom, 2018; 40 (8 ): 873–880. doi: 10.1007/s13258-018-0699-8 30047116
53 Zhai YN , Zhou SJ , Chen YJ , Wan X . The mitochondrial genome of a rare Chinese stag beetle Kirchnerius guangxii (Coleoptera: Lucanidae). Mitochondrial DNA B. 2020; 5 (2 ): 1633–1635. doi: 10.1080/23802359.2020.1745708
54 Xu MQ , Zhou SJ , Wan X . Phylogenetic Implication of Large Intergenic Spacers: Insights from a Mitogenomic Comparison of Prosopocoilus Stag Beetles (Coleoptera: Lucanidae). Animals. 2022; 12 (13 ): 1595. doi: 10.3390/ani12131595 35804495
55 Wu YY , Cao YY , Fang J , Wan X . The first complete mitochondrial genome of stag beetle from China, Prosopocoilus gracilis (Coleoptera, Lucanidae). Mitochondrial DNA A. 2016; 27 (4 ): 2633–2634. doi: 10.3109/19401736.2015.1041129 26024142
56 Liu J , Zhou SJ , Chen YJ , Wan X . Mitogenome of the Monotypic Genus Rhaetus (Coleoptera: Scarabaeidae: Lucanidae). J Entomol Sci. 2018; 53 (4 ): 503–513. doi: 10.18474/JES17-122.1
57 Lee JM , Park JH , Xi H , Park JS . Comprehensive analyses of the complete mitochondrial genome of Figulus binodulus (Coleoptera: Lucanidae). J Insect Sci. 2020; 20 (5 ): 1–13. doi: 10.1093/jisesa/ieaa090 32976575
58 Wang Q , Liu J , Lin ZQ , Wan X . The complete mitochondrial genome of Odontolabis fallaciosa (Coleoptera: Lucanidae) with its phylogenetic implications. Zool Syst. 2018; 43 (3 ): 268–275. doi: 10.11865/zs.201831
59 Wolstenholme DR . Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992 141 : 173–216. doi: 10.1016/s0074-7696(08)62066-5 1452431
60 Wei ZH . The complete mitochondrial genomes of five Agrilinae (Coleoptera, Buprestidae) species and phylogenetic implications. Zookeys. 2022; 1092 : 195–212. doi: 10.3897/zookeys.1092.80993 35586288
61 Sheffield NC , Song H , Cameron SL , Whiting MF . A Comparative Analysis of Mitochondrial Genomes in Coleoptera (Arthropoda: Insecta) and Genome Descriptions of Six New Beetles. Mol Biol Evol. 2008; 25 : 2499–2509. doi: 10.1093/molbev/msn198 18779259
62 Wu C , Zhou Y , Tian T , Li TJ , Chen B . First report of complete mitochondrial genome in the subfamily Alleculinae and mitochondrial genome-based phylogenetics in Tenebrionidae (Coleoptera: Tenebrionoidea). Insect Sci. 2022; 29 (4 ): 1226–1238. doi: 10.1111/1744-7917.12983 34791791
63 Huang XY , Chen B , Wei ZH , Shi AM . First Report of Complete Mitochondrial Genome in the Tribes Coomaniellini and Dicercini (Coleoptera: Buprestidae) and Phylogenetic Implications. Genes. 2022; 13 (6 ): 1074. doi: 10.3390/genes13061074 35741836
64 Herms DA , McCullough DG . Emerald ash borer invasion of North America: history, biology, ecology, impacts, and management. Annu Rev Entomol. 2014; 59 : 13–30. doi: 10.1146/annurev-ento-011613-162051 24112110
65 Cameron SL , Whiting MF . The complete mitochondrial genome of the tobacco hornworm, Manduca sexta, (Insecta: Lepidoptera: Sphingidae), and an examination of mitochondrial gene variability within butterflies and moths. Gene. 2008; 408 (1–2 ): 112–123. doi: 10.1016/j.gene.2007.10.023 18065166
66 Nie RE , Yang XK . Research progress in mitochondrial genomes of Coleoptera. Acta Entomol Sin. 2014; 57 (7 ): 860–868. doi: 10.16380/j.kcxb.2014.07.015 (In Chinese)
67 Gillespie JJ , Johnston JS , Cannone JJ , Gutell RR . Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements. Insect Mol Biol. 2006; 15 (5 ): 657–686. doi: 10.1111/j.1365-2583.2006.00689.x 17069639
68 Chai HN , Du YZ . The complete mitochondrial genome of the pink stem borer, Sesamia inferens, in comparison with four other Noctuid moths. Int J Mol Sci. 2012; 13 (8 ): 10236–10256. doi: 10.3390/ijms130810236 22949858
69 Zhang RY , Li J , Geng S , Yang J , Zhang X , An YX , et al . The first mitochondrial genome for Phaudidae (Lepidoptera) with phylogenetic analyses of Zygaenoidea. Int J Biol Macromol. 2020; 149 : 951–961. doi: 10.1016/j.ijbiomac.2020.01.307 32018010
70 Pang YT , Feng SQ , Wu Y , Zhang SF , Li ZH . Sequencing and analysis of the mitochondrial genome of Trogoderma variabile (Coleoptera: Dermestidae). Plant Quarantine. 2019; 33 (1 ): 8–13. doi: 10.19662/j.cnki.issn1005-2755.2019.01.002 (In Chinese)
71 Zeng LY , Pang YT , Feng SQ , Wang YN , Stejskal V , Aulicky R , et al . Comparative mitochondrial genomics of five Dermestid beetles (Coleoptera: Dermestidae) and its implications for phylogeny. Genomics. 2021; 113 (1 Pt 2 ): 927–934. doi: 10.1016/j.ygeno.2020.10.026 33131762
72 Zhang DX , Szymura JM , Hewitt GM . Evolution and structural conservation of the control region of insect mitochondrial DNA. J Mol Evol. 1995; 40 : 382–391. doi: 10.1007/BF00164024 7769615
73 Zhang DX , Hewitt GM . Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997; 25 (2 ): 99–120. doi: 10.1016/S0305-1978(96)00042-7
74 Clary DO , Wolstenholme DR . The mitochondrial DNA molecular of Drosophila yakuba: nucleotide sequence, gene organization, and genetic code. J Mol Evol. 1985; 22 : 252–271. doi: 10.1007/BF02099755 3001325
75 Dietz L , Seidel M , Eberle J , Misof B , Pacheco TL , Podsiadlowski L , et al . A transcriptome-based phylogeny of Scarabaeoidea confirms the sister group relationship of dung beetles and phytophagous pleurostict scarabs (Coleoptera). Syst Entomol. 2023; 48 (4 ): 672–686. doi: 10.1111/syen.12602
76 Šípek P , Fabrizi S , Eberle J , Ahrens D . A molecular phylogeny of rose chafers (Coleoptera: Scarabaeidae: Cetoniinae) reveals a complex and concerted morphological evolution related to their flight mode. Mol Phylogenet Evol. 2016; 101 : 163–175. doi: 10.1016/j.ympev.2016.05.012 27165937
77 Zhang SQ , Che LH , Li Y , Liang D , Pang H , Ślipiński A , et al . Evolutionary history of Coleoptera revealed by extensive sampling of genes and species. Nat Commun. 2018; 9 (1 ), 205. doi: 10.1038/s41467-017-02644-4 29335414
78 Mckenna DD , Wild AL , Kanda K , Bellamy CL , Beutel RG , Caterino MS , et al . The Beetle Tree of Life Reveals that Coleoptera Survived End Permian Mass Extinction to Diversify During the Cretaceous Terrestrial Revolution. Syst Entomol. 2015; 40 (4 ): 835–880. doi: 10.1111/syen.12132
