Cell Types of the Human Retina and Its Organoids at Single-Cell Resolution.

Human organoids recapitulating the cell-type diversity and function of their target organ are valuable for basic and translational research. We developed light-sensitive human retinal organoids with multiple nuclear and synaptic layers and functional synapses. We sequenced the RNA of 285,441 single cells from these organoids at seven developmental time points and from the periphery, fovea, pigment epithelium and choroid of light-responsive adult human retinas, and performed histochemistry. Cell types in organoids matured in vitro to a stable "developed" state at a rate similar to human retina development in vivo. Transcriptomes of organoid cell types converged toward the transcriptomes of adult peripheral retinal cell types. Expression of disease-associated genes was cell-type-specific in adult retina, and cell-type specificity was retained in organoids. We implicate unexpected cell types in diseases such as macular degeneration. This resource identifies cellular targets for studying disease mechanisms in organoids and for targeted repair in human retinas.

Self-organization and symmetry breaking in intestinal organoid development.

Intestinal organoids are complex three-dimensional structures that mimic the cell-type composition and tissue organization of the intestine by recapitulating the self-organizing ability of cell populations derived from a single intestinal stem cell. Crucial in this process is a first symmetry-breaking event, in which only a fraction of identical cells in a symmetrical sphere differentiate into Paneth cells, which generate the stem-cell niche and lead to asymmetric structures such as the crypts and villi. Here we combine single-cell quantitative genomic and imaging approaches to characterize the development of intestinal organoids from single cells. We show that their development follows a regeneration process that is driven by transient activation of the transcriptional regulator YAP1. Cell-to-cell variability in YAP1, emerging in symmetrical spheres, initiates Notch and DLL1 activation, and drives the symmetry-breaking event and formation of the first Paneth cell. Our findings reveal how single cells exposed to a uniform growth-promoting environment have the intrinsic ability to generate emergent, self-organized behaviour that results in the formation of complex multicellular asymmetric structures.

Loss of Hepatic Leucine-Rich Repeat-Containing G-Protein Coupled Receptors 4 and 5 Promotes Nonalcoholic Fatty Liver Disease.

The roof plate-specific spondin-leucine-rich repeat-containing G-protein coupled receptor 4/5 (LGR4/5)-zinc and ring finger 3 (ZNRF3)/ring finger protein 43 (RNF43) module is a master regulator of hepatic Wnt/ß-catenin signaling and metabolic zonation. However, its impact on nonalcoholic fatty liver disease (NAFLD) remains unclear. The current study investigated whether hepatic epithelial cell-specific loss of the Wnt/ß-catenin modulator Lgr4/5 promoted NAFLD. The 3- and 6-month-old mice with hepatic epithelial cell-specific deletion of both receptors Lgr4/5 (Lgr4/5dLKO) were compared with control mice fed with normal diet (ND) or high-fat diet (HFD). Six-month-old HFD-fed Lgr4/5dLKO mice developed hepatic steatosis and fibrosis but the control mice did not. Serum cholesterol-high-density lipoprotein and total cholesterol levels in 3- and 6-month-old HFD-fed Lgr4/5dLKO mice were decreased compared with those in control mice. An ex vivo primary hepatocyte culture assay and a comprehensive bile acid (BA) characterization in liver, plasma, bile, and feces demonstrated that ND-fed Lgr4/5dLKO mice had impaired BA secretion, predisposing them to develop cholestatic characteristics. Lipidome and RNA-sequencing analyses demonstrated severe alterations in several lipid species and pathways controlling lipid metabolism in the livers of Lgr4/5dLKO mice. In conclusion, loss of hepatic Wnt/ß-catenin activity by Lgr4/5 deletion led to loss of BA secretion, cholestatic features, altered lipid homeostasis, and deregulation of lipoprotein pathways. Both BA and intrinsic lipid alterations contributed to the onset of NAFLD.

A single-cell atlas of the airway epithelium reveals the CFTR-rich pulmonary ionocyte.

