A high-resolution map of human RNA translation.

Translated small open reading frames (smORFs) can have important regulatory roles and encode microproteins, yet their genome-wide identification has been challenging. We determined the ribosome locations across six primary human cell types and five tissues and detected 7,767 smORFs with translational profiles matching those of known proteins. The human genome was found to contain highly cell-type- and tissue-specific smORFs and a subset that encodes highly conserved amino acid sequences. Changes in the translational efficiency of upstream-encoded smORFs (uORFs) and the corresponding main ORFs predominantly occur in the same direction. Integration with 456 mass-spectrometry datasets confirms the presence of 603 small peptides at the protein level in humans and provides insights into the subcellular localization of these small proteins. This study provides a comprehensive atlas of high-confidence translated smORFs derived from primary human cells and tissues in order to provide a more complete understanding of the translated human genome.

Transient naive reprogramming corrects hiPS cells functionally and epigenetically.

Cells undergo a major epigenome reconfiguration when reprogrammed to human induced pluripotent stem cells (hiPS cells). However, the epigenomes of hiPS cells and human embryonic stem (hES) cells differ significantly, which affects hiPS cell function1-8. These differences include epigenetic memory and aberrations that emerge during reprogramming, for which the mechanisms remain unknown. Here we characterized the persistence and emergence of these epigenetic differences by performing genome-wide DNA methylation profiling throughout primed and naive reprogramming of human somatic cells to hiPS cells. We found that reprogramming-induced epigenetic aberrations emerge midway through primed reprogramming, whereas DNA demethylation begins early in naive reprogramming. Using this knowledge, we developed a transient-naive-treatment (TNT) reprogramming strategy that emulates the embryonic epigenetic reset. We show that the epigenetic memory in hiPS cells is concentrated in cell of origin-dependent repressive chromatin marked by H3K9me3, lamin-B1 and aberrant CpH methylation. TNT reprogramming reconfigures these domains to a hES cell-like state and does not disrupt genomic imprinting. Using an isogenic system, we demonstrate that TNT reprogramming can correct the transposable element overexpression and differential gene expression seen in conventional hiPS cells, and that TNT-reprogrammed hiPS and hES cells show similar differentiation efficiencies. Moreover, TNT reprogramming enhances the differentiation of hiPS cells derived from multiple cell types. Thus, TNT reprogramming corrects epigenetic memory and aberrations, producing hiPS cells that are molecularly and functionally more similar to hES cells than conventional hiPS cells. We foresee TNT reprogramming becoming a new standard for biomedical and therapeutic applications and providing a novel system for studying epigenetic memory.

Modelling human blastocysts by reprogramming fibroblasts into iBlastoids.

Human pluripotent and trophoblast stem cells have been essential alternatives to blastocysts for understanding early human development1-4. However, these simple culture systems lack the complexity to adequately model the spatiotemporal cellular and molecular dynamics that occur during early embryonic development. Here we describe the reprogramming of fibroblasts into in vitro three-dimensional models of the human blastocyst, termed iBlastoids. Characterization of iBlastoids shows that they model the overall architecture of blastocysts, presenting an inner cell mass-like structure, with epiblast- and primitive endoderm-like cells, a blastocoel-like cavity and a trophectoderm-like outer layer of cells. Single-cell transcriptomics further confirmed the presence of epiblast-, primitive endoderm-, and trophectoderm-like cells. Moreover, iBlastoids can give rise to pluripotent and trophoblast stem cells and are capable of modelling, in vitro, several aspects of the early stage of implantation. In summary, we have developed a scalable and tractable system to model human blastocyst biology; we envision that this will facilitate the study of early human development and the effects of gene mutations and toxins during early embryogenesis, as well as aiding in the development of new therapies associated with in vitro fertilization.

Reprogramming roadmap reveals route to human induced trophoblast stem cells.

