Search

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

SARS-CoV-2 infection of human cells is initiated by the binding of the viral Spike protein to its cell-surface receptor ACE2. We conducted a targeted CRISPRi screen to uncover druggable pathways controlling Spike protein binding to human cells. Here we show that the protein BRD2 is required for ACE2 transcription in human lung epithelial cells and cardiomyocytes, and BRD2 inhibitors currently evaluated in clinical trials potently block endogenous ACE2 expression and SARS-CoV-2 infection of human cells, including those of human nasal epithelia. Moreover, pharmacological BRD2 inhibition with the drug ABBV-744 inhibited SARS-CoV-2 replication in Syrian hamsters. We also found that BRD2 controls transcription of several other genes induced upon SARS-CoV-2 infection, including the interferon response, which in turn regulates the antiviral response. Together, our results pinpoint BRD2 as a potent and essential regulator of the host response to SARS-CoV-2 infection and highlight the potential of BRD2 as a therapeutic target for COVID-19.

Although coronavirus disease 2019 (COVID-19) causes cardiac dysfunction in up to 25% of patients, its pathogenesis remains unclear. Exposure of human induced pluripotent stem cell (iPSC)-derived heart cells to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) revealed productive infection and robust transcriptomic and morphological signatures of damage, particularly in cardiomyocytes. Transcriptomic disruption of structural genes corroborates adverse morphologic features, which included a distinct pattern of myofibrillar fragmentation and nuclear disruption. Human autopsy specimens from patients with COVID-19 reflected similar alterations, particularly sarcomeric fragmentation. These notable cytopathic features in cardiomyocytes provide insights into SARS-CoV-2-induced cardiac damage, offer a platform for discovery of potential therapeutics, and raise concerns about the long-term consequences of COVID-19 in asymptomatic and severe cases.

Abstract Background Many complementary solutions are available for the identifier mapping problem. This creates an opportunity for bioinformatics tool developers. Tools can be made to flexibly support multiple mapping services or mapping services could be combined to get broader coverage. This approach requires an interface layer between tools and mapping services. Results Here we present BridgeDb, a software framework for gene, protein and metabolite identifier mapping. This framework provides a standardized interface layer through which bioinformatics tools can be connected to different identifier mapping services. This approach makes it easier for tool developers to support identifier mapping. Mapping services can be combined or merged to support multi-omics experiments or to integrate custom microarray annotations. BridgeDb provides its own ready-to-go mapping services, both in webservice and local database forms. However, the framework is intended for customization and adaptation to any identifier mapping service. BridgeDb has already been integrated into several bioinformatics applications. Conclusion By uncoupling bioinformatics tools from mapping services, BridgeDb improves capability and flexibility of those tools. All described software is open source and available at http://www.bridgedb.org.

Abstract Background Biological pathways are a useful abstraction of biological concepts, and software tools to deal with pathway diagrams can help biological research. PathVisio is a new visualization tool for biological pathways that mimics the popular GenMAPP tool with a completely new Java implementation that allows better integration with other open source projects. The GenMAPP MAPP file format is replaced by GPML, a new XML file format that provides seamless exchange of graphical pathway information among multiple programs. Results PathVisio can be combined with other bioinformatics tools to open up three possible uses: visual compilation of biological knowledge, interpretation of high-throughput expression datasets, and computational augmentation of pathways with interaction information. PathVisio is open source software and available at http://www.pathvisio.org. Conclusion PathVisio is a graphical editor for biological pathways, with flexibility and ease of use as primary goals.

Mutations in the cardiac splicing factor RBM20 lead to malignant dilated cardiomyopathy (DCM). To understand the mechanism of RBM20-associated DCM, we engineered isogenic iPSCs with DCM-associated missense mutations in RBM20 as well as RBM20 knockout (KO) iPSCs. iPSC-derived engineered heart tissues made from these cell lines recapitulate contractile dysfunction of RBM20-associated DCM and reveal greater dysfunction with missense mutations than KO. Analysis of RBM20 RNA binding by eCLIP reveals a gain-of-function preference of mutant RBM20 for 3' UTR sequences that are shared with amyotrophic lateral sclerosis (ALS) and processing-body associated RNA binding proteins (FUS, DDX6). Deep RNA sequencing reveals that the RBM20 R636S mutant has unique gene, splicing, polyadenylation and circular RNA defects that differ from RBM20 KO. Super-resolution microscopy verifies that mutant RBM20 maintains very limited nuclear localization potential; rather, the mutant protein associates with cytoplasmic processing bodies (DDX6) under basal conditions, and with stress granules (G3BP1) following acute stress. Taken together, our results highlight a pathogenic mechanism in cardiac disease through splicing-dependent and -independent pathways.

