Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks

Multi-omics analysis in primary T cells elucidates mechanisms behind disease-associated genetic loci

BACKGROUND

Genome-wide association studies (GWAS) have uncovered the genetic basis behind many diseases and conditions. However, most of these genetic loci affect regulatory regions, making the interpretation challenging. Chromatin conformation has a fundamental role in gene regulation and is frequently used to associate potential target genes to regulatory regions. However, previous studies mostly used small sample sizes and immortalized cell lines instead of primary cells.

RESULTS

Here we present the most extensive dataset of chromatin conformation with matching gene expression and chromatin accessibility from primary CD4+ and CD8+ T cells to date, isolated from psoriatic arthritis patients and healthy controls. We generated 108 Hi-C libraries (49 billion reads), 128 RNA-seq libraries and 126 ATAC-seq libraries. These data enhance our understanding of the mechanisms by which GWAS variants impact gene regulation, revealing how genetic variation alters chromatin accessibility and structure in primary cells at an unprecedented scale. We refine the mapping of GWAS loci to implicated regulatory elements, such as CTCF binding sites and other enhancer elements, aiding gene assignment. We uncover BCL2L11 as the probable causal gene within the rheumatoid arthritis (RA) locus rs13396472, despite the GWAS variants' intronic positioning relative to ACOXL, and we identify mechanisms involving SESN3 dysregulation in the RA locus rs4409785.

CONCLUSIONS

Given these genes' significant role in T cell development and maturation, our work deepens our comprehension of autoimmune disease pathogenesis, suggesting potential treatment targets. In addition, our dataset provides a valuable resource for the investigation of immune-mediated diseases and gene regulatory mechanisms.

Global chromatin reorganization and regulation of genes with specific evolutionary ages during differentiation and cancer

Cancer cells are highly plastic, favoring adaptation to changing conditions. Genes related to basic cellular processes evolved in ancient species, while more specialized genes appeared later with multicellularity (metazoan genes) or even after mammals evolved. Transcriptomic analyses have shown that ancient genes are up-regulated in cancer, while metazoan-origin genes are inactivated. Despite the importance of these observations, the underlying mechanisms remain unexplored. Here, we study local and global epigenomic mechanisms that may regulate genes from specific evolutionary periods. Using evolutionary gene age data, we characterize the epigenomic landscape, gene expression regulation, and chromatin organization in several cell types: human embryonic stem cells, normal primary B-cells, primary chronic lymphocytic leukemia malignant B-cells, and primary colorectal cancer samples. We identify topological changes in chromatin organization during differentiation observing patterns in Polycomb repression and RNA polymerase II pausing, which are reversed during oncogenesis. Beyond the non-random organization of genes and chromatin features in the 3D epigenome, we suggest that these patterns lead to preferential interactions among ancient, intermediate, and recent genes, mediated by RNA polymerase II, Polycomb, and the lamina, respectively. Our findings shed light on gene regulation according to evolutionary age and suggest this organization changes across differentiation and oncogenesis.

Identification of a distal enhancer of Ucp1 essential for thermogenesis and mitochondrial function in brown fat

Uncoupling protein 1 (UCP1) is a crucial protein located in the mitochondrial inner membrane that mediates nonshivering thermogenesis. However, the molecular mechanisms by which enhancer-promoter chromatin interactions control Ucp1 transcriptional regulation in brown adipose tissue (BAT) are unclear. Here, we employed circularized chromosome conformation capture coupled with next-generation sequencing (4C-seq) to generate high-resolution chromatin interaction profiles of Ucp1 in interscapular brown adipose tissue (iBAT) and epididymal white adipose tissue (eWAT) and revealed marked changes in Ucp1 chromatin interaction between iBAT and eWAT. Next, we identified four iBAT-specific active enhancers of Ucp1, and three of them were activated by cold stimulation. Transcriptional repression of the Ucp1-En4 or Ucp1-En6 region significantly downregulated Ucp1 and impaired mitochondrial function in brown adipocytes. Furthermore, depletion of the cohesin subunit RAD21 decreased the interaction intensity between Ucp1-En4 and the Ucp1 promoter and downregulated Ucp1. EBF2 cooperated with the acetyltransferase CBP to regulate Ucp1-En4 activity and increase Ucp1 transcriptional activity. In vivo, lentivirus-mediated repression of Ucp1-En4 was injected into iBAT, resulting in impacted iBAT thermogenic capacity and impaired iBAT mitochondrial function under cold acclimation conditions. Studying the functional enhancers regulating Ucp1 expression in iBAT will provide important insights into the regulatory mechanisms of BAT activity.

LncRNA DANCR suppresses acute myocardial infarction in mice via mediating p-RXRA/TRAF2/NIK/IKK/NF-κB signaling pathway

OBJECTIVES

This study aimed to investigate the role of LncRNA differentiation antagonizing non-protein coding RNA (DANCR) in acute myocardial infarction (AMI).

METHODS

A mouse model of AMI was established, and the cardiac contractile function was detected. Next, cardiomyocytes treated with oxygen-glucose deprivation (OGD) were used for gain- and loss-function experiments in vitro. RIP analysis was used to verify the binding of DANCR and Retinoid X receptor alpha (RXRA), and Co-IP assay was used to measure the binding of phosphorylated RXRA to TNF receptor associated factor 2 (TRAF2).

RESULTS

The expression of DANCR in myocardial tissues of AMI mice were significantly downregulated. Overexpression of DANCR decreased myocardial infarct size and enhanced cardiac contractile function in AMI mice. Moreover, overexpression of DANCR promoted proliferation and inhibited apoptosis in OGD-induced cardiomyocytes. Mechanism studies demonstrated that DANCR interacted with RXRA and promoted glycogen synthase kinase 3beta (GSK3β)-mediated phosphorylation of RXRA, and phosphorylated RXRA interacted with TRAF2 protein to downregulate TRAF2 protein level. Bexarotene (Bex), an activator of RXRA, inhibited TRAF2 protein expression, while RXRA overexpression had no effect on TRAF2 protein expression. Bex treatment or silencing TRAF2 promoted proliferation and inhibited apoptosis in cardiomyocytes. In addition, silencing DANCR inhibited cardiomyocyte proliferation and induced apoptosis by activating the NIK/IKK/NF-κB pathway, while B022, an inhibitor of NIK, counteracted the effects of DANCR silencing on cardiomyocytes.

