Mammalian Mitochondrial Genes Research Articles

Structural basis for shape-selective recognition and aminoacylation of a D-armless human mitochondrial tRNA

Human mitochondrial gene expression relies on the specific recognition and aminoacylation of mitochondrial tRNAs (mtRNAs) by nuclear-encoded mitochondrial aminoacyl-tRNA synthetases (mt-aaRSs). Despite their essential role in cellular energy homeostasis, strong mutation pressure and genetic drift have led to an unparalleled sequence erosion of animal mtRNAs. The structural and functional consequences of this erosion are not understood. Here, we present cryo-EM structures of the human mitochondrial seryl-tRNA synthetase (mSerRS) in complex with mtRNASer(GCU). These structures reveal a unique mechanism of substrate recognition and aminoacylation. The mtRNASer(GCU) is highly degenerated, having lost the entire D-arm, tertiary core, and stable L-shaped fold that define canonical tRNAs. Instead, mtRNASer(GCU) evolved unique structural innovations, including a radically altered T-arm topology that serves as critical identity determinant in an unusual shape-selective readout mechanism by mSerRS. Our results provide a molecular framework to understand the principles of mito-nuclear co-evolution and specialized mechanisms of tRNA recognition in mammalian mitochondrial gene expression.

Applying frequency chaos game representation with perceptual image hashing to gene sequence phylogenetic analyses

As a very important research direction in the field of bioinformatics, sequence alignment plays a vital role in the research and development of biology. Converting genome sequence to graph by using frequency chaos game representation (FCGR) is an excellent gene sequence mapping technology, which can store rich genetic information into FCGR graphics. To each FCGR image, we construct its perceptual image hashing (PIH) matrix using the bicubic interpolation zooming. The difference of the perceptual hash matrix of each two images is calculated, and the clustering distance of the corresponding two gene sequences is represented by the differentials of the perceptual hash matrix. In this paper, we aligned and analyzed several typical genome sequence datasets including mammalian mitochondrial genes, human immunodeficiency virus 1 (HIV-1) and hepatitis E virus (HEV) to build their evolutionary trees. Experimental results showed that our PIH combining FCGR method (FCGR-PIH) has similar classification accuracy to the classical Clustal W sequence alignment method. Furthermore, 25 complete mitochondrial DNA sequences of cichlid fishes and 27 Escherichia coli/Shigella full genome sequences were selected from the AFproject test platform for tests. The performance benchmark rankings demonstrate the effectiveness of the FCGR-PIH algorithm and its potential for large-scale genome sequence analysis.

Next-generation sequencing profiling of mitochondrial genomes in gout

BackgroundAccumulating evidence implicates mitochondrial DNA (mtDNA) alleles, which are independent of the nuclear genome, in disease, especially in human metabolic diseases. However, this area of investigation has lagged behind in researching the nuclear alleles in complex traits, for example, in gout.MethodsNext-generation sequencing was utilized to investigate the relationship between mtDNA alleles and phenotypic variations in 52 male patients with gout and 104 age-matched male non-gout controls from the Taiwan Biobank whole-genome sequencing samples. Differences from a reference sequence (GRCh38) were identified. The sequence kernel association test (SKAT) was applied to identify gout-associated alleles in mitochondrial genes. The tools Polymorphism Phenotyping, Sorting Intolerant From Tolerant (SIFT), Predict the pathology of Mutations (PMUT), Human Mitochondrial Genome Database (mtDB), Multiple Alignment using Fast Fourier Transform (MAFFT), and Mammalian Mitochondrial tRNA Genes (Mamit-tRNA) were used to evaluate pathogenicity of alleles. Validation of selected alleles by quantitative polymerase chain reaction of single nucleotide polymorphisms (qPCR SNPs) was also performed.ResultsWe identified 456 alleles in patients with gout and 640 alleles in non-gout controls with 274 alleles shared by both. Mitochondrial genes were associated with gout, with MT-CO3, MT-TA, MT-TC, and MT-TT containing potentially pathogenic gout-associated alleles and displaying evidence of gene-gene interactions. All heteroplasmy levels of potentially pathogenic alleles exceeded metabolic thresholds for pathogenicity. Validation assays confirmed the next-generation sequencing results of selected alleles. Among them, potentially pathogenic MT-CO3 alleles correlated with high-density lipoprotein (HDL) levels (P = 0.034).ConclusionThis study provided two scientific insights. First, this was the most extensive mitochondrial genomic profiling associated with gout. Second, our results supported the roles of mitochondria in gout and HDL, and this comprehensive analysis framework can be applied to other diseases in which mitochondrial dysfunction has been implicated.

