Alternative splicing is an important mechanism for increasing protein diversity. However, its functional effects are largely unknown. Here, we present our new software workflow composed of the open-source application AltAnalyze and the Cytoscape plugin DomainGraph. Both programs provide an intuitive and comprehensive end-to-end solution for the analysis and visualization of alternative splicing data from Affymetrix Exon and Gene Arrays at the level of proteins, domains, microRNA binding sites, molecular interactions and pathways. Our software tools include easy-to-use graphical user interfaces, rigorous statistical methods (FIRMA, MiDAS and DABG filtering) and do not require prior knowledge of exon array analysis or programming. They provide new methods for automatic interpretation and visualization of the effects of alternative exon inclusion on protein domain composition and microRNA binding sites. These data can be visualized together with affected pathways and gene or protein interaction networks, allowing a straightforward identification of potential biological effects due to alternative splicing at different levels of granularity. Our programs are available at http://www.altanalyze.org and http://www.domaingraph.de. These websites also include extensive documentation, tutorials and sample data.

SNPLogic (http://www.snplogic.org) brings together single nucleotide polymorphism (SNP) information from numerous sources to provide a comprehensive SNP selection, annotation and prioritization system for design and analysis of genotyping projects. SNPLogic integrates information about the genetic context of SNPs (gene, chromosomal region, functional location, haplotypes tags and overlap with transcription factor binding sites, splicing sites, miRNAs and evolutionarily conserved regions), genotypic data (allele frequencies per population and validation method), coverage of commercial arrays (ParAllele, Affymetrix and Illumina), functional predictions (modeled on structure and sequence) and connections or established associations (biological pathways, gene ontology terms and OMIM disease terms). The SNPLogic web interface facilitates construction and annotation of user-defined SNP lists that can be saved, shared and exported. Thus, SNPLogic can be used to identify and prioritize candidate SNPs, assess custom and commercial arrays panels and annotate new SNP data with publicly available information. We have found integration of SNP annotation in the context of pathway information and functional prediction scores to be a powerful approach to the analysis and interpretation of SNP-disease association data.

Long-QT syndrome mutations can cause syncope and sudden death by prolonging the cardiac action potential (AP). Ion channels affected by mutations are various, and the influences of cellular calcium cycling on LQTS cardiac events are unknown. To better understand LQTS arrhythmias, we performed current-clamp and intracellular calcium ([Ca(2+)]i) measurements on cardiomyocytes differentiated from patient-derived induced pluripotent stem cells (iPS-CM). In myocytes carrying an LQT2 mutation (HERG-A422T), APs and [Ca(2+)]i transients were prolonged in parallel. APs were abbreviated by nifedipine exposure and further lengthened upon releasing intracellularly stored Ca(2+). Validating this model, control iPS-CM treated with HERG-blocking drugs recapitulated the LQT2 phenotype. In LQT3 iPS-CM, expressing NaV1.5-N406K, APs and [Ca(2+)]i transients were markedly prolonged. AP prolongation was sensitive to tetrodotoxin and to inhibiting Na(+)-Ca(2+) exchange. These results suggest that LQTS mutations act partly on cytosolic Ca(2+) cycling, potentially providing a basis for functionally targeted interventions regardless of the specific mutation site.

The objective of this study was to determine the role of A-Kinase Anchoring Protein (AKAP)-Lbc in the development of heart failure, by investigating AKAP-Lbc-protein kinase D1 (PKD1) signaling in vivo in cardiac hypertrophy. Using a gene-trap mouse expressing a truncated version of AKAP-Lbc (due to disruption of the endogenous AKAP-Lbc gene), that abolishes PKD1 interaction with AKAP-Lbc (AKAP-Lbc-ΔPKD), we studied two mouse models of pathological hypertrophy: i) angiotensin (AT-II) and phenylephrine (PE) infusion and ii) transverse aortic constriction (TAC)-induced pressure overload. Our results indicate that AKAP-Lbc-ΔPKD mice exhibit an accelerated progression to cardiac dysfunction in response to AT-II/PE treatment and TAC. AKAP-Lbc-ΔPKD mice display attenuated compensatory cardiac hypertrophy, increased collagen deposition and apoptosis, compared to wild-type (WT) control littermates. Mechanistically, reduced levels of PKD1 activation are observed in AKAP-Lbc-ΔPKD mice compared to WT mice, resulting in diminished phosphorylation of histone deacetylase 5 (HDAC5) and decreased hypertrophic gene expression. This is consistent with a reduced compensatory hypertrophy phenotype leading to progression of heart failure in AKAP-Lbc-ΔPKD mice. Overall, our data demonstrates a critical in vivo role for AKAP-Lbc-PKD1 signaling in the development of compensatory hypertrophy to enhance cardiac performance in response to TAC-induced pressure overload and neurohumoral stimulation by AT-II/PE treatment.

