Search

Abstract Background Schizophrenia is associated with rare copy-number (CN) mutations. Screening for such alleles genome-wide, though comprehensive, cannot study in-depth the causality of particular loci, therefore cannot provide the functional interpretation for the disease etiology. We hypothesized that CN mutations in the SELENBP1 locus could associate with the disorder and that these mutations could alter the gene product's activity in patients. Methods We analyzed SELENBP1 CN variation (CNV) in blood DNA from 49 schizophrenia patients and 49 controls (cohort A). Since CN of genes may vary among tissues, we investigated SELENBP1 CN in age- sex- and postmortem interval-matched cerebellar DNA samples from 14 patients and 14 controls (cohort B). Since CNV may either be de-novo or inherited we analyzed CNV of the SELENBP1 locus in blood DNA from 26 trios of schizophrenia probands and their healthy parents (cohort C). SELENBP1 mRNA levels were measured by real-time PCR. Results In cohort A reduced CN of the SELENBP1 locus was found in four patients but in none of the controls. In cohort B we found reduced CN of the SELENBP1 locus in two patients but in none of the controls. In cohort C three patients exhibited drastic CN reduction, not present in their parents, indicating de-novo mutation. A reduction in SELENBP1 mRNA levels in the postmortem cerebellar samples of schizophrenia patients was found. Conclusions We report a focused study of CN mutations in the selenium binding-protein1 (SELENBP1) locus previously linked with schizophrenia. We provide evidence for recurrence of decreased CN of the SELENBP1 locus in three unrelated patients' cohorts but not in controls, raising the possibility of functional involvement of these mutations in the etiology of the disease.

Abstract Background Although mitochondrial dysfunction is consistently manifested in patients with Type 2 Diabetes mellitus (T2DM), the association of mitochondrial DNA (mtDNA) sequence variants with T2DM varies among populations. These differences might stem from differing environmental influences among populations. However, other potentially important considerations emanate from the very nature of mitochondrial genetics, namely the notable high degree of partitioning in the distribution of human mtDNA variants among populations, as well as the interaction of mtDNA and nuclear DNA-encoded factors working in concert to govern mitochondrial function. We hypothesized that association of mtDNA genetic variants with T2DM could be revealed while controlling for the effect of additional inherited factors, reflected in family history information. Methods To test this hypothesis we set out to investigate whether mtDNA genetic variants will be differentially associated with T2DM depending on the diabetes status of the parents. To this end, association of mtDNA genetic backgrounds (haplogroups) with T2DM was assessed in 1055 Jewish patients with and without T2DM parents ('DP' and 'HP', respectively). Results Haplogroup J1 was found to be 2.4 fold under-represented in the 'HP' patients (p = 0.0035). These results are consistent with a previous observation made in Finnish T2DM patients. Moreover, assessing the haplogroup distribution in 'DP' versus 'HP' patients having diabetic siblings revealed that haplogroup J1 was virtually absent in the 'HP' group. Conclusion These results imply the involvement of inherited factors, which modulate the susceptibility of haplogroup J1 to T2DM.

Abstract Background For survival, scorpions depend on a wide array of short neurotoxic polypeptides. The venoms of scorpions from the most studied group, the Buthida, are a rich source of small, 23–78 amino acid-long peptides, well packed by either three or four disulfide bridges that affect ion channel function in excitable and non-excitable cells. Results In this work, by constructing a toxin transcripts data set from the venom gland of the scorpion Buthus occitanus israelis, we were able to follow the evolutionary path leading to mature toxin diversification and suggest a mechanism for leader peptide hyper-conservation. Toxins from each family were more closely related to one another than to toxins from other species, implying that fixation of duplicated genes followed speciation, suggesting early gene conversion events. Upon fixation, the mature toxin-coding domain was subjected to diversifying selection resulting in a significantly higher substitution rate that can be explained solely by diversifying selection. In contrast to the mature peptide, the leader peptide sequence was hyper-conserved and characterized by an atypical sub-neutral synonymous substitution rate. We interpret this as resulting from purifying selection acting on both the peptide and, as reported here for the first time, the DNA sequence, to create a toxin family-specific codon bias. Conclusion We thus propose that scorpion toxin genes were shaped by selective forces acting at three levels, namely (1) diversifying the mature toxin, (2) conserving the leader peptide amino acid sequence and intriguingly, (3) conserving the leader DNA sequences.

Abstract Background Recent genome-wide association studies searching for candidate susceptibility loci for common complex diseases such as type 2 diabetes mellitus (T2DM) and its common complications have uncovered novel disease-associated genes. Nevertheless these large-scale population screens often overlook the tremendous variation in the mitochondrial genome (mtDNA) and its involvement in complex disorders. Results We have analyzed the mitochondrial DNA (mtDNA) genetic variability in Ashkenazi (Ash), Sephardic (Seph) and North African (NAF) Jewish populations (total n = 1179). Our analysis showed significant differences (p < 0.001) in the distribution of mtDNA genetic backgrounds (haplogroups) among the studied populations. To test whether these differences alter the pattern of disease susceptibility, we have screened our three Jewish populations for an association of mtDNA genetic haplogroups with T2DM complications. Our results identified population-specific susceptibility factors of which the best example is the Ashkenazi Jewish specific haplogroup N1b1, having an apparent protective effect against T2DM complications in Ash (p = 0.006), being absent in the NAF population and under-represented in the Seph population. We have generated and analyzed whole mtDNA sequences from the disease associated haplogroups revealing mutations in highly conserved positions that are good candidates to explain the phenotypic effect of these genetic backgrounds. Conclusion Our findings support the possibility that recent bottleneck events leading to over-representation of minor mtDNA alleles in specific genetic isolates, could result in population-specific susceptibility loci to complex disorders.

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation., (© 2023. The Author(s).)

