Ribosome Levels Selectively Regulate Translation and Lineage Commitment in Human Hematopoiesis.

Blood cell formation is classically thought to occur through a hierarchical differentiation process, although recent studies have shown that lineage commitment may occur earlier in hematopoietic stem and progenitor cells (HSPCs). The relevance to human blood diseases and the underlying regulation of these refined models remain poorly understood. By studying a genetic blood disorder, Diamond-Blackfan anemia (DBA), where the majority of mutations affect ribosomal proteins and the erythroid lineage is selectively perturbed, we are able to gain mechanistic insight into how lineage commitment is programmed normally and disrupted in disease. We show that in DBA, the pool of available ribosomes is limited, while ribosome composition remains constant. Surprisingly, this global reduction in ribosome levels more profoundly alters translation of a select subset of transcripts. We show how the reduced translation of select transcripts in HSPCs can impair erythroid lineage commitment, illuminating a regulatory role for ribosome levels in cellular differentiation.

Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens.

Genetic screens help infer gene function in mammalian cells, but it has remained difficult to assay complex phenotypes-such as transcriptional profiles-at scale. Here, we develop Perturb-seq, combining single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-based perturbations to perform many such assays in a pool. We demonstrate Perturb-seq by analyzing 200,000 cells in immune cells and cell lines, focusing on transcription factors regulating the response of dendritic cells to lipopolysaccharide (LPS). Perturb-seq accurately identifies individual gene targets, gene signatures, and cell states affected by individual perturbations and their genetic interactions. We posit new functions for regulators of differentiation, the anti-viral response, and mitochondrial function during immune activation. By decomposing many high content measurements into the effects of perturbations, their interactions, and diverse cell metadata, Perturb-seq dramatically increases the scope of pooled genomic assays.

Epigenetic Memory Underlies Cell-Autonomous Heterogeneous Behavior of Hematopoietic Stem Cells.

Stem cells determine homeostasis and repair of many tissues and are increasingly recognized as functionally heterogeneous. To define the extent of-and molecular basis for-heterogeneity, we overlaid functional, transcriptional, and epigenetic attributes of hematopoietic stem cells (HSCs) at a clonal level using endogenous fluorescent tagging. Endogenous HSC had clone-specific functional attributes over time in vivo. The intra-clonal behaviors were highly stereotypic, conserved under the stress of transplantation, inflammation, and genotoxic injury, and associated with distinctive transcriptional, DNA methylation, and chromatin accessibility patterns. Further, HSC function corresponded to epigenetic configuration but not always to transcriptional state. Therefore, hematopoiesis under homeostatic and stress conditions represents the integrated action of highly heterogeneous clones of HSC with epigenetically scripted behaviors. This high degree of epigenetically driven cell autonomy among HSCs implies that refinement of the concepts of stem cell plasticity and of the stem cell niche is warranted.

Lactate activates the mitochondrial electron transport chain independently of its metabolism.

Lactate has long been considered a cellular waste product. However, we found that as extracellular lactate accumulates, it also enters the mitochondrial matrix and stimulates mitochondrial electron transport chain (ETC) activity. The resulting increase in mitochondrial ATP synthesis suppresses glycolysis and increases the utilization of pyruvate and/or alternative respiratory substrates. The ability of lactate to increase oxidative phosphorylation does not depend on its metabolism. Both L- and D-lactate are effective at enhancing ETC activity and suppressing glycolysis. Furthermore, the selective induction of mitochondrial oxidative phosphorylation by unmetabolized D-lactate reversibly suppressed aerobic glycolysis in both cancer cell lines and proliferating primary cells in an ATP-dependent manner and enabled cell growth on respiratory-dependent bioenergetic substrates. In primary T cells, D-lactate enhanced cell proliferation and effector function. Together, these findings demonstrate that lactate is a critical regulator of the ability of mitochondrial oxidative phosphorylation to suppress glucose fermentation.

Resilient anatomy and local plasticity of naive and stress haematopoiesis.

