Transcriptomic entropy benchmarks stem cell-derived cardiomyocyte maturation against endogenous tissue at single cell level

Protocol for quantifying stem-cell-derived cardiomyocyte maturity using transcriptomic entropy score

The inability to quantify cardiomyocyte (CM) maturation remains a significant barrier to evaluating the effects of ongoing efforts to produce adult-like CMs from pluripotent stem cells (PSCs). Here, we present a protocol to quantify stem-cell-derived CM maturity using a single-cell RNA sequencing-based metric "entropy score." We describe steps for generating an entropy score using customized R code. This tool can be used to quantify maturation levels of PSC-CMs and potentially other cell types. For complete details on the use and execution of this protocol, please refer to Kannan et al.1.

Estimating Transcriptome Diversity and Specialization in Capsicum annuum L

Chili pepper fruits of the genus Capsicum represent excellent experimental models to study the growth, development, and ripening processes in a non-climacteric species at the physiological, biochemical, and molecular levels. Fruit growth, development, and ripening involve a complex, harmonious, and finely controlled regulation of gene expression. The purpose of this study was to estimate the changes in transcriptome diversity and specialization, as well as gene specificities during fruit development in this crop, and to illustrate the advantages of estimating these parameters. To achieve these aims, we programmed and made publicly available an R package. In this study, we applied these methods to a set of 179 RNA-Seq libraries from a factorial experiment that includes 12 different genotypes at various stages of fruit development. We found that the diversity of the transcriptome decreases linearly from the flower to the mature fruit, while its specialization follows a complex and non-linear behavior during this process. Additionally, by defining sets of genes with different degrees of specialization and applying Gene Ontology enrichment analysis, we identified processes, functions, and components that play a central role in particular fruit development stages. In conclusion, the estimation of diversity, specialization, and specificity summarizes the global properties of the transcriptomes, providing insights that are difficult to achieve by other means.

Characterizing heterogeneity along EMT and metabolic axes in colorectal cancer reveals underlying consensus molecular subtype-specific trends

Colorectal cancer (CRC) is highly heterogeneous with variable survival outcomes and therapeutic vulnerabilities. A commonly used classification system in CRC is the Consensus Molecular Subtypes (CMS) based on gene expression patterns. However, how these CMS categories connect to axes of phenotypic plasticity and heterogeneity remains unclear. Here, in our analysis of CMS-specific TCGA data and 101 bulk transcriptomic datasets, we found the epithelial phenotype score to be consistently positively correlated with scores of glycolysis, OXPHOS and FAO pathways, while PD-L1 activity scores positively correlated with mesenchymal phenotype scoring, revealing possible interconnections among plasticity axes. Single-cell RNA-sequencing analysis of patient samples revealed that that CMS2 and CMS3 subtype samples were relatively more epithelial as compared to CMS1 and CMS4. CMS1 revealed two subpopulations: one close to CMS4 (more mesenchymal) and the other closer to CMS2 or CMS3 (more epithelial), indicating a partial EMT-like behavior. Consistent observations were made in single-cell analysis of metabolic axes and PD-L1 activity scores. Together, our results quantify the patterns of two functional interconnected axes of phenotypic heterogeneity - EMT and metabolic reprogramming - in a CMS-specific manner in CRC.

