A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development.

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

CRISPR Activation Reverses Haploinsufficiency and Functional Deficits Caused by TTN Truncation Variants.

BACKGROUND: TTN truncation variants (TTNtvs) are the most common genetic lesion identified in individuals with dilated cardiomyopathy, a disease with high morbidity and mortality rates. TTNtvs reduce normal TTN (titin) protein levels, produce truncated proteins, and impair sarcomere content and function. Therapeutics targeting TTNtvs have been elusive because of the immense size of TTN, the rarity of specific TTNtvs, and incomplete knowledge of TTNtv pathogenicity. METHODS: We adapted CRISPR activation using dCas9-VPR to functionally interrogate TTNtv pathogenicity and develop a therapeutic in human cardiomyocytes and 3-dimensional cardiac microtissues engineered from induced pluripotent stem cell models harboring a dilated cardiomyopathy-associated TTNtv. We performed guide RNA screening with custom TTN reporter assays, agarose gel electrophoresis to quantify TTN protein levels and isoforms, and RNA sequencing to identify molecular consequences of TTN activation. Cardiomyocyte epigenetic assays were also used to nominate DNA regulatory elements to enable cardiomyocyte-specific TTN activation. RESULTS: CRISPR activation of TTN using single guide RNAs targeting either the TTN promoter or regulatory elements in spatial proximity to the TTN promoter through 3-dimensional chromatin interactions rescued TTN protein deficits disturbed by TTNtvs. Increasing TTN protein levels normalized sarcomere content and contractile function despite increasing truncated TTN protein. In addition to TTN transcripts, CRISPR activation also increased levels of myofibril assembly-related and sarcomere-related transcripts. CONCLUSIONS: TTN CRISPR activation rescued TTNtv-related functional deficits despite increasing truncated TTN levels, which provides evidence to support haploinsufficiency as a relevant genetic mechanism underlying heterozygous TTNtvs. CRISPR activation could be developed as a therapeutic to treat a large proportion of TTNtvs.

Novel Pan-ERR Agonists Ameliorate Heart Failure Through Enhancing Cardiac Fatty Acid Metabolism and Mitochondrial Function.

BACKGROUND: Cardiac metabolic dysfunction is a hallmark of heart failure (HF). Estrogen-related receptors ERRα and ERRγ are essential regulators of cardiac metabolism. Therefore, activation of ERR could be a potential therapeutic intervention for HF. However, in vivo studies demonstrating the potential usefulness of ERR agonist for HF treatment are lacking, because compounds with pharmacokinetics appropriate for in vivo use have not been available. METHODS: Using a structure-based design approach, we designed and synthesized 2 structurally distinct pan-ERR agonists, SLU-PP-332 and SLU-PP-915. We investigated the effect of ERR agonist on cardiac function in a pressure overload-induced HF model in vivo. We conducted comprehensive functional, multi-omics (RNA sequencing and metabolomics studies), and genetic dependency studies both in vivo and in vitro to dissect the molecular mechanism, ERR isoform dependency, and target specificity. RESULTS: Both SLU-PP-332 and SLU-PP-915 significantly improved ejection fraction, ameliorated fibrosis, and increased survival associated with pressure overload-induced HF without affecting cardiac hypertrophy. A broad spectrum of metabolic genes was transcriptionally activated by ERR agonists, particularly genes involved in fatty acid metabolism and mitochondrial function. Metabolomics analysis showed substantial normalization of metabolic profiles in fatty acid/lipid and tricarboxylic acid/oxidative phosphorylation metabolites in the mouse heart with 6-week pressure overload. ERR agonists increase mitochondria oxidative capacity and fatty acid use in vitro and in vivo. Using both in vitro and in vivo genetic dependency experiments, we show that ERRγ is the main mediator of ERR agonism-induced transcriptional regulation and cardioprotection and definitively demonstrated target specificity. ERR agonism also led to downregulation of cell cycle and development pathways, which was partially mediated by E2F1 in cardiomyocytes. CONCLUSIONS: ERR agonists maintain oxidative metabolism, which confers cardiac protection against pressure overload-induced HF in vivo. Our results provide direct pharmacologic evidence supporting the further development of ERR agonists as novel HF therapeutics.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes.

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Epigenomic and Transcriptomic Dynamics During Human Heart Organogenesis.

