Characterization of mitochondrial proteomes of nonbilaterian animals.

Mitochondria require ~1,500 proteins for their maintenance and proper functionality, which constitute the mitochondrial proteome (mt-proteome). Although a few of these proteins, mostly subunits of the electron transport chain complexes, are encoded in mitochondrial DNA (mtDNA), the vast majority are encoded in the nuclear genome and imported to the organelle. Previous studies have shown a continuous and complex evolution of mt-proteome among eukaryotes. However, there was less attention paid to mt-proteome evolution within Metazoa, presumably because animal mtDNA and, by extension, animal mitochondria are often considered to be uniform. In this analysis, two bioinformatic approaches (Orthologue-detection and Mitochondrial Targeting Sequence prediction) were used to identify mt-proteins in 23 species from four nonbilaterian phyla: Cnidaria, Ctenophora, Placozoa, and Porifera, as well as two choanoflagellates, the closest animal relatives. Our results revealed a large variation in mt-proteome in nonbilaterian animals in size and composition. Myxozoans, highly reduced cnidarian parasites, possessed the smallest inferred mitochondrial proteomes, while calcareous sponges possessed the largest. About 513 mitochondrial orthologous groups were present in all nonbilaterian phyla and human. Interestingly, 42 human mitochondrial proteins were not identified in any nonbilaterian species studied and represent putative innovations along the bilaterian branch. Several of these proteins were involved in apoptosis and innate immunity, two processes known to evolve within Metazoa. Conversely, several proteins identified as mitochondrial in nonbilaterian phyla and animal outgroups were absent in human, representing cases of possible loss. Finally, a few human cytosolic proteins, such as histones and cytosolic ribosomal proteins, were predicted to be targeted to mitochondria in nonbilaterian animals. Overall, our analysis provides the first step in characterization of mt-proteomes in nonbilaterian animals and understanding evolution of animal mt-proteome. © 2018 IUBMB Life, 70(12):1289-1301, 2018.

Ancestral state reconstruction infers phytopathogenic origins of sooty blotch and flyspeck fungi on apple.

Members of the sooty blotch and flyspeck (SBFS) complex are epiphytic fungi in the Ascomycota that cause economically damaging blemishes of apples worldwide. SBFS fungi are polyphyletic, but approx. 96% of SBFS species are in the Capnodiales. Evolutionary origins of SBFS fungi remain unclear, so we attempted to infer their origins by means of ancestral state reconstruction on a phylogenetic tree built utilizing genes for the nuc 28S rDNA (approx. 830 bp from near the 59 end) and the second largest subunit of RNA polymerase II (RPB2). The analyzed taxa included the well-known genera of SBFS as well as non-SBFS fungi from seven families within the Capnodiales. The non-SBFS taxa were selected based on their distinct ecological niches, including plant-parasitic and saprophytic species. The phylogenetic analyses revealed that most SBFS species in the Capnodiales are closely related to plant-parasitic fungi. Ancestral state reconstruction provided strong evidence that plant-parasitic fungi were the ancestors of the major SBFS lineages. Knowledge gained from this study may help to better understand the ecology and evolution of epiphytic fungi.

Mitochondrial DNA of Clathrina clathrus (Calcarea, Calcinea): six linear chromosomes, fragmented rRNAs, tRNA editing, and a novel genetic code.