Creating cell-specific computational models of stem cell-derived cardiomyocytes using optical experiments

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have gained traction as a powerful model in cardiac disease and therapeutics research, since iPSCs are self-renewing and can be derived from healthy and diseased patients without invasive surgery. However, current iPSC-CM differentiation methods produce cardiomyocytes with immature, fetal-like electrophysiological phenotypes, and the variety of maturation protocols in the literature results in phenotypic differences between labs. Heterogeneity of iPSC donor genetic backgrounds contributes to additional phenotypic variability. Several mathematical models of iPSC-CM electrophysiology have been developed to help understand the ionic underpinnings of, and to simulate, various cell responses, but these models individually do not capture the phenotypic variability observed in iPSC-CMs. Here, we tackle these limitations by developing a computational pipeline to calibrate cell preparation-specific iPSC-CM electrophysiological parameters. We used the genetic algorithm (GA), a heuristic parameter calibration method, to tune ion channel parameters in a mathematical model of iPSC-CM physiology. To systematically optimize an experimental protocol that generates sufficient data for parameter calibration, we created simulated datasets by applying various protocols to a population of in silico cells with known conductance variations, and we fitted to those datasets. We found that calibrating models to voltage and calcium transient data under 3 varied experimental conditions, including electrical pacing combined with ion channel blockade and changing buffer ion concentrations, improved model parameter estimates and model predictions of unseen channel block responses. This observation held regardless of whether the fitted data were normalized, suggesting that normalized fluorescence recordings, which are more accessible and higher throughput than patch clamp recordings, could sufficiently inform conductance parameters. Therefore, this computational pipeline can be applied to different iPSC-CM preparations to determine cell line-specific ion channel properties and understand the mechanisms behind variability in perturbation responses.

Integrated modeling framework reveals co-regulation of transcription factors, miRNAs and lncRNAs on cardiac developmental dynamics

AIMS

Dissecting complex interactions among transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) are central for understanding heart development and function. Although computational approaches and platforms have been described to infer relationships among regulatory factors and genes, current approaches do not adequately account for how highly diverse, interacting regulators that include noncoding RNAs (ncRNAs) control cardiac gene expression dynamics over time.

METHODS

To overcome this limitation, we devised an integrated framework, cardiac gene regulatory modeling (CGRM) that integrates LogicTRN and regulatory component analysis bioinformatics modeling platforms to infer complex regulatory mechanisms. We then used CGRM to identify and compare the TF-ncRNA gene regulatory networks that govern early- and late-stage cardiomyocytes (CMs) generated by in vitro differentiation of human pluripotent stem cells (hPSC) and ventricular and atrial CMs isolated during in vivo human cardiac development.

RESULTS

Comparisons of in vitro versus in vivo derived CMs revealed conserved regulatory networks among TFs and ncRNAs in early cells that significantly diverged in late staged cells. We report that cardiac genes ("heart targets") expressed in early-stage hPSC-CMs are primarily regulated by MESP1, miR-1, miR-23, lncRNAs NEAT1 and MALAT1, while GATA6, HAND2, miR-200c, NEAT1 and MALAT1 are critical for late hPSC-CMs. The inferred TF-miRNA-lncRNA networks regulating heart development and contraction were similar among early-stage CMs, among individual hPSC-CM datasets and between in vitro and in vivo samples. However, genes related to apoptosis, cell cycle and proliferation, and transmembrane transport showed a high degree of divergence between in vitro and in vivo derived late-stage CMs. Overall, late-, but not early-stage CMs diverged greatly in the expression of "heart target" transcripts and their regulatory mechanisms.

CONCLUSIONS

In conclusion, we find that hPSC-CMs are regulated in a cell autonomous manner during early development that diverges significantly as a function of time when compared to in vivo derived CMs. These findings demonstrate the feasibility of using CGRM to reveal dynamic and complex transcriptional and posttranscriptional regulatory interactions that underlie cell directed versus environment-dependent CM development. These results with in vitro versus in vivo derived CMs thus establish this approach for detailed analyses of heart disease and for the analysis of cell regulatory systems in other biomedical fields.

Immune Infiltration in Tumor and Adjacent Non-Neoplastic Regions Codetermines Patient Clinical Outcomes in Early-Stage Lung Cancer

INTRODUCTION

In recent years, the proportion of patients with NSCLC diagnosed at an early stage has increased continuously.

METHODS

In this study, we analyzed samples and data collected from 119 samples from 67 early stage patients with NSCLC, including 52 pairs of tumor and adjacent non-neoplastic samples, and performed RNA-sequencing analysis with high sequencing depth.