RATIONALE: There is growing evidence that common variants and rare sequence alterations in regulatory sequences can result in birth defects or predisposition to disease. Congenital heart defects are the most common birth defect and have a clear genetic component, yet only a third of cases can be attributed to structural variation in the genome or a mutation in a gene. The remaining unknown cases could be caused by alterations in regulatory sequences. OBJECTIVE: Identify regulatory sequences and gene expression networks that are active during organogenesis of the human heart. Determine whether these sites and networks are enriched for disease-relevant genes and associated genetic variation. METHODS AND RESULTS: We characterized ChromHMM (chromatin state) and gene expression dynamics during human heart organogenesis. We profiled 7 histone modifications in embryonic hearts from each of 9 distinct Carnegie stages (13-14, 16-21, and 23), annotated chromatin states, and compared these maps to over 100 human tissues and cell types. We also generated RNA-sequencing data, performed differential expression, and constructed weighted gene coexpression networks. We identified 177 412 heart enhancers; 12 395 had not been previously annotated as strong enhancers. We identified 92% of all functionally validated heart-positive enhancers (n=281; 7.5× enrichment; P<2.2×10-16). Integration of these data demonstrated novel heart enhancers are enriched near genes expressed more strongly in cardiac tissue and are enriched for variants associated with ECG measures and atrial fibrillation. Our gene expression network analysis identified gene modules strongly enriched for heart-related functions, regulatory control by heart-specific enhancers, and putative disease genes. CONCLUSIONS: Well-connected hub genes with heart-specific expression targeted by embryonic heart-specific enhancers are likely disease candidates. Our functional annotations will allow for better interpretation of whole genome sequencing data in the large number of patients affected by congenital heart defects.

A bipartite boundary element restricts UBE3A imprinting to mature neurons.

Angelman syndrome (AS) is a severe neurodevelopmental disorder caused by the loss of function from the maternal allele of UBE3A, a gene encoding an E3 ubiquitin ligase. UBE3A is only expressed from the maternally inherited allele in mature human neurons due to tissue-specific genomic imprinting. Imprinted expression of UBE3A is restricted to neurons by expression of UBE3A antisense transcript (UBE3A-ATS) from the paternally inherited allele, which silences the paternal allele of UBE3A in cis However, the mechanism restricting UBE3A-ATS expression and UBE3A imprinting to neurons is not understood. We used CRISPR/Cas9-mediated genome editing to functionally define a bipartite boundary element critical for neuron-specific expression of UBE3A-ATS in humans. Removal of this element led to up-regulation of UBE3A-ATS without repressing paternal UBE3A However, increasing expression of UBE3A-ATS in the absence of the boundary element resulted in full repression of paternal UBE3A, demonstrating that UBE3A imprinting requires both the loss of function from the boundary element as well as the up-regulation of UBE3A-ATS These results suggest that manipulation of the competition between UBE3A-ATS and UBE3A may provide a potential therapeutic approach for AS.

GPCRomics: GPCR Expression in Cancer Cells and Tumors Identifies New, Potential Biomarkers and Therapeutic Targets.

G protein-coupled receptors (GPCRs), the largest family of targets for approved drugs, are rarely targeted for cancer treatment, except for certain endocrine and hormone-responsive tumors. Limited knowledge regarding GPCR expression in cancer cells likely has contributed to this lack of use of GPCR-targeted drugs as cancer therapeutics. We thus undertook GPCRomic studies to define the expression of endoGPCRs (which respond to endogenous molecules such as hormones, neurotransmitters and metabolites) in multiple types of cancer cells. Using TaqMan qPCR arrays to quantify the mRNA expression of ∼340 such GPCRs, we found that human chronic lymphocytic leukemia (CLL) cells/stromal cells associated with CLL, breast cancer cell lines, colon cancer cell lines, pancreatic ductal adenocarcinoma (PDAC) cells, cancer associated fibroblasts (CAFs), and PDAC tumors express 50 to >100 GPCRs, including many orphan GPCRs (which lack known physiologic agonists). Limited prior data exist regarding the expression or function of most of the highly expressed GPCRs in these cancer cells and tumors. Independent results from public cancer gene expression databases confirm the expression of such GPCRs. We propose that highly expressed GPCRs in cancer cells (for example, GPRC5A in PDAC and colon cancer cells and GPR68 in PDAC CAFs) may contribute to the malignant phenotype, serve as biomarkers and/or may be novel therapeutic targets for the treatment of cancer.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development.

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

The Gq signalling pathway inhibits brown and beige adipose tissue.

Brown adipose tissue (BAT) dissipates nutritional energy as heat via the uncoupling protein-1 (UCP1) and BAT activity correlates with leanness in human adults. Here we profile G protein-coupled receptors (GPCRs) in brown adipocytes to identify druggable regulators of BAT. Twenty-one per cent of the GPCRs link to the Gq family, and inhibition of Gq signalling enhances differentiation of human and murine brown adipocytes. In contrast, activation of Gq signalling abrogates brown adipogenesis. We further identify the endothelin/Ednra pathway as an autocrine activator of Gq signalling in brown adipocytes. Expression of a constitutively active Gq protein in mice reduces UCP1 expression in BAT, whole-body energy expenditure and the number of brown-like/beige cells in white adipose tissue (WAT). Furthermore, expression of Gq in human WAT inversely correlates with UCP1 expression. Thus, our data indicate that Gq signalling regulates brown/beige adipocytes and inhibition of Gq signalling may be a novel therapeutic approach to combat obesity.