Post-fertilisation sperm mitophagy: the tale of Mitochondrial Eve and Steve.

Preformationist William Harvey's proclamation of everything live coming from an egg still holds true for mammalian mitochondria and mitochondrial genes. At fertilisation, mitochondria carried into the oocyte cytoplasm by the spermatozoon are sought out and destroyed, leaving only oocyte mitochondria to propagate their mitochondrial (mt) DNA to offspring. This clonal inheritance mode, the 'mitochondrial Eve' paradigm, is mediated by oocytes' resident proteolytic, organelle-targeting mechanisms, including the substrate-specific ubiquitin proteasome system and the autophagic machinery for bulk protein and organelle degradation. Ubiquitination of sperm mitochondria within the cytoplasm of the fertilised oocyte was initially discovered in mammals. More recent studies in Drosophila and Caenorhabditis elegans implicated the ubiquitin-binding autophagy protein sequestosome 1 (SQSTM1) as the early adaptor channelling ubiquitinated sperm mitochondria towards the autophagic machinery. Downstream receptors include microtubule-associated protein 1 light chain 3α (LC3) and GABA type A receptor-associated protein (GABARAP). Among mammals, the domestic pig is the ideal mammalian model of mitochondrial inheritance because of rapid sperm mitophagy at the 1-cell stage of embryo development. Primary recognition of sperm mitochondria by SQSTM1 inside the porcine zygote is followed by GABARAP-containing autophagophore formation, and contributed to by valosin-containing protein (VCP), a 26S proteasome-presenting protein dislocase. Consequently, coinhibition of SQSTM1-GABARAP and VCP activities in the porcine zygotes, resulting in 2- to 4-cell embryos carrying intact sperm mitochondrial sheaths, revived the moniker of 'Mitochondrial Steve'. Further work will identify the determinants of species specificity of sperm mitophagy and explain the interplay and possible consequences of a mismatch between clonal mitochondrial genome and biparentally inherited chromosomal genes encoding for structural mitochondrial proteins and transcription factors. By better understanding sperm mitophagy and its potential failure, we may be able to alleviate mitochondrial disease and early pregnancy loss in livestock and improve their fitness, reproduction and ability to pass favourable production traits to offspring.

Complete mitochondrial genome of the stoat (Mustela erminea) and New Zealand fur seal (Arctocephalus forsteri) and their significance for mammalian phylogeny

The complete mitochondrial genome of three mustelid species, stoats (Mustela erminea), weasels (Mustela nivalis) and ferrets (Mustela furo), and the New Zealand fur seal (Arctocephalus forsteri) were sequenced using direct mitochondrial DNA extraction and overlapping long PCRs. The usual 37 mammalian mitochondrial genes (13 protein coding genes, 22 t-RNA and 2 r-RNA) were identified in all four mitogenomes. The divergence of stoats from other members of the sub-family Mustelinae was dated 4.5 million years ago. The mitogenomic data were consistent with a bear-like origin of seals.

Purifying Selection in Mammalian Mitochondrial Protein-Coding Genes Is Highly Effective and Congruent with Evolution of Nuclear Genes