Abstract Background Microarray technologies have evolved rapidly, enabling biologists to quantify genome-wide levels of gene expression, alternative splicing, and sequence variations for a variety of species. Analyzing and displaying these data present a significant challenge. Pathway-based approaches for analyzing microarray data have proven useful for presenting data and for generating testable hypotheses. Results To address the growing needs of the microarray community we have released version 2 of Gene Map Annotator and Pathway Profiler (GenMAPP), a new GenMAPP database schema, and integrated resources for pathway analysis. We have redesigned the GenMAPP database to support multiple gene annotations and species as well as custom species database creation for a potentially unlimited number of species. We have expanded our pathway resources by utilizing homology information to translate pathway content between species and extending existing pathways with data derived from conserved protein interactions and coexpression. We have implemented a new mode of data visualization to support analysis of complex data, including time-course, single nucleotide polymorphism (SNP), and splicing. GenMAPP version 2 also offers innovative ways to display and share data by incorporating HTML export of analyses for entire sets of pathways as organized web pages. Conclusion GenMAPP version 2 provides a means to rapidly interrogate complex experimental data for pathway-level changes in a diverse range of organisms.

Abstract Background The physiological regulation of G protein-coupled receptors, through desensitization and internalization, modulates the length of the receptor signal and may influence the development of tolerance and dependence in response to chronic drug treatment. To explore the importance of receptor regulation, we engineered a series of Gi-coupled receptors that differ in signal length, degree of agonist-induced internalization, and ability to induce adenylyl cyclase superactivation. All of these receptors, based on the kappa opioid receptor, were modified to be receptors activated solely by synthetic ligands (RASSLs). This modification allows us to compare receptors that have the same ligands and effectors, but differ only in desensitization and internalization. Results Removal of phosphorylation sites in the C-terminus of the RASSL resulted in a mutant that was resistant to internalization and less prone to desensitization. Replacement of the C-terminus of the RASSL with the corresponding portion of the mu opioid receptor eliminated the induction of AC superactivation, without disrupting agonist-induced desensitization or internalization. Surprisingly, removal of phosphorylation sites from this chimera resulted in a receptor that is constitutively internalized, even in the absence of agonist. However, the receptor still signals and desensitizes in response to agonist, indicating normal G-protein coupling and partial membrane expression. Conclusions These studies reveal that internalization, desensitization and adenylyl cyclase superactivation, all processes that decrease chronic Gi-receptor signals, are independently regulated. Furthermore, specific mutations can radically alter superactivation or internalization without affecting the efficacy of acute Gi signaling. These mutant RASSLs will be useful for further elucidating the temporal dynamics of the signaling of G protein-coupled receptors in vitro and in vivo.

Although COVID-19 causes cardiac dysfunction in up to 25% of patients, its pathogenesis remains unclear. Exposure of human iPSC-derived heart cells to SARS-CoV-2 revealed productive infection and robust transcriptomic and morphological signatures of damage, particularly in cardiomyocytes. Transcriptomic disruption of structural proteins corroborated adverse morphologic features, which included a distinct pattern of myofibrillar fragmentation and numerous iPSC-cardiomyocytes lacking nuclear DNA. Human autopsy specimens from COVID-19 patients displayed similar sarcomeric disruption, as well as cardiomyocytes without DNA staining. These striking cytopathic features provide new insights into SARS-CoV-2 induced cardiac damage, offer a platform for discovery of potential therapeutics, and raise serious concerns about the long-term consequences of COVID-19.

The excision of genomic sequences using paired CRISPR-Cas nucleases is a powerful tool to study gene function, create disease models and holds promise for therapeutic gene editing. However, our understanding of the factors that favor efficient excision is limited by the lack of a rapid, accurate measurement of DNA excision outcomes that is free of amplification bias. Here, we introduce ddXR (droplet digital PCR eXcision Reporter), a method that enables the accurate and sensitive detection of excisions and inversions independent of length. The method can be completed in a few hours without the need for next-generation sequencing. The ddXR method uncovered unexpectedly high rates of large (> 20 kb) excisions and inversions, while also revealing a surprisingly low dependence on linear distance, up to 170 kb. We further modified the method to measure precise repair of excision junctions and allele-specific excision, with important implications for disease modeling and therapeutic gene editing.

Previous studies have suggested that the loss of the translation initiation factor eIF4G1 homolog NAT1 induces excessive self-renewability of naive pluripotent stem cells (PSCs); yet the role of NAT1 in the self-renewal and differentiation of primed PSCs is still unclear. Here, we generate a conditional knockout of NAT1 in primed PSCs and use the cells for the functional analyses of NAT1. Our results show that NAT1 is required for the self-renewal and neural differentiation of primed PSCs. In contrast, NAT1 deficiency in naive pluripotency attenuates the differentiation to all cell types. We also find that NAT1 is involved in efficient protein expression of an RNA uridyltransferase, TUT7. TUT7 is involved in the neural differentiation of primed PSCs via the regulation of human endogenous retrovirus accumulation. These data demonstrate the essential roles of NAT1 and TUT7 in the precise transition of stem cell fate.

