Skip to main content
Download PDF

The mitochondrial genome of the 'twisted-wing parasite' Mengenilla australiensis (Insecta, Strepsiptera): a comparative study

BMC Genomics volume 10, Article number: 603 (2009) Cite this article

Abstract

Background

Strepsiptera are an unusual group of sexually dimorphic, entomophagous parasitoids whose evolutionary origins remain elusive. The lineage leading to Mengenilla australiensis (Family Mengenillidae) is the sister group to all remaining extant strepsipterans. It is unique in that members of this family have retained a less derived condition, where females are free-living from pupation onwards, and are structurally much less simplified. We sequenced almost the entire mitochondrial genome of M. australiensis as an important comparative data point to the already available genome of its distant relative Xenos vesparum (Family Xenidae). This study represents the first in-depth comparative mitochondrial genomic analysis of Strepsiptera.

Results

The partial genome of M. australiensis is presented as a 13421 bp fragment, across which all 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes and 18 transfer RNA (tRNA) sequences are identified. Two tRNA translocations disrupt an otherwise ancestral insect mitochondrial genome order. A+T content is measured at 84.3%, C-content is also very skewed. Compared with M. australiensis, codon bias in X. vesparum is more balanced. Interestingly, the size of the protein coding genome is truncated in both strepsipterans, especially in X. vesparum which, uniquely, has 4.3% fewer amino acids than the average holometabolan complement. A revised assessment of mitochondrial rRNA secondary structure based on comparative structural considerations is presented for M. australiensis and X. vesparum.

Conclusions

The mitochondrial genome of X. vesparum has undergone a series of alterations which are probably related to an extremely derived lifestyle. Although M. australiensis shares some of these attributes; it has retained greater signal from the hypothetical most recent common ancestor (MRCA) of Strepsiptera, inviting the possibility that a shift in the mitochondrial selective environment might be related to the specialization accompanying the evolution of a small, morphologically simplified completely host-dependent lifestyle. These results provide useful insights into the nature of the evolutionary transitions that accompanied the emergence of Strepsiptera, but we emphasize the need for adequate sampling across the order in future investigations concerning the extraordinary developmental and evolutionary origins of this group.

Background

Strepsiptera are an unusual group of obligate endoparasitoid insects [1]. They occur as a small (approx. 600 spp.) monophyletic insect order with uncertain evolutionary origins, inclusive of any clear understanding over the group's nearest extant relative. Strepsiptera parasitize 7 orders of insects, including silverfish (Thysanura); cockroaches (Blattaria); mantids (Mantodea); crickets and grasshoppers (Orthoptera); bugs (Hemiptera); wasps, ants and bees (Hymenoptera) and flies (Diptera). Understandably, the very uniqueness of Strepsiptera makes their placement within insects using morphological taxonomic methods a contentious task. Strepsiptera have been placed as the sister-group to myriad different insect groups, from beetles to true flies, even being placed outside of Holometabola: each hypothesis being founded on one or two 'key' characteristics that at one time or another have come under question [2–15]. Confident assertions of classification have especially been restricted because intermediate forms have largely gone extinct and are unrecorded in the fossil record (see [16] for a notable exception). Simplification of gross morphology during strepsipteran evolutionary specialization can also be viewed as a significant component of this problem: observable morphological variation is low and unevenly distributed between extremely dimorphic sexes.

Females are especially strongly simplified (lacking in most typical adult characters such as wings, legs or mouthparts), and in most strepsipteran species, they remain in the living host until the end of the reproductive cycle. Conversely, males metamorphose in a typical holometabolan fashion; developing wings and a usual suite of adult insect characters (such as antennae, mouthparts and compound eyes) only to leave the host immediately in search of a mate; usually in the form of a female-containing host [17–19]. Males deliver sperm through a brood canal opening in the cephalothorax (the modified "head") of the female, who as a reproductive adult is only partially exposed to the environment, as an extrusion between the tergites or sternites of living hosts. After fertilization, females are capable of producing many hundreds of thousands of active 1st instar larvae who emerge from the cephalothorax.

Members of the strepsipteran lineage, Mengenillidae, are the sister-group to Stylopidia, a clade that includes all other extant Strepsiptera [[20], unpublished data]. This family represents a transitional phase in the evolutionary specialization of Strepsiptera, whereby both sexes leave the host before pupation, and females do not reproduce or release progeny in the unusual fashion outlined above. The level of simplification in free living mengenillid females is much less extreme than in Stylopidia (e.g. X. vesparum); legs, mouthparts and compound eyes are all still present, although strongly reduced. The invention of a completely endoparasitic female was probably one of the most important novelties leading to the radiation of this unique group [21]. Evolutionary studies attempting to unravel the developmental and ecological phenomena that make Strepsiptera so biologically interesting therefore require that species from before this important transition occupy a central role in research.

Insect mitochondrial genomes provide a useful medium to deepen and connect comprehension of microevolutionary forces of populations, like neutral drift and selective sweeps, to macroevolutionary events affecting species and/or deeper levels of divergence. The abundance of mitochondria in most metazoan cell and tissue types makes mtDNA an easily obtainable, universally plentiful marker, where a lack of introns or duplicate genes and non-coding variable spacer regions make amplification of mtDNA relatively uncomplicated [22] (although pseudogenes co-opted by the nucleus from the mitochondrial genome can create separate difficulties [23–26]). The near-complete sequence of the mt genome of Xenos vesparum is available, and as a member of the more specialized clade Stylopidia, it is important to place it into context. The mitochondrial genome of Mengenilla australiensis (Family Mengenillidae) was therefore chosen for sequencing. Genome arrangement, nucleotide content, codon usage and the secondary structure of ribosomal RNA genes are each comparatively assessed between two distant Strepsiptera, and more widely across Holometabola (a major division of insects defined by the process of complete metamorphosis, including the 4 very species-rich orders; Lepidoptera, Diptera, Hymenoptera, Coleoptera and 6 smaller orders).

Results and discussion

Genome composition

The near-entire mt genome fragment of M. australiensis is presented as a 13421 bp sequence, across which 13 protein-coding genes (PCGs), a large (16S) and a small (12S) subunit rRNA gene (rrnL and rrnS respectively) and 18 tRNA sequences can be identified (Figure 1). PCGs are in expected positions, but two tRNA translocations disrupt an otherwise ancestral genome arrangement. Serine1 (S1: AGN) moves to a position between Alanine (A) and Arginine (R), whereas Valine (V) is either lost, or transferred to a position in a region that spans the flanking region, potentially anywhere between tRNA-Ile and nad2. All start codons across the PCGs (nad2 is only partial) begin with the typical M- or I- residue [27], and end in TAA, TAG, TA, or T. The control region itself, including up to 4 flanking tRNA genes could not be amplified. A variety of approaches were explored, none of which were successful in spanning the presumed A+T region. Similar difficulties were encountered during the amplification of the mt genome of X. vesparum [28]. In our study, genome-specific primers designed to span the A+T region also performed well when reverse compliment versions were implemented. It is hypothesized that the region is unusually long and/or too problematic for even advanced Taq polymerases to amplify, as a result of extreme repetitiveness or secondary structural folding issues (or both).

