Skip to main content
Download PDF

Exceptionally high cumulative percentage of NUMTs originating from linear mitochondrial DNA molecules in the Hydra magnipapillata genome

BMC Genomics volume 14, Article number: 447 (2013) Cite this article

Abstract

Background

In contrast to most animal genomes, mitochondrial genomes in species belonging to the phylum Cnidaria show distinct variations in genome structure, including the mtDNA structure (linear or circular) and the presence or absence of introns in protein-coding genes. Therefore, the analysis of nuclear insertions of mitochondrial sequences (NUMTs) in cnidarians allows us to compare the NUMT content in animals with different mitochondrial genome structures.

Results

NUMT identification in the Hydra magnipapillata, Nematostella vectensis and Acropora digitifera genomes showed that the NUMT density in the H. magnipapillata genome clearly exceeds that in other two cnidarians with circular mitochondrial genomes. We found that H. magnipapillata is an exceptional ancestral metazoan with a high NUMT cumulative percentage but a large genome, and its mitochondrial genome linearisation might be responsible for the NUMT enrichment. We also detected the co-transposition of exonic and intronic fragments within NUMTs in N. vectensis and provided direct evidence that mitochondrial sequences can be transposed into the nuclear genome through DNA-mediated fragment transfer. In addition, NUMT expression analyses showed that NUMTs are co-expressed with adjacent protein-coding genes, suggesting the relevance of their biological function.

Conclusions

Taken together, our results provide valuable information for understanding the impact of mitochondrial genome structure on the interaction of mitochondrial molecules and nuclear genomes.

Background

In eukaryotes, mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs) [1]. Du Buy and Riley identified NUMTs in a wide range of species, including plants, yeast, alveolates, nematodes, insects, and vertebrates [2–5]. Although NUMTs are present in numerous species, the numbers of these sequences widely vary. Some species (Monodelphis domestica) possess up to 2.04 Mbp of NUMTs, whereas no NUMTs have been detected in other species (e.g., Anopheles gambiae, Branchiostoma floridae, Danio rerio, Ciona savignyi) [4, 5]. To date, the highest cumulative percentage of NUMTs is 0.0861% in the Apis mellifera genome [5, 6]. NUMTs accumulate in the genomes over a continuous evolutionary process [7–11]. In general, each NUMT originates in one of two ways: the independent insertion of DNA from mitochondrial DNA into nuclear chromosomes and duplication after the insertion into nuclear DNA [9, 12–15]. Previous studies in human and honeybee genomes have suggested that only one-third of NUMTs were integrated as independent mitochondrial sequences, whereas the remaining two-thirds of NUMTs arose from duplications after insertion into the nuclear genome [6, 8]. Transposable elements or short-dispersed repeats have been associated with the on-going integration of mtDNA sequences into the nuclear genome and their subsequent duplication [16–18].

A large number of NUMTs are considered “dead on arrival” because these sequences are non-functional pseudo-genes, as evidenced through the presence of stop codons, frameshifts, and/or indels in their coding sequences [19, 20]. However, evidence of functional NUMTs has been observed in a few species. For instance, five NUMTs within gene-coding regions have been identified in expressed sequence tags (ESTs) in the honeybee genome [9]. In addition, Ricchetti and coworkers identified 22 out of 28 human-specific NUMTs inserted in known or predicted introns, 1 NUMT in an exon, 1 NUMT in a promoter region, and 4 NUMTs in intergenic regions. These authors also suggested that NUMTs would preferentially integrate into coding or regulatory sequences and cause insertions associated with human diseases and those induced through environmental insults [21].

Cnidaria is one of the earliest branches in the animal tree of life, as evidenced through fossil records dating approximately 600 million years ago [22]. The mtDNA in the phylum Cnidaria represents a “hot spot” of mitochondrial genomic diversity in animals due to variation in both the gene content and mtDNA genome architecture [23]. The cnidarian mitochondrial genomes have two unique characteristics. Among the four traditionally recognised cnidarian classes, species in the Scyphozoan, Cubozoan, and Hydrozoan classes display linear mitochondrial genomes, while those in the Anthozoa class exhibit circular mitochondrial genomes [24, 25]. The H. magnipapillata mitochondrial genome consists of two separate molecules: mitochondrial chromosome 1 (mt-Chr 1) and mitochondrial chromosome 2 (mt-Chr 2) [26]. Despite its fragmentary molecules, the length of the H. magnipapillata mitochondrial genome is 15,880 bp (8,194 bp and 7,686 bp, respectively), showing a typical metazoan mitochondrial genome size [23]. Identical inverted terminal repeats (ITR) occur on both H. magnipapillata mitochondrial chromosomes, similar to those in the Aurelia aurita and H. oligactis linear mitochondrial genome, and both mitochondrial chromosomes possess identical oriented sequences at the 5’ and 3’ ends (5’ and 3’ IOS) adjacent to the ITR [27, 28]. Another unique characteristic of cnidarian mitochondrial genomes is that introns have been observed in several classes [29, 30]. The N. vectensis and A. digitifera mitochondrial genome comprises molecular DNA containing a group I intron in the ND5 gene. However, the genes identified in group I are different [30, 31].

Our knowledge of mitochondrial DNA integrated into animal nuclear genomes is primarily limited to animals with circular mitochondrial genomes without introns. NUMT studies are not available in animals with linear mitochondrial genomes or those with intron-containing mitochondria. Therefore, it remains to be determined whether NUMT transposition differences exist between circular and linear mtDNAs [4, 5]. The completion of the mitochondrial genome sequence and recent availability of the genome draft sequences of Hydra magnipapillata, Nematostella vectensis, and Acropora digitifera have led to the description of NUMT features in the nuclear genomes of these three cnidarian genomes [26, 31–34].

Whole-genome shotgun sequencing was also used to sequence these three cnidarian genomes. Compared with the cost-effective data obtained from Illumina sequencing, the longer sequencing reads from Sanger sequencing or 454 pyro-sequencing were produced to assemble these genomes. Therefore, instead of the popular sequencing-by-hybridisation algorithm based on the k-mer content, overlap-layout-consensus and greedy algorithms were adopted for the assembly of the three cnidarian genomes. The advantage of genome assembly from long sequencing reads is that it provides a more accurate estimation of the NUMT content compared with genome assembly from short sequencing reads [32–34]. A description of NUMT features might provide information concerning the characteristics of mitochondrial pseudo-genes in these three species from the earliest branches in animal evolution, and it might also offer a comparison between NUMTs originating from circular mtDNA molecules and those originating from linear mtDNA molecules. In addition, the present study provides a detailed investigation on the genome-wide identification of NUMTs in three cnidarian species with distinct mitochondrial structures and explores the NUMT landscape in a species with linear mtDNA molecules.

Results

Exceptionally high frequency of NUMT insertions in the H. magnipapillatagenome

We determined the homology relationships between the three cnidarian nuclear genomes and their corresponding mitochondrial genomes using BLASTN searches. A total of 704, 24 and 1 NUMTs were inferred in Hydra magnipapillata, Nematostella vectensis, and Acropora digitifera, respectively (Figure 1, Additional file 1: Table S1, Table S2, Additional file 2: Table S3, and Additional file 3: Table S4). The length distribution of the BLASTN hits was 51 to 7,684 bp (mean: 753 bp) for H. magnipapillata, 76 to 1,836 bp (mean: 768 bp) for N. vectensis, and 86 bp for A. digitifera (Figure 2). No obvious differences were observed in the length of NUMTs between H. magnipapillata and N. vectensis (mean lengths: 753 bp and 768 bp, respectively; Student’s t test: P = 0.8536), while their mean lengths were much longer than those observed in many other metazoans (e.g., Amphimedon queenslandica: 222 bp; Drosophila sechellia: 240 bp; Homo sapiens: 309 bp; Nasonia vitripenni: 565 bp) [5, 35]. The largest proportion of NUMTs was between 100 and 200 bp in H. magnipapillata and between 800 and 900 bp in N. vectensis (Figure 2). The sequence similarity detected using BLASTN between NUMTs and homologous mtDNA sequences was 79.34 to 100%, 92.11 to 100%, and 95.35% in H. magnipapillata, N. vectensis, and A. digitifera, respectively. The mtDNA fraction covered by NUMTs was lower in N. vectensis (72.68%) and A. digitifera (0.47%) compared with H. magnipapillata (100%) (Table 1). NUMTs accounted for approximately 0.0470% (529,934 bp) of the H. magnipapillata genome (Table 1 and Figure 1C), which was much higher than the 0.0052% in N. vectensis (Table 1 and Figure 1B) and 0.00002% (86 bp) in A. digitifera (Table 1 and Figure 1A). The cumulative percentage was used as a measure of NUMT content in the genomes examined in this study, suggesting that the estimation of the NUMT content was not affected by differences in genome assembly fragmentation (Additional file 4: Table S5). Thus, although H. magnipapillata, N. vectensis, and A. digitifera belong to the same phylum, the NUMT compositions of the genomes of these three species are quite different. The number of NUMTs in the H. magnipapillata genome clearly exceeds those of the N. vectensis and A. digitifera genomes.