RNA-DNA differences (RDD) have previously been identified in the human mitochondrial RNA (mt-RNA) transcripts, yet their functional impact is poorly understood. By analyzing 4928 RNA-seq samples from 23 body sites, we found that mtDNA gene expression negatively correlated with the levels of both m 1 A 947 16 S rRNA modification (mtDNA position 2617) and the m 1 A 1812 ND5 mRNA modification (mtDNA position 13,710) in 15 and 14 body sites, respectively. Such correlation was not evident in all tested brain tissues, thus suggesting a tissue-specific impact of these modifications on mtDNA gene expression. To assess the response of the tested modifications to environmental cues, we analyzed pairs of skin samples that were either exposed to the sun or not. We found that the correlations of mtDNA gene expression with both mt-RNA modifications were compromised upon sun exposure. As a first step to explore the underlying mechanism, we analyzed RNA-seq data from keratinocytes that were exposed to increasing doses of UV irradiation. Similar to sun exposure, we found a significant decrease in mtDNA gene expression upon increase in UV dosage. In contrast, there was a significant increase in the m 1 A 947 16 S rRNA modification levels upon elevation in UV dose. Finally, we identified candidate modulators of such responses. Taken together, our results indicate that mt-RNA modifications functionally correlate with mtDNA gene expression, and responds to environmental cues, hence supporting their physiological importance., (© 2022. The Author(s), under exclusive licence to European Society of Human Genetics.)

Background: Mitochondrial DNA (mtDNA) is present at high copy numbers in animal cells, and though characterized by a single haplotype in each individual due to maternal germline inheritance, deleterious mutations and intact mtDNA molecules frequently co-exist (heteroplasmy). A number of factors, such as replicative segregation, mitochondrial bottlenecks, and selection, may modulate the exitance of heteroplasmic mutations. Since such mutations may have pathological consequences, they likely survive and are inherited due to functional complementation via the intracellular mitochondrial network. Here, we hypothesized that compromised mitochondrial fusion would hamper such complementation, thereby affecting heteroplasmy inheritance., Results: We assessed heteroplasmy levels in three Caenorhabditis elegans strains carrying different heteroplasmic mtDNA deletions (ΔmtDNA) in the background of mutant mitofusin (fzo-1). Animals displayed severe embryonic lethality and developmental delay. Strikingly, observed phenotypes were relieved during subsequent generations in association with complete loss of ΔmtDNA molecules. Moreover, deletion loss rates were negatively correlated with the size of mtDNA deletions, suggesting that mitochondrial fusion is essential and sensitive to the nature of the heteroplasmic mtDNA mutations. Introducing the ΔmtDNA into a fzo-1;pdr-1;+/ΔmtDNA (PARKIN ortholog) double mutant resulted in a skewed Mendelian progeny distribution, in contrast to the normal distribution in the fzo-1;+/ΔmtDNA mutant, and severely reduced brood size. Notably, the ΔmtDNA was lost across generations in association with improved phenotypes., Conclusions: Taken together, our findings show that when mitochondrial fusion is compromised, deleterious heteroplasmic mutations cannot evade natural selection while inherited through generations. Moreover, our findings underline the importance of cross-talk between mitochondrial fusion and mitophagy in modulating the inheritance of mtDNA heteroplasmy., (© 2022. The Author(s).)

Mitochondria are pivotal for bioenergetics, as well as in cellular response to viral infections. Nevertheless, their role in COVID-19 was largely overlooked. Here, we analyzed available bulk RNA-seq datasets from COVID-19 patients and corresponding healthy controls (three blood datasets, N = 48 healthy, 119 patients; two respiratory tract datasets, N = 157 healthy, 524 patients). We found significantly reduced mtDNA gene expression in blood, but not in respiratory tract samples from patients. Next, analysis of eight single-cells RNA-seq datasets from peripheral blood mononuclear cells, nasopharyngeal samples, and Bronchoalveolar lavage fluid (N = 1,192,243 cells), revealed significantly reduced mtDNA gene expression especially in immune system cells from patients. This is associated with elevated expression of nuclear DNA-encoded OXPHOS subunits, suggesting compromised mitochondrial-nuclear co-regulation. This, together with elevated expression of ROS-response genes and glycolysis enzymes in patients, suggest rewiring toward glycolysis, thus generating beneficial conditions for SARS-CoV-2 replication. Our findings underline the centrality of mitochondrial dysfunction in COVID-19., Competing Interests: The authors declare no competing interests. The manuscript is based on analyses of publicly available data. The source of each dataset is mentioned within the STAR Methods section., (© 2021 The Authors.)

Out of many intracellular bacteria, only the mitochondria and chloroplasts abandoned their independence billions of years ago and became endosymbionts within the host eukaryotic cell. Consequently, one cannot grow eukaryotic cells without their mitochondria, and the mitochondria cannot divide outside of the cell, thus reflecting interdependence. Here, we argue that such interdependence underlies the fundamental role of mitochondrial activities in the emergence of metazoans. Several lines of evidence support our hypothesis: ( a ) Differentiation and embryogenesis rely on mitochondrial function; ( b ) mitochondrial metabolites are primary precursors for epigenetic modifications (such as methyl and acetyl), which are critical for chromatin remodeling and gene expression, particularly during differentiation and embryogenesis; and ( c ) mitonuclear coregulation adapted to accommodate both housekeeping and tissue-dependent metabolic needs. We discuss the evolution of the unique mitochondrial genetic system, mitochondrial metabolites, mitonuclear coregulation, and their critical roles in the emergence of metazoans and in human disorders.