The functions of epithelial tissues are dictated by the types, abundance and distribution of the differentiated cells they contain. Attempts to restore tissue function after damage require knowledge of how physiological tasks are distributed among cell types, and how cell states vary between homeostasis, injury-repair and disease. In the conducting airway, a heterogeneous basal cell population gives rise to specialized luminal cells that perform mucociliary clearance1. Here we perform single-cell profiling of human bronchial epithelial cells and mouse tracheal epithelial cells to obtain a comprehensive census of cell types in the conducting airway and their behaviour in homeostasis and regeneration. Our analysis reveals cell states that represent known and novel cell populations, delineates their heterogeneity and identifies distinct differentiation trajectories during homeostasis and tissue repair. Finally, we identified a novel, rare cell type that we call the 'pulmonary ionocyte', which co-expresses FOXI1, multiple subunits of the vacuolar-type H+-ATPase (V-ATPase) and CFTR, the gene that is mutated in cystic fibrosis. Using immunofluorescence, modulation of signalling pathways and electrophysiology, we show that Notch signalling is necessary and FOXI1 expression is sufficient to drive the production of the pulmonary ionocyte, and that the pulmonary ionocyte is a major source of CFTR activity in the conducting airway epithelium.

Molecular characterization of chronic cutaneous wounds reveals subregion- and wound type-specific differential gene expression.

A limited understanding of the pathology underlying chronic wounds has hindered the development of effective diagnostic markers and pharmaceutical interventions. This study aimed to elucidate the molecular composition of various common chronic ulcer types to facilitate drug discovery strategies. We conducted a comprehensive analysis of leg ulcers (LUs), encompassing venous and arterial ulcers, foot ulcers (FUs), pressure ulcers (PUs), and compared them with surgical wound healing complications (WHCs). To explore the pathophysiological mechanisms and identify similarities or differences within wounds, we dissected wounds into distinct subregions, including the wound bed, border, and peri-wound areas, and compared them against intact skin. By correlating histopathology, RNA sequencing (RNA-Seq), and immunohistochemistry (IHC), we identified unique genes, pathways, and cell type abundance patterns in each wound type and subregion. These correlations aim to aid clinicians in selecting targeted treatment options and informing the design of future preclinical and clinical studies in wound healing. Notably, specific genes, such as PITX1 and UPP1, exhibited exclusive upregulation in LUs and FUs, potentially offering significant benefits to specialists in limb preservation and clinical treatment decisions. In contrast, comparisons between different wound subregions, regardless of wound type, revealed distinct expression profiles. The pleiotropic chemokine-like ligand GPR15L (C10orf99) and transmembrane serine proteases TMPRSS11A/D were significantly upregulated in wound border subregions. Interestingly, WHCs exhibited a nearly identical transcriptome to PUs, indicating clinical relevance. Histological examination revealed blood vessel occlusions with impaired angiogenesis in chronic wounds, alongside elevated expression of genes and immunoreactive markers related to blood vessel and lymphatic epithelial cells in wound bed subregions. Additionally, inflammatory and epithelial markers indicated heightened inflammatory responses in wound bed and border subregions and reduced wound bed epithelialization. In summary, chronic wounds from diverse anatomical sites share common aspects of wound pathophysiology but also exhibit distinct molecular differences. These unique molecular characteristics present promising opportunities for drug discovery and treatment, particularly for patients suffering from chronic wounds. The identified diagnostic markers hold the potential to enhance preclinical and clinical trials in the field of wound healing.

BET bromodomain inhibitors regulate keratinocyte plasticity.