The reprogramming of human somatic cells to primed or naive induced pluripotent stem cells recapitulates the stages of early embryonic development1-6. The molecular mechanism that underpins these reprogramming processes remains largely unexplored, which impedes our understanding and limits rational improvements to reprogramming protocols. Here, to address these issues, we reconstruct molecular reprogramming trajectories of human dermal fibroblasts using single-cell transcriptomics. This revealed that reprogramming into primed and naive pluripotency follows diverging and distinct trajectories. Moreover, genome-wide analyses of accessible chromatin showed key changes in the regulatory elements of core pluripotency genes, and orchestrated global changes in chromatin accessibility over time. Integrated analysis of these datasets revealed a role for transcription factors associated with the trophectoderm lineage, and the existence of a subpopulation of cells that enter a trophectoderm-like state during reprogramming. Furthermore, this trophectoderm-like state could be captured, which enabled the derivation of induced trophoblast stem cells. Induced trophoblast stem cells are molecularly and functionally similar to trophoblast stem cells derived from human blastocysts or first-trimester placentas7. Our results provide a high-resolution roadmap for the transcription-factor-mediated reprogramming of human somatic cells, indicate a role for the trophectoderm-lineage-specific regulatory program during this process, and facilitate the direct reprogramming of somatic cells into induced trophoblast stem cells.

A single-cell atlas identifies pretreatment features of primary imatinib resistance in chronic myeloid leukemia.

Primary resistance to tyrosine kinase inhibitors (TKIs) is a significant barrier to optimal outcomes in chronic myeloid leukemia (CML), but factors contributing to response heterogeneity remain unclear. Using single-cell RNA (scRNA) sequencing, we identified 8 statistically significant features in pretreatment bone marrow, which correlated with either sensitivity (major molecular response or MMR) or extreme resistance to imatinib (eventual blast crisis [BC] transformation). Employing machine-learning, we identified leukemic stem cell (LSC) and natural killer (NK) cell gene expression profiles predicting imatinib response with >80% accuracy, including no false positives for predicting BC. A canonical erythroid-specifying (TAL1/KLF1/GATA1) regulon was a hallmark of LSCs from patients with MMR and was associated with erythroid progenitor [ERP] expansion in vivo (P < .05), and a 2- to 10-fold (6.3-fold in group A vs 1.09-fold in group C) erythroid over myeloid bias in vitro. Notably, ERPs demonstrated exquisite TKI sensitivity compared with myeloid progenitors (P < .001). These LSC features were lost with progressive resistance, and MYC- and IRF1-driven inflammatory regulons were evident in patients who progressed to transformation. Patients with MMR also exhibited a 56-fold expansion (P < .01) of a normally rare subset of hyperfunctional adaptive-like NK cells, which diminished with progressive resistance, whereas patients destined for BC accumulated inhibitory NKG2A+ NK cells favoring NK cell tolerance. Finally, we developed antibody panels to validate our scRNA-seq findings. These panels may be useful for prospective studies of primary resistance, and in assessing the contribution of predetermined vs acquired factors in TKI response heterogeneity.

Striatin plays a major role in angiotensin II-induced cardiomyocyte and cardiac hypertrophy in mice in vivo.

The three striatins (STRN, STRN3, STRN4) form the core of STRiatin-Interacting Phosphatase and Kinase (STRIPAK) complexes. These place protein phosphatase 2A (PP2A) in proximity to protein kinases thereby restraining kinase activity and regulating key cellular processes. Our aim was to establish if striatins play a significant role in cardiac remodelling associated with cardiac hypertrophy and heart failure. All striatins were expressed in control human hearts, with up-regulation of STRN and STRN3 in failing hearts. We used mice with global heterozygote gene deletion to assess the roles of STRN and STRN3 in cardiac remodelling induced by angiotensin II (AngII; 7 days). Using echocardiography, we detected no differences in baseline cardiac function or dimensions in STRN+/- or STRN3+/- male mice (8 weeks) compared with wild-type littermates. Heterozygous gene deletion did not affect cardiac function in mice treated with AngII, but the increase in left ventricle mass induced by AngII was inhibited in STRN+/- (but not STRN3+/-) mice. Histological staining indicated that cardiomyocyte hypertrophy was inhibited. To assess the role of STRN in cardiomyocytes, we converted the STRN knockout line for inducible cardiomyocyte-specific gene deletion. There was no effect of cardiomyocyte STRN knockout on cardiac function or dimensions, but the increase in left ventricle mass induced by AngII was inhibited. This resulted from inhibition of cardiomyocyte hypertrophy and cardiac fibrosis. The data indicate that cardiomyocyte striatin is required for early remodelling of the heart by AngII and identify the striatin-based STRIPAK system as a signalling paradigm in the development of pathological cardiac hypertrophy.