The bone marrow adjusts blood cell production to meet physiological demands in response to insults. The spatial organization of normal and stress responses are unknown owing to the lack of methods to visualize most steps of blood production. Here we develop strategies to image multipotent haematopoiesis, erythropoiesis and lymphopoiesis in mice. We combine these with imaging of myelopoiesis1 to define the anatomy of normal and stress haematopoiesis. In the steady state, across the skeleton, single stem cells and multipotent progenitors distribute through the marrow enriched near megakaryocytes. Lineage-committed progenitors are recruited to blood vessels, where they contribute to lineage-specific microanatomical structures composed of progenitors and immature cells, which function as the production sites for each major blood lineage. This overall anatomy is resilient to insults, as it was maintained after haemorrhage, systemic bacterial infection and granulocyte colony-stimulating factor (G-CSF) treatment, and during ageing. Production sites enable haematopoietic plasticity as they differentially and selectively modulate their numbers and output in response to insults. We found that stress responses are variable across the skeleton: the tibia and the sternum respond in opposite ways to G-CSF, and the skull does not increase erythropoiesis after haemorrhage. Our studies enable in situ analyses of haematopoiesis, define the anatomy of normal and stress responses, identify discrete microanatomical production sites that confer plasticity to haematopoiesis, and uncover unprecedented heterogeneity of stress responses across the skeleton.

Circulating tumor cell clusters are oligoclonal precursors of breast cancer metastasis.

Circulating tumor cell clusters (CTC clusters) are present in the blood of patients with cancer but their contribution to metastasis is not well defined. Using mouse models with tagged mammary tumors, we demonstrate that CTC clusters arise from oligoclonal tumor cell groupings and not from intravascular aggregation events. Although rare in the circulation compared with single CTCs, CTC clusters have 23- to 50-fold increased metastatic potential. In patients with breast cancer, single-cell resolution RNA sequencing of CTC clusters and single CTCs, matched within individual blood samples, identifies the cell junction component plakoglobin as highly differentially expressed. In mouse models, knockdown of plakoglobin abrogates CTC cluster formation and suppresses lung metastases. In breast cancer patients, both abundance of CTC clusters and high tumor plakoglobin levels denote adverse outcomes. Thus, CTC clusters are derived from multicellular groupings of primary tumor cells held together through plakoglobin-dependent intercellular adhesion, and though rare, they greatly contribute to the metastatic spread of cancer.

Compatibility rules of human enhancer and promoter sequences.

Gene regulation in the human genome is controlled by distal enhancers that activate specific nearby promoters1. A proposed model for this specificity is that promoters have sequence-encoded preferences for certain enhancers, for example, mediated by interacting sets of transcription factors or cofactors2. This 'biochemical compatibility' model has been supported by observations at individual human promoters and by genome-wide measurements in Drosophila3-9. However, the degree to which human enhancers and promoters are intrinsically compatible has not yet been systematically measured, and how their activities combine to control RNA expression remains unclear. Here we design a high-throughput reporter assay called enhancer × promoter self-transcribing active regulatory region sequencing (ExP STARR-seq) and applied it to examine the combinatorial compatibilities of 1,000 enhancer and 1,000 promoter sequences in human K562 cells. We identify simple rules for enhancer-promoter compatibility, whereby most enhancers activate all promoters by similar amounts, and intrinsic enhancer and promoter activities multiplicatively combine to determine RNA output (R2 = 0.82). In addition, two classes of enhancers and promoters show subtle preferential effects. Promoters of housekeeping genes contain built-in activating motifs for factors such as GABPA and YY1, which decrease the responsiveness of promoters to distal enhancers. Promoters of variably expressed genes lack these motifs and show stronger responsiveness to enhancers. Together, this systematic assessment of enhancer-promoter compatibility suggests a multiplicative model tuned by enhancer and promoter class to control gene transcription in the human genome.

A persistent ultraviolet outflow from an accreting neutron star binary transient.

All disc-accreting astrophysical objects produce powerful disc winds. In compact binaries containing neutron stars or black holes, accretion often takes place during violent outbursts. The main disc wind signatures during these eruptions are blue-shifted X-ray absorption lines, which are preferentially seen in disc-dominated 'soft states'1,2. By contrast, optical wind-formed lines have recently been detected in 'hard states', when a hot corona dominates the luminosity3. The relationship between these signatures is unknown, and no erupting system has as yet revealed wind-formed lines between the X-ray and optical bands, despite the many strong resonance transitions in this ultraviolet (UV) region4. Here we report that the transient neutron star binary Swift J1858.6-0814 exhibits wind-formed, blue-shifted absorption lines associated with C IV, N V and He II in time-resolved UV spectroscopy during a luminous hard state, which we interpret as a warm, moderately ionized outflow component in this state. Simultaneously observed optical lines also display transient blue-shifted absorption. Decomposing the UV data into constant and variable components, the blue-shifted absorption is associated with the former. This implies that the outflow is not associated with the luminous flares in the data. The joint presence of UV and optical wind features reveals a multi-phase and/or spatially stratified evaporative outflow from the outer disc5. This type of persistent mass loss across all accretion states has been predicted by radiation-hydrodynamic simulations6 and helps to explain the shorter-than-expected duration of outbursts7.