Although most acute skin wounds heal rapidly, non-healing skin ulcers represent an increasing and substantial unmet medical need that urgently requires effective therapeutics. Keratinocytes resurface wounds to re-establish the epidermal barrier by transitioning to an activated, migratory state, but this ability is lost in dysfunctional chronic wounds. Small-molecule regulators of keratinocyte plasticity with the potential to reverse keratinocyte malfunction in situ could offer a novel therapeutic approach in skin wound healing. Utilizing high-throughput phenotypic screening of primary keratinocytes, we identify such small molecules, including bromodomain and extra-terminal domain (BET) protein family inhibitors (BETi). BETi induce a sustained activated, migratory state in keratinocytes in vitro, increase activation markers in human epidermis ex vivo and enhance skin wound healing in vivo. Our findings suggest potential clinical utility of BETi in promoting keratinocyte re-epithelialization of skin wounds. Importantly, this novel property of BETi is exclusively observed after transient low-dose exposure, revealing new potential for this compound class.

Systematic identification of novel cancer genes through analysis of deep shRNA perturbation screens.

Systematic perturbation screens provide comprehensive resources for the elucidation of cancer driver genes. The perturbation of many genes in relatively few cell lines in such functional screens necessitates the development of specialized computational tools with sufficient statistical power. Here we developed APSiC (Analysis of Perturbation Screens for identifying novel Cancer genes) to identify genetic drivers and effectors in perturbation screens even with few samples. Applying APSiC to the shRNA screen Project DRIVE, APSiC identified well-known and novel putative mutational and amplified cancer genes across all cancer types and in specific cancer types. Additionally, APSiC discovered tumor-promoting and tumor-suppressive effectors, respectively, for individual cancer types, including genes involved in cell cycle control, Wnt/ß-catenin and hippo signalling pathways. We functionally demonstrated that LRRC4B, a putative novel tumor-suppressive effector, suppresses proliferation by delaying cell cycle and modulates apoptosis in breast cancer. We demonstrate APSiC is a robust statistical framework for discovery of novel cancer genes through analysis of large-scale perturbation screens. The analysis of DRIVE using APSiC is provided as a web portal and represents a valuable resource for the discovery of novel cancer genes.

CPSF3-dependent pre-mRNA processing as a druggable node in AML and Ewing's sarcoma.

The post-genomic era has seen many advances in our understanding of cancer pathways, yet resistance and tumor heterogeneity necessitate multiple approaches to target even monogenic tumors. Here, we combine phenotypic screening with chemical genetics to identify pre-messenger RNA endonuclease cleavage and polyadenylation specificity factor 3 (CPSF3) as the target of JTE-607, a small molecule with previously unknown target. We show that CPSF3 represents a synthetic lethal node in a subset of acute myeloid leukemia (AML) and Ewing's sarcoma cancer cell lines. Inhibition of CPSF3 by JTE-607 alters expression of known downstream effectors in AML and Ewing's sarcoma lines, upregulates apoptosis and causes tumor-selective stasis in mouse xenografts. Mechanistically, it prevents the release of newly synthesized pre-mRNAs, resulting in read-through transcription and the formation of DNA-RNA hybrid R-loop structures. This study implicates pre-mRNA processing, and specifically CPSF3, as a druggable target providing an avenue to therapeutic intervention in cancer.

Author Correction: CPSF3-dependent pre-mRNA processing as a druggable node in AML and Ewing's sarcoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Transcriptional profiling of hepatocytes infected with the replicative form of the malaria parasite Plasmodium cynomolgi.

BACKGROUND: The zoonotic simian parasite Plasmodium cynomolgi develops into replicating schizonts and dormant hypnozoites during the infection of hepatocytes and is used as a model organism to study relapsing malaria. The transcriptional profiling of P. cynomolgi liver stages was previously reported and revealed many important biological features of the parasite but left out the host response to malaria infection. METHODS: Previously published RNA sequencing data were used to quantify the expression of host genes in rhesus macaque hepatocytes infected with P. cynomolgi in comparison to either cells from uninfected samples or uninfected bystander cells. RESULTS: Although the dataset could not be used to resolve the transcriptional profile of hypnozoite-infected hepatocytes, it provided a snapshot of the host response to liver stage schizonts at 9-10 day post-infection and identified specific host pathways that are modulated during the exo-erythrocytic stage of P. cynomolgi. CONCLUSIONS: This study constitutes a valuable resource characterizing the hepatocyte response to P. cynomolgi infection and provides a framework to build on future research that aims at understanding hepatocyte-parasite interactions during relapsing malaria infection.