Figure 1
figure 1

NUMT distribution in three cnidarian mitochondrial genomes. (A) NUMT distribution in the A. digitifera mitochondrial genome. The mtDNA genome contains the genes for 13 energy pathway proteins, two rDNAs and two tRNAs. A single intron occurs in the ND5 gene. Note that the intron occurs within the ND5 gene and contains the only copies of the tRNA-Trp, ATP8, COX1, tRNA-Met, and l-rRNA genes. The black line represents the only NUMT in the A. digitifera genome. (B) NUMT distribution in the N. vectensis mitochondrial genome. A single intron occurs in the ND5 gene. Note that the intron occurs within the ND5 gene and contains the only copies of the ND1 and ND3 genes. Each black line represents a NUMT insertion in the reference genome. (C) NUMT distribution in the H. magnipapillata mitochondrial genome. For clarity, the NUMTs whose length is longer than 500 bp are indicated with black lines.

Full size image
Figure 2
figure 2

NUMT lengths distribution.

Full size image
Table 1 Sizes of mtDNA, nuclear genomes, and NUMTs detected through a BLASTN search
Full size table

NUMTs have been observed in several metazoans, but they are rare or even absent in the oldest branch of animal evolution [5, 35]. Unexpectedly, we detected NUMT proliferation in the H. magnipapillata genome. BLASTN searches within the H. magnipapillata genome assembly using two mtDNA molecules, mt-Chr 1 and mt-Chr 2 as queries yielded 353 and 351 NUMT sequences, respectively (Additional file 1: Table S1 and Table S2). As shown in Figure 1C, these NUMTs covered almost the entire H. magnipapillata mitochondrial genome. The similarities between NUMTs and their corresponding mitochondrial counterparts varied from 79.34 to 100% in mt-Chr 1 and 79.35 to 100% in mt-Chr 2. NUMTs from mt-Chr 1 and mt-Chr 2 showed extensive differences in length, ranging from 51 to 7,684 bp and 51 to 7,421 bp, respectively. The mean lengths of the mt-Chr 2-specific NUMTs were longer than those of the mt-Chr 1-specific NUMTs (mt-Chr 1-specific mean = 694 bp and mt-Chr 2-specific mean = 811 bp; Student’s t test: P = 0.0760; Figure 2). The median number of NUMT length in mt-Chr 1-NUMT, mt-Chr 2-NUMT and N. vectensis-NUMT was 464 bp (its variance: 572,384), 519 bp (its variance: 956,478), and 861 bp (its variance: 142,365), respectively. In the H. magnipapillata genome, 529,934 bp (0.0470%) comprised sequences corresponding to a mitochondrial origin (Table 1), which was approximately 33.37 times longer than the entire original mitochondrial genome. We did not detect any insertion preference for certain mitochondrial genes in the H. magnipapillata genome. Approximately equal amounts of NUMTs in both mitochondrial genomes indicated that there was no obvious bias in the transferred region of both mitochondrial genomes.

Characterisation of NUMTs identified in the three cnidarian genomes

In the three cnidarian genomes, mitochondrial protein-coding genes were transferred into the nuclear genome at different frequencies. Of the 17 mitochondrial genes with NUMT insertions, 17 counterparts were identified in the H. magnipapillata genome, 15 counterparts were identified in the N. vectensis genome, and only 1 counterpart was identified in the A. digitifera genome. Mitochondrial genes with relatively complete structures were observed in both the H. magnipapillata and N. vectensis genomes.

Protein-coding genes, ribosomal genes, and tRNA genes were identified in H. magnipapillata NUMTs, showing varying insertion frequencies. In total, we identified 1,133 mitochondrial gene fragment counterparts in H. magnipapillata, including 257 relatively complete mitochondrial genes located in the 138 NUMTs. Among these, 13 mitochondrial protein-coding genes, 3 tRNAs, and 2 rRNA genes were detected. In contrast, of the 7 NUMTs inferred from the N. vectensis genome, only 9 relatively complete mitochondrial genes were detected, including 2 copies of tRNA-Met genes and 7 mitochondrial protein-coding genes. Unlike those of H. magnipapillata and N. vectensis, the NUMTs of A. digitifera lacked complete mitochondrial genes. The only NUMT identified in the A. digitifera genome was homologous to the ND3 gene fragment. The H. magnipapillata mitochondrial genes were present in multiple copies in the nuclear genome, and the copy numbers of these genes ranged from 19 to 46. Most NUMTs detected in H. magnipapillata were non-functional pseudo-genes, as evidenced by the mutations and short insertions in the coding regions of these genes. However, few protein-coding genes in NUMTs could be translated using mitochondrial and universal codons (e.g., ND4 L, Additional file 5: Figure S1A), and several NUMT tRNAs could be folded into perfect structures (e.g., tRNA-Trp, Additional file 5: Figure S1B).

NUMT duplications in the H. magnipapillatagenome

As in other animals (e.g., human, honeybee), NUMT duplications have occurred in the H. magnipapillata genome [8, 9]. Of all NUMTs in H. magnipapillata (n = 704), 39 mtDNA insertions were duplicated at least once in the nuclear genome. Of these 39 cases, 34 insertions were duplicated twice, 4 insertions were duplicated 3 times, and 1 insertion was duplicated 6 times (Additional file 6: Table S6). The duplicated NUMTs ranged in size from 53 to 2,000 bp, with mtDNA sequence identity ranging from 79.34 to 100%. However, NUMT duplication events were absent in the N. vectensis and A. digitifera genomes.

We identified at least 46 NUMTs that originated from neighbouring mtDNA regions but were located away from each other in the H. magnipapillata nuclear genome. We compared the sequences between pairs of NUMTs manually or using the RepeatMasker program. These sequences were classified into the following two groups: (1) low-complexity regions containing AT-rich or repetitive sequences; and (2) sequences containing transposable elements (Additional file 7: Figure S2).

Insertion of NUMTs in nuclear genes

NUMTs in Cnidaria are located in non-genic or low-gene-density regions, consistent with most bilaterians [19, 36]. We inferred the presence of 23, 1, and 0 intronic NUMTs in H. magnipapillata, N. vectensis, and A. digitifera, respectively. The genomic locations of NUMTs further revealed that most H. magnipapillata NUMTs (94.46%, 665 of 704) were located in non-genic regions (74.74%, 526 of 665) or low-gene-density regions (19.7%, 139 of 665). For N. vectensis, 87.5% (21 of 24) of the NUMTs were inserted in non-genic regions (52.38%, 11 of 24) or low-gene-density regions (47.62%, 10 of 24), while 12.5% (3 out of 24) of the NUMTs were located in regions of high gene density. However, we only detected one NUMT located in a non-genic region in A. digitifera.