Figure 1
figure 1

Summary schematic of the mt genome of M. australiensis. Genes transcribed on the leading (α) strand are given in blue (PCGs) or light blue (tRNAs) and run clockwise. Genes of the opposite polarity are given in red (PCGs) or orange (tRNAs). Abbreviations are described in the text.

Full size image

Figure 2 compares the genome organization of M. australiensis and X. vesparum against the inferred ancestral arrangement [29, 30]; the former having just two gene order rearrangements, compared to 4 in the latter. In the context of variation across a dataset containing 68 mitochondrial genomes taken from 6 major holometabolan groups compiled from Genbank (see methods), this would be consistent with a mostly ancestral arrangement being retained between the origins of the major holometabolan radiations, since every order contains representative ancestral (or near-ancestral as in M. australiensis) genome arrangements. Within certain orders, especially Hymenoptera (and certain hemimetabolous groups not discussed here [31]), lineages have undergone significant independent modifications (see [32] for a full review of hymenopteran mt genome evolution). Some synapomorphic alterations to this ancestral arrangement also exist in certain orders. These include a tRNA arrangement shift from IQM to MIQ across Lepidoptera [33] (outside of the fragment shown in Figure 2), and a likely shift from WCY to CWY in Neuroptera - a modification that is not shared by the other neuropterid orders [34]. But largely, the 68 genomes corroborate a hypothesis positing that the holometabolan orders emerged during a (relatively) brief period of ancient rapid radiation. Or in the context of phylogenetic tree shape, as an initial phase of short internal nodes, during which there would be insufficient time for changes in mitochondrial genome organization to become fixed between orders, followed by a longer period of external node expansion [35] in which enough time would pass for modifications to emerge independently along multiple branches within and between orders. With the inclusion of M. australiensis, whose genome does not share any of the alterations of X. vesparum, and is more representative of an ancestral insect arrangement, this also appears consistent for Strepsiptera: although greater sampling is required for a general picture of mt genome structural evolution to emerge.

Figure 2
figure 2

Linear comparison of partial mt genome structure in M. australiensis and X. vesparum. Elements in red denote deviations from the ancestral holometabolan arrangement [29, 30]. Readers are referred to the primary research for a complete catalogue of genome rearrangements in Holometabola: Hymenoptera [32]; Coleoptera [88–93]; Neuropterida [34, 94, 95], a shift from WCY to CWY is probably synapomorphic for Neuroptera but not for other neuropterid orders; Lepidoptera [33, 75–81] are characterized by a synapomorphic tRNA rearrangement MIQ, from IQM, 5' of nad2 (not shown); Diptera [60–74].

Full size image

The lengths of PCGs are truncated in both strepsipteran species, and especially so in X. vesparum. Across the 68 holometabolan dataset, M. australiensis has a reduction in mean content of 6.2 amino acids per gene. In X. vesparum, this reaches 11.8 amino acid deletions per gene. X. vesparum has the shortest nad2, cox1, cox2, atp8, atp6, cox3, nad3 and nad4 gene (Additional file 1), representing a 4.3% loss in amino acids from the average holometabolan complement. Its total coding genome is shorter by 154 amino acids, the next shortest is M. australiensis (81) followed by Bombus ignitus (67) and M. bicolor (45) (Figure 3). Given that substantial loss of gene content might be expected to severely compromise gene functionality and efficiency, it would be of considerable interest to investigate whether this kind of genomic streamlining is related in any way to the peculiar lifestyle of small endoparasitoid insects. In particular, does bottlenecking in strepsipteran populations lead to slightly deleterious mutations, like codon deletions, being fixed through random drift? Strepsipteran metapopulations could be experiencing the extreme extinction-recolonization population dynamics required for low effective population sizes to become influential, but there are no a priori reasons to suspect these should be so different from other hymenopteran host-parasitoid systems that depend on similar insect host groups (e.g. Evania appendigaster, Venturia canescens; both sampled in this study). Alternatively, bottlenecking could result from the population dynamics of mitochondria themselves, through very low numbers of these organelles being passed via the germ line of strepsipteran eggs, which are known to be extremely small [36].

Figure 3
figure 3

Comparison of the total mt coding genome size (in amino acids) across Holometabola. Black: Thysanura; Green: Diptera; Orange: Lepidoptera; Red: Hymenoptera; Blue: Coleoptera; Turquoise: Neuropterida; Burgundy: Strepsiptera. Each bar represents one species' genome. The superimposed line indicates the mean coding genome size across all 68 taxa.

Full size image

Analyses to date do not reveal general patterns of mitochondrial genome evolution across holometabolan parasitic lineages [37, 38], but these have concentrated largely on the comparative analysis of genome organization (gene order/orientation) itself. Within orders like Hymenoptera, although a correlation between the extent of mt genome modification and parasitism appears to be lacking [39, 40] (with an additional implication that the position of mt genes might largely be neutral [32]), the truncation of PCGs as presented in this report, probably represents a separate issue. Determining the extent to which PCG truncation occurs within Strepsiptera must be a target for future investigations. Unfortunately, direct measures of insect mitochondrial effective population size have not yet been investigated.

Overall nucleotide composition is typically A+T rich in M. australiensis (84.3%; 6th highest amongst holometabolan insects), and like in X. vesparum, C-skew is high (+0.27). A scatter of the variation in skew and A+T% across the 68 holometabolan genomes (plus 3 thysanuran genomes) is shown in Figures 4A and 4B. Notably, in Figure 4B, of the most C-skewed and G+C% poor genomes (occupying the lower right quarter of this graph), 9 (of 12) data points are hymenopteran, 2 (of 2) are strepsipteran and 1 (of 12) is lepidopteran. A similar graph that including hemimetabolous insects [41] shows that only 1 other genome (Schizaphis graminums; Aphidoidea, 1 of 11 hemipteran genomes) is so C-skewed and A+T% rich. Conversely, the thysanuran genomes and a subset of Coleoptera are very A-skewed and distinctly more A+T% balanced (Figure 4A). This pattern is discussed in more detail in the following section.

Figure 4
figure 4

Graphical summary of nucleotide content across 68 holometabolan mt genomes. Measured in bp percentage (Y-axis) and level of nucleotide skew (X-axis) (described in methods, following [41]). Green: Diptera; Orange: Lepidoptera; Red: Hymenoptera; Blue: Coleoptera; Turquoise: Neuropterida; Burgundy: Strepsiptera. 1 and 2 refer to M. australiensis and X. vesparum respectively. A) A+T% vs A-skew B) G+C% vs C-skew. 3 genomes from Thysanura (in black), the host group for M. australiensis are also shown (see methods).

Full size image

Codon Usage