The mammalian mitochondrial genomes differ from the nuclear genomes by maternal inheritance, absence of recombination, and higher mutation rate. All these differences decrease the effective population size of mitochondrial genome and make it more susceptible to accumulation of slightly deleterious mutations. It was hypothesized that mitochondrial genes, especially in species with low effective population size, irreversibly degrade leading to decrease of organismal fitness and even to extinction of species through the mutational meltdown. To interrogate this hypothesis, we compared the purifying selections acting on the representative set of mitochondrial (potentially degrading) and nuclear (potentially not degrading) protein-coding genes in species with different effective population size. For 21 mammalian species, we calculated the ratios of accumulation of slightly deleterious mutations approximated by Kn/Ks separately for mitochondrial and nuclear genomes. The 75% of variation in Kn/Ks is explained by two independent variables: type of a genome (mitochondrial or nuclear) and effective population size of species approximated by generation time. First, we observed that purifying selection is more effective in mitochondria than in the nucleus that implies strong evolutionary constraints of mitochondrial genome. Mitochondrial de novo nonsynonymous mutations have at least 5-fold more harmful effect when compared with nuclear. Second, Kn/Ks of mitochondrial and nuclear genomes is positively correlated with generation time of species, indicating relaxation of purifying selection with decrease of species-specific effective population size. Most importantly, the linear regression lines of mitochondrial and nuclear Kn/Ks's from generation times of species are parallel, indicating congruent relaxation of purifying selection in both genomes. Thus, our results reveal that the distribution of selection coefficients of de novo nonsynonymous mitochondrial mutations has a similar shape with the distribution of de novo nonsynonymous nuclear mutations, but its mean is five times smaller. The harmful effect of mitochondrial de novo nonsynonymous mutations triggers highly effective purifying selection, which maintains the fitness of the mammalian mitochondrial genome.

Highlights from the Literature

Highlights from the Literature

The role of mammalian PPR domain proteins in the regulation of mitochondrial gene expression

Pentatricopeptide repeat (PPR) domain proteins are a large family of RNA-binding proteins that are involved in the maturation and translation of organelle transcripts in eukaryotes. They were first identified in plant organelles and their important role in mammalian mitochondrial gene regulation is now emerging. Mammalian PPR proteins, like their plant counterparts, have diverse roles in mitochondrial transcription, RNA metabolism and translation and consequently are important for mitochondrial function and cell health. Here we discuss the current knowledge about the seven mammalian PPR proteins identified to date and their roles in the regulation of mitochondrial gene expression. Furthermore we discuss the mitochondrial RNA targets of the mammalian PPR proteins and methods to investigate the RNA targets of these mitochondrial RNA-binding proteins. This article is part of a Special Issue entitled: Mitochondrial Gene Expression.

Mitochondrial transformation: time for concerted action

We are frustratingly devoid of robust methods to genetically transform mitochondria in nearly all species. This extended Opinion explores the attempts to date, and their successes or failures.

Bayesian analysis of amino acid substitution models

Models of amino acid substitution present challenges beyond those often faced with the analysis of DNA sequences. The alignments of amino acid sequences are often small, whereas the number of parameters to be estimated is potentially large when compared with the number of free parameters for nucleotide substitution models. Most approaches to the analysis of amino acid alignments have focused on the use of fixed amino acid models in which all of the potentially free parameters are fixed to values estimated from a large number of sequences. Often, these fixed amino acid models are specific to a gene or taxonomic group (e.g. the Mtmam model, which has parameters that are specific to mammalian mitochondrial gene sequences). Although the fixed amino acid models succeed in reducing the number of free parameters to be estimated--indeed, they reduce the number of free parameters from approximately 200 to 0--it is possible that none of the currently available fixed amino acid models is appropriate for a specific alignment. Here, we present four approaches to the analysis of amino acid sequences. First, we explore the use of a general time reversible model of amino acid substitution using a Dirichlet prior probability distribution on the 190 exchangeability parameters. Second, we then explore the behaviour of prior probability distributions that are'centred' on the rates specified by the fixed amino acid model. Third, we consider a mixture of fixed amino acid models. Finally, we consider constraints on the exchangeability parameters as partitions,similar to how nucleotide substitution models are specified, and place a Dirichlet process prior model on all the possible partitioning schemes.