Mechanisms of cyclic AMP/protein kinase A- and glucocorticoid-mediated apoptosis using S49 lymphoma cells as a model system.

Cyclic AMP/protein kinase A (cAMP/PKA) and glucocorticoids promote the death of many cell types, including cells of hematopoietic origin. In wild-type (WT) S49 T-lymphoma cells, signaling by cAMP and glucocorticoids converges on the induction of the proapoptotic B-cell lymphoma-family protein Bim to produce mitochondria-dependent apoptosis. Kin(-), a clonal variant of WT S49 cells, lacks PKA catalytic (PKA-Cα) activity and is resistant to cAMP-mediated apoptosis. Using sorbitol density gradient fractionation, we show here that in kin(-) S49 cells PKA-Cα is not only depleted but the residual PKA-Cα mislocalizes to heavier cell fractions and is not phosphorylated at two conserved residues (Ser(338) or Thr(197)). In WT S49 cells, PKA-regulatory subunit I (RI) and Bim coimmunoprecipitate upon treatment with cAMP analogs and forskolin (which increases endogenous cAMP concentrations). By contrast, in kin(-) cells, expression of PKA-RIα and Bim is prominently decreased, and increases in cAMP do not increase Bim expression. Even so, kin(-) cells undergo apoptosis in response to treatment with the glucocorticoid dexamethasone (Dex). In WT cells, glucorticoid-mediated apoptosis involves an increase in Bim, but in kin(-) cells, Dex-promoted cell death appears to occur by a caspase 3-independent apoptosis-inducing factor pathway. Thus, although cAMP/PKA-Cα and PKA-R1α/Bim mediate apoptotic cell death in WT S49 cells, kin(-) cells resist this response because of lower levels of PKA-Cα and PKA-RIα subunits as well as Bim. The findings for Dex-promoted apoptosis imply that these lymphoma cells have adapted to selective pressure that promotes cell death by altering canonical signaling pathways.

A second trigeminal CGRP receptor: function and expression of the AMY1 receptor.

OBJECTIVE: The trigeminovascular system plays a central role in migraine, a condition in need of new treatments. The neuropeptide, calcitonin gene-related peptide (CGRP), is proposed as causative in migraine and is the subject of intensive drug discovery efforts. This study explores the expression and functionality of two CGRP receptor candidates in the sensory trigeminal system. METHODS: Receptor expression was determined using Taqman G protein-coupled receptor arrays and immunohistochemistry in trigeminal ganglia (TG) and the spinal trigeminal complex of the brainstem in rat and human. Receptor pharmacology was quantified using sensitive signaling assays in primary rat TG neurons. RESULTS: mRNA and histological expression analysis in rat and human samples revealed the presence of two CGRP-responsive receptors (AMY1: calcitonin receptor/receptor activity-modifying protein 1 [RAMP1]) and the CGRP receptor (calcitonin receptor-like receptor/RAMP1). In support of this finding, quantification of agonist and antagonist potencies revealed a dual population of functional CGRP-responsive receptors in primary rat TG neurons. INTERPRETATION: The unexpected presence of a functional non-canonical CGRP receptor (AMY1) at neural sites important for craniofacial pain has important implications for targeting the CGRP axis in migraine.

Proteomic and Metabolic Analyses of S49 Lymphoma Cells Reveal Novel Regulation of Mitochondria by cAMP and Protein Kinase A.

Cyclic AMP (cAMP), acting via protein kinase A (PKA), regulates many cellular responses, but the role of mitochondria in such responses is poorly understood. To define such roles, we used quantitative proteomic analysis of mitochondria-enriched fractions and performed functional and morphologic studies of wild-type (WT) and kin(-) (PKA-null) murine S49 lymphoma cells. Basally, 75 proteins significantly differed in abundance between WT and kin(-) S49 cells. WT, but not kin(-), S49 cells incubated with the cAMP analog 8-(4-chlorophenylthio)adenosine cAMP (CPT-cAMP) for 16 h have (a) increased expression of mitochondria-related genes and proteins, including ones in pathways of branched-chain amino acid and fatty acid metabolism and (b) increased maximal capacity of respiration on branched-chain keto acids and fatty acids. CPT-cAMP also regulates the cellular rate of ATP-utilization, as the rates of both ATP-linked respiration and proton efflux are decreased in WT but not kin(-) cells. CPT-cAMP protected WT S49 cells from glucose or glutamine deprivation, In contrast, CPT-cAMP did not protect kin(-) cells or WT cells treated with the PKA inhibitor H89 from glutamine deprivation. Under basal conditions, the mitochondrial structure of WT and kin(-) S49 cells is similar. Treatment with CPT-cAMP produced apoptotic changes (i.e. decreased mitochondrial density and size and loss of cristae) in WT, but not kin(-) cells. Together, these findings show that cAMP acts via PKA to regulate multiple aspects of mitochondrial function and structure. Mitochondrial perturbation thus likely contributes to cAMP/PKA-mediated cellular responses.