RESULTS

We found that immune-related genes were highly enriched among the differentially expressed genes and observed significantly higher inferred immune infiltration levels in adjacent non-neoplastic samples than in tumor samples. In survival analysis, the infiltration of certain immune cell types in tumor, but not adjacent non-neoplastic, samples were associated with overall patient survival, and excitingly, the differential infiltration between paired samples (tumor minus non-neoplastic) was more prognostic than expression in either non-neoplastic or tumor tissues. We also performed B cell receptor (BCR) and T cell receptor (TCR) repertoire analysis and observed more BCR/TCR clonotypes and increased BCR clonality in tumor than in non-neoplastic samples. Finally, we carefully quantified the fraction of the five histologic subtypes in our adenocarcinoma samples and found that higher histologic pattern complexity was associated with higher immune infiltration and low TCR clonality in the tumor-proximal regions.

CONCLUSIONS

Our results indicated significantly differential immune characteristics between tumor and adjacent non-neoplastic samples and suggested that the two regions provided complementary prognostic values in early-stage NSCLCs.

Adapting to a new environment: postnatal maturation of the human cardiomyocyte

The postnatal mammalian heart undergoes remarkable developmental changes, which are stimulated by the transition from the intrauterine to extrauterine environment. With birth, increased oxygen levels promote metabolic, structural and biophysical maturation of cardiomyocytes, resulting in mature muscle with increased efficiency, contractility and electrical conduction. In this Topical Review article, we highlight key studies that inform our current understanding of human cardiomyocyte maturation. Collectively, these studies suggest that human atrial and ventricular myocytes evolve quickly within the first year but might not reach a fully mature adult phenotype until nearly the first decade of life. However, it is important to note that fetal, neonatal and paediatric cardiac physiology studies are hindered by a number of limitations, including the scarcity of human tissue, small sample size and a heavy reliance on diseased tissue samples, often without age-matched healthy controls. Future developmental studies are warranted to expand our understanding of normal cardiac physiology/pathophysiology and inform age-appropriate treatment strategies for cardiac disease.

Dynamic changes in P300 enhancers and enhancer-promoter contacts control mouse cardiomyocyte maturation

Cardiomyocyte differentiation continues throughout murine gestation and into the postnatal period, driven by temporally regulated expression changes in the transcriptome. The mechanisms that regulate these developmental changes remain incompletely defined. Here, we used cardiomyocyte-specific ChIP-seq of the activate enhancer marker P300 to identify 54,920 cardiomyocyte enhancers at seven stages of murine heart development. These data were matched to cardiomyocyte gene expression profiles at the same stages and to Hi-C and H3K27ac HiChIP chromatin conformation data at fetal, neonatal, and adult stages. Regions with dynamic P300 occupancy exhibited developmentally regulated enhancer activity, as measured by massively parallel reporter assays in cardiomyocytes in vivo, and identified key transcription factor-binding motifs. These dynamic enhancers interacted with temporal changes of the 3D genome architecture to specify developmentally regulated cardiomyocyte gene expressions. Our work provides a 3D genome-mediated enhancer activity landscape of murine cardiomyocyte development.

Trajectory reconstruction identifies dysregulation of perinatal maturation programs in pluripotent stem cell-derived cardiomyocytes

A limitation in the application of pluripotent stem cell-derived cardiomyocytes (PSC-CMs) is the failure of these cells to achieve full functional maturity. The mechanisms by which directed differentiation differs from endogenous development, leading to consequent PSC-CM maturation arrest, remain unclear. Here, we generate a single-cell RNA sequencing (scRNA-seq) reference of mouse in vivo CM maturation with extensive sampling of previously difficult-to-isolate perinatal time periods. We subsequently generate isogenic embryonic stem cells to create an in vitro scRNA-seq reference of PSC-CM-directed differentiation. Through trajectory reconstruction, we identify an endogenous perinatal maturation program that is poorly recapitulated in vitro. By comparison with published human datasets, we identify a network of nine transcription factors (TFs) whose targets are consistently dysregulated in PSC-CMs across species. Notably, these TFs are only partially activated in common ex vivo approaches to engineer PSC-CM maturation. Our study can be leveraged toward improving the clinical viability of PSC-CMs.