Sponges (phylum Porifera) are a large and ancient group of morphologically simple but ecologically important aquatic animals. Although their body plan and lifestyle are relatively uniform, sponges show extensive molecular and genetic diversity. In particular, mitochondrial genomes from three of the four previously studied classes of Porifera (Demospongiae, Hexactinellida, and Homoscleromorpha) have distinct gene contents, genome organizations, and evolutionary rates. Here, we report the mitochondrial genome of Clathrina clathrus (Calcinea, Clathrinidae), a representative of the fourth poriferan class, the Calcarea, which proves to be the most unusual. Clathrina clathrus mitochondrial DNA (mtDNA) consists of six linear chromosomes 7.6-9.4 kb in size and encodes at least 37 genes: 13 protein codings, 2 ribosomal RNAs (rRNAs), and 24 transfer RNAs (tRNAs). Protein genes include atp9, which has now been found in all major sponge lineages, but no atp8. Our analyses further reveal the presence of a novel genetic code that involves unique reassignments of the UAG codons from termination to tyrosine and of the CGN codons from arginine to glycine. Clathrina clathrus mitochondrial rRNAs are encoded in three (srRNA) and ≥6 (lrRNA) fragments distributed out of order and on several chromosomes. The encoded tRNAs contain multiple mismatches in the aminoacyl acceptor stems that are repaired posttranscriptionally by 3'-end RNA editing. Although our analysis does not resolve the phylogenetic position of calcareous sponges, likely due to their high rates of mitochondrial sequence evolution, it confirms mtDNA as a promising marker for population studies in this group. The combination of unusual mitochondrial features in C. clathrus redefines the extremes of mtDNA evolution in animals and further argues against the idea of a "typical animal mtDNA."

Characterization of the Mitochondrial Proteome in the Ctenophore Mnemiopsis leidyi Using MitoPredictor.

Mitochondrial proteomes have been experimentally characterized for only a handful of animal species. However, the increasing availability of genomic and transcriptomic data allows one to infer mitochondrial proteins using computational tools. MitoPredictor is a novel random forest classifier, which utilizes orthology search, mitochondrial targeting signal (MTS) identification, and protein domain content to infer mitochondrial proteins in animals. MitoPredictor's output also includes an easy-to-use R Shiny applet for the visualization and analysis of the results. In this article, we provide a guide for predicting and analyzing the mitochondrial proteome of the ctenophore Mnemiopsis leidyi using MitoPredictor.

Mitochondrial DNA of the Demosponge Acanthella acuta: Linear Architecture and Other Unique Features.

While Acanthella acuta Schmidt 1862, a common demosponge found in the Mediterranean Sea and Atlantic Ocean, is morphologically similar to other sponges, its mitochondrial DNA (mtDNA) is unique within the class. In contrast to all other studied demosponges, the mtDNA of A. acuta is inferred to be linear and displays several unusual features such as inverted terminal repeats, group II introns in three mitochondrial genes, and two unique open reading frames (ORFs): one of which (ORF1535) combines a DNA polymerase domain with a DNA-directed RNA polymerase domain, while the second bears no discernible similarity to any reported sequences. The group II intron within the cox2 gene is the first such intron reported in an animal. Our phylogenetic analyses indicate that the cox1 intron is related to similar introns found in other demosponges, while the cox2 intron is likely not of animal origin. The two domains found within ORF1535 do not share a common origin and, along with the cox2 intron, were likely acquired by horizontal gene transfer. The findings of this paper open new avenues of exploration in the understanding of mtDNA linearization within Metazoa.

Expansion of the MutS Gene Family in Plants.

The MutS gene family is distributed across the tree of life and is involved in recombination, DNA repair, and protein translation. Multiple evolutionary processes have expanded the set of MutS genes in plants relative to other eukaryotes. Here, we investigate the origins and functions of these plant-specific genes. Land plants, green algae, red algae, and glaucophytes share cyanobacterial-like MutS1 and MutS2 genes that presumably were gained via plastid endosymbiotic gene transfer. MutS1 was subsequently lost in some taxa, including seed plants, whereas MutS2 was duplicated in Viridiplantae (i.e., land plants and green algae) with widespread retention of both resulting paralogs. Viridiplantae also have two anciently duplicated copies of the eukaryotic MSH6 gene (i.e., MSH6 and MSH7) and acquired MSH1 via horizontal gene transfer - potentially from a nucleocytovirus. Despite sharing the same name, "plant MSH1" is not directly related to the gene known as MSH1 in some fungi and animals, which may be an ancestral eukaryotic gene acquired via mitochondrial endosymbiosis and subsequently lost in most eukaryotic lineages. There has been substantial progress in understanding the functions of MSH1 and MSH6/MSH7 in plants, but the roles of the cyanobacterial-like MutS1 and MutS2 genes remain uncharacterized. Known functions of bacterial homologs and predicted protein structures, including fusions to diverse nuclease domains, provide hypotheses about potential molecular mechanisms. Because most plant-specific MutS proteins are targeted to the mitochondria and/or plastids, the expansion of this family appears to have played a large role in shaping plant organelle genetics.