DISCOVER-seq (discovery of in situ Cas off-targets and verification by sequencing) is a broadly applicable approach for unbiased CRISPR-Cas off-target identification in cells and tissues. It leverages the recruitment of DNA repair factors to double-strand breaks (DSBs) after genome editing with CRISPR nucleases. Here, we describe a detailed experimental protocol and analysis pipeline with which to perform DISCOVER-seq. The principle of this method is to track the precise recruitment of MRE11 to DSBs by chromatin immunoprecipitation followed by next-generation sequencing. A customized open-source bioinformatics pipeline, BLENDER (blunt end finder), then identifies off-target sequences genome wide. DISCOVER-seq is capable of finding and measuring off-targets in primary cells and in situ. The two main advantages of DISCOVER-seq are (i) low false-positive rates because DNA repair enzyme binding is required for genome edits to occur and (ii) its applicability to a wide variety of systems, including patient-derived cells and animal models. The whole protocol, including the analysis, can be completed within 2 weeks.

Repair of double strand DNA breaks (DSBs) can result in gene disruption or gene modification via homology directed repair (HDR) from donor DNA. Altering cellular responses to DSBs may rebalance editing outcomes towards HDR and away from other repair outcomes. Here, we utilize a pooled CRISPR screen to define host cell involvement in HDR between a Cas9 DSB and a plasmid double stranded donor DNA (dsDonor). We find that the Fanconi Anemia (FA) pathway is required for dsDonor HDR and that other genes act to repress HDR. Small molecule inhibition of one of these repressors, CDC7, by XL413 and other inhibitors increases the efficiency of HDR by up to 3.5 fold in many contexts, including primary T cells. XL413 stimulates HDR during a reversible slowing of S-phase that is unexplored for Cas9-induced HDR. We anticipate that XL413 and other such rationally developed inhibitors will be useful tools for gene modification.

Cardiac tissues are able to adjust their contractile behavior to adapt to the local mechanical environment. Nonuniformity of the native tissue mechanical properties contributes to the development of heart dysfunctions, yet the current in vitro cardiac tissue models often fail to recapitulate the mechanical nonuniformity. To address this issue, a 3D cardiac microtissue model is developed with engineered mechanical nonuniformity, enabled by 3D-printed hybrid matrices composed of fibers with different diameters. When escalating the complexity of tissue mechanical environments, cardiac microtissues start to develop maladaptive hypercontractile phenotypes, demonstrated in both contractile motion analysis and force-power analysis. This novel hybrid system could potentially facilitate the establishment of "pathologically-inspired" cardiac microtissue models for deeper understanding of heart pathology due to nonuniformity of the tissue mechanical environment.

The CRISPR/Cas system is a highly specific genome editing tool capable of distinguishing alleles differing by even a single base pair. Target sites might carry genetic variations that are not distinguishable by sgRNA designing tools based on one reference genome. AlleleAnalyzer is an open-source software that incorporates single-nucleotide variants and short insertions and deletions to design sgRNAs for precisely editing 1 or multiple haplotypes of a sequenced genome, currently supporting 11 Cas proteins. It also leverages patterns of shared genetic variation to optimize sgRNA design for different human populations. AlleleAnalyzer is available at https://github.com/keoughkath/AlleleAnalyzer .

Human pluripotent stem cells (hPSCs) have the intrinsic ability to self-organize into complex multicellular organoids that recapitulate many aspects of tissue development. However, robustly directing morphogenesis of hPSC-derived organoids requires novel approaches to accurately control self-directed pattern formation. Here, we combined genetic engineering with computational modeling, machine learning, and mathematical pattern optimization to create a data-driven approach to control hPSC self-organization by knock down of genes previously shown to affect stem cell colony organization, CDH1 and ROCK1. Computational replication of the in vitro system in silico using an extended cellular Potts model enabled machine learning-driven optimization of parameters that yielded emergence of desired patterns. Furthermore, in vitro the predicted experimental parameters quantitatively recapitulated the in silico patterns. These results demonstrate that morphogenic dynamics can be accurately predicted through model-driven exploration of hPSC behaviors via machine learning, thereby enabling spatial control of multicellular patterning to engineer human organoids and tissues. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

For long QT syndrome (LQTS), recent progress in genome-sequencing technologies enabled the identification of rare genomic variants with diagnostic, prognostic, and therapeutic implications. However, pathogenic stratification of the identified variants remains challenging, especially in variants of uncertain significance. This study aimed to propose a phenotypic cell-based diagnostic assay for identifying LQTS to recognize pathogenic variants in a high-throughput manner suitable for screening. We investigated the response of LQT2-induced pluripotent stem cell (iPSC)-derived cardiomyocytes (iPSC-CMs) following IKr blockade using a multi-electrode array, finding that the response to IKr blockade was significantly smaller than in Control-iPSC-CMs. Furthermore, we found that LQT1-iPSC-CMs and LQT3-iPSC-CMs could be distinguished from Control-iPSC-CMs by IKs blockade and INa blockade, respectively. This strategy might be helpful in compensating for the shortcomings of genetic testing of LQTS patients.

