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.

Identification of LY3522348: A Highly Selective and Orally Efficacious Ketohexokinase Inhibitor.

The identification of clinical candidate LY3522348 (compound 23) is described. LY3522348 is a highly selective, oral dual inhibitor of human ketohexokinase isoforms C and A (hKHK-C, hKHK-A). Optimization began with highly efficient (S)-2-(2-methylazetidin-1-yl)-6-(1H-pyrazol-4-yl)-4-(trifluoromethyl)nicotinonitrile (3). Efforts focused on developing absorption, distribution, metabolism, potency, and in vitro safety profiles to support oral QD dosing in patients. Structure-based design leveraged vectors for substitution of the pyrazole ring, which provided an opportunity to interact with several different proximal amino acid residues in the protein. LY3522348 displayed a robust pharmacodynamic response in a mouse model of fructose metabolism and was advanced into clinical trials.

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.

Hypoxia extends lifespan and neurological function in a mouse model of aging.

There is widespread interest in identifying interventions that extend healthy lifespan. Chronic continuous hypoxia delays the onset of replicative senescence in cultured cells and extends lifespan in yeast, nematodes, and fruit flies. Here, we asked whether chronic continuous hypoxia is beneficial in mammalian aging. We utilized the Ercc1 Δ/- mouse model of accelerated aging given that these mice are born developmentally normal but exhibit anatomic, physiological, and biochemical features of aging across multiple organs. Importantly, they exhibit a shortened lifespan that is extended by dietary restriction, the most potent aging intervention across many organisms. We report that chronic continuous 11% oxygen commenced at 4 weeks of age extends lifespan by 50% and delays the onset of neurological debility in Ercc1 Δ/- mice. Chronic continuous hypoxia did not impact food intake and did not significantly affect markers of DNA damage or senescence, suggesting that hypoxia did not simply alleviate the proximal effects of the Ercc1 mutation, but rather acted downstream via unknown mechanisms. To the best of our knowledge, this is the first study to demonstrate that "oxygen restriction" can extend lifespan in a mammalian model of aging.

Matrix prior for data transfer between single cell data types in latent Dirichlet allocation.

Single cell ATAC-seq (scATAC-seq) enables the mapping of regulatory elements in fine-grained cell types. Despite this advance, analysis of the resulting data is challenging, and large scale scATAC-seq data are difficult to obtain and expensive to generate. This motivates a method to leverage information from previously generated large scale scATAC-seq or scRNA-seq data to guide our analysis of new scATAC-seq datasets. We analyze scATAC-seq data using latent Dirichlet allocation (LDA), a Bayesian algorithm that was developed to model text corpora, summarizing documents as mixtures of topics defined based on the words that distinguish the documents. When applied to scATAC-seq, LDA treats cells as documents and their accessible sites as words, identifying "topics" based on the cell type-specific accessible sites in those cells. Previous work used uniform symmetric priors in LDA, but we hypothesized that nonuniform matrix priors generated from LDA models trained on existing data sets may enable improved detection of cell types in new data sets, especially if they have relatively few cells. In this work, we test this hypothesis in scATAC-seq data from whole C. elegans nematodes and SHARE-seq data from mouse skin cells. We show that nonsymmetric matrix priors for LDA improve our ability to capture cell type information from small scATAC-seq datasets.

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

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.

Structure-Based Design of Active-Site-Directed, Highly Potent, Selective, and Orally Bioavailable Low-Molecular-Weight Protein Tyrosine Phosphatase Inhibitors.

Protein tyrosine phosphatases constitute an important class of drug targets whose potential has been limited by the paucity of drug-like small-molecule inhibitors. We recently described a class of active-site-directed, moderately selective, and potent inhibitors of the low-molecular-weight protein tyrosine phosphatase (LMW-PTP). Here, we report our extensive structure-based design and optimization effort that afforded inhibitors with vastly improved potency and specificity. The leading compound inhibits LMW-PTP potently and selectively (Ki = 1.2 nM, >8000-fold selectivity). Many compounds exhibit favorable drug-like properties, such as low molecular weight, weak cytochrome P450 inhibition, high metabolic stability, moderate to high cell permeability (Papp > 0.2 nm/s), and moderate to good oral bioavailability (% F from 23 to 50% in mice), and therefore can be used as in vivo chemical probes to further dissect the complex biological as well as pathophysiological roles of LMW-PTP and for the development of therapeutics targeting LMW-PTP.

Comprehensive characterization of tissue-specific chromatin accessibility in L2 Caenorhabditis elegans nematodes.

