Association of Human iPSC Gene Signatures and X Chromosome Dosage with Two Distinct Cardiac Differentiation Trajectories

Sex in cardiovascular disease: Why this biological variable should be considered in in vitro models

Cardiovascular disease (CVD), the world's leading cause of death, exhibits notable epidemiological, clinical, and pathophysiological differences between sexes. Many such differences can be linked back to cardiovascular sexual dimorphism, yet sex-specific in vitro models are still not the norm. A lack of sex reporting and apparent male bias raises the question of whether in vitro CVD models faithfully recapitulate the biology of intended treatment recipients. To ensure equitable treatment for the overlooked female patient population, sex as a biological variable (SABV) inclusion must become commonplace in CVD preclinical research. Here, we discuss the role of sex in CVD and underlying cardiovascular (patho)physiology. We review shortcomings in current SABV practices, describe the relevance of sex, and highlight emerging strategies for SABV inclusion in three major in vitro model types: primary cell, stem cell, and three-dimensional models. Last, we identify key barriers to inclusive design and suggest techniques for overcoming them.

IFN-γ activates an immune-like regulatory network in the cardiac vascular endothelium

The regulatory mechanisms underlying the response to pro-inflammatory cytokines during myocarditis are poorly understood. Here, we use iPSC-derived cardiovascular progenitor cells (CVPCs) to model the response to interferon gamma (IFN-γ) during myocarditis. We generate RNA-seq and ATAC-seq for four CVPCs that were treated with IFN-γ and compare them with paired untreated controls. Transcriptional differences after treatment show that IFN-γ initiates an innate immune cell-like response in the vascular cardiac endothelium. IFN-γ treatment also shifts the CVPC transcriptome towards the adult coronary artery and aorta profiles and expands the relative endothelial cell population in all four CVPC lines. Analysis of the accessible chromatin shows that IFN-γ is a potent chromatin remodeler and establishes an IRF-STAT immune-cell like regulatory network. Our findings reveal insights into the endothelial-specific protective mechanisms during myocarditis.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

Noninvasive analysis of contractility during identical maturations revealed two phenotypes in ventricular but not in atrial iPSC-CM

Patient-derived induced pluripotent stem cells (iPSCs) can be differentiated into atrial and ventricular cardiomyocytes to allow for personalized drug screening. A hallmark of differentiation is the manifestation of spontaneous beating in a two-dimensional (2-D) cell culture. However, an outstanding observation is the high variability in this maturation process. We valued that contractile parameters change during differentiation serving as an indicator of maturation. Consequently, we recorded noninvasively spontaneous motion activity during the differentiation of male iPSC toward iPSC cardiomyocytes (iPSC-CMs) to further analyze similar maturated iPSC-CMs. Surprisingly, our results show that identical differentiations into ventricular iPSC-CMs are variable with respect to contractile parameters resulting in two distinct subpopulations of ventricular-like cells. In contrast, differentiation into atrial iPSC-CMs resulted in only one phenotype. We propose that the noninvasive and cost-effective recording of contractile activity during maturation using a smartphone device may help to reduce the variability in results frequently reported in studies on ventricular iPSC-CMs.NEW & NOTEWORTHY Differentiation of induced pluripotent stem cells (iPSCs) into iPSC-derived cardiomyocytes (iPSC-CMs) exhibits a high variability in mature parameters. Here, we monitored noninvasively contractile parameters of iPSC-CM during full-time differentiation using a smartphone device. Our results show that parallel maturations of iPSCs into ventricular iPSC-CMs, but not into atrial iPSC-CMs, resulted in two distinct subpopulations of iPSC-CMs. These findings suggest that our cost-effective method may help to compare iPSC-CMs at the same maturation level.

Complex regulatory networks influence pluripotent cell state transitions in human iPSCs

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

eQTL mapping in fetal-like pancreatic progenitor cells reveals early developmental insights into diabetes risk

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

APLNR marks a cardiac progenitor derived with human induced pluripotent stem cells

