Research Article |
Corresponding author: Zhou Yafu ( zyf820207@126.com ) Corresponding author: Shi Fuming ( shif_m@126.com ) Academic editor: Fernando Montealegre-Z
© 2018 Mao Shaoli, Yuan Hao, Lu Chao, Zhou Yafu, Shi Fuming, Wang Yuchao.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Shaoli M, Hao Y, Chao L, Yafu Z, Fuming S, Yuchao W (2018) The complete mitochondrial genome of Xizicus (Haploxizicus) maculatus revealed by next-generation sequencing and phylogenetic implication (Orthoptera, Meconematinae). ZooKeys 773: 57-67. https://doi.org/10.3897/zookeys.773.24156
|
Xizicus Gorochov, 1993, the quiet-calling katydid, is a diverse genus with 68 species in world, which includes more than 45 species in China, has undergone numerous taxonomic revisions with contradicting conclusions. In this study the complete mitochondrial genome of Xizicus (Haploxizicus) maculatus collected from Hainan for the first time was sequenced using the Next-Generation Sequencing (NGS) technology. The length of whole mitogenome is 16,358 bp and contains the typical gene arrangement, base composition, and codon usage found in other related species. The overall base composition of the mitochondrial genome is 37.0 % A, 32.2 % T, 20.2 % C, and 10.6 % G. All 13 protein-coding genes (PCGs) began with typical ATN initiation codon. Nine of the 13 PCGs have a complete termination codon, but the remaining four genes (COI, COIII, ND5, and ND4) terminate with an incomplete T. Phylogenetic analyses are carried out based on the concatenated dataset of 13 PCGs and two rRNAs of Tettigoniidae species available in GenBank. Both Bayesian inference and Maximum Likelihood analyses recovered each subfamily as a monophyletic group. Regardless of the position of Lipotactinae, the relationships among the subfamilies of Tettigoniidae were as follows: ((((Tettigoniinae, Bradyporinae) Meconematinae) Conocephalinae) Hexacentrinae). The topological structure of the phylogeny trees showed that the Xizicus (Haploxizicus) maculatus is closer to Xizicus (Xizicus) fascipes than Xizicus (Eoxizicus) howardi.
mitochondrial genome, Next-Generation Sequencing, phylogenetic relationship, Xizicus (Haploxizicus) maculatus
Insect mitochondrial genome (mitogenome) occurs as a small (15–20kb), circular, and double-stranded DNA molecules, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes (
Xizicus Gorochov, 1993 is a diverse genus in Meconematinae with 68 species in the world, currently divided into seven subgenera by traditional taxonomy based on comparative morphology: Xizicus s. str., Axizicus, Eoxizicus, Furcixizicus, Haploxizicus, Paraxizicus, Zangxizicus according to the OSF website (
To date, the mitogenomes of Xizicus (Xizicus) fascipes and Xizicus (Eoxizicus) howardi have been sequenced (
The specimen of X. (H.) maculatus was collected at Jianfengling in Hainan, China in June 2017 and stored in 100 % ethanol at −4 °C. Genomic DNA was extracted from the leg muscle tissue of a single adult male species of using a DNeasy Blood & Tissue Kit (Qiagen, USA) according to the manufacturer’s instruction and sent to a company for library prep and sequencing (Genesky Biotechnologies Inc., Shanghai). The library was prepared using a TruSeq DNA sample Preparation kit (Vanzyme, China) and sequenced with 150 bp pair-end reads on the Illumina Hiseq 2500 sequencing platform (Illumina, USA).
13,799,778 raw reads were sequenced by the Illumina Hiseq 2500 platform. The raw paired-end reads were filtered to obtain high-quality clean reads by using CLC Genomics Workbench 8 (CLC Bio, Aarhus, Denmark) with default parameters. Then the filtered reads were aligned to the mitochondrial genome of X. (X.) fascipes (JQ326212) as a reference using MITObim v1.8 (Hahn et al. 2013) and Mira 4.0.2 (
Phylogenetic analyses were performed on the concatenated datasets of PCGs and rRNAs of the newly sequenced mitogenome and 24 Tettigoniidae species downloaded from GenBank, with two Phaneropteridae taxa (Ducetia japonica and Phyllomimus sinicus) selected as the outgroups. Alignment of each protein-coding gene inferred from the amino acid alignment was performed using MEGA v7.0 (
The complete mitogenome of X. (H.) maculatus is 16,358 bp in length and has been deposited in GenBank under accession no. MG779499. Xizicus (H.) maculatus mtDNA is larger than that observed in other species in Meconematinae, which typically ranged from 16,044 bp (
The mitochondrial genome structure is detailed in Table
Gene/region | Position | Size | Direction | Initiation codon | Termination codon | anticodon |
---|---|---|---|---|---|---|
tRNA Ile | 1-66(-3) | 66 | F | GAT | ||
tRNA Gln | 64-132(+8) | 69 | R | TTG | ||
tRNA Met | 141-204 | 64 | F | CTA | ||
ND2 | 205-1233(-2) | 1029 | F | ATT | TAA | |
tRNATrp | 1232-1297(-8) | 66 | F | TCA | ||
tRNA Cys | 1290-1353 | 64 | R | GCA | ||
tRNA Tyr | 1354-1419(-8) | 66 | R | GTA | ||
COI | 1412-2951 | 1540 | F | ATT | T | |
tRNA Leu(UUR) | 2452-3016(+3) | 65 | F | TAA | ||
COII | 3020-3703(+1) | 684 | F | ATG | TAA | |
tRNA Lys | 3705-3774(-1) | 70 | F | CTT | ||
tRNA Asp | 3774-3839 | 66 | F | GTC | ||
ATP8 | 3840-4004(-7) | 165 | F | ATT | TAA | |
ATP6 | 3998-4675(-1) | 678 | F | ATG | TAA | |
COIII | 4675-5461 | 787 | F | ATG | T | |
tRNAGly | 5462-5526 | 65 | F | TCC | ||
ND3 | 5527-5880(+2) | 354 | F | ATC | TAA | |
tRNA Ala | 5883-5946(-1) | 64 | F | TGC | ||
tRNAArg | 5946-6008(+16) | 63 | F | TCG | ||
tRNA Asn | 6025-6090(+2) | 66 | F | GTT | ||
tRNA Ser(AGN) | 6093-6159 | 67 | F | GCT | ||
tRNAGlu | 6160-6226(+14) | 67 | F | TTC | ||
tRNA Phe | 6305-6421 | 65 | R | GAA | ||
ND5 | 6306-8037 | 1732 | R | ATT | T | |
tRNA His | 8038-8101 | 64 | R | GTG | ||
ND4 | 8102-9440(-7) | 1339 | R | ATG | T | |
ND4L | 9434-9730(+1) | 297 | R | ATG | TAA | |
tRNA Thr | 9732-9802(-1) | 71 | F | TGT | ||
tRNA Pro | 9802-9867(+1) | 66 | R | TGG | ||
ND6 | 9869-10396(-1) | 528 | F | ATA | TAA | |
Cytb | 10396-11532(-2) | 1137 | F | ATG | TAG | |
tRNA Ser(UCN) | 11531-11599(+17) | 69 | F | TGA | ||
ND1 | 11617-12570(-6) | 954 | R | ATA | TAG | |
tRNA Leu(CUN) | 12565-12628 | 64 | R | TAG | ||
lrRNA | 12629-13932 | 1304 | R | |||
tRNA Val | 13933-14003 | 71 | R | TAC | ||
srRNA | 14004-14788 | 785 | R | |||
Control region | 14789-16358 | 1570 | F |
The total length of all 13 PCGs was 11,224 bp, and the overall A+T content of X. (H.) maculatusPCGs was 68.3 %. The initiation codons of all PCGs were typical with ATN (COII, ATP6, COIII, ND4, ND4L, and Cytb with ATG, ND2, COI, ATP8, ND5 with ATT; ND6, ND1 with ATA, Only ND3 with ATC). The ATN codon was prevalently regarded as initiation codons for mitogenome PCGs in insects. Seven genes (ND2, COII, ATP8, ATP6, ND3, ND4L, ND6) used TAA as the termination codons, and two genes (Cytb, ND1) were stopped with TAG. COI, COII, ND5, and ND4 had an incomplete termination codon T (Table
The four most-used amino acids in X. (H.) maculatus were Leu (16.1 %), Ser (8.9 %), Phe (8.7 %), and Ile (8.4 %), whose proportions were similar to those observed in other Tettigoniidae species. All codons were present in the protein-coding genes of this mitogenome. Excluding incomplete termination codons, there were 3,735 codons in the X. (H.) maculatus protein-coding genes. The codon usage in X. (H.) maculatus appeared to be typical of other insect mitochondrial sequences. The RSCU analysis indicated that codons including A or T at the third position were always overused compared with other synonymous codons in Xizicus (Fig.
The length of tRNA genes ranged from 63 to 71 bp and the relative locations for each tRNA are shown in Table
The lrRNA and srRNA were 1304bp and 785bp in length, respectively. They were located between tRNALeu(CUN) and A+T-rich region, being separated by tRNAVal.
The control region was 1570 bp in length and located between srRNA and tRNAIle in X. (H.) maculatus mitogenome, and was composed of 64.4 % A and T nucleotides (Fig.
Bayesian analyses and maximum likelihood produced identical topologies using the best-fit partitioning scheme and site-homogeneous models, excepting the location of Lipotactes tripyrga (Lipotactinae) and genera relationships in Meconematinae (Fig.
In present study, regardless of the position of Lipotactinae, the relationships among the subfamilies of Tettigoniidae were as follows: ((((Tettigoniinae, Bradyporinae) Meconematinae) Conocephalinae) Hexacentrinae), which was congruent with the phylogenetic results using site-homogeneous both in ML and BI analyses (
The generic relationships within the subfamily Meconematinae were not identical in the BI and ML trees. The relationships between the four genera were as follows: (((Pseudocosmetura, Decma) Pseudokuzicus) Xizicus) in the BI tree (Fig.
This study is supported by the National Natural Science Foundation of China (No. 31601887, 31672259), Scientific Research Program of Shaanxi Province (No. 2015KJXX-90).