Single-cell analysis reveals context-dependent, cell-level selection of mtDNA.

Heteroplasmy occurs when wild-type and mutant mitochondrial DNA (mtDNA) molecules co-exist in single cells1. Heteroplasmy levels change dynamically in development, disease and ageing2,3, but it is unclear whether these shifts are caused by selection or drift, and whether they occur at the level of cells or intracellularly. Here we investigate heteroplasmy dynamics in dividing cells by combining precise mtDNA base editing (DdCBE)4 with a new method, SCI-LITE (single-cell combinatorial indexing leveraged to interrogate targeted expression), which tracks single-cell heteroplasmy with ultra-high throughput. We engineered cells to have synonymous or nonsynonymous complex I mtDNA mutations and found that cell populations in standard culture conditions purge nonsynonymous mtDNA variants, whereas synonymous variants are maintained. This suggests that selection dominates over simple drift in shaping population heteroplasmy. We simultaneously tracked single-cell mtDNA heteroplasmy and ancestry, and found that, although the population heteroplasmy shifts, the heteroplasmy of individual cell lineages remains stable, arguing that selection acts at the level of cell fitness in dividing cells. Using these insights, we show that we can force cells to accumulate high levels of truncating complex I mtDNA heteroplasmy by placing them in environments where loss of biochemical complex I activity has been reported to benefit cell fitness. We conclude that in dividing cells, a given nonsynonymous mtDNA heteroplasmy can be harmful, neutral or even beneficial to cell fitness, but that the 'sign' of the effect is wholly dependent on the environment.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans.

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

CRISPR-free base editors with enhanced activity and expanded targeting scope in mitochondrial and nuclear DNA.

The all-protein cytosine base editor DdCBE uses TALE proteins and a double-stranded DNA-specific cytidine deaminase (DddA) to mediate targeted C•G-to-T•A editing. To improve editing efficiency and overcome the strict TC sequence-context constraint of DddA, we used phage-assisted non-continuous and continuous evolution to evolve DddA variants with improved activity and expanded targeting scope. Compared to canonical DdCBEs, base editors with evolved DddA6 improved mitochondrial DNA (mtDNA) editing efficiencies at TC by 3.3-fold on average. DdCBEs containing evolved DddA11 offered a broadened HC (H = A, C or T) sequence compatibility for both mitochondrial and nuclear base editing, increasing average editing efficiencies at AC and CC targets from less than 10% for canonical DdCBE to 15-30% and up to 50% in cell populations sorted to express both halves of DdCBE. We used these evolved DdCBEs to efficiently install disease-associated mtDNA mutations in human cells at non-TC target sites. DddA6 and DddA11 substantially increase the effectiveness and applicability of all-protein base editing.

High-Throughput Measurement of Mitochondrial RNA Turnover in Human Cultured Cells.

RNA turnover is an essential part of the gene expression pathway, and there are several experimental approaches for its determination. High-throughput measurement of global RNA turnover rates can provide valuable information about conditions or proteins that impact gene expression. Here, we present a protocol for mitochondrial RNA turnover analysis which involves metabolic labeling of RNA coupled with quantitative high-throughput fluorescent microscopy. This approach gives an excellent opportunity to discover new factors involved in mitochondrial gene regulation when combined with loss-of-function screening strategy.

MitoCarta3.0: an updated mitochondrial proteome now with sub-organelle localization and pathway annotations.

The mammalian mitochondrial proteome is under dual genomic control, with 99% of proteins encoded by the nuclear genome and 13 originating from the mitochondrial DNA (mtDNA). We previously developed MitoCarta, a catalogue of over 1000 genes encoding the mammalian mitochondrial proteome. This catalogue was compiled using a Bayesian integration of multiple sequence features and experimental datasets, notably protein mass spectrometry of mitochondria isolated from fourteen murine tissues. Here, we introduce MitoCarta3.0. Beginning with the MitoCarta2.0 inventory, we performed manual review to remove 100 genes and introduce 78 additional genes, arriving at an updated inventory of 1136 human genes. We now include manually curated annotations of sub-mitochondrial localization (matrix, inner membrane, intermembrane space, outer membrane) as well as assignment to 149 hierarchical 'MitoPathways' spanning seven broad functional categories relevant to mitochondria. MitoCarta3.0, including sub-mitochondrial localization and MitoPathway annotations, is freely available at http://www.broadinstitute.org/mitocarta and should serve as a continued community resource for mitochondrial biology and medicine.

A bacterial cytidine deaminase toxin enables CRISPR-free mitochondrial base editing.