The human nuclear factor related to kappa-B-binding protein (NFRKB) is a 1299-residue protein that is a component of the metazoan INO80 complex involved in chromatin remodeling, transcription regulation, DNA replication and DNA repair. Although full length NFRKB is predicted to be around 65% disordered, comparative sequence analysis identified several potentially structured sections in the N-terminal region of the protein. These regions were targeted for crystallographic studies, and the structure of one of these regions spanning residues 370-495 was determined using the JCSG high-throughput structure determination pipeline. The structure reveals a novel, mostly helical domain reminiscent of the winged-helix fold typically involved in DNA binding. However, further analysis shows that this domain does not bind DNA, suggesting it may belong to a small group of winged-helix domains involved in protein-protein interactions.

The molecular and functional diversity of G protein-coupled receptors is essential to many physiological processes. However, this diversity presents a significant challenge to understanding the G protein-mediated signaling events that underlie a specific physiological response. To increase our understanding of these processes, we sought to gain control of the timing and specificity of G(s) signaling in vivo. We used naturally occurring human mutations to develop two G(s)-coupled engineered receptors that respond solely to a synthetic ligand (RASSLs). Our G(s)-coupled RASSLs are based on the melanocortin-4 receptor, a centrally expressed receptor that plays an important role in the regulation of body weight. These RASSLs are not activated by the endogenous hormone alpha-melanocyte-stimulating hormone but respond potently to a selective synthetic ligand, tetrahydroisoquinoline. The RASSL variants reported here differ in their intrinsic basal activities, allowing the separation of the effects of basal signaling from ligand-mediated activation of the G(s) pathway in vivo. These RASSLs can be used to activate G(s) signaling in any tissue, but would be particularly useful for analyzing downstream events that mediate body weight regulation in mice. Our study also demonstrates the use of human genetic variation for protein engineering.

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.

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.

Precise editing of human genomes in pluripotent stem cells by homology-driven repair of targeted nuclease-induced cleavage has been hindered by the difficulty of isolating rare clones. We developed an efficient method to capture rare mutational events, enabling isolation of mutant lines with single-base substitutions without antibiotic selection. This method facilitates efficient induction or reversion of mutations associated with human disease in isogenic human induced pluripotent stem cells.

WikiPathways provides a collaborative platform for creating, updating, and sharing pathway diagrams and serves as an example of content curation by the biology community.