Lipid Droplet-Derived Monounsaturated Fatty Acids Traffic via PLIN5 to Allosterically Activate SIRT1.

Lipid droplets (LDs) provide a reservoir for triacylglycerol storage and are a central hub for fatty acid trafficking and signaling in cells. Lipolysis promotes mitochondrial biogenesis and oxidative metabolism via a SIRT1/PGC-1α/PPARα-dependent pathway through an unknown mechanism. Herein, we identify that monounsaturated fatty acids (MUFAs) allosterically activate SIRT1 toward select peptide-substrates such as PGC-1α. MUFAs enhance PGC-1α/PPARα signaling and promote oxidative metabolism in cells and animal models in a SIRT1-dependent manner. Moreover, we characterize the LD protein perilipin 5 (PLIN5), which is known to enhance mitochondrial biogenesis and function, to be a fatty-acid-binding protein that preferentially binds LD-derived monounsaturated fatty acids and traffics them to the nucleus following cAMP/PKA-mediated lipolytic stimulation. Thus, these studies identify the first-known endogenous allosteric modulators of SIRT1 and characterize a LD-nuclear signaling axis that underlies the known metabolic benefits of MUFAs and PLIN5.

Cycling cancer persister cells arise from lineages with distinct programs.

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure1, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment2. Although most cancer persisters remain arrested in the presence of the drug, a rare subset can re-enter the cell cycle under constitutive drug treatment. Little is known about the non-genetic mechanisms that enable cancer persisters to maintain proliferative capacity in the presence of drugs. To study this rare, transiently resistant, proliferative persister population, we developed Watermelon, a high-complexity expressed barcode lentiviral library for simultaneous tracing of each cell's clonal origin and proliferative and transcriptional states. Here we show that cycling and non-cycling persisters arise from different cell lineages with distinct transcriptional and metabolic programs. Upregulation of antioxidant gene programs and a metabolic shift to fatty acid oxidation are associated with persister proliferative capacity across multiple cancer types. Impeding oxidative stress or metabolic reprogramming alters the fraction of cycling persisters. In human tumours, programs associated with cycling persisters are induced in minimal residual disease in response to multiple targeted therapies. The Watermelon system enabled the identification of rare persister lineages that are preferentially poised to proliferate under drug pressure, thus exposing new vulnerabilities that can be targeted to delay or even prevent disease recurrence.

Genome-wide enhancer maps link risk variants to disease genes.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease1. Many of the underlying causal variants may affect enhancers2,3, but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types4. Here we apply this ABC model to create enhancer-gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions.

MAFG-driven astrocytes promote CNS inflammation.

Multiple sclerosis is a chronic inflammatory disease of the CNS1. Astrocytes contribute to the pathogenesis of multiple sclerosis2, but little is known about the heterogeneity of astrocytes and its regulation. Here we report the analysis of astrocytes in multiple sclerosis and its preclinical model experimental autoimmune encephalomyelitis (EAE) by single-cell RNA sequencing in combination with cell-specific Ribotag RNA profiling, assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), genome-wide analysis of DNA methylation and in vivo CRISPR-Cas9-based genetic perturbations. We identified astrocytes in EAE and multiple sclerosis that were characterized by decreased expression of NRF2 and increased expression of MAFG, which cooperates with MAT2α to promote DNA methylation and represses antioxidant and anti-inflammatory transcriptional programs. Granulocyte-macrophage colony-stimulating factor (GM-CSF) signalling in astrocytes drives the expression of MAFG and MAT2α and pro-inflammatory transcriptional modules, contributing to CNS pathology in EAE and, potentially, multiple sclerosis. Our results identify candidate therapeutic targets in multiple sclerosis.

Live-animal imaging of native haematopoietic stem and progenitor cells.