Cnidarian phylogenetic relationships as revealed by mitogenomics.

BACKGROUND: Cnidaria (corals, sea anemones, hydroids, jellyfish) is a phylum of relatively simple aquatic animals characterized by the presence of the cnidocyst: a cell containing a giant capsular organelle with an eversible tubule (cnida). Species within Cnidaria have life cycles that involve one or both of the two distinct body forms, a typically benthic polyp, which may or may not be colonial, and a typically pelagic mostly solitary medusa. The currently accepted taxonomic scheme subdivides Cnidaria into two main assemblages: Anthozoa (Hexacorallia + Octocorallia) - cnidarians with a reproductive polyp and the absence of a medusa stage - and Medusozoa (Cubozoa, Hydrozoa, Scyphozoa, Staurozoa) - cnidarians that usually possess a reproductive medusa stage. Hypothesized relationships among these taxa greatly impact interpretations of cnidarian character evolution. RESULTS: We expanded the sampling of cnidarian mitochondrial genomes, particularly from Medusozoa, to reevaluate phylogenetic relationships within Cnidaria. Our phylogenetic analyses based on a mitochogenomic dataset support many prior hypotheses, including monophyly of Hexacorallia, Octocorallia, Medusozoa, Cubozoa, Staurozoa, Hydrozoa, Carybdeida, Chirodropida, and Hydroidolina, but reject the monophyly of Anthozoa, indicating that the Octocorallia + Medusozoa relationship is not the result of sampling bias, as proposed earlier. Further, our analyses contradict Scyphozoa [Discomedusae + Coronatae], Acraspeda [Cubozoa + Scyphozoa], as well as the hypothesis that Staurozoa is the sister group to all the other medusozoans. CONCLUSIONS: Cnidarian mitochondrial genomic data contain phylogenetic signal informative for understanding the evolutionary history of this phylum. Mitogenome-based phylogenies, which reject the monophyly of Anthozoa, provide further evidence for the polyp-first hypothesis. By rejecting the traditional Acraspeda and Scyphozoa hypotheses, these analyses suggest that the shared morphological characters in these groups are plesiomorphies, originated in the branch leading to Medusozoa. The expansion of mitogenomic data along with improvements in phylogenetic inference methods and use of additional nuclear markers will further enhance our understanding of the phylogenetic relationships and character evolution within Cnidaria.

Phylomitogenomics bolsters the high-level classification of Demospongiae (phylum Porifera).

Class Demospongiae is the largest in the phylum Porifera (Sponges) and encompasses nearly 8,000 accepted species in three subclasses: Keratosa, Verongimorpha, and Heteroscleromorpha. Subclass Heteroscleromorpha contains ∼90% of demosponge species and is subdivided into 17 orders. The higher level classification of demosponges underwent major revision as the result of nearly three decades of molecular studies. However, because most of the previous molecular work only utilized partial data from a small number of nuclear and mitochondrial (mt) genes, this classification scheme needs to be tested by larger datasets. Here we compiled a mt dataset for 136 demosponge species-including 64 complete or nearly complete and six partial mt-genome sequences determined or assembled for this study-and used it to test phylogenetic relationships among Demospongiae in general and Heteroscleromorpha in particular. We also investigated the phylogenetic position of Myceliospongia araneosa, a highly unusual demosponge without spicules and spongin fibers, currently classified as Demospongiae incertae sedis, for which molecular data were not available. Our results support the previously inferred sister-group relationship between Heteroscleromorpha and Keratosa + Verongimorpha and suggest five main clades within Heteroscleromorpha: Clade C0 composed of order Haplosclerida; Clade C1 composed of Scopalinida, Sphaerocladina, and Spongillida; Clade C2 composed of Axinellida, Biemnida, Bubarida; Clade C3 composed of Tetractinellida; and Clade C4 composed of Agelasida, Clionaida, Desmacellida, Merliida, Suberitida, Poecilosclerida, Polymastiida, and Tethyida. The inferred relationships among these clades were (C0(C1(C2(C3+C4)))). Analysis of molecular data from M. araneosa placed it in the C3 clade as a sister taxon to the highly skeletonized tetractinellids Microscleroderma sp. and Leiodermatium sp. Molecular clock analysis dated divergences among the major clades in Heteroscleromorpha from the Cambrian to the Early Silurian, the origins of most heteroscleromorph orders in the middle Paleozoic, and the most basal splits within these orders around the Paleozoic to Mesozoic transition. Overall, the results of this study are mostly congruent with the accepted classification of Heteroscleromorpha, but add temporal perspective and new resolution to phylogenetic relationships within this subclass.