G protein-coupled receptors (GPCRs) signal through a limited number of G-protein pathways and play crucial roles in many biological processes. Studies of their in vivo functions have been hampered by the molecular and functional diversity of GPCRs and the paucity of ligands with specific signaling effects. To better compare the effects of activating different G-protein signaling pathways through ligand-induced or constitutive signaling, we developed a new series of RASSLs (receptors activated solely by synthetic ligands) that activate different G-protein signaling pathways. These RASSLs are based on the human 5-HT(4b) receptor, a GPCR with high constitutive G(s) signaling and strong ligand-induced G-protein activation of the G(s) and G(s/q) pathways. The first receptor in this series, 5-HT(4)-D(100)A or Rs1 (RASSL serotonin 1), is not activated by its endogenous agonist, serotonin, but is selectively activated by the small synthetic molecules GR113808, GR125487, and RO110-0235. All agonists potently induced G(s) signaling, but only a few (e.g., zacopride) also induced signaling via the G(q) pathway. Zacopride-induced G(q) signaling was enhanced by replacing the C-terminus of Rs1 with the C-terminus of the human 5-HT(2C) receptor. Additional point mutations (D(66)A and D(66)N) blocked constitutive G(s) signaling and lowered ligand-induced G(q) signaling. Replacing the third intracellular loop of Rs1 with that of human 5-HT(1A) conferred ligand-mediated G(i) signaling. This G(i)-coupled RASSL, Rs1.3, exhibited no measurable signaling to the G(s) or G(q) pathway. These findings show that the signaling repertoire of Rs1 can be expanded and controlled by receptor engineering and drug selection.

Age-related macular degeneration and genetic forms of blindness such as Best Disease and Retinitis Pigmentosa can be caused by degeneration of the Retinal Pigment Epithelium (RPE). RPE generated from patient-derived induced pluripotent stem cells (iPSCs) is valuable for both the study of disease mechanisms and development of therapeutic strategies. However, protocols to produce iPSC-derived RPE in vitro are often inefficient, labor-intensive, low-throughput, and highly variable between cell lines and within batches. Here, we report a robust, scalable method to generate iPSC-RPE using doxycycline-inducible expression of eye field transcription factors OTX2, PAX6 and MITF paired with RPE-permissive culture media. Doxycycline addition induces exogenous expression of these transcription factors in Best Disease patient- and wildtype iPSCs to efficiently produce monolayers of RPE with characteristic morphology and gene expression. Further, these RPE monolayers display functionality features including light absorption via pigmentation, polarity-driven fluid transport, and phagocytosis. With this method, we achieve a highly efficient and easily scalable differentiation without the need for mechanical isolation or enrichment methods, generating RPE cultures applicable for in vitro studies.

G protein-coupled receptor (GPCR) signaling in osteoblasts (OBs) is an important regulator of bone formation. We previously described a mouse model expressing Rs1, an engineered constitutively active Gs-coupled GPCR, under the control of the 2.3 kb Col I promoter. These mice showed a dramatic age-dependent increase in trabecular bone of femurs. Here, we further evaluated the effects of enhanced Gs signaling in OBs on intramembranous bone formation by examining calvariae of 1- and 9-week-old Col1(2.3)/Rs1 mice and characterized the in vivo gene expression specifically occurring in osteoblasts with activated Gs G protein-coupled receptor signaling, at the cellular level rather than in a whole bone. Rs1 calvariae displayed a dramatic increase in bone volume with partial loss of cortical structure. By immunohistochemistry, Osterix was detected in cells throughout the inter-trabecular space while Osteocalcin was expressed predominantly in cells along bone surfaces, suggesting the role of paracrine mediators secreted from OBs driven by 2.3 kb Col I promoter could influence early OB commitment, differentiation, and/or proliferation. Gene expression analysis of calvarial OBs revealed that genes affected by Rs1 signaling include those encoding proteins important for cell differentiation, cytokines and growth factors, angiogenesis, coagulation, and energy metabolism. The set of Gs-GPCRs and other GPCRs that may contribute to the observed skeletal phenotype and candidate paracrine mediators of the effect of Gs signaling in OBs were also determined. Our results identify novel detailed in vivo cellular changes of the anabolic response of the skeleton to Gs signaling in mature OBs.

