Polymorphic NumtS trace human population relationships

The human genome is constantly subjected to evolutionary forces which shape its architecture. Insertions of mitochondrial DNA sequences into nuclear genome (NumtS) have been described in several eukaryotic species, including Homo sapiens and other primates. The ongoing process of the generation of NumtS has made them valuable markers in primate phylogenetic studies, as well as potentially informative loci for reconstructing the genetic history of modern humans. Here, we report the identification of 53 human-specific NumtS by inspection of the UCSC genome browser, showing that they may be direct insertions of mitochondrial DNA into the human nuclear DNA after the human-chimpanzee split. In silico analyses allowed us to identify 14 NumtS which are polymorphic in terms of their presence/absence within the human genome in individuals of different ancestry. The allele frequencies of these polymorphic NumtS were calculated for 1000 Genomes Project sequence data from 13 populations worldwide, and principal components analysis and hierarchical clustering methods allowed the detection of strong signals of geographical structure related to the genetic diversity of these loci. All identified polymorphic human-specific NumtS together with a tandemly duplicated NumtS have also been validated by PCR amplification on a panel of 60 samples belonging to five native populations worldwide, confirming the expected NumtS variability. On the basis of these findings, we have succeeded in depicting the landscape of variation of a series of NumtS in several ethnic groups, making an advance in their identification as useful markers in the study on human population genetics.     

家马核基因组中线粒体核内插入序列分析

在各种真核生物核基因组中,存在一些由线粒体基因组转移进入核基因组中的DNA片段,这些被认为是分子化石的片段叫做线粒体核内插入序列(Numt)。由于Numt与真实的线粒体序列高度相似,因此它的存在必然会成为PCR扩增线粒体DNA的不利因素。利用已经公布的家马(Equus caballus)基因组序列(2007年9月公布,GenBank登录号为NC_009144-NC_009175)对家马Numt进行了深入分析,共发现200个可能的Numt,长度范围为29到3727bp,其中有10个的长度大于800bp。分析结果显示由于不存在线粒体控制区域的疑似Numt,因此对基于此区域的群体遗传学研究不会产生影响。本研究还发现在家马进化过程中,第1号和27号染色体更倾向于接受线粒体序列的转移。以上结果将为今后马科动物的研究提供重要的参考信息,有助于避免在线粒体DNA研究中由于Numt污染的存在而得出错误的实验结果。     

Identification of three populations ofOphiostoma ulmi (aggressive subgroup) by mitochondrial DNA restriction-site mapping and nuclear DNA fingerprinting

Three genetically distinct populations of the Dutch elm pathogenOphistoma ulmi within the aggressive subgroup were defined by the hybridization of a human minisatellite DNA sequence (HVR 33.6) to polymorphic dispersed DNA sequences within theO. ulmi nuclear genomes. For the 10 isolates examined there was a close correlation between nuclear DNA fingerprints and mitochondrial (mt) DNA restriction patterns. A restriction-site map was constructed for the mitochondrial genomes for each of these populations. The three mt DNA maps corresponded to genome sizes of 49.1 (Type I), 49.9 (Type II), and 53.9 (Type III) kilobase pairs (kbp) of DNA. The Type I and Type II mt genomes differed from the Type III mt genome by discrete length mutations of 4.8 and 4.0 kbp, respectively. It is unknown whether these length mutations resulted from insertions into or deletions from a progenitor mitochondrial genome. There was no correlation between the mitochondrial or nuclear genotypes and the geographical source of the isolates.     

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.     

Mitochondrial pseudogenes: evolution's misplaced witnesses

Nuclear copies of mitochondrial DNA (mtDNA) have contaminated PCR-based mitochondrial studies of over 64 different animal species. Since the last review of these nuclear mitochondrial pseudogenes (Numts) in animals, Numts have been found in 53 of the species studied. The recent evidence suggests that Numts are not equally abundant in all species, for example they are more common in plants than in animals, and also more numerous in humans than in Drosophila. Methods for avoiding Numts have now been tested, and several recent studies demonstrate the potential utility of Numt DNA sequences in evolutionary studies. As relics of ancient mtDNA, these pseudogenes can be used to infer ancestral states or root mitochondrial phylogenies. Where they are numerous and selectively unconstrained, Numts are ideal for the study of spontaneous mutation in nuclear genomes.     

Whole Mitochondrial and Plastid Genome SNP Analysis of Nine Date Palm Cultivars Reveals Plastid Heteroplasmy and Close Phylogenetic Relationships among Cultivars