CLIPreg: constructing translational regulatory networks from CLIP-, Ribo- and RNA-seq.

MOTIVATION: The creation and analysis of gene regulatory networks have been the focus of bioinformatics research and underpins much of what is known about gene regulation. However, as a result of a bias in the availability of data types that are collected, the vast majority of gene regulatory network resources and tools have focused on either transcriptional regulation or protein-protein interactions. This has left other areas of regulation, for instance, translational regulation, vastly underrepresented despite them having been shown to play a critical role in both health and disease. RESULTS: In order to address this, we have developed CLIPreg, a package that integrates RNA, Ribo and CLIP- sequencing data in order to construct translational regulatory networks coordinated by RNA-binding proteins and micro-RNAs. This is the first tool of its type to be created, allowing for detailed investigation into a previously unseen layer of regulation. AVAILABILITY AND IMPLEMENTATION: CLIPreg is available at https://github.com/SGDDNB/CLIPreg. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Discovering microproteins: making the most of ribosome profiling data.

Building a reference set of protein-coding open reading frames (ORFs) has revolutionized biological process discovery and understanding. Traditionally, gene models have been confirmed using cDNA sequencing and encoded translated regions inferred using sequence-based detection of start and stop combinations longer than 100 amino-acids to prevent false positives. This has led to small ORFs (smORFs) and their encoded proteins left un-annotated. Ribo-seq allows deciphering translated regions from untranslated irrespective of the length. In this review, we describe the power of Ribo-seq data in detection of smORFs while discussing the major challenge posed by data-quality, -depth and -sparseness in identifying the start and end of smORF translation. In particular, we outline smORF cataloguing efforts in humans and the large differences that have arisen due to variation in data, methods and assumptions. Although current versions of smORF reference sets can already be used as a powerful tool for hypothesis generation, we recommend that future editions should consider these data limitations and adopt unified processing for the community to establish a canonical catalogue of translated smORFs.

An atlas of human long non-coding RNAs with accurate 5' ends.

Long non-coding RNAs (lncRNAs) are largely heterogeneous and functionally uncharacterized. Here, using FANTOM5 cap analysis of gene expression (CAGE) data, we integrate multiple transcript collections to generate a comprehensive atlas of 27,919 human lncRNA genes with high-confidence 5' ends and expression profiles across 1,829 samples from the major human primary cell types and tissues. Genomic and epigenomic classification of these lncRNAs reveals that most intergenic lncRNAs originate from enhancers rather than from promoters. Incorporating genetic and expression data, we show that lncRNAs overlapping trait-associated single nucleotide polymorphisms are specifically expressed in cell types relevant to the traits, implicating these lncRNAs in multiple diseases. We further demonstrate that lncRNAs overlapping expression quantitative trait loci (eQTL)-associated single nucleotide polymorphisms of messenger RNAs are co-expressed with the corresponding messenger RNAs, suggesting their potential roles in transcriptional regulation. Combining these findings with conservation data, we identify 19,175 potentially functional lncRNAs in the human genome.

IL-11 is a crucial determinant of cardiovascular fibrosis.