Discovery of a ZIP7 inhibitor from a Notch pathway screen.

The identification of activating mutations in NOTCH1 in 50% of T cell acute lymphoblastic leukemia has generated interest in elucidating how these mutations contribute to oncogenic transformation and in targeting the pathway. A phenotypic screen identified compounds that interfere with trafficking of Notch and induce apoptosis via an endoplasmic reticulum (ER) stress mechanism. Target identification approaches revealed a role for SLC39A7 (ZIP7), a zinc transport family member, in governing Notch trafficking and signaling. Generation and sequencing of a compound-resistant cell line identified a V430E mutation in ZIP7 that confers transferable resistance to the compound NVS-ZP7-4. NVS-ZP7-4 altered zinc in the ER, and an analog of the compound photoaffinity labeled ZIP7 in cells, suggesting a direct interaction between the compound and ZIP7. NVS-ZP7-4 is the first reported chemical tool to probe the impact of modulating ER zinc levels and investigate ZIP7 as a novel druggable node in the Notch pathway.

mTORC1 signaling suppresses Wnt/ß-catenin signaling through DVL-dependent regulation of Wnt receptor FZD level.

Wnt/ß-catenin signaling plays pivotal roles in cell proliferation and tissue homeostasis by maintaining somatic stem cell functions. The mammalian target of rapamycin (mTOR) signaling functions as an integrative rheostat that orchestrates various cellular and metabolic activities that shape tissue homeostasis. Whether these two fundamental signaling pathways couple to exert physiological functions still remains mysterious. Using a genome-wide CRISPR-Cas9 screening, we discover that mTOR complex 1 (mTORC1) signaling suppresses canonical Wnt/ß-catenin signaling. Deficiency in tuberous sclerosis complex 1/2 (TSC1/2), core negative regulators of mTORC1 activity, represses Wnt/ß-catenin target gene expression, which can be rescued by RAD001. Mechanistically, mTORC1 signaling regulates the cell surface level of Wnt receptor Frizzled (FZD) in a Dishevelled (DVL)-dependent manner by influencing the association of DVL and clathrin AP-2 adaptor. Sustained mTORC1 activation impairs Wnt/ß-catenin signaling and causes loss of stemness in intestinal organoids ex vivo and primitive intestinal progenitors in vivo. Wnt/ß-catenin-dependent liver metabolic zonation gene expression program is also down-regulated by mTORC1 activation. Our study provides a paradigm that mTORC1 signaling cell autonomously regulates Wnt/ß-catenin pathway to influence stem cell maintenance.

Genome-wide CRISPR screen for PARKIN regulators reveals transcriptional repression as a determinant of mitophagy.

PARKIN, an E3 ligase mutated in familial Parkinson's disease, promotes mitophagy by ubiquitinating mitochondrial proteins for efficient engagement of the autophagy machinery. Specifically, PARKIN-synthesized ubiquitin chains represent targets for the PINK1 kinase generating phosphoS65-ubiquitin (pUb), which constitutes the mitophagy signal. Physiological regulation of PARKIN abundance, however, and the impact on pUb accumulation are poorly understood. Using cells designed to discover physiological regulators of PARKIN abundance, we performed a pooled genome-wide CRISPR/Cas9 knockout screen. Testing identified genes individually resulted in a list of 53 positive and negative regulators. A transcriptional repressor network including THAP11 was identified and negatively regulates endogenous PARKIN abundance. RNAseq analysis revealed the PARKIN-encoding locus as a prime THAP11 target, and THAP11 CRISPR knockout in multiple cell types enhanced pUb accumulation. Thus, our work demonstrates the critical role of PARKIN abundance, identifies regulating genes, and reveals a link between transcriptional repression and mitophagy, which is also apparent in human induced pluripotent stem cell-derived neurons, a disease-relevant cell type.

