A combination treatment based on drug repurposing demonstrates mutation-agnostic efficacy in pre-clinical retinopathy models.

Inherited retinopathies are devastating diseases that in most cases lack treatment options. Disease-modifying therapies that mitigate pathophysiology regardless of the underlying genetic lesion are desirable due to the diversity of mutations found in such diseases. We tested a systems pharmacology-based strategy that suppresses intracellular cAMP and Ca2+ activity via G protein-coupled receptor (GPCR) modulation using tamsulosin, metoprolol, and bromocriptine coadministration. The treatment improves cone photoreceptor function and slows degeneration in Pde6ßrd10 and RhoP23H/WT retinitis pigmentosa mice. Cone degeneration is modestly mitigated after a 7-month-long drug infusion in PDE6A-/- dogs. The treatment also improves rod pathway function in an Rpe65-/- mouse model of Leber congenital amaurosis but does not protect from cone degeneration. RNA-sequencing analyses indicate improved metabolic function in drug-treated Rpe65-/- and rd10 mice. Our data show that catecholaminergic GPCR drug combinations that modify second messenger levels via multiple receptor actions provide a potential disease-modifying therapy against retinal degeneration.

Conditional deletion of miR-204 and miR-211 in murine retinal pigment epithelium results in retinal degeneration.

MicroRNAs (miRs) are short, evolutionarily conserved noncoding RNAs that canonically downregulate expression of target genes. The miR family composed of miR-204 and miR-211 is among the most highly expressed miRs in the retinal pigment epithelium (RPE) in both mouse and human and also retains high sequence identity. To assess the role of this miR family in the developed mouse eye, we generated two floxed conditional KO mouse lines crossed to the RPE65-ERT2-Cre driver mouse line to perform an RPE-specific conditional KO of this miR family in adult mice. After Cre-mediated deletion, we observed retinal structural changes by optical coherence tomography; dysfunction and loss of photoreceptors by retinal imaging; and retinal inflammation marked by subretinal infiltration of immune cells by imaging and immunostaining. Single-cell RNA sequencing of diseased RPE and retinas showed potential miR-regulated target genes, as well as changes in noncoding RNAs in the RPE, rod photoreceptors, and Müller glia. This work thus highlights the role of miR-204 and miR-211 in maintaining RPE function and how the loss of miRs in the RPE exerts effects on the neural retina, leading to inflammation and retinal degeneration.

Restoring retinal polyunsaturated fatty acid balance and retina function by targeting ceramide in AdipoR1-deficient mice.

Mutations in the adiponectin receptor 1 gene (AdipoR1) lead to retinitis pigmentosa and are associated with age-related macular degeneration. This study explores the effects of AdipoR1 gene deficiency in mice, revealing a striking decline in ω3 polyunsaturated fatty acids (PUFA), an increase in ω6 fatty acids, and elevated ceramides in the retina. The AdipoR1 deficiency impairs peroxisome proliferator-activated receptor α signaling, which is crucial for FA metabolism, particularly affecting proteins associated with FA transport and oxidation in the retina and retinal pigmented epithelium. Our lipidomic and proteomic analyses indicate changes that could affect membrane composition and viscosity through altered ω3 PUFA transport and synthesis, suggesting a potential influence of AdipoR1 on these properties. Furthermore, we noted a reduction in the Bardet-Biedl syndrome proteins, which are crucial for forming and maintaining photoreceptor outer segments that are PUFA-enriched ciliary structures. Diminution in Bardet-Biedl syndrome-proteins content combined with our electron microscopic observations raises the possibility that AdipoR1 deficiency might impair ciliary function. Treatment with inhibitors of ceramide synthesis led to substantial elevation of ω3 LC-PUFAs, alleviating photoreceptor degeneration and improving retinal function. These results serve as the proof of concept for a ceramide-targeted strategy to treat retinopathies linked to PUFA deficiency, including age-related macular degeneration.

Distinct mouse models of Stargardt disease display differences in pharmacological targeting of ceramides and inflammatory responses.