Impaired Reorganization of Centrosome Structure Underlies Human Infantile Dilated Cardiomyopathy

BACKGROUND

During cardiomyocyte maturation, the centrosome, which functions as a microtubule organizing center in cardiomyocytes, undergoes dramatic structural reorganization where its components reorganize from being localized at the centriole to the nuclear envelope. This developmentally programmed process, referred to as centrosome reduction, has been previously associated with cell cycle exit. However, understanding of how this process influences cardiomyocyte cell biology, and whether its disruption results in human cardiac disease, remains unknown. We studied this phenomenon in an infant with a rare case of infantile dilated cardiomyopathy (iDCM) who presented with left ventricular ejection fraction of 18% and disrupted sarcomere and mitochondria structure.

METHODS

We performed an analysis beginning with an infant who presented with a rare case of iDCM. We derived induced pluripotent stem cells from the patient to model iDCM in vitro. We performed whole exome sequencing on the patient and his parents for causal gene analysis. CRISPR/Cas9-mediated gene knockout and correction in vitro were used to confirm whole exome sequencing results. Zebrafish and Drosophila models were used for in vivo validation of the causal gene. Matrigel mattress technology and single-cell RNA sequencing were used to characterize iDCM cardiomyocytes further.

RESULTS

Whole exome sequencing and CRISPR/Cas9 gene knockout/correction identified RTTN, the gene encoding the centrosomal protein RTTN (rotatin), as the causal gene underlying the patient's condition, representing the first time a centrosome defect has been implicated in a nonsyndromic dilated cardiomyopathy. Genetic knockdowns in zebrafish and Drosophila confirmed an evolutionarily conserved requirement of RTTN for cardiac structure and function. Single-cell RNA sequencing of iDCM cardiomyocytes showed impaired maturation of iDCM cardiomyocytes, which underlie the observed cardiomyocyte structural and functional deficits. We also observed persistent localization of the centrosome at the centriole, contrasting with expected programmed perinuclear reorganization, which led to subsequent global microtubule network defects. In addition, we identified a small molecule that restored centrosome reorganization and improved the structure and contractility of iDCM cardiomyocytes.

CONCLUSIONS

This study is the first to demonstrate a case of human disease caused by a defect in centrosome reduction. We also uncovered a novel role for RTTN in perinatal cardiac development and identified a potential therapeutic strategy for centrosome-related iDCM. Future study aimed at identifying variants in centrosome components may uncover additional contributors to human cardiac disease.

Single-cell RNA sequencing reveals maturation trajectory in human pluripotent stem cell-derived cardiomyocytes in engineered tissues

Cardiac in vitro models have become increasingly obtainable and affordable with the optimization of human pluripotent stem cell-derived cardiomyocyte (hPSC-CM) differentiation. However, these CMs are immature compared to their in vivo counterparts. Here we study the cellular phenotype of hPSC-CMs by comparing their single-cell gene expression and functional profiles in three engineered cardiac tissue configurations: human ventricular (hv) cardiac anisotropic sheet, cardiac tissue strip, and cardiac organoid chamber (hvCOC), with spontaneously aggregated 3D cardiac spheroids (CS) as control. The CM maturity was found to increase with increasing levels of complexity of the engineered tissues from CS to hvCOC. The contractile components are the first function to mature, followed by electrophysiology and oxidative metabolism. Notably, the 2D tissue constructs show a higher cellular organization whereas metabolic maturity preferentially increases in the 3D constructs. We conclude that the tissue engineering models resembling configurations of native tissues may be reliable for drug screening or disease modeling.

