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.

Circulating N-lactoyl-amino acids and N-formyl-methionine reflect mitochondrial dysfunction and predict mortality in septic shock.

INTRODUCTION: Sepsis is a highly morbid condition characterized by multi-organ dysfunction resulting from dysregulated inflammation in response to acute infection. Mitochondrial dysfunction may contribute to sepsis pathogenesis, but quantifying mitochondrial dysfunction remains challenging. OBJECTIVE: To assess the extent to which circulating markers of mitochondrial dysfunction are increased in septic shock, and their relationship to severity and mortality. METHODS: We performed both full-scan and targeted (known markers of genetic mitochondrial disease) metabolomics on plasma to determine markers of mitochondrial dysfunction which distinguish subjects with septic shock (n = 42) from cardiogenic shock without infection (n = 19), bacteremia without sepsis (n = 18), and ambulatory controls (n = 19) - the latter three being conditions in which mitochondrial function, proxied by peripheral oxygen consumption, is presumed intact. RESULTS: Nine metabolites were significantly increased in septic shock compared to all three comparator groups. This list includes N-formyl-L-methionine (f-Met), a marker of dysregulated mitochondrial protein translation, and N-lactoyl-phenylalanine (lac-Phe), representative of the N-lactoyl-amino acids (lac-AAs), which are elevated in plasma of patients with monogenic mitochondrial disease. Compared to lactate, the clinical biomarker used to define septic shock, there was greater separation between survivors and non-survivors of septic shock for both f-Met and the lac-AAs measured within 24 h of ICU admission. Additionally, tryptophan was the one metabolite significantly decreased in septic shock compared to all other groups, while its breakdown product kynurenate was one of the 9 significantly increased. CONCLUSION: Future studies which validate the measurement of lac-AAs and f-Met in conjunction with lactate could define a sepsis subtype characterized by mitochondrial dysfunction.

Lipoylation is dependent on the ferredoxin FDX1 and dispensable under hypoxia in human cells.

Iron-sulfur clusters (ISC) are essential cofactors that participate in electron transfer, environmental sensing, and catalysis. Amongst the most ancient ISC-containing proteins are the ferredoxin (FDX) family of electron carriers. Humans have two FDXs- FDX1 and FDX2, both of which are localized to mitochondria, and the latter of which is itself important for ISC synthesis. We have previously shown that hypoxia can eliminate the requirement for some components of the ISC biosynthetic pathway, but FDXs were not included in that study. Here, we report that FDX1, but not FDX2, is dispensable under 1% O2 in cultured human cells. We find that FDX1 is essential for production of the lipoic acid cofactor, which is synthesized by the ISC-containing enzyme lipoyl synthase. While hypoxia can rescue the growth phenotype of either FDX1 or lipoyl synthase KO cells, lipoylation in these same cells is not rescued, arguing against an alternative biosynthetic route or salvage pathway for lipoate in hypoxia. Our work reveals the divergent roles of FDX1 and FDX2 in mitochondria, identifies a role for FDX1 in lipoate synthesis, and suggests that loss of lipoic acid can be tolerated under low oxygen tensions in cell culture.

Semi-deconvolution of bulk and single-cell RNA-seq data with application to metastatic progression in breast cancer.

MOTIVATION: Identifying cell types and their abundances and how these evolve during tumor progression is critical to understanding the mechanisms of metastasis and identifying predictors of metastatic potential that can guide the development of new diagnostics or therapeutics. Single-cell RNA sequencing (scRNA-seq) has been especially promising in resolving heterogeneity of expression programs at the single-cell level, but is not always feasible, e.g. for large cohort studies or longitudinal analysis of archived samples. In such cases, clonal subpopulations may still be inferred via genomic deconvolution, but deconvolution methods have limited ability to resolve fine clonal structure and may require reference cell type profiles that are missing or imprecise. Prior methods can eliminate the need for reference profiles but show unstable performance when few bulk samples are available. RESULTS: In this work, we develop a new method using reference scRNA-seq to interpret sample collections for which only bulk RNA-seq is available for some samples, e.g. clonally resolving archived primary tissues using scRNA-seq from metastases. By integrating such information in a Quadratic Programming framework, our method can recover more accurate cell types and corresponding cell type abundances in bulk samples. Application to a breast tumor bone metastases dataset confirms the power of scRNA-seq data to improve cell type inference and quantification in same-patient bulk samples. AVAILABILITY AND IMPLEMENTATION: Source code is available on Github at https://github.com/CMUSchwartzLab/RADs.

Head tremor in cervical dystonia: Quantifying severity with computer vision.

BACKGROUND: Head tremor (HT) is a common feature of cervical dystonia (CD), usually quantified by subjective observation. Technological developments offer alternatives for measuring HT severity that are objective and amenable to automation. OBJECTIVES: Our objectives were to develop CMOR (Computational Motor Objective Rater; a computer vision-based software system) to quantify oscillatory and directional aspects of HT from video recordings during a clinical examination and to test its convergent validity with clinical rating scales. METHODS: For 93 participants with isolated CD and HT enrolled by the Dystonia Coalition, we analyzed video recordings from an examination segment in which participants were instructed to let their head drift to its most comfortable dystonic position. We evaluated peak power, frequency, and directional dominance, and used Spearman's correlation to measure the agreement between CMOR and clinical ratings. RESULTS: Power averaged 0.90 (SD 1.80) deg2/Hz, and peak frequency 1.95 (SD 0.94) Hz. The dominant HT axis was pitch (antero/retrocollis) for 50%, roll (laterocollis) for 6%, and yaw (torticollis) for 44% of participants. One-sided t-tests showed substantial contributions from the secondary (t = 18.17, p < 0.0001) and tertiary (t = 12.89, p < 0.0001) HT axes. CMOR's HT severity measure positively correlated with the HT item on the Toronto Western Spasmodic Torticollis Rating Scale-2 (Spearman's rho = 0.54, p < 0.001). CONCLUSIONS: We demonstrate a new objective method to measure HT severity that requires only conventional video recordings, quantifies the complexities of HT in CD, and exhibits convergent validity with clinical severity ratings.