The complete mitochondrial genome of Hemathlophorus brevigenatus Wei, 2005 (Hymenoptera: Tenthredinidae) with phylogenetic analysis

Abstract The mitochondrial genome of Hemathlophorus brevigenatus Wei, 2005 collected from Huanggang Mountain of China is described using the NGS approach. The length of the sequence is 15,452 bp containing 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region. The overall A + T content is 79.5%. tRNA rearrangements occur in the MQI cluster. Phylogenetic analysis of H. brevigenatus resolved it in a clade with Allantus togatus in Allantinae which provides new evidence for the phylogeny of Tenthredinidae.

erected the tribe Athlophorini for Hemathlophorus Malaise, 1945 and Athlophorus Burmeister, 1847 under the subfamily Sioblinae of Tenthredinidae. Although Hemathlophorus and Athlophorus are two peculiar genera among the Tenthredinidae as shown by the fore wing with a long vein R þ M and a more or less basally constricted abdomen, Lacourt's peculiar systematic arrangement of the genera and the tribe was not adopted by sawfly researchers, such as Wei and Nie (1998), Saini (2006), andTaeger et al. (2010). In the three later systems, Hemathlophorus and Athlophorus were placed into the Allantinae of Tenthredinidae. In this study, we sequenced the mitochondrial genome of Hemathlophorus brevigenatus Wei, 2005 and inferred its phylogenetic history with other mitochondrial genomes of Symphytan species to clarify the phylogenetic position of Hemathlophorus within the Tenthredinidae.
Specimens of H. brevigenatus were collected from Huanggangshan Mountain, Fujian Province, China (27.81 N,117.95 E) in April 2019. The specimen was deposited at the Asia Sawfly Museum, Nanchang (ASMN) (Meicai Wei, weimc@ 126.com) under the voucher number CSCS-Hym-MC0177. Whole genomic DNA was extracted from the thorax muscle of a female adult using the DNeasy Blood &Tissue Kit (Qiagen, Valencia, CA). Genomic DNA was analyzed with the high-throughput Illumina Hiseq 4000 platform with 150 bp paired-end reads. DNA sequences were assembled using MitoZ (Meng et al. 2019) and Geneious Prime 2019.2.1 (Biomatters Ltd., Auckland, New Zealand). A total of 43,205,384 raw reads were assembled using MitoZ and resulted in a contig 15,108 bp in length with 37 genes. The control region was assembled using trnM, trnI from the Xenapatidea procincta (MW487928) as a reference. The contig was extended 3 0 by 117 bp and 5 0 by 582 bp using the "Map to Reference" function in Geneious Prime. This extension that resulted overlapped with the 15,108 bp contig and closed the gap. Thus, a 482 bp long control region was obtained. The verification is conducted by assembly using X. procincta as the reference. The annotations of tRNAs were generated using the MITOS web server (Bernt et al. 2013). Protein-coding genes (PCGs) were annotated by the open reading frames between the flanking tRNAs, and then defined by comparative analyses (Cameron 2014).
The complete mitochondrial genome sequences were aligned with ClustalW using default settings (Thompson et al. 1994) and concatenated with SequenceMatrix v1.7.8 (Vaidya et al. 2011). The phylogenetic tree was constructed using IQ-TREE (maximum-likelihood) (Nguyen et al. 2015) with the GTR þ MTART model and 1000 bootstrap replicates. To investigate the phylogenetic position of H. brevigenatus, 10 unsaturated nucleotide sequences (atp8, nad2, nad6 are excluded) of 37 Hymenopterans are aligned separately then concatenated, resulting in an alignment of 9836 bp.

Disclosure statement
No potential conflict of interest was reported by the authors. The authors alone are responsible for the content and writing of the article.

Data availability statement
The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at [https://www.ncbi.nlm.nih.gov] (https://www.ncbi.nlm.nih.gov/) under the accession number MW632125. The associated BioProject, SRA, and BioSample numbers are PRJNA714776, SRR14233978, and SAMN18397095 respectively.