Cardiomyocytes can be readily derived from human induced pluripotent stem cell (hiPSC) lines, yet its efficacy varies across different batches of the same and different hiPSC lines. To unravel the inconsistencies of in vitro cardiac differentiation, we utilized single cell transcriptomics on hiPSCs undergoing cardiac differentiation and identified cardiac and extra-cardiac lineages throughout differentiation. We further identified APLNR as a surface marker for in vitro cardiac progenitors and immunomagnetically isolated them. Differentiation of isolated in vitro APLNR+ cardiac progenitors derived from multiple hiPSC lines resulted in predominantly cardiomyocytes accompanied with cardiac mesenchyme. Transcriptomic analysis of differentiating in vitro APLNR+ cardiac progenitors revealed transient expression of cardiac progenitor markers before further commitment into cardiomyocyte and cardiac mesenchyme. Analysis of in vivo human and mouse embryo single cell transcriptomic datasets have identified APLNR expression in early cardiac progenitors of multiple lineages. This platform enables generation of in vitro cardiac progenitors from multiple hiPSC lines without genetic manipulation, which has potential applications in studying cardiac development, disease modelling and cardiac regeneration.

Sex as a Biological Variable in Tissue Engineering and Regenerative Medicine

Although sex differences have been noted in cellular function and behavior, therapy efficacy, and disease incidence and outcomes, the adoption of sex as a biological variable in tissue engineering and regenerative medicine remains limited. Furthering the development of personalized, precision medicine requires considering biological sex at the bench and in the clinic. This review provides the basis for considering biological sex when designing tissue-engineered constructs and regenerative therapies by contextualizing sex as a biological variable within the tissue engineering triad of cells, matrices, and signals. To achieve equity in biological sex within medicine requires a cultural shift in science and engineering research, with active engagement by researchers, clinicians, companies, policymakers, and funding agencies.

The Foundational Data Initiative for Parkinson Disease: Enabling efficient translation from genetic maps to mechanism

The Foundational Data Initiative for Parkinson Disease (FOUNDIN-PD) is an international collaboration producing fundamental resources for Parkinson disease (PD). FOUNDIN-PD generated a multi-layered molecular dataset in a cohort of induced pluripotent stem cell (iPSC) lines differentiated to dopaminergic (DA) neurons, a major affected cell type in PD. The lines were derived from the Parkinson's Progression Markers Initiative study, which included participants with PD carrying monogenic PD variants, variants with intermediate effects, and variants identified by genome-wide association studies and unaffected individuals. We generated genetic, epigenetic, regulatory, transcriptomic, and longitudinal cellular imaging data from iPSC-derived DA neurons to understand molecular relationships between disease-associated genetic variation and proximate molecular events. These data reveal that iPSC-derived DA neurons provide a valuable cellular context and foundational atlas for modeling PD genetic risk. We have integrated these data into a FOUNDIN-PD data browser as a resource for understanding the molecular pathogenesis of PD.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Action potential variability in human pluripotent stem cell-derived cardiomyocytes obtained from healthy donors

Human pluripotent stem cells (PSC) have been used for disease modelling, after differentiation into the desired cell type. Electrophysiologic properties of cardiomyocytes derived from pluripotent stem cells are extensively used to model cardiac arrhythmias, in cardiomyopathies and channelopathies. This requires strict control of the multiple variables that can influence the electrical properties of these cells. In this article, we report the action potential variability of 780 cardiomyocytes derived from pluripotent stem cells obtained from six healthy donors. We analyze the overall distribution of action potential (AP) data, the distribution of action potential data per cell line, per differentiation protocol and batch. This analysis indicates that even using the same cell line and differentiation protocol, the differentiation batch still affects the results. This variability has important implications in modeling arrhythmias and imputing pathogenicity to variants encountered in patients with arrhythmic diseases. We conclude that even when using isogenic cell lines to ascertain pathogenicity to variants associated to arrythmias one should use cardiomyocytes derived from pluripotent stem cells using the same differentiation protocol and batch and pace the cells or use only cells that have very similar spontaneous beat rates. Otherwise, one may find phenotypic variability that is not attributable to pathogenic variants.

Cell type determination for cardiac differentiation occurs soon after seeding of human-induced pluripotent stem cells

BACKGROUND

Cardiac differentiation of human-induced pluripotent stem (hiPS) cells consistently produces a mixed population of cardiomyocytes and non-cardiac cell types, even when using well-characterized protocols. We sought to determine whether different cell types might result from intrinsic differences in hiPS cells prior to the onset of differentiation.