Mutations in gene regulatory elements have been associated with a wide range of complex neuropsychiatric disorders. However, due to their cell-type specificity and difficulties in characterizing their regulatory targets, the ability to identify causal genetic variants has remained limited. To address these constraints, we perform an integrative analysis of chromatin interactions, open chromatin regions and transcriptomes using promoter capture Hi-C, assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing, respectively, in four functionally distinct neural cell types: induced pluripotent stem cell (iPSC)-induced excitatory neurons and lower motor neurons, iPSC-derived hippocampal dentate gyrus-like neurons and primary astrocytes. We identify hundreds of thousands of long-range cis-interactions between promoters and distal promoter-interacting regions, enabling us to link regulatory elements to their target genes and reveal putative processes that are dysregulated in disease. Finally, we validate several promoter-interacting regions by using clustered regularly interspaced short palindromic repeats (CRISPR) techniques in human excitatory neurons, demonstrating that CDK5RAP3, STRAP and DRD2 are transcriptionally regulated by physically linked enhancers.

CRISPR-Cas genome editing induces targeted DNA damage but can also affect off-target sites. Current off-target discovery methods work using purified DNA or specific cellular models but are incapable of direct detection in vivo. We developed DISCOVER-Seq (discovery of in situ Cas off-targets and verification by sequencing), a universally applicable approach for unbiased off-target identification that leverages the recruitment of DNA repair factors in cells and organisms. Tracking the precise recruitment of MRE11 uncovers the molecular nature of Cas activity in cells with single-base resolution. DISCOVER-Seq works with multiple guide RNA formats and types of Cas enzymes, allowing characterization of new editing tools. Off-targets can be identified in cell lines and patient-derived induced pluripotent stem cells and during adenoviral editing of mice, paving the way for in situ off-target discovery within individual patient genotypes during therapeutic genome editing.

MESP1 is a key transcription factor in development of early cardiovascular tissue and it is required for induction of the cardiomyocyte (CM) gene expression program, but its role in vascular development is unclear. Here, we used inducible CRISPRi knock-down of MESP1 to analyze the molecular processes of the early differentiation stages of human induced pluripotent stem cells into mesoderm and subsequently vascular progenitor cells. We found that expression of the mesodermal marker, BRACHYURY (encoded by T) was unaffected in MESP1 knock-down cells as compared to wild type cells suggesting timely movement through the primitive streak whereas another mesodermal marker MIXL1 was slightly, but significantly decreased. In contrast, the expression of the vascular cell surface marker KDR was decreased and CD31 and CD34 expression were substantially reduced in MESP1 knock-down cells supporting inhibition or delay of vascular specification. In addition, mRNA microarray data revealed several other altered gene expressions including the EMT regulating transcription factors SNAI1 and TWIST1, which were both significantly decreased indicating that MESP1 knock-down cells are less likely to undergo EMT during vascular progenitor differentiation. Our study demonstrates that while leaving primitive streak markers unaffected, MESP1 expression is required for timely vascular progenitor specification. Thus, MESP1 expression is essential for the molecular features of early CM, EC and VSMC lineage specification.

The integration of in vitro cardiac tissue models, human induced pluripotent stem cells (hiPSCs) and genome-editing tools allows for the enhanced interrogation of physiological phenotypes and recapitulation of disease pathologies. Here, using a cardiac tissue model consisting of filamentous three-dimensional matrices populated with cardiomyocytes derived from healthy wild-type (WT) hiPSCs (WT hiPSC-CMs) or isogenic hiPSCs deficient in the sarcomere protein cardiac myosin-binding protein C (MYBPC3-/- hiPSC-CMs), we show that the WT microtissues adapted to the mechanical environment with increased contraction force commensurate to matrix stiffness, whereas the MYBPC3-/- microtissues exhibited impaired force development kinetics regardless of matrix stiffness and deficient contraction force only when grown on matrices with high fibre stiffness. Under mechanical overload, the MYBPC3-/- microtissues had a higher degree of calcium transient abnormalities, and exhibited an accelerated decay of calcium dynamics as well as calcium desensitization, which accelerated when contracting against stiffer fibres. Our findings suggest that MYBPC3 deficiency and the presence of environmental stresses synergistically lead to contractile deficits in cardiac tissues.

Morphogenesis involves interactions of asymmetric cell populations to form complex multicellular patterns and structures comprised of distinct cell types. However, current methods to model morphogenic events lack control over cell-type co-emergence and offer little capability to selectively perturb specific cell subpopulations. Our in vitro system interrogates cell-cell interactions and multicellular organization within human induced pluripotent stem cell (hiPSC) colonies. We examined effects of induced mosaic knockdown of molecular regulators of cortical tension (ROCK1) and cell-cell adhesion (CDH1) with CRISPR interference. Mosaic knockdown of ROCK1 or CDH1 resulted in differential patterning within hiPSC colonies due to cellular self-organization, while retaining an epithelial pluripotent phenotype. Knockdown induction stimulates a transient wave of differential gene expression within the mixed populations that stabilized in coordination with observed self-organization. Mosaic patterning enables genetic interrogation of emergent multicellular properties, which can facilitate better understanding of the molecular pathways that regulate symmetry-breaking during morphogenesis.