Extensive and evolutionarily persistent mitochondrial tRNA editing in Velvet Worms (phylum Onychophora).

Mitochondrial genomes of onychophorans (velvet worms) present an interesting problem: Some previous studies reported them lacking several transfer RNA (tRNA) genes, whereas others found that all their tRNA genes were present but severely reduced. To resolve this discrepancy, we determined complete mitochondrial DNA (mtDNA) sequences of the onychophorans Oroperipatus sp. and Peripatoides sympatrica as well as cDNA sequences from 14 and 10 of their tRNAs, respectively. We show that tRNA genes in these genomes are indeed highly reduced and encode truncated molecules, which are restored to more conventional structures by extensive tRNA editing. During this editing process, up to 34 nucleotides are added to the tRNA sequences encoded in Oroperipatus sp. mtDNA, rebuilding the aminoacyl acceptor stem, the TΨC arm, and in some extreme cases, the variable arm and even a part of the anticodon stem. The editing is less extreme in P. sympatrica in which at least a part of the TΨC arm is always encoded in mtDNA. When the entire TΨC arm is added de novo in Oroperipatus sp., the sequence of this arm is either identical or similar among different tRNA species, yet the sequences show substantial variation for each tRNA. These observations suggest that the arm is rebuilt, at least in part, by a template-independent mechanism and argue against the alternative possibility that tRNA genes or their parts are imported from the nucleus. By contrast, the 3' end of the aminoacyl acceptor stem is likely restored by a template-dependent mechanism. The extreme tRNA editing reported here has been preserved for >140 My as it was found in both extant families of onychophorans. Furthermore, a similar type of tRNA editing may be present in several other groups of arthropods, which show a high degree of tRNA gene reduction in their mtDNA.

Comparison of Haemophilus parasuis reference strains and field isolates by using random amplified polymorphic DNA and protein profiles.

BACKGROUND: Haemophilus parasuis is the causative agent of Glässer's disease and is a pathogen of swine in high-health status herds. Reports on serotyping of field strains from outbreaks describe that approximately 30% of them are nontypeable and therefore cannot be traced. Molecular typing methods have been used as alternatives to serotyping. This study was done to compare random amplified polymorphic DNA (RAPD) profiles and whole cell protein (WCP) lysate profiles as methods for distinguishing H. parasuis reference strains and field isolates. RESULTS: The DNA and WCP lysate profiles of 15 reference strains and 31 field isolates of H. parasuis were analyzed using the Dice and neighbor joining algorithms. The results revealed unique and reproducible DNA and protein profiles among the reference strains and field isolates studied. Simpson's index of diversity showed significant discrimination between isolates when three 10 mer primers were combined for the RAPD method and also when both the RAPD and WCP lysate typing methods were combined. CONCLUSIONS: The RAPD profiles seen among the reference strains and field isolates did not appear to change over time which may reflect a lack of DNA mutations in the genes of the samples. The recent field isolates had different WCP lysate profiles than the reference strains, possibly because the number of passages of the type strains may affect their protein expression.

Large dataset of octocoral mitochondrial genomes provides new insights into mt-mutS evolution and function.