A large-scale screening of the normalized mammalian mitochondrial gene expression profiles

Mammalian mitochondrial genomes are organized in a conserved and extremely compact manner, encoding molecules that play a vital role in oxidative phosphorylation (OXPHOS) and carry out a number of other important biological functions. A large-scale screening of the normalized mitochondrial gene expression profiles generated from publicly available mammalian serial analysis of gene expression (SAGE) datasets (over 17.7 millions of tags) was performed in this study. Acquired SAGE libraries represent an extensive range of human, mouse, rat, bovine and swine cell and tissue samples (normal and pathological) in a variety of conditions. Using a straightforward in silico algorithm, variations in total mitochondrial gene expression, as well as in the expression of individual genes encoded by mitochondrial genomes are addressed, and common patterns in the species- and tissue-specific mitochondrial gene expression profiles are discussed.

Forces maintaining organellar genomes: is any as strong as genetic code disparity or hydrophobicity?

It remains controversial why mitochondria and chloroplasts retain the genes encoding a small subset of their constituent proteins, despite the transfer of so many other genes to the nucleus. Two candidate obstacles to gene transfer, suggested long ago, are that the genetic code of some mitochondrial genomes differs from the standard nuclear code, such that a transferred gene would encode an incorrect amino acid sequence, and that the proteins most frequently encoded in mitochondria are generally very hydrophobic, which may impede their import after synthesis in the cytosol. More recently it has been suggested that both these interpretations suffer from serious "false positives" and "false negatives": genes that they predict should be readily transferred but which have never (or seldom) been, and genes whose transfer has occurred often or early, even though this is predicted to be very difficult. Here I consider the full known range of ostensibly problematic such genes, with particular reference to the sequences of events that could have led to their present location. I show that this detailed analysis of these cases reveals that they are in fact wholly consistent with the hypothesis that code disparity and hydrophobicity are much more powerful barriers to functional gene transfer than any other. The popularity of the contrary view has led to the search for other barriers that might retain genes in organelles even more powerfully than code disparity or hydrophobicity; one proposal, concerning the role of proteins in redox processes, has received widespread support. I conclude that this abandonment of the original explanations for the retention of organellar genomes has been premature. Several other, relatively minor, obstacles to gene transfer certainly exist, contributing to the retention of relatively many organellar genes in most lineages compared to animal mtDNA, but there is no evidence for obstacles as severe as code disparity or hydrophobicity. One corollary of this conclusion is that there is currently no reason to suppose that engineering nuclear versions of the remaining mammalian mitochondrial genes, a feat that may have widespread biomedical relevance, should require anything other than sequence alterations obviating code disparity and causing modest reductions in hydrophobicity without loss of enzymatic function.

Phosphorylation of rat mitochondrial transcription termination factor (mTERF) is required for transcription termination but not for binding to DNA.

Despite the crucial importance of mitochondrial transcription, knowledge of its regulation is poor. Therefore, characterization of mammalian mitochondrial transcription termination factor (mTERF) functionality and regulation is of fundamental biological interest in order to understand the regulation of mitochondrial transcription. Here we report that mTERF is the first protein having a role in mammalian mitochondrial gene expression that appears to be controlled by phosphorylation. Recombinant mature rat mTERF protein has specific DNA-binding capacity for the sequence required for transcription termination. Furthermore, unlike recombinant human mTERF, the rat protein bound to its mitochondrial DNA binding site promotes the termination of transcription initiated with heterologous RNA polymerase. Interestingly, mTERF is a phosphoprotein with four phosphate groups, and while the DNA-binding activity of mTERF is unaffected by the phosphorylation/dephosphorylation state, only the phosphorylated form of the protein is active for termination activity. Moreover, natural human mTERF is also a phosphoprotein and its termination activity is inhibited by dephosphorylation. These data suggest that mTERF functioning in vivo is regulated by phosphorylation.

Root factors induce mitochondrial-related gene expression and fungal respiration during the developmental switch from asymbiosis to presymbiosis in the arbuscular mycorrhizal fungus Gigaspora rosea.