RESULTS

By associating individual differentiated cells that share a common hiPS cell precursor, we tested whether expression variability is predetermined from the hiPS cell state. In a single experiment, cells that shared a progenitor were more transcriptionally similar to each other than to other cells in the differentiated population. However, when the same hiPS cells were differentiated in parallel, we did not observe high transcriptional similarity across differentiations. Additionally, we found that substantial cell death occurs during differentiation in a manner that suggested all cells were equally likely to survive or die, suggesting that there is no intrinsic selection bias for cells descended from particular hiPS cell progenitors. We thus wondered how cells grow spatially during differentiation, so we labeled cells by expression of marker genes and found that cells expressing the same marker tended to occur in patches. Our results suggest that cell type determination across multiple cell types, once initiated, is maintained in a cell-autonomous manner for multiple divisions.

CONCLUSIONS

Altogether, our results show that while substantial heterogeneity exists in the initial hiPS cell population, it is not responsible for the variability observed in differentiated outcomes; instead, factors specifying the various cell types likely act during a window that begins shortly after the seeding of hiPS cells for differentiation.

In heart failure reactivation of RNA-binding proteins is associated with the expression of 1,523 fetal-specific isoforms

Reactivation of fetal-specific genes and isoforms occurs during heart failure. However, the underlying molecular mechanisms and the extent to which the fetal program switch occurs remains unclear. Limitations hindering transcriptome-wide analyses of alternative splicing differences (i.e. isoform switching) in cardiovascular system (CVS) tissues between fetal, healthy adult and heart failure have included both cellular heterogeneity across bulk RNA-seq samples and limited availability of fetal tissue for research. To overcome these limitations, we have deconvoluted the cellular compositions of 996 RNA-seq samples representing heart failure, healthy adult (heart and arteria), and fetal-like (iPSC-derived cardiovascular progenitor cells) CVS tissues. Comparison of the expression profiles revealed that reactivation of fetal-specific RNA-binding proteins (RBPs), and the accompanied re-expression of 1,523 fetal-specific isoforms, contribute to the transcriptome differences between heart failure and healthy adult heart. Of note, isoforms for 20 different RBPs were among those that reverted in heart failure to the fetal-like expression pattern. We determined that, compared with adult-specific isoforms, fetal-specific isoforms encode proteins that tend to have more functions, are more likely to harbor RBP binding sites, have canonical sequences at their splice sites, and contain typical upstream polypyrimidine tracts. Our study suggests that compared with healthy adult, fetal cardiac tissue requires stricter transcriptional regulation, and that during heart failure reversion to this stricter transcriptional regulation occurs. Furthermore, we provide a resource of cardiac developmental stage-specific and heart failure-associated genes and isoforms, which are largely unexplored and can be exploited to investigate novel therapeutics for heart failure.

Heart failure is a chronic condition in which the heart does not pump enough blood. It has been shown that in heart failure, the adult heart reverts to a fetal-like metabolic state and oxygen consumption. Additionally, genes and isoforms that are expressed in the heart only during fetal development (i.e. not in the healthy adult heart) are turned on in heart failure. However, the underlying molecular mechanisms and the extent to which the switch to a fetal gene program occurs remains unclear. In this study, we initially characterized the differences between the fetal and adult heart transcriptomes (entire set of expressed genes and isoforms). We found that RNA binding proteins (RBPs), a family of genes that regulate multiple aspects of a transcript’s maturation, including transcription, splicing and post-transcriptional modifications, play a central role in the differences between fetal and adult heart tissues. We observed that many RBPs that are only expressed in the heart during fetal development become reactivated in heart failure, resulting in the expression of 1,523 fetal-specific isoforms. These findings suggest that reactivation of fetal-specific RBPs in heart failure drives a transcriptome-wide switch to expression of fetal-specific isoforms; and hence that RBPs could potentially serve as novel therapeutic targets.

Ultra-Sharp Nanowire Arrays Natively Permeate, Record, and Stimulate Intracellular Activity in Neuronal and Cardiac Networks