We identified a fraction of NUMTs (3.25%, 23 of 704) in H. magnipapillata located in introns, similar to 1 NUMT from scaffold 1 of N. vectensis (Table 2). In H. magnipapillata, genes typically contained a single NUMT. For example, the s-rRNA gene, ND4, ND1, COX1, and COX2 fragments were transposed into introns of protein kinase C, delta, neuroendocrine convertase 1, sox10, and cytochrome b5 reductase 4, respectively. Genes containing more than 1 NUMT were identified in 3 cases (Table 2). The only intron of the mitoferrin-1 gene contained 8 NUMTs, intron 3 of the MAD homolog 4-interacting transcription coactivator 1 gene contained 4 NUMTs, and intron 2 of LOC100199452 contained 2 NUMTs. All NUMTs within a gene were in the same orientation (Additional file 8: Figure S3). Twelve mitochondrial genes had intronic NUMT counterparts. ψCOX1 was identified in 8 different intronic NUMTs, the s-rRNA gene was detected in 5 intronic NUMTs, and all other genes in intronic NUMTs were present at a low frequency. No observed NUMT insertions overlapped with exons.

Table 2 NUMTs in the predicted nuclear genes
Full size table

The expression of NUMTs suggests its functional relevance

Generally, NUMTs lose their functions and are considered “dead on arrival”. These sequences are likely located in intergenic regions or introns and display no transcriptional activity. However, we observed that NUMTs were co-expressed with adjacent protein-coding genes. In this study, all the H. magnipapillata ESTs were retrieved from GenBank to identify NUMT expression. The mitochondrial transcripts were filtered, and authentic ESTs with NUMT insertions were subjected to further analysis. The absence in the H. magnipapillata assembly might reflect the use of an individual with a unique NUMT insertion at a specific genomic locus in the cDNA library construction. This unique NUMT insertion event did not occur in the individuals involved in the Hydra genome-sequencing project. If the NUMT was expressed through neighbouring promoters, we will be able to identify the corresponding insertion in the EST data. Several of the NUMTs identified in H. magnipapillata overlapped with the ESTs. After removing contaminants, we obtained 7 ESTs corresponding to the NUMTs in the partial sequences. BLASTX results show that 4 ESTs (CV151845, CV284212, CX770377, DN243213) partially corresponded to the COX3 gene and that 1 EST partially corresponded to the ND6 gene (CA301932). None of these sequences could be perfectly translated due to the presence of mitochondrial stop codons. NUMTs were inserted into their 5’ untranslated regions (UTRs) in 2 other ESTs (CO538443, DN603666). CO538443 is transcribed from a transcript containing the complete exons of hydramacin-1, and DN603666 originated from the arminin 1b gene.

Discussion

Unique mitochondrial genome structure in cnidarians

In our study, the number of NUMTs in the Hydra magnipapillata genome was clearly higher than in Nematostella vectensis and Acropora digitifera. A total of 704 NUMTs were identified in H. magnipapillata, representing approximately 0.0470% of the genome (Table 1). The H. magnipapillata genome shows a higher cumulative NUMT percentage than the other invertebrate species investigated and a lower NUMT cumulative percentage than that of Apis mellifera (0.0861%) [4, 5, 9].

The structure of cnidarian mtDNA is variable, as both circular and linear mitochondrial molecules have been identified [24, 37, 38]. H. magnipapillata contains two linear mtDNAs, and the N. vectensis and A. digitifera mtDNAs are single circular molecules. Due to increased susceptibility to exonuclease activity, linear chromosomes are less stable than circular chromosomes and are more easily broken when cells experience damage or the degradation of abnormal mitochondria [23, 39]. In addition, in contrast to the two shorter linear DNA chromosomes in H. magnipapillata, the circular chromosomes in the other two species must be linearised before integration into the nuclear genome. The frequency of NUMTs might be associated with the number and stability of mitochondrial molecules [4, 5]. Therefore, when they enter the nucleus, the mtDNAs in H. magnipapillata provide more sources of integrated fragments than those in N. vectensis and A. digitifera. That is, compared with an animal with a circular mitochondrial chromosome, an animal with a linear mitochondrial chromosome might have more sources to generate mitochondrial fragments that integrate into the nucleus.

In the older phyletic lineage of animal including Placozoa, Nematomorpha, and Platyhelminthes, the available data suggest that the ratio of NUMTs to genome size ranges from 0.0002% in Trichoplax adhaerens to 0.0204% in Brugia malayi[5]. Compared with the species in the older phyletic lineages, H. magnipapillata contains a much higher NUMT cumulative percentage. The size of the H. magnipapillata genome is of the same order of magnitude as those in the species with the highest NUMT cumulative percentage (Figure 3A). Therefore, the H. magnipapillata genome has experienced an exceptionally high frequency of NUMT insertions throughout its evolutionary history. In addition, as shown in Figure 3B, our analysis of 49 animal genomes revealed that no correlation exists between the genome size and NUMT cumulative percentage (Pearson correlation test: P = 0.3532). Except in Monodelphis domestica, the NUMT cumulative percentage is relatively low in most species with large genomes, particularly in invertebrates. Therefore, H. magnipapillata is an exception in that it has a large genome and a high NUMT cumulative percentage.

Figure 3
figure 3

The ratio of NUMTs to the total size of the nuclear genome (%) in the species concerned and the NUMT cumulative percentage is not correlated to genome sizes. (A) The ratio of NUMTs to the total size of the nuclear genome (%) in the species concerned. The x-axis contains the name of each species (including species described in the research by Hazkani-Covos [5] and three Cnidaria species), while the y-axis numbers represent the ratio of NUMTs to the total size of the nuclear genome. “*” represents species with large nuclear genomes (lager than 1 Gbp). (B) The NUMT cumulative percentage is not correlated with genome size. A log-log scale graph shows the dependence between NUMT cumulative percentage and genome size. Genome size information was retrieved from the review of Einat Hazkani-Covo et al. [5], “*” indicates H. magnipapillata.

Full size image

The evolutionary trend of plant mitochondrial genomes is opposite to that of metazoan mitochondrial genomes, in which the genome has become smaller and more compact [40]. Despite their conserved function in eukaryotes, plant mitochondrial genomes have a variety of unique features compared with those of metazoan animals. In contrast to metazoan animals, which have experienced dramatic reduction of mitochondrial genome size after a prior mitochondria-to-nucleus rate acceleration, the transfer of mitochondrial genomes to the nuclear genomes is still on-going in plants. Because of different transfer rates among plant mitochondrial genomes, their gene repertoires and mitochondrial genome sizes vary considerably. In addition to the ability to transfer DNA fragments to the nucleus, plants also have a propensity to integrate DNA fragments from various foreign sources, including the chloroplast, plastid and nuclear genomes, via intracellular transfer, and from other species via horizontal transfer [41]. Multiple and on-going gene transfer events are responsible for the considerable variations in genome size among plant mitochondrial genomes. Therefore, a comparison of the NUMT content among plants is not feasible due to the complicated evolutionary characteristics of mitochondrial genomes, which are similar for most protist species. Protists have diverged so far from metazoan animals that their mitochondrial genomes display striking diversity in size and complexity and possess large, spacious, gene-rich sequences. The flux from the mitochondria to the nucleus is a widespread and on-going phenomenon in protist species [42]. The mitochondrial genome linearisation in the Hydra has split its mitochondrial genome into two linear pieces. The linear mitochondrial genomes in protist species consist of only a single type of DNA molecule, suggesting a distinct mechanism for mitochondrial genome linearisation [43]. Taken together, a comparison of NUMT content among metazoan animals, plants, and protists cannot avoid the bias that results from their distinct mitochondrial genome characteristics.