Recently developed single-cell technologies allow researchers to characterize cell states at ever greater resolution and scale. Caenorhabditis elegans is a particularly tractable system for studying development, and recent single-cell RNA-seq studies characterized the gene expression patterns for nearly every cell type in the embryo and at the second larval stage (L2). Gene expression patterns give insight about gene function and into the biochemical state of different cell types; recent advances in other single-cell genomics technologies can now also characterize the regulatory context of the genome that gives rise to these gene expression levels at a single-cell resolution. To explore the regulatory DNA of individual cell types in C. elegans, we collected single-cell chromatin accessibility data using the sci-ATAC-seq assay in L2 larvae to match the available single-cell RNA-seq data set. By using a novel implementation of the latent Dirichlet allocation algorithm, we identify 37 clusters of cells that correspond to different cell types in the worm, providing new maps of putative cell type-specific gene regulatory sites, with promise for better understanding of cellular differentiation and gene regulation.

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.

Avocado: a multi-scale deep tensor factorization method learns a latent representation of the human epigenome.

The human epigenome has been experimentally characterized by thousands of measurements for every basepair in the human genome. We propose a deep neural network tensor factorization method, Avocado, that compresses this epigenomic data into a dense, information-rich representation. We use this learned representation to impute epigenomic data more accurately than previous methods, and we show that machine learning models that exploit this representation outperform those trained directly on epigenomic data on a variety of genomics tasks. These tasks include predicting gene expression, promoter-enhancer interactions, replication timing, and an element of 3D chromatin architecture.

Inferring Causal Gene Regulatory Networks from Coupled Single-Cell Expression Dynamics Using Scribe.

Here, we present Scribe (https://github.com/aristoteleo/Scribe-py), a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs restricted directed information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime"-ordered single-cell data compared with true time-series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as "RNA velocity" restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses highlight a shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and suggest ways of overcoming it.

Characterization of a novel AICARFT inhibitor which potently elevates ZMP and has anti-tumor activity in murine models.

AICARFT is a folate dependent catalytic site within the ATIC gene, part of the purine biosynthetic pathway, a pathway frequently upregulated in cancers. LSN3213128 is a potent (16 nM) anti-folate inhibitor of AICARFT and selective relative to TS, SHMT1, MTHFD1, MTHFD2 and MTHFD2L. Increases in ZMP, accompanied by activation of AMPK and cell growth inhibition, were observed with treatment of LY3213128. These effects on ZMP and proliferation were dependent on folate levels. In human breast MDA-MB-231met2 and lung NCI-H460 cell lines, growth inhibition was rescued by hypoxanthine, but not in the A9 murine cell line which is deficient in purine salvage. In athymic nude mice, LSN3213128 robustly elevates ZMP in MDA-MB-231met2, NCI-H460 and A9 tumors in a time and dose dependent manner. Significant tumor growth inhibition in human breast MDA-MB231met2 and lung NCI-H460 xenografts and in the syngeneic A9 tumor model were observed with oral administration of LSN3213128. Strikingly, AMPK appeared activated within the tumors and did not change even at high levels of intratumoral ZMP after weeks of dosing. These results support the evaluation of LSN3213128 as an antineoplastic agent.

Functional characterization of enhancer evolution in the primate lineage.

BACKGROUND: Enhancers play an important role in morphological evolution and speciation by controlling the spatiotemporal expression of genes. Previous efforts to understand the evolution of enhancers in primates have typically studied many enhancers at low resolution, or single enhancers at high resolution. Although comparative genomic studies reveal large-scale turnover of enhancers, a specific understanding of the molecular steps by which mammalian or primate enhancers evolve remains elusive. RESULTS: We identified candidate hominoid-specific liver enhancers from H3K27ac ChIP-seq data. After locating orthologs in 11 primates spanning around 40 million years, we synthesized all orthologs as well as computational reconstructions of 9 ancestral sequences for 348 active tiles of 233 putative enhancers. We concurrently tested all sequences for regulatory activity with STARR-seq in HepG2 cells. We observe groups of enhancer tiles with coherent trajectories, most of which can be potentially explained by a single gain or loss-of-activity event per tile. We quantify the correlation between the number of mutations along a branch and the magnitude of change in functional activity. Finally, we identify 84 mutations that correlate with functional changes; these are enriched for cytosine deamination events within CpGs. CONCLUSIONS: We characterized the evolutionary-functional trajectories of hundreds of liver enhancers throughout the primate phylogeny. We observe subsets of regulatory sequences that appear to have gained or lost activity. We use these data to quantify the relationship between sequence and functional divergence, and to identify CpG deamination as a potentially important force in driving changes in enhancer activity during primate evolution.

PREDICTD PaRallel Epigenomics Data Imputation with Cloud-based Tensor Decomposition.

The Encyclopedia of DNA Elements (ENCODE) and the Roadmap Epigenomics Project seek to characterize the epigenome in diverse cell types using assays that identify, for example, genomic regions with modified histones or accessible chromatin. These efforts have produced thousands of datasets but cannot possibly measure each epigenomic factor in all cell types. To address this, we present a method, PaRallel Epigenomics Data Imputation with Cloud-based Tensor Decomposition (PREDICTD), to computationally impute missing experiments. PREDICTD leverages an elegant model called "tensor decomposition" to impute many experiments simultaneously. Compared with the current state-of-the-art method, ChromImpute, PREDICTD produces lower overall mean squared error, and combining the two methods yields further improvement. We show that PREDICTD data captures enhancer activity at noncoding human accelerated regions. PREDICTD provides reference imputed data and open-source software for investigating new cell types, and demonstrates the utility of tensor decomposition and cloud computing, both promising technologies for bioinformatics.

CMPF, a Metabolite Formed Upon Prescription Omega-3-Acid Ethyl Ester Supplementation, Prevents and Reverses Steatosis.

Prescription ω-3 fatty acid ethyl ester supplements are commonly used for the treatment of hypertriglyceridemia. However, the metabolic profile and effect of the metabolites formed by these treatments remain unknown. Here we utilized unbiased metabolomics to identify 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF) as a significant metabolite of the ω-3-acid ethyl ester prescription Lovaza™ in humans. Administration of CMPF to mice before or after high-fat diet feeding at exposures equivalent to those observed in humans increased whole-body lipid metabolism, improved insulin sensitivity, increased beta-oxidation, reduced lipogenic gene expression, and ameliorated steatosis. Mechanistically, we find that CMPF acutely inhibits ACC activity, and induces long-term loss of SREBP1c and ACC1/2 expression. This corresponds to an induction of FGF21, which is required for long-term steatosis protection, as FGF21KO mice are refractory to the improved metabolic effects. Thus, CMPF treatment in mice parallels the effects of human Lovaza™ supplementation, revealing that CMPF may contribute to the improved metabolic effects observed with ω-3 fatty acid prescriptions.