We report innovative scalable, vertical, ultra-sharp nanowire arrays that are individually addressable to enable long-term, native recordings of intracellular potentials. Stable amplitudes of intracellular potentials from 3D tissue-like networks of neurons and cardiomyocytes are obtained. Individual electrical addressability is necessary for high-fidelity intracellular electrophysiological recordings. This study paves the way toward predictive, high-throughput, and low-cost electrophysiological drug screening platforms.

Single-cell sequencing reveals lineage-specific dynamic genetic regulation of gene expression during human cardiomyocyte differentiation

Dynamic and temporally specific gene regulatory changes may underlie unexplained genetic associations with complex disease. During a dynamic process such as cellular differentiation, the overall cell type composition of a tissue (or an in vitro culture) and the gene regulatory profile of each cell can both experience significant changes over time. To identify these dynamic effects in high resolution, we collected single-cell RNA-sequencing data over a differentiation time course from induced pluripotent stem cells to cardiomyocytes, sampled at 7 unique time points in 19 human cell lines. We employed a flexible approach to map dynamic eQTLs whose effects vary significantly over the course of bifurcating differentiation trajectories, including many whose effects are specific to one of these two lineages. Our study design allowed us to distinguish true dynamic eQTLs affecting a specific cell lineage from expression changes driven by potentially non-genetic differences between cell lines such as cell composition. Additionally, we used the cell type profiles learned from single-cell data to deconvolve and re-analyze data from matched bulk RNA-seq samples. Using this approach, we were able to identify a large number of novel dynamic eQTLs in single cell data while also attributing dynamic effects in bulk to a particular lineage. Overall, we found that using single cell data to uncover dynamic eQTLs can provide new insight into the gene regulatory changes that occur among heterogeneous cell types during cardiomyocyte differentiation.

Capturing Pluripotency and Beyond

During the development of a multicellular organism, the specification of different cell lineages originates in a small group of pluripotent cells, the epiblasts, formed in the preimplantation embryo. The pluripotent epiblast is protected from premature differentiation until exposure to inductive cues in strictly controlled spatially and temporally organized patterns guiding fetus formation. Epiblasts cultured in vitro are embryonic stem cells (ESCs), which recapitulate the self-renewal and lineage specification properties of their endogenous counterparts. The characteristics of totipotency, although less understood than pluripotency, are becoming clearer. Recent studies have shown that a minor ESC subpopulation exhibits expanded developmental potential beyond pluripotency, displaying a characteristic reminiscent of two-cell embryo blastomeres (2CLCs). In addition, reprogramming both mouse and human ESCs in defined media can produce expanded/extended pluripotent stem cells (EPSCs) similar to but different from 2CLCs. Further, the molecular roadmaps driving the transition of various potency states have been clarified. These recent key findings will allow us to understand eutherian mammalian development by comparing the underlying differences between potency network components during development. Using the mouse as a paradigm and recent progress in human PSCs, we review the epiblast's identity acquisition during embryogenesis and their ESC counterparts regarding their pluripotent fates and beyond.

A library of induced pluripotent stem cells from clinically well-characterized, diverse healthy human individuals

A library of well-characterized human induced pluripotent stem cell (hiPSC) lines from clinically healthy human subjects could serve as a useful resource of normal controls for in vitro human development, disease modeling, genotype-phenotype association studies, and drug response evaluation. We report generation and extensive characterization of a gender-balanced, racially/ethnically diverse library of hiPSC lines from 40 clinically healthy human individuals who range in age from 22 to 61 years. The hiPSCs match the karyotype and short tandem repeat identities of their parental fibroblasts, and have a transcription profile characteristic of pluripotent stem cells. We provide whole-genome sequencing data for one hiPSC clone from each individual, genomic ancestry determination, and analysis of mendelian disease genes and risks. We document similar transcriptomic profiles, single-cell RNA-sequencing-derived cell clusters, and physiology of cardiomyocytes differentiated from multiple independent hiPSC lines. This extensive characterization makes this hiPSC library a valuable resource for many studies on human biology.

•

A library of induced pluripotent stem cells from 40 healthy human subjects

•

Racially/ethnically diverse subjects of clinically well-characterized health

•

Whole-genome sequencing identifies variants of mild common phenotypes or incomplete penetrance