Quantification of abnormal contractile motions of cardiac tissue has been a noteworthy challenge and significant limitation in assessing and classifying the drug-induced arrhythmias (i.e., Torsades de pointes). To overcome these challenges, researchers have taken advantage of computational image processing tools to measure contractile motion from cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs). However, the amplitude and frequency analysis of contractile motion waveforms does not produce sufficient information to objectively classify the degree of variations between two or more sets of cardiac contractile motions. In this paper, we generated contractile motion data from beating hiPSC-CMs using motion tracking software based on optical flow analysis, and then implemented a computational algorithm, phase space reconstruction (PSR), to derive parameters (embedding, regularity, and fractal dimensions) to further characterize the dynamic nature of the cardiac contractile motions. Application of drugs known to cause cardiac arrhythmia induced significant changes to these resultant dimensional parameters calculated from PSR analysis. Integrating this new computational algorithm with the existing analytical toolbox of cardiac contractile motions will allow us to expand current assessments of cardiac tissue physiology into an automated, high-throughput, and quantifiable manner which will allow more objective assessments of drug-induced proarrhythmias.

The creation of human induced pluripotent stem cells (hiPSCs) has provided an unprecedented opportunity to study tissue morphogenesis and organ development through 'organogenesis-in-a-dish'. Current approaches to cardiac organoid engineering rely on either direct cardiac differentiation from embryoid bodies (EBs) or generation of aligned cardiac tissues from predifferentiated cardiomyocytes from monolayer hiPSCs. To experimentally model early cardiac organogenesis in vitro, our protocol combines biomaterials-based cell patterning with stem cell organoid engineering. 3D cardiac microchambers are created from 2D hiPSC colonies; these microchambers approximate an early-development heart with distinct spatial organization and self-assembly. With proper training in photolithography microfabrication, maintenance of human pluripotent stem cells, and cardiac differentiation, a graduate student with guidance will likely be able to carry out this experimental protocol, which requires ∼3 weeks. We envisage that this in vitro model of human early heart development could serve as an embryotoxicity screening assay in drug discovery, regulation, and prescription for healthy fetal development. We anticipate that, when applied to hiPSC lines derived from patients with inherited diseases, this protocol can be used to study the disease mechanisms of cardiac malformations at an early stage of embryogenesis.

Genome editing holds great promise for experimental biology and potential clinical use. To successfully utilize genome editing, it is critical to sensitively detect and quantify its outcomes: homology-directed repair (HDR) and nonhomologous end joining (NHEJ). This has been difficult at endogenous gene loci and instead is frequently done using artificial reporter systems. Here, we describe a droplet digital PCR (ddPCR)-based method to simultaneously measure HDR and NHEJ at endogenous gene loci. This highly sensitive and quantitative method may significantly contribute to a better understanding of DNA repair mechanisms underlying genome editing and to the improvement of genome editing technology by allowing for efficient and systematic testing of many genome editing conditions in parallel.

Molecular chaperones regulate quality control in the human proteome, pathways that have been implicated in many diseases, including heart failure. Mutations in the BAG3 gene, which encodes a co-chaperone protein, have been associated with heart failure due to both inherited and sporadic dilated cardiomyopathy. Familial BAG3 mutations are autosomal dominant and frequently cause truncation of the coding sequence, suggesting a heterozygous loss-of-function mechanism. However, heterozygous knockout of the murine BAG3 gene did not cause a detectable phenotype. To model BAG3 cardiomyopathy in a human system, we generated an isogenic series of human induced pluripotent stem cells (iPSCs) with loss-of-function mutations in BAG3. Heterozygous BAG3 mutations reduced protein expression, disrupted myofibril structure, and compromised contractile function in iPSC-derived cardiomyocytes (iPS-CMs). BAG3-deficient iPS-CMs were particularly sensitive to further myofibril disruption and contractile dysfunction upon exposure to proteasome inhibitors known to cause cardiotoxicity. We performed affinity tagging of the endogenous BAG3 protein and mass spectrometry proteomics to further define the cardioprotective chaperone complex that BAG3 coordinates in the human heart. Our results establish a model for evaluating protein quality control pathways in human cardiomyocytes and their potential as therapeutic targets and susceptibility factors for cardiac drug toxicity.

RationaleDuring each beat, cardiac myocytes (CMs) generate the mechanical output necessary for heart function through contractile mechanisms that involve shortening of sarcomeres along myofibrils. Human-induced pluripotent stem cells (hiPSCs) can be differentiated into CMs (hiPSC-CMs) that model cardiac contractile mechanical output more robustly when micropatterned into physiological shapes. Quantifying the mechanical output of these cells enables us to assay cardiac activity in a dish.ObjectiveWe sought to develop a computational platform that integrates analytic approaches to quantify the mechanical output of single micropatterned hiPSC-CMs from microscopy videos.Methods and resultsWe micropatterned single hiPSC-CMs on deformable polyacrylamide substrates containing fluorescent microbeads. We acquired videos of single beating cells, of microbead displacement during contractions, and of fluorescently labeled myofibrils. These videos were independently analyzed to obtain parameters that capture the mechanical output of the imaged single cells. We also developed novel methods to quantify sarcomere length from videos of moving myofibrils and to analyze loss of synchronicity of beating in cells with contractile defects. We tested this computational platform by detecting variations in mechanical output induced by drugs and in cells expressing low levels of myosin-binding protein C.ConclusionsOur method can measure the cardiac function of single micropatterned hiPSC-CMs and determine contractile parameters that can be used to elucidate mechanisms that underlie variations in CM function. This platform will be amenable to future studies of the effects of mutations and drugs on cardiac function.