Mutations in many visual cycle enzymes in photoreceptors and retinal pigment epithelium (RPE) cells can lead to the chronic accumulation of toxic retinoid byproducts, which poison photoreceptors and the underlying RPE if left unchecked. Without a functional ATP-binding cassette, sub-family A, member 4 (ABCA4), there is an elevation of all-trans-retinal and prolonged buildup of all-trans-retinal adducts, resulting in a retinal degenerative disease known as Stargardt-1 disease. Even in this monogenic disorder, there is significant heterogeneity in the time to onset of symptoms among patients. Using a combination of molecular techniques, we studied Abca4 knockout (simulating human noncoding disease variants) and Abca4 knock-in mice (simulating human misfolded, catalytically inactive protein variants), which serve as models for Stargardt-1 disease. We compared the two strains to ascertain whether they exhibit differential responses to agents that affect cytokine signaling and/or ceramide metabolism, as alterations in either of these pathways can exacerbate retinal degenerative phenotypes. We found different degrees of responsiveness to maraviroc, a known immunomodulatory CCR5 antagonist, and to the ceramide-lowering agent AdipoRon, an agonist of the ADIPOR1 and ADIPOR2 receptors. The two strains also display different degrees of transcriptional deviation from matched WT controls. Our phenotypic comparison of the two distinct Abca4 mutant-mouse models sheds light on potential therapeutic avenues previously unexplored in the treatment of Stargardt disease and provides a surrogate assay for assessing the effectiveness for genome editing.

spinDrop: a droplet microfluidic platform to maximise single-cell sequencing information content.

Droplet microfluidic methods have massively increased the throughput of single-cell sequencing campaigns. The benefit of scale-up is, however, accompanied by increased background noise when processing challenging samples and the overall RNA capture efficiency is lower. These drawbacks stem from the lack of strategies to enrich for high-quality material or specific cell types at the moment of cell encapsulation and the absence of implementable multi-step enzymatic processes that increase capture. Here we alleviate both bottlenecks using fluorescence-activated droplet sorting to enrich for droplets that contain single viable cells, intact nuclei, fixed cells or target cell types and use reagent addition to droplets by picoinjection to perform multi-step lysis and reverse transcription. Our methodology increases gene detection rates fivefold, while reducing background noise by up to half. We harness these properties to deliver a high-quality molecular atlas of mouse brain development, despite starting with highly damaged input material, and provide an atlas of nascent RNA transcription during mouse organogenesis. Our method is broadly applicable to other droplet-based workflows to deliver sensitive and accurate single-cell profiling at a reduced cost.

Stress resilience-enhancing drugs preserve tissue structure and function in degenerating retina via phosphodiesterase inhibition.

Chronic, progressive retinal diseases, such as age-related macular degeneration (AMD), diabetic retinopathy, and retinitis pigmentosa, arise from genetic and environmental perturbations of cellular and tissue homeostasis. These disruptions accumulate with repeated exposures to stress over time, leading to progressive visual impairment and, in many cases, legal blindness. Despite decades of research, therapeutic options for the millions of patients suffering from these disorders remain severely limited, especially for treating earlier stages of pathogenesis when the opportunity to preserve the retinal structure and visual function is greatest. To address this urgent, unmet medical need, we employed a systems pharmacology platform for therapeutic development. Through integrative single-cell transcriptomics, proteomics, and phosphoproteomics, we identified universal molecular mechanisms across distinct models of age-related and inherited retinal degenerations, characterized by impaired physiological resilience to stress. Here, we report that selective, targeted pharmacological inhibition of cyclic nucleotide phosphodiesterases (PDEs), which serve as critical regulatory nodes that modulate intracellular second messenger signaling pathways, stabilized the transcriptome, proteome, and phosphoproteome through downstream activation of protective mechanisms coupled with synergistic inhibition of degenerative processes. This therapeutic intervention enhanced resilience to acute and chronic forms of stress in the degenerating retina, thus preserving tissue structure and function across various models of age-related and inherited retinal disease. Taken together, these findings exemplify a systems pharmacology approach to drug discovery and development, revealing a new class of therapeutics with potential clinical utility in the treatment or prevention of the most common causes of blindness.

Tim-3 adapter protein Bat3 acts as an endogenous regulator of tolerogenic dendritic cell function.