PCGs from Strepsiptera and the holometabolan dataset were imported into INCA v2.1 [42] to analyze codon usage. M. australiensis and X. vesparum are compared in Figure 5: codon bias is significantly relaxed in X. vesparum (2-tailed paired T-Test of ENC/MILC vectors (see methods); P < 1 × 10-5), and an apparent switch in threonine (T) preference from ACU to ACA has occurred. In Figure 6A, these levels are assessed in context, alongside the holometabolan dataset using the ENC and MILC measures of codon bias for individual PCGs (see methods; [43]). In X. vesparum, all 12 PCGs occupy a much more relaxed, spread apart region of ENC/MILC space (the higher the value, the less codon bias), with average gene ENC and MILC values of 36.26 and 0.71 respectively. In M. australiensis these are 29.63 and 0.56, and in Holometabola, they average 32.77 (ENC) and 0.62 (MILC) over all genes. It follows that the values of codon bias among PCGs in M. australiensis should be found in the densest cluster of background holometabolan PCGs. The thysanuran genomes (empty black circles), alongside three coleopteran genomes (forming the isolated cluster of empty blue circles with MILC values > 1) have particularly balanced codon usage. Although beetles generally follow this rule; Pyrophorus divergens, Tetraphalerus bruchi and Tribolium castaneum in particular, have unusually high MILC values. Interestingly, unlike the thysanurans, these do not result in correspondingly elevated ENC values, possibly because of differences in the way MILC and ENC estimate bias [43].

Figure 5
figure 5

Codon usage in M. australiensis (red) and X. vesparum (black). Codon frequencies are given (out of 1) for individual amino acids.

Full size image
Figure 6
figure 6

Evaluation of codon bias across Holometabola. A) Codon bias in M. australiensis (red filled circle) and X. vesparum (black filled circle) measured using ENC and MILC indices (described in methods), are compared against remaining Holometabola, whose data points are given as empty circles: the colour scheme follows taxonomic group as in Figure 3 and 4. Data points refer to individual PCGs. B) MILC values are compared against (G+C)3: which is G+C content at the 3rd codon position of each PCG. This serves as an approximate measure of background nucleotide composition across individual PCGs.

Full size image

These results are consistent with prevailing neutral mutational theories positing that genomic G+C content is the most significant factor in determining codon bias between organisms [44–47]. This may explain why X. vesparum has significantly relaxed ENC-MILC values, whose global G+C content is 5% higher than M. australiensis. For the thysanuran and beetle genomes, although G+C content in the 3rd codon position of individual PCGs is not necessarily elevated (Figure 6B), it is global G+C% content that appears to matter: it has been documented that codon bias can be accurately predicted from intergenic regions [44], and it is unsurprising that these genomes should also have the highest genome wide G+C content across the 68 Holometabola dataset (Additional file 2). It is noticeable that the coleopteran T. bruchi has a very A-skewed genome (Figure 4A), which could also partly explain some of the discrepancies appearing between the different methods.

Secondary structure of ribosomal and transfer RNAs

With the addition of M. australiensis, a comparative revision of strepsipteran mitochondrial rRNA secondary structure was possible. Figures 7 and 8 are complete models for the rrnL in M. australiensis and X. vesparum respectively. The addition of new comparative evidence (from M. australiensis and Apis mellifera [48]) since Carapelli et al. (2006, [28]) enables some improvement over the structural model for X. vesparum. These modifications are highlighted in blue. Regions in which sequence variability is still too high for good structural covariation are highlighted in red. Overall, M. australiensis predictions are largely consistent with current insect consensus models but substantial portions of secondary structure remain problematic in X. vesparum. For example, the helix at the base of domain IV in M. australiensis is supported by comparative evidence from Drosophila melanogaster and Apis mellifera. In X. vesparum, most of this primary sequence in this helix is highly divergent and covarying base-pairing cannot be identified. Where this helix terminates in the single-strand bulge however, the base pairing AUG-UAU is supported by M. australiensis (AAG-UUU) which provides additional corroborative support for the bps in D. melanogaster and A. mellifera, which have, AAG-UUU and AGG-UCU respectively.

Figure 7
figure 7

Predicted secondary structure of the rrnL in M. australiensis. Regions in red receive little comparative support and remain ambiguous. Canonical A-U and G-C base pairings are given as '-'. G-U, G-A, and other non-canonical bp combinations are denoted as, '•', '○' and '●' respectively. Predicted tertiary interactions are represented by a line. Missing data are coded as 'N'.

Full size image
Figure 8
figure 8

Revised secondary structural prediction of the rrnL in X. vesparum. Regions in red receive little comparative support and remain ambiguous; regions in blue are structural predictions which have improved with the inclusion of comparative data from M. australiensis. Structural annotations follow figure 7.

Full size image

rrnS secondary structural models for M. australiensis and X. vesparum are presented in Figure 9. Absence of sequence across domain I precludes a complete analysis, but partial models for domains II and III are generally conserved across Strepsiptera. Regions highlighted in red are more labile, and adequate consensus models for these structures are still lacking. The structural predictions for 18 tRNA genes are given in Figure 10. There is less scope for variation within extremely length-constrained tRNA genes. Despite this, the TψP stem is still too variable to be useful in a deep comparative framework: even between M. australiensis and X. vesparum, base-pairings demonstrate little covariation. The three remaining stem-loop structures are more conserved across Strepsiptera, and more useful in a comparative structural context. Within these regions, X. vesparum demonstrates more substantial modifications to typical insect tRNA structure than M. australiensis - consisting largely of single base or base-pair substitutions/compensations. For example, the proximal base-pairs of the acceptor stem of tRNA-Pro (P) in X. vesparum appears as UCAG-CUGA. In M. australiensis, Lepidoptera (Ochrogaster lunifer [41]) and Hymenoptera (Vanhornia eucnemidarum [39]) it is CAAA-UUUG. Similarly, the proximal DHU stem consists of a canonical A-U base-pairing in X. vesparum, where in M. australiensis and other insect groups it appears as a non-canonical G-U. Further, the base-pairing C-G, found across disparate insect lineages is substituted by A-U in the proximal stems of the tRNA-Lys (K) acceptor and the tRNA-Glu (E) DHU helices of X. vesparum.

Figure 9
figure 9

Secondary structural predictions for the rrnS in Strepsiptera. A) M. australiensis, B) X. vesparum. Regions in red receive little comparative support and remain ambiguous. Structural annotations follow figure 7.

Full size image
Figure 10
figure 10

Secondary structural predictions for M. australiensis tRNA genes. Pertinent aspects relating to comparative structural differences between X. vesparum [28] and other insect lineages [39, 41, 58] are discussed in the main text.

Full size image

Conclusions