Fibrosis is a common pathology in cardiovascular disease. In the heart, fibrosis causes mechanical and electrical dysfunction and in the kidney, it predicts the onset of renal failure. Transforming growth factor ß1 (TGFß1) is the principal pro-fibrotic factor, but its inhibition is associated with side effects due to its pleiotropic roles. We hypothesized that downstream effectors of TGFß1 in fibroblasts could be attractive therapeutic targets and lack upstream toxicity. Here we show, using integrated imaging-genomics analyses of primary human fibroblasts, that upregulation of interleukin-11 (IL-11) is the dominant transcriptional response to TGFß1 exposure and required for its pro-fibrotic effect. IL-11 and its receptor (IL11RA) are expressed specifically in fibroblasts, in which they drive non-canonical, ERK-dependent autocrine signalling that is required for fibrogenic protein synthesis. In mice, fibroblast-specific Il11 transgene expression or Il-11 injection causes heart and kidney fibrosis and organ failure, whereas genetic deletion of Il11ra1 protects against disease. Therefore, inhibition of IL-11 prevents fibroblast activation across organs and species in response to a range of important pro-fibrotic stimuli. These results reveal a central role of IL-11 in fibrosis and we propose that inhibition of IL-11 is a potential therapeutic strategy to treat fibrotic diseases.

PKN2 deficiency leads both to prenatal 'congenital' cardiomyopathy and defective angiotensin II stress responses.

The protein kinase PKN2 is required for embryonic development and PKN2 knockout mice die as a result of failure in the expansion of mesoderm, cardiac development and neural tube closure. In the adult, cardiomyocyte PKN2 and PKN1 (in combination) are required for cardiac adaptation to pressure-overload. The specific role of PKN2 in contractile cardiomyocytes during development and its role in the adult heart remain to be fully established. We used mice with cardiomyocyte-directed knockout of PKN2 or global PKN2 haploinsufficiency to assess cardiac development and function using high resolution episcopic microscopy, MRI, micro-CT and echocardiography. Biochemical and histological changes were also assessed. Cardiomyocyte-directed PKN2 knockout embryos displayed striking abnormalities in the compact myocardium, with frequent myocardial clefts and diverticula, ventricular septal defects and abnormal heart shape. The sub-Mendelian homozygous knockout survivors developed cardiac failure. RNASeq data showed up-regulation of PKN2 in patients with dilated cardiomyopathy, suggesting an involvement in adult heart disease. Given the rarity of homozygous survivors with cardiomyocyte-specific deletion of PKN2, the requirement for PKN2 in adult mice was explored using the constitutive heterozygous PKN2 knockout. Cardiac hypertrophy resulting from hypertension induced by angiotensin II was reduced in these haploinsufficient PKN2 mice relative to wild-type littermates, with suppression of cardiomyocyte hypertrophy and cardiac fibrosis. It is concluded that cardiomyocyte PKN2 is essential for heart development and the formation of compact myocardium and is also required for cardiac hypertrophy in hypertension. Thus, PKN signalling may offer therapeutic options for managing congenital and adult heart diseases.

Cardiomyocyte BRAF and type 1 RAF inhibitors promote cardiomyocyte and cardiac hypertrophy in mice in vivo.

The extracellular signal-regulated kinase 1/2 (ERK1/2) cascade promotes cardiomyocyte hypertrophy and is cardioprotective, with the three RAF kinases forming a node for signal integration. Our aims were to determine if BRAF is relevant for human heart failure, whether BRAF promotes cardiomyocyte hypertrophy, and if Type 1 RAF inhibitors developed for cancer (that paradoxically activate ERK1/2 at low concentrations: the 'RAF paradox') may have the same effect. BRAF was up-regulated in heart samples from patients with heart failure compared with normal controls. We assessed the effects of activated BRAF in the heart using mice with tamoxifen-activated Cre for cardiomyocyte-specific knock-in of the activating V600E mutation into the endogenous gene. We used echocardiography to measure cardiac dimensions/function. Cardiomyocyte BRAFV600E induced cardiac hypertrophy within 10 d, resulting in increased ejection fraction and fractional shortening over 6 weeks. This was associated with increased cardiomyocyte size without significant fibrosis, consistent with compensated hypertrophy. The experimental Type 1 RAF inhibitor, SB590885, and/or encorafenib (a RAF inhibitor used clinically) increased ERK1/2 phosphorylation in cardiomyocytes, and promoted hypertrophy, consistent with a 'RAF paradox' effect. Both promoted cardiac hypertrophy in mouse hearts in vivo, with increased cardiomyocyte size and no overt fibrosis. In conclusion, BRAF potentially plays an important role in human failing hearts, activation of BRAF is sufficient to induce hypertrophy, and Type 1 RAF inhibitors promote hypertrophy via the 'RAF paradox'. Cardiac hypertrophy resulting from these interventions was not associated with pathological features, suggesting that Type 1 RAF inhibitors may be useful to boost cardiomyocyte function.