Dendritic cells (DCs) sense environmental cues and adopt either an immune-stimulatory or regulatory phenotype, thereby fine-tuning immune responses. Identifying endogenous regulators that determine DC function can thus inform the development of therapeutic strategies for modulating the immune response in different disease contexts. Tim-3 plays an important role in regulating immune responses by inhibiting the activation status and the T cell priming ability of DC in the setting of cancer. Bat3 is an adaptor protein that binds to the tail of Tim-3; therefore, we studied its role in regulating the functional status of DCs. In murine models of autoimmunity (experimental autoimmune encephalomyelitis) and cancer (MC38-OVA-implanted tumor), lack of Bat3 expression in DCs alters the T cell compartment-it decreases TH1, TH17 and cytotoxic effector cells, increases regulatory T cells, and exhausted CD8+ tumor-infiltrating lymphocytes, resulting in the attenuation of autoimmunity and acceleration of tumor growth. We found that Bat3 expression levels were differentially regulated by activating versus inhibitory stimuli in DCs, indicating a role for Bat3 in the functional calibration of DC phenotypes. Mechanistically, loss of Bat3 in DCs led to hyperactive unfolded protein response and redirected acetyl-coenzyme A to increase cell intrinsic steroidogenesis. The enhanced steroidogenesis in Bat3-deficient DC suppressed T cell response in a paracrine manner. Our findings identified Bat3 as an endogenous regulator of DC function, which has implications for DC-based immunotherapies.

Inhibition of ceramide accumulation in AdipoR1-/- mice increases photoreceptor survival and improves vision.

Adiponectin receptor 1 (ADIPOR1) is a lipid and glucose metabolism regulator that possesses intrinsic ceramidase activity. Mutations of the ADIPOR1 gene have been associated with nonsyndromic and syndromic retinitis pigmentosa. Here, we show that the absence of AdipoR1 in mice leads to progressive photoreceptor degeneration, significant reduction of electroretinogram amplitudes, decreased retinoid content in the retina, and reduced cone opsin expression. Single-cell RNA-Seq results indicate that ADIPOR1 encoded the most abundantly expressed ceramidase in mice and one of the 2 most highly expressed ceramidases in the human retina, next to acid ceramidase ASAH1. We discovered an accumulation of ceramides in the AdipoR1-/- retina, likely due to insufficient ceramidase activity for healthy retina function, resulting in photoreceptor death. Combined treatment with desipramine/L-cycloserine (DC) lowered ceramide levels and exerted a protective effect on photoreceptors in AdipoR1-/- mice. Moreover, we observed improvement in cone-mediated retinal function in the DC-treated animals. Lastly, we found that prolonged DC treatment corrected the electrical responses of the primary visual cortex to visual stimuli, approaching near-normal levels for some parameters. These results highlight the importance of ADIPOR1 ceramidase in the retina and show that pharmacological inhibition of ceramide generation can provide a therapeutic strategy for ADIPOR1-related retinopathy.

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.

TIM-3 restrains anti-tumour immunity by regulating inflammasome activation.

T cell immunoglobulin and mucin-containing molecule 3 (TIM-3), first identified as a molecule expressed on interferon-γ producing T cells1, is emerging as an important immune-checkpoint molecule, with therapeutic blockade of TIM-3 being investigated in multiple human malignancies. Expression of TIM-3 on CD8+ T cells in the tumour microenvironment is considered a cardinal sign of T cell dysfunction; however, TIM-3 is also expressed on several other types of immune cell, confounding interpretation of results following blockade using anti-TIM-3 monoclonal antibodies. Here, using conditional knockouts of TIM-3 together with single-cell RNA sequencing, we demonstrate the singular importance of TIM-3 on dendritic cells (DCs), whereby loss of TIM-3 on DCs-but not on CD4+ or CD8+ T cells-promotes strong anti-tumour immunity. Loss of TIM-3 prevented DCs from expressing a regulatory program and facilitated the maintenance of CD8+ effector and stem-like T cells. Conditional deletion of TIM-3 in DCs led to increased accumulation of reactive oxygen species resulting in NLRP3 inflammasome activation. Inhibition of inflammasome activation, or downstream effector cytokines interleukin-1ß (IL-1ß) and IL-18, completely abrogated the protective anti-tumour immunity observed with TIM-3 deletion in DCs. Together, our findings reveal an important role for TIM-3 in regulating DC function and underscore the potential of TIM-3 blockade in promoting anti-tumour immunity by regulating inflammasome activation.