The mitochondrial genome of X. vesparum displays a number of characteristics that are not shared by its distant relative M. australiensis. In the latter, only 2 tRNA translocations disrupt an otherwise ancestral insect genome organization, whereas in X. vesparum, 3 tRNA translocations and at least 2 duplications have occurred (Figure 2). Codon bias in X. vesparum is also significantly relaxed: M. australiensis occupies a much more typical region of ENC/MILC space (Figures 5 and 6). Further, rRNA secondary structural model predictions in M. australiensis are much more consistent with current insect consensus structures (Figures 7, 8, 9 and 10). In X. vesparum, increased variability at the primary sequence level precludes the assembly of well corroborated models in several parts of the rRNA genes. Both strepsipteran genomes are unusually C-skewed and A+T% rich. Interestingly, PCGs are noticeably truncated in both strepsipteran taxa, but especially so in X. vesparum which has lost nearly twice as many codons (154) as M. australiensis (81), equating to 4.3% and 2.2% fewer codons than the average holometabolan genome respectively.

These results provide useful insight into the evolution of strepsipteran mitochondria spanning the shift from a largely free-living insect into an extremely unique and simplified endoparasitoid, since the emergence of Strepsiptera over 100 million years ago [49]. The modifications to the mitochondrial genomes of X. vesparum and M. australiensis appear consistent with organisms that have evolved extremely derived lifestyles. M. australiensis represents a transitional phase in the evolutionarily specialization of Strepsiptera, and the composition, architecture and structure of its genome reflects this. These observations raise important questions about the changing selective environment in (mitochondrial) genomes belonging to small, and almost entirely host-dependent endoparasitoids. Future investigations into the evolutionary and developmental origins of this unique biological system must ensure that strepsipteran species from Mengenillidae are adequately sampled: wide sampling across this group is crucial if the underlying ancestral signal is to be maximized.

Methods

Specimen description

DNA was used from a small set of free-flying males collected in light traps, leaving little possibility for contamination through thysanuran host DNA. Nuclear and mitochondrial gene alignments from a widely sampled Strepsiptera taxon set (unpublished) confirm this. We also include 3 thysanuran genomes in Figure 8 to emphasize the genome composition differences between M. australiensis and its target host group. Specimens were collected on the 16th March 2006 in Australia, Queensland, Blackdown Tablelands NP, South Mimosa Creek, 50 m dstr road, 794 mao. Coordinates as follows: S23 47.687'E 149 04. 195'. Light trap loc 34 (Collectors were N. Jönsson, T. Malm & D. Williams). Complete species name: Mengenilla australiensis Kifune & Hirashima 1983.

Genome isolation and sequencing

M. australiensis specimens were extracted for PCR either by macerating tissue and digesting overnight at 50°C with Chelex 100 grade resin (Bio-Rad Cat. no 142-1253, Hemel Hampstead, UK) and enzyme Proteinase K (BIOLINE, Cat. no BIO-37037, London, UK) or by employing a QIAGEN blood and tissue column extraction protocol. Amplification of the M. australiensis mitochondrion was carried out by amplifying three ~3 kb fragments, with intervening sections that did not overlap being amplified by standard PCR. Long PCR mixes containing 1 μl BIO-X-ACT long DNA polymerase (BIOLINE, Cat. no. BIO-21049), 1 μl dNTPs, 25 μl 2× Polymate additive (BIOLINE, Cat. no. BIO-37041) for A+T rich sequences, 5 μl 10× opti-buffer, 4 μl MgCl2 solution, 3 μl forward/reverse primers and 5-10 μl of genomic DNA were made up to 50 μl reactions with ultrapure (NANOpure Diamond™) water. Long PCR cycling parameters were as follows: 92°C for 2 minutes; 10 cycles of: 92°C for 10 seconds, 53°C for 30 seconds and 68°C for 13 minute; 28 cycles of: 92°C for 15 seconds, 53°C for 30 seconds and 68°C for 14 minute; with an extension step of 7 minutes at 68°C. Long PCR fragments were gel extracted (QIAGEN Cat no. 28704, Hemel Hampstead, UK) and used as templates in subsequent sequencing reactions. Standard PCR mixes contained 0.5 μl Taq polymerase, 1 μl dNTPs, 5 μl MgCl2 free buffer, 4 μl MgCl2solution, 3 μl forward/reverse primers and 2 μl of genomic DNA, made up to 50 μl reactions with ultrapure (NANOpure Diamond™) water. Standard PCR cycling parameters were as follows: 94°C for 2 minutes, then 35 cycles of 94°C for 1 minute, T°C for 45 seconds (T varied with primer pair) and 72°C for 1 minute, with an extension step of 5 minutes at 72°C. For difficult regions, Phusion high fidelity Taq was utilized, following the manufacturer's instructions. Standard PCR products were cleaned using 2 μl of Shrimp Alkaline Phosphatase (SAP) and 3 μl of Exonuclease I (in a 1:10 SAP buffer dilution) per 50 μl PCR reaction, or fragments were extracted from agarose gels and sequenced directly. Problematic fragments that did not sequence directly were inserted into pGEM-T Easy plasmid vectors (Promega, Southampton, UK), and multiple clones were sequenced for each fragment. Sequencing reactions were performed using the BigDye® Terminator v3.1 cycle sequencing kit with the following modifications: 1 μl of BigDye®, 1.5 μl 5× Buffer, 1 μl primer (3.2 pmol), 1-3 μl genomic DNA made up to 10 μl reactions with ultrapure (NANOpure Diamond™) water. Sequence reads were generated using an Applied Biosystems 3730xl DNA Analyzer. Sequences were imported into GridinSoft Notepad (Lite Edition; http://notepad.gridinsoft.com) and BioEdit (version 7.0.5.3; [50]) for analysis. Primer pairs employed are given in Additional file 3. The annotated genome fragment can be found in Genbank under the accession number GU188852. Voucher specimens and extractions are deposited in the Hope Entomological Collections, Oxford University Museum of Natural History, Oxford, UK, and the Swedish Museum of Natural History, Stockholm, Sweden.

Genome charactarization

PCGs and their boundaries were diagnosed via sequence comparison across alignments derived from the 68 holometabolan genome dataset. Start codons found to be in-frame and not overlapping with upstream genes on the same strand were usually identified as 5' gene boundaries. Stop codons were typically identified as TAA, TAG, TA, or T. A complete annotation of the M. australiensis genome is given in Additional file 4. One non-coding fragment, 24 bp in length, was identified between the stop codon of nad1 and the beginning of tRNA-Ser2. Secondary structural predictions for ribosomal and transfer RNAs were ascertained by implementing a comparative structural framework [51–55] (for an online tutorial see the jRNA web site: http://hymenoptera.tamu.edu/rna/), in conjunction with a thermodynamic-based RNA folding algorithm (mfold; [56]). mfold gave initial raw estimates of structure: a tool that is especially useful for variable domains with highly divergent primary sequence. In this approach, structures with the lowest thermodynamic stability values were most highly considered (although structures with only marginal thermodynamic differences were also evaluated) and robust comparative evidence where available, took priority. The following ribosomal RNA structural templates from the Comparative RNA Web (CRW) site database [57] were used: Escherichia coli, Drosophila melanogaster, Drosophila virilis, Chorthippus parallelus and Apis mellifera [48]. For tRNA genes, the following models were used as comparative data points: X. vesparum (Strepsiptera: [28]), O. lunifer (Lepidoptera: [41]), V. eucnemidarum (Hymenoptera: [39]) and Cryptopygus antarcticus (Collembola: [58]), structures were verified for accuracy in tRNA-Scan [59]. Gene order was assessed across Holometabola using the inferred ancestral insect arrangement as a comparative standard [29, 30]. As a measure of gene length fluctuation, total and individual gene length was counted across strepsipteran PCGs, and compared across the dataset of 68 holometabolan genomes. To include nad2 gene, the 5' portion (approx. 100 amino acids) was excluded from analysis so that M. australiensis could be assessed alongside the entire dataset.