ShinyCell: simple and sharable visualization of single-cell gene expression data.

MOTIVATION: As the generation of complex single-cell RNA sequencing datasets becomes more commonplace it is the responsibility of researchers to provide access to these data in a way that can be easily explored and shared. Whilst it is often the case that data is deposited for future bioinformatic analysis many studies do not release their data in a way that is easy to explore by non-computational researchers. RESULTS: In order to help address this we have developed ShinyCell, an R package that converts single-cell RNA sequencing datasets into explorable and shareable interactive interfaces. These interfaces can be easily customized in order to maximize their usability and can be easily uploaded to online platforms to facilitate wider access to published data. AVAILABILITY AND IMPLEMENTATION: ShinyCell is available at https://github.com/SGDDNB/ShinyCell and https://figshare.com/projects/ShinyCell/100439. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

GeneSwitches: ordering gene expression and functional events in single-cell experiments.

SUMMARY: Emerging single-cell RNA-sequencing data technologies has made it possible to capture and assess the gene expression of individual cells. Based on the similarity of gene expression profiles, many tools have been developed to generate an in silico ordering of cells in the form of pseudo-time trajectories. However, these tools do not provide a means to find the ordering of critical gene expression changes over pseudo-time. We present GeneSwitches, a tool that takes any single-cell pseudo-time trajectory and determines the precise order of gene expression and functional-event changes over time. GeneSwitches uses a statistical framework based on logistic regression to identify the order in which genes are either switched on or off along pseudo-time. With this information, users can identify the order in which surface markers appear, investigate how functional ontologies are gained or lost over time and compare the ordering of switching genes from two related pseudo-temporal processes. AVAILABILITY: GeneSwitches is available at https://geneswitches.ddnetbio.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Whole-Genome Sequencing of Finnish Type 1 Diabetic Siblings Discordant for Kidney Disease Reveals DNA Variants associated with Diabetic Nephropathy.

BACKGROUND: Several genetic susceptibility loci associated with diabetic nephropathy have been documented, but no causative variants implying novel pathogenetic mechanisms have been elucidated. METHODS: We carried out whole-genome sequencing of a discovery cohort of Finnish siblings with type 1 diabetes who were discordant for the presence (case) or absence (control) of diabetic nephropathy. Controls had diabetes without complications for 15-37 years. We analyzed and annotated variants at genome, gene, and single-nucleotide variant levels. We then replicated the associated variants, genes, and regions in a replication cohort from the Finnish Diabetic Nephropathy study that included 3531 unrelated Finns with type 1 diabetes. RESULTS: We observed protein-altering variants and an enrichment of variants in regions associated with the presence or absence of diabetic nephropathy. The replication cohort confirmed variants in both regulatory and protein-coding regions. We also observed that diabetic nephropathy-associated variants, when clustered at the gene level, are enriched in a core protein-interaction network representing proteins essential for podocyte function. These genes include protein kinases (protein kinase C isoforms ε and ι) and protein tyrosine kinase 2. CONCLUSIONS: Our comprehensive analysis of a diabetic nephropathy cohort of siblings with type 1 diabetes who were discordant for kidney disease points to variants and genes that are potentially causative or protective for diabetic nephropathy. This includes variants in two isoforms of the protein kinase C family not previously linked to diabetic nephropathy, adding support to previous hypotheses that the protein kinase C family members play a role in diabetic nephropathy and might be attractive therapeutic targets.

Widespread Translational Control of Fibrosis in the Human Heart by RNA-Binding Proteins.