Bacterial toxins represent a vast reservoir of biochemical diversity that can be repurposed for biomedical applications. Such proteins include a group of predicted interbacterial toxins of the deaminase superfamily, members of which have found application in gene-editing techniques1,2. Because previously described cytidine deaminases operate on single-stranded nucleic acids3, their use in base editing requires the unwinding of double-stranded DNA (dsDNA)-for example by a CRISPR-Cas9 system. Base editing within mitochondrial DNA (mtDNA), however, has thus far been hindered by challenges associated with the delivery of guide RNA into the mitochondria4. As a consequence, manipulation of mtDNA to date has been limited to the targeted destruction of the mitochondrial genome by designer nucleases9,10.Here we describe an interbacterial toxin, which we name DddA, that catalyses the deamination of cytidines within dsDNA. We engineered split-DddA halves that are non-toxic and inactive until brought together on target DNA by adjacently bound programmable DNA-binding proteins. Fusions of the split-DddA halves, transcription activator-like effector array proteins, and a uracil glycosylase inhibitor resulted in RNA-free DddA-derived cytosine base editors (DdCBEs) that catalyse C•G-to-T•A conversions in human mtDNA with high target specificity and product purity. We used DdCBEs to model a disease-associated mtDNA mutation in human cells, resulting in changes in respiration rates and oxidative phosphorylation. CRISPR-free DdCBEs enable the precise manipulation of mtDNA, rather than the elimination of mtDNA copies that results from its cleavage by targeted nucleases, with broad implications for the study and potential treatment of mitochondrial disorders.

Human REXO2 controls short mitochondrial RNAs generated by mtRNA processing and decay machinery to prevent accumulation of double-stranded RNA.

RNA decay is a key element of mitochondrial RNA metabolism. To date, the only well-documented machinery that plays a role in mtRNA decay in humans is the complex of polynucleotide phosphorylase (PNPase) and SUV3 helicase, forming the degradosome. REXO2, a homolog of prokaryotic oligoribonucleases present in humans both in mitochondria and the cytoplasm, was earlier shown to be crucial for maintaining mitochondrial homeostasis, but its function in mitochondria has not been fully elucidated. In the present study, we created a cellular model that enables the clear dissection of mitochondrial and non-mitochondrial functions of human REXO2. We identified a novel mitochondrial short RNA, referred to as ncH2, that massively accumulated upon REXO2 silencing. ncH2 degradation occurred independently of the mitochondrial degradosome, strongly supporting the hypothesis that ncH2 is a primary substrate of REXO2. We also investigated the global impact of REXO2 depletion on mtRNA, revealing the importance of the protein for maintaining low steady-state levels of mitochondrial antisense transcripts and double-stranded RNA. Our detailed biochemical and structural studies provide evidence of sequence specificity of the REXO2 oligoribonuclease. We postulate that REXO2 plays dual roles in human mitochondria, 'scavenging' nanoRNAs that are produced by the degradosome and clearing short RNAs that are generated by RNA processing.

Mitochondrial Gene Expression and Beyond-Novel Aspects of Cellular Physiology.

Mitochondria are peculiar organelles whose proper function depends on the crosstalk between two genomes, mitochondrial and nuclear. The human mitochondrial genome (mtDNA) encodes only 13 proteins; nevertheless, its proper expression is essential for cellular homeostasis, as mtDNA-encoded proteins are constituents of mitochondrial respiratory complexes. In addition, mtDNA expression results in the production of RNA molecules, which influence cell physiology once released from the mitochondria into the cytoplasm. As a result, dysfunctions of mtDNA expression may lead to pathologies in humans. Here, we review the mechanisms of mitochondrial gene expression with a focus on recent findings in the field. We summarize the complex turnover of mitochondrial transcripts and present an increasing body of evidence indicating new functions of mitochondrial transcripts. We discuss mitochondrial gene regulation in different cellular contexts, focusing on stress conditions. Finally, we highlight the importance of emerging aspects of mitochondrial gene regulation in human health and disease.

Quantitative proteomics revealed C6orf203/MTRES1 as a factor preventing stress-induced transcription deficiency in human mitochondria.