CONCLUSIONS

Studies demonstrated that DANCR suppressed AMI in mice by mediating p-RXRA/TRAF2/NIK/IKK/NF-κB pathway.

Global chromatin reorganization and regulation of genes with specific evolutionary ages during differentiation and cancer

Cancer cells are highly plastic, allowing them to adapt to changing conditions. Genes related to basic cellular processes evolved in ancient species, while more specialized genes appeared later with multicellularity (metazoan genes) or even after mammals evolved. Transcriptomic analyses have shown that ancient genes are up-regulated in cancer, while metazoan-origin genes are inactivated. Despite the importance of these observations, the underlying mechanisms remain unexplored. Here, we study local and global epigenomic mechanisms that may regulate genes from specific evolutionary periods. Using evolutionary gene age data, we characterize the epigenomic landscape, gene expression regulation, and chromatin organization in three cell types: human embryonic stem cells, normal B-cells, and primary cells from Chronic Lymphocytic Leukemia, a B-cell malignancy. We identify topological changes in chromatin organization during differentiation observing patterns in Polycomb repression and RNA Polymerase II pausing, which are reversed during oncogenesis. Beyond the non-random organization of genes and chromatin features in the 3D epigenome, we suggest that these patterns lead to preferential interactions among ancient, intermediate, and recent genes, mediated by RNA Polymerase II, Polycomb, and the lamina, respectively. Our findings shed light on gene regulation according to evolutionary age and suggest this organization changes across differentiation and oncogenesis.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

Priority index for critical Covid-19 identifies clinically actionable targets and drugs

While genome-wide studies have identified genomic loci in hosts associated with life-threatening Covid-19 (critical Covid-19), the challenge of resolving these loci hinders further identification of clinically actionable targets and drugs. Building upon our previous success, we here present a priority index solution designed to address this challenge, generating the target and drug resource that consists of two indexes: the target index and the drug index. The primary purpose of the target index is to identify clinically actionable targets by prioritising genes associated with Covid-19. We illustrate the validity of the target index by demonstrating its ability to identify pre-existing Covid-19 phase-III drug targets, with the majority of these targets being found at the leading prioritisation (leading targets). These leading targets have their evolutionary origins in Amniota ('four-leg vertebrates') and are predominantly involved in cytokine-cytokine receptor interactions and JAK-STAT signaling. The drug index highlights opportunities for repurposing clinically approved JAK-STAT inhibitors, either individually or in combination. This proposed strategic focus on the JAK-STAT pathway is supported by the active pursuit of therapeutic agents targeting this pathway in ongoing phase-II/III clinical trials for Covid-19.

In-Depth Genomic Analysis: The New Challenge in Congenital Heart Disease

The use of next-generation sequencing has provided new insights into the causes and mechanisms of congenital heart disease (CHD). Examinations of the whole exome sequence have detected detrimental gene variations modifying single or contiguous nucleotides, which are characterised as pathogenic based on statistical assessments of families and correlations with congenital heart disease, elevated expression during heart development, and reductions in harmful protein-coding mutations in the general population. Patients with CHD and extracardiac abnormalities are enriched for gene classes meeting these criteria, supporting a common set of pathways in the organogenesis of CHDs. Single-cell transcriptomics data have revealed the expression of genes associated with CHD in specific cell types, and emerging evidence suggests that genetic mutations disrupt multicellular genes essential for cardiogenesis. Metrics and units are being tracked in whole-genome sequencing studies.

Neuroplastin Expression in Male Mice Is Essential for Fertility, Mating, and Adult Testosterone Levels

Male reproduction depends on hormonally driven behaviors and numerous genes for testis development and spermatogenesis. Neuroplastin-deficient (Nptn-/-) male mice cannot sire offspring. By immunohistochemistry, we characterized neuroplastin expression in the testis. Breeding, mating behavior, hormonal regulation, testicular development, and spermatogenesis were analyzed in cell-type specific neuroplastin mutant mice. Leydig, Sertoli, peritubular myoid, and germ cells express Np, but spermatogenesis and sperm number are not affected in Nptn-/- males. Neuroplastin lack from CNS neurons or restricted to spermatogonia or Sertoli cells permitted reproduction. Normal luteinizing hormone (LH) and follicle-stimulating hormone (FSH) blood levels in Nptn-/- males support undisturbed hormonal regulation in the brain. However, Nptn-/- males lack mounting behavior accompanied by low testosterone blood levels. Testosterone rise from juvenile to adult blood levels is absent in Nptn-/- males. LH-receptor stimulation raising intracellular Ca2+ in Leydig cells triggers testosterone production. Reduced Plasma Membrane Ca2+ ATPase 1 (PMCA1) in Nptn-/- Leydig cells suggests that Nptn-/- Leydig cells produce sufficient testosterone for testis and sperm development, but a lack of PMCA-Np complexes prevents the increase from reaching adult blood levels. Behavioral immaturity with low testosterone blood levels underlies infertility of Nptn-/- males, revealing that Np is essential for reproduction.

Priority index for asthma (PIA): In silico discovery of shared and distinct drug targets for adult- and childhood-onset disease

Asthma is a chronic disease that is caused by a combination of genetic risks and environmental triggers and can affect both adults and children. Genome-wide association studies have revealed partly distinct genetic architectures for its two age-of-onset subtypes (namely, adult-onset and childhood-onset). We reason that identifying shared and distinct drug targets between these subtypes may inform the development of subtype-specific therapeutic strategies. In attempting this, we here introduce Priority Index for Asthma or PIA, a genetics-led and network-driven drug target prioritisation tool for asthma. We demonstrate the validity of the tool in improving drug target prioritisation for asthma compared to the status quo methods, as well as in capturing the underlying etiology and existing therapeutics for the disease. We also illustrate how PIA can be used to prioritise drug targets for adult- and childhood-onset asthma, as well as to identify shared and distinct pathway crosstalk genes. Shared crosstalk genes are mostly involved in JAK-STAT signaling, with clinical evidence supporting that targeting this pathway may be a promising drug repurposing opportunity for both subtypes. Crosstalk genes specific to childhood-onset asthma are enriched for PI3K-AKT-mTOR signaling, and we identify genes that are already targeted by licensed medications as repurposed drug candidates for this subtype. We make all our results accessible and reproducible at http://www.genetictargets.com/PIA. Collectively, our study has significant implications for asthma computational medicine research and can guide the future development of subtype-specific therapeutic strategies for the disease.