Considering that the species that have diverged in a short evolutionary time-scale possess similar mitochondrial genome characteristics, we also compared the NUMT content in phylogenetically close species to determine the influences of genome structures on NUMT integration in non-metazoan genomes. To determine whether the variations in the protist NUMT content are influenced by mitochondrial genome structures according to the trend observed in metazoan species, we compared the NUMT content in Plasmodium species, which display similar mitochondrial genome characteristics and different mitochondrial genome structures. The genome sequences for six Plasmodium species, including four species with linear mitochondrial genomes and two species with circular mitochondrial genomes, have been determined. Our measurements of NUMT content in these species showed that the NUMT content of three species (3 in 4) with linear mitochondrial genomes is more abundant than that of species with circular mitochondrial genomes, suggesting that this trend also follows that observed in metazoan species (Additional file 9: Table S7). Plasmodium gallinaceum (GC content: 20.6%; Genome size: 16,913,475 bp; NUMT cumulative percentage: 0.0020%) and Plasmodium berghei (GC content: 23.7%, Genome size: 17,952,627 bp; NUMT cumulative percentage: 0%) have similar genome sizes and GC content, but their NUMT cumulative percentage is different, suggesting that the genome size and GC content are not the primary factors to influence the NUMT cumulative percentage.

Another difference between these three species is the genome size. The genome sizes in N. vectensis (357 Mbp) and A. digitifera (420 Mbp) are of the same order of magnitude, whereas the genome of H. magnipapillata (1,240 Mbp) is approximately three times as large as those of the other two species. A comprehensive correlation analysis from multiple genomes shows that the NUMT content is strongly correlated with genome size both in plants and animals [5, 44]. Therefore, it is expected that the cumulative size of NUMTs in H. magnipapillata will be larger than that of the other two species. However, our results indicate that H. magnipapillata is an exception in that it has a large genome size but a high NUMT cumulative percentage. Thus, the NUMT content in H. magnipapillata cannot be fully explained by genome size when compared with the NUMT content of other animals. That is, the increased genome size of H. magnipapillata is not primarily responsible for the large number of NUMTs. Therefore, our results suggest that linearisation of the mitochondrial genome might primarily reflect the high number of NUMTs in the H. magnipapillata genome.

Cnidaria mtDNA displays variations not only in the genome structure but also in genomic components. Introns are frequently identified in cnidarian mitochondrial genomes and have only been reported in cnidarians and sponges among the Metazoa [29, 30, 45]. Both N. vectensis and A. digitifera mitochondrial genomes contain an ND5 intron. The mitochondrial sequences can be integrated in the nuclear genome through two candidate mechanisms: direct DNA transfer and a cDNA intermediate [1, 46, 47]. Based on examples from plants, prevailing views have focused on cDNA as a vehicle when mitochondrial genes enter the nucleus [47]. Evidence supporting this view includes the finding that integrated nuclear copies of genes that originate from the mitochondrial genome often lack introns. In our study, introns were present in both the N. vectensis and A. digitifera mitochondrial genomes. The simultaneous presence of both the intron and exons from ND5 in the NUMT (scaffold_2815 in N. vectensis) indicates that direct DNA transfer is one route of mtDNA transfer into the N. vectensis genome. Therefore, unlike previously studied genomes exhibiting an absence of introns in the mitochondrial genome, our results provide direct evidence that mtDNA integrates into the nuclear genome through direct DNA transfer.

NUMT content exhibits dramatic differences in three cnidarian species

An exceptionally high density of NUMTs was detected in the H. magnipapillata genome. The number of NUMTs in H. magnipapillata clearly exceeded that in N. vectensis and A. digitifera. All NUMTs exhibited high similarity with their mitochondrial counterparts (92.11 to 100% in N. vectensis, 79.34 to 100% in H. magnipapillata, and 95.35% in A. digitifera). No NUMT duplication events were detected in the N. vectensis genome, suggesting that all NUMTs were integrated from the mitochondrial genome. The high sequence similarity between NUMTs and their mitochondrial counterparts indicates that the NUMT integration occurred recently in the N. vectensis genome. In A. digitifera, only one NUMT was identified, and the scarcity of NUMTs suggests that the A. digitifera genome resists NUMT integration. Consequently, although the three species are in the same phylum, their NUMT content and characteristics differ dramatically.

Several early studies suggested that transposable elements or short-dispersed repeats were associated with the on-going integration of mtDNA sequences and their subsequent duplication within the nuclear genome [16–18]. NUMTs are derived from independent insertion events and subsequently undergo further genomic arrangements, resulting in the interruption of older NUMTs through the insertion of repetitive elements. Undoubtedly, duplication increases the NUMT content in H. magnipapillata, as evidenced by the absence of NUMT duplications in A. digitifera and N. vectensis. The transposable element content is much higher in H. magnipapillata than in N. vectensis and A. digitifera[32–34]. Transposable elements have burst three times and are still active in the H. magnipapillata genome [33]. Unlike the H. magnipapillata genome, transposable elements in N. vectensis comprise a small fractionn of the genome and are all relatively young, and NUMTs are more than 90% similar to their mitochondrial counterparts [32]. In addition, several NUMTs in H. magnipapillata were disrupted through transposable elements, while this disruption was not observed in N. vectensis. Thus, a low portion of transposable elements is correlated with a low number of NUMTs. In short, these results imply that the variations in the NUMT content in the three cnidarian species might be correlated with the content and activity of transposable elements in their genomes.

Insertion of NUMTs in nuclear genes and NUMT expression

In the present study, partial sequences in 7 ESTs corresponded to NUMTs. Notably, all NUMT insertions were absent in the original assembly of the H. magnipapillata genome. These differences are polymorphic sites between the genome sequencing and the EST sequencing individuals. The H. magnipapillata genome is estimated to contain approximately 0.7% single-nucleotide polymorphisms between alleles, indicating its substantial heterozygosity [33]. It is expected that several polymorphic NUMT insertions might have occurred in different individuals in relatively recent evolutionary times. Therefore, these NUMT insertions might not be fixed in the H. magnipapillata genome.

NUMTs are located in regions with low gene content [4, 19, 36]. This notion has been recently challenged through evidence that the integration of NUMTs in the fish and sponge genomes occurs primarily in known or predicted coding regions [35, 48]. In addition, 80% of human-specific NUMTs are integrated into known or predicted genes [21]. In our study, 23 NUMTs were inserted into the introns of H. magnipapillata, with 14 NUMTs belonging to known or predicted genes (Table 2). The co-expression of NUMTs with protein-coding genes suggested that NUMT expression might be functionally relevant and important for interactions with neighboring protein-coding genes. Cytochrome b5 reductase 4, toxin A, protein kinase C delta, and proteasome subunit beta type-2 are related to oxidative stress [49–53]; SOX10, proteasome subunit beta type-2, protein kinase C delta, and toxin A are associated with mitochondrial apoptosis or mitochondrial damage [52, 54–57]. Mitoferrin-1 functions as an essential importer of iron for the mitochondrial haeme and iron sulphur cluster in erythroblasts and is necessary for erythropoiesis [58]. Six of the 14 known or predicted genes have been associated with mitochondrial activation, suggesting that the proteins that mediate mitochondrial activation have more opportunities to contact the mitochondria. Genes associated with mitochondrial damage are targets for mitochondrial fragment insertion. These results are consistent with previous results suggesting that genetic or environmental factors that increase the frequency of chromosome breaks play a critical role in providing the impetus for the continued invasion of the human genome through mitochondrial DNA [21].

The primary mechanisms responsible for the appearance of new exons are gene duplication, exon shuffling, lateral gene transfer, and alternative splicing [59]. However, the expressed NUMTs detected in the present study suggest that NUMTs are another possible source of new exons. Four ESTs (CV151845, CV284212, CX770377, DN243213) are non-functional because they contain stop codons when translated with both the universal and mitochondrial genetic codes. Although these four transcripts do not encode proteins and their functions are unknown, they correspond to the same gene, COX3. One possible explanation is that NUMTs are transcribed as non-coding RNAs, as previously suggested [60]. Two other ESTs (CO538443, DN603666) contain NUMTs in their 5’UTRs. CO538443 is transcribed from a complete gene, hydramacin-1, while DN603666 is only transcribed from part of arminin 1b[61, 62]. These results suggest that the insertion of NUMT into a gene may have no effect on its transcription and gradually becomes part of the gene, but it might affect the expression of genes encoding proteins. Unlike the previous report that mitochondrial-derived protein-coding exons belong to various functional classes, our results show that NUMT exons belong to the same functional class [63]. Both hydramacin-1 and arminin 1b are novel antimicrobial proteins that were identified during the investigation of the epithelial defence of H. magnipapillata[61, 62]. Whether all exonised NUMTs in the H. magnipapillata genome distribute in immune-related genes requires further investigation. Although only 7 expressed NUMTs were detected in our study, we provide evidence that NUMTs could be expressed in the transcriptome, which might have the potential to affect co-transcriptional genes.