In the Australian redclaw crayfish, Cherax quadricarinatus (WZ/ZZ system), intersexuals, although exhibiting both male and female gonopores, are functional males bearing a female genotype (WZ males). Therefore, the occurrence of the unusual homogametic WW females in nature is plausible. We developed W/Z genomic sex markers and used them to investigate the genotypic structure of experimental and native C. quadricarinatus populations in Australia. We discovered, for the first time, the natural occurrence of WW females in crustacean populations. By modeling population dynamics, we found that intersexuals contribute to the growth rate of crayfish populations in the short term. Given the vastly fragmented C. quadricarinatus habitat, which is characterized by drought-flood cycles, we speculate that intersexuals contribute to the fitness of this species since they lead to occasional increment in the population growth rate which potentially supports crayfish population restoration and establishment under extinction threats or colonization events., Competing Interests: A patent regarding sex-specific genomic markers in the Australian redclaw crayfish is pending (International application number: PCT/IL2018/051046, International publication number: WO/2019/058371)., (© 2020 The Authors.)

The higher order organization of eukaryotic and prokaryotic genomes is pivotal in the regulation of gene expression. Specifically, chromatin accessibility in eukaryotes and nucleoid accessibility in bacteria are regulated by a cohort of proteins to alter gene expression in response to diverse physiological conditions. By contrast, prior studies have suggested that the mitochondrial genome (mtDNA) is coated solely by mitochondrial transcription factor A (TFAM), whose increased cellular concentration was proposed to be the major determinant of mtDNA packaging in the mitochondrial nucleoid. Nevertheless, recent analysis of DNase-seq and ATAC-seq experiments from multiple human and mouse samples suggest gradual increase in mtDNA occupancy during the course of embryonic development to generate a conserved footprinting pattern which correlate with sites that have low TFAM occupancy in vivo (ChIP-seq) and tend to adopt G-quadruplex structures. These findings, along with recent identification of mtDNA binding by known modulators of chromatin accessibility such as MOF, suggest that mtDNA higher order organization is generated by cross talk with the nuclear regulatory system, may have a role in mtDNA regulation, and is more complex than once thought., (Copyright © 2019 Mishmar, Levin, Naeem and Sondheimer.)

Mitochondrial complex I (CI) is the largest multi-subunit oxidative phosphorylation (OXPHOS) protein complex. Recent availability of a high-resolution human CI structure, and from two non-human mammals, enabled predicting the impact of mutations on interactions involving each of the 44 CI subunits. However, experimentally assessing the impact of the predicted interactions requires an easy and high-throughput method. Here, we created such a platform by cloning all 37 nuclear DNA (nDNA) and 7 mitochondrial DNA (mtDNA)-encoded human CI subunits into yeast expression vectors to serve as both 'prey' and 'bait' in the split murine dihydrofolate reductase (mDHFR) protein complementation assay (PCA). We first demonstrated the capacity of this approach and then used it to examine reported pathological OXPHOS CI mutations that occur at subunit interaction interfaces. Our results indicate that a pathological frame-shift mutation in the MT-ND2 gene, causing the replacement of 126 C-terminal residues by a stretch of only 30 amino acids, resulted in loss of specificity in ND2-based interactions involving these residues. Hence, the split mDHFR PCA is a powerful assay for assessing the impact of disease-causing mutations on pairwise protein-protein interactions in the context of a large protein complex, thus offering a possible mechanistic explanation for the underlying pathogenicity.

Unlike the nuclear genome, the mammalian mitochondrial genome (mtDNA) is thought to be coated solely by mitochondrial transcription factor A (TFAM), whose binding sequence preferences are debated. Therefore, higher-order mtDNA organization is considered much less regulated than both the bacterial nucleoid and the nuclear chromatin. However, our recently identified conserved DNase footprinting pattern in human mtDNA, which co-localizes with regulatory elements and responds to physiological conditions, likely reflects a structured higher-order mtDNA organization. We hypothesized that this pattern emerges during embryogenesis. To test this hypothesis, we analyzed assay for transposase-accessible chromatin sequencing (ATAC-seq) results collected during the course of mouse and human early embryogenesis. Our results reveal, for the first time, a gradual and dynamic emergence of the adult mtDNA footprinting pattern during embryogenesis of both mammals. Taken together, our findings suggest that the structured adult chromatin-like mtDNA organization is gradually formed during mammalian embryogenesis., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