Date palm is a very important crop in western Asia and northern Africa, and it is the oldest domesticated fruit tree with archaeological records dating back 5000 years. The huge economic value of this crop has generated considerable interest in breeding programs to enhance production of dates. One of the major limitations of these efforts is the uncertainty regarding the number of date palm cultivars, which are currently based on fruit shape, size, color, and taste. Whole mitochondrial and plastid genome sequences were utilized to examine single nucleotide polymorphisms (SNPs) of date palms to evaluate the efficacy of this approach for molecular characterization of cultivars. Mitochondrial and plastid genomes of nine Saudi Arabian cultivars were sequenced. For each species about 60 million 100 bp paired-end reads were generated from total genomic DNA using the Illumina HiSeq 2000 platform. For each cultivar, sequences were aligned separately to the published date palm plastid and mitochondrial reference genomes, and SNPs were identified. The results identified cultivar-specific SNPs for eight of the nine cultivars. Two previous SNP analyses of mitochondrial and plastid genomes identified substantial intra-cultivar ( = intra-varietal) polymorphisms in organellar genomes but these studies did not properly take into account the fact that nearly half of the plastid genome has been integrated into the mitochondrial genome. Filtering all sequencing reads that mapped to both organellar genomes nearly eliminated mitochondrial heteroplasmy but all plastid SNPs remained heteroplasmic. This investigation provides valuable insights into how to deal with interorganellar DNA transfer in performing SNP analyses from total genomic DNA. The results confirm recent suggestions that plastid heteroplasmy is much more common than previously thought. Finally, low levels of sequence variation in plastid and mitochondrial genomes argue for using nuclear SNPs for molecular characterization of date palm cultivars.     

Population‐specific dynamics and selection patterns of transposable element insertions in European natural populations

Transposable elements (TEs) are ubiquitous sequences in genomes of virtually all species. While TEs have been investigated for several decades, only recently we have the opportunity to study their genome‐wide population dynamics. Most of the studies so far have been restricted either to the analysis of the insertions annotated in the reference genome or to the analysis of a limited number of populations. Taking advantage of the European Drosophila population genomics consortium (DrosEU) sequencing data set, we have identified and measured the dynamics of TEs in a large sample of European Drosophila melanogaster natural populations. We showed that the mobilome landscape is population‐specific and highly diverse depending on the TE family. In contrast with previous studies based on SNP variants, no geographical structure was observed for TE abundance or TE divergence in European populations. We further identified de novo individual insertions using two available programs and, as expected, most of the insertions were present at low frequencies. Nevertheless, we identified a subset of TEs present at high frequencies and located in genomic regions with a high recombination rate. These TEs are candidates for being the target of positive selection, although neutral processes should be discarded before reaching any conclusion on the type of selection acting on them. Finally, parallel patterns of association between the frequency of TE insertions and several geographical and temporal variables were found between European and North American populations, suggesting that TEs can be potentially implicated in the adaptation of populations across continents.     

Continued colonization of the human genome by mitochondrial DNA

Integration of mitochondrial DNA fragments into nuclear chromosomes (giving rise to nuclear DNA sequences of mitochondrial origin, or NUMTs) is an ongoing process that shapes nuclear genomes. In yeast this process depends on double-strand-break repair. Since NUMTs lack amplification and specific integration mechanisms, they represent the prototype of exogenous insertions in the nucleus. From sequence analysis of the genome of Homo sapiens, followed by sampling humans from different ethnic backgrounds, and chimpanzees, we have identified 27 NUMTs that are specific to humans and must have colonized human chromosomes in the last 4–6 million years. Thus, we measured the fixation rate of NUMTs in the human genome. Six such NUMTs show insertion polymorphism and provide a useful set of DNA markers for human population genetics. We also found that during recent human evolution, Chromosomes 18 and Y have been more susceptible to colonization by NUMTs. Surprisingly, 23 out of 27 human-specific NUMTs are inserted in known or predicted genes, mainly in introns. Some individuals carry a NUMT insertion in a tumor-suppressor gene and in a putative angiogenesis inhibitor. Therefore in humans, but not in yeast, NUMT integrations preferentially target coding or regulatory sequences. This is indeed the case for novel insertions associated with human diseases and those driven by environmental insults. We thus propose a mutagenic phenomenon that may be responsible for a variety of genetic diseases in humans and suggest that genetic or environmental factors that increase the frequency of chromosome breaks provide the impetus for the continued colonization of the human genome by mitochondrial DNA.     

Extensive mitochondrial DNA transfer in a rapidly evolving rodent has been mediated by independent insertion events and by duplications

