- Research article
- Open Access
First record of translocation in Culicidae (Diptera) mitogenomes: evidence from the tribe Sabethini
BMC Genomics volume 20, Article number: 721 (2019)
The tribe Sabethini (Diptera: Culicidae) contains important vectors of the yellow fever virus and presents remarkable morphological and ecological diversity unequalled in other mosquito groups. However, there is limited information about mitochondrial genomes (mitogenomes) from these species. As mitochondrial genetics has been fundamental for posing evolutionary hypotheses and identifying taxonomical markers, in this study we sequenced the first sabethine mitogenomes: Sabethes undosus, Trichoprosopon pallidiventer, Runchomyia reversa, Limatus flavisetosus, and Wyeomyia confusa. In addition, we performed phylogenetic analyses of Sabethini within Culicidae and compared its mitogenomic architecture to that of other insects.
Similar to other insects, the Sabethini mitogenome contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region. However, the gene order is not the same as that in other mosquitoes; the tyrosine (Y) and cysteine (C) tRNA genes have translocated. In general, mitogenome rearrangements within insects are uncommon events; the translocation reported here is unparalleled among Culicidae and can be considered an autapomorphy for the Neotropical sabethines.
Our study provides clear evidence of gene rearrangements in the mitogenomes of these Neotropical genera in the tribe Sabethini. Gene order can be informative at the taxonomic level of tribe. The translocations found, along with the mitogenomic sequence data and other recently published findings, reinforce the status of Sabethini as a well-supported monophyletic taxon. Furthermore, T. pallidiventer was recovered as sister to R. reversa, and both were placed as sisters of other Sabethini genera (Sabethes, Wyeomyia, and Limatus).
Insect mitochondrial DNA (mtDNA or the mitogenome) consists of a compact circular molecule usually 15–18 kb in length. It encodes 37 genes: 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. All of these genes are involved in translation of the PCGs . This complement of 37 genes is highly conserved in metazoan animals, with rare exceptions . The relative positions of PCGs and rRNA genes show limited variation among different phyla compared with tRNA genes . For this reason, it has been proposed that tRNA gene rearrangements may be useful phylogenetic markers [1, 4]. The ancestral insect mitogenome differs from the ancestral arthropod mitogenome only by the position of one tRNA gene . Although notable deviations from the ancestral insect mitogenomes in structure, gene arrangement, and gene content have been reported within insects, it is clear by analysing the numerous mitogenomes currently available that these exceptions are found only in highly derived portions of the insect phylogeny .
Among insects, one of the most studied groups is the family Culicidae (mosquitoes) because of their medical and veterinary importance worldwide . A total of 3564 species of Culicidae are currently recognised , classified in two subfamilies and 113 genera. There are approximately 1 hundred mitogenomes available representing several mosquito species, mainly from medically important genera of pathogen-vectoring species such as Aedes, Anopheles and Culex . All currently available Culicidae mitogenomes have similar gene order and arrangement to the ancestral insect mitogenome [8,9,10,11], with the exception of an inversion reported by Beard et al. . Despite their importance as transmitters of the yellow fever virus, to the best of our knowledge, there is only one study addressing the phylogenetic relationships within the tribe Sabethini, and the relationship between sabethines and other culicids using mitogenomes . Although Aragão et al.  did not address gene arrangements within the mitogenome, there are some unique characteristics of the mitochondria of sabethines that separate them from other culicids. Members of this tribe are among the most morphologically and biologically diverse of all mosquitoes, and present greatest species diversity in the tropics . In general, immature Sabethini species use as their habitat phytotelma, which can be formed by leaves, flowers, bamboo internodes, bromeliads, pitcher plants, and many other plant structures [14, 15].
Although long recognized as a monophyletic group, the taxonomic placement of sabethines within the Culicidae has varied, as well as the genus relationships within the tribe (see Judd ). Some advances have been made in the classification of this tribe [16,17,18,19,20], but the Sabethini as a whole is far from being organized into natural groups . Results of the analyses corroborate the classification proposed by Judd in recognizing both the tribe Sabethini and Neotropical sabethines as monophyletic groups . Currently, there are two major groups within the Sabethini, “New World” (Neotropical) and “Old World” (mainly Oriental/Australian) groups [14, 22] with approximately 700 valid species assigned into 13 genera . Regarding the Sabethini taxonomic placement within Culicidae, there is considerable controversy in the literature. Harbach and Kitching  placed tribe Sabethini close to tribe Culicini, based on morphological characters. In contrast, Reidenbach et al. , using morphological and genetic data, recovered the Sabethini and Aedini tribes as closely related.
In light of these knowledge gaps and this controversial situation, we started investigating the phylogeny of Culicidae at the genus level based on mitogenomes (Lorenz et al. MS submitted). Results corroborated the monophyly of Sabethini and incidentally unveiled tRNA translocations in some samples of that tribe. Then, we decided to perform a new study, now at the species level, to further understand the information obtained from sabethine mitochondrial data. The goals of the current study were to: 1) verify the arrangement and order of genes in the mitogenomes of Sabethini (Sabethes undosus, Trichoprosopon pallidiventer, Runchomyia reversa, Limatus flavisetosus, and Wyeomyia confusa) and 2) establish phylogenetic relationships among Sabethini species and with other culicids. Our results showed that the gene order contains phylogenetic information that complements the monophyly of the tribe Sabethini.
General characteristics of the mitogenome of the five newly sequenced Sabethini species from this study are in Table 1. The mitochondrial genomes are similar to those of three Sabethes species studied by Aragão et al. , ranging from 15,302 to 16,037 bp in length. Length variation is mainly due to differences in the length of control regions in Sabethini species. In addition to this AT-rich control region, the mitogenomes also contained 37 genes: 13 PCGs, 22 tRNA genes, and two rRNA genes (Fig. 1).
Gene order in the mitogenomes in the Sabethini differed from other mosquito species. Tyrosine (Y) and cysteine (C) tRNAs are translocated (Fig. 2). In the available mitogenomes from other groups, Y and C are located between COI and tRNA-W genes, whereas in the Sabethini, they are between the tRNA-I gene and the AT control region. Our results suggest that this translocation is an autapomorphy of the Neotropical sabethines.
Results of phylogenetic analyses showed that the representative species of Sabethini are monophyletic (Fig. 3). The phylogenetic analysis based on BI and ML recovered identical tree topologies, with variations in the statistical support of some branches. The three Sabethes species available in GenBank clustered with S. undosus. The strongly supported clade composed of T. pallidiventer and R. reversa was sister to the remaining sabethines. All clades corroborate the taxonomic classification of Sabethini proposed by Reidenbach et al.  based on morphology and sequence data.
Regarding the translocation found in Sabethini, we searched the literature to determine whether such translocation events are common in other insect orders and mapped them onto mtDNA following Cameron  (Fig. 4). We represent only a simplification of few insect orders to exemplify rearrangements occurrence; Embioptera, Mantodea, and Dermaptera also present this phenomenon. There are many groups in which a given rearrangement has been observed in only a single species, but the taxonomic range of these rearrangements is currently unknown. However, it is apparent that mitogenome rearrangements within insects, with some exceptions such as hemipteroids (Phthiraptera, Psocoptera and Thysanoptera) and hymenopterans, are uncommon events.
The tribe Sabethini has remarkable morphological and ecological diversity, unparalleled by other Culicidae tribes. Although widely recognized as a distinct group, the phylogenetic relationships within Sabethini are poorly known . Here, we report for the first time the mitogenomes of five Sabethini species. Knowledge of the mitogenomes of these species can be useful in elucidating phylogenetic relationships within Sabethini. Mitogenomes for the newly sequenced species are similar to those of other Sabethes species studied by Aragão et al. . The control region shows greater variability in length than other regions of the mitogenome, as has been observed in other Culicidae species . Similar to other Culicidae [9, 24], the third codon positions has higher AT content than the first and second codons. The transcription hypothesis of codon usage proposed by Sun et al.  says that “the high availability of ATPs, along with the lack of other NTPs, leads to the maximization of the use of adenines in the third codon position, increasing the efficiency of the transcription”. Furthermore, purifying selection against deleterious mutations is less severe on the third codon position . Thus, the higher AT content in the third codon position is likely associated with biased usage within synonymous codons.
Regarding gene order, the translocation reported here is likely an autapomorphy of the Neotropical sabethines. This evolutionary event occurred in the Sabethini ancestor at least 75 million years ago based on the dates inferred by Reidenbach et al. . Translocation events within tRNA gene regions have been found in Collembola , but very few gene rearrangements have been examined in invertebrates, thus it is difficult to hypothesize about the mechanisms involved in Sabethini mitogenome translocations. Tandem duplications have been proposed to produce translocations . Alternatively, translocations may be produced through the illicit priming of mitochondrial replication by a tRNA molecule (originally proposed by Cantatore et al. ). Specifically, after mitochondrial replication is initiated, failure to cleave the tRNA primer from the nascent DNA strand could lead to the ultimate incorporation of a tRNA gene into the mitochondrial genome. It is unclear which model is more likely to have occurred in Sabethini.
Sabethini phylogenetic relationships were robustly resolved by analysis of mitogenomes (Fig. 3). The lineage composed of Trichoprosopon pallidiventer and Runchomyia reversa was sister to the group of Sabethes, Wyeomyia, and Limatus, which formed a highly supported clade. The phylogenetic tree recovered here supports those found by Reidenbach et al.  using six nuclear genes and 80 morphological characters. Our results are also in agreement with those of Aragão et al. , where the species of Sabethes were monophyletic.
Other contemporary studies of our team lead one to believe in the monophyly of Sabethini based on mitogenomes (Lorenz et al., MS submitted) and pointed to putative autapomorphic wing shape characters (Lorenz & Suesdek, MS submitted). In addition, our study also addressed the phylogenetic placement of the genus Runchomyia using DNA sequence data for the first time. This genus is the most difficult to characterize of the sabethines studied here because no clear autapomorphies are known for this genus from any life stage . Along with Runchomyia, the genus Trichoprosopon represent the most divergent branch within Sabethini, as proposed by Edwards . This group has both morphological synapomorphies with Neotropical taxa and retains many of plesiomorphic characteristics of Oriental and Australian groups . Regarding Sabethini phylogenetic placement within Culicidae, our findings suggest that this tribe is closer to both Culicini and Aedini than to the remaining groups studied. This finding differ from those obtained by Harbach and Kitching  and Reidenbach et al.  using distinct data sets.
All species studied here are Neotropical, and previous studies strongly support these taxa as a monophyletic group [14, 30]. Species of the genera Trichoprosopon and Sabethes have been found to be competent to transmit arboviruses to humans and other primates . The genus Wyeomyia has also been implicated in the transmission of yellow fever virus . The geographical range of Sabethes, Trichoprosopon and Wyeomyia is restricted to the Neotropics. Furthermore, Sabethini from the Oriental/Australian region are not vectors of arboviruses. This suggests that vector competence in Neotropical Sabethini for arboviruses may have arisen after Gondwanan separation of the Americas and Africa . The present study also indicates that vector competence to arboviruses may have arisen more than once in Culicidae, in agreement with the hypothesis proposed by Judd . An alternative and parsimonious explanation is that all culicids share vectorial competence but it is expressed unequally across species.
Currently, with representative mitogenomes sequenced from each insect order, no gene rearrangements have been found to be shared between orders . However, we found that synapomorphic mtDNA rearrangements can distinguish groups at taxonomic scales below the ordinal level. Rearrangements occur throughout the mitogenome, with inversions or local rearrangements being most common, whereas translocations are rare. Most gene translocations found in insect mitogenomes are useful for high-level phylogenetic reconstructions [5, 33]. Nevertheless, the results found in the present study suggest that gene order can be informative at lower taxonomic ranks than previously hypothesized, and reinforce the findings from sequence based phylogenetic analysis. The remaining genes in sabethine mitogenomes have the same order as the ancestral hexapodan gene arrangement . Further studies addressing the mitogenomes of Sabethini will provide new insights and additional tools to investigate relationships within a group of mosquitoes of public health importance. In addition, the translocation of two tRNA genes is a remarkable phenomenon that merits investigation with broader taxon sampling.
In this study, we report the complete mitochondrial genomes of five Sabethini species and the occurrence of a tRNA translocation. The length of these mitogenomes was similar to that of other Sabethes species, and the AT-rich control region was the most length variable region. PCG third codon positions had a higher AT content than the first or the second codons, similar to other Culicidae species or animals in general. Regarding gene order, the translocation reported here is unprecedented among Culicidae and represents an autapomorphy for the Neotropical sabethines. Our analyses suggest that gene order can be informative at low taxonomic levels, such as tribes, reinforcing the results from sequence based phylogenetic analysis. The mitogenomes also corroborated the monophyly of Sabethini and provided robust phylogenetic support for the lineage of Trichoprosopon plus Runchomyia as sisters to a clade composed of Sabethes, Wyeomyia, and Limatus. We suggest that new insights on the phylogeny of the group could be addressed further by using larger taxon sampling, including Oriental and Australian taxa.
Sample collection, data used, and DNA extraction
Mosquito species used in this study share similar habitats, and were captured in the same region of Atlantic forest, São Paulo state, Brazil (24°53′S, 47°51′W), except for T. pallidiventer, which was captured in the Amazon forest, Amazonas, Brazil (3°05′S, 60°00′W). Field-collected specimens were identified to species level using morphological keys of Forattini  and Lane . Adult females were collected and stored at − 80 °C until DNA extraction. Genomic DNA was extracted from each mosquito individually using a Qiagen DNeasy® Blood and Tissue Kit (Qiagen Ltd., Crawley, UK) following the same procedure as Foster et al. . DNA was examined after electrophoresis on 1% agarose gels and used for PCR amplification. A list of all species used in this study is in Table 2.
PCR amplification and sequencing
The mitochondrial genome of each individual was amplified using the one-step long-range polymerase chain reaction method . Reactions were prepared as follows: 25 μL of GoTaq Long PCR Master Mix 2x (Promega, WI), 0.3 μM of each primer (forward and reverse), 1–2 μL of extracted DNA, and sterile water to make up final volume of 50 μL. The primers used for mitochondrial genome amplification follow Hwang et al. : HPK16Saa for the forward strand and HPK16Sbb for the reverse strand. The size of the amplified fragment was approximately 15,300 bp. Thermocycler conditions consisted of: initial denaturation at 94 °C for 2 min, followed by 39 cycles of 94 °C for 15 s, annealing at 65 °C for 20 s, and 65 °C for 15 min, and one final extension cycle at 72 °C for 10 min. The only step that varied between species was annealing temperature: 65 °C for Wyeomyia confusa and Trichoprosopon pallidiventer, 62 °C for Limatus flavisetosus, and 55 °C for Sabethes undosus and Runchomyia reversa. Amplicons were purified using a DNA Clean & Concentrator kit (Zymo Research, CA) and quantified using a Qubit 2.0 fluorometer (Life Technologies, OR).
Mitogenome sequence assembly and analysis
Next-generation sequencing was employed to obtain mitochondrial DNA sequences of five Sabethini species. Barcoded libraries were constructed from the PCR products using the Nextera XT DNA Sample Preparation Kit (Illumina, IL) and sequenced on the Illumina MiSeq platform with paired-end 250 bp read chemistry. Raw nucleotide sequences were checked using FastQC software . Mitochondrial genomes were assembled using Mira v4  and Newbler v2.9 software, and visualized in Tablet . Mitochondrial genes were annotated using MITOS . DNA and amino acid (AA) sequence alignments were performed with MEGA 6.06 software  and the ClustalW algorithm. Phylogenetic inferences were made using the 13 mitochondrial PCGs from the five Sabethini species sequenced here and other Culicidae available on GenBank (see Table 2). To perform a more comprehensive phylogenetic analysis of the tribe Sabethini, in addition to the newly sequenced samples from our study, we also included in our analysis three species of Sabethes previously sequenced by Aragão et al. : Sabethes chloropterus (MF957172), Sabethes glaucodaemon (MF957173), and Sabethes belisarioi (MF957171). We tested all of the following dataset combinations: amino acids, all DNA sites, or 1st and 2nd codon positions alone. All tests produced the same topology, differing only nodal support. Sequences were concatenated using FASconCAT-G v. 1.04 . All phylogenetic analyses were partitioned based on gene, with gamma-distributed site heterogeneity, and substitution models estimated independently for each partition. Phylogenetic analysis was performed using Maximum Likelihood (ML) with bootstraping = 1000 in RAxML v8.2.11 software . Bayesian phylogenetic inference (BI) was conducted using MrBayes v. 3.2.6  with four runs of eight chains each, run for at least 5 million generations. The tree generated was visualized and customized in FigTree v1.4.3 software .
Availability of data and materials
All the nucleotide data generated in this study were submitted to the NCBI sequence nucleotide database under the accession numbers: MK575482, MK575487, MK575488, MK575490, and MK575492.
Cameron SL. Insect mitochondrial genomics: implications for evolution and phylogeny. Annu Rev Entomol. 2014;59:95–117.
Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27:1767–80.
Pääbo S, Thomas WK, Whitfield KM, Kumazawa Y, Wilson AC. Rearrangements of mitochondrial transfer RNA genes in marsupials. J Mol Evol. 1991;33:426–30.
Macey JR, Larson A, Ananjeva NB, Fang Z, Papenfuss TJ. Two novel gene orders and the role of light-strand replication in rearrangement of the vertebrate mitochondrial genome. Mol Biol Evol. 1997;14:91–104.
Boore JL, Lavrov DV, Brown WM. Gene translocation links insects and crustaceans. Nature. 1998;392:660–7.
Foster WA, Walker ED. Mosquitoes (Culicidae). In Medical and veterinary Entomology: Academic Press; 2019.
Mosquito Taxonomic Inventory. 2019. http://mosquito-taxonomic-inventory.info/. Accessed 20 May 2019.
GenBank in National Center for Biotechnology Information. 2018. https://www.ncbi.nlm.nih.gov/. Accessed 13 Dec 2018.
Demari-Silva B, Foster PG, de Oliveira TM, Bergo ES, Sanabani SS, Pessôa R, Sallum MAM. Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera: Culicidae), members of the coronator group. BMC Genomics. 2015;16:831.
Oliveira TMP, Foster PG, Bergo ES, Nagaki SS, Sanabani SS, Marinotti O, Sallum MAM. Mitochondrial genomes of Anopheles (Kerteszia)(Diptera: Culicidae) from the Atlantic forest, Brazil. J Med Entomol. 2016;53:790–7.
Hao YJ, Zou YL, Ding YR, Xu WY, Yan ZT, Li XD, Chen B. Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes. Sci Rep. 2017;7(1):7666.
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.
Aragão A, Neto JPN, Cruz ACR, Casseb SMM, Cardoso JF, da Silva SP, Ishikawa EAY. Description and phylogeny of the mitochondrial genome of Sabethes chloropterus, Sabethes glaucodaemon and Sabethes belisarioi (Diptera: Culicidae). Genomics. 2018;6:1–8.
Judd D. Review of the systematics and phylogenetic relationships of the Sabethini (Diptera: Culicidae). Syst Entomol. 1996;21:129–50.
Frank JH, Lounibos LP. Phytotelmata: swamps or islands? Fla Entomol. 1987;70:14–20.
Zavortink TJ. Mosquito studies (Diptera, Culicidae). XXXV. The new sabethine genus Johnbelkinia and a preliminary reclassification of the composite genus Trichoprosopon. Contr Am Ent Inst. 1979;12:1-61.
Mattingly PF. Medical entomology studies the subgenera Rachiontomyia Tricholeptomyia and Tripteroide in the oriental region (Diptera: Culicidae). Contrib Am Ent Inst. 1981;17:1–147.
Harbach RE, Peyton EL. A new subgenus of Wyeomyia (Diptera: Culicidae), with the reclassification and redescription of Wyeomyia (Davismyia) arborea, Wyeomyia (Dendromyia) tarsata and Sabethes (Sabethes) carrilloi. Mosq Syst. 1991;23:92-109.
Harbach RE, Peyton EL. Morphology and evolution of the larval maxilla and its importance in the classification of the Sabethini (Diptera: Culicidae). Mosq Syst. 1993;25:1-16.
Motta MA, Lourenço-de-Oliveira R. Wyeomyia luteoventralis Theobald, the type species of the subgenus Dendromyia Theobald (Diptera: Culicidae). Mem Inst Osw Cruz. 1995;90:375–85.
Harbach RE, Kitching IJ. Phylogeny and classification of the Culicidae (Diptera). System Entomol. 1998;23:327–70.
Ward RA. Third supplement to ‘A catalog of the mosquitoes of the World' (Diptera: Culicidae). Mosq Syst. 1992;24:177-230.
Reidenbach KR, Cook S, Bertone MA, Harbach RE, Wiegmann BM, Besansky NJ. Phylogenetic analysis and temporal diversification of mosquitoes (Diptera: Culicidae) based on nuclear genes and morphology. BMC Evol Biol. 2009;9:298–302.
Krzywinski J, Li C, Morris M, Conn JE, Lima JB, Povoa MM, Wilkerson RC. Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex. Mol Phylogenet Evol. 2011;58:469–77.
Sun WY, Xu DL, Chen HX, Shi W, Sundberg P, Strand M, Sun SC. Complete mitochondrial genome sequences of two parasitic/commensal nemerteans, Gononemertes parasita and Nemertopsis tetraclitophila (Nemertea: Hoplonemertea). Parasit Vectors. 2014;7:273.
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.
Dowton M, Austin AD. Evolutionary dynamics of a mitochondrial rearrangement" hot spot" in the hymenoptera. Mol Biol Evol. 1999;16:298–309.
Cantatore P, Gadaleta MN, Roberti M, Saccone C, Wilson AC. Duplication and remoulding of tRNA genes during the evolutionary rearrangement of mitochondrial genomes. Nature. 1987;329:853.
Edwards FW. Diptera: Fam. Culicidae. Gen. Ins. 1932;194:30-67.
Harbach RE. The Culicidae (Diptera): a review of taxonomy, classification and phylogeny. Zootaxa. 1668;2007:591–638.
Mattingly PF. Tr.(Tripteroides). Mosq Syst. 1974;6:231.
Galindo P, Trapido H. Forest mosquitoes associated with sylvan yellow fever in Nicaragua. Am J Trop Med Hyg. 1957;6:145–52.
Boore JL, Brown WM. Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998;8:668–74.
Crease TJ. The complete sequence of the mitochondrial genome of Daphnia pulex (Cladocera: Crustacea). Gene. 1999;233:89–99.
Forattini OP. Culicidologia médica: identificaçäo, biologia e epidemiologia. In: Culicidologia Médica: Identificaçäo, Biologia E Epidemiologia, vol. 2. In press; 2002.
Lane J. Neotropical Culicidae, vol. 1. Sao Paulo, In press: University of Sao Paulo; 1953.
Foster PG, Bergo ES, Bourke BP, Oliveira TM, Nagaki SS, Sant’Ana DC, Sallum MAM. Phylogenetic analysis and DNA-based species confirmation in Anopheles (Nyssorhynchus). PLoS One. 2013;8:e54063.
Matthews BJ, Dudchenko O, Kingan SB, Koren S, Antoshechkin I, Crawford JE, Weedall GD. Improved reference genome of Aedes aegypti informs arbovirus vector control. Nature. 2018;563:501–9.
Zhang H, Xing D, Wang G, Li C, Zhao T. Sequencing and analysis of the complete mitochondrial genome of Aedes albopictus (Diptera: Culicidae) in China. Mitochondr Dna Part A. 2016;27:2787–8.
Luo QC, Hao YJ, Meng F, Li TJ, Ding YR, Hua YQ, Chen B. The mitochondrial genomes of Culex tritaeniorhynchus and Culex pipiens pallens (Diptera: Culicidae) and comparison analysis with two other Culex species. Parasit Vectors. 2016;9:406.
Hardy CM, Court LN, Morgan MJ. The complete mitochondrial DNA genome of Aedes vigilax (Diptera: Culicidae). Mitochondr Dna Part A. 2016;27:2552–3.
Lemos PDS, Monteiro HADO, Castro FC, Lima CPSD, Silva DEAD, Vasconcelos JMD, Nunes MRT. Characterization of mitochondrial genome of Haemagogus janthinomys (Diptera: Culicidae). Mitochondr Dna Part A. 2017;28:50–1.
Moreno M, Marinotti O, Krzywinski J, Tadei WP, James AA, Achee NL, Conn JE. Complete mtDNA genomes of Anopheles darlingi and an approach to anopheline divergence time. Malaria J. 2010;9:127.
Foster PG, de Oliveira TMP, Bergo ES, Conn JE, Sant’Ana DC, Nagaki SS, Sallum MAM. Phylogeny of Anophelinae using mitochondrial protein coding genes. Roy Soc Open Sci. 2017;4:170758.
Hwang UW, Park CJ, Yong TS, Kim W. One-step PCR amplification of complete arthropod mitochondrial genomes. Mol Phylogenet Evol. 2001;19:345–52.
Andrews S. FastQC: quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc. Accessed 20 Oct 2018.
Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. Comput Sci Biol Proc. 1999;99:45–56.
Milne I, Bayer M, Cardle L, Shaw P, Stephen G, Wright F, Marshall D. Tablet - next generation sequence assembly visualization. Bioinformatics. 2009;26:401–2.
Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, Stadler PF. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69:313–9.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Kück P, Meusemann K. FASconCAT: convenient handling of data matrices. Mol Phylogenet Evol. 2010;56:1115–8.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Höhna S, Huelsenbeck JP. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Rambaut A. FigTree version 1.4.0. 2016. http://tree.bio.ed.ac.uk/software/figtree. Accessed 10 Oct 2018.
LS has been fellow of Conselho Nacional de Desenvolvimento Científico e Tecnológico (CAPES) grants 311805/2014–0 and 311984/2018–5. CL was granted by fellowships no. 2013/05521–9 and 2015/12784–1 from São Paulo Research Foundation (FAPESP). MAMS is financially supported by São Paulo Research Foundation (FAPESP) grant no. 2014/26229–9 and CNPq no. 301877/2016–5. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.