•

Similar physiology of cardiomyocytes from independent hiPSC clones and individuals

A framework for developing sex-specific engineered heart models

The convergence of tissue engineering and patient-specific stem cell biology has enabled the engineering of in vitro tissue models that allow the study of patient-tailored treatment modalities. However, sex-related disparities in health and disease, from systemic hormonal influences to cellular-level differences, are often overlooked in stem cell biology, tissue engineering and preclinical screening. The cardiovascular system, in particular, shows considerable sex-related differences, which need to be considered in cardiac tissue engineering. In this Review, we analyse sex-related properties of the heart muscle in the context of health and disease, and discuss a framework for including sex-based differences in human cardiac tissue engineering. We highlight how sex-based features can be implemented at the cellular and tissue levels, and how sex-specific cardiac models could advance the study of cardiovascular diseases. Finally, we define design criteria for sex-specific cardiac tissue engineering and provide an outlook to future research possibilities beyond the cardiovascular system.

SCReadCounts: estimation of cell-level SNVs expression from scRNA-seq data

BACKGROUND

Recent studies have demonstrated the utility of scRNA-seq SNVs to distinguish tumor from normal cells, characterize intra-tumoral heterogeneity, and define mutation-associated expression signatures. In addition to cancer studies, SNVs from single cells have been useful in studies of transcriptional burst kinetics, allelic expression, chromosome X inactivation, ploidy estimations, and haplotype inference.

RESULTS