Inherited myeloproliferative neoplasm risk affects haematopoietic stem cells.

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P < 5.0 × 10-8), 7 of which have not been previously reported. We find that there is a shared genetic architecture between MPN risk and several haematopoietic traits from distinct lineages; that there is an enrichment for MPN risk variants within accessible chromatin of HSCs; and that increased MPN risk is associated with longer telomere length in leukocytes and other clonal haematopoietic states-collectively suggesting that MPN risk is associated with the function and self-renewal of HSCs. We use gene mapping to identify modulators of HSC biology linked to MPN risk, and show through targeted variant-to-function assays that CHEK2 and GFI1B have roles in altering the function of HSCs to confer disease risk. Overall, our results reveal a previously unappreciated mechanism for inherited MPN risk through the modulation of HSC function.

Cumulus provides cloud-based data analysis for large-scale single-cell and single-nucleus RNA-seq.

Massively parallel single-cell and single-nucleus RNA sequencing has opened the way to systematic tissue atlases in health and disease, but as the scale of data generation is growing, so is the need for computational pipelines for scaled analysis. Here we developed Cumulus-a cloud-based framework for analyzing large-scale single-cell and single-nucleus RNA sequencing datasets. Cumulus combines the power of cloud computing with improvements in algorithm and implementation to achieve high scalability, low cost, user-friendliness and integrated support for a comprehensive set of features. We benchmark Cumulus on the Human Cell Atlas Census of Immune Cells dataset of bone marrow cells and show that it substantially improves efficiency over conventional frameworks, while maintaining or improving the quality of results, enabling large-scale studies.

Integrated scRNA-Seq Identifies Human Postnatal Thymus Seeding Progenitors and Regulatory Dynamics of Differentiating Immature Thymocytes.

During postnatal life, thymopoiesis depends on the continuous colonization of the thymus by bone-marrow-derived hematopoietic progenitors that migrate through the bloodstream. The current understanding of the nature of thymic immigrants is largely based on data from pre-clinical models. Here, we employed single-cell RNA sequencing (scRNA-seq) to examine the immature postnatal thymocyte population in humans. Integration of bone marrow and peripheral blood precursor datasets identified two putative thymus seeding progenitors that varied in expression of CD7; CD10; and the homing receptors CCR7, CCR9, and ITGB7. Whereas both precursors supported T cell development, only one contributed to intrathymic dendritic cell (DC) differentiation, predominantly of plasmacytoid dendritic cells. Trajectory inference delineated the transcriptional dynamics underlying early human T lineage development, enabling prediction of transcription factor (TF) modules that drive stage-specific steps of human T cell development. This comprehensive dataset defines the expression signature of immature human thymocytes and provides a resource for the further study of human thymopoiesis.

Decoding human fetal liver haematopoiesis.

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells and multipotent progenitors (HSC/MPPs) but remains poorly defined in humans. Here, using single-cell transcriptome profiling of approximately 140,000 liver and 74,000 skin, kidney and yolk sac cells, we identify the repertoire of human blood and immune cells during development. We infer differentiation trajectories from HSC/MPPs and evaluate the influence of the tissue microenvironment on blood and immune cell development. We reveal physiological erythropoiesis in fetal skin and the presence of mast cells, natural killer and innate lymphoid cell precursors in the yolk sac. We demonstrate a shift in the haemopoietic composition of fetal liver during gestation away from being predominantly erythroid, accompanied by a parallel change in differentiation potential of HSC/MPPs, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a reference for harnessing the therapeutic potential of HSC/MPPs.

Transcriptional Atlas of Intestinal Immune Cells Reveals that Neuropeptide α-CGRP Modulates Group 2 Innate Lymphoid Cell Responses.