All studied octocoral mitochondrial genomes (mt-genomes) contain a homologue of the Escherichia coli mutS gene, a member of a gene family encoding proteins involved in DNA mismatch repair, other types of DNA repair, meiotic recombination, and other functions. Although mutS homologues are found in all domains of life, as well as viruses, octocoral mt-mutS is the only such gene found in an organellar genome. While the function of mtMutS is not known, its domain architecture, conserved sequence, and presence of several characteristic residues suggest its involvement in mitochondrial DNA repair. This inference is supported by exceptionally low rates of mt-sequence evolution observed in octocorals. Previous studies of mt-mutS have been limited by the small number of octocoral mt-genomes available. We utilized sequence-capture data from the recent Quattrini et al. 2020 study [Nature Ecology & Evolution 4:1531-1538] to assemble complete mt-genomes for 94 species of octocorals. Combined with sequences publicly available in GenBank, this resulted in a dataset of 184 complete mt-genomes, which we used to re-analyze the conservation and evolution of mt-mutS. In our analysis, we discovered the first case of mt-mutS loss among octocorals in one of the two Pseudoanthomastus spp. assembled from Quattrini et al. data. This species displayed accelerated rate and changed patterns of nucleotide substitutions in mt-genome, which we argue provide additional evidence for the role of mtMutS in DNA repair. In addition, we found accelerated mt-sequence evolution in the presence of mt-mutS in several octocoral lineages. This accelerated evolution did not appear to be the result of relaxed selection pressure and did not entail changes in patterns of nucleotide substitutions. Overall, our results support previously reported patterns of conservation in mt-mutS and suggest that mtMutS is involved in DNA repair in octocoral mitochondria. They also indicate that the presence of mt-mutS contributes to, but does not fully explain, the low rates of sequence evolution in octocorals.

Rapid proliferation of repetitive palindromic elements in mtDNA of the endemic Baikalian sponge Lubomirskia baicalensis.

Animal mitochondrial DNA (mtDNA) is a remarkably compact molecule largely because of the scarcity of noncoding "selfish" DNA. Recently, however, we found that mitochondrial genomes of several phylogenetically diverse species of demosponges contain small repetitive palindromic sequences, interspersed within intergenic regions and fused in protein and ribosomal RNA genes. Here, I report and analyze the proliferation of such elements in the mitochondrial genome of the endemic sponge of Lake Baikal Lubomirskia baicalensis. Because Baikal sponges are closely related to the circumglobally distributed freshwater sponge Ephydatia muelleri with which they shared a common ancestor approximately 3-10 Ma, both the rate of single nucleotide substitutions and the rate of palindromic repeat insertions can be calculated in this system. I found the rate of nucleotide substitutions in mtDNA of freshwater sponges to be extremely low (0.5-1.6 x 10(-9) per site per year), more similar to that in plants than bilaterian animals. By contrast, the per/nucleotide rate of insertions of repetitive elements is at least four times higher. This rapid rate of proliferation combined with the broad phylogenetic distribution of hairpin elements can make them a defining force in the evolution of mitochondrial genomes of demosponges.

Parallel loss of nuclear-encoded mitochondrial aminoacyl-tRNA synthetases and mtDNA-encoded tRNAs in Cnidaria.

Unlike most animal mitochondrial (mt) genomes, which encode a set of 22 transfer RNAs (tRNAs) sufficient for mt protein synthesis, those of cnidarians have only retained one or two tRNA genes. Whether the missing cnidarian mt-tRNA genes relocated outside the main mt chromosome or were lost remains unclear. It is also unknown what impact the loss of tRNA genes had on other components of the mt translational machinery. Here, we explored the nuclear genome of the cnidarian Nematostella vectensis for the presence of mt-tRNA genes and their corresponding mt aminoacyl-tRNA synthetases (mt-aaRS). We detected no candidates for mt-tRNA genes and only two mt-aaRS orthologs. At the same time, we found that all but one cytosolic aaRS appear to be targeted to mitochondria. These results indicate that the loss of mt-tRNAs in Cnidaria is genuine and occurred in parallel with the loss of nuclear-encoded mt-aaRS. Our phylogenetic analyses of individual aaRS revealed that although the nearly total loss of mt-aaRS is rare, aaRS gene deletion and replacement have occurred throughout the evolution of Metazoa.