The human genome produces thousands of long noncoding RNAs (lncRNAs)-transcripts >200 nucleotides long that do not encode proteins. Although critical roles in normal biology and disease have been revealed for a subset of lncRNAs, the function of the vast majority remains untested. We developed a CRISPR interference (CRISPRi) platform targeting 16,401 lncRNA loci in seven diverse cell lines, including six transformed cell lines and human induced pluripotent stem cells (iPSCs). Large-scale screening identified 499 lncRNA loci required for robust cellular growth, of which 89% showed growth-modifying function exclusively in one cell type. We further found that lncRNA knockdown can perturb complex transcriptional networks in a cell type-specific manner. These data underscore the functional importance and cell type specificity of many lncRNAs.

The enteric nervous system (ENS) of the gastrointestinal tract controls many diverse functions, including motility and epithelial permeability. Perturbations in ENS development or function are common, yet there is no human model for studying ENS-intestinal biology and disease. We used a tissue-engineering approach with embryonic and induced pluripotent stem cells (PSCs) to generate human intestinal tissue containing a functional ENS. We recapitulated normal intestinal ENS development by combining human-PSC-derived neural crest cells (NCCs) and developing human intestinal organoids (HIOs). NCCs recombined with HIOs in vitro migrated into the mesenchyme, differentiated into neurons and glial cells and showed neuronal activity, as measured by rhythmic waves of calcium transients. ENS-containing HIOs grown in vivo formed neuroglial structures similar to a myenteric and submucosal plexus, had functional interstitial cells of Cajal and had an electromechanical coupling that regulated waves of propagating contraction. Finally, we used this system to investigate the cellular and molecular basis for Hirschsprung's disease caused by a mutation in the gene PHOX2B. This is, to the best of our knowledge, the first demonstration of human-PSC-derived intestinal tissue with a functional ENS and how this system can be used to study motility disorders of the human gastrointestinal tract.

Fibrodysplasia ossificans progressiva (FOP) patients carry a missense mutation in ACVR1 [617G > A (R206H)] that leads to hyperactivation of BMP-SMAD signaling. Contrary to a previous study, here we show that FOP fibroblasts showed an increased efficiency of induced pluripotent stem cell (iPSC) generation. This positive effect was attenuated by inhibitors of BMP-SMAD signaling (Dorsomorphin or LDN1931890) or transducing inhibitory SMADs (SMAD6 or SMAD7). In normal fibroblasts, the efficiency of iPSC generation was enhanced by transducing mutant ACVR1 (617G > A) or SMAD1 or adding BMP4 protein at early times during the reprogramming. In contrast, adding BMP4 at later times decreased iPSC generation. ID genes, transcriptional targets of BMP-SMAD signaling, were critical for iPSC generation. The BMP-SMAD-ID signaling axis suppressed p16/INK4A-mediated cell senescence, a major barrier to reprogramming. These results using patient cells carrying the ACVR1 R206H mutation reveal how cellular signaling and gene expression change during the reprogramming processes.

Tissue engineering approaches have the potential to increase the physiologic relevance of human iPS-derived cells, such as cardiomyocytes (iPS-CM). However, forming Engineered Heart Muscle (EHM) typically requires >1 million cells per tissue. Existing miniaturization strategies involve complex approaches not amenable to mass production, limiting the ability to use EHM for iPS-based disease modeling and drug screening. Micro-scale cardiospheres are easily produced, but do not facilitate assembly of elongated muscle or direct force measurements. Here we describe an approach that combines features of EHM and cardiospheres: Micro-Heart Muscle (μHM) arrays, in which elongated muscle fibers are formed in an easily fabricated template, with as few as 2,000 iPS-CM per individual tissue. Within μHM, iPS-CM exhibit uniaxial contractility and alignment, robust sarcomere assembly, and reduced variability and hypersensitivity in drug responsiveness, compared to monolayers with the same cellular composition. μHM mounted onto standard force measurement apparatus exhibited a robust Frank-Starling response to external stretch, and a dose-dependent inotropic response to the β-adrenergic agonist isoproterenol. Based on the ease of fabrication, the potential for mass production and the small number of cells required to form μHM, this system provides a potentially powerful tool to study cardiomyocyte maturation, disease and cardiotoxicology in vitro.