Conclusions

Our knowledge of mitochondrial DNA integrated into animal nuclear genomes is primarily limited to animals with circular mitochondrial genomes without introns. NUMT studies are not available in animals with linear mitochondrial genomes or those with intron-containing mitochondria. In this study, we provided the first report of NUMT in the phylum Cnidaria, containing mitochondrial genomes with distinct variations in genome structure, including the mtDNA structure (linear or circular) and the presence or absence of introns in protein-coding genes. Our analysis shows that mitochondrial genome linearisation might be responsible for the enrichment of NUMTs in the Hydra magnipapillata genome. In addition, co-transposition of exonic and intronic fragments within NUMTs in Nematostella vectensis provides direct evidence that mitochondrial sequences might be transposed into the nuclear genome through DNA-mediated fragment transfer. Furthermore, expressed NUMTs might be detected in the transcriptome, suggesting that these sequences might have biological relevance and the potential to affect co-transcriptional genes. Taken together, our results provide valuable information about the impact of different mitochondrial genome structures in NUMT evolution.

Methods

The mitochondrial genomes were retrieved from the NCBI database for H. magnipapillata (NC_011220, NC_011221) and N. vectensis (NC_008164) [26, 48]. The genome sequences were downloaded from the NCBI database (h7, 28-JAN-2009; ASM20922v1, 22-AUG-2007; Adig_1.0, 28-JUL-2011).

The mitochondrial genome of A. digitifera is not available yet. Following the strategies used to study H. magnipapillata mtDNAs, we identified the A. digitifera mtDNA sequences from the genome assembly [26]. This strategy for assembling the mtDNAs from genome sequencing data was also approved for the mtDNA identification of invertebrates [26, 64]. Briefly, we performed homology searches against the A. digitifera draft genome using the TBLASTN program with known mitochondrial proteins from closely related species. Subsequently, A. digitifera mitochondrial sequences were downloaded from the NCBI database. These sequences were BLASTN searched against the resulting genome sequences obtained from the TBLASTN searches to recover authentic mitochondrial sequences and to assemble mitochondrial genomes. The annotation of A. digitifera mtDNA was performed according to the annotation procedure of invertebrate mtDNAs, as described in our previous studies [65, 66]. Gene annotation for proteins and rRNAs was performed manually and using the DOGMA program [67]. The tRNA genes were identified using tRNAscan-SE, employing the cove only search mode and the invertebrate mitochondrial genetic code [68].

BLASTN searches of nuclear genome sequences were performed using their corresponding mitochondrial genome sequences, and the maximum expectation value was set to E = 10-4, and hits with lengths less than 50 nt were ignored [69]. The total number, size, and locations of NUMTs were determined from the BLASTN results using custom Perl scripts. BLASTN searches were also conducted locally using individual mitochondrial genes as query sequences in each of the H. magnipapillata, N. vectensis, and A. digitifera genome sequences [26, 32, 34, 69]. The sequences of mitochondrial genes were analysed using BioEdit software and custom Perl scripts. Secondary structures of mt-tRNAs and their corresponding NUMT-tRNAs were predicted using the tRNAscan-SE and Mfold programs with default values [70].