Methylation of histone H3 lysine 4 is required for maintenance of beta cell function in adult mice

AIMS/HYPOTHESIS

Beta cells control glucose homeostasis via regulated production and secretion of insulin. This function arises from a highly specialised gene expression programme that is established during development and then sustained, with limited flexibility, in terminally differentiated cells. Dysregulation of this programme is seen in type 2 diabetes but mechanisms that preserve gene expression or underlie its dysregulation in mature cells are not well resolved. This study investigated whether methylation of histone H3 lysine 4 (H3K4), a marker of gene promoters with unresolved functional importance, is necessary for the maintenance of mature beta cell function.

METHODS

Beta cell function, gene expression and chromatin modifications were analysed in conditional Dpy30 knockout mice, in which H3K4 methyltransferase activity is impaired, and in a mouse model of diabetes.

RESULTS

H3K4 methylation maintains expression of genes that are important for insulin biosynthesis and glucose responsiveness. Deficient methylation of H3K4 leads to a less active and more repressed epigenome profile that locally correlates with gene expression deficits but does not globally reduce gene expression. Instead, developmentally regulated genes and genes in weakly active or suppressed states particularly rely on H3K4 methylation. We further show that H3K4 trimethylation (H3K4me3) is reorganised in islets from the Lepr^db/db mouse model of diabetes in favour of weakly active and disallowed genes at the expense of terminal beta cell markers with broad H3K4me3 peaks.

CONCLUSIONS/INTERPRETATION

Sustained methylation of H3K4 is critical for the maintenance of beta cell function. Redistribution of H3K4me3 is linked to gene expression changes that are implicated in diabetes pathology.

In vivo-like Culture of Monophagous Animal Organ using Dietary Components

Animals depend on other species to live, with monophagy being an extreme mode. Monophagous animals depend on their diet not only for nutritients but also for developmental and reproductive controls. Thus, dietary components may be useful in culturing tissues from monophagous animals. We hypothesized that a dedifferentiated tissue from the monophagous silkworm, Bombyx mori, would re-differentiate when cultured in a medium containing an extract of mulberry (Morus alba) leaves, the only food of B. mori. Over 40 fat-body transcriptomes were sequenced, and we concluded that it is possible to establish in vivo-like silkworm tissue cultures using their diet.

Transcriptome Dedifferentiation Observed in Animal Primary Cultures is Essential to Plant Reprogramming

Tissue culture environment liberate cells from ordinary laws of multi-cellular organisms. This liberation enables cells several behaviors, such as growth, dedifferentiation, acquisition of pluripotency, immortalization and reprogramming. Each phenomenon is relating to each other and hardly to determine. Recently, dedifferentiation of animal cell was quantified as increasing liberality which is information entropy of transcriptome. The increasing liberality induced by tissue culture may reappear in plant cells too. Here we corroborated it. Measuring liberality during reprogramming of plant cells suggested that reprogramming is a combined phenomenon of dedifferentiation and re-differentiation.

Cross-Organ Transcriptomic Comparison Reveals Universal Factors During Maturation

Various cell types can be derived from stem cells. However, these cells are immature and do not match their adult counterparts in functional capabilities, limiting their use in disease modeling and cell therapies. Thus, it is crucial to understand the mechanisms of maturation in vivo. However, it is unknown if there are genes and pathways conserved across organs during maturation. To address this, we performed a time-series analysis of the transcriptome of the mouse heart, brain, liver, and kidney and analyzed their trajectories over time. In addition, gene regulatory networks were reconstructed to determine overlapping expression patterns. Based on these, we identified commonly upregulated and downregulated pathways across all four organs. Key upstream regulators were also predicted based on the temporal expression of downstream genes. These findings suggest the presence of universal regulators during organ maturation, which may help us develop a general strategy to mature stem cell-derived cells in vitro.

Using human induced pluripotent stem cell-derived cardiomyocytes to understand the mechanisms driving cardiomyocyte maturation