Developing technologies for efficient and scalable disruption of gene expression will provide powerful tools for studying gene function, developmental pathways, and disease mechanisms. Here, we develop clustered regularly interspaced short palindromic repeat interference (CRISPRi) to repress gene expression in human induced pluripotent stem cells (iPSCs). CRISPRi, in which a doxycycline-inducible deactivated Cas9 is fused to a KRAB repression domain, can specifically and reversibly inhibit gene expression in iPSCs and iPSC-derived cardiac progenitors, cardiomyocytes, and T lymphocytes. This gene repression system is tunable and has the potential to silence single alleles. Compared with CRISPR nuclease (CRISPRn), CRISPRi gene repression is more efficient and homogenous across cell populations. The CRISPRi system in iPSCs provides a powerful platform to perform genome-scale screens in a wide range of iPSC-derived cell types, dissect developmental pathways, and model disease.

Precise genome-editing relies on the repair of sequence-specific nuclease-induced DNA nicking or double-strand breaks (DSBs) by homology-directed repair (HDR). However, nonhomologous end-joining (NHEJ), an error-prone repair, acts concurrently, reducing the rate of high-fidelity edits. The identification of genome-editing conditions that favor HDR over NHEJ has been hindered by the lack of a simple method to measure HDR and NHEJ directly and simultaneously at endogenous loci. To overcome this challenge, we developed a novel, rapid, digital PCR-based assay that can simultaneously detect one HDR or NHEJ event out of 1,000 copies of the genome. Using this assay, we systematically monitored genome-editing outcomes of CRISPR-associated protein 9 (Cas9), Cas9 nickases, catalytically dead Cas9 fused to FokI, and transcription activator-like effector nuclease at three disease-associated endogenous gene loci in HEK293T cells, HeLa cells, and human induced pluripotent stem cells. Although it is widely thought that NHEJ generally occurs more often than HDR, we found that more HDR than NHEJ was induced under multiple conditions. Surprisingly, the HDR/NHEJ ratios were highly dependent on gene locus, nuclease platform, and cell type. The new assay system, and our findings based on it, will enable mechanistic studies of genome-editing and help improve genome-editing technology.

Neural crest cells (NCC) are stem cells that generate different lineages, including neuroendocrine, melanocytic, cartilage, and bone. The differentiation potential of NCC varies according to the level from which cells emerge along the neural tube. For example, only anterior "cranial" NCC form craniofacial bone, whereas solely posterior "trunk" NCC contribute to sympathoadrenal cells. Importantly, the isolation of human fetal NCC carries ethical and scientific challenges, as NCC induction typically occur before pregnancy is detectable. As a result, current knowledge of NCC biology derives primarily from non-human organisms. Important differences between human and non-human NCC, such as expression of HNK1 in human but not mouse NCC, suggest a need to study human NCC directly. Here, we demonstrate that current protocols to differentiate human pluripotent stem cells (PSC) to NCC are biased toward cranial NCC. Addition of retinoic acid drove trunk-related markers and HOX genes characteristic of a posterior identity. Subsequent treatment with bone morphogenetic proteins (BMPs) enhanced differentiation to sympathoadrenal cells. Our approach provides methodology for detailed studies of human NCC, and clarifies roles for retinoids and BMPs in the differentiation of human PSC to trunk NCC and to sympathoadrenal lineages.

Human pluripotent stem cells (hPS cells) are rapidly emerging as a powerful tool for biomedical discovery. The advent of human induced pluripotent stem cells (hiPS cells) with human embryonic stem (hES)-cell-like properties has led to hPS cells with disease-specific genetic backgrounds for in vitro disease modeling and drug discovery as well as mechanistic and developmental studies. To fully realize this potential, it will be necessary to modify the genome of hPS cells with precision and flexibility. Pioneering experiments utilizing site-specific double-strand break (DSB)-mediated genome engineering tools, including zinc finger nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs), have paved the way to genome engineering in previously recalcitrant systems such as hPS cells. However, these methods are technically cumbersome and require significant expertise, which has limited adoption. A major recent advance involving the clustered regularly interspaced short palindromic repeats (CRISPR) endonuclease has dramatically simplified the effort required for genome engineering and will likely be adopted widely as the most rapid and flexible system for genome editing in hPS cells. In this unit, we describe commonly practiced methods for CRISPR endonuclease genomic editing of hPS cells into cell lines containing genomes altered by insertion/deletion (indel) mutagenesis or insertion of recombinant genomic DNA.

Tissue morphogenesis and organ formation are the consequences of biochemical and biophysical cues that lead to cellular spatial patterning in development. To model such events in vitro, we use PEG-patterned substrates to geometrically confine human pluripotent stem cell colonies and spatially present mechanical stress. Modulation of the WNT/β-catenin pathway promotes spatial patterning via geometric confinement of the cell condensation process during epithelial-mesenchymal transition, forcing cells at the perimeter to express an OCT4+ annulus, which is coincident with a region of higher cell density and E-cadherin expression. The biochemical and biophysical cues synergistically induce self-organizing lineage specification and creation of a beating human cardiac microchamber confined by the pattern geometry. These highly defined human cardiac microchambers can be used to study aspects of embryonic spatial patterning, early cardiac development and drug-induced developmental toxicity.

