- Research article
- Open Access
Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera:Culicidae), members of the coronator group
© Demari-Silva et al. 2015
- Received: 21 April 2015
- Accepted: 23 September 2015
- Published: 21 October 2015
The Coronator Group currently encompasses six morphologically similar species (Culex camposi Dyar, Culex coronator Dyar and Knab, Culex covagarciai Forattini, Culex usquatus Dyar, Culex usquatissimus Dyar, and Culex ousqua Dyar). Culex coronator has been incriminated as a potential vector of West Nile Virus (WNV), Saint Louis Encephalitis Virus (SLEV), and Venezuelan Equine Encephalitis Virus (VEEV). The complete mitochondrial genome of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi was sequenced, annotated, and analyzed to provide genetic information about these species.
The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi varied from 15,573 base pairs in Cx. usquatus to 15,576 in Cx. coronator. They contained 37 genes (13 protein-encoding genes, 2 rRNA genes, and 22 tRNA genes) and the AT-rich control region. Comparative analyses of the 37 genes demonstrated the mitochondrial genomes to be composed of variable and conserved genes. Despite the small size, the ATP8, ATP6 plus NADH5 protein-encoding genes were polymorphic, whereas tRNAs and rRNAs were conserved. The control region contained some poly-T stretch. The Bayesian phylogenetic tree corroborated that both the Coronator Group and the Culex pipens complex are monophyletic taxa.
The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx. usquatissimus and Cx. camposi share the same gene composition and arrangement features that match to those reported for most Culicidae species. They are composed of the same 37 genes and the AT-rich control region, which contains poly-T stretches that may be involved in the functional role of the mitochondrial genome. Taken together, results of the dN/dS ratios, the sliding window analyses and the Bayesian phylogenetic analyses suggest that ATP6, ATP8 and NADH5 are promising genes to be employed in phylogenetic studies involving species of the Coronator Group, and probably other species groups of the subgenus Culex. Bayesian topology corroborated the morphological hypothesis of the Coronator Group as monophyletic lineage within the subgenus Culex.
- Coronator Group
- Mitochondrial genome
Culex camposi Dyar, Culex coronator Dyar and Knab, Culex covagarciai Forattini, Culex ousqua Dyar, Culex usquatus Dyar and Culex usquatissimus Dyar are members of the Coronator Group of the subgenus Culex . Morphological features of the fourth-instar larva, pupa and female cannot distinguish among species of the Coronator Group. However, traits of the male genitalia allow species identification. Species of the Coronator Group are largely distributed and some are sympatric in the Neotropics . In southeastern Brazil, for instance, the geographical distribution of Cx. usquatus, Cx. coronator and Cx. camposi overlaps . Phylogenetic relationships among the species remain unresolved, because results of studies employing DNA sequences of the cytochrome c oxidase subunit 1 (COX1) mitochondrial gene [3, 4] and of the internal transcribed spacer 2 (ITS2) of the ribosomal DNA  showed poorly supported clades. Likely because of the lack of studies, there is no evidence of hybridization among species of the Coronator Group.
Culex coronator is a potential vector of arboviruses. As such, this species may participate in the dynamics of the transmission of the Saint Louis Encephalitis Virus (VESL), the Venezuelan Equine Encephalitis Virus (VEEV), the Mucambo Virus (MV) [6, 7], and the West Nile Virus (WNV) . Experiments of susceptibility of Cx. coronator to the WNV showed that the population from Florida, USA, is competent in disseminating the WN-FL03p2-3 strain under some conditions; however, the transmission rate decreases at lower temperatures (0–17 % at 25 °C and 28–67 % at 28 °C) . The impossibility of accurately identifying a species can lead to inconclusive studies regarding the vector status of it. The problematic identification of species of the Coronator Group may be causing difficulties in defining the medical importance of this group.
The mitochondrial genome of the majority of the metazoan organisms is a small genome (ranging from 15 to 20 kb). It contains 37 genes, which encode 13 Protein-Coding Genes (PCGs) involved in oxidative phosphorylation, 2rRNA (rrnL and rrnS) and 22 tRNAs genes, which are necessary for the translation of the proteins encoded by the genes [9, 10]. The mitochondrial genome also encompasses a region that is rich in adenine (A) and thymine (T). This AT-rich region starts and controls the replication of the mitochondrial genome . Because of the lack of recombination and fast evolutionary processes in comparison to those that involve nuclear genes, the mitochondrial genome seems to be a source of information that can be employed to address evolutionary processes of both vertebrates and invertebrates [12–14].
Although being largely employed in molecular taxonomy and phylogeny, only Culex quinquefasciatus Say and Culex pipiens L. possess their mitochondrial genome available in the Genbank. Considering the potential taxonomic and phylogenetic content of genes of the mitochondrial genome for Culex species and the role of Cx. coronator as vector of arboviruses, four species of the Coronator Group were the focus of this study. Next generation sequencing technology was applied to sequence the mitochondrial genomes of Culex coronator, Culex usquatus, Culex camposi, and Culex usquatissimus. The objectives of the study were to: 1) describe the mitochondrial genome of four species of the Coronator Group; 2) to analyze the mitochondrial genome of four species of the Coronator Group; 3) to compare the genomes of species of the Coronator Group with other species of the Culicidae family available in the Genbank; 4) to assess the monophyly of the Coronator Group.
Mosquito sampling and DNA extraction
Species, specimen codes, gender, states and geographical coordinates of the collection of localities in Brazil
Rio Grande do Sul
Mato Grosso do Sul
PCR amplification and sequencing
Genomic DNA was extracted from each mosquito individually using QIAgen Dneasy® Blood and Tissue kit (QIAgen Ltd., Crawley, UK) following the same procedures adopted by Foster et al. . The mitochondrial genome of each individual was amplified in two PCRs. A long PCR amplification generated a fragment of ~12,000 base pairs (bp) and a second reaction targeted a fragment of ~4,500 bp. Two sets of primers were used as follows: 16Sa (5′ CGCCTGTTTATCAAAAACAT 3′) and LCO1490 (5′ GGTCAACAAATCATAAAGATATTGG 3′) to amplify ~12,000-bp fragment, and 16Sb (5′ CCGGTTGAACTCAGATCA 3′) and HCO2198 (5′ TAAACTTCAGGGTGACCAAAAAATCA 3′) to generate fragments of ~4,500-bp. See Folmer et al. 1994  for additional information about the LCO1490 and HCO 2198 primers, and Kambhampati and Smith (1995)  for 16Sa and 16Sb primers. Each PCR reaction was performed using GoTaq® Long PCR Master Mix (Promega, Wisconsin, USA) following the protocols proposed by the manufacturer. The long PCR products were purified in DNA Clean & Concentrator™ (Zymo Research, California, USA) and quantified with Qubit (LifeTechnologies, Oregon, USA) for library preparation. Individually barcoded PCR libraries were created for the 12,000 bp and 4,500 bp PCR products employing Nextera® XT DNA Sample Preparation Kit (Illumina, Illinois, USA) and sequenced on the Illumina MiSeq platform with paired-end 250 base pairs chemistry.
Genome assembly and annotation
Mira v4 software was used for sequence assembly and Tablet to view the assembly [18, 19]. The MITOS website  was used to delimit the genes. Annotation of protein coding genes was conducted using the Wise2 v2.2 package and the HMMER2 models  (HMMER v3.1b2; http://hmmer.janelia.org/).
Nucleotide composition values were calculated using the MEGA 6.06 software . The GC—and AT-skews were used to determine strand asymmetry . To measure these values, the following formulas were used: GC-skew (G-C) / (G + C) and AT-skew (A-T) / (A + T). The secondary structure of tRNAs was inferred in tRNAscan-SE 1.21 . A sliding window of 200 bp and steps of 25 bp was implemented to estimate nucleotide divergence. The analysis was performed using the five complete mitochondrial genome sequences in DnaSP v. 5.10.01 . Ratios between non-synonymous (dN) and synonymous (dS) substitution among protein-coding sequences of all pairs of species were calculated in KaKs calculator  using a modified Yang-Nielsen algorithm. This ratio is currently employed to estimate if genes are under negative (purifying) selection (dN/dS < 1), positive (adaptative) selection (dN/dS > 1), or neutral evolution.
Phylogenetic analysis was carried out using Bayesian inference from the concatenated sequences of all protein coding genes of the five newly sequenced Culex specimens, plus the mitochondrial genomes of four species available in the GenBank, as follows: Cx. quinquefasciatus [HQ724617], Cx. pipiens [HQ724616, HQ724615], Aedes albopicutus [NC_006817] and Anopheles darlingi [NC_014275]. The latter two species were employed as outgroups. Multiple alignments were generated in Clustal X  and then edited visually in BioEdit  to ensure that genes were in the frame. The best-fit model of nucleotide substitution for each gene was estimated in JModeltest 0.1.1 package , under Akaike Information Criterion (AIC). Bayesian analysis was performed in MrBayes v. 3.1.2  and the dataset was partitioned by genes. Two independent runs with one cold and three heated chains each, were implemented for 1,000,000 generations with trees sampled every 100 generations; the consensus topology was generated after a burnin of 25 % of the retained trees.
Mitochondrial genome organization and nucleotide composition
Structural feature of the mitochondrial genome of four species of Cx. coronator, Cx. usquatus, Cx. camposi, and Cx. usquatissimus
Mitochondrial gene structure in base pairs
Cx. usquatissimus AC
Cx. usquatissimus RO
Overall nucleotide composition and A + T content per structural feature, and by codon position in protein-coding genes
%A + T
%A + T in PCG
Cx. usquatissimus AC
Cx. usquatissimus RO
The protein-coding genes (PCGs)
Except for the COX1 gene in which no canonical ATN was found, the majority of the PCGs of Culex species possessed the start codon ATN (Additional file 1). The absence of the ATN start codon in the COX1 gene was found in Cx. quinquefasciatus, Cx. pipiens, Ae. aegytpti and in other insects [31–35] and therefore it has been widely discussed . In 1983, de Brujin  proposed a tetranucleotide (ATAA) positioned upstream in the proximity of the opening codon of COX1 (TCG) as a start codon for Drosophila melanogaster. The author hypothesized that if an ATA, a recognized start codon, overlaps a stop codon TAA in the frame −1, the ATA of the tetranucleotide ATAA could function as a highly conserved 0 frame and as an opening frame. Notwithstanding, in the mitochondrial genome of Culex species, the closest ATAA tetranucleotide is located within the trnL2 gene, overlapping 36 base pairs of the gene. Consequently, it is plausible to suppose that the hexanucleotide ATTTAA, which usually flanks the beginning of the COX1 in mosquitoes, is involved in the initiation of the translation signal of this gene, as it was proposed for Tetrodontophora bielanensis .
Eleven genes (ATP6, ATP8, COX3, CytB, NAD2, NAD3, NAD4, NAD4L, NAD5, NAD6) have complete TAA stop codons, except for COX1 and COX2 genes, which terminate with an incomplete codon T. This phenomenon has been extensively reported in other mitochondrial genomes, especially in COX genes, and so far it has been related to posttranscriptional polyadenylation, during which residual adenines are added in order to originate TAA terminators [34, 36, 41].
Relative synonymous codon usage (RSCU) and overall average of codon usage in the PCGs of Cx. coronator, Cx. usquatus, Cx. camposi, and Cx. usquatissimus
The set of 22 tRNA genes in all species varied in length from 64 bp in trnR gene to 72 bp in trnV gene (Additional file 1). The AT content ranged from 78.6 % to 79.1 % (Table 3). The anticodons of all tRNAs in Culex species were the same as that of other insects [31–35] (Additional file 1). The conventional cloverleaf-like structure of tRNAs was found in 20 (Additional file 2) of the 22 tRNAs described for Culex species. In contrast, the trnR and trnS1 genes could not be folded into cloverleaf-like structure. In the trnS1, several similar cases were reported in previously sequenced metazoan mitochondrial genomes because of the replacement of the DHU (dihydrouridine) arm by a DHU-loop [43, 44]. The absence of a secondary structure in the trnS1 has been reported for Collembola species .
Similar to the mitochondrial genome of other metazoan species [32–37], the rrnL is located between trnL1 and trnV ribosomal genes, and the rrnS is situated between the trnV gene and the AT-rich region. In all species, the rrnL and rrnS are 2,124 base pairs long (Additional file 1). The AT content was 82 %, a value that is slightly higher than the overall content, however, similar to that observed in the 3rd codon position (Table 3). Both tRNAs and rRNAs genes are more conserved than the PCGs, with only 10 variable sites in tRNAs (throughout 1,462pb - 1,463 bp) and 11 substitutions among 2,124 bp of the two rRNAs genes. Evolutionary studies relative to the tRNA and rRNA in Drosophila and Anopheles showed signatures of negative selection, because non-pairing substitutions such as G-C to G-U can affect the secondary structure of these genes [41, 45].
Intergenic spacer sequences
A total of 52 bp to 56 bp non-coding nucleotides were found in the mitochondrial sequences of the species included in this study (Table 2). The length of each spacer mostly varies from 1 to 5 bp, except for two spacers: one located between the trnk and the tranD genes (12 bp in Cx. coronator, and 11 bp in the remaining species); and the other between trnS2 and NAD1 genes with 19 bp in all samples. Although those spacers account for only 0.3 % of the whole mitochondrial genomes, they are composed almost exclusively of adenine and thymine, with an AT content that is higher than that observed in the AT-control region, reaching 94.5 % in Cx. usquatissimus from Acre (Table 3).
AT-rich control region
The size variation, ranging from 719 bp in Cx. usquatus to 724 bp in Cx. coronator (Table 2) is within the range variation observed in the mitochondrial genomes reported for other Culicidae (from 531 bp in Anopheles gambiae to 731 bp in Culex quinquefasciatus). Those differences in length were caused by several indels of one to three nucleotides distributed throughout the AT regions of the mitochondrial genomes. Considering both indels and nucleotide substitutions, 34 variable sites were observed in Culex species. Similar to nuclear noncoding DNA, the AT-control region is subjected to a high level of polymorphisms that lead to variation in size and base composition. However, in species of the Anopheles gambiae complex, there is evidence that this region change more slowly than the 3rd codon position of PCGs . Additionally, among different subgenera of Anopheles, there are conservative poly-T regions of 15 to 18 bp . Comparing the control region of the four Culex species studied, we observed 4 conservative T-stretch regions of 6 to 15 bp. This was likely caused because this region contains several regulatory elements, including the origin of replication and transcription of the major strand of the mitochondrial genome  and, consequently, has a functional role, which is conserved among genera.
Sliding window analysis
The rrnL and rrnS demonstrated to be the most conserved genes within the mitochondrial genomes of species of the Coronator Group, showing low nucleotide diversity. Although rrnL usually shows high variability among some groups, i.e., Taeniidae , Nematoda , and Octopodus , herein, the highest nucleotide diversity found in the rrnL was 0.02. Per se, the rrnL and rrnS genes seem not to be good markers for taxonomy and evolutionary studies focusing on species of the Coronator Group.
Considering both the epidemiological importance of Cx. coronator and the problematic identification of females based on morphological characters, the development of molecular markers can facilitate species identification as well as establish phylogenetic relationships within the group and other closely related species. Taken together, results of the sliding window analysis and the dN/dS ratios indicate that the NADH5, ATP8 and ATP6 genes represent potential markers for future investigations regarding species of the Coronator Group, primarily because they retain high variations among PCGs with an elevated proportion of non-synonymous substitutions. Notwithstanding, the low levels of diversity shown herein, demonstrates that the utility of sequences of the mitochondrial gene needs further evaluation in future studies employing a larger sample size and other species of the subgenus Culex as well as other Culex (Culex) species.
Average of nucleotide divergence under Maximum Composite Likelihood of 13 protein-coding genes of the mitochondrial genome of Cx. coronator, Cx. usquatus, Cx. camposi, and Cx. usquatissimus and four specimens obtained from GenBank. The averages are expressed as percentages; GenBank accession numbers are indicated between brackets
1. Cx. pipiens Greece [HQ724615]
2. Cx. pipiens Turkey [HQ724616]
3. Cx. quinquefasciatus USA [HQ724617]
4. Cx. usquatissimus AC
5. Cx. usquatissimus RO
6. Cx. usquatus
7. Cx. camposi
8. Cx. coronator
The complete mitochondrial genomes of Cx. usquatus, Cx. usquatissimus, Cx. coronator, and Cx. camposi share the same gene order and gene content of other Culicidae species with the mitochondrial genome reported to date. The overall length and the AT content of the four Culex species are within the range observed in species of the Anopheles albitarsis complex, Cx. quinquefasciatus, and Ae. aegypti. Except for trnR and trnS1, the tRNAs were folded into the typical cloverleaf-like structure. Similar to other metazoan and species of the family Culicidae, the control region presented some conservative poly T- stretches, which seems to be involved in the functional roles of mitochondrial genome. Results of the Bayesian phylogeny corroborated the monophyly of the Coronator Group. Nevertheless, new insights on the evolution of the group can be addressed using a larger species sampling, including morphologically similar species and species of other Culex subgenera, for example, Lutzia, Phenacomyia and Phytothelmatomyia species. Probably because of the functional roles played in the respiratory chain by the mitochondrial genome, PCGs showed high purifying selection and low nucleotide divergence among the species investigated here.
This work was supported by FAPESP (Grant 2014/26229–7) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Grant 301666/2011–3).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Forattini OP. Entomologia Médica. 2° Volume. Culicini: Culex, Aedes e Psorophora. São Paulo, Editora da Universidade de São Paulo; 1965. Pages: 121–129.Google Scholar
- Bram R. Classification of Culex subgenus Culex (Diptera: Culicidae). Proc US Nat Mus. 1967;120:1–123.View ArticleGoogle Scholar
- Demari-Silva B, Vesgueiro FT, Sallum MAM, Marrelli MT. Taxonomic and phylogenetic relationships between species of the genus Culex (Diptera: Culicidae) from Brazil inferred from the cytochrome c oxidase I mitochondrial gene. J Med Entomol. 2011;48:272–9.View ArticlePubMedGoogle Scholar
- Laurito M, de Oliveira TMP, Almirón WR, Sallum MAM. COI barcode versus morphological identification of Culex (Culex) (Diptera: Culicidae) species: a case study using samples from Argentina and Brazil. Mem Inst Oswaldo Cruz. 2013;108:110–22.PubMed CentralView ArticlePubMedGoogle Scholar
- Vesgueiro FT, Demari-Silva B, Malfronte RS, Sallum MAM, Marrelli MT. Intragenomic variation in the second internal transcribed spacer of the ribosomal DNA of species of the genera Culex and Lutzia (Diptera: Culicidae). Mem Inst Oswaldo Cruz. 2011;106:1–8.View ArticlePubMedGoogle Scholar
- Consoli RAGB, Lourenço-de-Oliveira R. Principais mosquitos de importância sanitária no Brasil. Rio de Janeiro: FIOCRUZ; 1998.Google Scholar
- Forattini OP. Culicidologia Médica: identificação, biologia e epidemiologia. São Paulo: EDUSP; 2002. pp. 702–703.Google Scholar
- Alto BW, Connelly CR, O’Meara GF, Hickman D, Karr N. Reproductive biology and susceptibility of Florida Culex coronator to infection with West Nile Virus. Vector Borne Zoonotic Dis. 2014;14:606–14.PubMed CentralView ArticlePubMedGoogle Scholar
- Saraste M. Oxidativephosphorylationat thefindesiècle. Science. 1999;283:1488–93.View ArticlePubMedGoogle Scholar
- Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27:1767–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Goddard JM, Wolstenholme DR. Origin and direction of replication in mitochondrial DNA molecules from Drosophila melanogaster. Proc Natl Acad Sci U S A. 1978;75:3886–90.PubMed CentralView ArticlePubMedGoogle Scholar
- Botero-Castro F, Tilak MK, Justy F, Catzeflis F, Delsuc F, Douzery EJ. Next-generation sequencing and phylogenetic signal of complete mitochondrial genomes for resolving the evolutionary history of leaf-nosed bats (Phyllostomidae). Mol Phylogenet Evol. 2013;693:728–39.View ArticleGoogle Scholar
- Schierwater B, Stadler P, Desalle R, Podsiadlowski L. Mitogenomic and metazoan evolution. Mol Phylogenet Evol. 2013;69:311–2.View ArticlePubMedGoogle Scholar
- Williams ST, Foster PG, Littlewood DT. The complete mitochondrial genome of a turbinid vetigastropod from MiSeq Illumina sequencing of genomic DNA and steps towards a resolved gastropod phylogeny. Gene. 2014;533:38–47.View ArticlePubMedGoogle Scholar
- Foster PG, Bergo ES, Bourke BP, Oliveira TM, Nagaki SS, Sant’Ana DC, et al. Phylogenetic analysis and DNA-based species confirmation in Anopheles (Nyssorhynchus). PLoS One. 2013;8:e54063.PubMed CentralView ArticlePubMedGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.PubMedGoogle Scholar
- Kambhampati S, Smith PT. PCR primers for amplification of four insect mitochondrial gene fragments. Insect Mol Biol. 1995;4:233–6.View ArticlePubMedGoogle Scholar
- Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. Comput Sci Biol Proc. 1999;99:45–56. German Conference on Bioinformatics GCB.Google Scholar
- Milne I, Stephen G, Bayer M, Cock PJ, Pritchard L, Cardle L, et al. Using Tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 2013;14:193–202.View ArticlePubMedGoogle Scholar
- Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, et al. MITOS: Improved de novo Metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69:313–9.View ArticlePubMedGoogle Scholar
- Birney E, Durbin R. Using Gene Wise in the Drosophila annotation experiment. Genome Res. 2000;10:547–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Phylogenet Evol. 2013;30:2725–9.View ArticleGoogle Scholar
- Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995;41:353–8.View ArticlePubMedGoogle Scholar
- Schattner P, Brooks AN, Lowe TM. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005;1:33.Google Scholar
- Librado P, Rozas J. DnaSPv5: a software for comprehensive analysis of DNA. Bioinformatics. 2009;25:1451–2.View ArticlePubMedGoogle Scholar
- Zhang Z, Li J, Zhao XQ, Wang J, Wong GK, Yu J. KaKs Calculator: Calculating Ka and Ks through model selection and model averaging. Genomics Proteomics Bioinformatics. 2006;4:259–63.View ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. ClustalW and ClustalX version2.0. Bioinformatics. 2007;23:2947–8.View ArticlePubMedGoogle Scholar
- BioEdit HT. a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser (Oxf). 1999;41:95–8.Google Scholar
- Posada D. jModelTest:phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.View ArticlePubMedGoogle Scholar
- Beard CB, Hamm DM, Collins FH. The mitochondrial genome of the mosquito Anopheles gambiae: DNA sequence, genome organization, and comparisons with mitochondrial sequences of other insects. Insect Mol Biol. 1993;2:103–24.View ArticlePubMedGoogle Scholar
- Behura SK, Lobo NF. Haas B, deBruyn B, Lovin DD, Shumway MF et al. Complete sequences of mitochondrial genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes. Insect Biochem Mol Biol. 2011;41:770–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Logue K, Chan ER, Phipps T, Small ST, Reimer L, Henry-Halldin C, et al. Mitochondrial genome sequences reveal deep divergences among Anopheles punctulatus sibling species in Papua New Guinea. Malar J. 2013;12:64.PubMed CentralView ArticlePubMedGoogle Scholar
- Krzywinski J, Li C, Morris M, Conn JE, Lima JB, Povoa MM, et al. Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex. Mol Phylogenet Evol. 2011;58:469–77.PubMed CentralView ArticlePubMedGoogle Scholar
- Lloyd R, Streeter SD, Foster PG, Littlewood DTJ, Huntley J, Beckham GT, et al. The complete mitochondrial genome of Limnoria quadripunctata Holthuis (Isopoda: Limnoriidae). Mitochondrial DNA. 2014;10:1–2.Google Scholar
- Chen YH, Huang DY, Wang Y, Zhu CD, Hao JS. The complete mitochondrial genome of the endangered Apollo butterfly, Parnassius apollo (Lepidoptera: Papilionidae) and its comparison to other Papilionidae species. J Asia-Pacific Entomology. 2014;17:663–71.View ArticleGoogle Scholar
- Sun WY, Xu DL, Chen HX, Shi W, Sundberg P, Sun SC. Complete mitochondrial genome sequences of two parasitic/commensal nemerteans, Gononemertes parasite and Nemertopsis tetraclitophila (Nemertea: Hoplonemertea). Parasit Vectors. 2014;19:273.View ArticleGoogle Scholar
- Sun Z, Wan DG, Murphy RW, Ma L, Zhang XS, Huang DW. Comparison of base composition and codon usage in Insects. Genes Genomics. 2009;31:65–71.View ArticleGoogle Scholar
- Nardi F, Carapelli A, Fanciulli PP, Dallai R, Frati F. The complete mitochondrial DNA sequence of the basal hexapod Tetrodontophora bielanensis: evidence for heteroplasmy and tRNA translocations. Mol Biol Evol. 2001;18:1293–304.View ArticlePubMedGoogle Scholar
- de Bruijn MH. Drosophila melanogaster mitochondrial DNA, a novel organization and genetic code. Nature. 1983;304:234–41.View ArticlePubMedGoogle Scholar
- Krzywinski J, Grushko OG, Besansky NJ. Analysis of the complete mitochondrial DNA from Anopheles funestus: an improved dipteran mitochondrial genome annotation and a temporal dimension of mosquito evolution. Mol Phylogenet Evol. 2006;39:417–23.View ArticlePubMedGoogle Scholar
- Salvato P, Simonato M, Battisti A, Negrisolo E. The complete mitochondrial genome of the bag-shelter moth Ochrogaster lunifer (Lepidoptera, Notodontidae). BMC Genomics. 2008;9:331.PubMed CentralView ArticlePubMedGoogle Scholar
- Montooth KL, Abt DN, Hofmann JW, Rand DM. Comparative genomics of Drosophila mtDNA: novel features of conservation and change across functional domains and lineages. J Mol Evol. 2009;69:94–114.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao F, Huang DY, Sun XY, Shi QH, Hao JS, Yang Q. The first mitochondrial genome for the butterfly family Riodinidae (Abisara fylloides) and its systematic implications. Zoological Res. 2013;34:109–19.Google Scholar
- Cook CE. The complete mitochondrial genome of the stomatopod crustacean Squilla mantis. BMC Genomics. 2005;6:105.PubMed CentralView ArticlePubMedGoogle Scholar
- Caccone A, Garcia BA, Powell JR. Evolution of the mitochondrial DNA control region in the Anopheles gambiae complex. Insect Mol Biol. 1996;5:51–9.View ArticlePubMedGoogle Scholar
- Moritz C, Dowling TE, Brown WM. Evolution of Animal Mitochondrial DNA: Relevance for Population Biology and Systematics. Annu Rev Ecol Syst. 1987;18:269–92.View ArticleGoogle Scholar
- Hebert PD, Cywinska A, Ball SL, de Waard JR. Biological identifications through DNA barcodes. Proc Biol Sci. 2003;270:313–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Jia WZ, Yan HB, Guo AJ, Zhu XQ, Wang YC, Shi WG, et al. Complete mitochondrial genomes of Taenia multiceps, T.hydatigena and T.pisiformis: additional molecular markers for tape worm genus of human and animal health significance. BMC Genomics. 2010;11:447.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie Y, Zhang Z, Niu L, Wang Q, Wang C, Lan J, et al. The mitochondrial genome of Baylisascaris procyonis. PLoS One. 2011;6:e27066.PubMed CentralView ArticlePubMedGoogle Scholar
- Jabbar A, Beveridge I, Mohandas N, Chilton NB, Littlewood DT, Jex AR, et al. Analyses of mitochondrial amino acid sequence datasets support the proposal that specimens of Hypodontus macropi from three species of macropodid hosts represent distinct species. BMC Evol Biol. 2013;13:259.PubMed CentralView ArticlePubMedGoogle Scholar
- Broughton RE, Milam JE, Roe BA. The complete sequence of the zebra fish (Daniorerio) mitochondrial genome and evolutionary patterns invertebrate mitochondrial DNA. Genome Res. 2001;11:1958–67.PubMed CentralPubMedGoogle Scholar
- Le TH, Pearson MS, Blair D, Dai N, Zhang LH, McManus DP. Complete mitochondrial genomes confirm the distinctiveness of the horse-dog and sheep-dog strains of Echinococcus granulosus. Parasitology. 2002;124:97–112.View ArticlePubMedGoogle Scholar
- Cheng R, Zheng X, Ma Y, Li Q. The complete mitochondrial genomes of two octopods Cistopus chinensis and Cistopus taiwanicus: Revealing the phylogenetic position of the genus Cistopus within the Order Octopoda. PLoS One. 2013;17:e842.Google Scholar
- Lloyd RE, Foster PG, Guille M, Littlewood DT. Next generation sequencing and comparative analyses of Xenopus mitogenomes. BMC Genomics. 2012;13:496.PubMed CentralView ArticlePubMedGoogle Scholar
- Laurito M, Almirón W. Phylogenetic relationships of Culex (Culex) species (Diptera, Culicidae) from Argentina based on morphological characters. Zootaxa. 2013;3652:117–55.View ArticlePubMedGoogle Scholar