Long QT syndrome (LQTS) is caused by functional alterations in cardiac ion channels and is associated with prolonged cardiac repolarization time and increased risk of ventricular arrhythmias. Inherited type 2 LQTS (LQT2) and drug-induced LQTS both result from altered function of the hERG channel. We investigated whether the electrophysiological characteristics of LQT2 can be recapitulated in vitro using induced pluripotent stem cell (iPSC) technology. Spontaneously beating cardiomyocytes were differentiated from two iPSC lines derived from an individual with LQT2 carrying the R176W mutation in the KCNH2 (HERG) gene. The individual had been asymptomatic except for occasional palpitations, but his sister and father had died suddenly at an early age. Electrophysiological properties of LQT2-specific cardiomyocytes were studied using microelectrode array and patch-clamp, and were compared with those of cardiomyocytes derived from control cells. The action potential duration of LQT2-specific cardiomyocytes was significantly longer than that of control cardiomyocytes, and the rapid delayed potassium channel (I(Kr)) density of the LQT2 cardiomyocytes was significantly reduced. Additionally, LQT2-derived cardiac cells were more sensitive than controls to potentially arrhythmogenic drugs, including sotalol, and demonstrated arrhythmogenic electrical activity. Consistent with clinical observations, the LQT2 cardiomyocytes demonstrated a more pronounced inverse correlation between the beating rate and repolarization time compared with control cells. Prolonged action potential is present in LQT2-specific cardiomyocytes derived from a mutation carrier and arrhythmias can be triggered by a commonly used drug. Thus, the iPSC-derived, disease-specific cardiomyocytes could serve as an important platform to study pathophysiological mechanisms and drug sensitivity in LQT2.

We introduce GO-Elite, a flexible and powerful pathway analysis tool for a wide array of species, identifiers (IDs), pathways, ontologies and gene sets. In addition to the Gene Ontology (GO), GO-Elite allows the user to perform over-representation analysis on any structured ontology annotations, pathway database or biological IDs (e.g. gene, protein or metabolite). GO-Elite exploits the structured nature of biological ontologies to report a minimal set of non-overlapping terms. The results can be visualized on WikiPathways or as networks. Built-in support is provided for over 60 species and 50 ID systems, covering gene, disease and phenotype ontologies, multiple pathway databases, biomarkers, and transcription factor and microRNA targets. GO-Elite is available as a web interface, GenMAPP-CS plugin and as a cross-platform application.

Lipolysis in adipocytes is regulated by phosphorylation of lipid droplet-associated proteins, including perilipin 1A and hormone-sensitive lipase (HSL). Perilipin 1A is potentially phosphorylated by cAMP(adenosine 3',5'-cyclic monophosphate)-dependent protein kinase (PKA) on several sites, including conserved C-terminal residues, serine 497 (PKA-site 5) and serine 522 (PKA-site 6). To characterize perilipin 1A phosphorylation, novel monoclonal antibodies were developed, which selectively recognize perilipin 1A phosphorylation at PKA-site 5 and PKA-site 6. Utilizing these novel antibodies, as well as antibodies selectively recognizing HSL phosphorylation at serine 563 or serine 660, we used high content analysis to examine the phosphorylation of perilipin 1A and HSL in adipocytes exposed to lipolytic agents. We found that perilipin PKA-site 5 and HSL-serine 660 were phosphorylated to a similar extent in response to forskolin (FSK) and L-γ-melanocyte stimulating hormone (L-γ-MSH). In contrast, perilipin PKA-site 6 and HSL-serine 563 were phosphorylated more slowly and L-γ-MSH was a stronger agonist for these sites compared to FSK. When a panel of lipolytic agents was tested, including multiple concentrations of isoproterenol, FSK, and L-γ-MSH, the pattern of results was virtually identical for perilipin PKA-site 5 and HSL-serine 660, whereas a distinct pattern was observed for perilipin PKA-site 6 and HSL-serine 563. Notably, perilipin PKA-site 5 and HSL-serine 660 feature two arginine residues upstream from the phospho-acceptor site, which confers high affinity for PKA, whereas perilipin PKA-site 6 and HSL-serine 563 feature only a single arginine. Thus, we suggest perilipin 1A and HSL are differentially phosphorylated in a similar manner at the initiation of lipolysis and arginine residues near the target serines may influence this process.

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.

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.

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.

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.