Signaling abnormalities in immune responses in the small intestine can trigger chronic type 2 inflammation involving interaction of multiple immune cell types. To systematically characterize this response, we analyzed 58,067 immune cells from the mouse small intestine by single-cell RNA sequencing (scRNA-seq) at steady state and after induction of a type 2 inflammatory reaction to ovalbumin (OVA). Computational analysis revealed broad shifts in both cell-type composition and cell programs in response to the inflammation, especially in group 2 innate lymphoid cells (ILC2s). Inflammation induced the expression of exon 5 of Calca, which encodes the alpha-calcitonin gene-related peptide (α-CGRP), in intestinal KLRG1+ ILC2s. α-CGRP antagonized KLRG1+ ILC2s proliferation but promoted IL-5 expression. Genetic perturbation of α-CGRP increased the proportion of intestinal KLRG1+ ILC2s. Our work highlights a model where α-CGRP-mediated neuronal signaling is critical for suppressing ILC2 expansion and maintaining homeostasis of the type 2 immune machinery.

A Cellular Taxonomy of the Bone Marrow Stroma in Homeostasis and Leukemia.

Stroma is a poorly defined non-parenchymal component of virtually every organ with key roles in organ development, homeostasis, and repair. Studies of the bone marrow stroma have defined individual populations in the stem cell niche regulating hematopoietic regeneration and capable of initiating leukemia. Here, we use single-cell RNA sequencing (scRNA-seq) to define a cellular taxonomy of the mouse bone marrow stroma and its perturbation by malignancy. We identified seventeen stromal subsets expressing distinct hematopoietic regulatory genes spanning new fibroblastic and osteoblastic subpopulations including distinct osteoblast differentiation trajectories. Emerging acute myeloid leukemia impaired mesenchymal osteogenic differentiation and reduced regulatory molecules necessary for normal hematopoiesis. These data suggest that tissue stroma responds to malignant cells by disadvantaging normal parenchymal cells. Our taxonomy of the stromal compartment provides a comprehensive bone marrow cell census and experimental support for cancer cell crosstalk with specific stromal elements to impair normal tissue function and thereby enable emergent cancer.

Optimal-Transport Analysis of Single-Cell Gene Expression Identifies Developmental Trajectories in Reprogramming.

Understanding the molecular programs that guide differentiation during development is a major challenge. Here, we introduce Waddington-OT, an approach for studying developmental time courses to infer ancestor-descendant fates and model the regulatory programs that underlie them. We apply the method to reconstruct the landscape of reprogramming from 315,000 single-cell RNA sequencing (scRNA-seq) profiles, collected at half-day intervals across 18 days. The results reveal a wider range of developmental programs than previously characterized. Cells gradually adopt either a terminal stromal state or a mesenchymal-to-epithelial transition state. The latter gives rise to populations related to pluripotent, extra-embryonic, and neural cells, with each harboring multiple finer subpopulations. The analysis predicts transcription factors and paracrine signals that affect fates and experiments validate that the TF Obox6 and the cytokine GDF9 enhance reprogramming efficiency. Our approach sheds light on the process and outcome of reprogramming and provides a framework applicable to diverse temporal processes in biology.

Genetic determinants of co-accessible chromatin regions in activated T cells across humans.

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4+ T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

Hill kinetics as a noise filter: the role of transcription factor autoregulation in gene cascades.

An intuition based on deterministic models of chemical kinetics is that population heterogeneity of transcription factor levels in cells is transmitted unchanged downstream to the target genes. We use a stochastic model of a two-gene cascade with a self-regulating upstream gene to show that, counter to the intuition, there is no simple mapping (bimodal to bimodal, unimodal to unimodal) between the shapes of the distributions of transcription factor numbers and target protein numbers in cells. Due to the presence of the two regulations, the system contains two nonlinear transfer functions, defined by the Hill kinetics of transcription factor binding. The transfer function of the regulator can "interfere" with the transfer function of the target, converting the bimodal input into a unimodal output or vice versa. We show that this effect can be predicted by a geometric construction. As an example application of the method, we present a case study of a system of several downstream genes of different sensitivities, controlled by a common transcription factor which also regulates its own transcription. We show that a single regulator can induce qualitatively different patterns (binary or graded) of responses to a signal in different downstream genes, depending on whether the sensitivity regions of the transfer functions of the upstream and downstream genes overlap or not. Alternatively, the same model can be interpreted as describing a single downstream gene that has different sensitivities in different cell lines due to mutations. Our model shows, therefore, a possible kinetic mechanism by which different genes can interpret the same biological signal in a different manner.