RNA interference in marine and freshwater sponges: actin knockdown in Tethya wilhelma and Ephydatia muelleri by ingested dsRNA expressing bacteria.

BACKGROUND: The marine sponge Tethya wilhelma and the freshwater sponge Ephydatia muelleri are emerging model organisms to study evolution, gene regulation, development, and physiology in non-bilaterian animal systems. Thus far, functional methods (i.e., loss or gain of function) for these organisms have not been available. RESULTS: We show that soaking developing freshwater sponges in double-stranded RNA and/or feeding marine and freshwater sponges bacteria expressing double-stranded RNA can lead to RNA interference and reduction of targeted transcript levels. These methods, first utilized in C. elegans, have been adapted for the development and feeding style of easily cultured marine and freshwater poriferans. We demonstrate phenotypic changes result from 'knocking down' expression of the actin gene. CONCLUSION: This technique provides an easy, efficient loss-of-function manipulation for developmental and gene regulatory studies in these important non-bilaterian animals.

Multiple Losses of MSH1, Gain of mtMutS, and Other Changes in the MutS Family of DNA Repair Proteins in Animals.

MutS is a key component of the mismatch repair (MMR) pathway. Members of the MutS protein family are present in prokaryotes, eukaryotes, and viruses. Six MutS homologs (MSH1-6) have been identified in yeast, of which three function in nuclear MMR, while MSH1 functions in mitochondrial DNA repair. MSH proteins are believed to be well conserved in animals, except for MSH1-which is thought to be lost. Two intriguing exceptions to this general picture have been found, both in the class Anthozoa within the phylum Cnidaria. First, an ortholog of the yeast-MSH1 was reported in one hexacoral species. Second, a MutS homolog (mtMutS) has been found in the mitochondrial genome of all octocorals. To understand the origin and potential functional implications of these exceptions, we investigated the evolution of the MutS family both in Cnidaria and in animals in general. Our study confirmed the acquisition of octocoral mtMutS by horizontal gene transfer from a giant virus. Surprisingly, we identified MSH1 in all hexacorals and several sponges and placozoans. By contrast, MSH1 orthologs were lacking in other cnidarians, ctenophores, and bilaterian animals. Furthermore, while we identified MSH2 and MSH6 in nearly all animals, MSH4, MSH5, and, especially, MSH3 were missing in multiple species. Overall, our analysis revealed a dynamic evolution of the MutS family in animals, with multiple losses of MSH1, MSH3, some losses of MSH4 and MSH5, and a gain of the octocoral mtMutS. We propose that octocoral mtMutS functionally replaced MSH1 that was present in the common ancestor of Anthozoa.

Causes and consequences of mitochondrial proteome size variation in animals.

Despite a conserved set of core mitochondrial functions, animal mitochondrial proteomes show a large variation in size. We analyzed putative mechanisms behind and functional significance of this variation by performing comparative analysis of the experimentally-verified mitochondrial proteomes of four bilaterian animals (human, mouse, Caenorhabditis elegans, and Drosophila melanogaster) and two non-animal outgroups (Acanthamoeba castellanii and Saccharomyces cerevisiae). We found that of several factors affecting mitochondrial proteome size, evolution of novel mitochondrial proteins in mammals and loss of ancestral proteins in protostomes were the main contributors. Interestingly, the gain and loss of the N-terminal mitochondrial targeting signal was not a major factor in the proteome size evolution.

MMPdb and MitoPredictor: Tools for facilitating comparative analysis of animal mitochondrial proteomes.

Data on experimentally-characterized animal mitochondrial proteomes (mt-proteomes) are limited to a few model organisms and are scattered across multiple databases, impeding a comparative analysis. We developed two resources to address these problems. First, we re-analyzed proteomic data from six species with experimentally characterized mt-proteomes: animals (Homo sapiens, Mus musculus, Caenorhabditis elegans, and Drosophila melanogaster), and outgroups (Acanthamoeba castellanii and Saccharomyces cerevisiae) and created the Metazoan Mitochondrial Proteome Database (MMPdb) to host the results. Second, we developed a novel pipeline, "MitoPredictor" that uses a Random Forest classifier to infer mitochondrial localization of proteins based on orthology, mitochondrial targeting signal prediction, and protein domain analyses. Both tools generate an R Shiny applet that can be used to visualize and interact with the results and can be used on a personal computer. MMPdb is also available online at https://mmpdb.eeob.iastate.edu/.