During spore germination, arbuscular mycorrhizal (AM) fungi show limited hyphal development in the absence of a host plant (asymbiotic). In the presence of root exudates, they switch to a new developmental stage (presymbiotic) characterized by extensive hyphal branching. Presymbiotic branching of the AM fungus Gigaspora rosea was induced in liquid medium by a semipurified exudate fraction from carrot (Daucus carota) root organ cultures. Changes in RNA accumulation patterns were monitored by differential display analysis. Differentially appearing cDNA fragments were cloned and further analyzed. Five cDNA fragments could be identified that show induced RNA accumulation 1 h after the addition of root exudate. Sequence similarities of two fragments to mammalian Nco4 and mitochondrial rRNA genes suggested that root exudates could influence fungal respiratory activity. To support this hypothesis, additional putative mitochondrial related-genes were shown to be induced by root exudates. These genes were identified after subtractive hybridization and putatively encode a pyruvate carboxylase and a mitochondrial ADP/ATP translocase. The gene GrosPyc1 for the pyruvate carboxylase was studied in more detail by cloning a cDNA and by quantifying its RNA accumulation. The hypothesis that respiratory activity of AM fungi is stimulated by root exudates was confirmed by physiological and cytological analyses in G. rosea and Glomus intraradices. Oxygen consumption and reducing activity of both fungi was induced after 3 and 2 h of exposition with the root factor, respectively, and the first respiration activation was detected in G. intraradices after approximately 90 min. In addition, changes in mitochondrial morphology, orientation, and overall biomass were detected in G. rosea after 4 h. In summary, the root-exuded factor rapidly induces the expression of certain fungal genes and, in turn, fungal respiratory activity before intense branching. This defines the developmental switch from asymbiosis to presymbiosis, first by gene activation (0.5-1 h), subsequently on the physiological level (1.5-3 h), and finally as a morphological response (after 5 h).

Partitioning the variation in mammalian substitution rates.

We have used analysis of variance to partition the variation in synonymous and amino acid substitution rates between three effects (gene, lineage, and a gene-by-lineage interaction) in mammalian nuclear and mitochondrial genes. We find that gene effects are stronger for amino acid substitution rates than for synonymous substitution rates and that lineage effects are stronger for synonymous substitution rates than for amino acid substitution rates. Gene-by-lineage interactions, equivalent to overdispersion corrected for lineage effects, are found in amino acid substitutions but not in synonymous substitutions. The variance in the ratio of amino acid and synonymous substitution rates is dominated by gene effects, but there is also a significant gene-by-lineage interaction.

Sources of Incongruence among Mammalian Mitochondrial Sequences: COII, COIII, and ND6 Genes Are Main Contributors

To investigate the origins of incongruence among mammalian mitochondrial protein-coding genes, we compiled a matrix that included 13 protein-coding-genes for 41 mammals from 14 different orders. This matrix was examined for congruence using different partitioning strategies. The incongruence length difference test showed significant incongruence among the 13 gene partitions used simultaneously, and the result was not affected by third codon or transversion weighting. In the pair-wise comparisons, significant incongruence was detected between NADH:ubiquinone oxidoreductase subunit 6 gene (ND6), cytochrome oxidase subunit II (COII), or cytochrome oxidase subunit III (COIII) gene partitioned individually against the rest of the genes. Omission of any of the 14 mammalian orders alone or in combinations from the matrix did not result in a statistically significant improvement of congruence, suggesting that taxonomic sampling will not improve congruence among the data sets. However, omission of the ND6, COII, and COIII significantly improved congruence in our data matrix. Possible origins of unusual phylogenetic properties of the three genes are discussed.

An araC-controlled Bacterialcre Expression System to Produce DNA Minicircle Vectors for Nuclear and Mitochondrial Gene Therapy

