Mitochondrial phylogenomics and genetic relationships of closely related pine moth (Lasiocampidae: Dendrolimus) species in China, using whole mitochondrial genomes
© Qin et al.; licensee BioMed Central. 2015
Received: 1 July 2014
Accepted: 23 April 2015
Published: 4 June 2015
Pine moths (Lepidoptera; Bombycoidea; Lasiocampidae: Dendrolimus spp.) are among the most serious insect pests of forests, especially in southern China. Although COI barcodes (a standardized portion of the mitochondrial cytochrome c oxidase subunit I gene) can distinguish some members of this genus, the evolutionary relationships of the three morphospecies Dendrolimus punctatus, D. tabulaeformis and D. spectabilis have remained largely unresolved. We sequenced whole mitochondrial genomes of eight specimens, including D. punctatus wenshanensis. This is an unambiguous subspecies of D. punctatus, and was used as a reference for inferring the relationships of the other two morphospecies of the D. punctatus complex. We constructed phylogenetic trees from this data, including twelve published mitochondrial genomes of other Bombycoidea species, and examined the relationships of the Dendrolimus taxa using these trees and the genomic features of the mitochondrial genome.
The eight fully sequenced mitochondrial genomes from the three morphospecies displayed similar genome structures as other Bombycoidea species in terms of gene content, base composition, level of overall AT-bias and codon usage. However, the Dendrolimus genomes possess a unique feature in the large ribosomal 16S RNA subunits (rrnL), which are more than 60 bp longer than other members of the superfamily and have a higher AC proportion. The eight mitochondrial genomes of Dendrolimus were highly conservative in many aspects, for example with identical stop codons and overlapping regions. But there were many differences in start codons, intergenic spacers, and numbers of mismatched base pairs of tRNA (transfer RNA genes).
Our results, based on phylogenetic trees, genetic distances, species delimitation and genomic features (such as intergenic spacers) of the mitochondrial genome, indicated that D. tabulaeformis is as close to D. punctatus as is D. punctatus wenshanensis, whereas D. spectabilis evolved independently from D. tabulaeformis and D. punctatus. Whole mitochondrial DNA phylogenies showed that D. spectabilis formed a well-supported monophyletic clade, with a clear species boundary separating it from the other congeners examined here. However, D. tabulaeformis often clustered with D. punctatus and with the subspecies D. punctatus wenshanensis. Genetic distance analyses showed that the distance between D. tabulaeformis and D. punctatus is generally less than the intraspecific distance of D. punctatus and its subspecies D. punctatus wenshanensis. In the species delimitation analysis of Poisson Tree Processes (PTP), D. tabulaeformis, D. punctatus and D. punctatus wenshanensis clustered into a putative species separated from D. spectabilis. In comparison with D. spectabilis, D. tabulaeformis and D. punctatus also exhibit a similar structure in intergenic spacer characterization. These different types of evidence suggest that D. tabulaeformis is very close to D. punctatus and its subspecies D. punctatus wenshanensis, and is likely to be another subspecies of D. punctatus.
Whole mitochondrial genomes possess relatively rich genetic information compared with the traditional use of single or multiple genes for phylogenetic purposes. They can be used to better infer phylogenetic relationships and degrees of relatedness of taxonomic groups, at least from the aspect of maternal lineage: caution should be taken due to the maternal-only inheritance of this genome. Our results indicate that D. spectabilis is an independent lineage, while D. tabulaeformis shows an extremely close relationship to D. punctatus.
KeywordsMitochondrial genome Dendrolimus punctatus D. tabulaeformis D. spectabilis Subspecies
Insect mitochondrial genomes are usually a closed-circular molecule 14-20 kilo base pairs (kbp) in size . They contain 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and a large non-coding region called the control region (CR, or the AT-rich region) which includes the essential regulatory elements for transcription and replication [1,2]. Due to their unique characters, which include small genome size, easily accessible nature, faster nucleotide substitution rates and the presence of strictly orthologous genes, mitochondrial genomes have been widely used as molecular markers for phylogenetic analyses and the investigation of questions concerning comparative and evolutionary genomics [3-10]. Whole mitochondrial genomes provide not only more general genetic information than shorter sequences of individual genes such as the COI gene, but also sets of genome-level characters, such as the relative position of different genes, structural genomic features and compositional features [1,2,11,12]. Whole mitochondrial genomes are widely used to infer phylogenetic relationships at different taxonomic levels [7,13-18]. Furthermore, mitochondrial genomes often evolve faster through higher mutation rates than nuclear genomes, especially in intergenic regions [19,20], and can be powerful markers for the inference of phylogenetic relationships among closely related taxa [9,21-24].
Dendrolimus species (Lepidoptera, Lasiocampidae) are the most serious phytophagous pests worldwide. Dendrolimus pini causes destructive damage of the Scots pine (Pinus sylvestris) in Europe , and Dendrolimus moths cause extensive forest damage in China [26-31]. The occurrence of different species of these pine moths and the existence of natural hybrids leads to heterosis and strong tolerance to pesticides, giving high survival rates and making them difficult to eradicate . The need for improved biological pest control means that it is very important to elucidate the genetic relationships of pine moth species. However, three (D. punctatus , D. tabulaeformis , D. spectabilis ) of the six commonly occurring Dendrolimus species in China cannot be readily discriminated, and their taxonomic status is uncertain. Previous studies revealed that these three can hybridize, and it has been suggested that D. tabulaeformis and D. spectabilis are subspecies of D. punctatus [26,27], although other studies treat all three as different species [29-31]. A recent extensive molecular systematic study, looking at mitochondrial COI and several nuclear genes, used different phylogenetic and DNA barcoding methods to assess this complex. The phylogenetic relationships of D. punctatus, D. tabulaeformis, D. spectabilis could not be fully resolved, although their close relationship was confirmed . This species group thus provides a good model for investigating the utility of the mitochondrial genome in exploring relationships of closely related species groups.
To clarify the taxonomic and phylogenetic relationships of these species, we report here eight complete mitochondrial genomes of the three morphospecies Dendrolimus punctatus, D. tabulaeformis and D. spectabilis. We compare their genome structures in detail and analyze the relationships between them.
Specimen collection, DNA extraction, PCR amplification, cloning and sequencing
Eight adults of three morphospecies (D. punctatus (2), D. punctatus wenshanensis (2), D. tabulaeformis (2) and D. spectabilis (2)) were collected from five locations (Additional file 1). A taxonomic expert (Professor Chun-sheng Wu, Institute of Zoology, Chinese Academy of Sciences) used traditional morphological approaches to identify these specimens. All specimens were preserved in 95% ethanol and maintained at 4°C for long term storage. DNA was isolated from thoracic muscle tissue using the DNeasy Blood and Tissue kit (QIAGEN) following the manufacturer’s protocol.
Mitochondrial genomes were PCR-amplified and sequenced by Sanger sequencing. Whole genomes of the eight specimens were mainly assembled from 14 overlapping PCR fragments (Additional file 2). Eight pairs of universal primers were selected and used for initial amplifications [3,4]. From these sequences, six species-specific primer sets were designed to amplify remaining sections. A specific fragment of DNA was amplified using universal primers by the following conditions: an initial denaturation at 94°C for 4 min, followed by 40 cycles of denaturation at 94°C for 30 s, annealing at 45–55°C for 30 s, elongation at 72°C for 45 s, and a final extension step of 72°C for 10 min. Species-specific primers were amplified with long PCR reaction conditions: an initial denaturation at 92°C for 2 min, followed by 40 cycles of denaturation at 92°C for 30 s, annealing at 45–55°C for 30 s, elongation at 60°C for 12 min, and a final extension step of 60°C for 20 min. All reactions were performed using Takara LA taq (TaKaRa Co., Dalian, China). PCR fragments containing the control region could not be sequenced directly; they were cloned into the pEASY-T3 Cloning Vector (Beijing TransGen Biotech Co., Ltd., Beijing, China) and then sequenced by M13-F (CGCCAGGGTTTTCCCAGTCACGAC) and M13-R (GAGCGGATAACAATTTCACACAGG) primers.
Sequence assembly and annotation
Raw sequences were checked manually by eye using BioEdit  and assembled using the SeqMan program from the laser gene package DNAStar (Madison, USA). The tRNA genes were identified by tRNAscan-SE Search Server v.1.21 , with source set for Mito/Chloroplast and Genetic Code for tRNA Isotype Prediction set for Invertebrate Mito. To detect overall tRNA genes, the cove cutoff score could be set lower (≥1). The locations of protein-coding and rRNA genes were determined by comparing homologous sequences with other published Lepidoptera mitochondrial sequences downloaded from GenBank.
Evolutionary relationships based on phylogenetic analysis
There has been much controversy concerning the evolutionary relationships of closely related species of Dendrolimus [26-31]. The three morphospecies D. spectabilis, D. punctatus and D. tabulaeformis can hybridize with each other. Furthermore, D. punctatus and D. tabulaeformis have highly similar diets and morphological features. Different geographical environments and host plants can result in varied morphology, and D. spectabilis and D. tabulaeformis have been treated as subspecies of D. punctatus [26-28]. We assessed their genetic relationships using whole mitochondrial genome data (including genomic features), with D. punctatus wenshanensis (an unambiguous subspecies of D. punctatus) as a baseline reference to infer relationships between the two ambiguous morphospecies (D. tabulaeformis and D. spectabilis) and D. punctatus.
List of taxa analyzed in this study
Lee et al., unpulished
Samia cynthia ricini
In this study
In this study
In this study
In this study
In this study
In this study
Dendrolimus punctatus wenshanensis03
In this study
Dendrolimus punctatus wenshanensis06
In this study
The ML analyses were conducted using RaxML  with the selected partition scheme (Additional file 3). Because different partition models cannot be set in RaxML, GTRGAMMAI was selected for both the 13PCGs dataset and the 37genes dataset. In RaxML, the rapid hill-climbing algorithm starting from 100 randomized maximum-parsimony trees was used for ML searches. Confidence values of the ML trees were evaluated via bootstrap tests with 1000 iterations.
The BI analyses were conducted using MrBayes ver.3.1  with the best partition scheme (Additional file 3). All BI analyses were conducted with two sets of four Markov chains (one cold and three hot chains). Three criteria were used to verify convergence of BI analyses in order to be confident about the reliability of the results: (1) the average standard deviation of split frequencies was less than 0.01, (2) the potential scale reduction factors (PSRF)  were close to 1.0 for all parameters, and (3) the effective sample size value in tracer software exceeded 200 . Each set was sampled with a burnin of 25%. Bayesian posterior probabilities were estimated as the confidence values of the BI tree (BBP).
Evolutionary relationships based on genome analysis and genetic distance calculations
Genome analysis was performed with the 20 mitochondrial sequences of Bombycoidea (Table 1). Nucleotide compositions, codon usage (excluding stop codons) and Relative Synonymous Codon Usage (RSCU) were calculated with MEGA 4.0 . Composition skew analysis used the formulae AT skew = [A-T]/[A + T] and GC skew = [G-C]/[G + C] .
A second investigation used the COI gene alone, but from much larger sample sizes ((D. punctatus (195), D. tabulaeformis (6), D. spectabilis (14), D. punctatus wenshanensis (35)) and different geographical populations. K2P genetic distances were calculated with a Perl script . Confidence intervals were calculated using R with t.test function in a stats package ( http://cran.r-project.org/).
Evolutionary relationships based on species delimitation
We applied the Poisson Tree Processes (PTP) model to delimit species on a rooted phylogenetic tree for our two datasets (13PCGs and 37gene) . The PTP program models speciation or coalescent events relative to number of substitutions represented by branch lengths, and uses heuristic algorithms to estimate the most likely classification of branches into population and species-level processes. It does not require ultrametrization of trees and outperforms the commonly used General Mixed Yule Coalescent (GMYC) model when species are closely related and evolutionary distances small [52,53]. The BI phylogenetic trees obtained with MrBayes 3.1 , described previously, were used as input trees. PTP analyses were run from Python using the ETE (Python Environment for Tree Exploration) package  for tree manipulation and visualization.
Mitochondrial genome organization
Eight complete mitochondrial genomes from Dendrolimus were sequenced, annotated and deposited in GenBank (Table 1). Genome organizations are presented in an additional file (Additional file 5 a-h). All genomes contained the same set of 37 genes (13 protein-coding genes, 22 tRNA genes, and two rRNA genes) and a putative control region . Genome sizes ranged from 15,407 bp to 15,419 bp, and all exhibited similar sequence characteristics. The gene order of these eight pine moths was identical to other ditrysian lepidopterans with the trnM gene location type (trnM-trnI-trnQ) [55-57]. All 13 PCGs start with typical ATN codons, except for cox1 which used CGA. Nine genes (nad2, atp8, atp6, cox3, nad5, nad4l, nad6, cob, nad1) share the same complete termination codon TAA, and four genes use incomplete stop codons (a single T for cox1, cox2 and nad4, TA for nad3) (Additional file 6). Incomplete stop codons are common in lepidopteran mitochondrial genomes and are presumed to be completed via post-transcriptional polyadenylation [58,59].
As for other Lepidoptera, 22 tRNA genes were detected in Dendrolimus mitochondrial genomes, and ranged in size from 64 to 71 bp. Fourteen tRNA genes were coded on the majority strand (J-strand) with eight coded on the minority strand (N-strand). All the tRNA genes have classic cloverleaf secondary structures except for the trnS1 (AGN) gene, where the dihydrouridine (DHU) arm is replaced by unpaired nucleotides. This feature is common to most Lepidopteran mitochondrial genomes except the tortricid Adoxophyes honmai, which has all tRNA genes with complete cloverleaf structures . These aberrant tRNA genes may be modified via RNA-editing mechanisms [61,62]. Many mismatches were found in each of the newly sequenced genomes resulting in some small differences (Additional file 7). There were 16 identical mismatched base pairs and G-U wobble pairs in tRNA structures for all eight samples.
Base composition, AT-skew and codon usage were calculated for the eight newly sequenced genomes and the 12 published Bombycidae mitochondrial genomes (Table 1). The base composition of the J-strand of the Dendrolinus species fell within the range of other Bombycoidea species, and their nucleotide compositions were significantly biased toward A and T (Additional file 8). The nucleotide skew statistics for the entire majority strand of eight individuals indicated a weak A skew and a moderate C skew. In the protein coding genes, both the second and the third positions had negative AT-skew and GC-skew, with the first position having positive AT-skew and GC-skew; overall, PCGs had a negative AT-skew and a positive GC-skew (Additional file 8). The nucleotide skew statistics for the tRNA genes of the eight Dendrolimus individuals indicated a weak A skew and a moderate C skew. The rrnL gene of Dendrolimus is more than 60 bp longer than that of other Bombycoidea species, and the nucleotide skew statistics for rrnL of the eight individuals showed higher A- and C-skews than other Bombycoidea.
Codon usage and RSCU results of our eight moths and the other Bombycoidea species are also given (Additional file 9). Codon usage bias has been detected in currently sequenced mitochondrial genomes. Excluding stop codons, there were 3726 codons in both D. spectabilis and D. tabulaeformis,. One additional codon was detected in D. punctatus. The codon families show very similar behavior among all 20 species considered (Additional file 9). In the eight newly sequenced genomes, eight codon families (Phe, Leu2, Ile, Met, Ser2, Tyr, Asn, Gly) had at least 50 codons per thousand codons (CDs), and two codon families (Leu2, Ile) had at least 100 codons per thousand CDs; other Bombycoidea species had three codon families (Phe, Leu2, Ile) with at least 100 codons per thousand CDs.
RSCU results of Bombycoidea showed that codons exhibit a strong AT-bias in the third codon position, and GC-rich codons are less preferred among the usage of both four- and twofold degenerate codons  (Additional file 9). In our eight specimens, many codons were absent, with AGG absent in D. punctatus wenshanensis03, D. punctatus wenshanensis 06, D. punctatus04, D. tabulaeformis06 and D. tabulaeformis38; UGC and AGG were absent in D. spectabilis02 and D. spectabilis13 (Additional file 9).
Evolutionary relationships based on phylogenetic analyses
Evolutionary relationships based on genetic distance: individual genes and regions versus whole mitochondrial genome analysis
In most genes and regions (such as intergenic spacers, tRNA genes, rRNA genes, whole genome, atp6, cox1, nd2, nd4, nd4l, nd5 and nd6), the average K2P distances of D. punctatus to D. tabulaeformis were less than the K2P intraspecies distances of D. punctatus to D. punctatus or the distances of D. punctatus to its subspecies D. punctatus wenshanensis (Additional file 4). But for the atp8, cox2, cox3, cytb, nd1 and nd3 genes, the average genetic distances of D. punctatus to D. tabulaeformis were higher than the average intraspecific distances of D. punctatus and distances of D. punctatus to D. punctatus wenshanensis (Figure 1). However, for most regions and genes, especially for the complete genome, genetic distances between D. punctatus and D. tabulaeformis were less than the intraspecific differences of D. punctatus. These findings suggest that D. punctatus and D. tabulaeformis are extremely similar. The average K2P interspecific distances of D. punctatus and D. spectabilis, based on whole genomes or genes or regions, were higher than the K2P intraspecific distances of D. punctatus and distances of D. punctatus to D. punctatus wenshanensis (Figure 1). These results therefore suggest that D. spectabilis is a relatively distinct lineage from D. punctatus and D. tabulaeformis.
However, there were some differences between the two datasets, the one using just eight samples but sequenced for the entire cox1 gene and the other examining many more samples but only for the COI region of cox1. In the mitochondrial genome dataset, the results from the cox1 gene show that the average genetic distance of D. punctatus to D. tabulaeformis was higher than the average genetic distance of D. punctatus to D. punctatus wenshanensis but less than the intraspecific distance of D. punctatus. In the larger COI dataset, the average interspecific distance of D. punctatus and D. tabulaeformis was higher than average intraspecific distance of D. punctatus but less than the genetic distance of D. punctatus to D. punctatus wenshanensis. There are two possible causes of this discrepancy: firstly, sample sizes are small for the cox1 dataset so that confidence intervals cannot be estimated and compared, and secondly, the cox1 gene in the mitochondrial genome dataset was 1,531 bp long, much longer than the length of the COI barcode sequence examined here (588 bp).
Evolutionary relationships based on species delimitation
Evolutionary relationships with a mitochondrial genome feature: intergenic spacers
Mitochondrial genomes possess many organizational characteristics. The intergenic regions, for example, do not encode genes and as a result are under more relaxed selection and evolve faster than coding genes [19,20]. Consequently, intergenic regions are expected to be more variable than coding genes, and therefore could be a useful feature when studying the relationships of closely related species.
There were seventeen intergenic spacers in the eight sequenced mitochondrial genomes (Additional file 12). Six were longer than 10 bp: trnM-nd2, trnY-cox1, atp6-cox3, trnA-trnR, trnN-trnS (AGN) and nd4-nd4l. Comparing these six regions of the eight individuals, many differences can be observed (Additional file 13). The intergenic spacer trnM-nd2 was the largest at 58 bp, and is common to most lepidopteran mitochondrial genomes. In this spacer, the two individuals of D. spectabilis had identical nucleotide patterns. Except for one specimen, D. punctatus, D. punctatus wenshanensis and D. tabulaeformis were also identical, but differed by two nucleotides from D. spectabilis. The exception was D. punctatus04, with two different bases. In the other five intergenic spacers, the two D. spectabilis always showed an identical nucleotide distribution, which differed from those of D. punctatus, D. punctatus wenshanensis and D. tabulaeformis. The latter taxa present a mainly similar nucleotide pattern, but D. punctatus04 and D. punctatus05 often showed some differences. In the trnY-cox1 region, D. punctatus05 has 11 nucleotide differences, and for the atp6-cox3, trnA-trnR, trnN-trnS(AGN) and nd4-nd4l regions, D. punctatus04 showed one to three nucleotides difference. D. tabulaeformis always showed an identical nucleotide pattern to D. punctatus wenshanensis. These findings further support the conclusions that D. tabulaeformis has as close a relationship with D. punctatus as has D. punctatus wenshanensis, and that D. spectabilis is distinct.
Conclusions and Discussion
In this study, whole mitochondrial genomes were sequenced using conventional Sanger sequencing with the Primer-Walking method being applied to long fragment sequencing. Sequences were corrected manually to improve the accuracy of sequence determination, and assembled using overlap regions to obtain full-length mitochondrial genomes. Even though earlier laboratory hybridization studies suggested that both the current morphospecies D. spectabilis and D. tabulaeformis are subspecies of D. punctatus [26,27], there is as yet no consensus on the exact relationships of the three taxa [29-31]. After assessing the whole mitochondrial genome data in multiple ways, we find that D. spectabilis forms an independent lineage from D. punctatus, while D. tabulaeformis has a very close relationship with D. punctatus. Like D. punctatus wenshanensis, D. tabulaeformis is likely to be a subspecies of D. punctatus, and this conclusion is also supported by an earlier DNA barcoding based study .
Geographically, D. spectabilis is distributed in northern China, and its range does not overlap with D. punctatus (central and southern China). Both genetic and distributional data strongly suggest that D. spectabilis has evolved as a separate lineage or species, perhaps recently. In contrast, D. tabulaeformis is mostly distributed in central China and geographically has not totally separated from D. punctatus. The eight specimens that formed the basis of this study were selected from different locations (Additional file 1): two individuals of D. punctatus from Baise, Guangxi, two D. punctatus wenshanensis from Wenshan Shilin, Yunnan, two D. tabulaeformis from Chengde, Hebei, D. spectabilis02 from Tongliaokulun, Neimenggu and D. spectabilis13 from Taian, Shandong. D. punctatus, D. punctatus wenshanensis and D. tabulaeformis came from different sampling sites, but still have very high sequence similarity, and genetically cannot be unambiguously distinguished. However, the two individuals of D. spectabilis, although from different sampling sites, showed a high degree of sequence identity to one another and significant differences from D. punctatus, D. punctatus wenshanensis and D. tabulaeformis.
The main goal of this study was to explore the phylogenetic relationships of closely related species in the genus Dendrolimus. The mitochondrial DNA genome was selected because of its relatively faster rate of evolution compared with nuclear markers. However, we note that our study does have several limitations. Firstly, the sample size of each species is relatively small compared with studies where a few molecular markers but many specimens are used. We only studied the whole mitochondrial genome from two specimens of each taxon due to high sequencing costs. An increased sample size would be desirable in future studies. However, we did analyse many more individuals using the COI barcoding region alone, coming to similar conclusions. Secondly, the mitochondrial genome only represents the evolution of maternal lineages, and misleading conclusions may be drawn when there are inconsistencies between mitochondrial and nuclear genealogy. Therefore, the introduction of many, not just one or two, nuclear markers in future studies is also highly desirable. Furthermore, with the development of next-generation sequencing technology, more and more mitochondrial genomes and transcriptomes are being examined to infer phylogenetic relationships [65-67]. We plan to include nuclear markers, including transcriptome analysis, in a future study of this taxon group.
This research was supported by China National Funds for Distinguished Young Scientists (to Zhang, Grant No. 31425023), by the Natural Science Foundation of China (to Zhang, Grant No. 31071963 and 31272340), by a Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (to Zhang, Grant No. PHR201107120), by the National High-tech Research and Development Project (863) of China (2012AA021601), the Program for Changjiang Scholars and Innovative Research Team in University (IRT13081), and a Science and Technology Foundation Project (2012FY110803).
- Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27(8):1767–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Taanman J-W. The mitochondrial genome: structure, transcription, translation and replication. Biochimica et Biophysica Acta (BBA)-BioenergeticS. 1999;1410(2):103–23.View ArticleGoogle Scholar
- Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach AT: Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annual Rev Ecol Evol Systemat. 2006;37:545-579.Google Scholar
- 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.View ArticleGoogle Scholar
- Salvato P, Simonato M, Battisti A, Negrisolo E. The complete mitochondrial genome of the bag-shelter moth Ochrogaster lunifer (Lepidoptera, Notodontidae). BMC Genomics. 2008;9(1):331.PubMed CentralPubMedView ArticleGoogle Scholar
- Reyes A, Gissi C, Pesole G, Saccone C. Asymmetrical directional mutation pressure in the mitochondrial genome of mammals. Mol Biol Evol. 1998;15(8):957–66.PubMedView ArticleGoogle Scholar
- Ingman M, Kaessmann H, PaÈaÈbo S, Gyllensten U. Mitochondrial genome variation and the origin of modern humans. Nature. 2000;408(6813):708–13.PubMedView ArticleGoogle Scholar
- Lin C-P, Danforth BN. How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets. Mol Phylogenet Evol. 2004;30(3):686–702.PubMedView ArticleGoogle Scholar
- Gissi C, Iannelli F, Pesole G. Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species. Heredity. 2008;101(4):301–20.PubMedView ArticleGoogle Scholar
- Elson JL, Lightowlers RN. Mitochondrial DNA clonality in the dock: can surveillance swing the case? Trends Genet. 2006;22(11):603–7.PubMedView ArticleGoogle Scholar
- Cameron SL. Insect mitochondrial genomics: implications for evolution and phylogeny. Annu Rev Entomol. 2014;59:95–117.PubMedView ArticleGoogle Scholar
- Clayton DA. Replication of animal mitochondrial DNA. Cell. 1982;28(4):693–705.PubMedView ArticleGoogle Scholar
- Cameron SL, Lambkin CL, Barker SC, Whiting MF. A mitochondrial genome phylogeny of Diptera: whole genome sequence data accurately resolve relationships over broad timescales with high precision. Syst Entomol. 2007;32(1):40–59.View ArticleGoogle Scholar
- Dowton M, Cameron SL, Austin AD, Whiting MF. Phylogenetic approaches for the analysis of mitochondrial genome sequence data in the Hymenoptera–a lineage with both rapidly and slowly evolving mitochondrial genomes. Mol Phylogenet Evol. 2009;52(2):512–9.PubMedView ArticleGoogle Scholar
- Fenn JD, Song H, Cameron SL, Whiting MF. A preliminary mitochondrial genome phylogeny of Orthoptera (Insecta) and approaches to maximizing phylogenetic signal found within mitochondrial genome data. Mol Phylogenet Evol. 2008;49(1):59–68.PubMedView ArticleGoogle Scholar
- Nardi F, Spinsanti G, Boore JL, Carapelli A, Dallai R, Frati F. Hexapod origins: monophyletic or paraphyletic? Science. 2003;299(5614):1887–9.PubMedView ArticleGoogle Scholar
- Boore JL, Brown WM. Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998;8(6):668–74.PubMedView ArticleGoogle Scholar
- Shao R, Barker S. Mitochondrial genomes of parasitic arthropods: implications for studies of population genetics and evolution. Parasitology. 2007;134(02):153–67.PubMedView ArticleGoogle Scholar
- Burger G, Gray MW, Franz Lang B. Mitochondrial genomes: anything goes. Trends Genet. 2003;19(12):709–16.PubMedView ArticleGoogle Scholar
- Ghikas DV, Kouvelis VN, Typas MA. Phylogenetic and biogeographic implications inferred by mitochondrial intergenic region analyses and ITS1-5.8 S-ITS2 of the entomopathogenic fungi Beauveria bassiana and B. brongniartii. BMC Microbiol. 2010;10(1):174.PubMed CentralPubMedView ArticleGoogle Scholar
- Torriani SF, Penselin D, Knogge W, Felder M, Taudien S, Platzer M, et al. Comparative analysis of mitochondrial genomes from closely related Rhynchosporium species reveals extensive intron invasion. Fungal Genet Biol. 2014;62:34–42.PubMedView ArticleGoogle Scholar
- Arctander P. Comparison of a mitochondrial gene and a corresponding nuclear pseudogene. Proc R Soc London, Ser B. 1995;262(1363):13–9.View ArticleGoogle Scholar
- Krzywinski J, Grushko OG, Besansky NJ. Analysis of the complete mitochondrial DNA from Anopheles funestus: An improved dipteran mitochondrial genome annotation and a temporal dimension of mosquito evolution. Mol Phylogenet Evol. 2006;39(2):417–23.PubMedView ArticleGoogle Scholar
- Lopez JV, Cevario S, O'Brien SJ. Complete Nucleotide Sequences of the Domestic Cat (Felis catus) Mitochondrial Genome and a Transposed mtDNA Tandem Repeat (Numt) in the Nuclear Genome. Genomics. 1996;33(2):229–46.PubMedView ArticleGoogle Scholar
- le Mellec A, Michalzik B. Impact of a pine lappet (Dendrolimus pini) mass outbreak on C and N fluxes to the forest floor and soil microbial properties in a Scots pine forest in Germany. Can J Forest Res. 2008;38(7):1829–41.View ArticleGoogle Scholar
- Zhao QS, Wu WB, Lv GP, Yuan X, Li SK. Hybridization among pine moths species. Acta Entomol Sin. 1992;35(1):28–32.Google Scholar
- Zhao QS, Wu WB, Lv GP, Chen TF, Lin QY. Hybridization among pine moths species and its genetic rule. Scientia Silvae Sinicae. 1999;35(4):45–50.Google Scholar
- Zhang AB, Wang ZJ, Tan SJ, Li DM. Monitoring the masson pine moth, Dendrolimus punctatus Walker Lepidoptera: Lasiocampidae with synthetic sex pheromone-baited traps in Qianshan County. China JPN J Appl Entomol Z. 2003;38(2):177–86.View ArticleGoogle Scholar
- Zhang AB, Kong XB, Li DM. DNA finger printing evidence for the phylogenetic relationship of eight species and subspecies of Dendrolimus (Lepidoptera: Lasiocampidae) in China. Acta Entomol Sin. 2004;47(2):236–42.Google Scholar
- Zhang AB, Zhang Z, Wang HB, Kong XB. Geographical distribution of Lasiocampidae in China and its relationship with environmental factors. J Beijing Forestry Uni. 2004;26(4):54–60.Google Scholar
- Zhang AB, Li DM, Chen J. Geohistories of Dendrolimus punctatus Walker and its host plant pine genus in China. Chinese Bulletin Entomol. 2004;2:146–50.Google Scholar
- Walker F, Gray JE. List of the specimens of lepidopterous insects in the collection of the British Museum: order of the Trustees. 1855.Google Scholar
- Tsai P, Liu Y. A study of the genus Dendrolimus (Lasiocampidae) of China with descriptions of two new species and one new subspecies. Acta Entomol Sin. 1962;11(3):237–52.Google Scholar
- Butler A. Descriptions of new species of Heterocera from Japan, Part 3, Geometries. Ann Mag Nat Hist. 1878;5:392–419.View ArticleGoogle Scholar
- Dai Q-Y, Gao Q, Wu C-S, Chesters D, Zhu C-D, Zhang A-B. Phylogenetic reconstruction and DNA barcoding for closely related pine moth species (Dendrolimus) in China with multiple gene markers. PLoS One. 2012;7(4), e32544.PubMed CentralPubMedView ArticleGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic acids symposium series. 1999;1999:95–8.Google Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):0955–64.View ArticleGoogle Scholar
- Yang X, Xue D, Han H. The complete mitochondrial genome of Biston panterinaria (Lepidoptera: Geometridae), with phylogenetic utility of mitochondrial genome in the Lepidoptera. Gene. 2013;515(2):349–58.PubMedView ArticleGoogle Scholar
- Yang L, Wei Z-J, Hong G-Y, Jiang S-T, Wen L-P. The complete nucleotide sequence of the mitochondrial genome of Phthonandria atrilineata (Lepidoptera: Geometridae). Mol Biol Rep. 2009;36(6):1441–9.PubMedView ArticleGoogle Scholar
- Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC bioinformatics. 2004;5(1):113.PubMed CentralPubMedView ArticleGoogle Scholar
- Vaidya G, Lohman DJ, Meier R. SequenceMatrix: concatenation software for the fast assembly of multi‐gene datasets with character set and codon information. Cladistics. 2011;27(2):171–80.View ArticleGoogle Scholar
- Lanfear R, Calcott B, Ho SY, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29(6):1695–701.PubMedView ArticleGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.PubMedView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572–4.PubMedView ArticleGoogle Scholar
- Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences. Stat Sci 1992;7:457-472Google Scholar
- Rambaut A, Drummond A: Tracer v1. 4. In.; 2007 [http://beast.bio.ed.ac.uk/Tracer].
- Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9.PubMedView ArticleGoogle Scholar
- Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995;41(3):353–8.PubMedView ArticleGoogle Scholar
- Hebert PD, Cywinska A, Ball SL. Biological identifications through DNA barcodes. Proc R Soc London, Ser B. 2003;270(1512):313–21.View ArticleGoogle Scholar
- Hebert PD, Ratnasingham S, de Waard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc London, Ser B. 2003;270 Suppl 1:S96–9.View ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29(22):2869–76.PubMed CentralPubMedView ArticleGoogle Scholar
- Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, et al. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006;55(4):595–609.PubMedView ArticleGoogle Scholar
- Huerta-Cepas J, Dopazo J, Gabaldón T. ETE: a python Environment for Tree Exploration. BMC bioinformatics. 2010;11(1):24.PubMed CentralPubMedView ArticleGoogle Scholar
- Negrisolo E, Babbucci M, Patarnello T. The mitochondrial genome of the ascalaphid owlfly Libelloides macaronius and comparative evolutionary mitochondriomics of neuropterid insects. BMC Genomics. 2011;12(1):221.PubMed CentralPubMedView ArticleGoogle Scholar
- 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):112–23.PubMedView ArticleGoogle Scholar
- Cao Y-Q, Ma C, Chen J-Y, Yang D-R. The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera. BMC Genomics. 2012;13(1):276.PubMed CentralPubMedView ArticleGoogle Scholar
- Ojala D, Montoya J, Attardi G. tRNA punctuation model of RNA processing in human mitochondria. 1981.Google Scholar
- Ojala D, Merkel C, Gelfand R, Attardi G. The tRNA genes punctuate the reading of genetic information in human mitochondrial DNA. Cell. 1980;22(2):393–403.PubMedView ArticleGoogle Scholar
- Lee E-S, Shin KS, Kim M-S, Park H, Cho S, Kim C-B. The mitochondrial genome of the smaller tea tortrix Adoxophyes honmai (Lepidoptera: Tortricidae). Gene. 2006;373:52–7.PubMedView ArticleGoogle Scholar
- Lavrov DV, Brown WM, Boore JL. A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc Natl Acad Sci. 2000;97(25):13738–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Masta SE, Boore JL. The complete mitochondrial genome sequence of the spider Habronattus oregonensis reveals rearranged and extremely truncated tRNAs. Mol Biol Evol. 2004;21(5):893–902.PubMedView ArticleGoogle Scholar
- Chai H-N, Du Y-Z, Zhai B-P. Characterization of the complete mitochondrial genomes of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae). Int J Biol Sci. 2012;8(4):561.PubMed CentralPubMedView ArticleGoogle Scholar
- Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111–20.PubMedView ArticleGoogle Scholar
- Kawahara AY, Breinholt JW. Phylogenomics provides strong evidence for relationships of butterflies and moths. Proc R Soc Lond B Biol Sci. 2014;281(1788):20140970.View ArticleGoogle Scholar
- Timmermans MJ, Lees DC, Simonsen TJ: Towards a mitogenomic phylogeny of Lepidoptera. Mol Phylogen Evol. 2014;79:169-178.Google Scholar
- Bazinet AL, Cummings MP, Mitter KT, Mitter CW. Can RNA-Seq Resolve the Rapid Radiation of Advanced Moths and Butterflies (Hexapoda: Lepidoptera: Apoditrysia)? An Exploratory Study. PLoS One. 2013;8(12), e82615.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.