Comparative study of human mitochondrial proteome reveals extensive protein subcellular relocalization after gene duplications.

BACKGROUND: Gene and genome duplication is the principle creative force in evolution. Recently, protein subcellular relocalization, or neolocalization was proposed as one of the mechanisms responsible for the retention of duplicated genes. This hypothesis received support from the analysis of yeast genomes, but has not been tested thoroughly on animal genomes. In order to evaluate the importance of subcellular relocalizations for retention of duplicated genes in animal genomes, we systematically analyzed nuclear encoded mitochondrial proteins in the human genome by reconstructing phylogenies of mitochondrial multigene families. RESULTS: The 456 human mitochondrial proteins selected for this study were clustered into 305 gene families including 92 multigene families. Among the multigene families, 59 (64%) consisted of both mitochondrial and cytosolic (non-mitochondrial) proteins (mt-cy families) while the remaining 33 (36%) were composed of mitochondrial proteins (mt-mt families). Phylogenetic analyses of mt-cy families revealed three different scenarios of their neolocalization following gene duplication: 1) relocalization from mitochondria to cytosol, 2) from cytosol to mitochondria and 3) multiple subcellular relocalizations. The neolocalizations were most commonly enabled by the gain or loss of N-terminal mitochondrial targeting signals. The majority of detected subcellular relocalization events occurred early in animal evolution, preceding the evolution of tetrapods. Mt-mt protein families showed a somewhat different pattern, where gene duplication occurred more evenly in time. However, for both types of protein families, most duplication events appear to roughly coincide with two rounds of genome duplications early in vertebrate evolution. Finally, we evaluated the effects of inaccurate and incomplete annotation of mitochondrial proteins and found that our conclusion of the importance of subcellular relocalization after gene duplication on the genomic scale was robust to potential gene misannotation. CONCLUSION: Our results suggest that protein subcellular relocalization is an important mechanism for the retention and gain of function of duplicated genes in animal genome evolution.

The mitochondrial genomes of sponges provide evidence for multiple invasions by Repetitive Hairpin-forming Elements (RHE).

BACKGROUND: The mitochondrial (mt) genomes of sponges possess a variety of features, which appear to be intermediate between those of Eumetazoa and non-metazoan opisthokonts. Among these features is the presence of long intergenic regions, which are common in other eukaryotes, but generally absent in Eumetazoa. Here we analyse poriferan mitochondrial intergenic regions, paying particular attention to repetitive sequences within them. In this context we introduce the mitochondrial genome of Ircinia strobilina (Lamarck, 1816; Demospongiae: Dictyoceratida) and compare it with mtDNA of other sponges. RESULTS: Mt genomes of dictyoceratid sponges are identical in gene order and content but display major differences in size and organization of intergenic regions. An even higher degree of diversity in the structure of intergenic regions was found among different orders of demosponges. One interesting observation made from such comparisons was of what appears to be recurrent invasions of sponge mitochondrial genomes by repetitive hairpin-forming elements, which cause large genome size differences even among closely related taxa. These repetitive hairpin-forming elements are structurally and compositionally divergent and display a scattered distribution throughout various groups of demosponges. CONCLUSION: Large intergenic regions of poriferan mt genomes are targets for insertions of repetitive hairpin- forming elements, similar to the ones found in non-metazoan opisthokonts. Such elements were likely present in some lineages early in animal mitochondrial genome evolution but were subsequently lost during the reduction of intergenic regions, which occurred in the Eumetazoa lineage after the split of Porifera. Porifera acquired their elements in several independent events. Patterns of their intra-genomic dispersal can be seen in the mt genome of Vaceletia sp.

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae).

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.