Genomic Evidence Supports the Recognition of Endometriosis as an Inflammatory Systemic Disease and Reveals Disease-Specific Therapeutic Potentials of Targeting Neutrophil Degranulation

Background

Endometriosis, classically viewed as a localized disease, is increasingly recognized as a systemic disease with multi-organ effects. This disease is highlighted by systemic inflammation in affected organs and by high comorbidity with immune-mediated diseases.

Results

We provide genomic evidence to support the recognition of endometriosis as an inflammatory systemic disease. This was achieved through our genomics-led target prioritization, called 'END', that leverages the value of multi-layered genomic datasets (including genome-wide associations in disease, regulatory genomics, and protein interactome). Our prioritization recovered existing proof-of-concept therapeutic targeting in endometriosis and outperformed competing prioritization approaches (Open Targets and Naïve prioritization). Target genes at the leading prioritization revealed molecular hallmarks (and possibly the cellular basis as well) that are consistent with systemic disease manifestations. Pathway crosstalk-based attack analysis identified the critical gene AKT1. In the context of this gene, we further identified genes that are already targeted by licensed medications in other diseases, such as ESR1. Such analysis was supported by current interests targeting the PI3K/AKT/mTOR pathway in endometriosis and by the fact that therapeutic agents targeting ESR1 are now under active clinical trials in disease. The construction of cross-disease prioritization map enabled the identification of shared and distinct targets between endometriosis and immune-mediated diseases. Shared target genes identified opportunities for repurposing existing immunomodulators, particularly disease-modifying anti-rheumatic drugs (such as TNF, IL6 and IL6R blockades, and JAK inhibitors). Genes highly prioritized only in endometriosis revealed disease-specific therapeutic potentials of targeting neutrophil degranulation - the exocytosis that can facilitate metastasis-like spread to distant organs causing inflammatory-like microenvironments.

Conclusion