The biology of haematopoietic stem cells (HSCs) has predominantly been studied under transplantation conditions1,2. It has been particularly challenging to study dynamic HSC behaviour, given that the visualization of HSCs in the native niche in live animals has not, to our knowledge, been achieved. Here we describe a dual genetic strategy in mice that restricts reporter labelling to a subset of the most quiescent long-term HSCs (LT-HSCs) and that is compatible with current intravital imaging approaches in the calvarial bone marrow3-5. We show that this subset of LT-HSCs resides close to both sinusoidal blood vessels and the endosteal surface. By contrast, multipotent progenitor cells (MPPs) show greater variation in distance from the endosteum and are more likely to be associated with transition zone vessels. LT-HSCs are not found in bone marrow niches with the deepest hypoxia and instead are found in hypoxic environments similar to those of MPPs. In vivo time-lapse imaging revealed that LT-HSCs at steady-state show limited motility. Activated LT-HSCs show heterogeneous responses, with some cells becoming highly motile and a fraction of HSCs expanding clonally within spatially restricted domains. These domains have defined characteristics, as HSC expansion is found almost exclusively in a subset of bone marrow cavities with bone-remodelling activity. By contrast, cavities with low bone-resorbing activity do not harbour expanding HSCs. These findings point to previously unknown heterogeneity within the bone marrow microenvironment, imposed by the stages of bone turnover. Our approach enables the direct visualization of HSC behaviours and dissection of heterogeneity in HSC niches.

Variability in the analysis of a single neuroimaging dataset by many teams.

Data analysis workflows in many scientific domains have become increasingly complex and flexible. Here we assess the effect of this flexibility on the results of functional magnetic resonance imaging by asking 70 independent teams to analyse the same dataset, testing the same 9 ex-ante hypotheses1. The flexibility of analytical approaches is exemplified by the fact that no two teams chose identical workflows to analyse the data. This flexibility resulted in sizeable variation in the results of hypothesis tests, even for teams whose statistical maps were highly correlated at intermediate stages of the analysis pipeline. Variation in reported results was related to several aspects of analysis methodology. Notably, a meta-analytical approach that aggregated information across teams yielded a significant consensus in activated regions. Furthermore, prediction markets of researchers in the field revealed an overestimation of the likelihood of significant findings, even by researchers with direct knowledge of the dataset2-5. Our findings show that analytical flexibility can have substantial effects on scientific conclusions, and identify factors that may be related to variability in the analysis of functional magnetic resonance imaging. The results emphasize the importance of validating and sharing complex analysis workflows, and demonstrate the need for performing and reporting multiple analyses of the same data. Potential approaches that could be used to mitigate issues related to analytical variability are discussed.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes.

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2-4 and coronary heart disease5-this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Expansion of the sagittal suture induces proliferation of skeletal stem cells and sustains endogenous calvarial bone regeneration.

In newborn humans, and up to approximately 2 y of age, calvarial bone defects can naturally regenerate. This remarkable regeneration potential is also found in newborn mice and is absent in adult mice. Since previous studies showed that the mouse calvarial sutures are reservoirs of calvarial skeletal stem cells (cSSCs), which are the cells responsible for calvarial bone regeneration, here we hypothesized that the regenerative potential of the newborn mouse calvaria is due to a significant amount of cSSCs present in the newborn expanding sutures. Thus, we tested whether such regenerative potential can be reverse engineered in adult mice by artificially inducing an increase of the cSSCs resident within the adult calvarial sutures. First, we analyzed the cellular composition of the calvarial sutures in newborn and in older mice, up to 14-mo-old mice, showing that the sutures of the younger mice are enriched in cSSCs. Then, we demonstrated that a controlled mechanical expansion of the functionally closed sagittal sutures of adult mice induces a significant increase of the cSSCs. Finally, we showed that if a calvarial critical size bone defect is created simultaneously to the mechanical expansion of the sagittal suture, it fully regenerates without the need for additional therapeutic aids. Using a genetic blockade system, we further demonstrate that this endogenous regeneration is mediated by the canonical Wnt signaling. This study shows that controlled mechanical forces can harness the cSSCs and induce calvarial bone regeneration. Similar harnessing strategies may be used to develop novel and more effective bone regeneration autotherapies.

Identification of a ß-arrestin-biased negative allosteric modulator for the ß2-adrenergic receptor.