To aid these types of studies, we have developed a tool, SCReadCounts, for cell-level tabulation of the sequencing read counts bearing SNV reference and variant alleles from barcoded scRNA-seq alignments. Provided genomic loci and expected alleles, SCReadCounts generates cell-SNV matrices with the absolute variant- and reference-harboring read counts, as well as cell-SNV matrices of expressed Variant Allele Fraction (VAFRNA) suitable for a variety of downstream applications. We demonstrate three different SCReadCounts applications on 59,884 cells from seven neuroblastoma samples: (1) estimation of cell-level expression of known somatic mutations and RNA-editing sites, (2) estimation of cell- level allele expression of biallelic SNVs, and (3) a discovery mode assessment of the reference and each of the three alternative nucleotides at genomic positions of interest that does not require prior SNV information. For the later, we applied SCReadCounts on the coding regions of KRAS, where it identified known and novel somatic mutations in a low-to-moderate proportion of cells. The SCReadCounts read counts module is benchmarked against the analogous modules of GATK and Samtools. SCReadCounts is freely available ( https://github.com/HorvathLab/NGS ) as 64-bit self-contained binary distributions for Linux and MacOS, in addition to Python source.

CONCLUSIONS

SCReadCounts supplies a fast and efficient solution for estimation of cell-level SNV expression from scRNA-seq data. SCReadCounts enables distinguishing cells with monoallelic reference expression from those with no gene expression and is applicable to assess SNVs present in only a small proportion of the cells, such as somatic mutations in cancer.

Sexual Dimorphism of the Heart: Genetics, Epigenetics, and Development

The democratization of genomic technologies has revealed profound sex biases in expression patterns in every adult tissue, even in organs with no conspicuous differences, such as the heart. With the increasing awareness of the disparities in cardiac pathophysiology between males and females, there are exciting opportunities to explore how sex differences in the heart are established developmentally. Although sexual dimorphism is traditionally attributed to hormonal influence, expression and epigenetic sex biases observed in early cardiac development can only be accounted for by the difference in sex chromosome composition, i.e., XX in females and XY in males. In fact, genes linked to the X and Y chromosomes, many of which encode regulatory factors, are expressed in cardiac progenitor cells and at every subsequent developmental stage. The effect of the sex chromosome composition may explain why many congenital heart defects originating before gonad formation exhibit sex biases in presentation, mortality, and morbidity. Some transcriptional and epigenetic sex biases established soon after fertilization persist in cardiac lineages, suggesting that early epigenetic events are perpetuated beyond early embryogenesis. Importantly, when sex hormones begin to circulate, they encounter a cardiac genome that is already functionally distinct between the sexes. Although there is a wealth of knowledge on the effects of sex hormones on cardiac function, we propose that sex chromosome-linked genes and their downstream targets also contribute to the differences between male and female hearts. Moreover, identifying how hormones influence sex chromosome effects, whether antagonistically or synergistically, will enhance our understanding of how sex disparities are established. We also explore the possibility that sexual dimorphism of the developing heart predicts sex-specific responses to environmental signals and foreshadows sex-biased health-related outcomes after birth.

Identification of X-chromosomal genes that drive sex differences in embryonic stem cells through a hierarchical CRISPR screening approach

BACKGROUND

X-chromosomal genes contribute to sex differences, in particular during early development, when both X chromosomes are active in females. Double X-dosage shifts female pluripotent cells towards the naive stem cell state by increasing pluripotency factor expression, inhibiting the differentiation-promoting MAP kinase (MAPK) signaling pathway, and delaying differentiation.

RESULTS

To identify the genetic basis of these sex differences, we use a two-step CRISPR screening approach to comprehensively identify X-linked genes that cause the female pluripotency phenotype in murine embryonic stem cells. A primary chromosome-wide CRISPR knockout screen and three secondary screens assaying for different aspects of the female pluripotency phenotype allow us to uncover multiple genes that act in concert and to disentangle their relative roles. Among them, we identify Dusp9 and Klhl13 as two central players. While Dusp9 mainly affects MAPK pathway intermediates, Klhl13 promotes pluripotency factor expression and delays differentiation, with both factors jointly repressing MAPK target gene expression.

CONCLUSIONS

Here, we elucidate the mechanisms that drive sex-induced differences in pluripotent cells and our approach serves as a blueprint to discover the genetic basis of the phenotypic consequences of other chromosomal effects.

Population-scale single-cell RNA-seq profiling across dopaminergic neuron differentiation

Studying the function of common genetic variants in primary human tissues and during development is challenging. To address this, we use an efficient multiplexing strategy to differentiate 215 human induced pluripotent stem cell (iPSC) lines toward a midbrain neural fate, including dopaminergic neurons, and use single-cell RNA sequencing (scRNA-seq) to profile over 1 million cells across three differentiation time points. The proportion of neurons produced by each cell line is highly reproducible and is predictable by robust molecular markers expressed in pluripotent cells. Expression quantitative trait loci (eQTL) were characterized at different stages of neuronal development and in response to rotenone-induced oxidative stress. Of these, 1,284 eQTL colocalize with known neurological trait risk loci, and 46% are not found in the Genotype-Tissue Expression (GTEx) catalog. Our study illustrates how coupling scRNA-seq with long-term iPSC differentiation enables mechanistic studies of human trait-associated genetic variants in otherwise inaccessible cell states.

Revealing Instability: Genetic Variation Underlies Variability in mESC Pluripotency

mESCs can self-renew indefinitely in vitro. However, depending on culture conditions some strains are more unstable than others. In this issue of Cell Stem Cell, Skelly et al. (2020) and Ortmann et al. (2020) shed light into the role genetic variation plays in control of ground state pluripotency.

In vitro Differentiation of Human iPSC-derived Cardiovascular Progenitor Cells (iPSC-CVPCs)

Induced pluripotent stem cell derived cardiovascular progenitor cells (iPSC-CVPCs) provide an unprecedented platform for examining the molecular underpinnings of cardiac development and disease etiology, but also have great potential to play pivotal roles in the future of regenerative medicine and pharmacogenomic studies. Biobanks like iPSCORE ( Stacey et al., 2013 ; Panopoulos et al., 2017 ), which contain iPSCs generated from hundreds of genetically and ethnically diverse individuals, are an invaluable resource for conducting these studies. Here, we present an optimized, cost-effective and highly standardized protocol for large-scale derivation of human iPSC-CVPCs using small molecules and purification using metabolic selection. We have successfully applied this protocol to derive iPSC-CVPCs from 154 different iPSCORE iPSC lines obtaining large quantities of highly pure cardiac cells. An important component of our protocol is Cell confluency estimates (ccEstimate), an automated methodology for estimating the time when an iPSC monolayer will reach 80% confluency, which is optimal for initiating iPSC-CVPC derivation, and enables the protocol to be readily used across iPSC lines with different growth rates. Moreover, we showed that cellular heterogeneity across iPSC-CVPCs is due to varying proportions of two distinct cardiac cell types: cardiomyocytes (CMs) and epicardium-derived cells (EPDCs), both of which have been shown to have a critical function in heart regeneration. This protocol eliminates the need of iPSC line-to-line optimization and can be easily adapted and scaled to high-throughput studies or to generate large quantities of cells suitable for regenerative medicine applications.

Variation in RARG increases susceptibility to doxorubicin-induced cardiotoxicity in patient specific induced pluripotent stem cell-derived cardiomyocytes

Doxorubicin is a potent anticancer drug used to treat a variety of cancer types. However, its use is limited by doxorubicin-induced cardiotoxicity (DIC). A missense variant in the RARG gene (S427L; rs2229774) has been implicated in susceptibility to DIC in a genome wide association study. The goal of this study was to investigate the functional role of this RARG variant in DIC. We used induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) from patients treated with doxorubicin. iPSC-CMs from individuals who experienced DIC (cases) showed significantly greater sensitivity to doxorubicin compared to iPSC-CMs from doxorubicin-treated individuals who did not develop DIC (controls) in cell viability and optical mapping experiments. Using CRISPR/Cas9, we generated isogenic cell lines that differed only at the RARG locus. Genetic correction of RARG-S427L to wild type resulted in reduced doxorubicin-induced double stranded DNA breaks, reactive oxygen species production, and cell death. Conversely, introduction of RARG-S427L increased susceptibility to doxorubicin. Finally, genetic disruption of the RARG gene resulted in protection from cell death due to doxorubicin treatment. Our findings suggest that the presence of RARG-S427L increases sensitivity to DIC, establishing a direct, causal role for this variant in DIC.