The presence of CpG motifs and their associated sequences in bacterial DNA causes an immunotoxic response following the delivery of these plasmid vectors into mammalian hosts. We describe a biotechnological approach to the elimination of this problem by the creation of a bacterial cre recombinase expression system, tightly controlled by the arabinose regulon. This permits the Cre-mediated and -directed excision of the entire bacterial vector sequences from plasmid constructs to create supercoiled gene expression minicircles for gene therapy. Minicircle yields using standard culture volumes are sufficient for most in vitro and in vivo applications whereas minicircle expression in vitro is significantly increased over standard plasmid transfection. By the simple expedient of removing the bacterial DNA complement, we significantly reduce the size and CpG content of these expression vectors, which should also reduce DNA-induced inflammatory responses in a dose-dependent manner. We further describe the generation of minicircle expression vectors for mammalian mitochondrial gene therapy, for which no other vector systems currently exist. The removal of bacterial vector sequences should permit appropriate transcription and correct transcriptional cleavage from the mitochondrial minicircle constructs in a mitochondrial environment and brings the realization of mitochondrial gene therapy a step closer.

Comments on Allard and Carpenter (1996), or the “Aquatic Ape” Hypothesis Revisited

The results of Allard and Carpenter's (Cladistics 12, 183–198, 1996) paper on weighting and congruence among mammalian mitochondrial genes are an artefact of errors in their data matrix; their “blue whale” ATPASE8 sequence is human, the actual blue whale sequence is assigned to the grey seal, and the “horse” sequence is that of the harbor seal. When these errors are corrected there is no evidence that the mitochondrial genes are incongruent.

The rate heterogeneity of nonsynonymous substitutions in mammalian mitochondrial genes.

Substitution rates at the three codon positions (r1, r2, and r3) of mammalian mitochondrial genes are in the order of r3 > r1 > r2, and the rate heterogeneity at the three positions, as measured by the shape parameter of the gamma distribution (alpha 1, alpha 2, and alpha 3), is in the order of alpha 3 > alpha 1 > alpha 2. The causes for the rate heterogeneity at the three codon positions remain unclear and, in particular, there has been no satisfactory explanation for the observation of alpha 1 > alpha 2. I attempted to dissect the causes of rate heterogeneity by studying the pattern of nonsynonymous substitutions with respect to codon positions in 10 mitochondrial genes from 19 mammalian species. Nonsynonymous substitutions involve more different amino acid replacements at the second than at the first codon position, which results in r1 > r2. The difference between r1 and r2 increases with the intensity of purifying selection, and so does the rate heterogeneity in nonsynonymous substitutions among sites at the same codon position. All mitochondrial genes appear to have functionally important and unimportant codons, with the latter having all three codon positions prone to nonsynonymous substitutions. Within the functionally important codons, the second codon position is much more conservative than the codon position. This explains why alpha 1 > alpha 2. The result suggests that overweighting of the second codon position in phylogenetic analysis may be a misguided practice.

Estimating the transition/transversion ratio from independent pairwise comparisons with an assumed phylogeny.

A method is presented for estimating the transition/transversion ratio (TI/TV), based on phylogenetically independent comparisons. TI/TV is a parameter of some models used in phylogeny estimation intended to reflect the fact that nucleotide substitutions are not all equally likely. Previous attempts to estimate TI/TV have commonly faced three problems: (1) few taxa; (2) nonindependence among pairwise comparisons; and (3) multiple hits make the apparent TI/TV between two sequences decrease over time since their divergence, giving a misleading impression of relative substitution probabilities. We have made use of the time dependency, modeling how the observed TI/TV changes over time and extrapolating to estimate the "instantaneous" TI/TV-the relevant parameter for phylogenetic inference. To illustrate our method, TI/TV was estimated for two mammalian mitochondrial genes. For 26 pairs of cytochrome b sequences, the estimate of TI/TV was 5.5; 16 pairs of 12s rRNA yielded an estimate of 9.5. These estimates are higher than those given by the maximum likelihood method and than those obtained by averaging all possible pairwise comparisons (with or without a two-parameter correction for multiple substitutions). We discuss strengths, weaknesses, and further uses of our method.