- Research article
- Open Access
Mitochondrial genome evolution in fire ants (Hymenoptera: Formicidae)
© Gotzek et al; licensee BioMed Central Ltd. 2010
Received: 29 October 2009
Accepted: 7 October 2010
Published: 7 October 2010
Complete mitochondrial genome sequences have become important tools for the study of genome architecture, phylogeny, and molecular evolution. Despite the rapid increase in available mitogenomes, the taxonomic sampling often poorly reflects phylogenetic diversity and is often also biased to represent deeper (family-level) evolutionary relationships.
We present the first fully sequenced ant (Hymenoptera: Formicidae) mitochondrial genomes. We sampled four mitogenomes from three species of fire ants, genus Solenopsis, which represent various evolutionary depths. Overall, ant mitogenomes appear to be typical of hymenopteran mitogenomes, displaying a general A+T-bias. The Solenopsis mitogenomes are slightly more compact than other hymentoperan mitogenomes (~15.5 kb), retaining all protein coding genes, ribosomal, and transfer RNAs. We also present evidence of recombination between the mitogenomes of the two conspecific Solenopsis mitogenomes. Finally, we discuss potential ways to improve the estimation of phylogenies using complete mitochondrial genome sequences.
The ant mitogenome presents an important addition to the continued efforts in studying hymenopteran mitogenome architecture, evolution, and phylogenetics. We provide further evidence that the sampling across many taxonomic levels (including conspecifics and congeners) is useful and important to gain detailed insights into mitogenome evolution. We also discuss ways that may help improve the use of mitogenomes in phylogenetic analyses by accounting for non-stationary and non-homogeneous evolution among branches.
There has been a rapid proliferation of whole mitochondrial genomes (mitogenomes) sequenced in recent years, no doubt driven in part by the increasing speed and decreasing cost of sequencing technologies. Whole mitogenomes are increasingly used in phylogenetic studies [1–6] and in analyses of genome rearrangements [7–11], which can also be used for phylogenetic inference [8, 12–16].
Hence, we here present the first complete mitogenomes of ants (Hymentopera: Formicidae; four mitogenomes from three fire ant species) and discuss their evolution in light of the other hymenopteran mitochondrial genomes. The four mitogenomes presented here represent important contributions to the ever expanding dataset of complete hymenopteran mitochondrial genomes in that they represent a previously unsampled, highly diverse, and ecologically dominant vespoid family (Formicidae). Our mitogenome data for three fire ant species belonging to two species groups [19, 20] also provides insights into mitogenome evolution at lower taxonomic levels. Such data generally are rare , with the notable exceptions in insects of a study employing comparable sampling for Nasonia  and much more extensive sampling in Drosophila (38 mitogenomes; GenBank Sept. 21, 2009).
Results and Discussion
For the present study we sequenced four complete mitogenomes from three fire ant species. These included two mitogenomes from Solenopsis invicta Buren, 1972 (referred to hereafter as "PMS" and "VMS" mitogenomes), one mitogenome from the closely related species S. richteri Buren, 1972 ("richteri" mitogenome) and one mitogenome from a more distantly related fire ant species S. geminata (Fabricius, 1804) ("geminata" mitogenome). These first complete mitogenome data for ants expand taxon selection of one of the largest and most diverse hymenopteran superfamilies (Vespoidea).
Summary of Solenopsis mitogenomes.
Gene content and order
Comparison of various methods to identify tRNAs in the Solenopsis mitogenomes.
The location of tRNAs differs from the hypothetical ancestral hymenopteran mitogenome , which is also a typical feature of hymenopteran mitogenome architecture [11, 25]. However, with only three apparent translocations, the Solenopsis mitogenome architecture appears to be less derived than that of the other hymenopteran mitogenomes [11, 25–32]. Most mitogenome rearrangements in Hymenoptera appear to be selectively neutral and involve tRNA translocations around the coxII - atp8 junction [29, 30] and the nd3 - nd5 junction [11, 25]. The Solenopsis mitogenomes show no variation in the coxII - atp8 junction and a translocation of tRNA-N from the nd3 - nd5 junction to the SrRNA - nd2 junction. The tRNA-V translocation also moved to the SrRNA - nd2 junction. Various types of gene order rearrangements are generally recognized, differing by their location (local vs. distant) and whether they retain their original orientation (inverted vs. not inverted) and generally can be classified into one of several categories: local inversions, local translocations (gene shuffling), translocations, and remote inversions (an inverted translocation) .
Visual inspection of tRNAs across species reveals that there is considerable amount of variation due to point mutations and indels in almost all tRNAs. This variation results in differences in ability of the three tRNA identification methods to correctly identify tRNAs across species (Table 2), suggesting that congeneric comparative studies of tRNAs are fruitful endeavors to studying tRNA evolution and mitgenome architecture .
Nucleotide and codon bias
Results of tests of recombination for the translocation of 500 bp fragment at nt8,950-9,483.
Type of test
Consistent with this scenario, six additional methods of analysis showed statistically non-significant yet congruent signals of recombination in the same location. One likely reason for the comparatively weak support for the hypothesized recombination event is that the translocation could have occurred sufficiently long ago that the signal of the event has decayed. Our results are unlikely an experimental artifact, since the recombinant sequence was recovered from separate PCR (and sequencing) reactions in which the initial PCR reactions utilized different conserved primers flanking this region. Further, the sequence data were high quality, the sequences differ slightly from the other three sequences, and the sequence traces do not show any evidence of conflicting base calls.
The phylogenetic tree recovered from all four Solenopsis mitogenomes is identical to a tree generated previously using sequence data from a much shorter region (coxI/tRNA-L2/coxII) and demonstrates that the S. invicta mitogenomes are not monophyletic: the PMS invicta mitogenome is more closely related to richteri than the VMS invicta mitogenome  (Figure 4). This phylogenetic pattern is upheld throughout the complete mitochondrial alignment, except for the short fragment shown above to have been transferred during the recombination event (Figure 6).
Our phylogenetic analysis of all hymenopteran mitogenomes encountered similar problems described previously [6, 50–52] in that some expected relationships were not recovered, e.g. the sister taxa relationship between the Proctotrupomorpha (Nasonia and Vanhornia; sensu Rasnitsyn ) and the remainder of the Apocrita [50, 54] or the sister relationship of the Ichneumonoidea (represented by Ichneumonidae in this study) with the Aculeata (Vespomorpha sensu Rasnitsyn ) [6, 51, 52, 54]. Also, our phylogenetic analyses were unable to support the monophyly of the Vespoidea (ants and vespids; Figure 4) due to the placement of Solenopsis outside of the Vespidae + Apoidea clade with high branch support. While ants generally are considered to belong to the Vespoidea , some doubts have been raised regarding the monophyly of this superfamily [54, 56] and our results mirror those of Dowton et al.  using only the mitochondrial 16 S rRNA locus.
Our homotachous (site specific rates do not change along tree) Bayesian analysis recovered a topology fully concordant with that of Dowton et al.  (Figure 4B). This was not surprising, since we tried to follow the suggestions made by Dowton et al.  regarding the best analysis parameters for hymentoperan mitogenomic data (i.e., using nucleotide data, exclusion of 3rd codon positions, sampling several outgroups, and using a partitioned Bayesian framework). This suggests that the phylogenetic inference was robust to additional taxon sampling (Solenopsis, Diadegma, Evania). However, this topology differed from our other analyses in the placement of Ichneumonidae. Whereas the homotachous Bayesian analysis recovered the Ichneumonidae as a sister group to Schlettererius (Stephanidae), the other analyses placed the Ichneumonidae basal to an Aculeata + Evania clade. The analyses also differed in the placement of Primeuchroeus. Both maximum likelihood methods placed Primeuchroeus with Solenopsis, whereas both Bayesian analyses left its position unresolved at the base of the other Aculeata. However, all of these placements were poorly supported regardless of the branch support measure. Our analyses suggest that optimality criteria and models of nucleotide evolution have a stronger influence on the resulting phylogenetic inference in our dataset than taxon sampling. The importance of optimality criterion was previously also noted by Dowton et al.  and suggests that the extraction of phylogenetic signal from hymenopteran mitogenomic data is not simple.
While Dowton et al.  showed that the accuracy of mitogenomic phylogenetic analysis is greatly improved by the implementation of a particular partitioning scheme in a Bayesian framework, their analyses failed to account for heterotachous (site specific rate of evolution changes along the tree) and non-stationarity (base frequencies change along the tree) substitution processes between branches during phylogenetic analyses [57–60], of which the latter is an especially well-known feature of the hymenopteran mitogenomes [[6, 44, 61, 62] this study]. Hence, we compared the phylogenies derived under homotachous and stationary models of nucleotide substitution with those allowing site specific rates of change and base frequencies to change between branches. We implemented the covarion model  in MrBayes to accommodate heterotachy and the model of Galtier and Gouy  to account for non-stationarity (and non-homogeneity) using nhPhyML. Applying Galtier and Gouy's  model under the maximum likelihood criterion significantly improved the fit of the model to the data (nhPhyML: -lnL = 159,198.79453; PhyML: -lnL = 160,821.701853; df = 60, LRT = 3245.8125, p < 0.0001 ), but it did not result in a different topology (Figure 4A). Use of the covarion model (which accommodates heterotachy by allowing sites to change their substitution rate along branches by being switched on or off) in the Bayesian framework did not result in an improvement in log-likelihood (covarion model: -lnL = 156,475.84; stationary model: -lnL = 156,228.88), but recovered a different topology almost identical to that recovered under the maximum likelihood optimality criterion (Figure 4B). This decrease in log-likelihood of the covarion model was surprising [but see ], but the exact implementation of heterotachy may be important [60, 63] so the evaluation of other heterotachous models (e.g., mixed branch length model  or the RERV model ) may be warranted. More thorough studies are undoubtedly required to tease apart the contribution and possible interaction between different models and optimality criteria [4, 6, 60, 63, 66].
It is evident that reconstructing the hymenopteran phylogeny using only mitogenome data is problematic due to a variety of factors, e.g, differences in GC-content and rates of evolution among branches  or an ancient rapid diversification event [67, 68] which is known to greatly complicate phylogenetic inference . To account for these confounding factors, one should attempt to break up long branches by increased taxon sampling wherever possible. Including several species per genus (Nasonia, Bombus, and Solenopsis) as we have done is not very effective. Hence we suggest adding more genera and previously unsampled higher level taxa to more effectively break up long branches. While increased taxon sampling will likely rapidly occur in the near future, we agree with earlier suggestions of the importance of not relying on mitogenomes as the sole source of data for inferring phylogenetic relationships . Mitochondria essentially provide one very specific evolutionary history of a lineage since they are maternally inherited as single markers independently from the nuclear genome, and their utility can be further compromised if their transmission is influenced by intracellular symbionts such as Wolbachia . Given the rapidly decreasing costs in high-through-put sequencing and the various rapidly increasing genomic resources for several apocritan lineages (Apis and Nasonia genomes [71–73], Solenopsis, Camponotus, Vespula, and Microctonus genomic libraries [74–78]), we are hopeful that future phylogenetic work will be based on a broader genetic basis employing many more molecular characters .
To fully clarify the evolution of mitochondrial genomes in the Hymenoptera will require much greater sampling at all taxonomic levels (i.e., from the superfamily to the species and possibly intraspecific level), which should provide much needed data to fully address the patterns and rate of evolution and genome organization in this organelle. While this information will yield an increased understanding of hymenopteran evolutionary history, mitogenome studies need to be complemented by phylogenetic studies utilizing other sources of data against which we can compare and contrast the information yielded by whole mitochondrial genome analyses. Given the recent rate of publications with full mitogenomes, this will no doubt occur in the near future, yet could be greatly optimized through concerted sampling efforts by the greater scientific community.
Generation of sequence data
Complete mitochondrial genomes were generated for four individuals from three closely related Solenopsis species. The two S. invicta samples were collected in Pascagoula, MS (PMS) and Vera, MS (VMS), S. geminata was collected in Tallahassee, FL, and S. richteri was collected in Rosario, Argentina. Even though the two S. invicta mitogenomes were sequenced from specimens collected in the invasive range, the ~1 kb coxI/tRNA-L2/coxII section of these mitogenomes are identical to haplotypes from the native range (PMS is H22 [GenBank: AY249114]; VMS is H5 [GenBank: AY249097]) and represent two widely divergent clades within S. invicta and two geographically and genetically distinct populations (Argentina and Brazil) [49, 79]. We sequenced the entire mitogenome of each individual using a primer walking approach by performing 34-40 separate PCRs of genomic portions of the genome of variable size (~400-1,000 bp). We designed primers such that substantial overlap occurs between the various amplicons (allowing independent verification of sequence data by increasing depth of coverage) and such that the combination of all PCR amplicons spanned the entire genome. We initially attempted to amplify portions of the genome using conserved primers published in Simon et al.  as well as by designing conserved primers for a subset of coding genes by aligning available mitogenome sequences for Apis mellifera [GenBank: NC001566] and Drosophila yakuba [GenBank: NC001322]. We subsequently designed additional primers spanning other regions with no or low coverage using sequence data generated for fire ants. All primers developed and used for this study are presented in Additional file 1.
Several lines of evidence suggest that our sequences generated specifically represent mitochondrial genomic DNA rather than nuclear mitochondrial-like sequences (numts), which appear to be common and are generally short and highly fragmented in ants and other Hymenoptera [81–83]: The redundant yet independent PCR amplification of essentially every genomic region, the absence of ambiguous base calls characteristic of heterozygotes, the PCR amplification of the complete genome, and the generation of a contig identical in sequence to PMS using 454 pyrosequencing technology (DDS unpublished data).
All PCR amplicons were sequenced in both directions and each strand was assembled into single contigs with overlapping ends, indicating that our mitogenome sequences contained no gaps. Leading and lagging strand for each mitogenome were then aligned and manually checked for indels or ambiguous base calls. Mitogenomes were deposited in the NCBI GenBank database [GenBank: HQ215537, HQ215538, HQ215539, HQ215540].
Mitogenomes were annotated using the DOGMA webserver , which uses BLASTX against a custom database to identify protein coding genes. We verified all annotations made with DOGMA: coding regions were checked against a S. invicta EST database  and tRNAs were validated using ARWEN 1.2  and tRNAscan-SE 1.21  since DOGMA only uses COVE  to identify tRNAs. Generally, tRNAscan-SE has very low false positive rates and thus rarely mispredicts tRNAs (COVE scores ≥ 20 are usually considered reliable ), whereas ARWEN has a low false negative rate and usually identifies all tRNAs . Generally, DOGMA identified significantly more tRNAs than either ARWEN or tRNAscan-SE, sometimes with quite high COVE scores. Two tRNAs in particular were not recovered, tRNA-S1 and tRNA-N. These, however could be folded manually.
Nucleotide sequences were aligned based on amino acid alignments using MUSCLE 3.6 . Models of nucleotide evolution were estimated for protein coding genes using jModeltest . DnaSP 4.50.3  was used to estimate codon usage bias and nucleotide frequency bias [91–93]. The CODEML program in the PAML4.2 package  was used to test for site-specific evidence of positive selection while correcting for nucleotide bias . We employed the following parameters: runmode = 0, omega and kappa estimated (from three different starting points), empirical codon frequencies from each codon position (codonfreq = 2).
Following the recommendations of Posada , we employed a suite of recombination detection programs offered in the program packages TOPALi 2.5  and RPD 3b32  and the RecombiTest website  http://www.lifesci.sussex.ac.uk/CSE/test/index.php to test for recombination in the Solenopsis mitogenomes (see Table 3 for specific tests used). When any of the recombination tests only utilized three sequences at a time (e.g., RDP), analyses were repeated with every possible sequence triplet combination and p-values were Bonferroni corrected. All settings were left at the software default for the initial analyses, except for the PDM and LRT, where we used flexible window sizes. The highest acceptable p-value was 0.05 (unless Bonferroni corrected). Loosely following Tsaousis et al.'s  criteria for evidence of recombination, we consider as good evidence for recombination when more than one test detected a recombination event (although without regard to the test being a global or a local method). The more tests recovering evidence for recombination the more confident we are that it represents a true recombination event. Although this classification is admittedly arbitrary, we agree with White et al.  that identifying instances of recombination is inherently difficult and requires the heuristic use of several methods to identify potential recombinants.
Phylogenetic analyses were conducted on protein coding genes of the hymentoperan mitogenomes and 9 outgroups (3 flies [GenBank: X03240, AF260826, AJ242872], 3 beetles [GenBank: AJ312413, DQ768215, AB267275], and 3 moths [GenBank: AF442957, AF149768, AY242996]). jModeltest  was used to estimate the most appropriate model of nucleotide evolution for each codon position at each locus separately. Following the suggestion of Dowton et al.  we used the Bayesian approach using nucleotide sequences and implemented the GTR+I+Γ model of sequence evolution across genes and codon positions since jModeltest usually indentified this model as the best fitting for each data partition. MrBayes 3.1.2  was then used to recover phylogenetic hypotheses. All parameters were unlinked between partitions. Two independent analyses were run for three million generations, each with three heated and one cold chain. Parameters were sampled every 1000th generation. Convergence between runs was assessed when log-likelihoods had plateaued, PRSF factors were ~1, and split frequencies had dropped < 0.01. Samples taken prior to convergence were removed before samples were summarized. The same analysis was repeated implementing the covariotide model of sequence evolution to account for heterotachy (changes in site-specific evolutionary rates across lineages) , which has been shown to effectively accommodate heterotachy [, but see ]. Since this analysis took longer to converge, 5 million generations were run.
Maximum likelihood analyses were implemented on the PhyML 3.0 webserver  http://www.atgc-montpellier.fr/phyml/. We implemented the GTR+I+Γ model of nucleotide substitution on the unpartitioned dataset, estimated proportion of invariable sites and gamma shape parameter using six substitution rate categories, and optimized equilibrium frequencies, branch lengths, and tree topology (using the nearest-neighbor interchange [NNI] and sub-tree pruning and regrafting [SPR]) on five random starting trees. In addition to running one hundred bootstrap replicates to estimate levels of branch support, we also implemented the SH-like aLRT, which assesses the likelihood gain of the presence of that branch . To accommodate non-stationarity (changes of base frequencies between branches) we implemented nhPhyML-Discrete  using default options and the topology recovered from the heterotachous Bayesian analysis as the starting tree.
The evolution of tRNA-N was studied using phylogenetic analyses as suggested by Saks et al.  and Dowton and Austin , which were conducted using 100 bootstrap replicates in PhyML using the same configuration as described above. Other relevant hymenopteran tRNAs (D, N, and V) were downloaded from GenBank and aligned using MUSCLE. Unlike other authors [29, 37, 38], unpaired loops and anticodons were not removed following the suggestions of Wong et al. . However, we would like to point out that the phylogenetic analysis should only be interpreted as a heuristic tool, since the alignment of many very short, evolutionary very old, and highly AT-biased sequences is not trivial, regardless of alignment method used or prior editing to remove problematic areas.
We would like to thank Kenneth Ross (UGA) and Seán Brady (NMNH, Smithsonian Institution) for critically reviewing the manuscript. We also thank four anonymous reviewers for their comments which greatly improved the manuscript. This study was supported in part by a grant from the United States Department of Agriculture NRICGP (2006-35302-16561; Agreement No. 58-6615-7-137).
- Arnason U, Adegoke JA, Gullberg A, Harley EH, Janke A, Kullberg M: Mitogenomic relationships of placental mammals and molecular estimates of their divergences. Gene. 2008, 421: 37-51. 10.1016/j.gene.2008.05.024.PubMedGoogle Scholar
- Dong S, Kumazawa Y: Complete mitochondrial DNA sequences of six snakes: phylogenetic relationships and molecular evolution of genomic features. J Mol Evol. 2005, 61: 12-22. 10.1007/s00239-004-0190-9.PubMedGoogle Scholar
- Cameron SL, Lambkin CL, Barker SC, Whiting MF: A mitochondrial genome phylogeny of Diptera: whole genome sequence data accurately resolve relationships over broad timescales with high precision. Syst Entomol. 2007, 32: 40-59. 10.1111/j.1365-3113.2006.00355.x.Google Scholar
- Sheffield NC, Song H, Cameron SL, Whiting MF: Nonstationary evolution and compositional heterogeneity in beetle mitochondrial phylogenomics. Syst Biol. 2009, 58: 381-394. 10.1093/sysbio/syp037.PubMedGoogle Scholar
- Fenn JD, Song H, Cameron SL, Whiting MF: A preliminary mitochondrial genome phylogeny of Orthoptera (Insecta) and approaches to maximizing phylogenetic signal found within mitochondrial genome data. Mol Phylogenet Evol. 2008, 49: 59-68. 10.1016/j.ympev.2008.07.004.PubMedGoogle Scholar
- Dowton M, Cameron SL, Austin AD, Whiting MF: Phylogenetic approaches for the analysis of mitochondrial genome sequence data in the Hymenoptera - A lineage with both rapidly and slowly evolving mitochondrial sequences. Mol Phyl Evol. 2009, 52: 512-519. 10.1016/j.ympev.2009.04.001.Google Scholar
- Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 392: 667-668. 10.1038/33577.PubMedGoogle Scholar
- Dowton M: Relationships among the cyclostomebraconid (Hymenoptera: Braconidae) subfamilies inferred from a mitochondrial tRNA gene rearrangement. Mol Phylogenet Evol. 1999, 11: 283-287. 10.1006/mpev.1998.0580.PubMedGoogle Scholar
- Cameron SL, Johnson KP, Whiting MF: The mitochondrial genome of the screamer louse Bothriometopus (Phthiraptera: Ischnocera): effects of extensive gene rearrangements on the evolution of the genome. J Mol Evol. 2007, 65: 589-604. 10.1007/s00239-007-9042-8.PubMedGoogle Scholar
- Beckenbach AT, Joy JB: Evolution of the mitochondrial genomes of gall midge (Diptera: Cecidomyiidae): Rearrangements and several truncation of the tRNA genes. Genome Biol Evol. 2009, 1: 278-287. 10.1093/gbe/evp027.PubMed CentralPubMedGoogle Scholar
- Dowton M, Cameron SL, Dowavic JI, Austin AD, Whiting MJ: 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.PubMedGoogle Scholar
- Boore JL, Brown WM: Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998, 8: 668-674. 10.1016/S0959-437X(98)80035-X.PubMedGoogle Scholar
- Smith MJ, Arndt A, Gorski S, Fajber E: The phylogeny of echinoderm classes based on mitochondrial gene arrangements. J Mol Evol. 1993, 36: 545-54. 10.1007/BF00556359.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Kumazawa Y, Nishida M: Variations in mitochondrial tRNA gene organization of reptiles as phylogenetic markers. Mol Biol Evol. 1995, 12: 759-772.PubMedGoogle Scholar
- Curole JP, Kocher TD: Mitogenomics: digging deeper with complete mitochondrial genomes. Trends Ecol Evol. 1999, 14: 394-398. 10.1016/S0169-5347(99)01660-2.PubMedGoogle Scholar
- Grissell EE: Hymenopteran biodiversity: some alien notions. Am Entomol. 1999, 45: 235-244.Google Scholar
- Gissi C, Iannelli F, Pesole G: Evolution of the mitochondrial genome of Metazoa as exemplified by the comparison of congeneric species. Heredity. 2008, 101: 301-320. 10.1038/hdy.2008.62.PubMedGoogle Scholar
- Trager JC: A revision of the fire ants, Solenopsis geminata group (Hymenoptera: Formicidae: Myrmicinae). J New York Entomol Soc. 1991, 99: 141-198.Google Scholar
- Pitts J, McHugh Ross KG: Cladistic analysis of the fire ants of the Solenopsis saevissima species-group (Hymenoptera: Formicidae). Zool Scripta. 2005, 34: 493-505. 10.1111/j.1463-6409.2005.00203.x.Google Scholar
- Oliveira DC, Raychoudhury R, Lavrov DV, Werren JH: Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp Nasonia (Hymenoptera: Pteromalidae). Mol Biol Evol. 2008, 25: 2167-2180. 10.1093/molbev/msn159.PubMed CentralPubMedGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148: 929-936.PubMed CentralPubMedGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen A-MK: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 341-449.Google Scholar
- Castro LR, Dowton M: The position of the Hymenoptera within the Holometabola as inferred from the mitochondrial genome of Perga condei (Hymeoptera: Symphyta: Pergidae). Mol Phyl Evol. 2005, 14: 469-479. 10.1016/j.ympev.2004.11.005.Google Scholar
- 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.PubMedGoogle Scholar
- Wei S-J, Shi M, He J-H, Sharkey M, Chen X-X: The complete mitochondrial genome of Diadegma semiclausum (Hymenoptera: Ichneumonidae) indicates extensive independent evolutionary events. Genome. 2009, 52: 308-319. 10.1139/G09-008.PubMedGoogle Scholar
- Crozier RH, Crozier YC: The mitochondrial genome of the honeybee Apis mellifera - complete sequence and genome organization. Genetics. 1993, 133: 97-117.PubMed CentralPubMedGoogle Scholar
- Dowton M, Castro LR, Campbell SL, Bargon SD, Austin AD: Frequent mitochondrial gene rearrangements at the hymenopteran nad3-nad5 junction. J Mol Evol. 2003, 56: 517-526. 10.1007/s00239-002-2420-3.PubMedGoogle Scholar
- Dowton M, Austin AD: Evolutionary dynamics of a mitochondrial rearrangement "hot spot" in the Hymenoptera. Mol Biol Evol. 1999, 16: 298-309.PubMedGoogle Scholar
- Castro LR, Ruberu K, Dowton M: Mitochondrial genomes of Vanhornia eucnemidarum (Apocrita: Vanhorniidae) and Primeuchroeus spp. (Aculeata: Chrysididae): evidence of rearranged mitochondrial genomes within the Apocrita (Insecta: Hymenoptera). Genome. 2006, 49: 752-766. 10.1139/G06-030.PubMedGoogle Scholar
- Cha SY, Yoon HJ, Lee EM, Yoon MH, Hwang JS, Jin BR, Han YS, Kim I: The complete nucleotide sequence and gene organization of the mitochondrial genome of the bumblebee, Bombus ignitus (Hymenoptera: Apidae). Gene. 2007, 392: 206-220. 10.1016/j.gene.2006.12.031.PubMedGoogle Scholar
- 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.Google Scholar
- Cantatore P, Gadaleta MN, Roberti M, Saccne C, Wilson AC: Duplication and remoulding of tRNA genes during the evolutionary rearrangement of mitochondrial genomes. Nature. 1987, 329: 853-855. 10.1038/329853a0.PubMedGoogle Scholar
- Giegé R, Sissler M, Florentz C: Universal rules and idiosyncratic features in tRNA identity. Nucl Acid Res. 1998, 26: 5017-5035. 10.1093/nar/26.22.5017.Google Scholar
- Saks ME, Sampson JR, Abelson J: Evolution of a transfer RNA gene through a point mutation in the anticodon. Science. 1998, 279: 1665-1670. 10.1126/science.279.5357.1665.PubMedGoogle Scholar
- Tocchini-Valentini G, Saks ME, Abelson J: tRNA Leucine identity and recognition sets. J Mol Biol. 2000, 298: 779-793. 10.1006/jmbi.2000.3694.PubMedGoogle Scholar
- Higgs PG, Jameson D, Jow H, Rattray M: The evolution of tRNA-Leu genes in animal mitochondrial genomes. J Mol Evol. 2003, 57: 435-445. 10.1007/s00239-003-2494-6.PubMedGoogle Scholar
- Rawlings TA, Collins TM, Bieler R: Changing identities: tRNA duplication and remolding within animal mitochondrial genomes. Proc Nat'l Acad Sci USA. 2003, 100: 15700-15705. 10.1073/pnas.2535036100.Google Scholar
- Xue H, Tong KL, Marck C, Grosjean H, Wong JT: Transfer RNA paralogs: evidence for genetic code-amino acid biosynthesis coevolution and an archeal root of life. Gene. 2003, 310: 59-66. 10.1016/S0378-1119(03)00552-3.PubMedGoogle Scholar
- Lavrov DV, Lang BF: Transfer RNA recruitment in mitochondrial DNA. Trends Genet. 2005, 21: 129-133. 10.1016/j.tig.2005.01.004.PubMedGoogle Scholar
- Tang DTP, Glazov EA, McWilliam SM, Barris WC: Dalrymple: Analysis of the complement and molecular evolution of tRNA genes in cow. BMC Genomics. 2009, 10: 188-10.1186/1471-2164-10-188.PubMed CentralPubMedGoogle Scholar
- Rogers HH, Bergman CM, Griffiths-Jones S: The evolution of tRNA genes in Drosophila. Genome Biol Evol. 2010, 2: 467-477. 10.1093/gbe/evq034.PubMed CentralPubMedGoogle Scholar
- Seligmann H: Undetected antisense tRNAs in mitochondrial genomes?. Biology Direct. 2010, 5: 39-10.1186/1745-6150-5-39.PubMed CentralPubMedGoogle Scholar
- Wei S-J, Tang P, Zheng L-H, Shi M, Chen X-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. 2010, 37: 1931-1942. 10.1007/s11033-009-9640-1.PubMed CentralPubMedGoogle Scholar
- Tsaousis AD, Martin DP, Ladoukakis ED, Posada D, Zouros E: Widespread recombination in published animal mtDNA sequences. Mol Biol Evol. 2005, 22: 925-933. 10.1093/molbev/msi084.PubMedGoogle Scholar
- Piganeau G, Gardner M, Eyre-Walker A: A broad survey of recombination in animal mitochondria. Mol Biol Evol. 2004, 21: 2319-2325. 10.1093/molbev/msh244.PubMedGoogle Scholar
- Ujvari B, M Dowton, T. Madsen: Mitochondrial DNA recombination in a free-ranging Australian lizard. Biol Lett. 2007, 3: 189-192. 10.1098/rsbl.2006.0587.PubMed CentralPubMedGoogle Scholar
- Posada D: Evaluation of methods for detecting recombination from DNA sequences: empirical data. Mol Biol Evol. 2002, 19: 708-717.PubMedGoogle Scholar
- Shoemaker DD, Ahrens ME, Ross KG: Molecular phylogeny of fire ants of the Solenopsis saevissima species-group based on mtDNA sequences. Mol Phylogenet Evol. 2006, 38: 200-215. 10.1016/j.ympev.2005.07.014.PubMedGoogle Scholar
- Castro LR, Dowton M: Molecular analyses of the Apocrita (Insecta: Hymenoptera) suggests that the Chalcidoidea are sister to the diaprioid complex. Invert Syst. 2006, 20: 603-614. 10.1071/IS06002.Google Scholar
- Dowton M, Austin AD: Molecular phylogeny of the insect order Hymenoptera: Apocritan relationships. Proc Natl Acad Sci USA. 1994, 91: 9911-9915. 10.1073/pnas.91.21.9911.PubMed CentralPubMedGoogle Scholar
- Dowton M, Austin AD, Dillon N: Bartowsky: Molecular phylogeny of the apocritan wasps: the Proctotrupomorpha and Evaniamorpha. Syst Entomol. 1997, 22: 245-255. 10.1046/j.1365-3113.1997.d01-42.x.Google Scholar
- Rasnitsyn AP: An outline of evolution of the hymenopterous insects (Order Vespida). Oriental Insects. 1988, 22: 115-145.Google Scholar
- Sharkey MJ: Phylogeny and classification of Hymenoptera. Zootaxa. 2007, 1668: 521-548.Google Scholar
- Brothers DJ: Phylogeny and evolution of wasps, ants and bees (Hymenoptera, Chrysidoidea, Vespoidea, and Apoidea). Zool Scripta. 1999, 28: 233-249. 10.1046/j.1463-6409.1999.00003.x.Google Scholar
- Pilgrim EM, von Dohlen CD, Pitts JP: Molecular phylogenetics of Vespoidea indicate paraphyly of the superfamily and novel relationships of its component families and subfamilies. Zool Scripta. 2008, 37: 539-560. 10.1111/j.1463-6409.2008.00340.x.Google Scholar
- Galtier N, Gouy M: Inferring pattern and process: maximum-likelihood implementation of a nonhomogeneous model of DNA sequence evolution for phylogenetic analysis. Mol Biol Evol. 1998, 15: 871-879.PubMedGoogle Scholar
- Boussau B, Gouy M: Efficient likelihood computations with non-reversible models of evolution. Syst Biol. 2006, 55: 756-768. 10.1080/10635150600975218.PubMedGoogle Scholar
- Huelsenbeck JP: Testing a covariotide model of DNA substitution. Mol Biol Evol. 2002, 19: 698-707.PubMedGoogle Scholar
- Kolaczkowski B, Thornton JW: A mixed branch length model of heterotachy improves phylogenetic accuracy. Mol Biol Evol. 2008, 25: 1054-1066. 10.1093/molbev/msn042.PubMed CentralPubMedGoogle Scholar
- Castro LR, Austin AD, Dowton M: Contrasting rates of mitochondrial molecular evolution in parasitic Diptera and Hymenoptera. Mol Biol Evol. 2002, 19: 1100-1113.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Wang H-C, Specer M, Susko E, Roger AJ: Testing for covarion-like evolution in protein sequences. Mol Biol Evol. 2006, 24: 294-305. 10.1093/molbev/msl155.PubMedGoogle Scholar
- Kolaczkowski B, Thornton JW: Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous. Nature. 2004, 431: 980-984. 10.1038/nature02917.PubMedGoogle Scholar
- Wu J, Susko E, Roger AJ: An independent heterotachy model and its implications for phylogeny and divergence time estimation. Mol Phyl Evol. 2008, 46: 801-806. 10.1016/j.ympev.2007.06.020.Google Scholar
- Zhou Y, Rodrigue N, Lartillot N, Philippe H: Evaluation of the model handling heterotachy in phylogenetic inference. BMC Evol Biol. 2007, 7: 206-10.1186/1471-2148-7-206.PubMed CentralPubMedGoogle Scholar
- Grimaldi D, Engel MS: Evolution of the Insects. 2005, Cambridge: Cambridge University PressGoogle Scholar
- Brady SG, Larkin L, Danforth BN: Bees, ants, and stinging wasps (Aculeata). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford: Oxford University Press, 264-269.Google Scholar
- Whitfield JB, Kjer KM: Ancient rapid radiations of insects: Challenges for phylogenetic analysis. Ann Rev Entomol. 2008, 53: 449-472. 10.1146/annurev.ento.53.103106.093304.Google Scholar
- Hurst GDD, Jiggins FM: Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. Proc R Soc B. 2005, 272: 1525-1534. 10.1098/rspb.2005.3056.PubMed CentralPubMedGoogle Scholar
- Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443: 931-49. 10.1038/nature05260.Google Scholar
- Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 2006, 16: 1334-1338. 10.1101/gr.5204306.PubMed CentralPubMedGoogle Scholar
- Nasonia Genome Working Group: Functional and Evolutionary Insights from the Genomes of Three Parasitoid Nasonia Species. Science. 2010, 327: 343-348. 10.1126/science.1178028.Google Scholar
- Wilfert L, Muñoz-Torres M, Reber-Funk C, Schmid-Hempel R, Tomkins J, Gadau J, Schmid-Hempel P: Construction and characterization of a BAC-library for a key pollinator, the bumblebee Bombus terrestris L. Insect Soc. 2009, 56: 44 - 48-10.1007/s00040-008-1034-1.Google Scholar
- Wang J, Jemielity S, Uva P, Wurm Y, Gräff J, Keller L: An annotated cDNA library and microarray for large-scale gene-expression studies in the ant, Solenopsis invicta. Genome Biol. 2007, 8: R9-10.1186/gb-2007-8-1-r9.PubMed CentralPubMedGoogle Scholar
- Goodisman MAD, Isoe J, Wheeler DE, Wells MA: Evolution of insect metamorphosis: A microarray-based study of larval and adult gene expression in the ant Camponotus festinatus. Evolution. 2005, 59: 858-870.PubMedGoogle Scholar
- Hoffman EA, Goodisman MAD: Gene expression and the evolution of phenotypic diversity in social wasps. BMC Biology. 2007, 5: 23-10.1186/1741-7007-5-23.PubMed CentralPubMedGoogle Scholar
- Crawford AM, Brauning R, Smolenski G, Ferguson C, Barton D, Wheeler TT, McCulloch A: The constituents of Microctonus sp. parasitoid venoms. Insect Mol Biol. 2008, 17: 313-324. 10.1111/j.1365-2583.2008.00802.x.PubMedGoogle Scholar
- Ross KG, Krieger MJB, Keller L, Shoemaker DD: Genetic variation and structure in native populations of the fire ant Solenopsis invicta: evolutionary and demographic implications. Biol J Linn Soc. 2007, 92: 541-560. 10.1111/j.1095-8312.2007.00853.x.Google Scholar
- Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994, 87: 651-701.Google Scholar
- Martins J, Solomon SE, Mikheyev AS, Mueller UG, Ortiz A, Bacci M: Nuclear mitochondrial-like sequences in ants: evidence from Atta cephalotes (Formicidae: Attini). Insect Mol Biol. 2007, 16: 777-784. 10.1111/j.1365-2583.2007.00771.x.PubMedGoogle Scholar
- Behura SK: Analysis of nuclear copies of mitochondrial sequences in honeybee (Apis mellifera) genome. Mol Evol Biol. 2007, 24: 1492-1505. 10.1093/molbev/msm068.Google Scholar
- Pamilo P, Viljakainen L, Vihavainen A: Exceptionally High Density of NUMTs in the Honeybee Genome. Mol Biol Evol. 2007, 24: 1340-1346. 10.1093/molbev/msm055.PubMedGoogle Scholar
- Wyman SK, Jansen RK, Boore JL: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20: 3252-3255. 10.1093/bioinformatics/bth352.PubMedGoogle Scholar
- Laslett D, Canbäck B: ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bioinformatics. 2008, 24: 172-175. 10.1093/bioinformatics/btm573.PubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964. 10.1093/nar/25.5.955.PubMed CentralPubMedGoogle Scholar
- Eddy SR, Durbin R: RNA sequence analysis using covariance models. Nucleic Acids Res. 1994, 22: 2079-2088. 10.1093/nar/22.11.2079.PubMed CentralPubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-97. 10.1093/nar/gkh340.PubMed CentralPubMedGoogle Scholar
- Posada D: jModeltest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedGoogle Scholar
- Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.PubMedGoogle Scholar
- Morton BR: Chloroplast DNA codon use: Evidence for selection at the psb A locus based on tRNA availability. J Mol Evol. 1993, 37: 273-280. 10.1007/BF00175504.PubMedGoogle Scholar
- Wright F: The "effective number of codons" used in a gene. Gene. 1990, 87: 23-29. 10.1016/0378-1119(90)90491-9.PubMedGoogle Scholar
- Shields DC, Sharp PM, Higgins DG, Wright F: "Silent"sites in Drosophila genes are not neutral: Evidence of selection among synonymous codons. Mol Biol Evol. 1988, 5: 704-716.PubMedGoogle Scholar
- Yang Z: PAML 4: a program package for phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedGoogle Scholar
- Albu M, Min XJ, Hickey D, Golding B: Uncorrected nucleotide bias in mtDNA can mimic the effects of positive Darwinian selection. Mol Biol Evol. 2008, 25: 2521-2524. 10.1093/molbev/msn224.PubMed CentralPubMedGoogle Scholar
- Milne I, Lindner D, Bayer M, Husmeier D, McGuire G, Marshall DF, Wright F: TOPALi v2: a rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops. Bioinformatics. 2009, 25: 126-127. 10.1093/bioinformatics/btn575.PubMed CentralPubMedGoogle Scholar
- Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics. 2005, 21: 260-262. 10.1093/bioinformatics/bth490.PubMedGoogle Scholar
- White DJ, Wolff JN, Pierson M, Gemmell NJ: Revealing the hidden complexities of mtDNA inheritance. Mol Ecol. 2008, 17: 4925-4942. 10.1111/j.1365-294X.2008.03982.x.PubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.PubMedGoogle Scholar
- Anisimova M, Gascuel O: Approximate Likelihood-Ratio Test for Branches: A Fast, Accurate, and Powerful Alternative. Syst Biol. 2006, 55: 539-552. 10.1080/10635150600755453.PubMedGoogle Scholar
- Wong KM, Suchard MA, Huelsenbeck JP: Alignment uncertainty and genomic analysis. Science. 2008, 319: 473-476. 10.1126/science.1151532.PubMedGoogle Scholar
- Martin D, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000, 16: 562-563. 10.1093/bioinformatics/16.6.562.PubMedGoogle Scholar
- Padidam M, Sawyer S, Fauquet CM: Possible emergence of new gemini viruses by frequent recombination. Virology. 1999, 265: 218-225. 10.1006/viro.1999.0056.PubMedGoogle Scholar
- Martin DP, Posada D, Crandall KA, Williamson C: A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses. 2005, 21: 98-102. 10.1089/aid.2005.21.98.PubMedGoogle Scholar
- Maynard Smith J: Analyzing the mosaic structure of genes. J Mol Evol. 1992, 34: 126-129.Google Scholar
- Posada D, Crandall KA: Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. Proc Natl Acad Sci USA. 2001, 98: 13757-13762. 10.1073/pnas.241370698.PubMed CentralPubMedGoogle Scholar
- Gibbs MJ, Armstrong JS, Gibbs AJ: Sister-Scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000, 16: 573-582. 10.1093/bioinformatics/16.7.573.PubMedGoogle Scholar
- Boni MF, Posada D, Feldman MW: An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007, 176: 1035-1047. 10.1534/genetics.106.068874.PubMed CentralPubMedGoogle Scholar
- Weiller GF: Phylogenetic profiles: a graphical method for detecting genetic recombinations in homologous sequences. Mol Biol Evol. 1998, 15: 326-335.PubMedGoogle Scholar
- Holmes EC, Worobey M, Rambaut A: Phylogenetic evidence for recombination in dengue virus. Mol Biol Evol. 1999, 16: 405-409.PubMedGoogle Scholar
- McGuire G, Wright F: TOPAL 2.0: Improved detection of mosaic sequences within multiple alignments. Bioinformatics. 2000, 16: 130-134. 10.1093/bioinformatics/16.2.130.PubMedGoogle Scholar
- Husmeier D, Wright F: Probabilistic Divergence Measures for Detecting Interspecies Recombination. Bioinformatics. 2001, 17 (Suppl 1): S123-S131.PubMedGoogle Scholar
- Husmeier D, McGuire G: Detecting Recombination in 4-Taxa DNA Sequence Alignments with Bayesian Hidden Markov Models and Markov Chain Monte Carlo. Mol Biol Evol. 2003, 20: 315-337. 10.1093/molbev/msg039.PubMedGoogle Scholar
- Sawyer SA: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6: 526-538.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.