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.

Alveolin proteins in the Toxoplasma inner membrane complex form a highly interconnected structure that maintains parasite shape and replication.

Apicomplexan parasites possess several specialized structures to invade their host cells and replicate successfully. One of these is the inner membrane complex (IMC), a peripheral membrane-cytoskeletal system underneath the plasma membrane. It is composed of a series of flattened, membrane-bound vesicles and a cytoskeletal subpellicular network (SPN) comprised of intermediate filament-like proteins called alveolins. While the alveolin proteins are conserved throughout the Apicomplexa and the broader Alveolata, their precise functions and interactions remain poorly understood. Here, we describe the function of one of these alveolin proteins in Toxoplasma, IMC6. Disruption of IMC6 resulted in striking morphological defects that led to aberrant invasion and replication but surprisingly minor effects on motility. Deletion analyses revealed that the alveolin domain alone is largely sufficient to restore localization and partially sufficient for function. As this highlights the importance of the IMC6 alveolin domain, we implemented unnatural amino acid photoreactive crosslinking to the alveolin domain and identified multiple binding interfaces between IMC6 and 2 other cytoskeletal IMC proteins-IMC3 and ILP1. This provides direct evidence of protein-protein interactions in the alveolin domain and supports the long-held hypothesis that the alveolin domain is responsible for filament formation. Collectively, our study features the conserved alveolin proteins as critical components that maintain the parasite's structural integrity and highlights the alveolin domain as a key mediator of SPN architecture.

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.

Ribosomal frameshifting selectively modulates the assembly, function, and pharmacological rescue of a misfolded CFTR variant.

The cotranslational misfolding of the cystic fibrosis transmembrane conductance regulator chloride channel (CFTR) plays a central role in the molecular basis of CF. The misfolding of the most common CF variant (ΔF508) remodels both the translational regulation and quality control of CFTR. Nevertheless, it is unclear how the misassembly of the nascent polypeptide may directly influence the activity of the translation machinery. In this work, we identify a structural motif within the CFTR transcript that stimulates efficient -1 ribosomal frameshifting and triggers the premature termination of translation. Though this motif does not appear to impact the interactome of wild-type CFTR, silent mutations that disrupt this RNA structure alter the association of nascent ΔF508 CFTR with numerous translation and quality control proteins. Moreover, disrupting this RNA structure enhances the functional gating of the ΔF508 CFTR channel at the plasma membrane and its pharmacological rescue by the CFTR modulators contained in the CF drug Trikafta. The effects of the RNA structure on ΔF508 CFTR appear to be attenuated in the absence of the ER membrane protein complex, which was previously found to modulate ribosome collisions during "preemptive quality control" of a misfolded CFTR homolog. Together, our results reveal that ribosomal frameshifting selectively modulates the assembly, function, and pharmacological rescue of a misfolded CFTR variant. These findings suggest that interactions between the nascent chain, quality control machinery, and ribosome may dynamically modulate ribosomal frameshifting in order to tune the processivity of translation in response to cotranslational misfolding.

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.

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.

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.

SMCHD1 activates the expression of genes required for the expansion of human myoblasts.

SMCHD1 is an epigenetic regulatory protein known to modulate the targeted repression of large chromatin domains. Diminished SMCHD1 function in muscle fibers causes Facioscapulohumeral Muscular Dystrophy (FSHD2) through derepression of the D4Z4 chromatin domain, an event which permits the aberrant expression of the disease-causing gene DUX4. Given that SMCHD1 plays a broader role in establishing the cellular epigenome, we examined whether loss of SMCHD1 function might affect muscle homeostasis through additional mechanisms. Here we show that acute depletion of SMCHD1 results in a DUX4-independent defect in myoblast proliferation. Genomic and transcriptomic experiments determined that SMCHD1 associates with enhancers of genes controlling cell cycle to activate their expression. Amongst these cell cycle regulatory genes, we identified LAP2 as a key target of SMCHD1 required for the expansion of myoblasts, where the ectopic expression of LAP2 rescues the proliferation defect of SMCHD1-depleted cells. Thus, the epigenetic regulator SMCHD1 can play the role of a transcriptional co-activator for maintaining the expression of genes required for muscle progenitor expansion. This DUX4-independent role for SMCHD1 in myoblasts suggests that the pathology of FSHD2 may be a consequence of defective muscle regeneration in addition to the muscle wasting caused by spurious DUX4 expression.

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.