The 68 holometabola mitochondrial genome dataset was compiled from previously published data across the 6 major holometabolan divisions: Diptera [60–74], Lepidoptera [33, 75–81], Hymenoptera [39, 82–87], Coleoptera [88–93], Neuropterida [34, 94, 95] and Strepsiptera [28].

Nucleotide content and Codon usage

A+T% and G+C% values were calculated for the α strand of the M. australiensis genome, which runs clockwise in the direction of transcription from nad2 (Figure 1). Genomic A- and G-skew measures were calculated using total nucleotide % values in the following manner: [A-T]/[A+T] and [G-C]/[G+C] following Perna & Kocher [96]. A complete list of nucleotide content values for the 68 holometabolan dataset is given in Additional file 3. Codon usage was investigated by importing the complete cDNA datasets of 68 available holometabolan mitochondrial genomes, including X. vesparum and the sequenced genome of M. australiensis into INCA v2.1 [42]. Two independent measures of codon bias: ENC (Effective Number of Codons used) and MILC (Measure Independent of Length and Composition) are implemented, whose statistical justification and behaviour under simulation differ markedly [43]. These two independent measures were combined as a singe vector statistic as follows: √(ENC2 + MILC2). The shortest gene, atp8 was also excluded from analyses due to artificially increased ENC/MILC bias in genes shorter than around 100 amino acids. nad4L was of borderline length (approx. 90-105 amino acids) and was retained, although nad4L from Drosophila sechellia and Melipona bicolor produced arbitrary ENC 'cut-off' values of 61, because ENC behaves less stably at borderline gene lengths; these data points were removed. Readers are referred to [43] for a detailed explanation.

Abbreviations

mt:

mitochondrial

PCG:

protein-coding gene

A+T region:

the putative control region

α strand:

leading strand of transcription

ENC:

Effective Number of Codons used

MILC:

Measure Independent of Length and Composition

rrnL and rrnS:

large (16S) and small (12S) subunit ribosomal RNA (rRNA) gene

tRNA:

genes are denoted as single letter amino acid IUPAC-IUB abbreviations

atp6 and atp8:

ATP synthase subunits 6 and 8

cob :

apocytochrome b

cox1-3:

cytochrome c oxidase subunits 1-3

nad1-6 and nad4L:

NADH dehydrogenase subunits 1-6 and 4L

MRCA:

Most Recent Common Ancestor.