Maintenance of mitochondrial gene expression is crucial for cellular homeostasis. Stress conditions may lead to a temporary reduction of mitochondrial genome copy number, raising the risk of insufficient expression of mitochondrial encoded genes. Little is known how compensatory mechanisms operate to maintain proper mitochondrial transcripts levels upon disturbed transcription and which proteins are involved in them. Here we performed a quantitative proteomic screen to search for proteins that sustain expression of mtDNA under stress conditions. Analysis of stress-induced changes of the human mitochondrial proteome led to the identification of several proteins with poorly defined functions among which we focused on C6orf203, which we named MTRES1 (Mitochondrial Transcription Rescue Factor 1). We found that the level of MTRES1 is elevated in cells under stress and we show that this upregulation of MTRES1 prevents mitochondrial transcript loss under perturbed mitochondrial gene expression. This protective effect depends on the RNA binding activity of MTRES1. Functional analysis revealed that MTRES1 associates with mitochondrial RNA polymerase POLRMT and acts by increasing mitochondrial transcription, without changing the stability of mitochondrial RNAs. We propose that MTRES1 is an example of a protein that protects the cell from mitochondrial RNA loss during stress.

Versatile approach for functional analysis of human proteins and efficient stable cell line generation using FLP-mediated recombination system.

Deciphering a function of a given protein requires investigating various biological aspects. Usually, the protein of interest is expressed with a fusion tag that aids or allows subsequent analyses. Additionally, downregulation or inactivation of the studied gene enables functional studies. Development of the CRISPR/Cas9 methodology opened many possibilities but in many cases it is restricted to non-essential genes. Recombinase-dependent gene integration methods, like the Flp-In system, are very good alternatives. The system is widely used in different research areas, which calls for the existence of compatible vectors and efficient protocols that ensure straightforward DNA cloning and generation of stable cell lines. We have created and validated a robust series of 52 vectors for streamlined generation of stable mammalian cell lines using the FLP recombinase-based methodology. Using the sequence-independent DNA cloning method all constructs for a given coding-sequence can be made with just three universal PCR primers. Our collection allows tetracycline-inducible expression of proteins with various tags suitable for protein localization, FRET, bimolecular fluorescence complementation (BiFC), protein dynamics studies (FRAP), co-immunoprecipitation, the RNA tethering assay and cell sorting. Some of the vectors contain a bidirectional promoter for concomitant expression of miRNA and mRNA, so that a gene can be silenced and its product replaced by a mutated miRNA-insensitive version. Our toolkit and protocols have allowed us to create more than 500 constructs with ease. We demonstrate the efficacy of our vectors by creating stable cell lines with various tagged proteins (numatrin, fibrillarin, coilin, centrin, THOC5, PCNA). We have analysed transgene expression over time to provide a guideline for future experiments and compared the effectiveness of commonly used inducers for tetracycline-responsive promoters. As proof of concept we examined the role of the exoribonuclease XRN2 in transcription termination by RNAseq.

Effect of Everolimus on Heterogenous Renal Cancer Cells Populations Including Renal Cancer Stem Cells.

The aim of this study was to compare effect of everolimus on growth of different renal cell carcinoma (RCC) populations and develop experimental design to measure the early response of everolimus in clear cell RCC (ccRCC) cell lines including renal cancer stem cells. Effect of everolimus on RCC cell lines which include primary (786-0) and metastatic (ACHN) RCC cell lines as well as heterogenous populations of tumor cells of different histological RCC subtypes (clear cell RCC and papillary RCC) was measured when treated with everolimus in the range of 1-9 µM. Gene expression profiling using microarray was performed to determine the early response to everolimus in ccRCC cell lines after optimizing concentration of drug. Gene Set Enrichment Analysis (GSEA) was done which mainly focused on basic genes related to mTOR, hormonal and metabolic pathways. Everolimus acts on RCC cells in a dose-dependent manner. In all examined cell lines IC50 dose was possible to calculate after the third day of treatment. In ccRCC lines (parental and stem cell) everolimus changes expression of mTOR complexes elements and elements of related pathways when treated with optimized doses of drug. Characteristic expression profile for ccRCC cells at an early exposure time to everolimus is to elucidate. Wevarie include some basic observations derived from data analysis in the context of mechanism of action of drug with a view to better understand biology of renal cancer cells.

The Dynamics of mRNA Turnover Revealed by Single-Molecule Imaging in Single Cells.

RNA degradation plays a fundamental role in regulating gene expression. In order to characterize the spatiotemporal dynamics of RNA turnover in single cells, we developed a fluorescent biosensor based on dual-color, single-molecule RNA imaging that allows intact transcripts to be distinguished from stabilized degradation intermediates. Using this method, we measured mRNA decay in single cells and found that individual degradation events occur independently within the cytosol and are not enriched within processing bodies. We show that slicing of an mRNA targeted for endonucleolytic cleavage by the RNA-induced silencing complex can be observed in real time in living cells. This methodology provides a framework for investigating the entire life history of individual mRNAs from birth to death in single cells.