The complete mitochondrial genome of the bag-shelter moth Ochrogaster lunifer (Lepidoptera, Notodontidae)
© Salvato et al; licensee BioMed Central Ltd. 2008
Received: 21 May 2008
Accepted: 15 July 2008
Published: 15 July 2008
Knowledge of animal mitochondrial genomes is very important to understand their molecular evolution as well as for phylogenetic and population genetic studies. The Lepidoptera encompasses more than 160,000 described species and is one of the largest insect orders. To date only nine lepidopteran mitochondrial DNAs have been fully and two others partly sequenced. Furthermore the taxon sampling is very scant. Thus advance of lepidopteran mitogenomics deeply requires new genomes derived from a broad taxon sampling. In present work we describe the mitochondrial genome of the moth Ochrogaster lunifer.
The mitochondrial genome of O. lunifer is a circular molecule 15593 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. It contains also 7 intergenic spacers. The gene order of the newly sequenced genome is that typical for Lepidoptera and differs from the insect ancestral type for the placement of trnM. The 77.84% A+T content of its α strand is the lowest among known lepidopteran genomes. The mitochondrial genome of O. lunifer exhibits one of the most marked C-skew among available insect Pterygota genomes. The protein-coding genes have typical mitochondrial start codons except for cox1 that present an unusual CGA. The O. lunifer genome exhibits the less biased synonymous codon usage among lepidopterans. Comparative genomics analysis study identified atp6, cox1, cox2 as cox3, cob, nad1, nad2, nad4, and nad5 as potential markers for population genetics/phylogenetics studies. A peculiar feature of O. lunifer mitochondrial genome it that the intergenic spacers are mostly made by repetitive sequences.
The mitochondrial genome of O. lunifer is the first representative of superfamily Noctuoidea that account for about 40% of all described Lepidoptera. New genome shares many features with other known lepidopteran genomes. It differs however for its low A+T content and marked C-skew. Compared to other lepidopteran genomes it is less biased in synonymous codon usage. Comparative evolutionary analysis of lepidopteran mitochondrial genomes allowed the identification of previously neglected coding genes as potential phylogenetic markers. Presence of repetitive elements in intergenic spacers of O. lunifer genome supports the role of DNA slippage as possible mechanism to produce spacers during replication.
List of taxa analyzed in present paper
Hang and Zhang, unpublished
Baumann and Baumann, unpublished
Lepidopsocid sp. RS-2001
An et al., unpublished
Liu et al., unpublished
Lee et al., unpublished
Yang et al., unpublished
Hong et al., unpublished
Ho et al., unpublished
Lobo et al., unpublished
Anopheles quadrimaculatus A
Cibrario et al., unpublished
Azeredo-Espin et al., unpublished
Lessinger et al., unpublished
Ye et al., unpublished
Yu et al., unpublished
Ye et al., unpublished
Ye et al., unpublished
Apis mellifera ligustica
Silvestre and Arias, unpublished
Results and discussion
Genome organization, structure and composition
The O. lunifer mtDNA has the typical lepidopteran gene order [8, 9] that differs from the ancestral gene order of insects  for the placement of trnM. In the ancestral type (e.g. Drosophila yakuba mtDNA) the order in the α strand is: A+T region, trnI, trnQ, trnM, nad2. In all lepidopteran mtDNAs, sequenced to date, the order is: A+T region, trnM, trnI, trnQ, nad2 which implies the translocation of trnM [5–11]. This placement of trnM is a molecular feature exclusive to lepidopteran mtDNAs. Further genome sequencing is necessary to establish if this feature is a mitochondrial signature of the whole order Lepidoptera.
The composition of the α strand of O. lunifer mtDNA is A = 6252 (40.09%), T = 5886 (37.75%), G = 1179 (7.56%) and C = 2276 (14.60%).
The average A+T% value for the analyzed mtDNAs set is 76. 63 ± 4.84. The highest A+T% values are shared by the mtDNAs of three bees (Apis mellifera, Bombus ignitus and Melipona bicolor) and two bugs (Aleurodicus dugesii and Schizaphis graminum). All lepidopteran mtDNAs but O. lunifer exhibit high A+T% values. The A+T content of O. lunifer mtDNA is 77.84% that represents the lowest value for lepidopteran complete mtDNAs [5–8, 10]. The lowest A+T contents are found in the termite mtDNAs (Reticulitermes spp.). Extreme A+T values are also shared by species having highly re-arranged gene order . However the possession of a re-arranged genome is not sufficient per se to have an A+T content drastically departing from the average (e.g. Aleurochiton aceris and Bemisia tabaci). The A+T values appear to be linked to taxonomic relatedness at low rank (i.e. genus, family) (e.g. species of Drosophila, species of Bactrocera, members of family Apidae). The relation is not true at higher ranks (i.e. superfamily; order) where patterns become inconsistent and the A+T content can be very different among species as exemplified by Hemiptera (A. dugesii vs.Triatoma dimidiata).
The average A-skew is 0.04214 ± 0.11350 and most of pterygote mtDNAs are slightly to moderately A-skewed with values ranging from 0.00287 (B. ignitus) to 0.18247 (Locusta migratoria). The lepidopteran A-skews vary from -0.04748 (C. raphaelis) to 0.05872 (B. mori) with the O. lunifer mtDNA exhibiting a slight A-skew (0.03015). The Reticulitermes mtDNA genomes, having the lowest A+T% content, exhibit a very pronounced A-skew. Most marked T-skews are observed in the mtDNA genomes of Campanulotes bidentatus and Trialeuroides vaporarium that have low A+T% content and gene-orders different than insect ancestral gene order [1, 14, 15]. Gene order re-arrangement is not necessarily linked to strong A/T-skew as proved by the highly rearranged, but low skew, genome of Heterodoxus macropus .
The average G+C% content is 23.37 ± 4.84. The G+C% pattern among various species is obviously opposite to the A+T% thus it does not require further comments. More composite is the G/C-skew distribution. The average G-skew is -0.16006 ± 0.138235. Most of pterygota mtDNAs are C-skewed with G-skew values ranging from -0.32827 (Vanhornia eucnemidarum) to -0.01250 (Heterodoxus macropus). The main exception is represented by the mtDNA of bugs, while the highest G-skewed genome is that of C. bidentatus. Most of lepidopteran mtDNAs share very similar G-skew values that are included within the bulk of mtDNAs. The notable exception is represented by the newly determined mtDNA of O. lunifer that exhibits the second most pronounced C-skew (G-skew = -0.31751) among analyzed genomes.
G-skew can be markedly different even in species belonging to the same genus and having a very similar G+C content as well exemplified by Reticulitermes santonensis and Reticulitermes virginicus mtDNAs. The same reasoning applies at high taxonomic rank to the Hemiptera. The mtDNA of C. bidentatus exhibits very high A-skew and G-skew. However, this feature is not a general rule and extreme A-skew and G-skew are not necessarily reciprocally linked, as proved by species of genus Reticulitermes that exhibit very strong A-skews but not G-skews.
The list of currently available mtDNAs reveals that there is a strong bias in term of taxon sampling both at low and high taxonomic ranks within Pterygota. A direct consequence is that present knowledge of base composition and A/G skews reflects such biases and addition of a single taxon can change our view on these features. This point is well exemplified by the O. lunifer mtDNA that exhibits a A+T percentage different than other lepidopteran mtDNAs that share high A+T contents [8, 9]. Thus a broad and more balanced taxon sampling appears to be a mandatory goal to investigate and identify general patterns for the parameters considered above.
The mtDNA of O. lunifer contains the full set of PCGs usually present in animal mtDNA. PCGs are arranged along the genome according to the standard order of Insects  (Figure 1). The putative start codons of PCGs are those previously known for animal mtDNA i.e. ATN, GTG, TTG, GTT  with the only exception represented by the CGA start codon of cox1 gene. This non-canonical putative start codon is found also in the butterfly A. melete and in the moths A. honmai, B. mori, B. mandarina, M. sexta and P. atrilineata [5–8]. In the butterfly C. raphaelis the tetranucleotide TTAG is the putative start codon  and the six nucleotide TATTAG has been suggested as putative start codon for the moths O. nubilalis and O. furnicalis . An unusual start codon for cox1 gene is known in various arthropod mtDNA [e.g. ].
The cox1, cox2, nad5, and nad4 genes of O. lunifer mtDNA have incomplete stop codons. The presence of incomplete stop codons is a feature shared with all lepidopteran mtDNAs sequenced to date [5–10] and more in general with many arthropod mtDNAs .
The atp8 and a atp6 of O. lunifer are the only PCGs having a seven nucleotides overlap (Figure 1). This feature is common to all lepidopteran mtDNA genomes known [5–10] and is found in many animal mtDNAs .
Total number of non-stop codons (CDs) used by the 12 analyzed mtDNAs is very similar ranging from 3695 of C. raphaelis to 3732 of O. lunifer. The codon families exhibit a very similar behavior among considered species. The eight codon families with at least 50 CDs per thousand CDs (Leu2, Ile, Phe, Met, Asn, Ser2, Gly, Tyr) encompass an average 65.82% ± 1.20% of all CDs. The three families with at least 100 CDs per thousand CDs (Leu2, Ile, Phe) account for an average 35.36% ± 0.98% of all CDs (Figure 3). The A+T rich CDs are favored over synonymous CDs with lower A+T content as proved by RSCU results (Figure 4). This point is well exemplified by the Leu2 family where the TTA codon accounts for the large majority of CDs in the family (see below). Invertebrate mitochondrial code includes 62 amino-acid encoding codons . Among the 12 analyzed genomes the total number of used codons results to be directly linked to the A+T content. The C. raphaelis mtDNA, having the highest A+T% content (see Figure 2) uses 52 codons, and never utilized the 10 G+C rich codons listed in Figure 2. Conversely, O. lunifer mtDNA, characterized by the lowest A+T% among considered lepidopteran genomes, uses all 62 codons. Differences in the number of used CDs are present between species of the same genus (e.g. B. mandarina vs. B. mori) even if the discrepancies appear circumscribed to G+C rich CDs with very limited use (e.g. GCG and CGC). The Leu1 (average = 11.73 ± 3. 82%) and Leu2 (average = 88.44 ± 3.89%) codon families are very differently represented in lepidopteran PCGs while Ser1 (average = 34.95 ± 3.67%) and Ser2 (average = 64.05 ± 1.09%) exhibit a more balanced composition.
Four amino acid residues (Leu, Ile, Phe and Ser) account for more than 44.50% (average = 45.68 ± 0.58%) of all residues forming the 13 mitochondrial proteins. The Leu and Ile amino acids share hydrophobic lateral chains, Phe is also hydrophobic and Ser exhibits an aliphatic behavior  thus their massive presence is striking but not surprising for membrane proteins.
Transfer and ribosomal RNA genes
The trnA, trnD, trnG, trnK, trnL1, trnL2, trnQ, and trnS2 of Ochrogaster mtDNA show mismatches in their stems. Mismatches are located mostly in the acceptor and anticodon stems with a single exception represented by trnD that exhibits the mismatch on the TΨC stem. Mismatches on tRNA stems are known also for the trnA, trnL1, trnL2, and trnQ, of C. raphaelis . Mismatches observed in tRNAs are corrected through RNA-editing mechanisms that are well known for arthropod mtDNA [e.g. ].
Preliminary analysis performed on rrnL and rrnS of O. lunifer revealed that these genes are capable of folding into structures (data not shown) similar to those already produced for lepidopteran mitochondrial ribosomal subunits [8, 25, 26]. Further studies, that extend the taxon sampling, are currently in progress in our lab to better define rrnL and rrnS structures within the Thaumetopoeinae subfamily that includes also O. lunifer.
Non coding regions
The s1 spacer, located between trnQ and nad2, appears to be the result of a duplicated segment (Figure 7). The s1 spacer is present in all 12 lepidopteran mtDNAs so far sequenced while it is absent in other insects . While the genomic location is constant the sequence divergence is high among species . Further investigation with a broad taxon sampling within the Lepidoptera is necessary to assess if the s1 spacer is a constant molecular signature of lepidopteran mtDNA.
The s2 spacer, placed between trnC and trnY, derives from the triplication of a six nucleotides motif with minor changes (Figure 7). An 11 bp spacer between trnC and trnY is found also in the mtDNA of A. melete and shares the ACAATT motif with the s2 spacer of O. lunifer. Because no other known lepidopteran mtDNA exhibits such a spacer its presence in A. melete and O. lunifer has to be interpreted as the result of independent events.
Spacer s3, located between nad3 and trnA, exhibits a partial duplicated segment and a poly-T motif within the first 30 nt. The second half of s3 spacer is characterized by two microsatellite repeats (CA)10(TA)12. Spacers having the same genomic location, and containing TA microsatellites are found also in B. mori and B. mandarina mtDNA genomes.
Spacer s4, inserted between trnE and trnF, contains a 5' microsatellite (TA)23, while the 3' half seems to be the triplication of a 10 nucleotides motif with some changes (Figure 7). A spacer characterized by a different motif (TATTA)31, but having the same genomic placement, is found in the A. honmai mtDNA genome.
The spacer s5, located between trnS2 and nad1, contains the ATACTAA motif which is conserved across the Lepidoptera order . This motif is possibly fundamental to site recognition by the transcription termination peptide (mtTERM protein) . Spacer s5 is present in most insect mtDNAs even if the nucleotide sequence can be quite divergent .
The s6 spacer is located between trnS2 and -rrnL and exhibits a di-nucleotide microsatellite (TA)19 directly in contact with the 3' end of rrnL gene. To date spacer s6 is known only for the mtDNA of O. lunifer.
The s7 spacer coincides with the A+T region. Several features common to the Lepidoptera A+T region  are present in the s7 spacer. The ORβ (origin of the β strand replication) is located 21 bp downstream from rrnS gene in B. mori . It contains the motif ATAGA followed by an 18 bp poly-T stretch. A very similar pattern occurs in O. lunifer where the ATAGA motif is located 17 bp downstream from rnnS gene and is followed by a 20 bp poly-T stretch (Figure 7). A microsatellite-like (AT)7(TA)3 element preceded by the ATTTA motif is present in the 3' third of O. lunifer s7 spacer. The presence of a microsatellite preceded by the ATTA motif is also a feature found in the A+T regions of other Lepidoptera . Finally a 10 bp poly-A is present immediately upstream trnM. This poly-T (in the β strand) element is still a common feature of the A+T region in Lepidoptera [8, 28]. No large repeated segments were detected in the A+T region of O. lunifer. This arrangement is consistent with other lepidopteran A+T regions while markedly contrasts with patterns observed in other insect orders [8, 29].
Intergenic spacers containing repeated elements are scattered all over the lepidopteran mtDNAs while repeated elements are restricted mostly to the A+T region in other insects . Most parts of spacers of O. lunifer are made by repeated motifs. Predominance of repeated elements suggest that mtDNA expansion can be achieved through a miss-pairing duplication mechanism, i.e. DNA slippage, during genome replication. Several intergenic spacers are restricted to a single butterfly/moth species and have not counterparts even within Lepidoptera. Thus it is plausible to suggest that spacers production occurs independently and recursively within Lepidoptera. It remains unknown while this feature is so prominent in moths and butterflies and apparently limited, reduced or absent in other insect mtDNAs sequenced to date. This behavior requires further investigation provided that mtDNA intergenic spacers are found in non-insect Arthropoda as well as other animal phyla [e.g. [18, 30]].
The mitochondrial genome of O. lunifer is the first sequenced mtDNA for a representative of the Noctuoidea a superfamily that includes about 40% of all described lepidopteran species. The newly determined genome shares the gene order, the presence of intergenic spacers, and other features with previously known lepidopteran genomes. The placement of trnM immediately after the A+T region results to be an exclusive molecular signature of all lepidopteran mtDNAs sequenced to date. Further genome sequencing will establish if this feature characterizes the whole order Lepidoptera. The mtDNA of O. lunifer exhibits a peculiar low A+T content and marked C-skew. Compared to other lepidopteran genomes it is less biased in synonymous codon usage. Comparative analysis on codon usage among lepidopteran mitochondrial genomes identified atp6, cox1, cox2, cox3, cob, nad1, nad2, nad4, and nad5 as potential markers for phylogenetic and population genetic studies. Most of the genes listed above have been previously neglected for the tasks suggested here. The massive presence of repetitive elements in intergenic spacers of O. lunifer genome lead us to suggest an important role of DNA slippage as possible mechanism to produce spacers during replication.
Sample origin and DNA extraction
An ethanol-preserved larva specimen of Ochrogaster lunifer collected in Australia (Suburb of Kenmore, Queensland, 25th February 2005) by Myron P. Zalucki (University of Queensland) was used as starting material for this study. Total DNA was extracted by applying a salting-out protocol . Quality of DNA was assessed through electrophoresis in a 1% agarose gel and staining with ethidium bromide.
PCR amplification and sequencing of Ochogaster lunifer mtDNA
PCR amplification was performed using a mix of insect universal primers [32, 33] and primers specifically designed on the O. lunifer sequences. For a full list of successful primers as well as PCR conditions see Additional file 1. The PCR products were visualized in electrophoresis in a 1% agarose gel and staining with ethidium bromide. Each PCR product represented by a single electrophoretic band was purified with the ExoSAP-IT kit (Amersham Biosciences) and directly sequenced. Sequencing of both strands was performed at the BMR Genomics service (Padova, Italy) on automated DNA sequencers mostly employing the primers used for PCR amplification.
Sequence assembly and annotation
The mtDNA final consensus sequence was assembled using the SeqMan II program from the Lasergene software package (DNAStar, Madison, WI). Genes and strands nomenclature used in this paper follows Negrisolo et al. .
Sequence analysis was performed as follows. Initially the mtDNA sequence was translated into putative proteins using the Transeq program available at the EBI web site. The true identity of these polypeptides was established using the BLAST program [34, 35] available at the NCBI web site. Gene boundaries were determined as follows. The 5' ends of PEGs were inferred to be at the first legitimate in-frame start codon (ATN, GTG, TTG, GTT; ) in the open reading frame (ORF) that was not located within the upstream gene encoded on the same strand. The only exception was atp6, which has been previously demonstrated to overlap with its upstream gene atp8 in many mtDNAs . The PCG terminus was inferred to be at the first in-frame stop codon encountered. When the stop codon was located within the sequence of a downstream gene encoded on the same strand, a truncated stop codon (T or TA) adjacent to the beginning of the downstream gene was designated as the termination codon. This codon was thought to be completed by polyadenylation to a complete TAA stop codon after transcript processing. Finally pair-wise comparisons with orthologous proteins were performed with ClustalW program  to better define the limits of PCGs.
The transfer RNA genes were identified using the tRNAscan-SE program  or recognized manually as sequences having the appropriate anticodon and capable of folding into the typical cloverleaf secondary structure .
The boundaries of the ribosomal rrnL gene were assumed to be delimited by the ends of the trnV-s6 pair. The 3' end of rrnS gene was assumed to be delimited by the start of trnV while the 5'end was determined through comparison with orthologous genes of other Lepidoptera so far sequenced.
Nucleotide composition was calculated with the EditSeq program included in the Lasergene software package. The GC-skew = (G-C)/(G+C) and AT-skew = (A-T)/(A+T) were used  to measure the base compositional difference between the different strands or between genes coded on the alternative strands. The Relative Synonymous Codon Usage (RSCU) values were calculated with MEGA 4 program .
The codon usage by analyzed genomes was investigated by calculating the two indices ENC (Effective Number of Codon used)  and MILC (Measure Independent of Length and Composition . ENC and MILC values were calculated with the INCA 2.1 program .
- atp6 and atp8:
ATP synthase subunits 6 and 8
- cob :
cytochrome c oxidase subunits 1–3
- nad1-6 and nad4L:
NADH dehydrogenase subunits 1–6 and 4L
- rrnS and rrnL:
small and large subunit ribosomal RNA (rRNA) genes
- trnX :
transfer RNA (tRNA) genes, where X is the one-letter abbreviation of the corresponding amino acid
mitochondrial genomic spacers
- A+T region:
the putative control region
protein coding gene
Relative Synonymous Codon Usage
- ENC, MILC:
Measure Independent of Length and Composition
We express our sincere thanks to Myron P. Zalucki (School of Integrative Biology, University of Queensland, Brisbane, Australia) who kindly provided the specimen of Ochrogaster lunifer used in present study. We thank Filippo Calore (Albignasego, Padova, Italy) who painted the icon of O. lunifer included in Figure 1, using as template a picture publically available at the CSIRO web site. Finally we thank two anonymous referees that provided very useful suggestions that helped to improve the manuscript.
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780.PubMedView ArticleGoogle Scholar
- Taanman JW: The mitochondrial genome: structure, transcription, translation and replication. Biochim Biophys Acta. 1999, 1410: 103-123.PubMedView ArticleGoogle Scholar
- Resh VH, Cardé RG: Insecta, Overview. Encyclopedia of Insects. Edited by: Resh VH, CArdé RG. 2003, Academic Press, Burlington MA, USA, 564-566. 1266pp.Google Scholar
- Powell JA: Lepidoptera (Moths, Butterflies). Encyclopedia of Insects. Edited by: Resh VH, CArdé RG. 2003, Academic Press, Burlington MA, USA, 631-663. 1266pp.Google Scholar
- Yukuhiro K, Sezutsu H, Itoh M, Shimizu K, Banno Y: Significant levels of sequence divergence and gene rearrangements have occurred between the mitochondrial genomes of the wild mulberry silkmoth, Bombyx mandarina and its close relative, the domesticated silkmoth, Bombyx mori. Mol Biol Evol. 2002, 19: 1385-1389.PubMedView ArticleGoogle Scholar
- Kim I, Lee EM, Seol KY, Yun EY, Lee YB, Hwang JS, Jin BR: The mitochondrial genome of the Korean hairstreak, Coreana raphaelis (Lepidoptera: Lycaenidae). Insect Mol Biol. 2006, 15 (2): 217-225.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-57.PubMedView 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: 112-123.PubMedView ArticleGoogle Scholar
- Hong MY, Lee EM, Jo YH, Park HC, Kim SR, Hwang JS, Jin BR, Kang PD, Kim KG, Han YS, Kim I: Complete nucleotide sequence and organization of the mitogenome of the silk moth Caligula boisduvalii (Lepidoptera: Saturniidae) and comparison with other lepidopteran insects. Gene. 2008, 413: 49-57.PubMedView ArticleGoogle Scholar
- Coates BS, Sumerford DV, Hellmich RL, Lewis LC: Partial mitochondrial genome sequences of Ostrinia nubilalis and Ostrinia furnicalis. Int J Biol Sci. 2005, 1: 13-18.PubMedView ArticleGoogle Scholar
- Taylor MFJ, McKechnie SW, Pierce N, Kreitman M: The lepidopteran mitochondrial control region: structure and evolution. Mol Biol Evol. 1993, 10: 1259-1272.PubMedGoogle Scholar
- Perna NT, Kocher TD: Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995, 41: 353-358.PubMedView ArticleGoogle Scholar
- Cameron SL, Whiting MF: Mitochondrial genomic comparisons of the subterranean termites from the genus Reticulitermes (Insecta: Isoptera: Rhinotermitidae). Genome. 2007, 50: 188-202.PubMedView ArticleGoogle Scholar
- Thao ML, Baumann L, Baumann P: Organization of the mitochondrial genomes of whiteflies, aphids, and psyllids (Hemiptera, Sternorrhyncha). BMC Evolutionary Biology. 2004, 4: 25-PubMedView ArticleGoogle Scholar
- Covacin C, Shao R, Cameron S, Barker SC: Extraordinary number of gene rearrangements in the mitochondrial genomes of lice (Phthiraptera: Insecta). Insect Mol Biol. 2006, 15: 63-68.PubMedView ArticleGoogle Scholar
- Shao R, Campbell NJ, Barker SC: Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera). Mol Biol Evol. 2001, 18: 858-865.PubMedView ArticleGoogle Scholar
- Wolstenholme DR: Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992, 141: 173-216.PubMedView ArticleGoogle Scholar
- Negrisolo E, Minelli A, Valle G: Extensive gene order rearrangement in the mitochondrial genome of the centipede Scutigera coleoptrata. J Mol Evol. 2004, 58: 413-423.PubMedView ArticleGoogle Scholar
- Sharp PM, Tuohy TMF, Mosurski KR: Codon usage in yeast: Cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Research. 1986, 14: 5125-5143.PubMedView ArticleGoogle Scholar
- Patthy L: Protein evolution. 2008, Blackwell, London, 374-2Google Scholar
- Wright F: The 'effective number of codons' used in a gene. Gene. 1990, 87: 23-29.PubMedView ArticleGoogle Scholar
- Supek F, Vlahovičekl K: Comparison of codon usage measures and their applicability in prediction of microbial gene expressivity. BMC Bioinformatics. 2005, 6: 182-PubMedView ArticleGoogle Scholar
- Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Hebert PDN: DNA barcodes distinguish species of tropical Lepidoptera. Proc Natl Acad Sci USA. 2006, 103: 968-971.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 USA. 2000, 97: 13738-13742.PubMedView ArticleGoogle Scholar
- Niehuis O, Naumann CM, Misof B: Identification of evolutionary conserved structural elements in the mt SSU Rrna of Zygaenoidea (Lepidoptera): a comparative sequence analysis. Org Divers Evol. 2006, 6: 17-32.View ArticleGoogle Scholar
- Niehuis O, Yen S-H, Naumann CM, Misof B: Higher phylogeny of zygaenid moths (Insecta: Lepidoptera) inferred from nuclear and mitochondrial sequence data and the evolution of larval cuticular cavities for chemical defence. Mol Phylogenet Evol. 2006, 39: 812-829.PubMedView ArticleGoogle Scholar
- Saito S, Tamura K, Aotsuka T: Replication origin of mitochondrial DNA in insects. Genetics. 2005, 171: 1695-1705.PubMedView ArticleGoogle Scholar
- Vila M, Björklund M: The utility of the neglected mitochondrial controlregion for evolutionary studies in Lepidoptera (Insecta). J Mol Evol. 2004, 58: 280-290.PubMedView ArticleGoogle Scholar
- Zhang DX, Hweitt GM: Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997, 25: 99-120.View ArticleGoogle Scholar
- Boore JL: The complete sequence of the mitochondrial genome of Nautilus macromphalus (Mollusca: Cephalopoda). BMC Genomics. 2006, 7: 182-PubMedView ArticleGoogle Scholar
- Patwary MU, Kenchington EL, Bird CJ, Zouros E: The use of random amplified polymorphic DNA markers in genetic studies of the sea scallop Placopecten magellanicus (Gmelin, 1791). J Shellfish Res. 1994, 13: 547-553.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: 651-704.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 Review of Ecology, Evolution, and Systematics. 2006, 37: 545-579.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Tatusova TA, Madden TL: BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174 (2): 247-250.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: Clustal-W – improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.PubMedView ArticleGoogle Scholar
- Smith AE, Marcker KA: N-formylmethionyl transfer RNA in mitochondria from yeast and rat liver. J Mol Biol. 1968, 38: 241-243.PubMedView ArticleGoogle Scholar
- Fearnley IM, Walker JE: Initiation codons in mammalian mitochondria: differences in genetic code in the organelle. Biochemistry. 1987, 26: 8247-8251.PubMedView ArticleGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.PubMedView ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol. 2007, 24: 1596-1599.PubMedView ArticleGoogle Scholar
- Supek F, Vlahovičekl K: INCA: synonymous codon usage analysis and clustering by means of self-organizing map. Bioinformatics. 2004, 20: 2329-2330.PubMedView ArticleGoogle Scholar
- Flook PK, Rowell CH, Gellissen G: The sequence, organization, and evolution of the Locusta migratoria mitochondrial genome. J Mol Evol. 1995, 41: 928-941.PubMedView ArticleGoogle Scholar
- Fenn JD, Cameron SL, Whiting MF: The complete mitochondrial genome sequence of the Mormon cricket (Anabrus simplex: Tettigoniidae: Orthoptera) and an analysis of control region variability. Insect Mol Biol. 2007, 16: 239-252.PubMedView ArticleGoogle Scholar
- Zhou Z, Huang Y, Shi F: The mitochondrial genome of Ruspolia dubia (Orthoptera: Conocephalidae) contains a short A+T-rich region of 70 bp in length. Genome. 2007, 50: 855-866.PubMedView ArticleGoogle Scholar
- Kim I, Cha SY, Yoon MH, Hwang JS, Lee SM, Sohn HD, Jin BR: The complete nucleotide sequence and gene organization of the mitochondrial genome of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae). Gene. 2005, 353: 155-168.PubMedView ArticleGoogle Scholar
- Cameron SL, Barker SC, Whiting MF: Mitochondrial genomics and the new insect order Mantophasmatodea. Mol Phylogenet Evol. 2006, 38 (1): 274-279.PubMedView ArticleGoogle Scholar
- Yamauchi MM, Miya MU, Nishida M: Use of a PCR-based approach for sequencing whole mitochondrial genomes of insects: two examples (cockroach and dragonfly) based on the method developed for decapod crustaceans. Insect Mol Biol. 2004, 13: 435-442.PubMedView ArticleGoogle Scholar
- Stewart JB, Beckenbach AT: Insect mitochondrial genomics 2: The complete mitochondrial genome sequence of a giant stonefly, Pteronarcys princeps, asymmetric directional mutation bias, and conserved plecopteran A+T-region elements. Genome. 2006, 49: 815-824.PubMedView ArticleGoogle Scholar
- Cameron SL, Johnson KP, Whiting MF: The mitochondrial genome of the screamer louse Bothriometopus (phthiraptera: ischnocera): effects of extensive gene rearrangements on the evolution of the genome. J Mol Evol. 2007, 65: 589-604.PubMedView ArticleGoogle Scholar
- Stewart JB, Beckenbach AT: Insect mitochondrial genomics: the complete mitochondrial genome sequence of the meadow spittlebug Philaenus spumarius (Hemiptera: Auchenorrhyncha: Cercopoidae). Genome. 2005, 48: 46-54.PubMedView ArticleGoogle Scholar
- Dotson EM, Beard CB: Sequence and organization of the mitochondrial genome of the Chagas disease vector, Triatoma dimidiata. Insect Mol Biol. 2001, 10: 205-215.PubMedView ArticleGoogle Scholar
- Shao R, Dowton M, Murrell A, Barker SC: Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of Insects. Mol Biol Evol. 2003, 20: 1612-1619.PubMedView ArticleGoogle Scholar
- Shao R, Barker SC: The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta: Thysanoptera): convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes. Mol Biol Evol. 2003, 20: 362-370.PubMedView ArticleGoogle Scholar
- Stewart JB, Beckenbach AT: Phylogenetic and genomic analysis of the complete mitochondrial DNA sequence of the spotted asparagus beetle Crioceris duodecimpunctata. Mol Phylogenet Evol. 2003, 26 (3): 513-526.PubMedView ArticleGoogle Scholar
- Arnoldi FG, Ogoh K, Ohmiya Y, Viviani VR: Mitochondrial genome sequence of the Brazilian luminescent click beetle Pyrophorus divergens (Coleoptera: Elateridae): mitochondrial genes utility to investigate the evolutionary history of Coleoptera and its bioluminescence. Gene. 2007, 405: 1-9.PubMedView ArticleGoogle Scholar
- Bae JS, Kim I, Sohn HD, Jin BR: The mitochondrial genome of the firefly, Pyrocoelia rufa: complete DNA sequence, genome organization, and phylogenetic analysis with other insects. Mol Phylogenet Evol. 2004, 32 (3): 978-985.PubMedView ArticleGoogle Scholar
- Friedrich M, Muqim N: Sequence and phylogenetic analysis of the complete mitochondrial genome of the flour beetle Tribolium castanaeum. Mol Phylogenet Evol. 2003, 26 (3): 502-512.PubMedView ArticleGoogle Scholar
- Beard CB, Hamm DM, Collins FH: The mitochondrial genome of the mosquito Anopheles gambiae: DNA sequence, genome organization, and comparisons with mitochondrial sequences of other insects. Insect Mol Biol. 1993, 2: 103-124.PubMedView ArticleGoogle Scholar
- Mitchell SE, Cockburn AF, Seawright JA: The mitochondrial genome of Anopheles quadrimaculatus species A: complete nucleotide sequence and gene organization. Genome. 1993, 36: 1058-1073.PubMedView ArticleGoogle Scholar
- Lessinger AC, Martins Junqueira AC, Lemos TA, Kemper EL, da Silva FR, Vettore AL, Arruda P, Azeredo-Espin AM: The mitochondrial genome of the primary screwworm fly Cochliomyia hominivorax (Diptera: Calliphoridae). Insect Mol Biol. 2000, 9: 521-529.PubMedView ArticleGoogle Scholar
- Junqueira AC, Lessinger AC, Torres TT, Da Silva FR, Vettore AL, Arruda P, Azeredo Espin AM: The mitochondrial genome of the blowfly Chrysomya chloropyga (Diptera: Calliphoridae). Gene. 2004, 339: 7-15.PubMedView ArticleGoogle Scholar
- Lewis DL, Farr CL, Kaguni LS: Drosophila melanogaster mitochondrial DNA: completion of the nucleotide sequence and evolutionary comparisons. Insect Mol Biol. 1995, 4: 263-278.PubMedView ArticleGoogle Scholar
- Ballard JW: Comparative genomics of mitochondrial DNA in members of the Drosophila melanogaster subgroup. J Mol Evol. 2000, 51: 48-63.PubMedGoogle Scholar
- 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.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: 40-59.View ArticleGoogle Scholar
- Spanos L, Koutroumbas G, Kotsyfakis M, Louis C: The mitochondrial genome of the mediterranean fruit fly, Ceratitis capitata. Insect Mol Biol. 2000, 9: 139-144.PubMedView ArticleGoogle Scholar
- Nardi F, Carapelli A, Dallai R, Frati F: The mitochondrial genome of the olive fly Bactrocera oleae: two haplotypes from distant geographical locations. Insect Mol Biol. 2003, 12: 605-611.PubMedView ArticleGoogle Scholar
- Crozier RH, Crozier YC: The mitochondrial genome of the honeybee Apis mellifera: complete sequence and genome organization. Genetics. 1993, 133: 97-117.PubMedGoogle Scholar
- Cha SY, Yoon HJ, Lee EM, Yoon MH, Hwang JS, Jin BR, Han YS, Kim I: The complete nucleotide sequence and gene organization of the mitochondrial genome of the bumblebee, Bombus ignitus (Hymenoptera: Apidae). Gene. 2007, 392: 206-220.PubMedView ArticleGoogle Scholar
- Castro LR, Ruberu K, Dowton M: Mitochondrial genomes of Vanhornia eucnemidarum (Apocrita: Vanhorniidae) and Primeuchroeus spp. (Aculeata: Chrysididae): Evidence of rearranged mitochondrial genomes within the Apocrita (Insecta: Hymenoptera). Genome. 2006, 49: 752-766.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.