Catecholamine-stimulated ß2-adrenergic receptor (ß2AR) signaling via the canonical Gs-adenylyl cyclase-cAMP-PKA pathway regulates numerous physiological functions, including the therapeutic effects of exogenous ß-agonists in the treatment of airway disease. ß2AR signaling is tightly regulated by GRKs and ß-arrestins, which together promote ß2AR desensitization and internalization as well as downstream signaling, often antithetical to the canonical pathway. Thus, the ability to bias ß2AR signaling toward the Gs pathway while avoiding ß-arrestin-mediated effects may provide a strategy to improve the functional consequences of ß2AR activation. Since attempts to develop Gs-biased agonists and allosteric modulators for the ß2AR have been largely unsuccessful, here we screened small molecule libraries for allosteric modulators that selectively inhibit ß-arrestin recruitment to the receptor. This screen identified several compounds that met this profile, and, of these, a difluorophenyl quinazoline (DFPQ) derivative was found to be a selective negative allosteric modulator of ß-arrestin recruitment to the ß2AR while having no effect on ß2AR coupling to Gs. DFPQ effectively inhibits agonist-promoted phosphorylation and internalization of the ß2AR and protects against the functional desensitization of ß-agonist mediated regulation in cell and tissue models. The effects of DFPQ were also specific to the ß2AR with minimal effects on the ß1AR. Modeling, mutagenesis, and medicinal chemistry studies support DFPQ derivatives binding to an intracellular membrane-facing region of the ß2AR, including residues within transmembrane domains 3 and 4 and intracellular loop 2. DFPQ thus represents a class of biased allosteric modulators that targets an allosteric site of the ß2AR.

Image-seq: spatially resolved single-cell sequencing guided by in situ and in vivo imaging.

Tissue function depends on cellular organization. While the properties of individual cells are increasingly being deciphered using powerful single-cell sequencing technologies, understanding their spatial organization and temporal evolution remains a major challenge. Here, we present Image-seq, a technology that provides single-cell transcriptional data on cells that are isolated from specific spatial locations under image guidance, thus preserving the spatial information of the target cells. It is compatible with in situ and in vivo imaging and can document the temporal and dynamic history of the cells being analyzed. Cell samples are isolated from intact tissue and processed with state-of-the-art library preparation protocols. The technique therefore combines spatial information with highly sensitive RNA sequencing readouts from individual, intact cells. We have used both high-throughput, droplet-based sequencing as well as SMARTseq-v4 library preparation to demonstrate its application to bone marrow and leukemia biology. We discovered that DPP4 is a highly upregulated gene during early progression of acute myeloid leukemia and that it marks a more proliferative subpopulation that is confined to specific bone marrow microenvironments. Furthermore, the ability of Image-seq to isolate viable, intact cells should make it compatible with a range of downstream single-cell analysis tools including multi-omics protocols.

Intellectual disability and autism in propionic acidemia: a biomarker-behavioral investigation implicating dysregulated mitochondrial biology.

Propionic acidemia (PA) is an autosomal recessive condition (OMIM #606054), wherein pathogenic variants in PCCA and PCCB impair the activity of propionyl-CoA carboxylase. PA is associated with neurodevelopmental disorders, including intellectual disability (ID) and autism spectrum disorder (ASD); however, the correlates and mechanisms of these outcomes remain unknown. Using data from a subset of participants with PA enrolled in a dedicated natural history study (n = 33), we explored associations between neurodevelopmental phenotypes and laboratory parameters. Twenty (61%) participants received an ID diagnosis, and 12 of the 31 (39%) who were fully evaluated received the diagnosis of ASD. A diagnosis of ID, lower full-scale IQ (sample mean = 65 ± 26), and lower adaptive behavior composite scores (sample mean = 67 ± 23) were associated with several biomarkers. Higher concentrations of plasma propionylcarnitine, plasma total 2-methylcitrate, serum erythropoietin, and mitochondrial biomarkers plasma FGF21 and GDF15 were associated with a more severe ID profile. Reduced 1-13C-propionate oxidative capacity and decreased levels of plasma and urinary glutamine were also associated with a more severe ID profile. Only two parameters, increased serum erythropoietin and decreased plasma glutamine, were associated with ASD. Plasma glycine, one of the defining features of PA, was not meaningfully associated with either ID or ASD. Thus, while both ID and ASD were commonly observed in our PA cohort, only ID was robustly associated with metabolic parameters. Our results suggest that disease severity and associated mitochondrial dysfunction may play a role in CNS complications of PA and identify potential biomarkers and candidate surrogate endpoints.