If sequences originated from closely spaced regions in the mtDNA but were located away from each other in the nuclear DNA, the fragment between the two NUMTs was extracted from the genome sequence. Because NUMTs are associated with repetitive elements, we used the RepeatMasker program (http://www.repeatmasker.org/) to identify repetitive elements in these extracted fragments. As reported in Apis mellifera, NUMTs originating from the same location in mtDNA were considered to be duplicated NUMTs in the nuclear genome [9]. The seq_gene.md and seq_gene.q files containing RefSeq transcript information were downloaded from the NCBI database (http://ftp://ftp.ncbi.nih.gov/genomes/Hydra_magnipapillata/ARCHIVE/BUILD.1.1/mapview/seq_gene.md.gz) and were used to analyse the positions of NUMTs in gene models for H. magnipapillata. The presence of predicted genes in a 10-kbp window around each NUMT was considered to determine whether a NUMT was located in a genic or non-genic region (empty, no genes; low gene density, 1–2 genes; high gene density, more than two genes). The location information of NUMTs was retrieved using the JGI N. vectensis v1.0 (Nemve1) Genome Browser (http://genome.jgi-psf.org/Nemve1/Nemve1.home.html) and the A. digitifera Genome (Ver 1.1) (http://marinegenomics.oist.jp/genomes/downloads?project_id=3).

In total, 164,325 H. magnipapillata ESTs were downloaded from GenBank to detect NUMT expression and were searched with mitochondrial genome sequences using BLASTN searches. The BLASTN matches that contained partial sequences not originating from mitochondria were considered to be expressed NUMTs and were further BLASTN searched against the genome sequences to locate their genomic positions.

Abbreviations

NUMT:

Nuclear insertion of mitochondrial sequence

mtDNA:

mitochondrial DNA

bp:

base pair(s)

mt-Chr 1:

mitochondrial Chromosome 1

mt-Chr 2:

mitochondrial Chromosome 2

EST:

Expressed sequence tag

tRNA:

transfer RNA

s-rRNA:

12S ribosomal RNA

l-rRNA:

16S ribosomal RNA

ND1-6:

and 4 L, NADH dehydrogenase subunits 1–6 and 4 L

Met:

Methionine

Trp:

Tryptophan

COX 1–3:

Cytochrome c Oxidase subunits 1–3

UTR:

Untranslated region

Mbp:

Megabase pairs

CYTB:

Cytochrome b

NCBI:

National center for biotechnology information

ψCOX1:

COX1 pseudo-gene.

References

  1. 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. J Mol Evol. 1994, 39: 174-190.

    CAS  PubMed  Google Scholar 

  2. du Buy HG, Riley FL: Hybridization between the nuclear and kinetoplast DNA’s of Leishmania enriettii and between nuclear and mitochondrial DNA’s of mouse liver. Proc Natl Acad Sci USA. 1967, 57: 790-797. 10.1073/pnas.57.3.790.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Bensasson D, Zhang D, Hartl DL, Hewitt GM: Mitochondrial pseudogenes: evolution’s misplaced witnesses. Trends Ecol Evol. 2001, 16: 314-321. 10.1016/S0169-5347(01)02151-6.

    Article  PubMed  Google Scholar 

  4. Richly E, Leister D: NUMTs in sequenced eukaryotic genomes. Mol Biol Evol. 2004, 21: 1081-1084. 10.1093/molbev/msh110.

    Article  CAS  PubMed  Google Scholar 

  5. Hazkani-Covo E, Zeller RM, Martin W: Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes. PLoS Genet. 2010, 6: e1000834-10.1371/journal.pgen.1000834.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Pamilo P, Viljakainen L, Vihavainen A: Exceptionally high density of NUMTs in the honeybee genome. Mol Biol Evol. 2007, 24: 1340-1346. 10.1093/molbev/msm055.

    Article  CAS  PubMed  Google Scholar 

  7. Woischnik M, Moraes CT: Pattern of organization of human mitochondrial pseudogenes in the nuclear genome. Genome Res. 2002, 12: 885-893.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Hazkani-Covo E, Sorek R, Graur D: Evolutionary dynamics of large numts in the human genome: rarity of independent insertions and abundance of post-insertion duplications. J Mol Evol. 2003, 56: 169-174. 10.1007/s00239-002-2390-5.

    Article  CAS  PubMed  Google Scholar 

  9. Behura SK: Analysis of nuclear copies of mitochondrial sequences in honeybee (Apis mellifera) genome. Mol Biol Evol. 2007, 24: 1492-1505. 10.1093/molbev/msm068.

    Article  CAS  PubMed  Google Scholar 

  10. Triant DA, DeWoody JA: Extensive mitochondrial DNA transfer in a rapidly evolving rodent has been mediated by independent insertion events and by duplications. Gene. 2007, 401: 61-70. 10.1016/j.gene.2007.07.003.

    Article  CAS  PubMed  Google Scholar 

  11. Viljakainen L, Oliveira DC, Werren JH, Behura SK: Transfers of mitochondrial DNA to the nuclear genome in the wasp Nasonia vitripennis. Insect Mol Biol. 2010, 19 (Suppl 1): 27-35.

    Article  CAS  PubMed  Google Scholar 

  12. Bensasson D, Feldman MW, Petrov DA: Rates of DNA duplication and mitochondrial DNA insertion in the human genome. J Mol Evol. 2003, 57: 343-354. 10.1007/s00239-003-2485-7.

    Article  CAS  PubMed  Google Scholar 

  13. Cheng Z, Ventura M, She X, Khaitovich P, Graves T, Osoegawa K, Church D, DeJong P, Wilson RK, Paabo S, Rocchi M, Eichler EE: A genome-wide comparison of recent chimpanzee and human segmental duplications. Nature. 2005, 437: 88-93. 10.1038/nature04000.

    Article  CAS  PubMed  Google Scholar 

  14. Bailey JA, Eichler EE: Primate segmental duplications: crucibles of evolution, diversity and disease. Nat Rev Genet. 2006, 7: 552-564.

    Article  CAS  PubMed  Google Scholar 

  15. Krampis K, Tyler BM, Boore JL: Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum. Molecular plant-microbe interactions: MPMI. 2006, 19: 1329-1336. 10.1094/MPMI-19-1329.

    Article  CAS  PubMed  Google Scholar 

  16. Gellissen G, Michaelis G: Gene transfer. Mitochondria to nucleus. Ann NY Acad Sci. 1987, 503: 391-401. 10.1111/j.1749-6632.1987.tb40625.x.

    Article  CAS  PubMed  Google Scholar 

  17. Mishmar D, Ruiz-Pesini E, Brandon M, Wallace DC: Mitochondrial DNA-like sequences in the nucleus (NUMTs): insights into our African origins and the mechanism of foreign DNA integration. Hum Mutat. 2004, 23: 125-133. 10.1002/humu.10304.

    Article  CAS  PubMed  Google Scholar 

  18. Ossorio PN, Sibley LD, Boothroyd JC: Mitochondrial-like DNA sequences flanked by direct and inverted repeats in the nuclear genome of Toxoplasma gondii. J Mol Biol. 1991, 222: 525-536. 10.1016/0022-2836(91)90494-Q.

    Article  CAS  PubMed  Google Scholar 

  19. Pereira SL, Baker AJ: Low number of mitochondrial pseudogenes in the chicken (Gallus gallus) nuclear genome: implications for molecular inference of population history and phylogenetics. BMC Evol Biol. 2004, 4: 17-10.1186/1471-2148-4-17.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Triant DA, DeWoody JA: The occurrence, detection, and avoidance of mitochondrial DNA translocations in mammalian systematics and phylogeography. J Mammal. 2007, 88: 908-920. 10.1644/06-MAMM-A-204R1.1.

    Article  Google Scholar 

  21. Ricchetti M, Tekaia F, Dujon B: Continued colonization of the human genome by mitochondrial DNA. PLoS Biol. 2004, 2: E273-10.1371/journal.pbio.0020273.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Coates MM: Visual ecology and functional morphology of cubozoa (cnidaria). Integr Comp Biol. 2003, 43: 542-548. 10.1093/icb/43.4.542.

    Article  PubMed  Google Scholar 

  23. Kayal E, Bentlage B, Collins AG, Kayal M, Pirro S, Lavrov DV: Evolution of linear mitochondrial genomes in medusozoan cnidarians. Genome Biol Evol. 2012, 4: 1-12. 10.1093/gbe/evr123.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Bridge D, Cunningham CW, Schierwater B, DeSalle R, Buss LW: Class-level relationships in the phylum Cnidaria: evidence from mitochondrial genome structure. Proc Natl Acad Sci USA. 1992, 89: 8750-8753. 10.1073/pnas.89.18.8750.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Bridge D, Cunningham CW, DeSalle R, Buss LW: Class-level relationships in the phylum Cnidaria: molecular and morphological evidence. Mol Biol Evol. 1995, 12: 679-689.

    CAS  PubMed  Google Scholar 

  26. Voigt O, Erpenbeck D, Worheide G: A fragmented metazoan organellar genome: the two mitochondrial chromosomes of Hydra magnipapillata. BMC Genomics. 2008, 9: 350-10.1186/1471-2164-9-350.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Shao Z, Graf S, Chaga OY, Lavrov DV: Mitochondrial genome of the moon jelly Aurelia aurita (Cnidaria, Scyphozoa): a linear DNA molecule encoding a putative DNA-dependent DNA polymerase. Gene. 2006, 381: 92-101.

    Article  CAS  PubMed  Google Scholar 

  28. Kayal E, Lavrov DV: The mitochondrial genome of Hydra oligactis (Cnidaria, Hydrozoa) sheds new light on animal mtDNA evolution and cnidarian phylogeny. Gene. 2008, 410: 177-186. 10.1016/j.gene.2007.12.002.

    Article  CAS  PubMed  Google Scholar 

  29. Beagley CT, Okimoto R, Wolstenholme DR: The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): introns, a paucity of tRNA genes, and a near-standard genetic code. Genetics. 1998, 148: 1091-1108.

    PubMed Central  CAS  PubMed  Google Scholar 

  30. van Oppen MJ, Catmull J, McDonald BJ, Hislop NR, Hagerman PJ, Miller DJ: The mitochondrial genome of Acropora tenuis (Cnidaria; Scleractinia) contains a large group I intron and a candidate control region. J Mol Evol. 2002, 55: 1-13. 10.1007/s00239-001-0075-0.

    Article  CAS  PubMed  Google Scholar 

  31. Medina M, Collins AG, Takaoka TL, Kuehl JV, Boore JL: Naked corals: skeleton loss in Scleractinia. Proc Natl Acad Sci USA. 2006, 103: 9096-9100. 10.1073/pnas.0602444103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, Jurka J, Genikhovich G, Grigoriev IV, Lucas SM, Steele RE, Finnerty JR, Technau U, Martindale MQ, Rokhsar DS: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317: 86-94. 10.1126/science.1139158.

    Article  CAS  PubMed  Google Scholar 

  33. Chapman JA, Kirkness EF, Simakov O, Hampson SE, Mitros T, Weinmaier T, Rattei T, Balasubramanian PG, Borman J, Busam D, Disbennett K, Pfannkoch C, Sumin N, Sutton GG, Viswanathan LD, Walenz B, Goodstein DM, Hellsten U, Kawashima T, Prochnik SE, Putnam NH, Shu S, Blumberg B, Dana CE, Gee L, Kibler DF, Law L, Lindgens D, Martinez DE, Peng J, Wigge PA, Bertulat B, Guder C, Nakamura Y, Ozbek S, Watanabe H, Khalturin K, Hemmrich G, Franke A, Augustin R, Fraune S, Hayakawa E, Hayakawa S, Hirose M, Hwang JS, Ikeo K, Nishimiya-Fujisawa C, Ogura A, Takahashi T, Steinmetz PR, Zhang X, Aufschnaiter R, Eder MK, Gorny AK, Salvenmoser W, Heimberg AM, Wheeler BM, Peterson KJ, Bottger A, Tischler P, Wolf A, Gojobori T, Remington KA, Strausberg RL, Venter JC, Technau U, Hobmayer B, Bosch TC, Holstein TW, Fujisawa T, Bode HR, David CN, Rokhsar DS, Steele RE: The dynamic genome of Hydra. Nature. 2010, 464: 592-596. 10.1038/nature08830.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Shinzato C, Shoguchi E, Kawashima T, Hamada M, Hisata K, Tanaka M, Fujie M, Fujiwara M, Koyanagi R, Ikuta T, Fujiyama A, Miller DJ, Satoh N: Using the Acropora digitifera genome to understand coral responses to environmental change. Nature. 2011, 476: 320-323. 10.1038/nature10249.

    Article  CAS  PubMed  Google Scholar 

  35. Erpenbeck D, Voigt O, Adamski M, Woodcroft BJ, Hooper JN, Worheide G, Degnan BM: NUMTs in the sponge genome reveal conserved transposition mechanisms in metazoans. Mol Biol Evol. 2011, 28: 1-5. 10.1093/molbev/msq217.

    Article  CAS  PubMed  Google Scholar 

  36. Leister D: Origin, evolution and genetic effects of nuclear insertions of organelle DNA. Trends in genetics: TIG. 2005, 21: 655-663. 10.1016/j.tig.2005.09.004.

    Article  CAS  PubMed  Google Scholar 

  37. WARRIOR R, Gall J: The mitochondrial DNA of Hydra attenuata and Hydra littoralis consists of two linear molecules. Arch Sci. 1985, 38: 439-445.

    Google Scholar 

  38. Ender A, Schierwater B: Placozoa are not derived cnidarians: evidence from molecular morphology. Mol Biol Evol. 2003, 20: 130-134. 10.1093/molbev/msg018.

    Article  CAS  PubMed  Google Scholar 

  39. Campbell CL, Thorsness PE: Escape of mitochondrial DNA to the nucleus in yme1 yeast is mediated by vacuolar-dependent turnover of abnormal mitochondrial compartments. J Cell Sci. 1998, 111 (Pt 16): 2455-2464.

    CAS  PubMed  Google Scholar 

  40. Burger G, Gray MW, Lang BF: Mitochondrial genomes: anything goes. Trends in genetics: TIG. 2003, 19: 709-716. 10.1016/j.tig.2003.10.012.

    Article  CAS  PubMed  Google Scholar 

  41. Mower JP, Sloan DB, Alverson AJ: Plant mitochondrial genome diversity: the genomics revolution. Plant Genome Diversity, Volume 1. Edited by: Wendel JF, Greilhuber J, Dolezel J, Leitch IJ. 2012, Vienna: Springer Vienna, 123-144.

    Chapter  Google Scholar 

  42. Gray MW, Lang BF, Cedergren R, Golding GB, Lemieux C, Sankoff D, Turmel M, Brossard N, Delage E, Littlejohn TG, Plante I, Rioux P, Saint-Louis D, Zhu Y, Burger G: Genome structure and gene content in protist mitochondrial DNAs. Nucleic Acids Res. 1998, 26: 865-878. 10.1093/nar/26.4.865.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Burger G, Forget L, Zhu Y, Gray MW, Lang BF: Unique mitochondrial genome architecture in unicellular relatives of animals. Proc Natl Acad Sci USA. 2003, 100: 892-897. 10.1073/pnas.0336115100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Smith DR, Crosby K, Lee RW: Plastids and gene transfer: correlation between nuclear plastid DNA abundance and plastid number supports the limited transfer window hypothesis. Genome Biol Evol. 2011, 3: 365-10.1093/gbe/evr001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Rot C, Goldfarb I, Ilan M, Huchon D: Putative cross-kingdom horizontal gene transfer in sponge (Porifera) mitochondria. BMC Evol Biol. 2006, 6: 71-10.1186/1471-2148-6-71.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Nugent JM, Palmer JD: RNA-mediated transfer of the gene coxII from the mitochondrion to the nucleus during flowering plant evolution. Cell. 1991, 66: 473-481. 10.1016/0092-8674(81)90011-8.

    Article  CAS  PubMed  Google Scholar 

  47. Henze K, Martin W: How do mitochondrial genes get into the nucleus?. Trends in genetics: TIG. 2001, 17: 383-387. 10.1016/S0168-9525(01)02312-5.

    Article  CAS  PubMed  Google Scholar 

  48. Antunes A, Ramos MJ: Discovery of a large number of previously unrecognized mitochondrial pseudogenes in fish genomes. Genomics. 2005, 86: 708-717. 10.1016/j.ygeno.2005.08.002.

    Article  CAS  PubMed  Google Scholar 

  49. He D, Hagen SJ, Pothoulakis C, Chen M, Medina ND, Warny M, LaMont JT: Clostridium difficile toxin a causes early damage to mitochondria in cultured cells. Gastroenterology. 2000, 119: 139-150. 10.1053/gast.2000.8526.

    Article  CAS  PubMed  Google Scholar 

  50. Majumder PK, Mishra NC, Sun X, Bharti A, Kharbanda S, Saxena S, Kufe D: Targeting of protein kinase C delta to mitochondria in the oxidative stress response. Cell growth & differentiation: the molecular biology journal of the American Association for Cancer Research. 2001, 12: 465-470.

    CAS  Google Scholar 

  51. Stolzing A, Grune T: The proteasome and its function in the ageing process. Clin Exp Dermatol. 2001, 26: 566-572. 10.1046/j.1365-2230.2001.00867.x.

    Article  CAS  PubMed  Google Scholar 

  52. Garcia-Fernandez LF, Losada A, Alcaide V, Alvarez AM, Cuadrado A, Gonzalez L, Nakayama K, Nakayama KI, Fernandez-Sousa JM, Munoz A, Sanchez-Puelles JM: Aplidin induces the mitochondrial apoptotic pathway via oxidative stress-mediated JNK and p38 activation and protein kinase C delta. Oncogene. 2002, 21: 7533-7544. 10.1038/sj.onc.1205972.

    Article  CAS  PubMed  Google Scholar 

  53. Lee CG, Kim YW, Kim EH, Meng Z, Huang W, Hwang SJ, Kim SG: Farnesoid X receptor protects hepatocytes from injury by repressing miR-199a-3p, which increases levels of LKB1. Gastroenterology. 2012, 142: 1206-1217. 10.1053/j.gastro.2012.01.007. e1207

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Schulze-Osthoff K, Bakker AC, Vanhaesebroeck B, Beyaert R, Jacob WA, Fiers W: Cytotoxic activity of tumor necrosis factor is mediated by early damage of mitochondrial functions. Evidence for the involvement of mitochondrial radical generation. J Biol Chem. 1992, 267: 5317-5323.

    CAS  PubMed  Google Scholar 

  55. Cross T, Griffiths G, Deacon E, Sallis R, Gough M, Watters D, Lord JM: PKC-delta is an apoptotic lamin kinase. Oncogene. 2000, 19: 2331-2337. 10.1038/sj.onc.1203555.

    Article  CAS  PubMed  Google Scholar 

  56. Tanaka Y, Engelender S, Igarashi S, Rao RK, Wanner T, Tanzi RE, Sawa A, LD V, Dawson TM, Ross CA: Inducible expression of mutant alpha-synuclein decreases proteasome activity and increases sensitivity to mitochondria-dependent apoptosis. Hum Mol Genet. 2001, 10: 919-926. 10.1093/hmg/10.9.919.

    Article  CAS  PubMed  Google Scholar 

  57. Honore SM, Aybar MJ, Mayor R: Sox10 is required for the early development of the prospective neural crest in Xenopus embryos. Dev Biol. 2003, 260: 79-96. 10.1016/S0012-1606(03)00247-1.

    Article  CAS  PubMed  Google Scholar 

  58. Shaw GC, Cope JJ, Li L, Corson K, Hersey C, Ackermann GE, Gwynn B, Lambert AJ, Wingert RA, Traver D, Trede NS, Barut BA, Zhou Y, Minet E, Donovan A, Brownlie A, Balzan R, Weiss MJ, Peters LL, Kaplan J, Zon LI, Paw BH: Mitoferrin is essential for erythroid iron assimilation. Nature. 2006, 440: 96-100. 10.1038/nature04512.

    Article  CAS  PubMed  Google Scholar 

  59. Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003, 4: 865-875.

    Article  CAS  PubMed  Google Scholar 

  60. Tsuji J: Immigrants to the Nucleus: Analysis of Mitochondrially Derived Nuclear Genomic Regions (NUMT). Master thesis. 2010, Tokyo University: Computational Biology Department

    Google Scholar 

  61. Augustin R, Anton-Erxleben F, Jungnickel S, Hemmrich G, Spudy B, Podschun R, Bosch TC: Activity of the novel peptide arminin against multiresistant human pathogens shows the considerable potential of phylogenetically ancient organisms as drug sources. Antimicrob Agents Chemother. 2009, 53: 5245-5250. 10.1128/AAC.00826-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Jung S, Dingley AJ, Augustin R, Anton-Erxleben F, Stanisak M, Gelhaus C, Gutsmann T, Hammer MU, Podschun R, Bonvin AM, Leippe M, Bosch TC, Grotzinger J: Hydramacin-1, structure and antibacterial activity of a protein from the basal metazoan Hydra. J Biol Chem. 2009, 284: 1896-1905.

    Article  CAS  PubMed  Google Scholar 

  63. Noutsos C, Kleine T, Armbruster U, DalCorso G, Leister D: Nuclear insertions of organellar DNA can create novel patches of functional exon sequences. Trends in genetics : TIG. 2007, 23: 597-601. 10.1016/j.tig.2007.08.016.

    Article  CAS  PubMed  Google Scholar 

  64. Feldmeyer B, Hoffmeier K, Pfenninger M: The complete mitochondrial genome of Radix balthica (Pulmonata, Basommatophora), obtained by low coverage shot gun next generation sequencing. Mol Phylogenet Evol. 2010, 57: 1329-1333. 10.1016/j.ympev.2010.09.012.

    Article  CAS  PubMed  Google Scholar 

  65. Ma C, Yang P, Jiang F, Chapuis MP, Shali Y, Sword GA, Kang L: Mitochondrial genomes reveal the global phylogeography and dispersal routes of the migratory locust. Mol Ecol. 2012, 21: 4344-4358. 10.1111/j.1365-294X.2012.05684.x.

    Article  PubMed  Google Scholar 

  66. Ren J, Liu X, Jiang F, Guo X, Liu B: Unusual conservation of mitochondrial gene order in Crassostrea oysters: evidence for recent speciation in Asia. BMC Evol Biol. 2010, 10: 394-10.1186/1471-2148-10-394.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  67. Wyman SK, Jansen RK, Boore JL: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20: 3252-3255. 10.1093/bioinformatics/bth352.

    Article  CAS  PubMed  Google Scholar 

  68. Schattner P, Brooks AN, Lowe TM: The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005, 33: W686-W689. 10.1093/nar/gki366.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.

    Article  CAS  PubMed  Google Scholar 

  70. Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This study was financially supported through funding from the Natural Science Foundation Key Project of Yunnan Province, China (No. 2007C0003Z), the National Natural Science Foundation of China (No.30660024), the Applied and Basic Research Foundation of Yunnan Province, China (No. 2006C0034M), and the Science Foundation of the Education Commission of Yunnan Province, China (No. 5Y0196B).