BACKGROUND: Fibrosis is a common pathology in many cardiac disorders and is driven by the activation of resident fibroblasts. The global posttranscriptional mechanisms underlying fibroblast-to-myofibroblast conversion in the heart have not been explored. METHODS: Genome-wide changes of RNA transcription and translation during human cardiac fibroblast activation were monitored with RNA sequencing and ribosome profiling. We then used RNA-binding protein-based analyses to identify translational regulators of fibrogenic genes. The integration with cardiac ribosome occupancy levels of 30 dilated cardiomyopathy patients demonstrates that these posttranscriptional mechanisms are also active in the diseased fibrotic human heart. RESULTS: We generated nucleotide-resolution translatome data during the transforming growth factor ß1-driven cellular transition of human cardiac fibroblasts to myofibroblasts. This identified dynamic changes of RNA transcription and translation at several time points during the fibrotic response, revealing transient and early-responder genes. Remarkably, about one-third of all changes in gene expression in activated fibroblasts are subject to translational regulation, and dynamic variation in ribosome occupancy affects protein abundance independent of RNA levels. Targets of RNA-binding proteins were strongly enriched in posttranscriptionally regulated genes, suggesting genes such as MBNL2 can act as translational activators or repressors. Ribosome occupancy in the hearts of patients with dilated cardiomyopathy suggested the same posttranscriptional regulatory network was underlying cardiac fibrosis. Key network hubs include RNA-binding proteins such as Pumilio RNA binding family member 2 (PUM2) and Quaking (QKI) that work in concert to regulate the translation of target transcripts in human diseased hearts. Furthermore, silencing of both PUM2 and QKI inhibits the transition of fibroblasts toward profibrotic myofibroblasts in response to transforming growth factor ß1. CONCLUSIONS: We reveal widespread translational effects of transforming growth factor ß1 and define novel posttranscriptional regulatory networks that control the fibroblast-to-myofibroblast transition. These networks are active in human heart disease, and silencing of hub genes limits fibroblast activation. Our findings show the central importance of translational control in fibrosis and highlight novel pathogenic mechanisms in heart failure.

An atlas of active enhancers across human cell types and tissues.

Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.

Update of the FANTOM web resource: high resolution transcriptome of diverse cell types in mammals.

Upon the first publication of the fifth iteration of the Functional Annotation of Mammalian Genomes collaborative project, FANTOM5, we gathered a series of primary data and database systems into the FANTOM web resource (http://fantom.gsc.riken.jp) to facilitate researchers to explore transcriptional regulation and cellular states. In the course of the collaboration, primary data and analysis results have been expanded, and functionalities of the database systems enhanced. We believe that our data and web systems are invaluable resources, and we think the scientific community will benefit for this recent update to deepen their understanding of mammalian cellular organization. We introduce the contents of FANTOM5 here, report recent updates in the web resource and provide future perspectives.

RNA polymerase II primes Polycomb-repressed developmental genes throughout terminal neuronal differentiation.

Polycomb repression in mouse embryonic stem cells (ESCs) is tightly associated with promoter co-occupancy of RNA polymerase II (RNAPII) which is thought to prime genes for activation during early development. However, it is unknown whether RNAPII poising is a general feature of Polycomb repression, or is lost during differentiation. Here, we map the genome-wide occupancy of RNAPII and Polycomb from pluripotent ESCs to non-dividing functional dopaminergic neurons. We find that poised RNAPII complexes are ubiquitously present at Polycomb-repressed genes at all stages of neuronal differentiation. We observe both loss and acquisition of RNAPII and Polycomb at specific groups of genes reflecting their silencing or activation. Strikingly, RNAPII remains poised at transcription factor genes which are silenced in neurons through Polycomb repression, and have major roles in specifying other, non-neuronal lineages. We conclude that RNAPII poising is intrinsically associated with Polycomb repression throughout differentiation. Our work suggests that the tight interplay between RNAPII poising and Polycomb repression not only instructs promoter state transitions, but also may enable promoter plasticity in differentiated cells.