Mitochondrial pseudogenes in the nuclear genome of Aedes aegypti mosquitoes: implications for past and future population genetic studies
© Hlaing et al; licensee BioMed Central Ltd. 2009
Received: 06 October 2008
Accepted: 06 March 2009
Published: 06 March 2009
Mitochondrial DNA (mtDNA) is widely used in population genetic and phylogenetic studies in animals. However, such studies can generate misleading results if the species concerned contain nuclear copies of mtDNA (Numts) as these may amplify in addition to, or even instead of, the authentic target mtDNA. The aim of this study was to determine if Numts are present in Aedes aegypti mosquitoes, to characterise any Numts detected, and to assess the utility of using mtDNA for population genetics studies in this species.
BLAST searches revealed large numbers of Numts in the Ae. aegypti nuclear genome on 146 supercontigs. Although the majority are short (80% < 300 bp), some Numts are almost full length mtDNA copies. These long Numts are not due to misassembly of the nuclear genome sequence as the Numt-nuclear genome junctions could be recovered by amplification and sequencing. Numt evolution appears to be a complex process in Ae. aegypti with ongoing genomic integration, fragmentation and mutation and the secondary movement of Numts within the nuclear genome.
The PCR amplification of the putative mtDNA nicotinamide adenine dinucleotide dehydrogenase subunit 4 (ND4) gene from 166 Southeast Asian Ae. aegypti mosquitoes generated a network with two highly divergent lineages (clade 1 and clade 2). Approximately 15% of the ND4 sequences were a composite of those from each clade indicating Numt amplification in addition to, or instead of, mtDNA. Clade 1 was shown to be composed at least partially of Numts by the removal of clade 1-specific bases from composite sequences following enrichment of the mtDNA. It is possible that all the clade 1 sequences in the network were Numts since the clade 2 sequences correspond to the known mitochondrial genome sequence and since all the individuals that produced clade 1 sequences were also found to contain clade 2 mtDNA-like sequences using clade 2-specific primers. However, either or both sets of clade sequences could have Numts since the BLAST searches revealed two long Numts that match clade 2 and one long Numt that matches clade 1. The substantial numbers of mutations in cloned ND4 PCR products also suggest there are both recently-derived clade 1 and clade 2 Numt sequences.
We conclude that Numts are prevalent in Ae. aegypti and that it is difficult to distinguish mtDNA sequences due to the presence of recently formed Numts. Given this, future population genetic or phylogenetic studies in Ae. aegypti should use nuclear, rather than mtDNA, markers.
Mitochondrial DNA (mtDNA) has been used extensively over the last three decades in population genetic and phylogenetic studies in a wide range of animals from Drosophila to humans (e.g. [1–4]). The advantages of mtDNA include its general lack of recombination that results in a single demographic history for the whole molecule, and high copy number which allows ease of amplification . Further, the relatively high mutation rate of mtDNA generates correspondingly high levels of polymorphism and divergence . This makes mtDNA particularly informative for the determination of genetic population structure and inference of population history within species as well as for deducing phylogenetic relationships between closely related species. MtDNA has been applied to many insect taxa including bees  and ants , in addition to medically important insects such as Anopheles [8–10] and Aedes mosquitoes [11, 12] where it is particularly important to estimate gene flow for vector control purposes [13, 14]. Most recently, the use of mtDNA sequences has been proposed for several DNA bar-coding initiatives for taxonomic identification and biodiversity assessment .
However, mtDNA also has its disadvantages. MtDNA data can be misleading in population and phylogenetic studies as it is particularly prone to selective sweeps ; it can introgress with relative ease between species ; and its population dynamics may be driven by intracellular symbionts . One problem that can be particularly difficult to detect is the presence of nuclear mitochondrial pseudogenes (Numts) [19, 20]. Numts result from the translocation of mitochondrial sequences from the mitochondrial genome into the nuclear genome and, once integrated, these non-functional sequences accumulate mutations freely. The potential for Numt amplification in addition to, or even instead of, the authentic target mtDNA sequence can seriously confound population genetic and phylogenetic analyses (reviewed in [21, 22]). For example, the mistaken inclusion of Numt sequences in an mtDNA phylogeny resulted in incorrect phylogenetic relationships being proposed within the South American bird genus Scytalopus [23, 24]. The high prevalence of Numts in gorillas has also been problematic as they initially obscured the presence of two genetically divergent groups of gorillas which has important implications for understanding their evolutionary history and future conservation [25–27].
In 2001, over 82 different eukaryotes, including 20 insect species, were known to have Numts . However, there have been numerous reports of Numts since then, for example, in gorillas , domestic cats , chickens , and several insects such as bees  and ants . Although taxonomically widespread, there is substantial variation among species in Numt copy number . Even within insects, Numt copy number varies greatly with high numbers in Tribolium flour beetles, honeybee and the brown mountain grasshopper [30, 33], but few or none in Drosophila and Anopheles mosquitoes . Although a positive correlation between Numt copy number and genome size is not clear cut [29, 32], genome size may partially explain taxonomic variation in Numt copy number . In this context it worth noting that the genome of Ae. aegypti is particularly large, 1.3 Gb . Numt fragments tend to be small, with the vast majority (>90%) less than 1 kb, and often less than 100 bp [30, 32, 35]. The longest Numts have been found in species with high Numt copy number, such as humans, mouse and rice . Controlling for size of mtDNA, the largest reported Numt to date is 14,654 bp in humans which corresponds to ~88% of the mitochondrial genome .
MtDNA has been widely used to study genetic population structure in Ae. aegypti mosquitoes, the organism of study here [37–41]. This mosquito is thought to have spread relatively recently from its ancestral home in Africa to the rest of the Tropics where they are important vectors of dengue and yellow fever [42, 43]. An understanding of genetic population structure and gene flow in Ae. aegypti is therefore of particular importance as this information is required for vector control [44–46]. The presence of Numts in Ae. aegypti could therefore seriously confound the interpretation of previous population genetic studies in Ae. aegypti as well as restrict the use of mtDNA in future studies.
Numt presence can manifest itself in several ways including: PCR ghost bands; extra bands in restriction profiles; sequence ambiguities; and unexpected phylogenetic placement ( and references therein). In our mitochondrial sequences of this species we found a significant proportion of sequence ambiguities suggestive of the presence of Numts. In general, previous authors have not reported sequence ambiguities in Ae. aegypti. The study by Paduan and Ribolla  is one exception but here the sequence ambiguities were attributed to heteroplasmy (where wild type and mutant mitochondrial genomes co-exist in a cell ). Another possible indication of the presence of Numts is the deep clade structure within many worldwide populations of Ae. aegypti. Two divergent mtDNA clades have been reported from: Mexico [37, 49]; Thailand ; Venezuela ; Australia ; and the Americas . The aims of this study are therefore to determine if Numts are present in Ae. aegypti mosquitoes, to characterise any Numts present, and to assess the utility of using mtDNA for population genetics studies in this species. Using a combined experimental and bioinformatics approach utilizing the recent genome sequence of Ae. aegypti (available at http://aaegypti.vectorbase.org/Genome/Home/)  we show here that Numts are prevalent in Ae. aegypti and that many previous studies have likely included Numt sequences in their mtDNA sequences. Finally, we discuss the implications of the use of mtDNA in future population genetic studies of this species.
Numt identification using BLAST searches
Size distribution of Ae. aegypti Numts (protein coding, tRNA and rRNA) detected by BLAST searches
Number of Numts
Although the majority of Numts are small in length, the penultimate column in Figure 1 shows that 16 supercontigs contain Numts that span two or more genes. Of these, 7 supercontigs have sequences corresponding to 6 or more contiguous mtDNA genes. The remaining supercontigs contain one to three smaller fragments of mtDNA-like sequence with each fragment spanning no more than three genes. For these longer Numts there appears to be some pattern across supercontigs in the genetic make up of the Numts. For example, supercontigs 1.348, 1.62 and 1.287 all contain a long Numt that spans the contiguous genes ND1, ND2, and CO1 as well as part of the non-contiguous ND4L gene. In addition, supercontigs 1.600, 1.765 and 1.929 span the contiguous genes ND2, CO1, CO2, ATP8, ATP6, CO3, ND3 and ND5. With the exception of the first two, this set of contiguous genes is also contained by the Numt on supercontig 1.357. Although these contiguous genes were detected by BLAST using each protein coding gene individually, alignments between the mtDNA and the supercontigs showed that in all cases the Numts were contiguous across all coding and non-coding regions with E values < 10-14. Despite their apparent similarity, the start and end positions of these long Numts differ in every case. There are also two supercontigs (1.363 and 1.593) that contain almost full length mtDNA copies.
Supercontigs 1.600 and 1.765 have even greater numbers of substitutions compared to the mtDNA sequence and this coupled with the presence of stop codons is very clear evidence that these are Numt sequences. In all five of these long Numts the majority of substitutions are synonymous. This is not inconsistent with their being Numts as many of the observed substitutions may have occurred in the mtDNA rather than the Numt [22, 24] as, based on Drosophila, the mitochondrial mutation rate is likely to be ~10 times higher than the nuclear mutation rate . It is interesting to note that the gene fragments at the start of the Numt on supercontig 1.363, and at the starts and ends of the Numts on supercontigs 1.600 and 1.765, have much higher numbers of substitutions (indicated in Figure 2) with respect to the mtDNA sequence than the remaining genes in these Numts. This pattern is unlikely to indicate misassembly since it is found on Numts with both high and low levels of differentiation from the mitochondrial sequence. This pattern more likely indicates a heterogenous ancestry for these Numts.
Recovered nuclear genomic sequences
Primers used for mtDNA specific amplification and for amplifying across Numt junctions
Fragment size (bp)
Clade 2 specific ND4 gene region
Start junction of Numt on 1.363
End junction of Numt on 1.363
Start junction of Numt on 1.593
Start junction of Numt on 1.1330
Deep clade structure of putative mitochondrial haplotypes
In addition to the 141 clear sequences we also obtained 25 ambiguous sequences, corresponding to ~15% of the total sample. The ambiguous sequences possessed double peaks at 9–17 sites all of which were sites that differentiated the two clades. Since the forward primer used here differed from those used in previous studies  we repeated the amplification and sequencing of 12 individuals with ambiguous sequences using the orginal primers that amplified a shorter 359 bp fragment corresponding to the 3' end of our amplicon. Even for these shorter fragments, the ambiguities were consistently present, in the same individuals and involving the same sites that differed between the two clades. This situation was not changed even by redesign of the reverse primers. These sequence ambiguities indicate that Numts are being amplified instead of, or together with, Ae. aegypti mtDNA. The core sequence of clade 2 has a 100% sequence identity with the mtDNA sequence of Ae. aegypti (GenBank accession number: EU352212) which indicates that clade 1 most likely comprises Numts.
Confirmation that clade 1 contains Numt sequences
If all the clade 1 sequences obtained from the Southeast Asian mosquitoes were Numts we would also expect an mtDNA sequence to be present in each individual. We therefore developed primers that amplify 282 bp of the ND4 gene from clade 2 (i.e. known mtDNA-like) sequences only (Table 2). These clade 2/mtDNA-specific primers had 3' termini at bases that distinguished clade 2 from clade 1. The PCR conditions under which these primers were specific for the clade 2/mtDNA sequence were determined using clade 1 and clade 2 PCR products as templates. Under these conditions the mtDNA-specific primers generated PCR products from all individuals that had previously given rise to unambiguous clade 1 sequences using the original primers. Nine of these amplicons were sequenced and all sequences were confirmed to fall into clade 2. For sequences of both clades to be present within an individual, at least one of the sequences must be a Numt. Based on the results of the mtDNA enrichment which showed clade 1 Numts and that the known mtDNA sequence corresponds to clade 2, the most parsimonious explanation is that all the clade 1 sequences are Numts. However, we cannot exclude the possibility that the clade 2/mtDNA specific primers are amplifying clade 2-like Numts such as those detected on supercontigs 1.593 and 1.363 by the BLAST searches.
Sequences from cloned PCR products
High numbers of Numts in the Ae. aegypti nuclear genome
The results of the BLAST searches demonstrate the presence of high numbers of Numts in the genome of Ae. aegypti with Numts detected on a total of 146 different supercontigs. This is comparable to other insects that are considered to have large numbers of Numts including the honeybee, Apis mellifera, and the Tribolium flour beetle which have Numts at an estimated 575 and 57 genomic locations, respectively . It is notable that other dipterans studied to date, D. melanogaster and An. gambiae, which have few or no Numts, have much smaller genomes than Ae. aegypti; 176 Mb and 278 Mb compared to 1.3 Gb for Ae. aegypti. Almost half the genome of Ae. aegypti is composed of transposable elements  which may indicate that this species is generally unable to rid itself of non functional DNA, including Numts.
In agreement with studies in other taxa [e.g. ], the majority of the Numts detected here were short. However, we also noted many long Numts with six being greater than 7 kb in length with the two longest on supercontigs 1.363 and 1.593 (Figure 1). Although these latter two differed very little from the mitochondrial sequence they appear to be real rather than the result of misassembly as there are some frameshift and stop mutations (albeit very few in number) and the presence of Numt-genome junction points was confirmed by amplification and sequencing. At 15.46 kb and 13.07 kb, the Numts on supercontigs 1.363 and 1.593, respectively, are amongst the longest reported to date. Other species with long Numts are humans , Arabidopsis, rice, house mice , domestic cats , and voles  that have Numts of 14.65 kb, 20.13 kb, 13.32 kb, 12.4 kb, 7.9 kb and 4 kb in length, respectively.
In general, as the Numts in Ae. aegypti get progressively shorter in length their similarity to the mitochondrial sequence decreases (Figure 2). This is consistent with the repeated movement of long mtDNA sequences to the nucleus which then decay through time becoming both fragmented and mutated. This is also consistent with our observation that, despite some variation among genes, all the mitochondrial genes are represented in Numts in our study. This is similar to the situation in humans [36, 54] and mouse and Ciona  but contrasts with several other studies which have noted a very marked preference for particular mtDNA genes to become Numts [29, 30]. The situation in which the Numt on supercontig 1.1330 was found to have a homologous Numt with the same end point but slightly different genomic flanking sequence has been found previously . Such Numt copies could have arisen from the secondary integration of Numts from one nuclear genomic location to another [28, 36, 55]. Genomic rearrangements involving Numts seem common . In accord with this, the pattern seen here of high substitution rates at the extreme ends of three of the long Numts indicates that these Numts could have been formed by the homologous recombination of long mtDNA sequences into genome positions that already contained Numts. However, it is also difficult to exclude the possibility of genome misassembly in giving rise to this pattern. Overall, our data lead to the conclusion that Numt evolution in Ae. aegypti is a very complex and dynamic process.
Origin of Clade 1 and Clade 2 type Numts
The genome sequence was obtained from a highly inbred strain of Ae. aegypti . It seems unfeasible that this strain could have maintained two mtDNA sequences, corresponding to clades 1 and 2, through successive bottlenecked generations. The evidence presented here demonstrates clearly that at least some, and possibly all, clade 1 sequences originate from Numts and not from mtDNA. Not only do the clade 1 sequences generated from ND4 PCR products (Figure 3) correspond to the sequences removed by mtDNA enrichment and the Numt on supercontig 1.1330, but clade 2 sequences were also found in all individuals that generated clade 1 sequences. The clade 2 sequences generated from ND4 PCR products (Figure 3) correspond to the reported mtDNA sequence. However, it is possible that not all of these clade 2 sequences have arisen from mtDNA as the data indicate that there are also very recently derived Numts: firstly, the putative Numts on supercontigs 1.363 and 1.593 are almost identical to the mtDNA sequences; and secondly, the clade 2 sequences of cloned PCR products have high numbers of mutations (Figure 5). Although it could be argued that these findings are the result of misassembly of the nuclear genome sequence and high Taq error rate, respectively, these do not appear to be the most parsimonious arguments.
The sites that differ between our clade 1 and clade 2 sequences are the same as those found in a large number of other studies [37–41, 56] using the putative ND4 mtDNA gene. In this and previous studies, the primers used were designed using mtDNA sequences from related taxa (Ae. albopictus). It has been noted that this strategy is inherently prone to the co-amplifcation of Numts . Many of the previous studies e.g. [37–39, 41] relied on single-strand conformation polymorphism (SSCP) coupled with sequencing of a few individuals. If co-amplifying Numt sequences were present and formed additional bands on the gene they may have been disregarded, as SSCP gels can contain background bands due to the same sequence adopting alternative conformations. This may explain why, besides this one, there has only been one other report  of heterogenous sequences using this marker. In this latter case the authors attributed the heterogenous sequences to heteroplasmy because the mixed sequences contained the variation found in their clade 1 and clade 2 sequences which they believed to all be true mtDNA.
There are at least two possible explanations for the origin of the clade 1 and clade 2 type Numts. Firstly, the Numts could have originated from the same species at different time points with the clade 1 sequences being entirely due to Numts and transferring at an earlier time. The deep clade structure would therefore have been formed by the movement of clade-1 like mtDNA into the nuclear genome, presumably in Africa of the order of one million years ago (making the assumption that divergence rates will have been similar to those of mtDNA estimated at 2.3% per million years ). This might at first appear inconsistent with the lack of amplification of other recently derived Numts, particularly of sequences intermediate between clades 1 and 2 that might be expected from continual mtDNA transfers to the nucleus. However, the ease with which clade 1 Numts are amplified could be explained by their being present in multiple copies in the genome, for example, due to secondary integration as suggested above. The one coding and 16 silent substitutions separating clade 1 and clade 2 sequences indicate the action of purifying selection. These substitutions would therefore have accumulated primarily along the clade 2 mtDNA lineage. This is consistent with the clade 1 sequences being multi-copy Numts as their independently acquired substitutions would not be apparent in sequenced PCR products. The star-like structure of the clade 1 sequences in the haplotype network is, however, hard to explain under this scenario.
A second possible explanation is that both clade 1 and clade 2 Numts are of recent origin but have originated from different source populations/species with divergent mtDNAs. A likely (but not only) candidate for a source of divergent mtDNA is the forest form of Ae. aegypti from Africa. The domestic form of Ae. aegypti studied here derives from Africa where there is known to be ongoing gene flow between it and the forest form . Putative mtDNA ND4 haplotypes of the forest form are shared with those of the domestic form in West Africa  and (from a comparison with our sequences) these comprise both clade 1 and clade 2 sequences. However, it is possible that the forest and domestic form have recently made secondary contact following a period of allopatry during which their mtDNA could have diverged. Clade 1 and clade 2 Numts could have arisen pre or post secondary contact and in either the forest and/or the domestic form.
These alternative hypotheses for the origin of the clade 1 and clade 2 Numts (single source at different times versus different sources and recent origin) could be distinguished by sequencing enriched mtDNA from a large number of individuals, ideally including those of the forest form. Under the second hypothesis at least some individuals should have clade 1 type mtDNA. This should be performed in conjunction with amplifications of the total genomic DNA using clade specific primers to determine if clade 1 or clade 2 Numt sequences are present.
Implications for inference of population history and genetic structure
The unknowing inclusion of Numts sequences in mtDNA based analyses could lead to mistaken inferences of population structure and population history. Since Ae. aegypti mosquitoes have spread from Africa to the rest of the Tropics relatively recently, the deep clade structure found in previous studies has often been interpreted to suggest colonization from at least two different source populations (and/or by many individuals) [38, 40, 46, 60, 61]. However, the deep clade structure could instead be due to the presence of Numts. In this case, the two deep clades could also mistakenly be interpreted to indicate the presence of two distinct species outside of Africa. This exemplifies the problem pointed out by Song et al.  that the use of mtDNA genes for DNA barcoding could lead to the overestimation of species numbers.
Population pairwise F ST values estimated using clade 1 and clade 2 sequences
northwest Thailand (Chiang Mai)
northeast Thailand (Ubon Ratchathani)
northwest Thailand (Chiang Mai)
Clade 2 (mtDNA-like)
northeast Thailand (Ubon Ratchathani)
Clade 2 (mtDNA-like)
Clade 2 (mtDNA-like)
There are a number of means by which Numt co-amplification can be limited such as mitochondrial purification before DNA extraction, long PCR, reverse transcriptase PCR, or using tissue rich in mtDNA (e.g. muscle) but none are guaranteed to absolutely overcome the problem ( and references therein). Perhaps the simplest solution to the problem providing that the Numts are monophyletic with respect to the mtDNA sequences, is their removal by restriction enzyme digestion or their avoidance by the use of mtDNA-specific primers . However, these possible solutions will not work here if, as the BLAST results and sequences of clones indicate, there are very recently derived Numt sequences. We therefore advocate that future studies of population structure in Ae. aegypti use other markers. With the genome sequence available there are a large number of nuclear markers to choose from.
Using a variety of approaches we conclude that Numts are prevalent in Ae. aegypti and that Numt evolution in Ae. aegypti is a very complex and dynamic process. These data also show that the deep clade structure we and others have been obtaining when performing supposedly mtDNA-based studies in Ae. aegypti may in fact be due to Numts. Some interpretations of the population and colonisation history of Ae. aegypti made to date may therefore be unreliable. We advocate that future studies of genetic population structure and gene flow in Ae. aegypti avoid using mtDNA data and use nuclear markers instead.
BLASTN searches of Ae. aegypti genome sequence
This study utilized the genomic sequence of Ae. aegypti (version-AaegL1) sequenced from the Liverpool LVP strain to approximately 8× coverage by the Broad Institute and The Institute for Genomic Research (TIGR). The genome sequence has over 1.31 billion bp and comprises 4,758 supercontigs but these have not yet been physically mapped to chromosomes http://aaegypti.vectorbase.org/Genome/Home/.
This genomic sequence was searched for Numts by BLASTN  on VectorBase for each of the 13 protein-coding genes, 22 tRNA genes and 2 rRNA genes from the Ae. aegypti mitochondrial genome sequence of 16,655 bp (GenBank accession number: EU352212). Threshold levels for the inference of Numts from BLASTN hits were taken as expectation values (E values) of 10-4 or 10-14, as used in other studies [30, 32, 64]. We determined the extent to which the hits against single genes were contiguous and due to single large Numts by carrying out alignments of the mtDNA genome with the regions of supercontigs containing multiple hits using ClustalX (1.83) [65, 66]. The aligned sequences were visualized in MEGA ver. 4  to determine the lengths of Numt sequences and to identify the Numt-supercontig junction points. The same software was also used for the translation of protein coding sequences to identify missense mutations and stop codons.
Extraction of genomic DNA, PCR amplication and direct sequencing
DNA was extracted from individual mosquitoes using a standard phenol/chloroform method . The final DNA pellet was suspended in 20 μl water (dH2O) and was diluted 1:20 to make a working solution. A fragment of the nicotinamide adenine dinucleotide dehydrogenase (NADH) subunit 4 mitochondrial (ND4) gene was amplified using the primers: ND4F (5'-GTTTAGATATARTTTCTTAYGG-3') and ND4R (5'-CTTCGDCTTCCWADWCGTTC-3'). The ND4R primer was the same as that used in previous studies in Mexico [37, 49], Thailand  and Venezuela . The ND4F primer used here was designed in Primer3 ver. 0.4.0  from an alignment of the ND4 region of the mtDNA of Aedes albopictus (GenBank accession number: AY072044), Anopheles gambiae (GenBank accession number: L20934) and Anopheles quadrimaculatus (GenBank accession number: L04272). The combination of primers used here therefore amplified a fragment that had the same 3' end as that of the previous studies but which was longer, being 763 bp rather than 359 bp in length.
The concentrations of the PCR reactants were 1 × ammonium buffer (CLP, Northampton, UK), 2.5 mM MgCl2, 200 μM dNTP and 0.4 μM primers. Amplifications were carried out in 50 μl volumes using 1 μl of template DNA. The PCR programme was modified from the previous Thailand and Mexico studies [37, 38, 49]. In an intial hot start, the reagents were heated to 95°C for 2 mins prior to the addition of 1.25 units of Thermoprime plus DNA polymerase (CLP, Northampton, UK). This was followed by 38 amplification cycles consisting of 30 sec at 92°C, 1 min at 50°C and 40 sec at 72°C, followed by a final extension for 5 mins at 72°C. PCR reactions were carried out on a GeneAmp® PCR System 9700 thermocycler (Applied Biosystems, Warrington, UK). PCR products were purified using Montage columns (Millipore, Billerica, MA, USA) and sequenced in both directions (Macrogen Inc., Seoul, Korea). DNA sequences were assembled using the Sequencher multiple sequence editor program ver. 4.5 (Gene Codes Corporation, Ann Arbor, USA) and checked manually. A total length of 763 bp of ND4 sequence was obtained from each individual.
Median joining networks were constructed using DNA Alignment ver. 188.8.131.52 and Network ver. 184.108.40.206 http://www.fluxus-technology.com. Genetic differentiation between pairs of populations (FST values) was estimated from the sequence data using analysis of molecular variance (AMOVA) in ARLEQUIN version 3.01 . Significance was estimated by 1000 permuations of the haplotypes between the populations.
Cloning and sequencing of PCR products
PCR products of the ND4 gene from individuals that yielded a mixed clade 1 and clade 2 sequence were cloned using the pGEM-T Vector System II kit according to the manufacturer's instruction (Promega Corporation, Madison, WI, USA). The same PCR primers were used to amplify the cloned ND4 gene fragments from colonies containing inserts. Colony material was prepared for PCR by heating a small portion of an individual colony in 20 μl dH2O to 95°C for 5 mins and 0.5 μl of this was used in PCR reaction volumes of 25 μl. The amplifications were performed using the same reactant concentrations as above. The PCR conditions were initial denaturation at 94°C for 5 mins followed by 30 cycles of 92°C for 30 sec, 52°C for 1 min and 72°C for 40 sec, with a final extension at 72°C for 7 mins. The PCR products were purified and sequenced and the same procedures as outlined above were used to check sequences and to construct a median joining network.
Fresh mtDNA separation and PCR amplification
Laboratory-reared freshly killed Ae. aegypti mosquitoes originating from Thailand were used for the isolation and extraction of enriched mtDNA using an alkaline lysis method . For each individual, DNA was extracted using a standard phenol/chloroform DNA extraction method  from both the starting homogenate (which contained both nuclear and mtDNA) and from the final mtDNA-rich supernatant. PCR amplification and sequencing of the ND4 region was carried out as above.
Additional PCR reactions
Primers designed and used for the specific amplification of clade 2 ND4 gene sequences and for the amplification of Numt junctions are given in Table 2 (see main text for details of their application). PCR was performed as before with the exception that no hot start was used and that the annealing temperature was increased to 55°C to achieve high specificity.
We are grateful to Mr. Sein Thaung (DMR-LM, Yangon, Myanmar), and all the entomology field staff from Thailand, Myanmar and Cambodia for their invaluable assistance and help in the field mosquito sample collections. We also thank Douda Bensasson and Casey Bergman (University of Manchester) for helpful discussions. This study was funded by the Special Programme for Research and Training in Tropical Diseases, World Health Organization (WHO/TDR) Collaborative Research Project Grant ID-A40198 and Research Training Grant (RTG) ID-A60987.
- Brown MW: Polymorphism in mitochondrial DNA of humans as revealed by restriction endonuclease analysis. Proceedings of the National Academy of Sciences of the United States of America. 1980, 77 (6): 3605-3609.PubMed CentralView ArticlePubMedGoogle Scholar
- Hale L, Singh R: Mitochondrial DNA variation and genetic structure in populations of Drosophila melanogaster. Molecular Biology Evolution. 1987, 4 (6): 622-637.PubMedGoogle Scholar
- Avise J, Ball R, Arnold J: Current versus historical population sizes in vertebrate species with high gene flow: a comparison based on mitochondrial DNA lineages and inbreeding theory for neutral mutations. Mol Biol Evol. 1988, 5 (4): 331-344.PubMedGoogle Scholar
- Guillen AKZ, Barrett GM, Takenaka O: Genetic diversity among African great apes based on mitochondrial DNA sequences. Biodiversity and Conservation. 2005, 14 (9): 2221-2233.View ArticleGoogle Scholar
- Avise JC, Arnold J, Ball RM, Bermingham E, Lamb T, Neigel JE, Reeb CA, Saunders NC: Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics. Annual Review of Ecology and Systematics. 1987, 18 (1): 489-522.View ArticleGoogle Scholar
- Meixner M, Arias M, Sheppard W: Mitochondrial DNA polymorphisms in honey bee subspecies from Kenya. Apidologie. 2000, 31: 181-190.View ArticleGoogle Scholar
- Shoemaker DD, Ahrens ME, Ross KG: Molecular phylogeny of fire ants of the Solenopsis saevissima species-group based on mtDNA sequences. Molecular Phylogenetics and Evolution. 2006, 38 (1): 200-215.View ArticlePubMedGoogle Scholar
- Michel AP, Ingrasci MJ, Schemerhorn BJ, Kern M, Goff G, Coetzee M, Elissa N, Fontenille D, Vulule J, Lehmann T, et al: Rangewide population genetic structure of the African malaria vector Anopheles funestus. Molecular Ecology. 2005, 14 (14): 4235-4248.View ArticlePubMedGoogle Scholar
- Besansky NJ, Lehmann T, Fahey GT, Fontenille D, Braack LEO, Hawley WA, Collins FH: Patterns of mitochondrial variation within and between African malaria vectors, Anopheles gambiae and Anopheles arabiensis, suggest extensive gene flow. Genetics. 1997, 147 (4): 1817-1828.PubMed CentralPubMedGoogle Scholar
- Walton C, Handley JM, Tun-Lin W, Collins FH, Harbach RE, Baimai V, Butlin RK: Population structure and population history of Anopheles dirus mosquitoes in Southeast Asia. Molecular Biology and Evolution. 2000, 17 (6): 962-974.View ArticlePubMedGoogle Scholar
- Birungi J, Munstermann LE: Genetic structure of Aedes albopictus (Diptera: Culicidae) populations based on mitochondrial ND5 sequences: evidence for an independent invasion into Brazil and United States. Annals of the Entomological Society of America. 2002, 95: 125-132.View ArticleGoogle Scholar
- Cook S, Diallo M, Sall AA, Cooper A, Holmes EC: Mitochondrial markers for molecular identification of Aedes mosquitoes (Diptera: Culicidae) involved in transmission of arboviral disease in West Africa. Journal of Medical Entomology. 2005, 42: 19-28.View ArticlePubMedGoogle Scholar
- Remme JHF, Blas E, Chitsulo L, Desjeux PMP, Engers HD, Kanyok TP, Kayondo JFK, Kioy DW, Kumaraswami V, Lazdins JK, et al: Strategic emphases for tropical diseases research: a TDR perspective. Trends in Parasitology. 2002, 18 (10): 421-426.View ArticlePubMedGoogle Scholar
- WHO: Basic and strategic research, molecular entomology committee workplan. 2004, [http://www.who.int/tdr/disease/dengue/direction.htm]Google Scholar
- Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. PLoS Biology. 2003, 270 (1512): 313-321.Google Scholar
- Bazin E, Glemin S, Galtier N: Population size does not influence mitochondrial genetic diversity in animals. Science. 2006, 312 (5773): 570-572.View ArticlePubMedGoogle Scholar
- Ballard JWO, Whitlock MC: The incomplete natural history of mitochondria. Molecular Ecology. 2004, 13 (4): 729-744.View ArticlePubMedGoogle Scholar
- Hurst GDD, Jiggins FM: Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. PLoS Biology. 2005, 272 (1572): 1525-1534.Google Scholar
- Lopez JV, Yuhki N, Masuda R, Modi W, O'Brien SJ: Numt, a recent transfer and tandem amplification of mitochondrial DNA to the nuclear genome of the domestic cat. Journal of Molecular Evolution. 1994, 39 (2): 174-190.PubMedGoogle Scholar
- Zhang D-X, Hewitt G: Nuclear integrations: challenges for mitochondrial DNA markers. Trends in Ecology and Evolution. 1996, 11: 247-251.View ArticlePubMedGoogle Scholar
- Sorenson M, Quinn T: Numts: a challenge for avian systematics and population biology. The Auk. 1998, 115: 214-221.View ArticleGoogle Scholar
- Bensasson D, Zhang D-X, Hartl DL, Hewitt GM: Mitochondrial pseudogenes: evolution's misplaced witnesses. Trends in Ecology & Evolution. 2001, 16 (6): 314-321.View ArticleGoogle Scholar
- Arctander P, Fjeldsa J: Andean tapaculos of the genus Scytalopus (Aves, Rhinocryptidae); a study of speciation using DNA sequence data. Experientia Basel Supplementum. 1994, 68: 205-225.Google Scholar
- Arctander P: Comparison of a mitochondrial gene and a corresponding nuclear pseudogene. Proceedings of the Royal Society of London B. 1995, 262: 13-19.View ArticleGoogle Scholar
- Jensen-Seaman MI, Esteban E Sarmiento, Amos S Deinard, Kenneth K Kidd: Nuclear integrations of mitochondrial DNA in gorillas. American Journal of Primatology. 2004, 63 (3): 139-147.View ArticlePubMedGoogle Scholar
- Thalmann O, Hebler J, Poinar HN, Paabo S, Vigilant L: Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes. Molecular Ecology. 2004, 13 (2): 321-335.View ArticlePubMedGoogle Scholar
- Thalmann O, Serre D, Hofreiter M, Lukas D, Eriksson J, Vigilant L: Nuclear insertions help and hinder inference of the evolutionary history of gorilla mtDNA. Molecular Ecology. 2005, 14 (1): 179-188.View ArticlePubMedGoogle Scholar
- Antunes A, Pontius J, Ramos MJ, O'Brien SJ, Johnson WE: Mitochondrial Introgressions into the nuclear genome of the domestic cat. J Hered. 2007, 98 (5): 414-420.View ArticlePubMedGoogle Scholar
- Pereira S, Baker A: Low number of mitochondrial pseudogenes in the chicken (Gallus gallus) nuclear genome: implications for molecular inference of population history and phylogenetics. BMC Evolutionary Biology. 2004, 4 (1): 17-PubMed CentralView ArticlePubMedGoogle Scholar
- Pamilo P, Viljakainen L, Vihavainen A: Exceptionally high density of Numts in the honeybee genome. Molecular Biology and Evolution. 2007, 24 (6): 1340-1346.View ArticlePubMedGoogle 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 Molecular Biology. 2007, 16 (6): 777-784.View ArticlePubMedGoogle Scholar
- Richly E, Leister D: Numts in sequenced eukaryotic genomes. Molecular Biology and Evolution. 2004, 21: 1081-1084.View ArticlePubMedGoogle Scholar
- Bensasson D, Petrov DA, Zhang D-X, Hartl DL, Hewitt GM: Genomic gigantism: DNA loss is slow in mountain grasshoppers. Mol Biol Evol. 2001, 18 (2): 246-253.View ArticlePubMedGoogle Scholar
- Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, et al: Genome Sequence of Aedes aegypti, a Major Arbovirus Vector. Science. 2007, 316 (5832): 1718-1723.View ArticlePubMedGoogle Scholar
- Woischnik M, Moraes C: Pattern of organization of human mitochondrial pseudogenes in the nuclear genome. Genome Research. 2002, 12 (6): 885-893.PubMed CentralView ArticlePubMedGoogle Scholar
- Tourmen Y, Baris O, Dessen P, Jacques C, Malthiery Y, Reynier P: Structure and chromosomal distribution of human mitochondrial pseudogenes. Genomics. 2002, 80 (1): 71-77.View ArticlePubMedGoogle Scholar
- Gorrochotegui-Escalante N, Gomez-Machorro C, Lozano-Fuentes S, Fernandez-Salas I, Munoz MD, Farfan-Ale JA, Garcia-Rejon J, Beaty BJ, Black WC: Breeding structure of Aedes aegypti populations in Mexico varies by region. American Journal of Tropical Medicine and Hygiene. 2002, 66 (2): 213-222.PubMedGoogle Scholar
- Bosio CF, Harrington LC, Jones JW, Sithiprasasna R, Norris DE, Scott TW: Genetic structure of Aedes aegypti populations in Thailand using mitochondrial DNA. American Journal of Tropical Medicine and Hygiene. 2005, 72 (4): 434-442.PubMedGoogle Scholar
- Costa-da-Silva AL, Capurro ML, Bracco JE: Genetic lineages in the yellow fever mosquito Aedes aegypti (Diptera: Culicidae) from Peru. Mem Inst Oswaldo Cruz, Rio de Janeiro. 2005, 100 (6): 639-644.View ArticleGoogle Scholar
- Bracco JE, Capurro ML, Lourenco-de-Oliveira R, Sallum MAM: Genetic variability of Aedes aegypti in the Americas using a mitochondrial gene: evidence of multiple introductions. Mem Inst Oswaldo Cruz, Rio de Janeiro. 2007, 102 (5): ISSN 0074-0276Google Scholar
- Herrera F, Urdaneta L, Rivero J, Zoghbi N, Ruiz J, Carrasquel G, Martínez JA, Pernalete M, Villegas P, Montoya A, et al: Population genetic structure of the dengue mosquito Aedes aegypti in Venezuela. Memórias do Instituto Oswaldo Cruz, Rio de Janeiro. 2006, 101: 625-633.View ArticleGoogle Scholar
- Smith CEG: The history of dengue in tropical Asia and its probable relationship to the mosquito Aedes aegypti. Journal of Tropical Medicine and Hygiene. 1956, 59: 243-251.PubMedGoogle Scholar
- Tabachnick WJ, Powell JR: Worldwide survey of genetic variation in the yellow fever mosquito, Aedes aegypti. Genetical Research. 1979, 34 (3): 215-229.View ArticlePubMedGoogle Scholar
- Tien TK, Vazeille-Falcoz M, Mousson L, Hoang TH, Rodhain F, Huong NT, Failloux AB: Aedes aegypti in Ho Chi Minh City (Viet Nam): susceptibility to dengue 2 virus and genetic differentiation. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1999, 93 (6): 581-586.View ArticleGoogle Scholar
- Mousson L, Vazeille M, Chawprom S, Prajakwong S, Rodhain F, Failloux AB: Genetic structure of Aedes aegypti populations in Chiang Mai (Thailand) and relation with dengue transmission. Tropical Medicine & International Health. 2002, 7 (10): 865-872.View ArticleGoogle Scholar
- Failloux AB, Vazeille M, Rodhain F: Geographic genetic variation in populations of the dengue virus vector Aedes aegypti. Journal of Molecular Evolution. 2002, 55 (6): 653-663.View ArticlePubMedGoogle Scholar
- Paduan KDS, Ribolla PEM: Mitochondrial DNA polymorphism and heteroplasmy in populations of Aedes aegypti in Brazil. Journal of Medical Entomology. 2008, 59-67.Google Scholar
- Parr R, Maki J, Reguly B, Dakubo G, Aguirre A, Wittock R, Robinson K, Jakupciak J, Thayer R: The pseudo-mitochondrial genome influences mistakes in heteroplasmy interpretation. BMC Genomics. 2006, 7 (1): 185-PubMed CentralView ArticlePubMedGoogle Scholar
- Gorrochotegui-Escalante N, Munoz MD, Fernandez-Salas I, Beaty BJ, Black WC: Genetic isolation by distance among Aedes aegypti populations along the northeastern coast of Mexico. American Journal of Tropical Medicine and Hygiene. 2000, 62 (2): 200-209.PubMedGoogle Scholar
- Beebe NW, Whelan PI, Hurk Avd, Ritchie SA, Cooper RD: Genetic diversity of the dengue vector Aedes aegypti in Australia and implications for future surveillance and mainland incursion monitoring. Commun Dis Intell. 2005, 29 (3): 299-304.Google Scholar
- Haag-Liautard C, Coffey N, Houle D, Lynch M, Charlesworth B, Keightley PD: Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster. PLoS Biology. 2008, 6 (8): e204-PubMed CentralView ArticlePubMedGoogle Scholar
- Kobayashi N, Tamura K, Aotsuka T: PCR error and molecular population genetics. Biochemical Genetics. 1999, 37 (9110): 317-321.View ArticlePubMedGoogle Scholar
- Triant D, DeWoody J: Molecular analyses of mitochondrial pseudogenes within the nuclear genome of arvicoline rodents. Genetica. 2008, 132 (1): 21-33.View ArticlePubMedGoogle Scholar
- Ricchetti M, Tekaia F, Dujon B: Continued colonization of the human genome by mitochondrial DNA. PLoS Biology. 2004, 2 (9): 1313-1323.View ArticleGoogle Scholar
- Hazkani-Covo E, Graur D: A comparative analysis of Numt evolution in human and chimpazee. Molecular Biology and Evolution. 2007, 24 (1): 13-18.View ArticlePubMedGoogle Scholar
- Urdaneta-Marquez L, Bosio C, Herrera F, Rubio-Palis Y, Salasek M, Black W: Genetic relationships among Aedes aegypti collections in Venezuela as determined by mitochondrial DNA variation and nuclear single nucleotide polymorphisms. American Journal of Tropical Medicine and Hygiene. 2008, 78 (3): 479-491.PubMedGoogle Scholar
- Sorenson M, Fleischer R: Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus. Proceedings of the National Academy of Sciences of the United States of America. 1996, 93: 15239-15243.PubMed CentralView ArticlePubMedGoogle Scholar
- Brower AV: Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from patterns of mitochondrial DNA evolution. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91 (14): 6491-6495.PubMed CentralView ArticlePubMedGoogle Scholar
- Paupy C, Brengues C, Kamgang B, Herve J-P, Fontenille D, Simard F: Gene flow between domestic and sylvan populations of Aedes aegypti (Diptera: Culicidae) in North Cameroon. Journal of Medical Entomology. 2008, 45 (3): 391-400.View ArticlePubMedGoogle Scholar
- Huber K, Le Loan L, Chantha N, Failloux AB: Human transportation influences Aedes aegypti gene flow in Southeast Asia. Acta Tropica. 2004, 90 (1): 23-29.View ArticlePubMedGoogle Scholar
- Ravel S, Monteny N, Olmos DV, Verdugo JE, Cuny G: A preliminary study of the population genetics of Aedes aegypti (Diptera: Culicidae) from Mexico using microsatellite and AFLP markers. Acta Tropica. 2001, 78 (3): 241-250.View ArticlePubMedGoogle Scholar
- Song H, Buhay JE, Whiting MF, Crandall KA: Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proceedings of the National Academy of Sciences. 2008, 105 (36): 13486-13491.View ArticleGoogle Scholar
- Altschul S, Gish W, Miller W, Myers E, Lipman D: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Behura S: Analysis of nuclear copies of mitochondrial sequences in honey bee Apis mellifera genome. Mol Biol and Evol. 2007, 24 (7): 1492-1505.View ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research. 1997, 25: 4876-4882.PubMed CentralView ArticlePubMedGoogle Scholar
- Jeanmougin F, Thompson J, Gouy M, Higgins DG, Gibson TJ: Multiple sequence alignment with Clustal X. Trends in Biochemical Sciences. 1998, 23: 403-405.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Molecular Biology and Evolution. 2007, 24: 1596-1599.View ArticlePubMedGoogle Scholar
- Sambrook J, Russell DW: Molecular cloning: a laboratory manual. 2001, New York: Cold Spring Harbor Laboratory Press, 3Google Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for General Users and for Biologist Programmers. Methods in Molecular Biology, Bioinformatics Methods and Protocols. Edited by: Misener S, Krawetz SA. 1999, Totowa, NJ: Humana Press Inc, 132: 365-386.View ArticleGoogle Scholar
- Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16 (1): 37-48.View ArticlePubMedGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes – application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMed CentralPubMedGoogle Scholar
- Tamura K, Aotsuka T: Rapid isolation method of animal mitochondrial DNA by the alkaline lysis procedure. Biochemical Genetics. 1988, 26 (11): 815-819.View ArticlePubMedGoogle 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.