Contractile motion is the simplest metric of cardiomyocyte health in vitro, but unbiased quantification is challenging. We describe a rapid automated method, requiring only standard video microscopy, to analyze the contractility of human-induced pluripotent stem cell-derived cardiomyocytes (iPS-CM). New algorithms for generating and filtering motion vectors combined with a newly developed isogenic iPSC line harboring genetically encoded calcium indicator, GCaMP6f, allow simultaneous user-independent measurement and analysis of the coupling between calcium flux and contractility. The relative performance of these algorithms, in terms of improving signal to noise, was tested. Applying these algorithms allowed analysis of contractility in iPS-CM cultured over multiple spatial scales from single cells to three-dimensional constructs. This open source software was validated with analysis of isoproterenol response in these cells, and can be applied in future studies comparing the drug responsiveness of iPS-CM cultured in different microenvironments in the context of tissue engineering.

NANOG (from Irish mythology Tír na nÓg) transcription factor plays a central role in maintaining pluripotency, cooperating with OCT4 (also known as POU5F1 or OCT3/4), SOX2, and other pluripotency factors. Although the physiological roles of the NANOG protein have been extensively explored, biochemical and biophysical properties in relation to its structural analysis are poorly understood. Here we determined the crystal structure of the human NANOG homeodomain (hNANOG HD) bound to an OCT4 promoter DNA, which revealed amino acid residues involved in DNA recognition that are likely to be functionally important. We generated a series of hNANOG HD alanine substitution mutants based on the protein-DNA interaction and evolutionary conservation and determined their biological activities. Some mutant proteins were less stable, resulting in loss or decreased affinity for DNA binding. Overexpression of the orthologous mouse NANOG (mNANOG) mutants failed to maintain self-renewal of mouse embryonic stem cells without leukemia inhibitory factor. These results suggest that these residues are critical for NANOG transcriptional activity. Interestingly, one mutant, hNANOG L122A, conversely enhanced protein stability and DNA-binding affinity. The mNANOG L122A, when overexpressed in mouse embryonic stem cells, maintained their expression of self-renewal markers even when retinoic acid was added to forcibly drive differentiation. When overexpressed in epiblast stem cells or human induced pluripotent stem cells, the L122A mutants enhanced reprogramming into ground-state pluripotency. These findings demonstrate that structural and biophysical information on key transcriptional factors provides insights into the manipulation of stem cell behaviors and a framework for rational protein engineering.

We present a non-invasive method to characterize the function of pluripotent stem-cell-derived cardiomyocytes based on video microscopy and image analysis. The platform, called Pulse, generates automated measurements of beating frequency, beat duration, amplitude, and beat-to-beat variation based on motion analysis of phase-contrast images captured at a fast frame rate. Using Pulse, we demonstrate recapitulation of drug effects in stem-cell-derived cardiomyocytes without the use of exogenous labels and show that our platform can be used for high-throughput cardiotoxicity drug screening and studying physiologically relevant phenotypes.

Drug discovery and development are hampered by high failure rates attributed to the reliance on non-human animal models employed during safety and efficacy testing. A fundamental problem in this inefficient process is that non-human animal models cannot adequately represent human biology. Thus, there is an urgent need for high-content in vitro systems that can better predict drug-induced toxicity. Systems that predict cardiotoxicity are of uppermost significance, as approximately one third of safety-based pharmaceutical withdrawals are due to cardiotoxicty. Here, we present a cardiac microphysiological system (MPS) with the attributes required for an ideal in vitro system to predict cardiotoxicity: i) cells with a human genetic background; ii) physiologically relevant tissue structure (e.g. aligned cells); iii) computationally predictable perfusion mimicking human vasculature; and, iv) multiple modes of analysis (e.g. biological, electrophysiological, and physiological). Our MPS is able to keep human induced pluripotent stem cell derived cardiac tissue viable and functional over multiple weeks. Pharmacological studies using the cardiac MPS show half maximal inhibitory/effective concentration values (IC₅₀/EC₅₀) that are more consistent with the data on tissue scale references compared to cellular scale studies. We anticipate the widespread adoption of MPSs for drug screening and disease modeling.

Astrocytes express a variety of G protein-coupled receptors and might influence cognitive functions, such as learning and memory. However, the roles of astrocytic Gs-coupled receptors in cognitive function are not known. We found that humans with Alzheimer's disease (AD) had increased levels of the Gs-coupled adenosine receptor A2A in astrocytes. Conditional genetic removal of these receptors enhanced long-term memory in young and aging mice and increased the levels of Arc (also known as Arg3.1), an immediate-early gene that is required for long-term memory. Chemogenetic activation of astrocytic Gs-coupled signaling reduced long-term memory in mice without affecting learning. Like humans with AD, aging mice expressing human amyloid precursor protein (hAPP) showed increased levels of astrocytic A2A receptors. Conditional genetic removal of these receptors enhanced memory in aging hAPP mice. Together, these findings establish a regulatory role for astrocytic Gs-coupled receptors in memory and suggest that AD-linked increases in astrocytic A2A receptor levels contribute to memory loss.