The bacterial heritage of mitochondria, as well as its independent genome [mitochondrial DNA (mtDNA)] and polycistronic transcripts, led to the view that mitochondrial transcriptional regulation relies on an evolutionarily conserved, prokaryotic-like system that is separated from the rest of the cell. Indeed, mtDNA transcription was previously thought to be governed by a few dedicated direct regulators, namely, the mitochondrial RNA polymerase (POLRMT), two transcription factors (TFAM and TF2BM), one transcription elongation (TEFM), and one known transcription termination factor (mTERF1). Recent findings have, however, revealed that known nuclear gene expression regulators are also involved in mtDNA transcription and have identified novel transcriptional features consistent with adaptation of the mitochondria to the regulatory environment of the precursor of the eukaryotic cell. Finally, whereas mammals follow the human mtDNA transcription pattern, other organisms notably diverge in terms of mtDNA transcriptional regulation. Hence, mtDNA transcriptional regulation is likely more evolutionary diverse than once thought., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Mitochondrial dysfunction has repeatedly been reported associated with type 2 diabetes mellitus (T2DM) and metabolic syndrome (MS), as have mitochondrial DNA (mtDNA) tRNA and duplication mutations and mtDNA haplogroup lineages. We identified 19 Taiwanese T2DM and MS pedigrees from Taiwan, with putative matrilineal transmission, one of which harbored the pathogenic mtDNA tRNA Leu(UUR) nucleotide (nt) 3243A>G mutation on the N9a3 haplogroup background. We then recruited three independent Taiwanese cohorts, two from Taipei (N = 498, mean age 52 and N = 1002, mean age 44) and one from a non-urban environment (N = 501, mean age 57). All three cohorts were assessed for an array of metabolic parameters, their mtDNA haplogroups determined, and the haplogroups correlated with T2DM/MS phenotypes. Logistic regression analysis revealed that mtDNA haplogroups D5, F4, and N9a conferred T2DM protection, while haplogroups F4 and N9a were risk factors for hypertension (HTN), and F4 was a risk factor for obesity (OB). Additionally, the 5263C>T (ND2 A165V) variant commonly associated with F4 was associated with hypertension (HTN). Cybrids were prepared with macro-haplogroup N (defined by variants m.ND3 10398A (114T) and m.ATP6 8701A (59T)) haplogroups B4 and F1 mtDNAs and from macro-haplogroup M (variants m.ND3 10398G (114A) and m.ATP6 8701G (59A)) haplogroup M9 mtDNAs. Additionally, haplogroup B4 and F1 cybrids were prepared with and without the mtDNA variant in ND1 3394T>C (Y30H) reported to be associated with T2DM. Assay of mitochondria complex I in these cybrids revealed that macro-haplogroup N cybrids had lower activity than M cybrids, that haplogroup F cybrids had lower activity than B4 cybrids, and that the ND1 3394T>C (Y30H) variant reduced complex I on both the B4 and F1 background but with very different cumulative effects. These data support the hypothesis that functional mtDNA variants may contribute to the risk of developing T2DM and MS., (Copyright © 2018. Published by Elsevier B.V.)

Human mitochondrial DNA (mtDNA) is believed to lack chromatin and histones. Instead, it is coated solely by the transcription factor TFAM. We asked whether mtDNA packaging is more regulated than once thought. To address this, we analyzed DNase-seq experiments in 324 human cell types and found, for the first time, a pattern of 29 mtDNA Genomic footprinting (mt-DGF) sites shared by ∼90% of the samples. Their syntenic conservation in mouse DNase-seq experiments reflect selective constraints. Colocalization with known mtDNA regulatory elements, with G-quadruplex structures, in TFAM-poor sites (in HeLa cells) and with transcription pausing sites, suggest a functional regulatory role for such mt-DGFs. Altered mt-DGF pattern in interleukin 3-treated CD34 + cells, certain tissue differences, and significant prevalence change in fetal versus nonfetal samples, offer first clues to their physiological importance. Taken together, human mtDNA has a conserved protein-DNA organization, which is likely involved in mtDNA regulation., (© 2018 Blumberg et al.; Published by Cold Spring Harbor Laboratory Press.)

Oxidative phosphorylation (OXPHOS), a fundamental energy source in all human tissues, requires interactions between mitochondrial (mtDNA)- and nuclear (nDNA)-encoded protein subunits. Although such interactions are fundamental to OXPHOS, bi-genomic coregulation is poorly understood. To address this question, we analyzed ∼8500 RNA-seq experiments from 48 human body sites. Despite well-known variation in mitochondrial activity, quantity, and morphology, we found overall positive mtDNA-nDNA OXPHOS genes' co-expression across human tissues. Nevertheless, negative mtDNA-nDNA gene expression correlation was identified in the hypothalamus, basal ganglia, and amygdala (subcortical brain regions, collectively termed the "primitive" brain). Single-cell RNA-seq analysis of mouse and human brains revealed that this phenomenon is evolutionarily conserved, and both are influenced by brain cell types (involving excitatory/inhibitory neurons and nonneuronal cells) and by their spatial brain location. As the "primitive" brain is highly oxidative, we hypothesized that such negative mtDNA-nDNA co-expression likely controls for the high mtDNA transcript levels, which enforce tight OXPHOS regulation, rather than rewiring toward glycolysis. Accordingly, we found "primitive" brain-specific up-regulation of lactate dehydrogenase B ( LDHB ), which associates with high OXPHOS activity, at the expense of LDHA , which promotes glycolysis. Analyses of co-expression, DNase-seq, and ChIP-seq experiments revealed candidate RNA-binding proteins and CEBPB as the best regulatory candidates to explain these phenomena. Finally, cross-tissue expression analysis unearthed tissue-dependent splice variants and OXPHOS subunit paralogs and allowed revising the list of canonical OXPHOS transcripts. Taken together, our analysis provides a comprehensive view of mito-nuclear gene co-expression across human tissues and provides overall insights into the bi-genomic regulation of mitochondrial activities., (© 2018 Barshad et al.; Published by Cold Spring Harbor Laboratory Press.)