Improved target prioritization, along with an atlas of in silico predicted targets and repurposed drugs (available at https://23verse.github.io/end), provides genomic insights into endometriosis, reveals disease-specific therapeutic potentials, and expands the existing theories on the origin of disease.

Epigenome rewiring in human pluripotent stem cells

The epigenome plays a crucial role in modulating the activity of regulatory elements, thereby orchestrating diverse transcriptional programs during embryonic development. Human (h)PSC stepwise differentiation provides an excellent platform for capturing dynamic epigenomic events during lineage transition in human development. Here we discuss how recent technological advances, from epigenomic mapping to targeted perturbation, are providing a more comprehensive appreciation of remodeling of the chromatin landscape during human development with implications for aberrant rewiring in disease. We predict that the continuous innovation of hPSC differentiation methods, epigenome mapping, and CRISPR (clustered regularly interspaced short palindromic repeats) perturbation technologies will allow researchers to build toward not only a comprehensive understanding of the epigenomic mechanisms governing development, but also a highly flexible way to model diseases with opportunities for translation.

Retrograde Amnesia - A Question of Disturbed Calcium Levels?

Retrograde amnesia is the inability to remember events or information. The successful acquisition and memory of information is required before retrograde amnesia may occur. Often, the trigger for retrograde amnesia is a traumatic event. Loss of memories may be caused in two ways: either by loss/erasure of the memory itself or by the inability to access the memory, which is still present. In general, memories and learning are associated with a positive connotation although the extinction of unpleasant experiences and memories of traumatic events may be highly welcome. In contrast to the many experimental models addressing learning deficits caused by anterograde amnesia, the incapability to acquire new information, retrograde amnesia could so far only be investigated sporadically in human patients and in a limited number of model systems. Apart from models and diseases in which neurodegeneration or dementia like Alzheimer’s disease result in loss of memory, retrograde amnesia can be elicited by various drugs of which alcohol is the most prominent one and exemplifies the non-specific effects and the variable duration. External or internal impacts like traumatic brain injury, stroke, or electroconvulsive treatments may similarly result in variable degrees of retrograde amnesia. In this review, I will discuss a new genetic approach to induce retrograde amnesia in a mouse model and raise the hypothesis that retrograde amnesia is caused by altered intracellular calcium homeostasis. Recently, we observed that neuronal loss of neuroplastin resulted in retrograde amnesia specifically for associative memories. Neuroplastin is tightly linked to the expression of the main Ca²⁺ extruding pumps, the plasma membrane calcium ATPases (PMCAs). Therefore, neuronal loss of neuroplastin may block the retrieval and storage of associative memories by interference with Ca²⁺ signaling cascades. The possibility to elicit retrograde amnesia in a controlled manner allows to investigate the underlying mechanisms and may provide a deeper understanding of the molecular and circuit processes of memory.

Genomic enhancers in cardiac development and disease

The Human Genome Project marked a major milestone in the scientific community as it unravelled the ~3 billion bases that are central to crucial aspects of human life. Despite this achievement, it only scratched the surface of understanding how each nucleotide matters, both individually and as part of a larger unit. Beyond the coding genome, which comprises only ~2% of the whole genome, scientists have realized that large portions of the genome, not known to code for any protein, were crucial for regulating the coding genes. These large portions of the genome comprise the 'non-coding genome'. The history of gene regulation mediated by proteins that bind to the regulatory non-coding genome dates back many decades to the 1960s. However, the original definition of 'enhancers' was first used in the early 1980s. In this Review, we summarize benchmark studies that have mapped the role of cardiac enhancers in disease and development. We highlight instances in which enhancer-localized genetic variants explain the missing link to cardiac pathogenesis. Finally, we inspire readers to consider the next phase of exploring enhancer-based gene therapy for cardiovascular disease.

A novel RNA-mediated mechanism causing down-regulation of insulating promoter interactions in human embryonic stem cells

The genome-wide promoter interactome is primarily maintained and regulated by architectural proteins such as CTCF and cohesin. However, some studies suggest a role for non-coding RNAs (ncRNAs) in this process. We aimed to characterise the regulatory role of RNA-mediated promoter interactions in the control of gene expression. We integrated genome-wide datasets of RNA-chromatin and promoter-genome interactions in human embryonic stem cells (hESCs) to identify putative RNA-mediated promoter interactions. We discovered that CTCF sites were enriched in RNA-PIRs (promoter interacting regions co-localising with RNA-chromatin interaction sites) and genes interacting with RNA-PIRs containing CTCF sites showed higher expression levels. One of the long noncoding RNAs (lncRNAs) expressed in hESCs, Syntaxin 18-Antisense 1 (STX18-AS1), appeared to be involved in an insulating promoter interaction with the neighbouring gene, MSX1. By knocking down STX18-AS1, the MSX1 promoter-PIR interaction was intensified and the target gene (MSX1) expression was down-regulated. Conversely, reduced MSX1 promoter-PIR interactions, resulting from CRISPR-Cas9 deletion of the PIR, increased the expression of MSX1. We conclude that STX18-AS1 RNA antagonised local CTCF-mediated insulating promoter interactions to augment gene expression. Such down-regulation of the insulating promoter interactions by this novel mechanism may explain the higher expression of genes interacting with RNA-PIRs linked to CTCF sites.

Restriction enzyme selection dictates detection range sensitivity in chromatin conformation capture-based variant-to-gene mapping approaches

Promoter-focused chromatin conformation techniques directly detect interactions between gene promoters and distal genomic sequences, providing structural information relevant to gene regulation without the excessive non-genic architectural data generated by full-scale Hi-C. 3D promoter 'interactome' maps are crucial for understanding how epigenomic features such as histone modifications and open chromatin, or genetic variants identified in genome-wide association studies (GWAS), contribute to biological function. However, variation in sensitivity between such promoter-focused methods, principally due to restriction enzyme selection, has not been systematically assessed. Here, we performed a head-to-head comparison of promoter capture datasets using 4 cutters (DpnII or MboI) versus the 6 cutter HindIII from the same five cell types. While HindIII generally produces a higher signal-to-noise ratio for significant interactions in comparison to 4-cutters, we show that DpnII/MboI detects more proximal interactions and shows little overlap with the HindIII detection range. Promoter-interacting genomic regions mapped by 4-cutters are more enriched for regulatory features and disease-associated genetic variation than 6-cutters maps, suggesting that high-resolution maps better capture gene regulatory architectures than do lower resolution approaches.

Environmental Alterations during Embryonic Development: Studying the Impact of Stressors on Pluripotent Stem Cell-Derived Cardiomyocytes

Non-communicable diseases (NCDs) sauch as diabetes, obesity and cardiovascular diseases are rising rapidly in all countries world-wide. Environmental maternal factors (e.g., diet, oxidative stress, drugs and many others), maternal illnesses and other stressors can predispose the newborn to develop diseases during different stages of life. The connection between environmental factors and NCDs was formulated by David Barker and colleagues as the Developmental Origins of Health and Disease (DOHaD) hypothesis. In this review, we describe the DOHaD concept and the effects of several environmental stressors on the health of the progeny, providing both animal and human evidence. We focus on cardiovascular diseases which represent the leading cause of death worldwide. The purpose of this review is to discuss how in vitro studies with pluripotent stem cells (PSCs), such as embryonic and induced pluripotent stem cells (ESC, iPSC), can underpin the research on non-genetic heart conditions. The PSCs could provide a tool to recapitulate aspects of embryonic development “in a dish”, studying the effects of environmental exposure during cardiomyocyte (CM) differentiation and maturation, establishing a link to molecular mechanism and epigenetics.

Leveraging three-dimensional chromatin architecture for effective reconstruction of enhancer-target gene regulatory interactions

A growing amount of evidence in literature suggests that germline sequence variants and somatic mutations in non-coding distal regulatory elements may be crucial for defining disease risk and prognostic stratification of patients, in genetic disorders as well as in cancer. Their functional interpretation is challenging because genome-wide enhancer-target gene (ETG) pairing is an open problem in genomics. The solutions proposed so far do not account for the hierarchy of structural domains which define chromatin three-dimensional (3D) architecture. Here we introduce a change of perspective based on the definition of multi-scale structural chromatin domains, integrated in a statistical framework to define ETG pairs. In this work (i) we develop a computational and statistical framework to reconstruct a comprehensive map of ETG pairs leveraging functional genomics data; (ii) we demonstrate that the incorporation of chromatin 3D architecture information improves ETG pairing accuracy and (iii) we use multiple experimental datasets to extensively benchmark our method against previous solutions for the genome-wide reconstruction of ETG pairs. This solution will facilitate the annotation and interpretation of sequence variants in distal non-coding regulatory elements. We expect this to be especially helpful in clinically oriented applications of whole genome sequencing in cancer and undiagnosed genetic diseases research.

Neuroplastin in Neuropsychiatric Diseases

Molecular mechanisms underlying neuropsychiatric and neurodegenerative diseases are insufficiently elucidated. A detailed understanding of these mechanisms may help to further improve medical intervention. Recently, intellectual abilities, creativity, and amnesia have been associated with neuroplastin, a cell recognition glycoprotein of the immunoglobulin superfamily that participates in synapse formation and function and calcium signaling. Data from animal models suggest a role for neuroplastin in pathways affected in neuropsychiatric and neurodegenerative diseases. Neuroplastin loss or disruption of molecular pathways related to neuronal processes has been linked to various neurological diseases, including dementia, schizophrenia, and Alzheimer's disease. Here, we review the molecular features of the cell recognition molecule neuroplastin, and its binding partners, which are related to neurological processes and involved in learning and memory. The emerging functions of neuroplastin may have implications for the treatment of diseases, particularly those of the nervous system.

Detecting chromosomal interactions in Capture Hi-C data with CHiCAGO and companion tools

Capture Hi-C is widely used to obtain high-resolution profiles of chromosomal interactions involving, at least on one end, regions of interest such as gene promoters. Signal detection in Capture Hi-C data is challenging and cannot be adequately accomplished with tools developed for other chromosome conformation capture methods, including standard Hi-C. Capture Hi-C Analysis of Genomic Organization (CHiCAGO) is a computational pipeline developed specifically for Capture Hi-C analysis. It implements a statistical model accounting for biological and technical background components, as well as bespoke normalization and multiple testing procedures for this data type. Here we provide a step-by-step guide to the CHiCAGO workflow that is aimed at users with basic experience of the command line and R. We also describe more advanced strategies for tuning the key parameters for custom experiments and provide guidance on data preprocessing and downstream analysis using companion tools. In a typical experiment, CHiCAGO takes ~2-3 h to run, although pre- and postprocessing steps may take much longer.

Chromatin Looping Links Target Genes with Genetic Risk Loci for Dermatological Traits

Chromatin looping between regulatory elements and gene promoters presents a potential mechanism whereby disease risk variants affect their target genes. In this study, we use H3K27ac HiChIP, a method for assaying the active chromatin interactome in two cell lines: keratinocytes and skin lymphoma-derived CD8+ T cells. We integrate public datasets for a lymphoblastoid cell line and primary CD4+ T cells and identify gene targets at risk loci for skin-related disorders. Interacting genes enrich for pathways of known importance in each trait, such as cytokine response (psoriatic arthritis and psoriasis) and replicative senescence (melanoma). We show examples of how our analysis can inform changes in the current understanding of multiple psoriasis-associated risk loci. For example, the variant rs10794648, which is generally assigned to IFNLR1, was linked to GRHL3, a gene essential in skin repair and development, in our dataset. Our findings, therefore, indicate a renewed importance of skin-related factors in the risk of disease.

Long-range linkage disequilibrium in French beef cattle breeds

BACKGROUND

Linkage disequilibrium (LD) is a key parameter to study the history of populations and to identify and fine map quantitative trait loci (QTL) and it has been studied for many years in animal populations. The advent of new genotyping technologies has allowed whole-genome LD studies in most cattle populations. However, to date, long-range LD (LRLD) between distant variants on the genome has not been investigated in detail in cattle. Here, we present the first comprehensive study of LRLD in French beef cattle by analysing data on 672 Charolais (CHA), 462 Limousine (LIM) and 326 Blonde d'Aquitaine (BLA) individuals that were genotyped on the Illumina BovineHD Beadchip. Furthermore, whole-genome LD and haplotype block structure were analysed in these three breeds.

RESULTS

We computed linkage disequilibrium (r²) values for 5.9, 5.6 and 6.0 billion pairs of SNPs on the 29 autosomes of CHA, LIM and BLA, respectively. Mean r² values drop to less than 0.1 for distances between SNPs greater than 120 kb. However, for the first time, we detected the existence of LRLD in the three main French beef breeds. In total, 598, 266, and 795 LRLD events (r² ≥ 0.6) were detected in CHA, LIM and BLA, respectively. Each breed had predominantly population-specific LRLD interactions, although shared LRLD events occurred in a number of regions (55 LRLD events were shared between two breeds and nine between the three breeds). Examples of possible functional gene interactions and QTL co-location were observed with some of these LRLD events, which suggests epistatic selection.

CONCLUSIONS

We identified long-range linkage disequilibrium for the first time in French beef cattle populations. Epistatic selection may be the main source of the observed LRLD events, but other forces may also be involved. LRLD information should be accounted for in genome-wide association studies.

Neuroplastin expression is essential for hearing and hair cell PMCA expression

Hearing deficits impact on the communication with the external world and severely compromise perception of the surrounding. Deafness can be caused by particular mutations in the neuroplastin (Nptn) gene, which encodes a transmembrane recognition molecule of the immunoglobulin (Ig) superfamily and plasma membrane Calcium ATPase (PMCA) accessory subunit. This study investigates whether the complete absence of neuroplastin or the loss of neuroplastin in the adult after normal development lead to hearing impairment in mice analyzed by behavioral, electrophysiological, and in vivo imaging measurements. Auditory brainstem recordings from adult neuroplastin-deficient mice (Nptn^-/-) show that these mice are deaf. With age, hair cells and spiral ganglion cells degenerate in Nptn^-/- mice. Adult Nptn^-/- mice fail to behaviorally respond to white noise and show reduced baseline blood flow in the auditory cortex (AC) as revealed by single-photon emission computed tomography (SPECT). In adult Nptn^-/- mice, tone-evoked cortical activity was not detectable within the primary auditory field (A1) of the AC, although we observed non-persistent tone-like evoked activities in electrophysiological recordings of some young Nptn^-/- mice. Conditional ablation of neuroplastin in Nptn^{lox/loxEmx1Cre} mice reveals that behavioral responses to simple tones or white noise do not require neuroplastin expression by central glutamatergic neurons. Loss of neuroplastin from hair cells in adult Nptn^{Δlox/loxPrCreERT} mice after normal development is correlated with increased hearing thresholds and only high prepulse intensities result in effective prepulse inhibition (PPI) of the startle response. Furthermore, we show that neuroplastin is required for the expression of PMCA 2 in outer hair cells. This suggests that altered Ca²⁺ homeostasis underlies the observed hearing impairments and leads to hair cell degeneration. Our results underline the importance of neuroplastin for the development and the maintenance of the auditory system.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

Using functional genomics to advance the understanding of psoriatic arthritis

Psoriatic arthritis (PsA) is a complex disease where susceptibility is determined by genetic and environmental risk factors. Clinically, PsA involves inflammation of the joints and the skin, and, if left untreated, results in irreversible joint damage. There is currently no cure and the few treatments available to alleviate symptoms do not work in all patients. Over the past decade, genome-wide association studies (GWAS) have uncovered a large number of disease-associated loci but translating these findings into functional mechanisms and novel targets for therapeutic use is not straightforward. Most variants have been predicted to affect primarily long-range regulatory regions such as enhancers. There is now compelling evidence to support the use of chromatin conformation analysis methods to discover novel genes that can be affected by disease-associated variants. Here, we will review the studies published in the field that have given us a novel understanding of gene regulation in the context of functional genomics and how this relates to the study of PsA and its underlying disease mechanism.

Three-dimensional chromatin organization in cardiac development and disease

Recent technological advancements in the field of chromatin biology have rewritten the textbook on nuclear organization. We now appreciate that the folding of chromatin in the three-dimensional space (i.e. its 3D "architecture") is non-random, hierarchical, and highly complex. While 3D chromatin structure is partially encoded in the primary sequence and thereby broadly conserved across cell types and states, a substantial portion of the genome seems to be dynamic during development or in disease. Moreover, there is growing evidence that at least some of the 3D structure of chromatin is functionally linked to gene regulation, both being modulated by and impacting on multiple nuclear processes (including DNA replication, transcription, and RNA splicing). In recent years, these new concepts have nourished several investigations about the functional role of 3D chromatin topology dynamics in the heart during development and disease. This review aims to provide a comprehensive overview of our current understanding in this field, and to discuss how this knowledge can inform further research as well as clinical practice.

A modified protocol of Capture-C allows affordable and flexible high-resolution promoter interactome analysis

Large-scale epigenomic projects have mapped hundreds of thousands of potential regulatory sites in the human genome, but only a small proportion of these elements are proximal to transcription start sites. It is believed that the majority of these sequences are remote promoter-activating genomic sites scattered within several hundreds of kilobases from their cognate promoters and referred to as enhancers. It is still unclear what principles, aside from relative closeness in the linear genome, determine which promoter(s) is controlled by a given enhancer; however, this understanding is of great fundamental and clinical relevance. In recent years, C-methods (chromosome conformation capture-based methods) have become a powerful tool for the identification of enhancer-promoter spatial contacts that, in most cases, reflect their functional link. Here, we describe a new hybridisation-based promoter Capture-C protocol that makes use of biotinylated dsDNA probes generated by PCR from a custom pool of long oligonucleotides. The described protocol allows high-resolution promoter interactome description, providing a flexible and cost-effective alternative to the existing promoter Capture-C modifications. Based on the obtained data, we propose several tips on probe design that could potentially improve the results of future experiments.

Epigenomes of Human Hearts Reveal New Genetic Variants Relevant for Cardiac Disease and Phenotype

RATIONALE

Identifying genetic markers for heterogeneous complex diseases such as heart failure is challenging and requires prohibitively large cohort sizes in genome-wide association studies to meet the stringent threshold of genome-wide statistical significance. On the other hand, chromatin quantitative trait loci, elucidated by direct epigenetic profiling of specific human tissues, may contribute toward prioritizing subthreshold variants for disease association.

OBJECTIVE

Here, we captured noncoding genetic variants by performing epigenetic profiling for enhancer H3K27ac chromatin immunoprecipitation followed by sequencing in 70 human control and end-stage failing hearts.

METHODS AND RESULTS

We have mapped a comprehensive catalog of 47 321 putative human heart enhancers and promoters. Three thousand eight hundred ninety-seven differential acetylation peaks (FDR [false discovery rate], 5%) pointed to pathways altered in heart failure. To identify cardiac histone acetylation quantitative trait loci (haQTLs), we regressed out confounding factors including heart failure disease status and used the G-SCI (Genotype-independent Signal Correlation and Imbalance) test1 to call out 1680 haQTLs (FDR, 10%). RNA sequencing performed on the same heart samples proved a subset of haQTLs to have significant association also to gene expression (expression quantitative trait loci), either in cis (180) or through long-range interactions (81), identified by Hi-C (high-throughput chromatin conformation assay) and HiChIP (high-throughput protein centric chromatin) performed on a subset of hearts. Furthermore, a concordant relationship between the gain or disruption of TF (transcription factor)-binding motifs, inferred from alternative alleles at the haQTLs, implied a surprising direct association between these specific TF and local histone acetylation in human hearts. Finally, 62 unique loci were identified by colocalization of haQTLs with the subthreshold loci of heart-related genome-wide association studies datasets.

CONCLUSIONS

Disease and phenotype association for 62 unique loci are now implicated. These loci may indeed mediate their effect through modification of enhancer H3K27 acetylation enrichment and their corresponding gene expression differences (bioRxiv: https://doi.org/10.1101/536763). Graphical Abstract: A graphical abstract is available for this article.

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation

Genome-wide association studies have uncovered over a 100 genetic loci associated with atrial fibrillation (AF), the most common arrhythmia. Many of the top AF-associated loci harbor key cardiac transcription factors, including PITX2, TBX5, PRRX1, and ZFHX3. Moreover, the vast majority of the AF-associated variants lie within noncoding regions of the genome where causal variants affect gene expression by altering the activity of transcription factors and the epigenetic state of chromatin. In this review, we discuss a transcriptional regulatory network model for AF defined by effector genes in Genome-wide association studies loci. We describe the current state of the field regarding the identification and function of AF-relevant gene regulatory networks, including variant regulatory elements, dose-sensitive transcription factor functionality, target genes, and epigenetic states. We illustrate how altered transcriptional networks may impact cardiomyocyte function and ionic currents that impact AF risk. Last, we identify the need for improved tools to identify and functionally test transcriptional components to define the links between genetic variation, epigenetic gene regulation, and atrial function.

Systems Genetics for Mechanistic Discovery in Heart Diseases

Cardiovascular diseases are the leading cause of death worldwide. Complex diseases with highly heterogenous disease progression among patient populations, cardiovascular diseases feature multifactorial contributions from both genetic and environmental stressors. Despite significant effort utilizing multiple approaches from molecular biology to genome-wide association studies, the genetic landscape of cardiovascular diseases, particularly for the nonfamilial forms of heart failure, is still poorly understood. In the past decade, systems-level approaches based on omics technologies have become an important approach for the study of complex traits in large populations. These advances create opportunities to integrate genetic variation with other biological layers to identify and prioritize candidate genes, understand pathogenic pathways, and elucidate gene-gene and gene-environment interactions. In this review, we will highlight some of the recent progress made using systems genetics approaches to uncover novel mechanisms and molecular bases of cardiovascular pathophysiological manifestations. The key technology and data analysis platforms necessary to implement systems genetics will be described, and the current major challenges and future directions will also be discussed. For complex cardiovascular diseases, such as heart failure, systems genetics represents a powerful strategy to obtain mechanistic insights and to develop individualized diagnostic and therapeutic regiments, paving the way for precision cardiovascular medicine.

Efficient differentiation of human pluripotent stem cells into cardiomyocytes on cell sorting thermoresponsive surface

The current differentiation process of human pluripotent stem cells (hPSCs) into cardiomyocytes to enhance the purity of hPSC-derived cardiomyocytes requires some purification processes, which are laborious processes. We developed cell sorting plates, which are prepared from coating thermoresponsive poly(N-isopropylacrylamide) and extracellular matrix proteins. After hPSCs were induced into cardiomyocytes on the thermoresponsive surface coated with laminin-521 for 15 days, the temperature of the cell culture plates was decreased to 8-9 °C to detach the cells partially from the thermoresponsive surface. The detached cells exhibited a higher cardiomyocyte marker of cTnT than the remaining cells on the thermoresponsive surface as well as the cardiomyocytes after purification using conventional cell selection. The detached cells expressed several cardiomyocyte markers, such as α-actinin, MLC2a and NKX2.5. This study suggested that the purification of hPSC-derived cardiomyocytes using cell sorting plates with the thermoresponsive surface is a promising method for the purification of hPSC-derived cardiomyocytes without conventional laborious processes.

Functional effects of variation in transcription factor binding highlight long-range gene regulation by epromoters

Identifying DNA cis-regulatory modules (CRMs) that control the expression of specific genes is crucial for deciphering the logic of transcriptional control. Natural genetic variation can point to the possible gene regulatory function of specific sequences through their allelic associations with gene expression. However, comprehensive identification of causal regulatory sequences in brute-force association testing without incorporating prior knowledge is challenging due to limited statistical power and effects of linkage disequilibrium. Sequence variants affecting transcription factor (TF) binding at CRMs have a strong potential to influence gene regulatory function, which provides a motivation for prioritizing such variants in association testing. Here, we generate an atlas of CRMs showing predicted allelic variation in TF binding affinity in human lymphoblastoid cell lines and test their association with the expression of their putative target genes inferred from Promoter Capture Hi-C and immediate linear proximity. We reveal >1300 CRM TF-binding variants associated with target gene expression, the majority of them undetected with standard association testing. A large proportion of CRMs showing associations with the expression of genes they contact in 3D localize to the promoter regions of other genes, supporting the notion of 'epromoters': dual-action CRMs with promoter and distal enhancer activity.

Enhancer-gene maps in the human and zebrafish genomes using evolutionary linkage conservation

The spatiotemporal expression of genes is controlled by enhancer sequences that bind transcription factors. Identifying the target genes of enhancers remains difficult because enhancers regulate gene expression over long genomic distances. To address this, we used an evolutionary approach to build two genome-wide maps of predicted enhancer-gene associations in the human and zebrafish genomes. Evolutionary conserved sequences were linked to their predicted target genes using PEGASUS, a bioinformatics method that relies on evolutionary conservation of synteny. The analysis of these maps revealed that the number of predicted enhancers linked to a gene correlate with its expression breadth. Comparison of both maps identified hundreds of putative vertebrate ancestral regulatory relationships from which we could determine that predicted enhancer-gene distances scale with genome size despite strong positional conservation. The two maps represent a resource for further studies, including the prioritization of sequence variants in whole genome sequence of patients affected by genetic diseases.

Genome-wide association study on coronary artery disease in type 1 diabetes suggests beta-defensin 127 as a risk locus

AIMS

Diabetes is a known risk factor for coronary artery disease (CAD). There is accumulating evidence that CAD pathogenesis differs for individuals with type 1 diabetes (T1D). However, the genetic background has not been extensively studied. We aimed to discover genetic loci increasing CAD susceptibility, especially in T1D, to examine the function of these discoveries and to study the role of the known risk loci in T1D.

METHODS AND RESULTS

We performed the largest genome-wide association study to date for CAD in T1D, comprising 4869 individuals with T1D (cases/controls: 941/3928). Two loci reached genome-wide significance, rs1970112 in CDKN2B-AS1 [odds ratio (OR) = 1.32, P = 1.50 × 10-8], and rs6055069 on DEFB127 promoter (OR = 4.17, P = 2.35 × 10-9), with consistent results in survival analysis. The CDKN2B-AS1 variant replicated (P = 0.04) when adjusted for diabetic kidney disease in three additional T1D cohorts (cases/controls: 434/3123). Furthermore, we explored the function of the lead discoveries with a cardio-phenome-wide analysis. Among the eight suggestive loci (P < 1 × 10-6), rs70962766 near B3GNT2 associated with central blood pressure, rs1344228 near CNTNAP5 with intima media thickness, and rs2112481 on GRAMD2B promoter with serum leucocyte concentration. Finally, we calculated genetic risk scores for individuals with T1D with the known susceptibility loci. General population risk variants were modestly but significantly associated with CAD also in T1D (P = 4.21 × 10-7).

CONCLUSION

While general population CAD risk loci had limited effect on the risk in T1D, for the first time, variants at the CDKN2B-AS1 locus were robustly associated with CAD in individuals with T1D. The novel finding on β-defensin DEFB127 promoter provides a link between diabetes, infection susceptibility, and CAD, although pending on future confirmation.

Efficient differentiation of human ES and iPS cells into cardiomyocytes on biomaterials under xeno-free conditions

Current xeno-free and chemically defined methods for the differentiation of hPSCs (human pluripotent stem cells) into cardiomyocytes are not efficient and are sometimes not reproducible. Therefore, it is necessary to develop reliable and efficient methods for the differentiation of hPSCs into cardiomyocytes for future use in cardiovascular research related to drug discovery, cardiotoxicity screening, and disease modeling. We evaluated two representative differentiation methods that were reported previously, and we further developed original, more efficient methods for the differentiation of hPSCs into cardiomyocytes under xeno-free, chemically defined conditions. The developed protocol successively differentiated hPSCs into cardiomyocytes, approximately 90-97% of which expressed the cardiac marker cTnT, with beating speeds and sarcomere lengths that were similar to those of a healthy adult human heart. The optimal cell culture biomaterials for the cardiac differentiation of hPSCs were also evaluated using extracellular matrix-mimetic material-coated dishes. Synthemax II-coated and Laminin-521-coated dishes were found to be the most effective and efficient biomaterials for the cardiac differentiation of hPSCs according to the observation of hPSC-derived cardiomyocytes with high survival ratios, high beating colony numbers, a similar beating frequency to that of a healthy adult human heart, high purity levels (high cTnT expression) and longer sarcomere lengths similar to those of a healthy adult human heart.

Identification of atrial fibrillation associated genes and functional non-coding variants

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

Integrative characterization of G-Quadruplexes in the three-dimensional chromatin structure

DNA molecules are highly compacted in the eukaryotic nucleus where distal regulatory elements reach their targets through three-dimensional chromosomal interactions. G-quadruplexes, stable four-stranded non-canonical DNA structures, can change local chromatin organization through the exclusion of nucleosomes. However, the relationship between G-quadruplexes and higher-order genome organization remains unknown. Here, we found that G-quadruplexes are significantly enriched at boundaries of topological associated domains (TADs). Architectural protein occupancy, which plays critical roles in the formation of TADs, was highly correlated with the content of G-quadruplexes at TAD boundaries. Moreover, adjacent boundaries containing G-quadruplexes frequently interacted with each other because of the high enrichment of architectural protein binding sites. Similar to CCCTC-binding factor (CTCF) binding sites, G-quadruplexes also showed strong insulation ability in the separation of adjacent regions. Additionally, the insulation ability of CTCF binding sites and TAD boundaries was significantly reinforced by G-quadruplexes. Furthermore, G-quadruplex motifs on different strands were associated with the orientation of CTCF binding sites. These findings suggest a potential role for G-quadruplexes in loop extrusion. The enrichment of transcription factor binding sites (TFBSs) around regulatory elements containing G-quadruplexes led to frequent interactions between regulatory elements containing G-quadruplexes. Intriguingly, more than 99% of G-quadruplexes overlapped with TFBSs. The binding sites of CTCF and cohesin proteins were preferentially located surrounding G-quadruplexes. Accordingly, we proposed a new mechanism of long-distance gene regulation in which G-quadruplexes are involved in distal interactions between enhancers and promoters.

The design of a thermoresponsive surface for the continuous culture of human pluripotent stem cells

Commonly, stem cell culture is based on batch-type culture, which is laborious and expensive. We continuously cultured human pluripotent stem cells (hPSCs) on thermoresponsive dish surfaces, where hPSCs were partially detached on the same thermoresponsive dish by decreasing the temperature of the thermoresponsive dish to be below the lower critical solution temperature for only 30 min. Then, the remaining cells were continuously cultured in fresh culture medium, and the detached stem cells were harvested in the exchanged culture medium. hPSCs were continuously cultured for ten cycles on the thermoresponsive dish surface, which was prepared by coating the surface with poly(N-isopropylacrylamide-co-styrene) and oligovitronectin-grafted poly(acrylic acid-co-styrene) or recombinant vitronectin for hPSC binding sites to maintain hPSC pluripotency. After ten cycles of continuous culture on the thermoresponsive dish surface, the detached cells expressed pluripotency proteins and had the ability to differentiate into cells derived from the three germ layers in vitro and in vivo. Furthermore, the detached cells differentiated into specific cell lineages, such as cardiomyocytes, with high efficiency.

Concise Review: Genetic and Epigenetic Regulation of Cardiac Differentiation from Human Pluripotent Stem Cells

Human pluripotent stem cells (hPSCs), including both embryonic stem cells and induced pluripotent stem cells, are the ideal cell sources for disease modeling, drug discovery, and regenerative medicine. In particular, regenerative therapy with hPSC-derived cardiomyocytes (CMs) is an unmet medical need for the treatment of severe heart failure. Cardiac differentiation protocols from hPSCs are made on the basis of cardiac development in vivo. However, current protocols have yet to yield 100% pure CMs, and their maturity is low. Cardiac development is regulated by the cardiac gene network, including transcription factors (TFs). According to our current understanding of cardiac development, cardiac TFs are sequentially expressed during cardiac commitment in hPSCs. Expression levels of each gene are strictly regulated by epigenetic modifications. DNA methylation, histone modification, and noncoding RNAs significantly influence cardiac differentiation. These complex circuits of genetic and epigenetic factors dynamically affect protein expression and metabolic changes in cardiac differentiation and maturation. Here, we review cardiac differentiation protocols and their molecular machinery, closing with a discussion of the future challenges for producing hPSC-derived CMs. Stem Cells 2019;37:992-1002.

Dynamics of genome reorganization during human cardiogenesis reveal an RBM20-dependent splicing factory

Functional changes in spatial genome organization during human development are poorly understood. Here we report a comprehensive profile of nuclear dynamics during human cardiogenesis from pluripotent stem cells by integrating Hi-C, RNA-seq and ATAC-seq. While chromatin accessibility and gene expression show complex on/off dynamics, large-scale genome architecture changes are mostly unidirectional. Many large cardiac genes transition from a repressive to an active compartment during differentiation, coincident with upregulation. We identify a network of such gene loci that increase their association inter-chromosomally, and are targets of the muscle-specific splicing factor RBM20. Genome editing studies show that TTN pre-mRNA, the main RBM20-regulated transcript in the heart, nucleates RBM20 foci that drive spatial proximity between the TTN locus and other inter-chromosomal RBM20 targets such as CACNA1C and CAMK2D. This mechanism promotes RBM20-dependent alternative splicing of the resulting transcripts, indicating the existence of a cardiac-specific trans-interacting chromatin domain (TID) functioning as a splicing factory.

GOPHER: Generator Of Probes for capture Hi-C Experiments at high Resolution

BACKGROUND

Target enrichment combined with chromosome conformation capturing methodologies such as capture Hi-C (CHC) can be used to investigate spatial layouts of genomic regions with high resolution and at scalable costs. A common application of CHC is the investigation of regulatory elements that are in contact with promoters, but CHC can be used for a range of other applications. Therefore, probe design for CHC needs to be adapted to experimental needs, but no flexible tool is currently available for this purpose.

RESULTS

We present a Java desktop application called GOPHER (Generator Of Probes for capture Hi-C Experiments at high Resolution) that implements three strategies for CHC probe design. GOPHER's simple approach is similar to the probe design of previous approaches that employ CHC to investigate all promoters, with one probe being placed at each margin of a single digest that overlaps the transcription start site (TSS) of each promoter. GOPHER's simple-patched approach extends this methodology with a heuristic that improves coverage of viewpoints in which the TSS is located near to one of the boundaries of the digest. GOPHER's extended approach is intended mainly for focused investigations of smaller gene sets. GOPHER can also be used to design probes for regions other than TSS such as GWAS hits or large blocks of genomic sequence. GOPHER additionally provides a number of features that allow users to visualize and edit viewpoints, and outputs a range of files useful for documentation, ordering probes, and downstream analysis.

CONCLUSION

GOPHER is an easy-to-use and robust desktop application for CHC probe design. Source code and a precompiled executable can be downloaded from the GOPHER GitHub page at https://github.com/TheJacksonLaboratory/Gopher .