G Protein-Coupled Receptor (GPCR) Expression in Native Cells: "Novel" endoGPCRs as Physiologic Regulators and Therapeutic Targets.

G protein-coupled receptors (GPCRs), the largest family of signaling receptors in the human genome, are also the largest class of targets of approved drugs. Are the optimal GPCRs (in terms of efficacy and safety) currently targeted therapeutically? Especially given the large number (∼ 120) of orphan GPCRs (which lack known physiologic agonists), it is likely that previously unrecognized GPCRs, especially orphan receptors, regulate cell function and can be therapeutic targets. Knowledge is limited regarding the diversity and identity of GPCRs that are activated by endogenous ligands and that native cells express. Here, we review approaches to define GPCR expression in tissues and cells and results from studies using these approaches. We identify problems with the available data and suggest future ways to identify and validate the physiologic and therapeutic roles of previously unrecognized GPCRs. We propose that a particularly useful approach to identify functionally important GPCRs with therapeutic potential will be to focus on receptors that show selective increases in expression in diseased cells from patients and experimental animals.

Quantitative proteomics analysis of the cAMP/protein kinase A signaling pathway.

To define the proteins whose expression is regulated by cAMP and protein kinase A (PKA), we used a quantitative proteomics approach in studies of wild-type (WT) and kin- (PKA-null) S49 murine T lymphoma cells. We also compared the impact of endogenous increases in the level of cAMP [by forskolin (Fsk) and the phosphodiesterase inhibitor isobutylmethylxanthine (IBMX)] or by a cAMP analogue (8-CPT-cAMP). We identified 1056 proteins in WT and kin- S49 cells and found that 8-CPT-cAMP and Fsk with IBMX produced differences in protein expression. WT S49 cells had a correlation coefficient of 0.41 between DNA microarray data and the proteomics analysis in cells incubated with 8-CPT-cAMP for 24 h and a correlation coefficient of 0.42 between the DNA microarray data obtained at 6 h and the changes in protein expression after incubation with 8-CPT-cAMP for 24 h. Glutathione reductase (Gsr) had a higher level of basal expression in kin- S49 cells than in WT cells. Consistent with this finding, kin- cells are less sensitive to cell killing and generation of malondialdehyde than are WT cells incubated with H(2)O(2). Cyclic AMP acting via PKA thus has a broad impact on protein expression in mammalian cells, including in the regulation of Gsr and oxidative stress.

Increased expression of the pro-apoptotic protein BIM, a mechanism for cAMP/protein kinase A (PKA)-induced apoptosis of immature T cells.

The second messenger cAMP is proapoptotic for numerous cell types, but the mechanism for this proapoptotic action is not defined. Here, we use murine CD4(+)/CD8(+) S49 lymphoma cells and isolated thymocytes to assess this mechanism. In WT S49 cells, cAMP acts via protein kinase A (PKA) to induce G(1) phase cell cycle arrest and apoptosis. Treatment of WT and cAMP-Deathless (D-) S49 cells, which lack cAMP-promoted apoptosis, with the PKA agonist 8-(4-chlorophenylthio)-cAMP (CPT-cAMP) differentially regulates transcripts for numerous proapoptotic and antiapoptotic proteins. In contrast, kin-S49 cells (which lack PKA) show no cAMP-promoted changes in transcript expression. In this study, we use knockdown and overexpression approaches to define the role in cAMP/PKA-promoted apoptosis of the proapoptotic factor BIM (Bcl-2 interacting mediator of cell death), whose expression prominently increases in response to CPT-cAMP treatment of WT but not D- or kin- S49 cells. Conditional expression of BimL, one of the three major forms of Bim, increases apoptosis of WT, D-, and kin-S49 cells, whereas inhibition of cAMP-mediated induction of Bim isoforms by shRNAi attenuates CPT-cAMP-mediated apoptosis of WT S49 cells. Bim protein levels increase in subpopulations of CPT-cAMP-treated cells that undergo apoptosis. Thymic CD4(+)/CD8(+) cells isolated from Bim(-/-) mice corroborated the requirement of Bim expression for cAMP-promoted apoptosis. Thus, up-regulation of Bim appears to be an important determinant of cAMP/PKA-mediated apoptosis in immature T cells and may be a mechanism for such apoptosis in other cell types as well.