The central role of the mitochondrion for cellular and organismal metabolism is well known, yet its functional role in evolution has rarely been featured in leading international conferences. Moreover, the contribution of mitochondrial genetics to complex disease phenotypes is particularly important, and although major advances have been made in the field of genomics, mitochondrial genomic data have in many cases been overlooked. Accumulating data and new knowledge support a major contribution of this maternally inherited genome, and its interactions with the nucleus, to both major evolutionary processes and diverse disease phenotypes. These advances encouraged us to assemble the first Mitochondrial Genomics and Evolution (MGE) meeting-an SMBE satellite and Israeli Science foundation international conference (Israel, September 2017). Here, we report the content and outcome of the MGE meeting (https://www.mge2017.com/; last accessed November 5, 2017)., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Mitochondrial DNA (mtDNA) genes are long known to be cotranscribed in polycistrones, yet it remains impossible to study nascent mtDNA transcripts quantitatively in vivo using existing tools. To this end, we used deep sequencing (GRO-seq and PRO-seq) and analyzed nascent mtDNA-encoded RNA transcripts in diverse human cell lines and metazoan organisms. Surprisingly, accurate detection of human mtDNA transcription initiation sites (TISs) in the heavy and light strands revealed a novel conserved transcription pausing site near the light-strand TIS. This pausing site correlated with the presence of a bacterial pausing sequence motif, with reduced SNP density, and with a DNase footprinting signal in all tested cells. Its location within conserved sequence block 3 (CSBIII), just upstream of the known transcription-replication transition point, suggests involvement in such transition. Analysis of nonhuman organisms enabled de novo mtDNA sequence assembly, as well as detection of previously unknown mtDNA TIS, pausing, and transcription termination sites with unprecedented accuracy. Whereas mammals ( Pan troglodytes , Macaca mulatta , Rattus norvegicus , and Mus musculus ) showed a human-like mtDNA transcription pattern, the invertebrate pattern ( Drosophila melanogaster and Caenorhabditis elegans ) profoundly diverged. Our approach paves the path toward in vivo, quantitative, reference sequence-free analysis of mtDNA transcription in all eukaryotes., (© 2017 Blumberg et al.; Published by Cold Spring Harbor Laboratory Press.)

Human mtDNA genetic variants have traditionally been considered markers for ancient population migrations. However, during the past three decades, these variants have been associated with altered susceptibility to various phenotypes, thus supporting their importance for human health. Nevertheless, mtDNA disease association has frequently been supported only in certain populations, due either to population stratification or differential epistatic compensations among populations. To partially overcome these obstacles, we performed meta-analysis of the multiple mtDNA association studies conducted until 2016, encompassing 53,975 patients and 63,323 controls. Our findings support the association of mtDNA haplogroups and recurrent variants with specific phenotypes such as Parkinson's disease, type 2 diabetes, longevity, and breast cancer. Strikingly, our assessment of mtDNA variants' involvement with multiple phenotypes revealed significant impact for Caucasian haplogroups H, J, and K. Therefore, ancient mtDNA variants could be divided into those that affect specific phenotypes, versus others with a general impact on phenotype combinations. We suggest that the mtDNA could serve as a model for phenotype specificity versus allele heterogeneity.

Many lineage-defining (nodal) mutations possess high functionality. However, differentiating adaptive nodal mutations from those that are functionally compensated remains challenging. To address this challenge, we identified functional nodal mutations (fNMs) in ~3,400 nuclear DNA (nDNA) and 4 mitochondrial DNA (mtDNA) protein structures from 91 and 1,003 species, respectively, representing the entire mammalian, bird and reptile phylogeny. A screen for candidate compensatory mutations among co-occurring amino acid changes in close structural proximity revealed that such compensated fNMs encompass 37% and 27% of the mtDNA and nDNA datasets, respectively. Analysis of the remaining (non-compensated) mutations, which are enriched for adaptive mutations, showed that birds and mammals share most such recurrent fNMs (N = 51). Among the latter, we discovered mutations in thermoregulation-related genes. These represent the best candidates to explain the molecular basis of convergent body thermoregulation in birds and mammals. Our analysis reveals the landscape of possible mutational compensation and convergence in amniote phylogeny.

Mitochondrial DNA (mtDNA) variants have been traditionally used as markers to trace ancient population migrations. Although experiments relying on model organisms and cytoplasmic hybrids, as well as disease association studies, have served to underline the functionality of certain mtDNA SNPs, only little is known of the regulatory impact of ancient mtDNA variants, especially in terms of gene expression. By analyzing RNA-seq data of 454 lymphoblast cell lines from the 1000 Genomes Project, we found that mtDNA variants defining the most common African genetic background, the L haplogroup, exhibit a distinct overall mtDNA gene expression pattern, which was independent of mtDNA copy numbers. Secondly, intra-population analysis revealed subtle, yet significant, expression differences in four tRNA genes. Strikingly, the more prominent African mtDNA gene expression pattern best correlated with the expression of nuclear DNA-encoded RNA-binding proteins, and with SNPs within the mitochondrial RNA-binding proteins PTCD1 and MRPS7. Our results thus support the concept of an ancient regulatory transition of mtDNA-encoded genes as humans left Africa to populate the rest of the world., Competing Interests: The authors have declared that no competing interests exist.

The mitochondrial ribosome, which translates all mitochondrial DNA (mtDNA)-encoded proteins, should be tightly regulated pre- and post-transcriptionally. Recently, we found RNA-DNA differences (RDDs) at human mitochondrial 16S (large) rRNA position 947 that were indicative of post-transcriptional modification. Here, we show that these 16S rRNA RDDs result from a 1-methyladenosine (m1A) modification introduced by TRMT61B, thus being the first vertebrate methyltransferase that modifies both tRNA and rRNAs. m1A947 is conserved in humans and all vertebrates having adenine at the corresponding mtDNA position (90% of vertebrates). However, this mtDNA base is a thymine in 10% of the vertebrates and a guanine in the 23S rRNA of 95% of bacteria, suggesting alternative evolutionary solutions. m1A, uridine, or guanine may stabilize the local structure of mitochondrial and bacterial ribosomes. Experimental assessment of genome-edited Escherichia coli showed that unmodified adenine caused impaired protein synthesis and growth. Our findings revealed a conserved mechanism of rRNA modification that has been selected instead of DNA mutations to enable proper mitochondrial ribosome function., Competing Interests: The authors have declared that no competing interests exist.

RNA-seq is becoming a preferred tool for genomics studies of model and non-model organisms. However, DNA-based analysis of organisms lacking sequenced genomes cannot rely on RNA-seq data alone to isolate most genes of interest, as DNA codes both exons and introns. With this in mind, we designed a novel tool, LEMONS, that exploits the evolutionary conservation of both exon/intron boundary positions and splice junction recognition signals to produce high throughput splice-junction predictions in the absence of a reference genome. When tested on multiple annotated vertebrate mRNA data, LEMONS accurately identified 87% (average) of the splice-junctions. LEMONS was then applied to our updated Mediterranean chameleon transcriptome, which lacks a reference genome, and predicted a total of 90,820 exon-exon junctions. We experimentally verified these splice-junction predictions by amplifying and sequencing twenty randomly selected genes from chameleon DNA templates. Exons and introns were detected in 19 of 20 of the positions predicted by LEMONS. To the best of our knowledge, LEMONS is currently the only experimentally verified tool that can accurately predict splice-junctions in organisms that lack a reference genome.

Compatibility between the nuclear (nDNA) and mitochondrial (mtDNA) genomes is important for organismal health. However, its significance for major evolutionary processes such as speciation is unclear, especially in vertebrates. We previously identified a sharp mtDNA-specific sequence divergence between morphologically indistinguishable chameleon populations (Chamaeleo chamaeleon recticrista) across an ancient Israeli marine barrier (Jezreel Valley). Because mtDNA introgression and gender-based dispersal were ruled out, we hypothesized that mtDNA spatial division was maintained by mito-nuclear functional compensation. Here, we studied RNA-seq generated from each of ten chameleons representing the north and south populations and identified candidate nonsynonymous substitutions (NSSs) matching the mtDNA spatial distribution. The most prominent NSS occurred in 14 nDNA-encoded mitochondrial proteins. Increased chameleon sample size (N = 70) confirmed the geographic differentiation in POLRMT, NDUFA5, ACO1, LYRM4, MARS2, and ACAD9. Structural and functionality evaluation of these NSSs revealed high functionality. Mathematical modeling suggested that this mito-nuclear spatial divergence is consistent with hybrid breakdown. We conclude that our presented evidence and mathematical model underline mito-nuclear interactions as a likely role player in incipient speciation in vertebrates., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Parkin, which is mutated in most recessive Parkinsonism, is a key player in the selective removal of damaged mitochondria via mitophagy. Damaged mitochondria may carry mitochondrial DNA (mtDNA) mutations, thus creating a mixed mtDNA population within cells (heteroplasmy). It was previously shown that Parkin over-expression reduced the level of heteroplasmic mutations that alter mitochondrial membrane potential in human cytoplasmic hybrids. However, it remained unclear whether Parkin serves a similar role at the entire living organism, and whether this role is evolutionarily conserved. Here, we show that mutation in the Caenorhabditis elegans orthologue of Parkin (pdr-1) modulates the level of a large heteroplasmic mtDNA truncation. Massive parallel sequencing revealed that the mtDNAs of C. elegans wild type and pdr-1(gk448) mutant strains were virtually deprived of heteroplasmy, thus reflecting strong negative selection against dysfunctional mitochondria. Therefore, our findings show that the role of Parkin in the modulation of heteroplasmy is conserved between human and worm and raise the interesting possibility that mitophagy modulates the striking lack of heteroplasmy in C. elegans., (Copyright © 2014 © Elsevier B.V. and Mitochondria Research Society. Published by Elsevier B.V. All rights reserved.)

In contrast to the nuclear genome, the mitochondrial DNA (mtDNA) is maternally inherited and resides in multiple cellular copies that may vary in sequence (heteroplasmy). Although the interaction between mtDNA and nuclear DNA-encoded factors (mito-nuclear interaction) is vital, the mtDNA accumulates mutations an order of magnitude faster than the nuclear genome both during evolution and during the lifetime of the individual, thus requiring tight mito-nuclear co-evolution. These unique features drew the attention of many to suggest a role for the mitochondria in ageing. Although an excess of mtDNA mutations has been found in aged humans and animal models, most of these mutations had minor functional potential. Moreover, there are mtDNA mutations that recur in aged humans, but do not have any clear functionality. Nevertheless, accumulation of recurrent private mutations with minor functionality in the fast-ageing, mtDNA polymerase mutated mice (Pol-gamma), suggested that these very mtDNA alterations participate in ageing. This introduces a paradox: how would either single or recurrent mutations with negligible functionality play a role in a major chronic phenotype such as ageing?Here, we propose a hypothesis to partially resolve this paradox: accumulation of mitochondrial mutations with subtle functionality, which was overlooked by the mechanisms of selection, supplemented by slightly affected fusion-fission cycles, will hamper mitochondrial functional complementation within cells, disrupt mito-nuclear interactions and lead to ageing. Since certain mito-nuclear genotypes are less functionally compatible than others, and since the mtDNA and the nuclear genome segregate independently among generations, mild functionality of mutations will have differential effect on individuals in the population thus explaining the large variability in the ageing phenotype even within ethnic groups. We emphasize the role of recurrent mtDNA mutations with functional potential during evolution and during the lifetime of the individual.

Most cell functions are carried out by interacting factors, thus underlying the functional importance of genetic interactions between genes, termed epistasis. Epistasis could be under strong selective pressures especially in conditions where the mutation rate of one of the interacting partners notably differs from the other. Accordingly, the order of magnitude higher mitochondrial DNA (mtDNA) mutation rate as compared to the nuclear DNA (nDNA) of all tested animals, should influence systems involving mitochondrial-nuclear (mito-nuclear) interactions. Such is the case of the energy producing oxidative phosphorylation (OXPHOS) and mitochondrial translational machineries which are comprised of factors encoded by both the mtDNA and the nDNA. Additionally, the mitochondrial RNA transcription and mtDNA replication systems are operated by nDNA-encoded proteins that bind mtDNA regulatory elements. As these systems are central to cell life there is strong selection toward mito-nuclear co-evolution to maintain their function. However, it is unclear whether (A) mito-nuclear co-evolution befalls only to retain mitochondrial functions during evolution or, also, (B) serves as an adaptive tool to adjust for the evolving energetic demands as species' complexity increases. As the first step to answer these questions we discuss evidence of both negative and adaptive (positive) selection acting on the mtDNA and nDNA-encoded genes and the effect of both types of selection on mito-nuclear interacting factors. Emphasis is given to the crucial role of recurrent ancient (nodal) mutations in such selective events. We apply this point-of-view to the three available types of mito-nuclear co-evolution: protein-protein (within the OXPHOS system), protein-RNA (mainly within the mitochondrial ribosome), and protein-DNA (at the mitochondrial replication and transcription machineries).

Transcription of mitochondrial DNA (mtDNA)-encoded genes is thought to be regulated by a handful of dedicated transcription factors (TFs), suggesting that mtDNA genes are separately regulated from the nucleus. However, several TFs, with known nuclear activities, were found to bind mtDNA and regulate mitochondrial transcription. Additionally, mtDNA transcriptional regulatory elements, which were proved important in vitro, were harbored by a deletion that normally segregated among healthy individuals. Hence, mtDNA transcriptional regulation is more complex than once thought. Here, by analyzing ENCODE chromatin immunoprecipitation sequencing (ChIP-seq) data, we identified strong binding sites of three bona fide nuclear TFs (c-Jun, Jun-D, and CEBPb) within human mtDNA protein-coding genes. We validated the binding of two TFs by ChIP-quantitative polymerase chain reaction (c-Jun and Jun-D) and showed their mitochondrial localization by electron microscopy and subcellular fractionation. As a step toward investigating the functionality of these TF-binding sites (TFBS), we assessed signatures of selection. By analyzing 9,868 human mtDNA sequences encompassing all major global populations, we recorded genetic variants in tips and nodes of mtDNA phylogeny within the TFBS. We next calculated the effects of variants on binding motif prediction scores. Finally, the mtDNA variation pattern in predicted TFBS, occurring within ChIP-seq negative-binding sites, was compared with ChIP-seq positive-TFBS (CPR). Motifs within CPRs of c-Jun, Jun-D, and CEBPb harbored either only tip variants or their nodal variants retained high motif prediction scores. This reflects negative selection within mtDNA CPRs, thus supporting their functionality. Hence, human mtDNA-coding sequences may have dual roles, namely coding for genes yet possibly also possessing regulatory potential., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The mutation rate of the mitochondrial DNA (mtDNA), which is higher by an order of magnitude as compared with the nuclear genome, enforces tight mitonuclear coevolution to maintain mitochondrial activities. Interruption of such coevolution plays a role in interpopulation hybrid breakdown, speciation events, and disease susceptibility. Previously, we found an elevated amino acid replacement rate and positive selection in the nuclear DNA-encoded oxidative phosphorylation (OXPHOS) complex I subunit NDUFC2, a phenomenon important for the direct interaction of NDUFC2 with the mtDNA-encoded complex I subunit ND4. This finding underlines the importance of mitonuclear coevolution to physical interactions between mtDNA and nuclear DNA-encoded factors. Nevertheless, it remains unclear whether this interaction is important for the stability and activity of complex I. Here, we show that siRNA silencing of NDUFC2 reduced growth of human D-407 retinal pigment epithelial cells, significantly diminished mitochondrial membrane potential, and interfered with complex I integrity. Moreover, site-directed mutagenesis of a positively selected amino acid in NDUFC2 significantly interfered with the interaction of NDUFC2 with its mtDNA-encoded partner ND4. Finally, we show that a genotype combination involving this amino acid (NDUFC2 residue 46) and the mtDNA haplogroup HV likely altered susceptibility to type 2 diabetes mellitus in Ashkenazi Jews. Therefore, mitonuclear coevolution is important for maintaining mitonuclear factor interactions, OXPHOS, and for human health., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

RNA transcripts are generally identical to the underlying DNA sequences. Nevertheless, RNA-DNA differences (RDDs) were found in the nuclear human genome and in plants and animals but not in human mitochondria. Here, by deep sequencing of human mitochondrial DNA (mtDNA) and RNA, we identified three RDD sites at mtDNA positions 295 (C-to-U), 13710 (A-to-U, A-to-G), and 2617 (A-to-U, A-to-G). Position 2617, within the 16S rRNA, harbored the most prevalent RDDs (>30% A-to-U and ∼15% A-to-G of the reads in all tested samples). The 2617 RDDs appeared already at the precursor polycistrone mitochondrial transcript. By using traditional Sanger sequencing, we identified the A-to-U RDD in six different cell lines and representative primates (Gorilla gorilla, Pongo pigmaeus, and Macaca mulatta), suggesting conservation of the mechanism generating such RDD. Phylogenetic analysis of more than 1700 vertebrate mtDNA sequences supported a thymine as the primate ancestral allele at position 2617, suggesting that the 2617 RDD recapitulates the ancestral 16S rRNA. Modeling U or G (the RDDs) at position 2617 stabilized the large ribosomal subunit structure in contrast to destabilization by an A (the pre-RDDs). Hence, these mitochondrial RDDs are likely functional.

Mutations frequently reoccur in the human mitochondrial DNA (mtDNA). However, it is unclear whether recurrent mtDNA nodal mutations (RNMs), that is, recurrent mutations in stems of unrelated phylogenetic nodes, are functional and hence selectively constrained. To answer this question, we performed comprehensive parsimony and maximum likelihood analyses of 9,868 publicly available whole human mtDNAs revealing 1,606 single nodal mutations (SNMs) and 679 RNMs. We then evaluated the potential functionality of synonymous, nonsynonymous and RNA SNMs and RNMs. For synonymous mutations, we have implemented the Codon Adaptation Index. For nonsynonymous mutations, we assessed evolutionary conservation, and employed previously described pathogenicity score assessment tools. For RNA genes' mutations, we designed a bioinformatic tool which compiled evolutionary conservation and potential effect on RNA structure. While comparing the functionality scores of nonsynonymous and RNA SNMs and RNMs with those of disease-causing mtDNA mutations, we found significant difference (P < 0.001). However, 24 RNMs and 67 SNMs had comparable values with disease-causing mutations reflecting their potential function thus being the best candidates to participate in adaptive events of unrelated lineages. Strikingly, some functional RNMs occurred in unrelated mtDNA lineages that independently altered susceptibility to the same diseases, thus suggesting common functionality. To our knowledge, this is the most comprehensive analysis of selective signatures in the mtDNA not only within proteins but also within RNA genes. For the first time, we discover virtually all positively selected RNMs in our phylogeny while emphasizing their dual role in past evolutionary events and in disease today.

Heteroplasmy, the mixture of mitochondrial genomes (mtDNA), varies among individuals and cells. Heteroplasmy levels alter the penetrance of pathological mtDNA mutations, and the susceptibility to age-related diseases such as Parkinson's disease. Although mitochondrial dysfunction occurs in age-related type 2 diabetes mellitus (T2DM), the involvement of heteroplasmy in diabetes is unclear. We hypothesized that the heteroplasmic mutational (HM) pattern may change in T2DM. To test this, we used next-generation sequencing, i.e. massive parallel sequencing (MPS), along with PCR-cloning-Sanger sequencing to analyze HM in blood and skeletal muscle DNA samples from monozygotic (MZ) twins either concordant or discordant for T2DM. Great variability was identified in the repertoires and amounts of HMs among individuals, with a tendency towards more mutations in skeletal muscle than in blood. Whereas many HMs were unique, many were either shared among twin pairs or among tissues of the same individual, regardless of their prevalence. This suggested a heritable influence on even low abundance HMs. We found no clear differences between T2DM and controls. However, we found ~5-fold increase of HMs in non-coding sequences implying the influence of negative selection (P < 0.001). This negative selection was evident both in moderate to highly abundant heteroplasmy (>5% of the molecules per sample) and in low abundance heteroplasmy (<5% of the molecules). Although our study found no evidence supporting the involvement of HMs in the etiology of T2DM, the twin study found clear evidence of a heritable influence on the accumulation of HMs as well as the signatures of selection in heteroplasmic mutations.

Factors required for mitochondrial function are encoded both by the nuclear and mitochondrial genomes. The order of magnitude higher mutation rate of animal mitochondrial DNA (mtDNA) enforces tight co-evolution of mtDNA and nuclear DNA encoded factors. In this essay we argue that such co evolution exists at the population and inter-specific levels and affect disease susceptibility. We also argue for the existence of three modes of co-evolution in the mitochondrial genetic system, which include the interaction of mtDNA and nuclear DNA encoded proteins, nuclear protein - mtDNA-encoded RNA interaction within the mitochondrial translation machinery and nuclear DNA encoded proteins-mtDNA binging sites interaction in the frame of the mtDNA replication and transcription machineries. These modes of co evolution require co-regulation of the interacting factors encoded by the two genomes. Thus co evolution plays an important role in modulating mitochondrial activity. This article is part of a Special Issue entitled: Mitochondrial Gene Expression., (Copyright © 2011 Elsevier B.V. All rights reserved.)

The Levant is an important migration bridge, harboring border-zones between Afrotropical and palearctic species. Accordingly, Chameleo chameleon, a common species throughout the Mediterranean basin, is morphologically divided in the southern Levant (Israel) into two subspecies, Chamaeleo chamaeleon recticrista (CCR) and C. c. musae (CCM). CCR mostly inhabits the Mediterranean climate (northern Israel), while CCM inhabits the sands of the north-western Negev Desert (southern Israel). AFLP analysis of 94 geographically well dispersed specimens indicated moderate genetic differentiation (PhiPT = 0.097), consistent with the classical division into the two subspecies, CCR and CCM. In contrast, sequence analysis of a 637 bp coding mitochondrial DNA (mtDNA) fragment revealed two distinct phylogenetic clusters which were not consistent with the morphological division: one mtDNA cluster consisted of CCR specimens collected in regions northern of the Jezreel Valley and another mtDNA cluster harboring specimens pertaining to both the CCR and CCM subspecies but collected southern of the Jezreel Valley. AMOVA indicated clear mtDNA differentiation between specimens collected northern and southern to the Jezreel Valley (PhiPT = 0.79), which was further supported by a very low coalescent-based estimate of effective migration rates. Whole chameleon mtDNA sequencing (∼17,400 bp) generated from 11 well dispersed geographic locations revealed 325 mutations sharply differentiating the two mtDNA clusters, suggesting a long allopatric history further supported by BEAST. This separation correlated temporally with the existence of an at least 1 million year old marine barrier at the Jezreel Valley exactly where the mtDNA clusters meet. We discuss possible involvement of gender-dependent life history differences in maintaining such mtDNA genetic differentiation and suggest that it reflects (ancient) local adaptation to mitochondrial-related traits.