Cardiovascular diseases are the leading cause of mortality and reduced quality of life globally. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) provide a personalized platform to study inherited heart diseases, drug-induced cardiac toxicity, and cardiac regenerative therapy. However, the immaturity of CMs obtained by current strategies is a major hurdle in utilizing hiPSC-CMs at their fullest potential. Here, the major findings and limitations of current maturation methodologies to enhance the utility of hiPSC-CMs in the battle against a major source of morbidity and mortality are reviewed. The most recent knowledge of the potential signaling pathways involved in the transition of fetal to adult CMs are assimilated. In particular, we take a deeper look on role of nutrient sensing signaling pathways and the potential role of cap-independent translation mediated by the modulation of mTOR pathway in the regulation of cardiac gap junctions and other yet to be identified aspects of CM maturation. Moreover, a relatively unexplored perspective on how our knowledge on the effects of preterm birth on cardiovascular development can be actually utilized to enhance the current understanding of CM maturation is examined. Furthermore, the interaction between the evolving neonatal human heart and brown adipose tissue as the major source of neonatal thermogenesis and its endocrine function on CM development is another discussed topic which is worthy of future investigation. Finally, the current knowledge regarding transcriptional mediators of CM maturation is still limited. The recent studies have produced the groundwork to better understand CM maturation in terms of providing some of the key factors involved in maturation and development of metrics for assessment of maturation which proves essential for future studies on in vitro PSC-CMs maturation.

FitDevo: accurate inference of single-cell developmental potential using sample-specific gene weight

The quantification of developmental potential is critical for determining developmental stages and identifying essential molecular signatures in single-cell studies. Here, we present FitDevo, a novel method for inferring developmental potential using scRNA-seq data. The main idea of FitDevo is first to generate sample-specific gene weight (SSGW) and then infer developmental potential by calculating the correlation between SSGW and gene expression. SSGW is generated using a generalized linear model that combines sample-specific information and gene weight learned from a training dataset covering scRNA-seq data of 17 previously published datasets. We have rigorously validated FitDevo’s effectiveness using a testing dataset with scRNA-seq data from 28 existing datasets and have also demonstrated its superiority over current methods. Furthermore, FitDevo’s broad application scope has been illustrated using three practical scenarios: deconvolution analysis of epidermis, spatial transcriptomic data analysis of hearts and intestines, and developmental potential analysis of breast cancer. The source code and related data are available at https://github.com/jumphone/fitdevo.

Functional genomics and the future of iPSCs in disease modeling

Induced pluripotent stem cells (iPSCs) are valuable in disease modeling because of their potential to expand and differentiate into virtually any cell type and recapitulate key aspects of human biology. Functional genomics are genome-wide studies that aim to discover genotype-phenotype relationships, thereby revealing the impact of human genetic diversity on normal and pathophysiology. In this review, we make the case that human iPSCs (hiPSCs) are a powerful tool for functional genomics, since they provide an in vitro platform for the study of population genetics. We describe cutting-edge tools and strategies now available to researchers, including multi-omics technologies, advances in hiPSC culture techniques, and innovations in drug development. Functional genomics approaches based on hiPSCs hold great promise for advancing drug discovery, disease etiology, and the impact of genetic variation on human biology.

Deciphering Cardiac Biology and Disease by Single-Cell Transcriptomic Profiling

By detecting minute molecular changes in hundreds to millions of single cells, single-cell RNA sequencing allows for the comprehensive characterization of the diversity and dynamics of cells in the heart. Our understanding of the heart has been transformed through the recognition of cellular heterogeneity, the construction of regulatory networks, the building of lineage trajectories, and the mapping of intercellular crosstalk. In this review, we introduce cardiac progenitors and their transcriptional regulation during embryonic development, highlight cellular heterogeneity and cell subtype functions in cardiac health and disease, and discuss insights gained from the study of pluripotent stem-cell-derived cardiomyocytes.