Mitochondrial DNA translocations to the nucleus (numt pseudogenes) are pervasive among eukaryotes, but copy number within the nuclear genome varies widely among taxa. As an increasing number of genomes are sequenced in their entirety, the origins, transfer mechanisms and insertion sites of numts are slowly being characterized. We investigated mitochondrial transfers within a genetically diverse rodent lineage and here report 15 numts totaling 21.8 kb that are harbored within the nuclear genome of the vole Microtus rossiaemeridionalis. The 15 numts total 21.8 kb and range from 0.39 to over 3.0 kb in length. Phylogenetic analyses revealed that these numts resulted from three independent insertions to the nucleus, two of which were followed by subsequent nuclear duplication events. The dates of the two translocations that led to subsequent duplications were estimated at 1.97 and 1.19 MYA, which coincide with the origin and radiation of the genus Microtus. Numt sequence data from five Microtus species were used to estimate an average rate of nucleotide substitution as 2.6x10(-8) subs/site/yr. This substitution rate is higher than in many other mammals, but is concordant with the elevated rate of mtDNA substitution in this lineage. Our data suggest that numt translocation in Microtus is more extensive than in either Mus or in Rattus, consistent with the elevated rate of speciation, karyotypic rearrangement, and mitochondrial DNA evolution in Microtus.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Assembly-free quantification of vagrant DNA inserts

Inserts of DNA from extranuclear sources, such as organelles and microbes, are common in eukaryote nuclear genomes. However, sequence similarity between the nuclear and extranuclear DNA, and a history of multiple insertions, make the assembly of these regions challenging. Consequently, the number, sequence and location of these vagrant DNAs cannot be reliably inferred from the genome assemblies of most organisms. We introduce two statistical methods to estimate the abundance of nuclear inserts even in the absence of a nuclear genome assembly. The first (intercept method) only requires low-coverage (<1×) sequencing data, as commonly generated for population studies of organellar and ribosomal DNAs. The second method additionally requires that a subset of the individuals carry extranuclear DNA with diverged genotypes. We validated our intercept method using simulations and by re-estimating the frequency of human NUMTs (nuclear mitochondrial inserts). We then applied it to the grasshopper Podisma pedestris, exceptional for both its large genome size and reports of numerous NUMT inserts, estimating that NUMTs make up 0.056% of the nuclear genome, equivalent to >500 times the mitochondrial genome size. We also re-analysed a museomics data set of the parrot Psephotellus varius, obtaining an estimate of only 0.0043%, in line with reports from other species of bird. Our study demonstrates the utility of low-coverage high-throughput sequencing data for the quantification of nuclear vagrant DNAs. Beyond quantifying organellar inserts, these methods could also be used on endosymbiont-derived sequences. We provide an R implementation of our methods called “vagrantDNA” and code to simulate test data sets.     

Characterization of full-length LINE-1 insertions in 154 genomes

Long interspersed nuclear elements (LINEs) are retrotransposons that contribute to genetic variation in the human genome. LINE-1 elements in larger-scale studies are challenging to identify using sequencing technologies due to cost and scalability. We developed an approach using optical mapping for detection of full-length LINE-1 insertions and 10× sequencing for confirmation. We found 51 true positive full-length LINE-1 insertions, of which 4 are novel insertions, in NA12878. Repeating our analysis on a larger sample set representing 26 populations, we identified 329 full-length LINE-1 elements, of which 123 are novel. 24.8% of these 329 LINE-1 insertions were shared amongst all 5 superpopulations (AFR, AMR, EUR, EAS, SAS). The African superpopulation has a higher percentage of population-specific LINE-1 insertions than any other superpopulation. These data indicate that our approach can provide high-speed, cost-effective, and increased accuracy for LINE-1 detection. These data also provide an insight into variations of LINE-1 elements between different populations.     

Fugu genome does not contain mitochondrial pseudogenes

Contrary to previous observations that fish genomes are devoid of nuclear mitochondrial pseudogenes, a genome-wide survey identified a large number of "recent" and "ancient" nuclear mitochondrial DNA fragments (Numts) in the whole-genome sequences of the fugu (Takifugu rubripes), Tetraodon nigroviridis, and zebrafish (Danio rerio). We have analyzed the latest assembly (v4.0) of the fugu genome and show that, like the Anopheles genome, the fugu nuclear genome does not contain mitochondrial pseudogenes. Fugu assembly v4.0 contains a single scaffold representing the near complete sequence of the fugu mitochondria. The "recent" Numts identified by the previous study in fugu assembly v2.0 are in fact shotgun sequences of mitochondrial DNA that were misassembled with the nuclear sequences, whereas the "ancient" Numts appear to be the result of spurious matches. It is likely that the Numts identified in the genomes of Tetraodon and zebrafish are also similar artifacts. Shotgun sequences of whole genomes often include some mitochondrial sequences. Therefore, any Numts identified in shotgun-sequence assemblies should be verified by Southern hybridization or PCR amplification.     

RADseq as a valuable tool for plants with large genomes—A case study in cycads