The ModERN Resource: Genome-Wide Binding Profiles for Hundreds of Drosophila and Caenorhabditis elegans Transcription Factors.

To develop a catalog of regulatory sites in two major model organisms, Drosophila melanogaster and Caenorhabditis elegans, the modERN (model organism Encyclopedia of Regulatory Networks) consortium has systematically assayed the binding sites of transcription factors (TFs). Combined with data produced by our predecessor, modENCODE (Model Organism ENCyclopedia Of DNA Elements), we now have data for 262 TFs identifying 1.23 M sites in the fly genome and 217 TFs identifying 0.67 M sites in the worm genome. Because sites from different TFs are often overlapping and tightly clustered, they fall into 91,011 and 59,150 regions in the fly and worm, respectively, and these binding sites span as little as 8.7 and 5.8 Mb in the two organisms. Clusters with large numbers of sites (so-called high occupancy target, or HOT regions) predominantly associate with broadly expressed genes, whereas clusters containing sites from just a few factors are associated with genes expressed in tissue-specific patterns. All of the strains expressing GFP-tagged TFs are available at the stock centers, and the chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center and also through a simple interface (http://epic.gs.washington.edu/modERN/) that facilitates rapid accessibility of processed data sets. These data will facilitate a vast number of scientific inquiries into the function of individual TFs in key developmental, metabolic, and defense and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks and globally across the life spans of these two key model organisms.

Optimization of Hydroxyethylamine Transition State Isosteres as Aspartic Protease Inhibitors by Exploiting Conformational Preferences.

NMR conformational analysis of a hydroxyethylamine peptide isostere developed as an aspartic protease inhibitor shows that it is a flexible architecture. Cyclization to form pyrrolidines, piperidines, or morpholines results in a preorganization of the whole system in solution. The resulting conformation is similar to the conformation of the inhibitor in the active site of BACE-1. This entropic gain results in increased affinity for the enzyme when compared with the acyclic system. For morpholines 27 and 29, the combination of steric and electronic factors is exploited to orient substituents toward S1, S1', and S2' pockets both in the solution and in the bound states. These highly preorganized molecules proved to be the most potent compounds of the series. Additionally, the morpholines, unlike the pyrrolidine and piperidine analogues, have been found to be brain penetrant BACE-1 inhibitors.