References

  1. Kathirithamby J: Host-parasitoid associations in Strepsiptera. Annu Rev Ent. 2009, 54: 227-249. 10.1146/annurev.ento.54.110807.090525.

    Article  CAS  Google Scholar 

  2. Crowson RA: The phylogeny of Coleoptera. Annu Rev Ent. 1960, 5: 111-134. 10.1146/annurev.en.05.010160.000551.

    Article  Google Scholar 

  3. Crowson RA: The Biology of the Coleoptera. 1981, Academic Press, New York

    Google Scholar 

  4. Whiting MF, Wheeler WC: Insect homeotic transformation. Nature. 1994, 368: 696-10.1038/368696a0.

    Article  Google Scholar 

  5. Kristensen NP: Forty years' insect phylogenetics systematics. Zool Beitr. 1995, 36: 83-124.

    Google Scholar 

  6. Whiting MF, Carpenter JC, Wheeler QD, Wheeler WC: The Strepsiptera Problem: Phylogeny of the Holometabolous Insect Orders Inferred From 18S and 28S Ribosomal DNA sequences and Morphology. Syst Biol. 1997, 46: 1-68.

    CAS  PubMed  Google Scholar 

  7. Huelsenbeck JP: Is the Felsenstein Zone a fly trap?. Syst Biol. 1997, 46: 69-74.

    Article  CAS  PubMed  Google Scholar 

  8. Huelsenbeck JP: Systematic bias in phylogenetic analysis: is the Sterpsiptera problem solved?. Syst Biol. 1998, 47: 519-537.

    CAS  PubMed  Google Scholar 

  9. Hwang UW, Kim W, Tautz D, Friedrich M: Molecular phylogenetics at the Felsenstein zone: approaching the Strepsiptera problem using 5.8S and 28S rDNA sequences. Mol Phyl Evol. 1998, 9: 470-480. 10.1006/mpev.1998.0518.

    Article  CAS  Google Scholar 

  10. Huelsenbeck JP: A Bayesian perspective of the Strepsiptera problem. Tidjschr Ent. 2001, 144: 165-178.

    Google Scholar 

  11. Whiting MF: Long branch distractions and the Strepsiptera. Syst Biol. 1998, 47: 134-138. 10.1080/106351598261076.

    Article  CAS  PubMed  Google Scholar 

  12. Hünefeld F, Beutel RG: The sperm pumps of Strepsiptera and Antliophora (Hexapoda). J Zool Syst Evol Res. 2005, 43: 297-306. 10.1111/j.1439-0469.2005.00327.x.

    Article  Google Scholar 

  13. Wiegmann BM, Trautwein MD, Kim JW, Cassel BK, Bertone MA, Winterton SL, Yeates DK: Single-copy nuclear genes resolve the phylogeny of the holometabolous insects. BMC Biology. 2009, 7: 34-10.1186/1741-7007-7-34.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Beutel RG, Pohl H: Endopterygote systematics - where do we stand and what is the goal (Hexapoda, Arthropoda)?. Syst Ent. 2006, 31: 202-219. 10.1111/j.1365-3113.2006.00341.x.

    Article  Google Scholar 

  15. Kinzelbach R: The systematic position of Strepsiptera (Insecta). Am Entomol. 1990, 35: 292-303.

    Article  Google Scholar 

  16. Pohl H, Beutel RG, Kinzelbach R: Protoxenidae fam. nov. (Insecta, Strepsiptera) from Baltic amber - a 'missing link' in Strepsipteran phylogeny. Zool Scr. 2005, 34: 57-69. 10.1111/j.1463-6409.2005.00173.x.

    Article  Google Scholar 

  17. Kinzelbach RK: Strepsiptera. Die Tierwelt Deutschlands. 1978, 65: 166-

    Google Scholar 

  18. Kathirithamby J: Review of the order Strepsiptera. Syst Ent. 1989, 14: 41-62. 10.1111/j.1365-3113.1989.tb00265.x.

    Article  Google Scholar 

  19. Kathirithamby J: Strepsiptera. The Insects of Australia: A textbook for students and research workers. Edited by: Naumann ID, Carne PB, Lawrence JF, Nielson ES, Spradberry JP, Taylor RW, Whitten MJ, Littlejohn MJ. 1993, Melbourne: Melbourne University Press, 684-695. 2

    Google Scholar 

  20. Pohl H, Beutel RG: The phylogeny of Strepsiptera. Cladistics. 2005, 21: 328-374. 10.1111/j.1096-0031.2005.00074.x.

    Article  Google Scholar 

  21. Pohl H, Beutel RG: The evolution of Strepsiptera. Zoology. 2007, 111: 318-338. 10.1016/j.zool.2007.06.008.

    Article  Google Scholar 

  22. Harrison RG: Animal mitochondrial DNA as a genetic marker in population and evolutionary biology. Trends Ecol Evol. 1989, 4: 6-11. 10.1016/0169-5347(89)90006-2.

    Article  CAS  PubMed  Google Scholar 

  23. Zhang D-X, Hewitt GM: Nuclear integrations: challenges for mitochondrial DNA markers. Trends Ecol Evol. 1996, 11: 247-251. 10.1016/0169-5347(96)10031-8.

    Article  CAS  PubMed  Google Scholar 

  24. Perna NT, Kocher TD: Mitochondrial DNA: molecular fossils in the nucleus. Curr Biol. 1996, 6: 128-129. 10.1016/S0960-9822(02)00441-4.

    Article  CAS  PubMed  Google Scholar 

  25. Pons J, Vogler AP: Complex Pattern of Coalescence and Fast Evolution of a Mitochondrial rRNA Pseudogene in a Recent Radiation of Tiger Beetles. Mol Biol Evol. 2005, 22: 991-1000. 10.1093/molbev/msi085.

    Article  CAS  PubMed  Google Scholar 

  26. Song H, Buhay JE, Whiting MF, Crandall KA: Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. P Natl Acad Sci USA. 2008, 105: 13486-13491. 10.1073/pnas.0803076105.

    Article  CAS  Google Scholar 

  27. Wolstenholme DR: Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992, 141: 173-216. 10.1016/S0074-7696(08)62066-5.

    Article  CAS  PubMed  Google Scholar 

  28. Carapelli A, Vannini L, Nardi F, Boore JL, Beani L, Dallai R, Frati F: The mitochondrial genome of the entomophagous endoparasite Xenos vesparum (Insecta: Strepsiptera). Gene. 2006, 376: 248-259. 10.1016/j.gene.2006.04.005.

    Article  CAS  PubMed  Google Scholar 

  29. Boore JL, Collins TM, Stanton D, Daehler LL, Brown WM: Deducing the pattern of arthropod phylogeny from mitochondrial DNA rearrangements. Nature. 1995, 376: 163-165. 10.1038/376163a0.

    Article  CAS  PubMed  Google Scholar 

  30. Boore JL, Lavrov D, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 393: 667-668. 10.1038/33577.

    Article  Google Scholar 

  31. 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. 10.1007/s00239-007-9042-8.

    Article  CAS  PubMed  Google Scholar 

  32. Dowton M, Cameron SL, Dowavic JI, Austin AD, Whiting MF: Characterization of 67 mitochondrial tRNA gene rearrangements in the Hymenoptera suggests that mitochondrial tRNA gene position is selectively neutral. Mol Biol Evol. 2009, 26: 1607-1617. 10.1093/molbev/msp072.

    Article  CAS  PubMed  Google Scholar 

  33. 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. 10.1016/j.gene.2007.10.023.

    Article  CAS  PubMed  Google Scholar 

  34. Cameron SL, Sullivan J, Song H, Miller KB, Whiting MF: A mitochondrial genome phylogeny of the Neuropterida (lace-wings, alderflies and snakeflies) and their relationship to the other holometabolous insect orders. Zool Scr. 2009, 38: 575-590. 10.1111/j.1463-6409.2009.00392.x.

    Article  Google Scholar 

  35. Rokas A, Carroll SB: Bushes in the Tree of Life. PLOS Biol. 2006, 4: 1899-1904. 10.1371/journal.pbio.0040352.

    Article  CAS  Google Scholar 

  36. Maeta Y, Takahashi K, Shimada N: Host body size as a factor determining the egg complement of Strepsiptera, an insect parasite. Int J Insect Morphol Embryol. 1998, 27: 27-37. 10.1016/S0020-7322(97)00033-0.

    Article  Google Scholar 

  37. Castro LR, Austin AD, Dowton M: Contrasting rates of mitochondrial molecular evolution in parasitic Diptera and Hymenoptera. Mol Biol Evol. 2002, 19: 1100-1113.

    Article  CAS  PubMed  Google Scholar 

  38. Eggleton P, Belshaw R: Insect parasitoids: an evolutionary overview. Philos Trans R Soc Lond. 1992, 337: 1-20. 10.1098/rstb.1992.0079.

    Article  Google Scholar 

  39. 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. 10.1139/G06-030.

    Article  CAS  PubMed  Google Scholar 

  40. Dowton M, Austin AD: Increased genetic diversity in mitochondrial genes is correlated with the evolution of parasitism in the hymenoptera. J Mol Evol. 1995, 41: 958-965. 10.1007/BF00173176.

    Article  CAS  PubMed  Google Scholar 

  41. Salvato P, Simonato M, Battisti A, Negrisolo E: The complete mitochondrial genome of the bag-shelter moth Ochrogaster lunifer (Lepidoptera, Notodontidae). BMC Genomics. 2008, 9: 331-10.1186/1471-2164-9-331.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Supek F, Vlahovicek K: INCA: synonymous codon usage analysis and clustering by means of self-organizing map. Bioinformatics. 2004, 20: 2329-2330. 10.1093/bioinformatics/bth238.

    Article  CAS  PubMed  Google Scholar 

  43. Supek F, Vlahovicek K: Comparison of codon usage measures and their applicability in prediction of microbial gene expressivity. BMC Bioinformatics. 2005, 6: 182-10.1186/1471-2105-6-182.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Chen SL, Lee W, Hottes AK, Shapiro L, McAdams HH: Codon usage between genomes is constrained by genome-wide mutational processes. Proc Natl Acad Sci USA. 2004, 101: 3480-3485. 10.1073/pnas.0307827100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Kanaya S, Kinouchi M, Abe T, Kudo Y, Yamada Y, Nishi T, Mori H, Ikemura T: Analysis of codon usage diversity of bacterial genes with a self-organizing map (SOM): characterization of horizontally transferred genes with emphasis on the E. coli O157 genome. Gene. 2001, 276: 89-99. 10.1016/S0378-1119(01)00673-4.

    Article  CAS  PubMed  Google Scholar 

  46. Knight RD, Freeland SJ, Landweber LF: A simple model based on mutation and selection explains trends in codon and amino-acid usage and GC composition within and across genomes. Genome Biol. 2001, 2: RESEARCH0010-

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Herschberg R, Petrov DA: Selection on codon bias. Annu Rev Genet. 2008, 42: 287-299. 10.1146/annurev.genet.42.110807.091442.

    Article  Google Scholar 

  48. Gillespie JJ, Johnston JS, Cannone JJ, Gutell RR: Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements. Insect Mol Biol. 2006, 15: 657-686. 10.1111/j.1365-2583.2006.00689.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Grimaldi D, Kathirithamby J, Schawaroch V: Strepsiptera and triungula in Cretaceous amber. Insect Syst Evol. 2005, 36: 1-20.

    Article  Google Scholar 

  50. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.

    CAS  Google Scholar 

  51. Kjer KM, Baldridge GD, Fallon AM: Mosquito large subunit ribosomal RNA: simultaneous alignment of primary and secondary structure. Biochim Biophys Acta. 1994, 1217: 147-155.

    Article  CAS  PubMed  Google Scholar 

  52. Kjer KM: Use of rRNA secondary structure in phylogenetic studies to identify homologous positions: an example of alignment and data presentation from the frogs. Mol Phylogenet Evol. 1995, 4: 314-330. 10.1006/mpev.1995.1028.

    Article  CAS  PubMed  Google Scholar 

  53. Gillespie JJ: Characterizing regions of ambiguous alignment caused by the expansion and contraction of hairpin-stem loops in ribosomal RNA molecules. Mol Phylogenet Evol. 2004, 33: 936-943. 10.1016/j.ympev.2004.08.004.

    Article  CAS  PubMed  Google Scholar 

  54. Gillespie JJ, McKenna CH, Yoder MJ, Gutell RR, Johnston JS, Kathirithamby J, Cognato AI: Assessing the odd secondary structural properties of nuclear small subunit ribosomal RNA sequences (18S) of the twisted-wing parasites (Insecta: Strepsiptera). Insect Mol Biol. 2005, 14: 625-643. 10.1111/j.1365-2583.2005.00591.x.

    Article  CAS  PubMed  Google Scholar 

  55. Gillespie JJ, Yoder MJ, Wharton RA: Predicted secondary structures for 28S and 18S rRNA from Ichneumonoidea (Insecta: Hymenoptera: Apocrita): Impact on sequence alignment and phylogeny estimation. J Mol Evol. 2005, 61: 114-137. 10.1007/s00239-004-0246-x.

    Article  CAS  PubMed  Google Scholar 

  56. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. Cannone JJ, Subramanian S, Schnare MN, Collett JR, D'Souza LM, Du Y, Feng B, Lin N, Madabusi LV, Müller KM, Pande N, Shang Z, Yu N, Gutell RR: The Comparative RNA Web (CRW) site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs. BMC Bioinformatics. 2002, 3: 2-10.1186/1471-2105-3-2. Correction: BMC Bioinformatics. 3:15

    Article  PubMed Central  PubMed  Google Scholar 

  58. Carapelli A, Comandi S, Convey P, Nardi F, Frati F: The complete mitochondrial genome of the Antarctic springtail Cryptopygus antarcticus (Hexapoda: Collembola). BMC Genomics. 2008, 9: 315-10.1186/1471-2164-9-315.

    Article  PubMed Central  PubMed  Google Scholar 

  59. 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. 10.1093/nar/25.5.955.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Lessinger AC, Junqueira ACM, Lemos TA, Kemper EL, da Silva FR, Vettore AL, Arruda P, Azeredo-Espin AML: The mitochondrial genome of the primary screwworm fly Cochliomyia hominivorax (Diptera: Calliphoridae). Insect Mol Biol. 2000, 9: 521-529. 10.1046/j.1365-2583.2000.00215.x.

    Article  CAS  PubMed  Google Scholar 

  61. Ballard JWO: Comparative genomics of mitochondrial DNA in Drosophila simulans. J Mol Evol. 2000, 51: 64-75.

    CAS  PubMed  Google Scholar 

  62. 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. 10.1046/j.1365-2583.2000.00165.x.

    Article  CAS  PubMed  Google Scholar 

  63. 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. 10.1111/j.1365-2583.1995.tb00032.x.

    Article  CAS  PubMed  Google Scholar 

  64. Oliveira MT, Barau JG, Junqueira ACM, Feijao PC, da Rosa AC, Abreu CF, Azeredo-Espin AML, Lessinger AC: Structure and evolution of the mitochondrial genomes of Haematobia irritans and Stomoxys calcitrans: The Muscidae (Diptera: Calyptratae) perspective. Mol Phylogenet Evol. 2008, 48: 850-857. 10.1016/j.ympev.2008.05.022.

    Article  CAS  PubMed  Google Scholar 

  65. Stevens JR, Westr H, Wall R: Mitochondrial genomes of the sheep blowfly, Lucilia sericata, and the secondary blowfly, Chrysomya megacephala. Med Vet Entomol. 2008, 22: 89-91. 10.1111/j.1365-2915.2008.00710.x.

    Article  CAS  PubMed  Google Scholar 

  66. Yu DJ, Xu L, Nardi F, Li JG, Zhang RJ: The complete nucleotide sequence of the mitochondrial genome of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae). Gene. 2007, 396: 66-74. 10.1016/j.gene.2007.02.023.

    Article  CAS  PubMed  Google Scholar 

  67. 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. 10.1046/j.1365-2583.2003.00445.x.

    Article  CAS  PubMed  Google Scholar 

  68. 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. 10.1111/j.1365-3113.2006.00355.x.

    Article  Google Scholar 

  69. 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. 10.1007/BF02099755.

    Article  CAS  PubMed  Google Scholar 

  70. Ballard JWO: Comparative genomics of mitochondrial DNA in members of the Drosophila melanogaster subgroup. J Mol Evol. 2000, 51: 48-63.

    CAS  PubMed  Google Scholar 

  71. 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. 10.1016/j.gene.2004.06.031.

    Article  CAS  PubMed  Google Scholar 

  72. 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. 10.1111/j.1365-2583.1993.tb00131.x.

    Article  CAS  PubMed  Google Scholar 

  73. 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. 10.1139/g93-141.

    Article  CAS  PubMed  Google Scholar 

  74. Matsumoto Y, Yanase T, Tsuda T, Noda H: Species-specific mitochondrial gene rearrangements in biting midges and vector species identification. Med Vet Entomol. 2009, 23: 47-55. 10.1111/j.1365-2915.2008.00789.x.

    Article  CAS  PubMed  Google Scholar 

  75. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. 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.

    Article  CAS  PubMed  Google Scholar 

  77. 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: 217-225. 10.1111/j.1365-2583.2006.00630.x.

    Article  CAS  PubMed  Google Scholar 

  78. Yang L, Wei ZJ, Hong GY, Jiang ST, Wen LP: The complete nucleotide sequence of the mitochondrial genome of Phthonandria atrilineata (Lepidoptera: Geometridae). Mol Biol Rep. 2009, 36: 1441-1449. 10.1007/s11033-008-9334-0.

    Article  CAS  PubMed  Google Scholar 

  79. 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. 10.1016/j.gene.2008.01.019.

    Article  CAS  PubMed  Google Scholar 

  80. Lee ES, Shin KS, Kim MS, Park H, Cho S, Kim CB: The mitochondrial genome of the smaller tea tortrix Adoxophyes honmai (Lepidoptera: Tortricidae). Gene. 2006, 373: 52-57. 10.1016/j.gene.2006.01.003.

    Article  CAS  PubMed  Google Scholar 

  81. Liu Y, Li Y, Pan M, Dai F, Zhu X, Lu C, Xiang Z: The complete mitochondrial genome of the Chinese oak silkmoth, Antheraea pernyi (Lepidoptera: Saturniidae). Acta Biochim Biophys Sin (Shanghai). 2008, 40: 693-703. 10.1111/j.1745-7270.2008.00449.x.

    Article  CAS  Google Scholar 

  82. Castro LR, Dowton M: The position of the Hymenoptera within the Holometabola as inferred from the mitochondrial genome of Perga condei (Hymenoptera, Symphyta, Pergidae). Perga condei. 2005, 34: 469-479.

    CAS  Google Scholar 

  83. Cameron SL, Dowton M, Castro LR, Ruberu K, Whiting MF, Austin AD, Diement K, Stevens J: Mitochondrial genome organization and phylogeny of two vespid wasps. Genome. 2008, 51: 800-808. 10.1139/G08-066.

    Article  CAS  PubMed  Google Scholar 

  84. 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. 10.1016/j.gene.2006.12.031.

    Article  CAS  PubMed  Google Scholar 

  85. Silvestre D, Dowton M, Arias MC: The mitochondrial genome of the stingless bee Melipona bicolor (Hymenoptera, Apidae, Meliponini): sequence, gene organization and a unique tRNA translocation event conserved across the tribe Meliponini. Genet Mol Biol. 2008, 31: 451-460. 10.1590/S1415-47572008000300010.

    Article  CAS  Google Scholar 

  86. Crozier RH, Crozier YC: The mitochondrial genome of the honeybee Apis mellifera: complete sequence and genome organization. Genetics. 1993, 133: 97-117.

    PubMed Central  CAS  PubMed  Google Scholar 

  87. Wei S, Tang P, Zheng L, Shi M, Chen X: The complete mitochondrial genome of Evania appendigaster (Hymenoptera: Evaniidae) has low A+T content and a long intergenic spacer between atp8 and atp6. Mol Biol Rep. 2009,

    Google Scholar 

  88. Arnoldi FGC, 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. 10.1016/j.gene.2007.07.035.

    Article  CAS  PubMed  Google Scholar 

  89. 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: 978-985. 10.1016/j.ympev.2004.03.009.

    Article  CAS  PubMed  Google Scholar 

  90. Friedrich M, Muqim N: Sequence and phylogenetic analysis of the complete mitochondrial genome of the flour beetle Tribolium castanaeum. Mol Phylogenet Evol. 2003, 26: 502-512. 10.1016/S1055-7903(02)00335-4.

    Article  CAS  PubMed  Google Scholar 

  91. Li X, Ogoh K, Ohba N, Liang X, Ohmiya Y: Mitochondrial genomes of two luminous beetles, Rhagophthalmus lufengensis and R. ohbai (Arthropoda, Insecta, Coleoptera). Gene. 2007, 392: 196-205. 10.1016/j.gene.2006.12.017.

    Article  CAS  PubMed  Google Scholar 

  92. 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: 513-526. 10.1016/S1055-7903(02)00421-9.

    Article  CAS  PubMed  Google Scholar 

  93. Sheffield NC, Song H, Cameron SL, Whiting MF: A comparative analysis of mitochondrial genomes in Coleoptera (Arthropoda: Insecta) and genome descriptions of six new beetles. Mol Biol Evol. 2008, 25: 2499-2509. 10.1093/molbev/msn198.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  94. Beckenbach AT, Stewart JB: Insect mitochondrial genomics 3: the complete mitochondrial genome sequences of representatives from two neuropteroid orders: a dobsonfly (order Megaloptera) and a giant lacewing and an owlfly (order Neuroptera). Genome. 2009, 52: 31-38. 10.1139/G08-098.

    Article  CAS  PubMed  Google Scholar 

  95. Hua J, Li M, Dong P, Xie Q, Bu W: The mitochondrial genome of Protohermes concolorus Yang et Yang 1988 (Insecta: Megaloptera: Corydalidae). Mol Biol Rep. 2009, 37: 1757-1765. 10.1007/s11033-008-9379-0.

    Article  Google Scholar 

  96. Perna NT, Kocher TD: Patterns of Nueleotide Composition at Fourfold Degenerate Sites of Animal Mitochondrial Genomes. J Mol Evol. 1995, 41: 353-358. 10.1007/BF01215182.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We owe gratitude to Niklas Jönsson, Tobias Malm and Dawn Williams for collection of material, to Clive Cook Director of the Queensland Parks and Wildlife Service - Northern Region for the Collecting Permit (# WITK03408905), and to B. Viklund (Swedish Museum of Natural History, Stockholm, Sweden) for donation of material. Thanks must also be extended to Peter Holland for useful discussions, the Leverhulme Trust (grant no. F/08 502G) for funding to JK, and the Elizabeth Hannah Jenkinson fund (University of Oxford) for an award to DPM.

Author information

Authors and Affiliations

  1. Department of Zoology, University of Oxford, The Tinbergen Building, South Parks Road, Oxford, OX1 3PS, UK

    Dino P McMahon, Alexander Hayward & Jeyaraney Kathirithamby

Authors
  1. Dino P McMahon
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Alexander Hayward
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jeyaraney Kathirithamby
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Dino P McMahon.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DPM and AH devised and completed laboratory work. DPM conducted analyses and early manuscript draughts. All authors contributed to the final version of the manuscript.

Electronic supplementary material

12864_2009_2487_MOESM1_ESM.XLS

Additional file 1: Coding gene lengths. Spreadsheet containing coding gene lengths (in amino acids) across the Holometabola dataset. (XLS 31 KB)

12864_2009_2487_MOESM2_ESM.XLS

Additional file 2: Nucleotide content. spreadsheet of genomic nucleotide content, A-skew and C-skew across the Holometabola dataset. (XLS 40 KB)

Additional file 3: Primers. List of primer pairs used in the study. (XLS 20 KB)

Additional file 4: Gene annotation. Detailed summary of genome fragment architecture. (XLS 26 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and Permissions

About this article

Cite this article

McMahon, D.P., Hayward, A. & Kathirithamby, J. The mitochondrial genome of the 'twisted-wing parasite' Mengenilla australiensis (Insecta, Strepsiptera): a comparative study. BMC Genomics 10, 603 (2009). https://doi.org/10.1186/1471-2164-10-603

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-603

Keywords

  • Mitochondrial Genome
  • Codon Usage
  • tRNA Gene
  • Codon Bias
  • Much Recent Common Ancestor

BMC Genomics

ISSN: 1471-2164

Contact us