Age-associated inflammation inhibits epidermal stem cell function.

Altered stem cell homeostasis is linked to organismal aging. However, the mechanisms involved remain poorly understood. Here we report novel alterations in hair follicle stem cells during skin aging, including increased numbers, decreased function, and an inability to tolerate stress. Performing high-throughput RNA sequencing on aging stem cells, cytokine arrays, and functional assays, we identify an age-associated imbalance in epidermal Jak-Stat signaling that inhibits stem cell function. Collectively, this study reveals a role for the aging epidermis in the disruption of cytokine and stem cell homeostasis, suggesting that stem cell decline during aging may be part of broader tumor-suppressive mechanisms.

The exon quantification pipeline (EQP): a comprehensive approach to the quantification of gene, exon and junction expression from RNA-seq data.

The quantification of transcriptomic features is the basis of the analysis of RNA-seq data. We present an integrated alignment workflow and a simple counting-based approach to derive estimates for gene, exon and exon-exon junction expression. In contrast to previous counting-based approaches, EQP takes into account only reads whose alignment pattern agrees with the splicing pattern of the features of interest. This leads to improved gene expression estimates as well as to the generation of exon counts that allow disambiguating reads between overlapping exons. Unlike other methods that quantify skipped introns, EQP offers a novel way to compute junction counts based on the agreement of the read alignments with the exons on both sides of the junction, thus providing a uniformly derived set of counts. We evaluated the performance of EQP on both simulated and real Illumina RNA-seq data and compared it with other quantification tools. Our results suggest that EQP provides superior gene expression estimates and we illustrate the advantages of EQP's exon and junction counts. The provision of uniformly derived high-quality counts makes EQP an ideal quantification tool for differential expression and differential splicing studies. EQP is freely available for download at https://github.com/Novartis/EQP-cluster.

Integration of multi-omics data of a genome-reduced bacterium: Prevalence of post-transcriptional regulation and its correlation with protein abundances.

We developed a comprehensive resource for the genome-reduced bacterium Mycoplasma pneumoniae comprising 1748 consistently generated '-omics' data sets, and used it to quantify the power of antisense non-coding RNAs (ncRNAs), lysine acetylation, and protein phosphorylation in predicting protein abundance (11%, 24% and 8%, respectively). These factors taken together are four times more predictive of the proteome abundance than of mRNA abundance. In bacteria, post-translational modifications (PTMs) and ncRNA transcription were both found to increase with decreasing genomic GC-content and genome size. Thus, the evolutionary forces constraining genome size and GC-content modify the relative contributions of the different regulatory layers to proteome homeostasis, and impact more genomic and genetic features than previously appreciated. Indeed, these scaling principles will enable us to develop more informed approaches when engineering minimal synthetic genomes.

YAP1 Exerts Its Transcriptional Control via TEAD-Mediated Activation of Enhancers.

YAP1 is a major effector of the Hippo pathway and a well-established oncogene. Elevated YAP1 activity due to mutations in Hippo pathway components or YAP1 amplification is observed in several types of human cancers. Here we investigated its genomic binding landscape in YAP1-activated cancer cells, as well as in non-transformed cells. We demonstrate that TEAD transcription factors mediate YAP1 chromatin-binding genome-wide, further explaining their dominant role as primary mediators of YAP1-transcriptional activity. Moreover, we show that YAP1 largely exerts its transcriptional control via distal enhancers that are marked by H3K27 acetylation and that YAP1 is necessary for this chromatin mark at bound enhancers and the activity of the associated genes. This work establishes YAP1-mediated transcriptional regulation at distal enhancers and provides an expanded set of target genes resulting in a fundamental source to study YAP1 function in a normal and cancer setting.

A comprehensive assessment of RNA-seq protocols for degraded and low-quantity samples.