Author information

Authors and Affiliations

  1. Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming, Yunnan, 650201, China

    Shen Song & Yongwang Miao

  2. Institute of Zoology, Chinese Academy of Sciences, Beijing, 100101, China

    Feng Jiang & Wei Guo

  3. Beijing Institutes of Life Science, Chinese Academy of Sciences, Beijing, 100101, China

    Feng Jiang

  4. Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China

    Jianbo Yuan

Authors
  1. Shen Song
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Feng Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jianbo Yuan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Wei Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yongwang Miao
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Feng Jiang or Yongwang Miao.

Additional information

Competing interests

The authors declare that there are no competing interests.

Authors’ contributions

YM and FJ conceived and designed the study. SS, FJ, and YM performed the data analyses and drafted the manuscript. JY and WG assisted in bioinformatic analyses and interpreted the results. All authors have read and approved the manuscript.

Electronic supplementary material

12864_2013_5175_MOESM1_ESM.xls

Additional file 1: Table S1 and Table S2: NUMTs detected in the H. magnipapillata genome, and parameters of alignments returned on BLASTN searches.Start and end indicates the positions of the alignments in the H. magnipapillata mitochondrial (mt-Chr 1 and mt-Chr 2) and nuclear genomes. Orientation corresponds to whether integration in the nuclear genome is 5’ > 3’ (+) or 3’ > 5’ (−). The E-value and identity with H. magnipapillata mtDNA are indicated as the expected value and % of identity, respectively, for each returned alignment. (XLS 158 KB)

