<title>Background</title> <p>Animal mitochondrial genomes (mtDNAs) are usually circular molecules spanning 16–20 kbp that contain 13 protein-coding genes (PCGs), 2 ribosomal RNA and 22 transfer (tRNA) genes [##REF##10101183##1##]. Non-coding control elements, that regulate the transcription and replication of the genome, are also present in mtDNAs [##REF##10101183##1##,##REF##10076021##2##]. Mitochondrial genomes are very important subject for different scientific disciplines including animal health, comparative and evolutionary genomics, molecular evolution, phylogenetics and population genetics. However, current knowledge on mtDNAs is very uneven as well exemplified by sequences available in GenBank that were obtained mostly from vertebrate taxa. Insects constitute the most species-rich class among animals with almost a million of taxa described to date [##UREF##0##3##]. Within the insects, the Lepidoptera (butterflies plus moths) order accounts for more than 160,000 species [##UREF##1##4##]. Despite this huge taxonomic diversity the existing information on lepidopteran mtDNA is very limited. Complete sequences have been determined for the two butterflies <italic>Coreana raphaelis </italic>and <italic>Artogeia melete</italic>, and for the seven moths <italic>Adoxophyes honmai, Antheraea pernyi, Bombyx mori, Bombyx mandarina </italic>and <italic>Manduca sexta, Phthonandria atrilineata </italic>and <italic>Saturnia boisduvalii </italic>[##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9####18337026##9##] while near complete sequences exist for <italic>Ostrinia furnacalis </italic>and <italic>Ostrinia nubilalis </italic>[##REF##15951845##10##] (Table ##TAB##0##1##). Current genomic knowledge of Lepidoptera is very scanty and the covered taxon-sampling is extremely poor and limited to six superfamiles among the 45–48 known, and to 9 families of the recognized 120 [##UREF##1##4##]. A better understanding of the lepidopteran mtDNA requires an expansion of taxon and genome samplings. We were able to fully sequence the mitochondrial genome of the bag-shelter moth <italic>Ochrogaster lunifer</italic>. The newly determined mtDNA is the first complete sequence for the Superfamily Noctuoidea, a very large assemblage that accounts for about 40% of all described Lepidoptera. [##UREF##1##4##]. In the present paper the <italic>Ochrogaster </italic>genome is described and compared with mtDNAs of other lepidopterans as well as pterygote Insecta.</p>
<title>Methods</title> <title>Sample origin and DNA extraction</title> <p>An ethanol-preserved larva specimen of <italic>Ochrogaster lunifer </italic>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 [##UREF##5##31##]. Quality of DNA was assessed through electrophoresis in a 1% agarose gel and staining with ethidium bromide.</p> <title>PCR amplification and sequencing of <italic>Ochogaster lunifer </italic>mtDNA</title> <p>PCR amplification was performed using a mix of insect universal primers [##UREF##6##32##,##UREF##7##33##] and primers specifically designed on the <italic>O. lunifer </italic>sequences. For a full list of successful primers as well as PCR conditions see Additional file ##SUPPL##0##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.</p> <title>Sequence assembly and annotation</title> <p>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. [##REF##15114420##18##].</p> <p>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 [##REF##2231712##34##,##REF##10339815##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; [##REF##1452431##17##]) in the open reading frame (ORF) that was not located within the upstream gene encoded on the same strand. The only exception was <italic>atp6</italic>, which has been previously demonstrated to overlap with its upstream gene <italic>atp8 </italic>in many mtDNAs [##REF##1452431##17##]. 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 [##REF##7984417##36##] to better define the limits of PCGs.</p> <p>Irrespectively of the real initiation codon, a formyl-Met was assumed to be the starting amino acid for all the proteins as previously proved for other mitochondrial genomes [##REF##5760639##37##,##REF##2964865##38##].</p> <p>The transfer RNA genes were identified using the tRNAscan-SE program [##REF##9023104##39##] or recognized manually as sequences having the appropriate anticodon and capable of folding into the typical cloverleaf secondary structure [##REF##1452431##17##].</p> <p>The boundaries of the ribosomal <italic>rrnL </italic>gene were assumed to be delimited by the ends of the <italic>trnV</italic>-s6 pair. The 3' end of <italic>rrnS </italic>gene was assumed to be delimited by the start of <italic>trnV </italic>while the 5'end was determined through comparison with orthologous genes of other Lepidoptera so far sequenced.</p> <title>Genomic analysis</title> <p>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 [##REF##7563121##12##] 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 [##REF##17488738##40##].</p> <p>The codon usage by analyzed genomes was investigated by calculating the two indices ENC (Effective Number of Codon used) [##REF##2110097##21##] and MILC (Measure Independent of Length and Composition [##REF##16029499##22##]. ENC and MILC values were calculated with the INCA 2.1 program [##REF##15059815##41##].</p>
<title>Results and discussion</title> <title>Genome organization, structure and composition</title> <p>The mtDNA genome of <italic>O. lunifer </italic>is a circular molecule 15593 bp long. It includes the entire set of 37 genes usually present in animal mtDNAs [##REF##10101183##1##], i.e., 13 PCGs, 22 tRNA genes, and 2 ribosomal genes (Figure ##FIG##0##1##). The mtDNA genome of <italic>O. lunifer </italic>contains also 7 intergenic spacers (s1–s7), spanning at least 15 bp, described in a paragraph below. Genes on the same strand are overlapped (e.g. <italic>trnM </italic>vs. <italic>trnI</italic>; <italic>atp8 </italic>vs. <italic>atp6</italic>), contiguous, separated by few nucleotides or by intergenic spacers (e.g. <italic>nad3 </italic>vs. <italic>trnA</italic>; <italic>trnC </italic>vs. <italic>trnY</italic>). Genes on opposite strands exhibit a similar behavior (Figure ##FIG##0##1##).</p> <p>The <italic>O. lunifer </italic>mtDNA has the typical lepidopteran gene order [##REF##18065166##8##,##REF##18337026##9##] that differs from the ancestral gene order of insects [##REF##10101183##1##] for the placement of <italic>trnM</italic>. In the ancestral type (e.g. <italic>Drosophila yakuba </italic>mtDNA) the order in the α strand is: A+T region, <italic>trnI</italic>, <italic>trnQ</italic>, <italic>trnM</italic>, <italic>nad2</italic>. In all lepidopteran mtDNAs, sequenced to date, the order is: A+T region, <italic>trnM</italic>, <italic>trnI</italic>, <italic>trnQ</italic>, <italic>nad2 </italic>which implies the translocation of <italic>trnM </italic>[##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9##, ##REF##15951845##10##, ##REF##8277854##11####8277854##11##]. This placement of <italic>trnM </italic>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.</p> <p>The composition of the α strand of <italic>O. lunifer </italic>mtDNA is A = 6252 (40.09%), T = 5886 (37.75%), G = 1179 (7.56%) and C = 2276 (14.60%).</p> <p>The A+T% and G+C% values for the α strand as well as the A- and G-skews [##REF##7563121##12##] were calculated for all available complete mtDNA genomes of Pterygota and are presented in the scatter plots of Figure ##FIG##1##2##.</p> <p>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 (<italic>Apis mellifera</italic>, <italic>Bombus ignitus </italic>and <italic>Melipona bicolor</italic>) and two bugs (<italic>Aleurodicus dugesii </italic>and <italic>Schizaphis graminum</italic>). All lepidopteran mtDNAs but <italic>O. lunifer </italic>exhibit high A+T% values. The A+T content of <italic>O. lunifer </italic>mtDNA is 77.84% that represents the lowest value for lepidopteran complete mtDNAs [##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8####18065166##8##,##REF##15951845##10##]. The lowest A+T contents are found in the termite mtDNAs (<italic>Reticulitermes </italic>spp.). Extreme A+T values are also shared by species having highly re-arranged gene order [##REF##17546084##13##]. However the possession of a re-arranged genome is not sufficient <italic>per se </italic>to have an A+T content drastically departing from the average (e.g. <italic>Aleurochiton aceris </italic>and <italic>Bemisia tabaci</italic>). The A+T values appear to be linked to taxonomic relatedness at low rank (i.e. genus, family) (e.g. species of <italic>Drosophila</italic>, species of <italic>Bactrocera</italic>, members of family <italic>Apidae</italic>). 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 (<italic>A. dugesii </italic>vs.<italic>Triatoma dimidiata)</italic>.</p> <p>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 (<italic>B. ignitus</italic>) to 0.18247 (<italic>Locusta migratoria</italic>). The lepidopteran A-skews vary from -0.04748 (<italic>C. raphaelis</italic>) to 0.05872 (<italic>B. mori</italic>) with the <italic>O. lunifer </italic>mtDNA exhibiting a slight A-skew (0.03015). The <italic>Reticulitermes </italic>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 <italic>Campanulotes bidentatus </italic>and <italic>Trialeuroides vaporarium </italic>that have low A+T% content and gene-orders different than insect ancestral gene order [##REF##10101183##1##,##REF##15291971##14##,##REF##16469069##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 <italic>Heterodoxus macropus </italic>[##REF##11319269##16##].</p> <p>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 (<italic>Vanhornia eucnemidarum</italic>) to -0.01250 (<italic>Heterodoxus macropus</italic>). The main exception is represented by the mtDNA of bugs, while the highest G-skewed genome is that of <italic>C. bidentatus</italic>. 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 <italic>O. lunifer </italic>that exhibits the second most pronounced C-skew (G-skew = -0.31751) among analyzed genomes.</p> <p>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 <italic>Reticulitermes santonensis </italic>and <italic>Reticulitermes virginicus </italic>mtDNAs. The same reasoning applies at high taxonomic rank to the Hemiptera. The mtDNA of <italic>C. bidentatus </italic>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 <italic>Reticulitermes </italic>that exhibit very strong A-skews but not G-skews.</p> <p>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 <italic>O. lunifer </italic>mtDNA that exhibits a A+T percentage different than other lepidopteran mtDNAs that share high A+T contents [##REF##18065166##8##,##REF##18337026##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.</p> <title>Protein-coding genes</title> <p>The mtDNA of <italic>O. lunifer </italic>contains the full set of PCGs usually present in animal mtDNA. PCGs are arranged along the genome according to the standard order of Insects [##REF##10101183##1##] (Figure ##FIG##0##1##). The putative start codons of PCGs are those previously known for animal mtDNA i.e. ATN, GTG, TTG, GTT [##REF##1452431##17##] with the only exception represented by the CGA start codon of <italic>cox1 </italic>gene. This non-canonical putative start codon is found also in the butterfly <italic>A. melete </italic>and in the moths <italic>A. honmai, B. mori, B. mandarina</italic>, <italic>M. sexta </italic>and <italic>P. atrilineata </italic>[##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8####18065166##8##]. In the butterfly <italic>C. raphaelis </italic>the tetranucleotide TTAG is the putative start codon [##REF##16640732##6##] and the six nucleotide TATTAG has been suggested as putative start codon for the moths <italic>O. nubilalis </italic>and <italic>O. furnicalis </italic>[##REF##15951845##10##]. An unusual start codon for <italic>cox1 </italic>gene is known in various arthropod mtDNA [e.g. [##REF##15114420##18##]].</p> <p>The <italic>cox1</italic>, <italic>cox2</italic>, <italic>nad5</italic>, and <italic>nad4 </italic>genes of <italic>O. lunifer </italic>mtDNA have incomplete stop codons. The presence of incomplete stop codons is a feature shared with all lepidopteran mtDNAs sequenced to date [##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9##, ##REF##15951845##10####15951845##10##] and more in general with many arthropod mtDNAs [##REF##10101183##1##].</p> <p>The <italic>atp8 </italic>and a <italic>atp6 </italic>of <italic>O. lunifer </italic>are the only PCGs having a seven nucleotides overlap (Figure ##FIG##0##1##). This feature is common to all lepidopteran mtDNA genomes known [##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9##, ##REF##15951845##10####15951845##10##] and is found in many animal mtDNAs [##REF##10101183##1##].</p> <p>The abundance of codon families and Relative Synonymous Codon Usage (RSCU) [##REF##3526280##19##] in PCGs were investigated for all available lepidopteran mtDNAs and the results are summarized in Figures ##FIG##2##3## and ##FIG##3##4##. All first codons as well as stop codons, complete and incomplete, were excluded from the analysis to avoid biases due to unusual putative start codons and incomplete stop codons.</p> <p>Total number of non-stop codons (CDs) used by the 12 analyzed mtDNAs is very similar ranging from 3695 of <italic>C. raphaelis </italic>to 3732 of <italic>O. lunifer</italic>. 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 ##FIG##2##3##). The A+T rich CDs are favored over synonymous CDs with lower A+T content as proved by RSCU results (Figure ##FIG##3##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 [##REF##10101183##1##]. Among the 12 analyzed genomes the total number of used codons results to be directly linked to the A+T content. The <italic>C. raphaelis </italic>mtDNA, having the highest A+T% content (see Figure ##FIG##1##2##) uses 52 codons, and never utilized the 10 G+C rich codons listed in Figure ##FIG##1##2##. Conversely, <italic>O. lunifer </italic>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. <italic>B. mandarina </italic>vs. <italic>B. mori</italic>) 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.</p> <p>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 [##UREF##2##20##] thus their massive presence is striking but not surprising for membrane proteins.</p> <p>Codon usage by single PCGs was investigated by calculating the two indices ENC (Effective number of codon used) [##REF##2110097##21##] and MILC (Measure Independent of Length and Composition) [##REF##16029499##22##]. Both indices, based on different approaches [##REF##2110097##21##,##REF##16029499##22##] provide a measure of codon variability of PCGs. The ENC and MILC estimate the codon variability in a way that allows comparison among sequences having different lengths as is the case of various PCGs. Genes exhibiting a higher diversity in codon usage have generally a higher number of variable sites, a prerequisite to be potential phylogenetic markers. Thus the use of ENC and MILC scores, according to the new approach presented in this paper, is a way to study PCG sequences variability on a codon perspective. The best scores of both indices should allow to identify the more diverse PCGs in a approach complementary to the usual method based on evolutionary distances among orthologous sequences (e.g. [##REF##18065166##8##]). The assessment of genetic variability is an interesting point. Indeed some PCGs are standard marker for species recognition [##REF##16418261##23##] or have been extensively used as phylogenetic markers in Lepidoptera while others have received so far limited or no attention. Understanding the genetic diversity of each PCG is a prerequisite to determine its phylogenetic usefulness. The ENC and MILC values were calculated for all PCGs but <italic>atp8 </italic>that contains too a few codons to get reliable ENC/MILC estimations [##REF##16029499##22##]. Calculations were extended also to all 13 PCGs pooled as well as to the pooled PCGs belonging to α and β strands respectively. The scatter plot analysis is provided in Figure ##FIG##4##5##. As expected the greatest diversity in codon usage is found when all codons are considered. Good codon diversity is found also when all PCGs of α or β strands are considered. More interesting is the behavior of single genes. In this latter case sequences well established as phylogenetic markers (i.e. <italic>cox1</italic>, <italic>cob</italic>, <italic>nad5</italic>, and <italic>cox2</italic>) are intermixed with PCGs poorly or not considered by researchers (e.g. <italic>cox3</italic>, <italic>nad4</italic>, <italic>nad1</italic>, <italic>nad2</italic>). Our results suggest that the neglected PCGs should be considered as potential markers thus extending the number of mtDNA PCGs sampled for population as well as phylogenetic markers. Findings, based on codon diversity, must be integrated with direct comparisons of sequences [##REF##18065166##8##] that allow to better define the optimal task that each gene can perform i.e. to be used at low taxonomic level or at high taxonomic level.</p> <title>Transfer and ribosomal RNA genes</title> <p><italic>Ochrogaster </italic>genome has the characteristic 22 tRNAs set (Figure ##FIG##5##6##) present in most of animal mtDNAs [##REF##10101183##1##]. All tRNAs present the typical clover leaf secondary structure but <italic>trnS1 </italic>lacks the DHU stem. This feature is shared with the <italic>C. raphaelis </italic>mtDNA [##REF##16640732##6##] but is not a general feature of lepidopteran mtDNA as proved by <italic>A. honmai </italic>that has all tRNAs with a complete clover leaf structure [##REF##16488090##7##]. In general, the lack of DHU arm in <italic>trnS1 </italic>is a common condition in metazoan mtDNAs [##REF##11095730##24##].</p> <p>The <italic>trnA</italic>, <italic>trnD</italic>, <italic>trnG</italic>, <italic>trnK</italic>, <italic>trnL1</italic>, <italic>trnL2</italic>, <italic>trnQ</italic>, and <italic>trnS2 </italic>of <italic>Ochrogaster </italic>mtDNA show mismatches in their stems. Mismatches are located mostly in the acceptor and anticodon stems with a single exception represented by <italic>trnD </italic>that exhibits the mismatch on the TΨC stem. Mismatches on tRNA stems are known also for the <italic>trnA</italic>, <italic>trnL1</italic>, <italic>trnL2</italic>, and <italic>trnQ</italic>, of <italic>C. raphaelis </italic>[##REF##16640732##6##]. Mismatches observed in tRNAs are corrected through RNA-editing mechanisms that are well known for arthropod mtDNA [e.g. [##REF##11095730##24##]].</p> <p>Preliminary analysis performed on <italic>rrnL </italic>and <italic>rrnS </italic>of <italic>O. lunifer </italic>revealed that these genes are capable of folding into structures (data not shown) similar to those already produced for lepidopteran mitochondrial ribosomal subunits [##REF##18065166##8##,##UREF##3##25##,##REF##16483803##26##]. Further studies, that extend the taxon sampling, are currently in progress in our lab to better define <italic>rrnL </italic>and <italic>rrnS </italic>structures within the Thaumetopoeinae subfamily that includes also <italic>O. lunifer</italic>.</p> <title>Non coding regions</title> <p>The mtDNA genome of <italic>O. lunifer </italic>contains 7 intergenic spacers (s1–s7) spanning at least 15 bp (Figures ##FIG##0##1## and ##FIG##6##7##). The features of s1–s7 spacers are presented below with reference to the α strand for orientation and sequence motifs description.</p> <p>The s1 spacer, located between <italic>trnQ </italic>and <italic>nad2</italic>, appears to be the result of a duplicated segment (Figure ##FIG##6##7##). The s1 spacer is present in all 12 lepidopteran mtDNAs so far sequenced while it is absent in other insects [##REF##18065166##8##]. While the genomic location is constant the sequence divergence is high among species [##REF##18065166##8##]. 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.</p> <p>The s2 spacer, placed between <italic>trnC </italic>and <italic>trnY</italic>, derives from the triplication of a six nucleotides motif with minor changes (Figure ##FIG##6##7##). An 11 bp spacer between <italic>trnC </italic>and <italic>trnY </italic>is found also in the mtDNA of <italic>A. melete </italic>and shares the ACAATT motif with the s2 spacer of <italic>O. lunifer</italic>. Because no other known lepidopteran mtDNA exhibits such a spacer its presence in <italic>A. melete </italic>and <italic>O. lunifer </italic>has to be interpreted as the result of independent events.</p> <p>Spacer s3, located between <italic>nad3 </italic>and <italic>trnA</italic>, 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)<sub>10</sub>(TA)<sub>12</sub>. Spacers having the same genomic location, and containing TA microsatellites are found also in <italic>B. mori </italic>and <italic>B. mandarina </italic>mtDNA genomes.</p> <p>Spacer s4, inserted between <italic>trnE </italic>and <italic>trnF</italic>, contains a 5' microsatellite (TA)<sub>23</sub>, while the 3' half seems to be the triplication of a 10 nucleotides motif with some changes (Figure ##FIG##6##7##). A spacer characterized by a different motif (TATTA)<sub>31</sub>, but having the same genomic placement, is found in the <italic>A. honmai </italic>mtDNA genome.</p> <p>The spacer s5, located between <italic>trnS2 </italic>and <italic>nad1</italic>, contains the ATACTAA motif which is conserved across the Lepidoptera order [##REF##18065166##8##]. This motif is possibly fundamental to site recognition by the transcription termination peptide (mtTERM protein) [##REF##10076021##2##]. Spacer s5 is present in most insect mtDNAs even if the nucleotide sequence can be quite divergent [##REF##18065166##8##].</p> <p>The s6 spacer is located between <italic>trnS2 </italic>and -<italic>rrnL </italic>and exhibits a di-nucleotide microsatellite (TA)<sub>19 </sub>directly in contact with the 3' end of <italic>rrnL </italic>gene. To date spacer s6 is known only for the mtDNA of <italic>O. lunifer</italic>.</p> <p>The s7 spacer coincides with the A+T region. Several features common to the Lepidoptera A+T region [##REF##18065166##8##] are present in the s7 spacer. The O<sub>R</sub>β (origin of the β strand replication) is located 21 bp downstream from <italic>rrnS </italic>gene in <italic>B. mori </italic>[##REF##16118189##27##]. It contains the motif ATAGA followed by an 18 bp poly-T stretch. A very similar pattern occurs in <italic>O. lunifer </italic>where the ATAGA motif is located 17 bp downstream from <italic>rnnS </italic>gene and is followed by a 20 bp poly-T stretch (Figure ##FIG##6##7##). A microsatellite-like (AT)<sub>7</sub>(TA)<sub>3 </sub>element preceded by the ATTTA motif is present in the 3' third of <italic>O. lunifer </italic>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 [##REF##18065166##8##]. Finally a 10 bp poly-A is present immediately upstream <italic>trnM</italic>. This poly-T (in the β strand) element is still a common feature of the A+T region in Lepidoptera [##REF##18065166##8##,##REF##15045483##28##]. No large repeated segments were detected in the A+T region of <italic>O. lunifer</italic>. This arrangement is consistent with other lepidopteran A+T regions while markedly contrasts with patterns observed in other insect orders [##REF##18065166##8##,##UREF##4##29##].</p> <p>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 [##REF##18065166##8##]. Most parts of spacers of <italic>O. lunifer </italic>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. [##REF##15114420##18##,##REF##16854241##30##]].</p>
<title>Results and discussion</title> <title>Genome organization, structure and composition</title> <p>The mtDNA genome of <italic>O. lunifer </italic>is a circular molecule 15593 bp long. It includes the entire set of 37 genes usually present in animal mtDNAs [##REF##10101183##1##], i.e., 13 PCGs, 22 tRNA genes, and 2 ribosomal genes (Figure ##FIG##0##1##). The mtDNA genome of <italic>O. lunifer </italic>contains also 7 intergenic spacers (s1–s7), spanning at least 15 bp, described in a paragraph below. Genes on the same strand are overlapped (e.g. <italic>trnM </italic>vs. <italic>trnI</italic>; <italic>atp8 </italic>vs. <italic>atp6</italic>), contiguous, separated by few nucleotides or by intergenic spacers (e.g. <italic>nad3 </italic>vs. <italic>trnA</italic>; <italic>trnC </italic>vs. <italic>trnY</italic>). Genes on opposite strands exhibit a similar behavior (Figure ##FIG##0##1##).</p> <p>The <italic>O. lunifer </italic>mtDNA has the typical lepidopteran gene order [##REF##18065166##8##,##REF##18337026##9##] that differs from the ancestral gene order of insects [##REF##10101183##1##] for the placement of <italic>trnM</italic>. In the ancestral type (e.g. <italic>Drosophila yakuba </italic>mtDNA) the order in the α strand is: A+T region, <italic>trnI</italic>, <italic>trnQ</italic>, <italic>trnM</italic>, <italic>nad2</italic>. In all lepidopteran mtDNAs, sequenced to date, the order is: A+T region, <italic>trnM</italic>, <italic>trnI</italic>, <italic>trnQ</italic>, <italic>nad2 </italic>which implies the translocation of <italic>trnM </italic>[##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9##, ##REF##15951845##10##, ##REF##8277854##11####8277854##11##]. This placement of <italic>trnM </italic>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.</p> <p>The composition of the α strand of <italic>O. lunifer </italic>mtDNA is A = 6252 (40.09%), T = 5886 (37.75%), G = 1179 (7.56%) and C = 2276 (14.60%).</p> <p>The A+T% and G+C% values for the α strand as well as the A- and G-skews [##REF##7563121##12##] were calculated for all available complete mtDNA genomes of Pterygota and are presented in the scatter plots of Figure ##FIG##1##2##.</p> <p>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 (<italic>Apis mellifera</italic>, <italic>Bombus ignitus </italic>and <italic>Melipona bicolor</italic>) and two bugs (<italic>Aleurodicus dugesii </italic>and <italic>Schizaphis graminum</italic>). All lepidopteran mtDNAs but <italic>O. lunifer </italic>exhibit high A+T% values. The A+T content of <italic>O. lunifer </italic>mtDNA is 77.84% that represents the lowest value for lepidopteran complete mtDNAs [##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8####18065166##8##,##REF##15951845##10##]. The lowest A+T contents are found in the termite mtDNAs (<italic>Reticulitermes </italic>spp.). Extreme A+T values are also shared by species having highly re-arranged gene order [##REF##17546084##13##]. However the possession of a re-arranged genome is not sufficient <italic>per se </italic>to have an A+T content drastically departing from the average (e.g. <italic>Aleurochiton aceris </italic>and <italic>Bemisia tabaci</italic>). The A+T values appear to be linked to taxonomic relatedness at low rank (i.e. genus, family) (e.g. species of <italic>Drosophila</italic>, species of <italic>Bactrocera</italic>, members of family <italic>Apidae</italic>). 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 (<italic>A. dugesii </italic>vs.<italic>Triatoma dimidiata)</italic>.</p> <p>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 (<italic>B. ignitus</italic>) to 0.18247 (<italic>Locusta migratoria</italic>). The lepidopteran A-skews vary from -0.04748 (<italic>C. raphaelis</italic>) to 0.05872 (<italic>B. mori</italic>) with the <italic>O. lunifer </italic>mtDNA exhibiting a slight A-skew (0.03015). The <italic>Reticulitermes </italic>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 <italic>Campanulotes bidentatus </italic>and <italic>Trialeuroides vaporarium </italic>that have low A+T% content and gene-orders different than insect ancestral gene order [##REF##10101183##1##,##REF##15291971##14##,##REF##16469069##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 <italic>Heterodoxus macropus </italic>[##REF##11319269##16##].</p> <p>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 (<italic>Vanhornia eucnemidarum</italic>) to -0.01250 (<italic>Heterodoxus macropus</italic>). The main exception is represented by the mtDNA of bugs, while the highest G-skewed genome is that of <italic>C. bidentatus</italic>. 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 <italic>O. lunifer </italic>that exhibits the second most pronounced C-skew (G-skew = -0.31751) among analyzed genomes.</p> <p>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 <italic>Reticulitermes santonensis </italic>and <italic>Reticulitermes virginicus </italic>mtDNAs. The same reasoning applies at high taxonomic rank to the Hemiptera. The mtDNA of <italic>C. bidentatus </italic>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 <italic>Reticulitermes </italic>that exhibit very strong A-skews but not G-skews.</p> <p>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 <italic>O. lunifer </italic>mtDNA that exhibits a A+T percentage different than other lepidopteran mtDNAs that share high A+T contents [##REF##18065166##8##,##REF##18337026##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.</p> <title>Protein-coding genes</title> <p>The mtDNA of <italic>O. lunifer </italic>contains the full set of PCGs usually present in animal mtDNA. PCGs are arranged along the genome according to the standard order of Insects [##REF##10101183##1##] (Figure ##FIG##0##1##). The putative start codons of PCGs are those previously known for animal mtDNA i.e. ATN, GTG, TTG, GTT [##REF##1452431##17##] with the only exception represented by the CGA start codon of <italic>cox1 </italic>gene. This non-canonical putative start codon is found also in the butterfly <italic>A. melete </italic>and in the moths <italic>A. honmai, B. mori, B. mandarina</italic>, <italic>M. sexta </italic>and <italic>P. atrilineata </italic>[##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8####18065166##8##]. In the butterfly <italic>C. raphaelis </italic>the tetranucleotide TTAG is the putative start codon [##REF##16640732##6##] and the six nucleotide TATTAG has been suggested as putative start codon for the moths <italic>O. nubilalis </italic>and <italic>O. furnicalis </italic>[##REF##15951845##10##]. An unusual start codon for <italic>cox1 </italic>gene is known in various arthropod mtDNA [e.g. [##REF##15114420##18##]].</p> <p>The <italic>cox1</italic>, <italic>cox2</italic>, <italic>nad5</italic>, and <italic>nad4 </italic>genes of <italic>O. lunifer </italic>mtDNA have incomplete stop codons. The presence of incomplete stop codons is a feature shared with all lepidopteran mtDNAs sequenced to date [##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9##, ##REF##15951845##10####15951845##10##] and more in general with many arthropod mtDNAs [##REF##10101183##1##].</p> <p>The <italic>atp8 </italic>and a <italic>atp6 </italic>of <italic>O. lunifer </italic>are the only PCGs having a seven nucleotides overlap (Figure ##FIG##0##1##). This feature is common to all lepidopteran mtDNA genomes known [##REF##12140251##5##, ####REF##16640732##6##, ##REF##16488090##7##, ##REF##18065166##8##, ##REF##18337026##9##, ##REF##15951845##10####15951845##10##] and is found in many animal mtDNAs [##REF##10101183##1##].</p> <p>The abundance of codon families and Relative Synonymous Codon Usage (RSCU) [##REF##3526280##19##] in PCGs were investigated for all available lepidopteran mtDNAs and the results are summarized in Figures ##FIG##2##3## and ##FIG##3##4##. All first codons as well as stop codons, complete and incomplete, were excluded from the analysis to avoid biases due to unusual putative start codons and incomplete stop codons.</p> <p>Total number of non-stop codons (CDs) used by the 12 analyzed mtDNAs is very similar ranging from 3695 of <italic>C. raphaelis </italic>to 3732 of <italic>O. lunifer</italic>. 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 ##FIG##2##3##). The A+T rich CDs are favored over synonymous CDs with lower A+T content as proved by RSCU results (Figure ##FIG##3##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 [##REF##10101183##1##]. Among the 12 analyzed genomes the total number of used codons results to be directly linked to the A+T content. The <italic>C. raphaelis </italic>mtDNA, having the highest A+T% content (see Figure ##FIG##1##2##) uses 52 codons, and never utilized the 10 G+C rich codons listed in Figure ##FIG##1##2##. Conversely, <italic>O. lunifer </italic>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. <italic>B. mandarina </italic>vs. <italic>B. mori</italic>) 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.</p> <p>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 [##UREF##2##20##] thus their massive presence is striking but not surprising for membrane proteins.</p> <p>Codon usage by single PCGs was investigated by calculating the two indices ENC (Effective number of codon used) [##REF##2110097##21##] and MILC (Measure Independent of Length and Composition) [##REF##16029499##22##]. Both indices, based on different approaches [##REF##2110097##21##,##REF##16029499##22##] provide a measure of codon variability of PCGs. The ENC and MILC estimate the codon variability in a way that allows comparison among sequences having different lengths as is the case of various PCGs. Genes exhibiting a higher diversity in codon usage have generally a higher number of variable sites, a prerequisite to be potential phylogenetic markers. Thus the use of ENC and MILC scores, according to the new approach presented in this paper, is a way to study PCG sequences variability on a codon perspective. The best scores of both indices should allow to identify the more diverse PCGs in a approach complementary to the usual method based on evolutionary distances among orthologous sequences (e.g. [##REF##18065166##8##]). The assessment of genetic variability is an interesting point. Indeed some PCGs are standard marker for species recognition [##REF##16418261##23##] or have been extensively used as phylogenetic markers in Lepidoptera while others have received so far limited or no attention. Understanding the genetic diversity of each PCG is a prerequisite to determine its phylogenetic usefulness. The ENC and MILC values were calculated for all PCGs but <italic>atp8 </italic>that contains too a few codons to get reliable ENC/MILC estimations [##REF##16029499##22##]. Calculations were extended also to all 13 PCGs pooled as well as to the pooled PCGs belonging to α and β strands respectively. The scatter plot analysis is provided in Figure ##FIG##4##5##. As expected the greatest diversity in codon usage is found when all codons are considered. Good codon diversity is found also when all PCGs of α or β strands are considered. More interesting is the behavior of single genes. In this latter case sequences well established as phylogenetic markers (i.e. <italic>cox1</italic>, <italic>cob</italic>, <italic>nad5</italic>, and <italic>cox2</italic>) are intermixed with PCGs poorly or not considered by researchers (e.g. <italic>cox3</italic>, <italic>nad4</italic>, <italic>nad1</italic>, <italic>nad2</italic>). Our results suggest that the neglected PCGs should be considered as potential markers thus extending the number of mtDNA PCGs sampled for population as well as phylogenetic markers. Findings, based on codon diversity, must be integrated with direct comparisons of sequences [##REF##18065166##8##] that allow to better define the optimal task that each gene can perform i.e. to be used at low taxonomic level or at high taxonomic level.</p> <title>Transfer and ribosomal RNA genes</title> <p><italic>Ochrogaster </italic>genome has the characteristic 22 tRNAs set (Figure ##FIG##5##6##) present in most of animal mtDNAs [##REF##10101183##1##]. All tRNAs present the typical clover leaf secondary structure but <italic>trnS1 </italic>lacks the DHU stem. This feature is shared with the <italic>C. raphaelis </italic>mtDNA [##REF##16640732##6##] but is not a general feature of lepidopteran mtDNA as proved by <italic>A. honmai </italic>that has all tRNAs with a complete clover leaf structure [##REF##16488090##7##]. In general, the lack of DHU arm in <italic>trnS1 </italic>is a common condition in metazoan mtDNAs [##REF##11095730##24##].</p> <p>The <italic>trnA</italic>, <italic>trnD</italic>, <italic>trnG</italic>, <italic>trnK</italic>, <italic>trnL1</italic>, <italic>trnL2</italic>, <italic>trnQ</italic>, and <italic>trnS2 </italic>of <italic>Ochrogaster </italic>mtDNA show mismatches in their stems. Mismatches are located mostly in the acceptor and anticodon stems with a single exception represented by <italic>trnD </italic>that exhibits the mismatch on the TΨC stem. Mismatches on tRNA stems are known also for the <italic>trnA</italic>, <italic>trnL1</italic>, <italic>trnL2</italic>, and <italic>trnQ</italic>, of <italic>C. raphaelis </italic>[##REF##16640732##6##]. Mismatches observed in tRNAs are corrected through RNA-editing mechanisms that are well known for arthropod mtDNA [e.g. [##REF##11095730##24##]].</p> <p>Preliminary analysis performed on <italic>rrnL </italic>and <italic>rrnS </italic>of <italic>O. lunifer </italic>revealed that these genes are capable of folding into structures (data not shown) similar to those already produced for lepidopteran mitochondrial ribosomal subunits [##REF##18065166##8##,##UREF##3##25##,##REF##16483803##26##]. Further studies, that extend the taxon sampling, are currently in progress in our lab to better define <italic>rrnL </italic>and <italic>rrnS </italic>structures within the Thaumetopoeinae subfamily that includes also <italic>O. lunifer</italic>.</p> <title>Non coding regions</title> <p>The mtDNA genome of <italic>O. lunifer </italic>contains 7 intergenic spacers (s1–s7) spanning at least 15 bp (Figures ##FIG##0##1## and ##FIG##6##7##). The features of s1–s7 spacers are presented below with reference to the α strand for orientation and sequence motifs description.</p> <p>The s1 spacer, located between <italic>trnQ </italic>and <italic>nad2</italic>, appears to be the result of a duplicated segment (Figure ##FIG##6##7##). The s1 spacer is present in all 12 lepidopteran mtDNAs so far sequenced while it is absent in other insects [##REF##18065166##8##]. While the genomic location is constant the sequence divergence is high among species [##REF##18065166##8##]. 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.</p> <p>The s2 spacer, placed between <italic>trnC </italic>and <italic>trnY</italic>, derives from the triplication of a six nucleotides motif with minor changes (Figure ##FIG##6##7##). An 11 bp spacer between <italic>trnC </italic>and <italic>trnY </italic>is found also in the mtDNA of <italic>A. melete </italic>and shares the ACAATT motif with the s2 spacer of <italic>O. lunifer</italic>. Because no other known lepidopteran mtDNA exhibits such a spacer its presence in <italic>A. melete </italic>and <italic>O. lunifer </italic>has to be interpreted as the result of independent events.</p> <p>Spacer s3, located between <italic>nad3 </italic>and <italic>trnA</italic>, 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)<sub>10</sub>(TA)<sub>12</sub>. Spacers having the same genomic location, and containing TA microsatellites are found also in <italic>B. mori </italic>and <italic>B. mandarina </italic>mtDNA genomes.</p> <p>Spacer s4, inserted between <italic>trnE </italic>and <italic>trnF</italic>, contains a 5' microsatellite (TA)<sub>23</sub>, while the 3' half seems to be the triplication of a 10 nucleotides motif with some changes (Figure ##FIG##6##7##). A spacer characterized by a different motif (TATTA)<sub>31</sub>, but having the same genomic placement, is found in the <italic>A. honmai </italic>mtDNA genome.</p> <p>The spacer s5, located between <italic>trnS2 </italic>and <italic>nad1</italic>, contains the ATACTAA motif which is conserved across the Lepidoptera order [##REF##18065166##8##]. This motif is possibly fundamental to site recognition by the transcription termination peptide (mtTERM protein) [##REF##10076021##2##]. Spacer s5 is present in most insect mtDNAs even if the nucleotide sequence can be quite divergent [##REF##18065166##8##].</p> <p>The s6 spacer is located between <italic>trnS2 </italic>and -<italic>rrnL </italic>and exhibits a di-nucleotide microsatellite (TA)<sub>19 </sub>directly in contact with the 3' end of <italic>rrnL </italic>gene. To date spacer s6 is known only for the mtDNA of <italic>O. lunifer</italic>.</p> <p>The s7 spacer coincides with the A+T region. Several features common to the Lepidoptera A+T region [##REF##18065166##8##] are present in the s7 spacer. The O<sub>R</sub>β (origin of the β strand replication) is located 21 bp downstream from <italic>rrnS </italic>gene in <italic>B. mori </italic>[##REF##16118189##27##]. It contains the motif ATAGA followed by an 18 bp poly-T stretch. A very similar pattern occurs in <italic>O. lunifer </italic>where the ATAGA motif is located 17 bp downstream from <italic>rnnS </italic>gene and is followed by a 20 bp poly-T stretch (Figure ##FIG##6##7##). A microsatellite-like (AT)<sub>7</sub>(TA)<sub>3 </sub>element preceded by the ATTTA motif is present in the 3' third of <italic>O. lunifer </italic>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 [##REF##18065166##8##]. Finally a 10 bp poly-A is present immediately upstream <italic>trnM</italic>. This poly-T (in the β strand) element is still a common feature of the A+T region in Lepidoptera [##REF##18065166##8##,##REF##15045483##28##]. No large repeated segments were detected in the A+T region of <italic>O. lunifer</italic>. This arrangement is consistent with other lepidopteran A+T regions while markedly contrasts with patterns observed in other insect orders [##REF##18065166##8##,##UREF##4##29##].</p> <p>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 [##REF##18065166##8##]. Most parts of spacers of <italic>O. lunifer </italic>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. [##REF##15114420##18##,##REF##16854241##30##]].</p>
<title>Conclusion</title> <p>The mitochondrial genome of <italic>O. lunifer </italic>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 <italic>trnM </italic>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 <italic>O. lunifer </italic>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 <italic>atp6</italic>, <italic>cox1</italic>, <italic>cox2, cox3</italic>, <italic>cob</italic>, <italic>nad1</italic>, <italic>nad2, nad4</italic>, and <italic>nad5 </italic>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 <italic>O. lunifer </italic>genome lead us to suggest an important role of DNA slippage as possible mechanism to produce spacers during replication.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>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 <italic>Ochrogaster lunifer</italic>.</p> <title>Results</title> <p>The mitochondrial genome of <italic>O. lunifer </italic>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 <italic>trnM</italic>. The 77.84% A+T content of its α strand is the lowest among known lepidopteran genomes. The mitochondrial genome of <italic>O. lunifer </italic>exhibits one of the most marked C-skew among available insect Pterygota genomes. The protein-coding genes have typical mitochondrial start codons except for <italic>cox1 </italic>that present an unusual CGA. The <italic>O. lunifer </italic>genome exhibits the less biased synonymous codon usage among lepidopterans. Comparative genomics analysis study identified <italic>atp6</italic>, <italic>cox1</italic>, <italic>cox2 </italic>as <italic>cox3</italic>, <italic>cob</italic>, <italic>nad1</italic>, <italic>nad2, nad4</italic>, and <italic>nad5 </italic>as potential markers for population genetics/phylogenetics studies. A peculiar feature of <italic>O. lunifer </italic>mitochondrial genome it that the intergenic spacers are mostly made by repetitive sequences.</p> <title>Conclusion</title> <p>The mitochondrial genome of <italic>O. lunifer </italic>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 <italic>O. lunifer </italic>genome supports the role of DNA slippage as possible mechanism to produce spacers during replication.</p>
<title>Abbreviations</title> <p>mtDNA: mitochondrial DNA; <italic>atp6 </italic>and <italic>atp8</italic>: ATP synthase subunits 6 and 8; <italic>cob</italic>: apocytochrome b; <italic>cox1</italic>-<italic>3</italic>: cytochrome c oxidase subunits 1–3; <italic>nad1</italic>-<italic>6 </italic>and <italic>nad4L</italic>: NADH dehydrogenase subunits 1–6 and 4L; <italic>rrnS </italic>and <italic>rrnL</italic>: small and large subunit ribosomal RNA (rRNA) genes; <italic>trnX</italic>: transfer RNA (tRNA) genes, where X is the one-letter abbreviation of the corresponding amino acid; s1–s7: mitochondrial genomic spacers; A+T region: the putative control region; PCG: protein coding gene; RSCU: Relative Synonymous Codon Usage; ENC, MILC: Measure Independent of Length and Composition; aa: amino acids; nt: nucleotides; bp: base pairs.</p> <title>Authors' contributions</title> <p>PS and MS carried out the molecular experiments. AB and EN designed and coordinated all experiments. EN performed the genomic analyses. All authors contributed to the manuscript and then read and approved the final version.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We express our sincere thanks to Myron P. Zalucki (School of Integrative Biology, University of Queensland, Brisbane, Australia) who kindly provided the specimen of <italic>Ochrogaster lunifer </italic>used in present study. We thank Filippo Calore (Albignasego, Padova, Italy) who painted the icon of <italic>O. lunifer </italic>included in Figure ##FIG##0##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.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Map of the mitochondrial genome of <italic>O. lunifer</italic></bold>. Genes coded in the α strand (clockwise orientation) are blue or cyan colored. Genes coded in the β strand (anti-clockwise orientation) are red or orange colored. Alternation of colors was applied for clarity. Start, first position along α strand; end, last position along α strand; size, size of the sequence; inc, intergenic nucleotides; fcd, first codon; scd, stop codon. Incomplete stop codons are presented with parentheses. Negative inc values refer to overlapping nucleotides for genes located in the same or different strands. Gene names are the standard abbreviations used in this paper; tRNA genes are indicated by the single letter IUPAC-IUB abbreviation for their corresponding amino acid in the draw. s1–s7, intergenic spacers.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>A-skew vs. A+T% and G-skew vs. G+C% in the Pteryogota mtDNAs</bold>. Values were calculated on α strands for full lengh mtDNA genomes. The X axis provides the skews values, while the Y axis provides the A+T/G+C values. Named of species are colored according to their taxonomic placement at Order level (see Table 1).</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Codon distribution in lepidopteran mtDNAs</bold>. Numbers to the left refer to the total number of codon. CDspT, codons per thousands codons. Codon Families are provided on the x axis.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Relative Synonymous Codon Usage (RSCU) in lepidopteran mtDNAs</bold>. Codon Families are provided on the x axis. Red-colored codon, codon not present in the genome. Codon Families are provided on the x axis.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Scatter plot graphic of MILC vs. ENC calculated for PCGs of lepidopteran mtDNAs</bold>. Dots correspond to average values calculated for different genes. PCGs on α strand are blue-colored, PCGs on β strand are red colored. All pooled PCGs are presented as a green dot plot. Genes nomenclature as in main text.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>Secondary structures of transfer tRNAs in <italic>O. lunifer </italic>mtDNA.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Genomic spacers in the mtDNA of O. lunifer</bold>. The sequences of spacers are those present in the α strand.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>List of taxa analyzed in present paper</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Order</bold></td><td align="left"><bold>Family</bold></td><td align="left"><bold>Species</bold></td><td align="left"><bold>Acc. number</bold></td><td align="left"><bold>Reference</bold></td></tr></thead><tbody><tr><td align="left"><bold>Orthoptera</bold></td><td align="left">Acrididae</td><td align="left"><italic>Locusta migratoria</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001712">NC_001712</ext-link></td><td align="left">[##REF##8587138##42##]</td></tr><tr><td/><td align="left">Acrididae</td><td align="left"><italic>Oxya chinensis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_010219">NC_010219</ext-link></td><td align="left">Hang and Zhang, unpublished</td></tr><tr><td/><td align="left">Tettigonidae</td><td align="left"><italic>Anabrus simplex</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009967">NC_009967</ext-link></td><td align="left">[##REF##17316330##43##]</td></tr><tr><td/><td align="left">Tettigonidae</td><td align="left"><italic>Ruspolia dubia</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009876">NC_009876</ext-link></td><td align="left">[##REF##17893726##44##]</td></tr><tr><td/><td align="left">Gryllotalpidae</td><td align="left"><italic>Gryllotalpa orientalis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006678">NC_006678</ext-link></td><td align="left">[##REF##15950403##45##]</td></tr><tr><td align="left"><bold>Isoptera</bold></td><td align="left">Rhinotermitidae</td><td align="left"><italic>Reticulitermes flavipes</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009498">NC_009498</ext-link></td><td align="left">[##REF##17546084##13##]</td></tr><tr><td/><td align="left">Rhinotermitidae</td><td align="left"><italic>Reticulitermes hageni</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009501">NC_009501</ext-link></td><td align="left">[##REF##17546084##13##]</td></tr><tr><td/><td align="left">Rhinotermitidae</td><td align="left"><italic>Reticulitermes virginicus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009500">NC_009500</ext-link></td><td align="left">[##REF##17546084##13##]</td></tr><tr><td/><td align="left">Rhinotermitidae</td><td align="left"><italic>Reticulitermes santonensis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009499">NC_009499</ext-link></td><td align="left">[##REF##17546084##13##]</td></tr><tr><td align="left"><bold>Manthophasmatodea</bold></td><td align="left">Mantophasmatidae</td><td align="left"><italic>Sclerophasma paresisense</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007701">NC_007701</ext-link></td><td align="left">[##REF##16321547##46##]</td></tr><tr><td align="left"><bold>Mantodea</bold></td><td align="left">Mantidae</td><td align="left"><italic>Tamolanica tamolana</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007702">NC_007702</ext-link></td><td align="left">[##REF##16321547##46##]</td></tr><tr><td align="left"><bold>Blattaria</bold></td><td align="left">Blattidae</td><td align="left"><italic>Periplaneta fuliginosa</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006076">NC_006076</ext-link></td><td align="left">[##REF##15271216##47##]</td></tr><tr><td align="left"><bold>Plecoptera</bold></td><td align="left">Pteronarcyidae</td><td align="left"><italic>Pteronarcys princeps</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006133">NC_006133</ext-link></td><td align="left">[##REF##16936790##48##]</td></tr><tr><td align="left"><bold>Phthiraptera</bold></td><td align="left">Boopidae</td><td align="left"><italic>Heterodoxus macropus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002651">NC_002651</ext-link></td><td align="left">[##REF##11319269##16##]</td></tr><tr><td/><td align="left">Philopteridae</td><td align="left"><italic>Bothriometopus macrocnemis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009983">NC_009983</ext-link></td><td align="left">[##REF##17925995##49##]</td></tr><tr><td/><td align="left">Philopteridae</td><td align="left"><italic>Campanulotes bidentatus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007884">NC_007884</ext-link></td><td align="left">[##REF##16469069##15##]</td></tr><tr><td align="left"><bold>Hemiptera</bold></td><td align="left">Aleyrodidae</td><td align="left"><italic>Aleurochiton aceris</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006160">NC_006160</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aleyrodidae</td><td align="left"><italic>Aleurodicus dugesii</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005939">NC_005939</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aleyrodidae</td><td align="left"><italic>Bemisia tabaci</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006279">NC_006279</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aleyrodidae</td><td align="left"><italic>Neomaskellia andropogonis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006159">NC_006159</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aleyrodidae</td><td align="left"><italic>Tetraleurodes acaciae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006292">NC_006292</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aleyrodidae</td><td align="left"><italic>Trialeurodes vaporariorum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006280">NC_006280</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aphididae</td><td align="left"><italic>Schizaphis graminum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006158">NC_006158</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Cicadellidae</td><td align="left"><italic>Homalodisca coagulata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006899">NC_006899</ext-link></td><td align="left">Baumann and Baumann, unpublished</td></tr><tr><td/><td align="left">Psyllidae</td><td align="left"><italic>Pachypsylla venusta</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006157">NC_006157</ext-link></td><td align="left">[##REF##15291971##14##]</td></tr><tr><td/><td align="left">Aphrophoridae</td><td align="left"><italic>Philaenus spumarius</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005944">NC_005944</ext-link></td><td align="left">[##REF##15729396##50##]</td></tr><tr><td/><td align="left">Reduviidae</td><td align="left"><italic>Triatoma dimidiata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002609">NC_002609</ext-link></td><td align="left">[##REF##11437912##51##]</td></tr><tr><td align="left"><bold>Psocoptera</bold></td><td align="left">Lepidopsocidae</td><td align="left">Lepidopsocid sp. RS-2001</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004816">NC_004816</ext-link></td><td align="left">[##REF##12832626##52##]</td></tr><tr><td align="left"><bold>Thysanoptera</bold></td><td align="left">Thripidae</td><td align="left"><italic>Thrips imaginis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004371">NC_004371</ext-link></td><td align="left">[##REF##12644556##53##]</td></tr><tr><td align="left"><bold>Coleoptera</bold></td><td align="left">Cerambycidae</td><td align="left"><italic>Anoplophora glabripennis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008221">NC_008221</ext-link></td><td align="left">An et al., unpublished</td></tr><tr><td/><td align="left">Chrysomelidae</td><td align="left"><italic>Crioceris duodecimpunctata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003372">NC_003372</ext-link></td><td align="left">[##REF##12644408##54##]</td></tr><tr><td/><td align="left">Elateridae</td><td align="left"><italic>Pyrophorus divergens</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009964">NC_009964</ext-link></td><td align="left">[##REF##17942246##55##]</td></tr><tr><td/><td align="left">Lampyridae</td><td align="left"><italic>Pyrocoelia rufa</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003970">NC_003970</ext-link></td><td align="left">[##REF##15288070##56##]</td></tr><tr><td/><td align="left">Tenebrionidae</td><td align="left"><italic>Tribolium castaneum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003081">NC_003081</ext-link></td><td align="left">[##REF##12644407##57##]</td></tr><tr><td align="left"><bold>Lepidoptera</bold></td><td/><td/><td/><td/></tr><tr><td align="right">Tortricoidea</td><td align="left">Tortricidae</td><td align="left"><italic>Adoxophyes honmai</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008141">NC_008141</ext-link></td><td align="left">[##REF##16488090##7##]</td></tr><tr><td align="right">Bombycoidea</td><td align="left">Saturniidae</td><td align="left"><italic>Antheraea pernyi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004622">NC_004622</ext-link></td><td align="left">Liu et al., unpublished</td></tr><tr><td align="right">Bombycoidea</td><td align="left">Bombycidae</td><td align="left"><italic>Bombyx mandarina</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003395">NC_003395</ext-link></td><td align="left">[##REF##12140251##5##]</td></tr><tr><td align="right">Bombycoidea</td><td align="left">Bombycidae</td><td align="left"><italic>Bombyx mori</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002355">NC_002355</ext-link></td><td align="left">Lee et al., unpublished</td></tr><tr><td align="right">Bombycoidea</td><td align="left">Saturniidae</td><td align="left"><italic>Saturnia boisduvalii</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_010613">NC_010613</ext-link></td><td align="left">[##REF##18337026##9##]</td></tr><tr><td align="right">Geometroidea</td><td align="left">Geometridae</td><td align="left"><italic>Phthonandria atrilineata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_010522">NC_010522</ext-link></td><td align="left">Yang et al., unpublished</td></tr><tr><td align="right">Papilionoidea</td><td align="left">Pieridae</td><td align="left"><italic>Artogeia melete</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_010568">NC_010568</ext-link></td><td align="left">Hong et al., unpublished</td></tr><tr><td align="right">Papilionoidea</td><td align="left">Lycaenidae</td><td align="left"><italic>Coreana raphaelis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007976">NC_007976</ext-link></td><td align="left">[##REF##16640732##6##]</td></tr><tr><td align="right">Sphingoidea</td><td align="left">Sphingidae</td><td align="left"><italic>Manduca sexta</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="EU286785">EU286785</ext-link></td><td align="left">[##REF##18065166##8##]</td></tr><tr><td align="right"><bold>Noctuoidea</bold></td><td align="left"><bold>Notodontidae</bold></td><td align="left"><bold><italic>Ochrogaster lunifer</italic></bold></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AM946601">AM946601</ext-link></td><td align="left"><bold>This paper</bold></td></tr><tr><td align="right">Pyraloidea</td><td align="left">Crambidae</td><td align="left"><italic>Ostrinia furnacalis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003368">NC_003368</ext-link></td><td align="left">[##REF##15951845##10##]</td></tr><tr><td align="right">Pyraloidea</td><td align="left">Crambidae</td><td align="left"><italic>Ostrinia nubilalis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003367">NC_003367</ext-link></td><td align="left">[##REF##15951845##10##]</td></tr><tr><td align="left"><bold>Diptera</bold></td><td align="left">Ceratopogonidae</td><td align="left"><italic>Culicoides arakawae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009809">NC_009809</ext-link></td><td align="left">Matsumoto, unpublished</td></tr><tr><td/><td align="left">Culicidae</td><td align="left"><italic>Aedes albopictus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006817">NC_006817</ext-link></td><td align="left">Ho et al., unpublished</td></tr><tr><td/><td align="left">Culicidae</td><td align="left"><italic>Aedes aegypti</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_010241">NC_010241</ext-link></td><td align="left">Lobo et al., unpublished</td></tr><tr><td/><td align="left">Culicidae</td><td align="left"><italic>Anopheles gambiae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002084">NC_002084</ext-link></td><td align="left">[##REF##9087549##58##]</td></tr><tr><td/><td align="left">Culicidae</td><td align="left"><italic>Anopheles quadrimaculatus </italic>A</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_000875">NC_000875</ext-link></td><td align="left">[##REF##8112570##59##]</td></tr><tr><td/><td align="left">Calliphoridae</td><td align="left"><italic>Cochliomyia hominivorax</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002660">NC_002660</ext-link></td><td align="left">[##REF##11029671##60##]</td></tr><tr><td/><td align="left">Calliphoridae</td><td align="left"><italic>Lucilia sericata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009733">NC_009733</ext-link></td><td align="left">Cibrario et al., unpublished</td></tr><tr><td/><td align="left">Calliphoridae</td><td align="left"><italic>Chrysomya putoria</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002697">NC_002697</ext-link></td><td align="left">[##REF##15363841##61##]</td></tr><tr><td/><td align="left">Drosophilidae</td><td align="left"><italic>Drosophila melanogaster</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001709">NC_001709</ext-link></td><td align="left">[##REF##8825764##62##]</td></tr><tr><td/><td align="left">Drosophilidae</td><td align="left"><italic>Drosophila mauritiana</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005779">NC_005779</ext-link></td><td align="left">[##REF##10903372##63##]</td></tr><tr><td/><td align="left">Drosophilidae</td><td align="left"><italic>Drosophila sechellia</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005780">NC_005780</ext-link></td><td align="left">[##REF##10903372##63##]</td></tr><tr><td/><td align="left">Drosophilidae</td><td align="left"><italic>Drosophila simulans</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005781">NC_005781</ext-link></td><td align="left">[##REF##10903372##63##]</td></tr><tr><td/><td align="left">Drosophilidae</td><td align="left"><italic>Drosophila yakuba</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001322">NC_001322</ext-link></td><td align="left">[##REF##3001325##64##]</td></tr><tr><td/><td align="left">Oestridae</td><td align="left"><italic>Dermatobia hominis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006378">NC_006378</ext-link></td><td align="left">Azeredo-Espin et al., unpublished</td></tr><tr><td/><td align="left">Muscidae</td><td align="left"><italic>Haematobia irritans</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007102">NC_007102</ext-link></td><td align="left">Lessinger et al., unpublished</td></tr><tr><td/><td align="left">Nemestrinidae</td><td align="left"><italic>Trichophthalma punctata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008755">NC_008755</ext-link></td><td align="left">[##UREF##8##65##]</td></tr><tr><td/><td align="left">Syrphidae</td><td align="left"><italic>Simosyrphus grandicornis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008754">NC_008754</ext-link></td><td align="left">[##UREF##8##65##]</td></tr><tr><td/><td align="left">Tabanidae</td><td align="left"><italic>Cydistomyia duplonotata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008756">NC_008756</ext-link></td><td align="left">[##UREF##8##65##]</td></tr><tr><td/><td align="left">Tephritidae</td><td align="left"><italic>Ceratitis capitata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_000857">NC_000857</ext-link></td><td align="left">[##REF##10762421##66##]</td></tr><tr><td/><td align="left">Tephritidae</td><td align="left"><italic>Bactrocera carambolae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009772">NC_009772</ext-link></td><td align="left">Ye et al., unpublished</td></tr><tr><td/><td align="left">Tephritidae</td><td align="left"><italic>Bactrocera dorsalis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008748">NC_008748</ext-link></td><td align="left">Yu et al., unpublished</td></tr><tr><td/><td align="left">Tephritidae</td><td align="left"><italic>Bactrocera oleae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005333">NC_005333</ext-link></td><td align="left">[##REF##14986921##67##]</td></tr><tr><td/><td align="left">Tephritidae</td><td align="left"><italic>Bactrocera papayae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009770">NC_009770</ext-link></td><td align="left">Ye et al., unpublished</td></tr><tr><td/><td align="left">Tephritidae</td><td align="left"><italic>Bactrocera philippinensis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_009771">NC_009771</ext-link></td><td align="left">Ye et al., unpublished</td></tr><tr><td align="left"><bold>Hymenoptera</bold></td><td align="left">Apidae</td><td align="left"><italic>Apis mellifera ligustica</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001566">NC_001566</ext-link></td><td align="left">[##REF##8417993##68##]</td></tr><tr><td/><td align="left">Apidae</td><td align="left"><italic>Bombus ignitus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ870926">DQ870926</ext-link></td><td align="left">[##REF##17321076##69##]</td></tr><tr><td/><td align="left">Apidae</td><td align="left"><italic>Melipona bicolor</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004529">NC_004529</ext-link></td><td align="left">Silvestre and Arias, unpublished</td></tr><tr><td/><td align="left">Vanhornidae</td><td align="left"><italic>Vanhornia eucnemidarum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_008323">NC_008323</ext-link></td><td align="left">[##REF##16936784##70##]</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Additional file ##SUPPL##0##1##. List of primers and PCR conditions used in the sequencing of <italic>Ochogaster lunifer </italic>mtDNA.</p></caption></supplementary-material>
<graphic xlink:href="1471-2164-9-331-1"/> <graphic xlink:href="1471-2164-9-331-2"/> <graphic xlink:href="1471-2164-9-331-3"/> <graphic xlink:href="1471-2164-9-331-4"/> <graphic xlink:href="1471-2164-9-331-5"/> <graphic xlink:href="1471-2164-9-331-6"/> <graphic xlink:href="1471-2164-9-331-7"/>
<media xlink:href="1471-2164-9-331-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Resh", "Card\u00e9", "Resh VH, CArd\u00e9 RG"], "given-names": ["VH", "RG"], "article-title": ["Insecta, Overview"], "source": ["Encyclopedia of Insects"], "year": ["2003"], "publisher-name": ["Academic Press, Burlington MA, USA"], "fpage": ["564"], "lpage": ["566"], "comment": ["1266pp."]}, {"surname": ["Powell", "Resh VH, CArd\u00e9 RG"], "given-names": ["JA"], "article-title": ["Lepidoptera (Moths, Butterflies)"], "source": ["Encyclopedia of Insects"], "year": ["2003"], "publisher-name": ["Academic Press, Burlington MA, USA"], "fpage": ["631"], "lpage": ["663"], "comment": ["1266pp."]}, {"surname": ["Patthy"], "given-names": ["L"], "source": ["Protein evolution"], "year": ["2008"], "edition": ["2"], "publisher-name": ["Blackwell, London"], "fpage": ["374"]}, {"surname": ["Niehuis", "Naumann", "Misof"], "given-names": ["O", "CM", "B"], "article-title": ["Identification of evolutionary conserved structural elements in the mt SSU Rrna of Zygaenoidea (Lepidoptera): a comparative sequence analysis"], "source": ["Org Divers Evol"], "year": ["2006"], "volume": ["6"], "fpage": ["17"], "lpage": ["32"]}, {"surname": ["Zhang", "Hweitt"], "given-names": ["DX", "GM"], "article-title": ["Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies"], "source": ["Biochem Syst Ecol"], "year": ["1997"], "volume": ["25"], "fpage": ["99"], "lpage": ["120"]}, {"surname": ["Patwary", "Kenchington", "Bird", "Zouros"], "given-names": ["MU", "EL", "CJ", "E"], "article-title": ["The use of random amplified polymorphic DNA markers in genetic studies of the sea scallop "], "italic": ["Placopecten magellanicus "], "source": ["J Shellfish Res"], "year": ["1994"], "volume": ["13"], "fpage": ["547"], "lpage": ["553"]}, {"surname": ["Simon", "Frati", "Beckenbach", "Crespi", "Liu", "Flook"], "given-names": ["C", "F", "A", "B", "H", "P"], "article-title": ["Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers"], "source": ["Ann Entomol Soc Am"], "year": ["1994"], "volume": ["87"], "fpage": ["651"], "lpage": ["704"]}, {"surname": ["Simon", "Buckley", "Frati", "Stewart", "Beckenbach"], "given-names": ["C", "TR", "F", "JB", "AT"], "article-title": ["Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA"], "source": ["Annual Review of Ecology, Evolution, and Systematics"], "year": ["2006"], "volume": ["37"], "fpage": ["545"], "lpage": ["579"]}, {"surname": ["Cameron", "Lambkin", "Barker", "Whiting"], "given-names": ["SL", "CL", "SC", "MF"], "article-title": ["A mitochondrial genome phylogeny of Diptera: whole genome sequence data accurately resolve relationships over broad timescales with high precision"], "source": ["Syst Entomol"], "year": ["2007"], "volume": ["32"], "fpage": ["40"], "lpage": ["59"]}]
BMC Genomics. 2008 Jul 15; 9:331