BACKGROUND: RNA-sequencing (RNA-seq) has emerged as one of the most sensitive tool for gene expression analysis. Among the library preparation methods available, the standard poly(A) + enrichment provides a comprehensive, detailed, and accurate view of polyadenylated RNAs. However, on samples of suboptimal quality ribosomal RNA depletion and exon capture methods have recently been reported as better alternatives. METHODS: We compared for the first time three commercial Illumina library preparation kits (TruSeq Stranded mRNA, TruSeq Ribo-Zero rRNA Removal, and TruSeq RNA Access) as representatives of these three different approaches using well-established human reference RNA samples from the MAQC/SEQC consortium on a wide range of input amounts (from 100 ng down to 1 ng) and degradation levels (intact, degraded, and highly degraded). RESULTS: We assessed the accuracy of the generated expression values by comparison to gold standard TaqMan qPCR measurements and gained unprecedented insight into the limits of applicability in terms of input quantity and sample quality of each protocol. We found that each protocol generates highly reproducible results (R 2 > 0.92) on intact RNA samples down to input amounts of 10 ng. For degraded RNA samples, Ribo-Zero showed clear performance advantages over the other two protocols as it generated more accurate and better reproducible gene expression results even at very low input amounts such as 1 ng and 2 ng. For highly degraded RNA samples, RNA Access performed best generating reliable data down to 5 ng input. CONCLUSIONS: We found that the ribosomal RNA depletion protocol from Illumina works very well at amounts far below recommendation and over a good range of intact and degraded material. We also infer that the exome-capture protocol (RNA Access, Illumina) performs better than other methods on highly degraded and low amount samples.

Decatransin, a new natural product inhibiting protein translocation at the Sec61/SecYEG translocon.

A new cyclic decadepsipeptide was isolated from Chaetosphaeria tulasneorum with potent bioactivity on mammalian and yeast cells. Chemogenomic profiling in S. cerevisiae indicated that the Sec61 translocon complex, the machinery for protein translocation and membrane insertion at the endoplasmic reticulum, is the target. The profiles were similar to those of cyclic heptadepsipeptides of a distinct chemotype (including HUN-7293 and cotransin) that had previously been shown to inhibit cotranslational translocation at the mammalian Sec61 translocon. Unbiased, genome-wide mutagenesis followed by full-genome sequencing in both fungal and mammalian cells identified dominant mutations in Sec61p (yeast) or Sec61α1 (mammals) that conferred resistance. Most, but not all, of these mutations affected inhibition by both chemotypes, despite an absence of structural similarity. Biochemical analysis confirmed inhibition of protein translocation into the endoplasmic reticulum of both co- and post-translationally translocated substrates by both chemotypes, demonstrating a mechanism independent of a translating ribosome. Most interestingly, both chemotypes were found to also inhibit SecYEG, the bacterial Sec61 translocon homolog. We suggest 'decatransin' as the name for this new decadepsipeptide translocation inhibitor.

SMN2 splice modulators enhance U1-pre-mRNA association and rescue SMA mice.

Spinal muscular atrophy (SMA), which results from the loss of expression of the survival of motor neuron-1 (SMN1) gene, represents the most common genetic cause of pediatric mortality. A duplicate copy (SMN2) is inefficiently spliced, producing a truncated and unstable protein. We describe herein a potent, orally active, small-molecule enhancer of SMN2 splicing that elevates full-length SMN protein and extends survival in a severe SMA mouse model. We demonstrate that the molecular mechanism of action is via stabilization of the transient double-strand RNA structure formed by the SMN2 pre-mRNA and U1 small nuclear ribonucleic protein (snRNP) complex. The binding affinity of U1 snRNP to the 5' splice site is increased in a sequence-selective manner, discrete from constitutive recognition. This new mechanism demonstrates the feasibility of small molecule-mediated, sequence-selective splice modulation and the potential for leveraging this strategy in other splicing diseases.