Practical adoption of state-of-the-art hiPSC-cardiomyocyte differentiation techniques

Human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes are a valuable resource for cardiac therapeutic development; however, generation of these cells in large numbers and high purity is a limitation in widespread adoption. Here, design of experiments (DOE) is used to investigate the cardiac differentiation space of three hiPSC lines when varying CHIR99027 concentration and cell seeding density, and a novel image analysis is developed to evaluate plate coverage when initiating differentiation. Metabolic selection via lactate purifies hiPSC-cardiomyocyte populations, and the bioenergetic phenotype and engineered tissue mechanics of purified and unpurified hiPSC-cardiomyocytes are compared. Findings demonstrate that when initiating differentiation one day after hiPSC plating, low (3 μM) Chiron and 72 x 10³ cells/cm² seeding density result in peak cardiac purity (50–90%) for all three hiPSC lines. Our results confirm that metabolic selection with lactate shifts hiPSC-cardiomyocyte metabolism towards oxidative phosphorylation, but this more “mature” metabolic phenotype does not by itself result in a more mature contractile phenotype in engineered cardiac tissues at one week of culture in 3D tissues. This study provides widely adaptable methods including novel image analysis code and parameters for refining hiPSC-cardiomyocyte differentiation and describes the practical implications of metabolic selection of cardiomyocytes for downstream tissue engineering applications.

Cyclooxygenases Inhibitors Efficiently Induce Cardiomyogenesis in Human Pluripotent Stem Cells

Application of human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) is limited by the challenges in their efficient differentiation. Recently, the Wingless (Wnt) signaling pathway has emerged as the key regulator of cardiomyogenesis. In this study, we evaluated the effects of cyclooxygenase inhibitors on cardiac differentiation of hPSCs. Cardiac differentiation was performed by adherent monolayer based method using 4 hPSC lines (HES3, H9, IMR90, and ES4SKIN). The efficiency of cardiac differentiation was evaluated by flow cytometry and RT-qPCR. Generated hPSC-CMs were characterised using immunocytochemistry, electrophysiology, electron microscopy, and calcium transient measurements. Our data show that the COX inhibitors Sulindac and Diclofenac in combination with CHIR99021 (GSK-3 inhibitor) efficiently induce cardiac differentiation of hPSCs. In addition, inhibition of COX using siRNAs targeted towards COX-1 and/or COX-2 showed that inhibition of COX-2 alone or COX-1 and COX-2 in combination induce cardiomyogenesis in hPSCs within 12 days. Using IMR90-Wnt reporter line, we showed that inhibition of COX-2 led to downregulation of Wnt signalling activity in hPSCs. In conclusion, this study demonstrates that COX inhibition efficiently induced cardiogenesis via modulation of COX and Wnt pathway and the generated cardiomyocytes express cardiac-specific structural markers as well as exhibit typical calcium transients and action potentials. These cardiomyocytes also responded to cardiotoxicants and can be relevant as an in vitro cardiotoxicity screening model.