Full genome sequencing of organisms with large and complex genomes is intractable and cost ineffective under most research budgets. Cycads (Cycadales) represent one of the oldest lineages of the extant seed plants and, partly due to their age, have incredibly large genomes up to ~60 Gbp. Restriction site‐associated DNA sequencing (RADseq) offers an approach to find genome‐wide informative markers and has proven to be effective with both model and nonmodel organisms. We tested the application of RADseq using ezRAD across all 10 genera of the Cycadales including an example data set of Cycas calcicola representing 72 samples from natural populations. Using previously available plastid and mitochondrial genomes as references, reads were mapped recovering plastid and mitochondrial genome regions and nuclear markers for all of the genera. De novo assembly generated up to 138,407 high‐depth clusters and up to 1,705 phylogenetically informative loci for the genera, and 4,421 loci for the example assembly of C. calcicola. The number of loci recovered by de novo assembly was lower than previous RADseq studies, yet still sufficient for downstream analysis. However, the number of markers could be increased by relaxing our assembly parameters, especially for the C. calcicola data set. Our results demonstrate the successful application of RADseq across the Cycadales to generate a large number of markers for all genomic compartments, despite the large number of plastids present in a typical plant cell. Our modified protocol was adapted to be applied to cycads and other organisms with large genomes to yield many informative genome‐wide markers.     

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747‐Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA‐sequencing, RAD‐sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA‐sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO‐lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage‐specific selection on mitochondria‐related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.     

Promiscuous mitochondrial group II intron sequences in plant nuclear genomes

Gene translocations from the organelles to the nucleus are postulated by the endosymbiont hypothesis. We here report evidence for sequence insertions in the nuclear genomes of plants that are derived from noncoding regions of the mitochondrial genome. Fragments of mitochondrial group II introns are identified in the nuclear genomes of tobacco and a bean species. The duplicated intron sequences of 75–140 bp are derived from cis- and trans-splicing introns of genes encoding subunits 1 and 5 of the NADH dehydrogenase. The mitochondrial sequences are inserted in the vicinities of a lectin gene, different glucanase genes and a gene encoding a subunit of photosystem II. Sequence similarities between the nuclear and mitochondrial copies are in the range of 80 to 97%, suggesting recent transfer events that occurred in the basic glucanase genes before and in the lectin gene after the gene duplications in the evolution of the nuclear gene families. Overlapping regions of the same introns are in two instances also involved in intramitochondrial sequence duplications. Correspondence to: V. Knoop     

Genome‐wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high‐throughput sequencing have promoted the collection of reference genomes and genome‐wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single‐nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population‐specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.     

Sequencing of pooled DNA samples (Pool-Seq) uncovers complex dynamics of transposable element insertions in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that parasitize genomes by semi-autonomously increasing their own copy number within the host genome. While TEs are important for genome evolution, appropriate methods for performing unbiased genome-wide surveys of TE variation in natural populations have been lacking. Here, we describe a novel and cost-effective approach for estimating population frequencies of TE insertions using paired-end Illumina reads from a pooled population sample. Importantly, the method treats insertions present in and absent from the reference genome identically, allowing unbiased TE population frequency estimates. We apply this method to data from a natural Drosophila melanogaster population from Portugal. Consistent with previous reports, we show that low recombining genomic regions harbor more TE insertions and maintain insertions at higher frequencies than do high recombining regions. We conservatively estimate that there are almost twice as many "novel" TE insertion sites as sites known from the reference sequence in our population sample (6,824 novel versus 3,639 reference sites, with on average a 31-fold coverage per insertion site). Different families of transposable elements show large differences in their insertion densities and population frequencies. Our analyses suggest that the history of TE activity significantly contributes to this pattern, with recently active families segregating at lower frequencies than those active in the more distant past. Finally, using our high-resolution TE abundance measurements, we identified 13 candidate positively selected TE insertions based on their high population frequencies and on low Tajima's D values in their neighborhoods.     

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA‐originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last . A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415 bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein‐coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.     

Hybridisation,paternal leakage and mitochondrial DNA linearization in three anomalous fish (Scombridae)

Using mitochondrial DNA for species identification and population studies assumes that the genome is maternally inherited, circular, located in the cytoplasm and lacks recombination. This study explores the mitochondrial genomes of three anomalous mackerel. Complete mitochondrial genome sequencing plus nuclear microsatellite genotyping of these fish identified them as Scomberomorus munroi (spotted mackerel). Unlike normal S. munroi, these three fish also contained different linear, mitochondrial genomes of Scomberomorus semifasciatus (grey mackerel). The results are best explained by hybridisation, paternal leakage and mitochondrial DNA linearization. This unusual observation may provide an explanation for mtDNA outliers in animal population studies.