12864_2013_5175_MOESM2_ESM.xls

Additional file 2: Table S3: NUMTs detected in the N. vectensis genome, and parameters of alignments returned on BLASTN searches.Start and end indicate the positions of the alignments in the N. vectensis mitochondrial (mtDNA) and nuclear genomes. Orientation corresponds to whether integration in the nuclear genome is 5’ > 3’ (+) or 3’ > 5’ (−). The E-value and identity with N. vectensis mtDNA are indicated as the expected value and % of identity, respectively, for each returned alignment. (XLS 27 KB)

12864_2013_5175_MOESM3_ESM.xls

Additional file 3: Table S4: NUMTs detected in the A. digitifera genome, and parameters of alignments returned on BLASTN searches. Start and end indicate the positions of the alignments in the A. digitifera mitochondrial (mtDNA) and nuclear genomes. Orientation corresponds to whether integration in the nuclear genome is 5’ > 3’ (+) or 3’ > 5’ (−). The E-value and identity with A. digitifera mtDNA are indicated as the expected value and % of identity, respectively, for each returned alignment. (XLS 24 KB)

Additional file 4: Table S5: Sequencing and assembly statistics for three cnidarian genomes. (XLS 24 KB)

12864_2013_5175_MOESM5_ESM.pdf

Additional file 5: Figure S1: Protein-coding genes in NUMTs could be translated perfectly by mitochondrial and universal codons, and several NUMT tRNAs which could be folded to perfect structures. (A) Protein-coding genes (ND4 L) in NUMTs could be perfectly translated using mitochondrial and universal codons. (B) NUMT tRNAs (e.g., tRNA-Trp) could be folded into perfect structures through simulations using the Mfold web server. (PDF 339 KB)

Additional file 6: Table S6: A list of NUMT duplications identified in the H. magnipapillata genome. (XLS 32 KB)

12864_2013_5175_MOESM7_ESM.pdf

Additional file 7: Figure S2: Examples of H. magnipapillata NUMTs interrupted by repetitive elements and transposable elements. (A) Blue boxes depict NUMTs, and transparent boxes depict repetitive elements. (B) Blue boxes depict NUMTs, and transparent boxes depict transposable elements. The numbers (in bold) on top of the blue boxes indicate the beginning and end positions of NUMTs on the scaffold. (PDF 305 KB)

12864_2013_5175_MOESM8_ESM.pdf

Additional file 8: Figure S3: Schematic view of predicted genes harbouring NUMTs within introns. (A) Schematic view of selected H. magnipapillata predicted genes harbouring NUMTs within introns. White boxes depict exons, black boxes depict NUMTs located in the introns, and grey boxes depict UTRs. The white triangle indicates the orientation with respect to the gene. (B) Schematic view of the predicted N. vectensis genes harbouring NUMTs within introns. (PDF 375 KB)

12864_2013_5175_MOESM9_ESM.xls

Additional file 9: Table S7: NUMT content in protist nuclear genomes. Note: MT size, mitochondrial genome size. NUMT P, NUMT cumulative. (XLS 31 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Song, S., Jiang, F., Yuan, J. et al. Exceptionally high cumulative percentage of NUMTs originating from linear mitochondrial DNA molecules in the Hydra magnipapillata genome. BMC Genomics 14, 447 (2013). https://doi.org/10.1186/1471-2164-14-447

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-14-447

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us