Estimating the Allele-Specific Expression of SNVs From 10× Genomics Single-Cell RNA-Sequencing Data

With the recent advances in single-cell RNA-sequencing (scRNA-seq) technologies, the estimation of allele expression from single cells is becoming increasingly reliable. Allele expression is both quantitative and dynamic and is an essential component of the genomic interactome. Here, we systematically estimate the allele expression from heterozygous single nucleotide variant (SNV) loci using scRNA-seq data generated on the 10×Genomics Chromium platform. We analyzed 26,640 human adipose-derived mesenchymal stem cells (from three healthy donors), sequenced to an average of 150K sequencing reads per cell (more than 4 billion scRNA-seq reads in total). High-quality SNV calls assessed in our study contained approximately 15% exonic and >50% intronic loci. To analyze the allele expression, we estimated the expressed variant allele fraction (VAF_RNA) from SNV-aware alignments and analyzed its variance and distribution (mono- and bi-allelic) at different minimum sequencing read thresholds. Our analysis shows that when assessing positions covered by a minimum of three unique sequencing reads, over 50% of the heterozygous SNVs show bi-allelic expression, while at a threshold of 10 reads, nearly 90% of the SNVs are bi-allelic. In addition, our analysis demonstrates the feasibility of scVAF_RNA estimation from current scRNA-seq datasets and shows that the 3′-based library generation protocol of 10×Genomics scRNA-seq data can be informative in SNV-based studies, including analyses of transcriptional kinetics.

RNA-Based Strategies for Cardiac Reprogramming of Human Mesenchymal Stromal Cells

Multipotent adult mesenchymal stromal cells (MSCs) could represent an elegant source for the generation of patient-specific cardiomyocytes needed for regenerative medicine, cardiovascular research, and pharmacological studies. However, the differentiation of adult MSC into a cardiac lineage is challenging compared to embryonic stem cells or induced pluripotent stem cells. Here we used non-integrative methods, including microRNA and mRNA, for cardiac reprogramming of adult MSC derived from bone marrow, dental follicle, and adipose tissue. We found that MSC derived from adipose tissue can partly be reprogrammed into the cardiac lineage by transient overexpression of GATA4, TBX5, MEF2C, and MESP1, while cells isolated from bone marrow, and dental follicle exhibit only weak reprogramming efficiency. qRT-PCR and transcriptomic analysis revealed activation of a cardiac-specific gene program and up-regulation of genes known to promote cardiac development. Although we did not observe the formation of fully mature cardiomyocytes, our data suggests that adult MSC have the capability to acquire a cardiac-like phenotype when treated with mRNA coding for transcription factors that regulate heart development. Yet, further optimization of the reprogramming process is mandatory to increase the reprogramming efficiency.

Systematic genetic analysis of the MHC region reveals mechanistic underpinnings of HLA type associations with disease

The MHC region is highly associated with autoimmune and infectious diseases. Here we conduct an in-depth interrogation of associations between genetic variation, gene expression and disease. We create a comprehensive map of regulatory variation in the MHC region using WGS from 419 individuals to call eight-digit HLA types and RNA-seq data from matched iPSCs. Building on this regulatory map, we explored GWAS signals for 4083 traits, detecting colocalization for 180 disease loci with eQTLs. We show that eQTL analyses taking HLA type haplotypes into account have substantially greater power compared with only using single variants. We examined the association between the 8.1 ancestral haplotype and delayed colonization in Cystic Fibrosis, postulating that downregulation of RNF5 expression is the likely causal mechanism. Our study provides insights into the genetic architecture of the MHC region and pinpoints disease associations that are due to differential expression of HLA genes and non-HLA genes.

Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits

The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes.